1
Цель работы: исследование вычислительных свойств метода ускоренной верхней релаксации при решении уравнения Пуассона на плоскости
1. Конечно-разностная аппроксимация задачи Дирихле для уравнения Пуассона на плоскости.
2. Метод ускоренной верхней релаксации для решения системы линейных уравнений в методе конечных разностей для уравнения Пуассона
№ |
U(0,y) |
U(1,y) |
U(x,0) |
U(x,1) |
f(x,y) |
Точное решение |
1 |
0 |
y2+y |
0 |
x2+x |
2(x+y) |
xy(y+x) |
ЗАДАЧА 1. Для заданного варианта тестовой задачи Дирихле для уравнения Пуассона Uxx+Uyy=f(x,y) в прямоугольной области на плоскости 0≤x≤1, 0≤y≤1 выполнить следующие действия:
1) На основании программного блока предыдущей лабораторной работы создать новый программный блок, реализовав в нем расчет по методу ускоренной верхней релаксации с использованием ускоряющего множителя ω (добавив его в список параметров).
2) Составить отдельный программный блок, в котором реализовать вызов программного блока для уравнения Пуассона внутри цикла по значениям параметра ω=1.00, 1.01 .. 1.99 с шагом 0.01 при h=0.05, ε=10-5. Результатом программного блока являются вектор значений {ωi}и вектор количества итераций для получения решения {iti}.
3) Построить график зависимости числа итераций от значения ω.
4) На основании полученных результатов определить наилучшее значение ускоряющего множителя при решении задачи, при котором количество итераций минимально.
ПРИЛОЖЕНИЕ. ВАРИАНТЫ ТЕСТОВОЙ ЗАДАЧИ ДИРИХЛЕ
2
2 |
y2 |
(y+1)2 |
x2 |
(x+1)2 |
4 |
(x+y)2 |
3 |
sin(πy) |
- sin(πy) |
0 |
0 |
- 2π2cos(πx)sin(π y) |
cos(πx)sin(πy) |
4 |
0 |
6y-4y2-2 |
x2-3x |
0 |
4y2-6y+2+4x212x |
(x2-3x)(2y2-3y+1) |
5 |
y2y/2+5/16 |
y2y/2+5/16 |
x2x+5/16 |
x2x+13/16 |
4 |
(x-1/2)2+(y-1/4)2 |
6 |
y |
9y2-5y+1 |
x |
4x2+1 |
18x+8y |
y(2x-1)2+x(3y-1)2 |
7 |
-y2 |
2y2+5y-1 |
-x |
7x-1 |
6x-2 |
x(5y-1)+y2(3x-1) |
8 |
-y2 |
1+y-2y2 |
x2 |
2x2-x-1 |
2(y-x) |
x2(y+1)-y2(x+1) |
9 |
-y2 |
2y2+4y-1 |
-x2 |
3x2+3x-1 |
6x+8y-4 |
y2(3x-1)+x2(4y-1) |
10 |
-1-y |
2y2-1 |
x2-x-1 |
2x2+x-2 |
2+2y+4x |
(x2- 1)(y+1)+x(2y2-1) |
11 |
-y |
5y2+2y-1 |
-x |
3x2+4x-1 |
6y+10x |
y(3x2-1)+x(5y2-1) |
12 |
y2 |
3y2-2 |
-x-x2 |
x2-x+1 |
4(x2+y2) |
y2(2x2-3x+1)+ +x(3y2-1)-x2 |
13 |
y2 |
(y+1)2 |
x2 |
(x+1)2 |
4 |
(x+y)2 |
14 |
0 |
y2+y |
0 |
x2+x |
2(x+y) |
xy(y+x) |
15 |
y2 |
(y+1)2 |
x2 |
(x+1)2 |
4 |
(x+y)2 |
16 |
Sin(πy) |
- sin(πy) |
0 |
0 |
- 2π2cos(πx)sin(π y) |
cos(πx)sin(πy) |
17 |
0 |
6y-4y2-2 |
x2-3x |
0 |
4y2-6y+2+4x212x |
(x2-3x)(2y2-3y+1) |
18 |
y2y/2+5/16 |
y2y/2+5/16 |
x2x+5/16 |
x2x+13/16 |
4 |
(x-1/2)2+(y-1/4)2 |
19 |
y |
9y2-5y+1 |
x |
4x2+1 |
18x+8y |
y(2x-1)2+x(3y-1)2 |
20 |
-y2 |
2y2+5y-1 |
-x |
7x-1 |
6x-2 |
x(5y-1)+y2(3x-1) |
21 |
-y2 |
1+y-2y2 |
x2 |
2x2-x-1 |
2(y-x) |
x2(y+1)-y2(x+1) |
22 |
-y2 |
2y2+4y-1 |
-x2 |
3x2+3x-1 |
6x+8y-4 |
y2(3x-1)+x2(4y-1) |
23 |
-1-y |
2y2-1 |
x2-x-1 |
2x2+x-2 |
2+2y+4x |
(x2- 1)(y+1)+x(2y2-1) |
3
24 |
-y |
5y2+2y-1 |
-x |
3x2+4x-1 |
6y+10x |
y(3x2-1)+x(5y2-1) |
25 |
y2 |
3y2-2 |
-x-x2 |
x2-x+1 |
4(x2+y2) |
y2(2x2-3x+1)+ +x(3y2-1)-x2 |
26 |
y2 |
(y+1)2 |
X2 |
(x+1)2 |
4 |
(x+y)2 |
Составил: 2010-2014 v1.3
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.