При численном решении дифференциальных уравнений
неизбежно возникают ошибки, приводящие к некоторому отличию результатов численного
решения от истинных значений искомой функции, причём существует два типа источников этих
ошибок. Как уже отмечалось в разделе "Понятие порядка аппроксимации", одним источником
ошибок является сама аппроксимация, т.е. замена производных в исходных дифференциальных
уравнениях конечными разностями. Другой источник ошибок связан с погрешностью вычислений.
Ошибки могут возникать при неточности вычислений как самой разностной схемы, так и
начальных и граничных условий. В зависимости от особенностей вычислительного алгоритма эти
ошибки в процессе расчёта могут затухать или возрастать. Если они не возрастают, то говорят,
что разностная схема устойчива, если же возрастают, то - неустойчива.
Тогда разностную схему кратко можно представить в следующем виде:
Теорема 1. Разностная схема (3.2) с линейным оператором Lh устойчива, если при любом ![]() ![]()
Теорема 2. Пусть разностная схема (3.2) аппроксимирует дифференциальное уравнение (3.1) с порядком аппроксимации ![]()
Таким образом, сходимость решения разностной схемы к решению исходного уравнения имеет место только при выполнении двух требований: 1. разностная схема должна аппроксимировать исходное уравнение; 2. разностная схема должна быть устойчива. |