Аппроксимация
Параметры схемы (неопределенные пока коэффициентыa1, …,aK;c1, …,cK; b11, …,b1K, …,bK1, …,bKK), конечно, не произвольны. Их (все или часть) находят, прежде всего, из условий аппроксимации, получаемых из разложения (21) в ряд в точкеt=tn илиt =tn+1 с учетом (18) и его следствий
vt = f(t, v),
(т. е. рассматривая аппроксимацию на решениях (18)).
Поскольку в (21) при rkстоит множительτ, то в разложении правой части (21) нужно удерживать на один член ряда меньше, чем в левой. Индексnбудем опускать. Итак, левая часть (21):
Далее:
Если коэффициенты akвыбирать таким образом, чтобы
, (23)
то Подставляя эти разложения в (21) и группируя члены при одинаковых степеняхτ, получим (после деления на τ )
Очевидно, что при выполнении условия (вместе с (23))
(24)
обеспечивается 1-й порядок аппроксимации, при
(25)
(вместе с (23), (24)) — 2-й порядок аппроксимации. Условия 3-го порядка точности (вместе с (23) – (25)) есть
, (26)
а условия 4-го порядка (вместе с (23) – (26)):
,
(27)
.
Уравнения (23) – (27) составляют относительно неопределенных коэффициентов некоторую нелинейную систему. Очевидно, что привлекаемое число условий аппроксимации (т. е. порядок точности схемы рпри выполнении соответствующих условий гладкости) должно быть меньше числа отличных от нуля коэффициентов в (21), и соответствующая система должна быть разрешима. В частности, для явных схем, а для общих неявных схем.
- Численное интегрирование жестких системобыкновенных дифференциальных уравнений (оду)
- Жесткие оду
- Линейные однородные уравнения 1-го порядка
- Системы линейных однородных уравнений
- Пример: задача Коши для линейного однородного уравнения второго порядка
- Нелинейные жесткие уравнения
- Пример: сингулярно-возмущённая нелинейная система второго порядка
- Произвольная система нелинейных уравнений
- Примеры простейших разностных схем для жестких оду
- Способы построения схем
- Требования к численным методам решения жёстких систем оду
- Одношаговые методы типа Рунге–Кутты
- Алгоритм
- Аппроксимация
- Устойчивость
- Примеры схем Рунге–Кутты
- Линейные многошаговые схемы (методы типа Адамса)
- Алгоритм и аппроксимация
- Устойчивость
- Примеры линейных многошаговых схем
- Схемы для продолженных систем (схемы Обрешкова)
- Алгоритм и аппроксимация
- Устойчивость
- Контрольные вопросы
- Общие вопросы к лабораторным работам 1–3
- Схемы Рунге–Кутты (работа №1)
- Уравнение Ван-дер-Поля
- Система Ван-дер-Поля и траектории-утки
- Суточные колебания озона в атмосфере
- Уравнение Бонгоффера–Ван-дер-Поля
- Сингулярно-возмущенная система — модель двухлампового генератора Фрюгауфа
- Простейшая модель гликолиза
- Модель химических реакций Робертсона
- Модель дифференциации растительной ткани
- Задача e5
- Уравнение Релея
- Экогенетические модели
- Список литературы