29. Задача численного решения обыкновенных дифференциальных уравнений. Задача Коши и граничные задачи.

Рассмотрим решение обыкновенного дифференциального уравнения первого порядка:

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

Решение дифференциального уравнения ищется на интервале:

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

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

Ясно, что при n=1 можно говорить только о задачи Коши.

Примеры постановки задачи Коши:

Примеры краевых задач:

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

Численные методы решения задачи Коши для ОДУ первого порядка

Постановка задачи. Найти решение ОДУ первого порядка

 на отрезке  при условии 

При нахождении приближенного решения будем считать, что вычисления проводятся с расчетным шагом , расчетными узлами служат точки  промежутка [x0, xn].

Целью является построение таблицы

xi	x0	x1	…	xn
yi	y0	y1	…	yn

т.е. ищутся приближенные значения y в узлах сетки.

Интегрируя уравнение на отрезке , получим

Вполне естественным (но не единственным) путем получения численного решения является замена в нем интеграла какой–либо квадратурной формулой численного интегрирования. Если воспользоваться простейшей формулой левых прямоугольников первого порядка

,

30. Решение задачи Коши с помощью формулы Тейлора.

Чтобы построить ряд Тейлора в окрестности точки, нужно найти значения производных функции в этих точках:

и т.д.

В результате решения на интервале можно представить в виде конечного числа слагаемых ряда Тейлора:

Недостаток этого метода:

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

Большую группу методов, свободную от этих недостатков, образуют методы Рунга-Кутта. К этим методам относятся методы Эйлера и Адама.

Метод Эйлера представляет собой аппроксимации производной разделенной разностью первого порядка.

Отсюда получим разностное уравнение:

– интервал интегрирования дифференциального уравнения.

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

Дифференциальное уравнение может быть представлено эквивалентным интегральным уравнением:

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

Простейшая формула – формула прямоугольника.

Формула прямоугольника на интервале может быть записана в виде:

В этом случае имеет место формула Эйлера.

Точность формулы Эйлера имеет величину, пропорциональную квадрату шага дискретизации.

Более точная формула получится с помощью формулы трапеции:

Тогда расчетная формула принимает вид:

Эта формула является нелинейной по искомому значению функции в конце интервала.

Точность этой формулы имеет величину второго порядка .

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

цикл

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

Алгоритм записывается:

цикл

На практике:

Аналогичным образом строится формула Адамса, основывается на методе Рунга-Кутта.

При этом используются значения функции не только в концах интервала, но и в точках, лежащих внутри интервала.