Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
POS-KSC .DOC
Скачиваний:
60
Добавлен:
19.11.2018
Размер:
3.27 Mб
Скачать

7.4.1. Разностная схема уравнений

Разностную схему рассмотрим на примере уравнения Пуассона в прямоугольнике, используя для аппроксимации второй производной конечные разности второго порядка точности (7.1.4). Вводя сетку , получаем ,

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

(7.4.1)

Данная неявная схема охватывает все внутренние точки области , их количество . Таково же число уравнений и неизвестных в СЛАУ, построенной на основе (7.4.1).

Пусть для простоты для всех , т.е. для всех внутренних точек , , а граничные условия таковы: внизу , слева и справа , и только наверху задана отличная от нуля функция . Зададим , . Тогда , а число внутренних точек и уравнений . Матрица СЛАУ для данной задачи задается по следующему закону (на языке пакета Mathcad):

Здесь , - это индексы матрицы , они связаны с другими, ранее введенными индексами для узлов сетки. Сеточная функция , , выражается через найденный в результате решения СЛАУ вектор решения , следующим образом: .

Пятидиагональная матрица имеет следующее строение:

Вектор правой части СЛАУ задается в данной задаче по закону:

,

где - число, определяющее в индексах вектора решения начало последнего слоя внутренних узлов по оси , на которых учитывается заданное граничное условие. Заданная функция из уравнения, в данной задаче являющаяся константой, , входит в правую часть СЛАУ в виде слагаемого . Значения вектора правой части в данной задаче:

На рис.3 показано распределение функции решения аналогичной краевой задачи в двумерной области при порядке СЛАУ . Решение получено комбинированным методом Зейделя-ОСП при оптимальном параметре за итераций с относительной точностью решения в . Обычный метод Зейделя сходится здесь лишь за итераций и сопоставим по времени решения с прямым методом. Еще большее число требуемых итераций показывают в данной задаче метод ОСП с матрицей (1.2) - .

Отметим, что матрица задачи при и заданном постоянна и не зависит от краевых условий и источников , которые входят в правую часть СЛАУ. Соответственно задачи с различными краевыми условиями и источниками могут решаться с тем же самым оптимальным параметром, найденным один раз для данной сетки.

рис.3

    1. Лабораторные задания к теме «Численное решение уравнений в частных производных»

Лабораторные работы по теме могут быть выполнены с помощью математических пакетов программ Mathcad или Matlab. В результате работы должна быть представлена искомая сеточная функция в виде матрицы значений в узлах сетки либо в виде послойного по времени распределения значений сеточных векторов. Следует также привести графическое представление результатов.

7.5.1. Гиперболические уравнения

Варианты заданий для одномерного волнового уравнения с граничными и начальными условиями (см. 7.2).

№ п/п

Метод

1.1

1

1

0

0

2

0.1, 0.05

1

1.2

1

2

0

0

1

0.1, 0.05

1

1.3

1

0

0

2

0.1, 0.05

1

1.4

1

2

0

0

1

0.1, 0.05

1

1.5

1

0

0

2

0.1, 0.05

1

1.6

1

2

0

0

1

0.1, 0.05

1

1.7

1

0

0

2

0.1, 0.05

1

1.8

1

2

0

0

1

0.1, 0.05

1

1.9

1

0

0

2

0.1, 0.05

1

1.10

1

2

0

0

1

0.1, 0.05

1

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]