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

Страницы работы

Содержание работы

                                                                                                                                              1

ЛАБОРАТОРНАЯ РАБОТА ПММ-2-2-13 ИССЛЕДОВАНИЕ МЕТОДА УСКОРЕННОЙ ВЕРХНЕЙ  РЕЛАКСАЦИИ ПРИ РЕШЕНИИ УРАВНЕНИЯ ПУАССОНА НА ПЛОСКОСТИ

Цель работы: исследование вычислительных свойств метода ускоренной верхней релаксации при решении уравнения Пуассона на плоскости

ВОПРОСЫ ДЛЯ ИЗУЧЕНИЯ ПО КОНСПЕКТУ ЛЕКЦИЙ

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

-

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

-

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

Похожие материалы

Информация о работе