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

Исключая из этих двух уравнений http://ad.cctpu.edu.ru/Math_method/math/12_files/image066.gif, найдем

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

Выразим теперь отсюда d="Рисунок 265" src="images/referats/3126/image035.png" alt="http://ad.cctpu.edu.ru/Math_method/math/12_files/image070.gif">:

http://ad.cctpu.edu.ru/Math_method/math/12_files/image072.gif(2.35)

Но, согласно формуле (2.34),

http://ad.cctpu.edu.ru/Math_method/math/12_files/image074.gif (2.36)

Сравнивая теперь (2.35) и (2.36), найдем, что

http://ad.cctpu.edu.ru/Math_method/math/12_files/image076.gif (2.37)

Пусть теперь i >0, то есть i=1, 2, ., n2. Выражая http://ad.cctpu.edu.ru/Math_method/math/12_files/image078.gif по формуле (2.34), получим:

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

Подставляя это в формулу (2.33), будем иметь

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

Разрешая полученное уравнение относительноhttp://ad.cctpu.edu.ru/Math_method/math/12_files/image084.gif, находим

http://ad.cctpu.edu.ru/Math_method/math/12_files/image086.gif, или

http://ad.cctpu.edu.ru/Math_method/math/12_files/image088.gif. (2.38)

Отсюда, сравнивая формулы (2.34) и (2.38), получаем для коэффициентов http://ad.cctpu.edu.ru/Math_method/math/12_files/image060.gifи http://ad.cctpu.edu.ru/Math_method/math/12_files/image062.gifрекуррентные формулы:

http://ad.cctpu.edu.ru/Math_method/math/12_files/image092.gif (2.39)

Так как http://ad.cctpu.edu.ru/Math_method/math/12_files/image094.gif и http://ad.cctpu.edu.ru/Math_method/math/12_files/image096.gif уже определены по формулам (2.37), то, используя формулы (2.39), можно последовательно определить коэффициенты http://ad.cctpu.edu.ru/Math_method/math/12_files/image060.gif и http://ad.cctpu.edu.ru/Math_method/math/12_files/image062.gif до http://ad.cctpu.edu.ru/Math_method/math/12_files/image099.gif и http://ad.cctpu.edu.ru/Math_method/math/12_files/image101.gif включительно. Эти вычисления называются прямым ходом метода прогонки.

Из формулы (2.33) при i=n2 и второго краевого условия (2.32) получаем

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

Разрешая эту систему относительноhttp://ad.cctpu.edu.ru/Math_method/math/12_files/image105.gif, будем иметь

http://ad.cctpu.edu.ru/Math_method/math/12_files/image107.gif. (2.40)

Теперь, используя (2.34) и первое краевое условие (2.32), мы можем последовательно найти http://ad.cctpu.edu.ru/Math_method/math/12_files/image109.gif. Это − обратный ход метода прогонки.

Итак, получаем следующую цепочку:

http://ad.cctpu.edu.ru/Math_method/math/12_files/image111.gif (2.41)

Для простейших краевых условий http://ad.cctpu.edu.ru/Math_method/math/12_files/image113.gif

формулы для http://ad.cctpu.edu.ru/Math_method/math/12_files/image115.gifи http://ad.cctpu.edu.ru/Math_method/math/12_files/image105.gif упрощаются. Полагая в этом случае http://ad.cctpu.edu.ru/Math_method/math/12_files/image117.gifиз формул (2.37), (2.40), (2.41) будем иметь

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

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

1) Существует ли решение алгебраической системы типа (2.31)?

2) Как фактически находить это решение?

3) Сходится ли разностное решение к точному при стремлении шага сетки к 0?

Можно доказать, что если краевая задача имеет вид

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

причем р(x)>0, то решение системы (2.31), (2.32) существует и единственно. Фактическое отыскание решения можно провести, например, методом прогонки. На третий вопрос дает ответ следующая

Теорема

Если http://ad.cctpu.edu.ru/Math_method/math/12_files/image008.gif и http://ad.cctpu.edu.ru/Math_method/math/12_files/image012.gif дважды непрерывно дифференцируемы, то разностное решение, соответствующее схеме с заменой

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


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

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

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

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