Лекция_11_ИЭТ 2011 Задача Коши

Методы
вычислений. ИЭТ. 4-й семестр. Лекция 11.
Стр. 7

МЕТОДЫ ВЫЧИСЛЕНИЙ

ИЭТ, 4-й семестр,
22

Лекция 11

Численные методы
решения задачи Коши

для систем
обыкновенных дифференциальных уравнений

Метод Эйлера, методы Рунге-Кутты.
Многошаговые методы Адамса

Основные понятия
и обозначения.

Обыкновенным
дифференциальным уравнением первого
порядка называется уравнение вида
,
где— известная функция трех переменных,
определенная в области,x
— независимая переменная,
,
неизвестная функция,— ее производная. Функцияназываетсярешением
дифференциального уравнения,
если
она непрерывно дифференцируема на
и
при всех
.
График
решения дифференциального уравнения
называют интегральной
кривой

дифференциального уравнения. В дальнейшем
будем рассматривать обыкновенные
дифференциальные уравнения в нормальной
форме
, т.е.
уравнения, разрешенные относительно
производной.:.(1)

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

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

Справедлива
следующая теорема существования и
единственности

Теорема. Если
функция
и ее частная производнаянепрерывны в области,,
то на некотором интервалесуществует единственное решение
уравнения,
удовлетворяющее начальному условию.

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

Обыкновенным
дифференциальным уравнением
-го
порядка
называется
уравнение вида

,
где
— известная функцияпеременных, определенная в области,— независимая переменная,
неизвестная функция,— порядок уравнения. Функцияназываетсярешением
дифференциального
уравнения,
если она
раз непрерывно
дифференцируема на
и
для всех
.
В дальнейшем
будем рассматривать обыкновенные
дифференциальные уравнения в нормальной
форме:.

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

,
,
…,,
при определенных ограничениях на правую
часть уравнения имеет единственное
решение. Справедлива следующая теорема.

Теорема.
Если функция
,,
и ее частные производные,непрерывны в области,,
то на некотором интервалесуществует единственное решение
уравнения,
удовлетворяющее начальным условиям:

,
,
…,,.

Задача Коши для
дифференциального уравнения
-го
порядка

,

,,
…,

может быть сведена
к задаче Коши для системы
го
порядка
(
системы n
дифференциальных уравнений 1-го порядка),
которая в векторной форме имеет вид
,,
где,.

Обусловленность
задачи Коши.

В силу
приведенных выше теорем существования
и единственности решения задачи Коши
и известной из курса дифференциальных
уравнений непрерывной зависимости
решения задачи Коши от начальных данных,
эта задача корректна. Однако применять
численные методы решения задачи имеет
смысл только к хорошо
обусловленным
задачам.
В инженерных и научных расчетах достаточно
часто возникают плохо
обусловленные задачи.

Пример плохо
обусловленной задачи Коши.
Рассмотрим
задачу Коши
,,.
Легко видеть, что общее решение задачи,
а решение задачи Кошии.
Однако, если начальные данные содержат
малую погрешность,,
тои.
Здесь важно понимать, что плохая
обусловленность задачи проявляется на
относительно длинном промежутке
интегрирования уравнения.

Модельная задача.
Полезно рассмотреть модельную
задачу


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

Задача Коши на
конечном отрезке.

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

Вычтем из уравнения
(1) уравнение (2) и воспользуемся формулой
конечных приращений Лагранжа:

где
некоторое промежуточное значение. В
результате получим, что погрешностьудовлетворяет дифференциальному
уравнению:

и начальному
условию
.
Решение этой задачи Коши выражается
формулой

.Таким
образом, величинаиграет в задаче Коши роль коэффициента
роста ошибки. Заметим, что знак производнойоказывает существенное влияние на
поведение погрешности.
Если,
то величина,
а вместе с ней и модуль погрешности
монотонно возрастают. При этом
соответствующие интегральные кривые
расходятся. ( см. рис.) Иначе ведет себя
погрешность в случае.
Здесьис ростоммонотонно убывают, а соответствующие
интегральные кривые сближаются. Ошибка,
внесенная в начальное значение, имеет
тенденцию к затуханию. В случае, когда
производнаянезнакопостоянна, поведение погрешности
может быть более сложным.

Постановка задачи
численного решения задачи Коши.

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

Численное решение
задачи Коши состоит в построении таблицы
приближенных значений
решенияв узлах сетки,.
Если,,
то сетка называетсяравномерной,
шаг сетки.

Как уже говорилось,
задача Коши для любого дифференциального
уравнения
-го
порядка

сводится к задаче
Коши для системы
дифференциальных уравнений 1-го порядка.Численное
решение задачи Коши для системы

дифференциальных уравнений состоит в
построении таблицы приближенных значений
компонентвектора решенияв точках.

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

Классификация
методов численного решения задачи Коши
.
Численный
метод решения задачи Коши называется
одношаговым,
если для вычисления решения в точке
используется информация о решении
только в точке.
В противном случае метод называетсямногошаговым.

Если при этом
,
т.е. решение в точкеявно выражается через значения в
«предыдущих» точках,то
метод относится кявным
методам.
Если же
,
то методнеявный.

Для оценки
погрешности метода на одном шаге сетки
обычно точное решение раскладывают по
формуле Тейлора в окрестности узла
:

.

Если значение
в точкеизвестно, то в силу равенствазначение производнойтакже можно считать известным. Для того
чтобы вычислить производныеболее высокого порядка, входящие в
формулу, продифференцируем равенствопоt
, используя правило дифференцирования
сложной функции. Получим

