Студопедия Главная Случайная страница Обратная связь

Разделы: Автомобили Астрономия Биология География Дом и сад Другие языки Другое Информатика История Культура Литература Логика Математика Медицина Металлургия Механика Образование Охрана труда Педагогика Политика Право Психология Религия Риторика Социология Спорт Строительство Технология Туризм Физика Философия Финансы Химия Черчение Экология Экономика Электроника

Методы Рунге-Кутта





Например, в одном из усовершенствований метода Эйлера, который также называют методом Рунге-Кутта второго порядка, переход осуществляют в два этапа по формулам:

Zi = Yi +h/2*f(Xi,Yi) Yi+1=Yi+h*f(Xi+h/2,Zi).

При этом получается погрешность порядка h2.

А самый распространенный на практике метод - метод Рунге-Кутта четвертого порядка, в котором точность имеет порядок величины h4, а переход к следующей точке осуществляется с помощью четырех вспомогательных величин:

k1= h*f(Xi,Yi)

k2= h*f(Xi+h/2,Yi+k1/2)

k3= h*f(Xi+h/2,Yi+k2/2)

k4= h*f(Xi+h,Yi+k3)

После вычисления этих вспомогательных величин переход от точки (Xi,Yi) к следующей точке осуществляется по формуле Yi+1=Yi+1/6*(k1+2k2+2k3+k4).

Упражнение 4.5. Выяснить геометрический смысл перехода к следующей точке по формулам усовершенствованного метода Эйлера.

Оценка точности в методах Рунге-Кутта второго и четвертого порядков на практике производится с помощью метода двойного счета, сформулированного в предыдущем параграфе.

Упражнение 4.6. Выписать и объяснить формулы оценки точности в методах Рунге-Кутта второго и четвертого порядков.

Поясним происхождение формул в методах Рунге-Кутта. Для получения закона вычисления значения Y(x) в каждой следующей точке поступают приблизительно так: выписывают разложение неизвестной функции в ряд Тейлора в точке Xi, как мы проделывали это выше, затем берут несколько первых членов этого разложения, и преобразуют полученную формулу Тейлора. После подставления Xi вместо переменной X и получают окончательное правило перехода к следующей точке.

Легко убедиться, что при выписывании разложения в ряд Тейлора только до линейного члена и подстановки значения X мы получим формулу метода Эйлера, т.е. и он является частным случаем общих методов Рунге-Кутта.

Пример 4.3. Применим формулы разобранных методов для нахождения значения функции Y(x) в точке 1, если эта функция удовлетворяет уравнению y'=y с начальным условием y(0)=1.

При решении методом Эйлера (n=2) имеем:

Y1=Y0+hf(X0,Y0)=1+0.5f(0,1)=1.5, Y(b)=Y2=Y1+hf(X1,Y1)=1.5+0.5f(0.5,1.5)=2.25.

При решении методом Рунге-Кутта 2-го порядка (n=1) имеем:

Z0=Y0+h/2*f(X0,Y0)=1+0.5f(0,1)=1.5, Y(b)=Y1=Y0+hf(X0+h/2,Z0)=1+f(0.5,1.5)=2.5.

При решении методом Рунге-Кутта 4-го порядка (n=1) имеем:

k1=h*f(X0,Y0)= 1*f(0,1)= 1

k2=h*f(X0+h/2,Y0+k1/2)= 1*f(0.5,1.5)= 1.5

k3=h*f(X0+h/2,Y0+k2/2)=1*f(0.5,1.75)=1.75

k4=h*f(X0+h,Y0+k3)= 1*f(1,2.75)=2.75

Y(b)=Y1=Y0+1/6*(k1+2k2+2k3+k4)=1+1/6*(1+3+3.5+2.75)=65/24

Упражнение 4.7. Найти точный ответ и сравнить погрешности, полученные при решении тремя различными методами.

Контрольные вопросы.

1. Какие типы приближенных методов решения обыкновенных дифференциальных уравнений Вы знаете? Назовите по одному примеру каждого типа.

2. В чем суть метода Пикара? Объясните происхождение рекуррентной формулы метода.

3. В чем суть метода разложения функции Y(x) в ряд?

4. В чем суть метода Эйлера? Поясните графически.

