Точность и устойчивость

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

Вывод основного уравнения. Пусть уравнение (1) разрешено относительно у `, имеем задачу

. (5.2)

Все численные методы решения задачи (2) в качестве основного блока имеют некую процедуру вычисления значения у_i+1, если известно у_i и некоторые другие параметры. Т.е. вначале задаётся точка (х₀, у₀), по ней вычисляем точку (х₁, у₁), затем тем же способом вычисляем (х₂, у₂), и т.д. до конца таблицы.

Пусть у(х) – решение задачи (2). Подставим его в формулу (2) – получим тождество y`(x) º f(x, y(x)). Проинтегрируем его в пределах от х_i до x_i+1 = х_i + h:

. Отсюда имеем: . Или

. (5.3)

Формула (3) – интегральная форма задачи (2). Из этого уравнения надо найти у_i+1. Для построения численного метода необходимо воспользоваться какой-нибудь аппроксимацией интеграла (3) (квадратурной формулой). Получим основное уравнение

у_i+1– у_i » Ф(х_i, h, y_i+1–k,…, y_i, y_i+1), (5.4)

где Ф(…) – функция, являющаяся аппроксимацией интеграла в правой части формулы (3). Тип конкретного численного метода в основном определяется

а) выбранной аппроксимацией, т.е. видом функции Ф;

б) способом решения алгебраического уравнения (4).

Основные классы численных методов

1. Если в (4) k =1, то имеем уравнение у_i+1– у_i » Ф(х_i, h, y_i, y_i+1), т.е. не используются более ранние значения функции, чем y_i. В этом случае метод называется одношаговым. Если k > 1, т.е. надо задать y_i–1, y_i–2 и т.д., то метод будет многошаговым.

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

2. Если Ф не зависит от y_i+1, то решить уравнение (4) очень просто:

у_i+1 » у_i + Ф(х_i, h, y_i+1–k,…, y_i, y_i+1),

т.е. получим явную формулу для y_i+1. Соответствующие методы называются явными. В противоположность им, методы, в которых функция Ф зависит от у_i+1, называются неявными. При их реализации при каждом i возникает необходимость решения относительно у_i+1 некоторого нелинейного уравнения.

Три простейших метода. Получим формулы трёх простейших численных методов, используя простые (не составные) формулы численного интегрирования для аппроксимации интеграла (5.3).

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

= Ф(х_i, h, у_i).

Функция Ф не зависит от у_i+1 и в ней k = 1. Следовательно, получаем явный одношаговый метод

у_i+1 = у_i + h f(х_i, у_i), (5.5)

называемый явным методом Эйлера.

Пример 5.3. Требуется решить задачу Коши

.

Здесь f(x, y) = y² + 1. Аналитическое решение задачи y = tg x. Для этого уравнения формула (5) имеет вид

у_i+1 = у_i + h [(у_i)² + 1].

2) Воспользуемся формулой трапеций (4.1), согласно которой интеграл (5.3) имеет следующий вид:

.

При такой аппроксимации интеграла правая часть равенства (4) зависит от у_i+1 и поэтому имеем неявный одношаговый метод вида

у_i+1 – у_i = [ f (x_i, у_i) + f (x_i+1, у_i+1) ], (5.6)

называемый неявным методом Эйлера. Согласно этому методу, в примере 5.3 для искомого значения у_i+1 будем иметь квадратное уравнение,

,

решив которое получим расчетную формулу:

.

В общем случае уравнение (6) может не иметь аналитического решения. Тогда для организации расчетов по методу (6) можно воспользоваться каким-нибудь численным алгоритмом решения уравнений для определения значения у_i+1 при каждом I (например, методом обратной интерполяции).

3) Другой способ реализации метода (6) без непосредственного решения уравнения относительно у_i+1 состоит в соединении явного и неявного методов. Вместо точного значения у_i+1 в правую часть равенства (6) подставляется его приближенное значение, найденное по формуле явного метода (5), а затем применяется формула (6). В результате получаем следующие расчетные формулы:

= у_i + h f(х_i, у_i); у_i+1 = у_i + [ f(x_i, у_i) + f(x_i+1, ) ]. (5.7)

Такая идея соединения явного и неявного методов носит название схемы прогноза и коррекции. На первом шаге по явной формуле получается величина , называемая прогнозом, которая затем корректируется с использованием формулы неявного метода. Полученное значение у_i+1 считается окончательным искомым значением функции у(х) в i + 1-й точке. Метод (7) является одним из простейших вариантов такой схемы и называется модифицированным методом Эйлера или методом Эйлера-Коши.

Подставив левую формулу (7) в правую, получим единую расчетную формулу, которая для примера 5.3 будет иметь вид

.

Дата добавления: 2014-01-07