IV. Метод сеток

Пусть - сетка в некоторой областиG; - линейное пространство сеточных функций, заданных на;- линейное пространство гладких функций;- норма в;- норма в. Предполагается, что

во-первых, существует оператор проектирования такой, что

для любого,

во-вторых, нормы исогласованы, т.е.

.

Рассмотрим некоторый дифференциальный оператор L, заданный в , и оператор, преобразующий сеточную функциюв сеточную функцию, заданную на.

Погрешностью аппроксимации оператора L разностным оператором называется сеточная функция:

,

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

Если при, то говорят, что разностный оператораппроксимирует дифференциальный операторL. Если при этом

то разностный оператор аппроксимирует дифференциальный операторL с порядком

Порядок выполнения работы

Классификация краевой задачи.

Линейное дифференциальное уравнение второго порядка такого вида относится к уравнениям гиперболического типа. Данное уравнение описывает колебание струны в зависимости от времени. Здесь переменная у имеет физический смысл времени, поэтому в дальнейшем будет рассматриваться такая задача:

гдеи.

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

В данной задаче начальные и граничные условия имеют следующий смысл:

Начальные условия (в начальный момент времени t=0):

задаёт форму струны.

задаёт распределение скоростей.

Граничные условия:

закон движения левого конца струны.

- скорость на правом конце струны.

Натяжение струны и плотность материала струны равны 1.

- внешняя сила отсутствует.

Явная разностная схема

j+1 i i+1

j

j-1

Рис.1

Для сведения задачи к явной разностной схеме используем шаблон «крест», изображённый на рисунке 1.

Конечно-разностная схема для данной задачи имеет вид:

Обозначим и выразимчерез остальные значения сеточной функции, входящие в уравнение:

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

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

Вычтя из первого выражения второе получим формулу для начального условия второго рода с использованием мнимого узла.

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

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

Порядок аппроксимации данной разностной схемы равен .

Для нахождения решений используем уравнения.

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

Исследование на сходимость явной разностной схемы

При исследовании сходимости решения используем следующий метод: возьмём начальную сетку, вычислим значение Uв узлах этой сетки. Затем, увеличивая размерность матрицыU,при этом, оставляя первоначальные узлы, сравним максимальное расхождение двух матриц. Если эта величина будет монотонно убывать, то решение сходится и схема устойчива.

Результаты исследования сходимости матрицы решения следующие:

Размерность U : 100x100 ; 200x200;400x400.

Таблица норм (максимальных расхождений)

Размерность U	100х100	200х200	200х200	400х400	400х400	800х800
Норма (U2-U1)	0.0129		0.0064		0.0032

То есть решение сходится, если выполнено условие: .

Таблица решения

x\t	0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
0*	0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000
0,05*	0.1564 0.2562 0.3557 0.4547 0.5533 0.6516 0.7495 0.8470 0.9441 1.0409
0,1*	0.3090 0.4086 0.5075 0.6056 0.7029 0.7994 0.8952 0.9903 1.0847 1.1783
0,15*	0.4540 0.5534 0.6517 0.7489 0.8450 0.9399 1.0338 1.1265 1.2182 1.3089
0,2*	0.5878 0.6871 0.7849 0.8812 0.9761 1.0695 1.1616 1.2522 1.3415 1.4294
0,25*	0.7071 0.8062 0.9036 0.9992 1.0930 1.1852 1.2756 1.3643 1.4514 1.5368
0,3*	0.8090 0.9080 1.0050 1.0999 1.1929 1.2839 1.3729 1.4601 1.5453 1.6286
0,35*	0.8910 0.9899 1.0866 1.1810 1.2733 1.3634 1.4513 1.5371 1.6208 1.7025
0,4*	0.9511 1.0499 1.1463 1.2404 1.3321 1.4215 1.5087 1.5935 1.6761 1.7565
0,45*	0.9877 1.0865 1.1828 1.2766 1.3680 1.4570 1.5437 1.6279 1.7099 1.7895

Форма струны в разные моменты времени.

График, иллюстрирующий решение

Срезы решения