Файл: 17.06.вкр.docx

Она аппроксимирует первую производную в точке

В разложении Тейлора функции можно заменитьна –и получить левую разностную производную

.

Вычитаяизполучаем центральную разностную производную

[]

Если в ряде Тейлора оставить третий член ряда, то можно получить центральную разностную производную для аппроксимации

[]

Если исходить из разложения Тейлора функции двух переменных

+…

+…

то можно получить следующие аппроксимации частных производных

[]

[]

Разностные операторы, соответствующие дифференциальному уравнению, записываются во внутренних узлах сетки. Разностные операторы, соответствующие граничным условиям, записываются в граничных узлах сетки. В результате получаем систему алгебраических уравнений, число которых пропорционально числу внутренних узлов сеточной области. Для получения численного решения требуется решить эту систему уравнений.

Удобно, если предполагается использование ЭВМ для реализации вычислений, перейти к следующим обозначениям:

, =,,

=,,

[],(x,y) =[],

= [],(x,y) =[].

Пример 11.1. Рассмотрим задачу теплопроводности в стержне, начальная температура которого равна нулю. Пусть температура левого конца фиксирована, а на правом конце происходит теплообмен с окружающей средой, так что тепловой поток пропорционален разности температур конца стержня и среды. Пусть температура среды определяется g(t). Другими словами, решим задачу:

Решение.

Построим в плоскости прямоугольную сетку , узлыкоторой определяются формулами:,,

Значения на левой и нижней сторонах сетки известны их граничных и начальных условий. Наша задача состоит в отыскании остальных значений.

Для решения задачи заменим частные производные в уравнении теплопроводности их конечно-разностными аппроксимациями

[].

(x, y)[]

Подставим эти выражения в наше уравнение и разрешим получившееся уравнение относительно значений функций на верхнем временном слое. Имеем:

[] = =[],

Отсюда

= = +

Получившаяся формула выражает решение в данный момент времени через решение в предыдущий момент времени (индекс i относится к временной переменной).

Аппроксимируем производную в граничном условиина правом конце, заменивлевой разностной производной, поскольку правая разностная производная требует значений функции за пределами сетки:

,

Отсюда находим:

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

Шаг 1. Находим решение на сеточном слое t= , используя явную формулу

Шаг 2. Величину находим по формуле

Выполнив шаги, получаем решение для. Повторив эти шаги, получаем решение прии т.д.

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

ПримерРешим задачу:

Решение.

Положим– шаг изменения пространственной переменной. Заменим частные производные в волновом уравнении конечно-разностными аппроксимациями

(x, y) = [],(x, y) =[].

Получим:

[] =].

Отсюда

+

Получили явную разностную схему, которая будет устойчивой если . Отсюда k. Выберем

Построим алгоритм решения задачи:

Шаг 1.Вводим сетку:.Создаем нулевой массив значенийразмера m.

Шаг 2.Задаем значения

Шаг 3.Заполняем первую и вторую строки массиваграничными условиями,(нулевой начальной скорости соответствует совпадение значений (смещений) в первом и втором столбцах)

Шаг 4.Заполняем первый и последний столбец массива граничными условиями(на концах струны смещение равно нулю в любой момент времени).

Шаг 5. Находим решение , используя разностную схему

+