5. Какова общая схема численных методов решения дифференциальных уравнений первого порядка?

6. Каков порядок точности при решении дифференциальных уравнений методами Эйлера, Рунге-Кутта второго и четвертого порядков?

7. Каким образом на практике следят за точностью при решении дифференциальных уравнений методами Рунге-Кутта второго и четвертого порядков? Обоснуйте этот способ.

Содержание лабораторной работы.

Постановка задачи: По заданному обыкновенному дифференциальному уравнению на фиксированном отрезке и значению искомой функции в левом конце определить значение в правом конце с требуемой точностью.

Предварительная работа.

1. Для своего уравнения найти дома точное решение в заданных точках.

2. Найти методом Пикара третье приближение к решению своего уравнения, подставить заданные точки и найти погрешность.

3. С помощью метода разложения в ряд найти для своей задачи ответы с точностью 0.01.

4. Нарисовать для своей задачи с помощью метода Эйлера 5 звеньев ломаной, дающей представление об интегральной кривой.

Порядок работы:

1.Ответить на вопросы контролирующей программы.

2.Ввести в ЭВМ и отладить программы для вычисления ответа тремя способами численного решения уравнений: методами Рунге-Кутта 1-го, 2-го и 4-го порядков. Отладку производить на уравнении y'=y с начальным условием y(0)=1 и правым концом отрезка, равным 1.

3.Исполнить программу для своего варианта и записать ответы.

4.Дополнить программу вычисления по формуле Рунге-Кутта 4-го порядка так, чтобы по введенному e она с помощью метода двойного счета выдавала результат с требуемой точностью.

5.Оформить и сдать работу.

ОТЧЕТ должен содержать

1. название и цель работы,

2. домашнее исследование своей задачи методами Пикара, Эйлера и разложения в ряд.

3. тексты программ для всех трех методов,

4. ответы для своего варианта, точное аналитическое решение задачи.

 

 







Дата добавления: 2015-08-27; просмотров: 1242. Нарушение авторских прав; Мы поможем в написании вашей работы!




Функция спроса населения на данный товар Функция спроса населения на данный товар: Qd=7-Р. Функция предложения: Qs= -5+2Р,где...


Аальтернативная стоимость. Кривая производственных возможностей В экономике Буридании есть 100 ед. труда с производительностью 4 м ткани или 2 кг мяса...


Вычисление основной дактилоскопической формулы Вычислением основной дактоформулы обычно занимается следователь. Для этого все десять пальцев разбиваются на пять пар...


Расчетные и графические задания Равновесный объем - это объем, определяемый равенством спроса и предложения...

Тема 5. Анализ количественного и качественного состава персонала Персонал является одним из важнейших факторов в организации. Его состояние и эффективное использование прямо влияет на конечные результаты хозяйственной деятельности организации.

Билет №7 (1 вопрос) Язык как средство общения и форма существования национальной культуры. Русский литературный язык как нормированная и обработанная форма общенародного языка Важнейшая функция языка - коммуникативная функция, т.е. функция общения Язык представлен в двух своих разновидностях...

Патристика и схоластика как этап в средневековой философии Основной задачей теологии является толкование Священного писания, доказательство существования Бога и формулировка догматов Церкви...

Типовые ситуационные задачи. Задача 1.У больного А., 20 лет, с детства отмечается повышенное АД, уровень которого в настоящее время составляет 180-200/110-120 мм рт Задача 1.У больного А., 20 лет, с детства отмечается повышенное АД, уровень которого в настоящее время составляет 180-200/110-120 мм рт. ст. Влияние психоэмоциональных факторов отсутствует. Колебаний АД практически нет. Головной боли нет. Нормализовать...

Эндоскопическая диагностика язвенной болезни желудка, гастрита, опухоли Хронический гастрит - понятие клинико-анатомическое, характеризующееся определенными патоморфологическими изменениями слизистой оболочки желудка - неспецифическим воспалительным процессом...

Признаки классификации безопасности Можно выделить следующие признаки классификации безопасности. 1. По признаку масштабности принято различать следующие относительно самостоятельные геополитические уровни и виды безопасности. 1.1. Международная безопасность (глобальная и...

Studopedia.info - Студопедия - 2014-2025 год . (0.011 сек.) русская версия | украинская версия