Метод конечных разностей или метод сеток

http://ad.cctpu.edu.ru/Math_method/math/12_files/image123.gif

равномерно сходится к точному с погрешностью http://ad.cctpu.edu.ru/Math_method/math/12_files/image125.gif при c="images/referats/3126/image062.png" alt="http://ad.cctpu.edu.ru/Math_method/math/12_files/image127.gif">

Таким образом, схема (2.28), (2.29) дает приближенное решение краевой задачи, но точность ее весьма мала. Это связано с тем, что аппроксимация производной

http://ad.cctpu.edu.ru/Math_method/math/12_files/image129.gif

имеет низкий порядок точности − погрешность этой аппроксимации

http://ad.cctpu.edu.ru/Math_method/math/12_files/image131.gif

Более точную разностную схему можно получить, если при переходе от линейной краевой задачи к конечно-разностным уравнениям воспользоваться центральными формулами для производных:

http://ad.cctpu.edu.ru/Math_method/math/12_files/image133.gif, (2.42)

http://ad.cctpu.edu.ru/Math_method/math/12_files/image135.gif, (2.43)

i=1, 2, ., n.

Погрешность формулы (2.42) выражается так:

http://ad.cctpu.edu.ru/Math_method/math/12_files/image137.gif

то есть формула (2.42) имеет второй порядок точности относительно шага сетки h. Подставляя выражения (2.42), (2.43) в задачу (2.24), (2.25) и выполняя некоторые преобразования, получим следующую систему:

http://ad.cctpu.edu.ru/Math_method/math/12_files/image139.gif(2.44)

Где http://ad.cctpu.edu.ru/Math_method/math/12_files/image141.gif.

Система (2.44) снова трехдиагональная и ее решение также можно получить методом прогонки. Его алгоритм здесь будет выглядеть так. Сначала находят коэффициенты

http://ad.cctpu.edu.ru/Math_method/math/12_files/image143.gif (2.45)

Затем определяют коэффициенты http://ad.cctpu.edu.ru/Math_method/math/12_files/image145.gif по следующим рекуррентным формулам:

http://ad.cctpu.edu.ru/Math_method/math/12_files/image147.gif (2.46)

Обратный ход начинается с нахождения http://ad.cctpu.edu.ru/Math_method/math/12_files/image105.gif:

http://ad.cctpu.edu.ru/Math_method/math/12_files/image149.gif (2.47)

После этого находим http://ad.cctpu.edu.ru/Math_method/math/12_files/image151.gifпо формулам:

http://ad.cctpu.edu.ru/Math_method/math/12_files/image153.gif, (2.48)

http://ad.cctpu.edu.ru/Math_method/math/12_files/image155.gif. (2.49)

Относительно схемы (2.44) можно также доказать, что она имеет единственное решение при

http://ad.cctpu.edu.ru/Math_method/math/12_files/image157.gifи http://ad.cctpu.edu.ru/Math_method/math/12_files/image159.gif,

и это решение может быть найдено описанным методом прогонки. Кроме того, для схемы (2.44) имеет место

Теорема

Пусть решение граничной задачи (2.24), (2.25) единственно и непрерывно дифференцируемо на [a, b] до четвертого порядка точности включительно. Если выполняются условия

http://ad.cctpu.edu.ru/Math_method/math/12_files/image161.gif, http://ad.cctpu.edu.ru/Math_method/math/12_files/image159.gif, http://ad.cctpu.edu.ru/Math_method/math/12_files/image164.gif

то схема (2.44) будет равномерно сходиться к решению задачи (2.24), (2.25) с погрешностью http://ad.cctpu.edu.ru/Math_method/math/12_files/image166.gif.

Заметим, что условия, приводимые в теоремах, являются достаточными, а отнюдь не необходимыми. Поэтому в практике численных расчетов нарушение этих условий обычно не вызывает заметного ухудшения расчетных схем.

Страница:  1  2  3 


Другие рефераты на тему «Математика»:

Поиск рефератов

Последние рефераты раздела

Copyright © 2010-2024 - www.refsru.com - рефераты, курсовые и дипломные работы