и т.д.

Если расчетные
формулы численного метода согласуются
с разложением по формуле Тейлора до
членов порядка
,
то числоназываетсяпорядком
метода.

Построение
численных методов решения задачи Коши.
Основной
этап численного решения задачи Коши
состоит в построении дискретной
задачи Коши.

Прежде чем описывать
общий подход, построим простейшую
дискретную задачу.

Пусть
.
Обозначими запишем уравнениево всех узлах сетки,
заменяя точное значение производнойзначением
разностной производной.
Подставив разностную производную в
уравнение, получим системуN
линейных
уравнений
,относительно
N
неизвестных
(значение
известно из начальных условий). Эта
система уравнений и представляет собой
один из вариантовдискретной
задачи Коши.
Положив
,,
получимформулу
Эйлера
численного
решения
задачи Коши
приближенное решение
,
вычисленноеметодом
Эйлера.

Итак, формулы
— расчетные формулы метода Эйлера
приближенного решения задачи Коши.

Метод
Эйлера допускает простую геометрическую
интерпретацию. Пусть известна точка
интегральной кривой уравнения
.Касательная
к интегральной кривой уравнения,
проходящая через эту точку определяется
уравнением
.
Следовательно, вычисленная методом
Эйлера точка
,,,
лежит на этой касательной. Таким образом,
после выполненияN
шагов неизвестная интегральная кривая
заменяется ломаной линией( ломаной
Эйлера), для которой угловой коэффициент
звена равен значению.

Метод
Эйлера представляет явный одношаговый
метод. Для него погрешность аппроксимации
имеет вид
.

Оценка погрешности.
Для метода
Эйлера справедлива такая оценка
глобальной погрешности:.

Влияние
вычислительной погрешности.
Оценивая
метод Эйлера необходимо учитывать, что
при его реализации на ЭВМ неизбежно
возникнут ошибки округления. В результате
фактически вычисляемые значения
будут удовлетворять соотношению.
Величиныучитывают вклад погрешностей округления.
Оценка погрешности фактически вычисляемых
значений:.
Оказывается, что полная погрешность
убывает только лишь при уменьшении шагадо некоторого значения.

Вернемся к построению
дискретной задачи Коши. Построение
дискретной задачи заключается в
следующем. Сначала строим на
сетку.
В каждом узле сетки запишем дискретный
аналог уравнения, заменяя производнуюкакой-либо разностной производной —,
а правую часть заменяем какой-либо ее
аппроксимацией —,
построенной на тех же узлах.

Полученная система
N уравнений
относительно N
неизвестных
значений
и естьдискретная
задача Коши.

В частности,
дискретная задача Коши для явного
одношагового метода имеет вид

,
или, что то же самое
.

Приближенное
решение задачи Коши в узлах сетки
полагаем равным вычисленным значениям
:.

Основной
характеристикой дискретной задачи Коши

является степень
аппроксимации.

Например, об
одношаговой явной дискретной задаче
Коши говорят, что дискретное уравнение
аппроксимирует
дифференциальное уравнение
с порядком p
,
если
,
константаC
не зависит
от h.

Величина
называетсяпогрешностью
аппроксимации
.

Говорят, что
численный метод решения задачи Коши
сходится с
порядком p
,
если
константаC
не зависит
от h.
Величина
называетсялокальной
погрешностью метода.

Численный метод
решения задачи Коши называется
устойчивым,
если для всех достаточно малых h
справедливо
неравенство
.
Здесь константане зависит отh,
а звездочкой помечены приближенные
решения возмущенной задачи дискретной
задачи Коши.

Доказано, что для
устойчивой дискретной задачи Коши
порядок аппроксимации и порядок
сходимости совпадают. Этот порядок
называется порядком
метода.

Легко видеть, что
описанный выше метод Эйлера является
устойчивым методом первого порядка, а
локальная погрешность метода Эйлера —
величина
.

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

Решение задачи
Коши методом Рунге-Кутты

Методом Рунге—Кутты
обычно называют одношаговый метод
четвертого порядка, относящийся к
широкому классу методов Рунге—Кутты.
В этом методе величины
вычисляются по следующим формулам:

Построить
расчетные формулы метода Рунге-Кутты
можно различными способами. Например,
можно использовать следующий прием.
Положив
,
имеем.
Затем, заменив интеграл в правой части
специально построенной квадратурной
формулой (вычислить интеграл иначе
невозможно, поскольку он зависит от
неизвестного искомого решенияy(x)).

Если заменить
интеграл формулой левых прямоугольников
,получим метод
Эйлера.

Если заменить
интеграл формулой трапеций ,получим
.
Этот метод имеет второй порядок точности,
но является неявным. Заменим значениена «предсказываемое» методом Эйлера.
Получаем метод, который называют методом
Эйлера-Коши:

.

Метод Эйлера-Коши
относятся к классу методов
прогноза и коррекции
.

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

Метод Рунге-Кутты
устойчив. Локальная погрешность метода
равна
.

Практически оценить
величину
достаточно сложно. При оценке погрешности
обычно используютправило
Рунге
. Для
этого сначала проводят вычисления с
шагом
,
а затем — с шагом.
Если— приближение, вычисленное с шагом,
а— с шагом,
то справедлива оценка.

За оценку погрешности
с шагом
принимают величину.

Используя оценку
по правилу Рунге, можно строить адаптивные
алгоритмы Рунге-Кутты. В адаптивных
алгоритмах на каждом очередном шаге
производят половинное деление шага до
тех пор, пока
не
станет меньше заданной погрешности
решения.

Leave a Comment