Выражение ошибки численного решения через изменения коэффициентов дифференциального уравнения


 

Сравнивая точное решение (5) с приближенным численным решением (7) можно видеть, что процедура вычисления на шаге h для них одинакова: предыдущее значение X(t) умножается слева на некоторую матрицу. Однако в точном решении X(t) умножается на матрицу полной экспоненты eAh, а в приближенном - на матрицу усеченной экспоненты еAh.

Зададимся вопросом: для какого, измененного по отношению к (1), уравнения приближенное численное решение X(t+h) из(7) является точным?

Измененное уравнение (1) запишем следующим образом:

 

(2.1.10)

 

Его точное решение на шаге h

 

X(t+h) = e(A+dA)hX(t) (2.1.11)

 

Приравняв выражения (2) и (11), можно получить :

 

e(A+dA)h= eAh (2.1.12)

 

Это уравнение определяет матрицу dA из уравнения (10).

Матричную экспоненту (12) представим в виде:

 

e(A+dA)h= eAhedAh

 

Отсюда и из (6) можно получить:

 

edAh– I = - e-AhdeAh (2.1.13)

 

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

dAh - (Ah)r+1/(r+1)! (2.1.14)

 

Соотношения (13) и (14) позволяют вычислить матрицу dA, входящую в измененное уравнение (10). Точное решение этого уравнения совпадает с приближенным численным решением уравнения (1) методом r-того порядка. Таким образом, матрицу dA можно использовать в качестве характеристики погрешностей численного решения через изменение коэффициентов дифференциального уравнения (1).

Используя (14) можно выразить локальную ошибку численного решения (9) через dA:

 

dX(t+h) = - dAhX(t)

 

Если по условиям решаемой задачи можно сформулировать допустимые изменения коэффициентов dA, то соотношения (13) и (14) позволят осуществить предварительный выбор порядка метода r и шага дискретизации h.



Дата добавления: 2021-10-28; просмотров: 342;


Поиск по сайту:

Воспользовавшись поиском можно найти нужную информацию на сайте.

Поделитесь с друзьями:

Считаете данную информацию полезной, тогда расскажите друзьям в соц. сетях.
Poznayka.org - Познайка.Орг - 2016-2024 год. Материал предоставляется для ознакомительных и учебных целей.
Генерация страницы за: 0.008 сек.