Скачиваний:
43
Добавлен:
10.05.2014
Размер:
111.27 Кб
Скачать

Министерство образования и науки Российской Федерации

Федеральное государственное бюджетное образовательное учреждение Высшего профессионального образования

НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ЯДЕРНЫЙ УНИВЕРСИТЕТ

Отчёт по лабораторной работе № 4

«ЧИСЛЕННОЕ РЕШЕНИЕ СМЕШАННОЙ КРАЕВОЙ ЗАДАЧИ ДЛЯ ВОЛНОВОГО УРАВНЕНИЯ»

Программа/

Результат

Теоретический материал

Оценка

Группа: Ф6-0

Курс: Численные методы

Выполнил:

Преподаватель: Рябов Павел Николаевич

Москва 2012 год

1. Полная постановка задачи:

2. График точного решения для различных моментов времени :

3. Результаты, полученные при выполнении лабораторной работы:

Узлы пространственной сетки,

Численное решение задачи в узлах сетки в момент времени

Точное решение задачи в узлах сетки,

Модуль разности между точным и приближённым решением в узлах сетки,

Аппроксимация 1 порядка

Аппроксимация 2 порядка

Аппроксимация 1 порядка

Аппроксимация 2 порядка

0.0

0.05

0.10

0.15

0.20

0.25

0.30

0.35

0.40

0.45

0.50

0.55

0.60

0.65

0.70

0.75

0.80

0.85

0.90

0.95

1.00

-1.13983

-1.03998

-0.93572

-0.82723

-0.71471

-0.59873

-0.47916

-0.35650

-0.23187

-0.10619

0.01988

0.14568

0.27062

0.39396

0.51465

0.63156

0.74502

0.85561

0.96218

1.06475

1.16305

-1.156245

-1.055745

-0.950904

-0.841939

-0.729267

-0.613178

-0.494183

-0.372864

-0.249750

-0.125404

-0.000433

0.124540

0.248892

0.372009

0.493296

0.612188

0.728199

0.840899

0.949858

1.054765

1.155339

-1.155293

-1.054748

-0.949872

-0.840939

-0.728282

-0.612297

-0.493438

-0.372213

-0.249170

-0.124896

0.000000

0.124896

0.249170

0.372213

0.493438

0.612297

0.728282

0.840939

0.949872

1.054748

1.155293

1.55E-02

1.48E-02

1.42E-02

1.37E-02

1.36E-02

1.36E-02

1.43E-02

1.57E-02

1.73E-02

1.87E-02

1.99E-02

2.08E-02

2.14E-02

2.18E-02

2.12E-02

1.93E-02

1.67E-02

1.47E-02

1.23E-02

1.00E-02

7.76E-03

9.52E-04

9.97E-04

1.03E-03

1.00E-03

9.86E-04

8.81E-04

7.45E-04

6.51E-04

5.80E-04

5.08E-04

4.33E-04

3.55E-04

2.78E-04

2.03E-04

1.43E-04

1.09E-04

8.22E-05

3.98E-05

1.43E-05

1.73E-05

4.66E-05

Максимальный модуль разности

при t = 0.5

2.18e-2

1.03e-3

Узлы пространственной сетки,

Численное решение задачи в узлах сетки в момент времени

Точное решение задачи в узлах сетки,

Модуль разности между точным и приближенным решением в узлах сетки,

Аппроксимация 1 порядка

Аппроксимация 2 порядка

Аппроксимация 1 порядка

Аппроксимация 2 порядка

0.0

0.05

0.10

0.15

0.20

0.25

0.30

0.35

0.40

0.45

0.50

0.55

0.60

0.65

0.70

0.75

0.80

0.85

0.90

0.95

1.00

-1.8691

-1.8131

-1.7533

-1.6895

-1.6216

-1.5494

-1.4730

-1.3925

-1.3085

-1.2204

-1.1272

-1.0290

-0.9266

-0.8192

-0.7071

-0.5916

-0.4730

-0.3520

-0.2291

-0.1048

0.0202

-1.90450

-1.85007

-1.79143

-1.72845

-1.66092

-1.58878

-1.51186

-1.43018

-1.34367

-1.25236

-1.15636

-1.05577

-0.95079

-0.84173

-0.72898

-0.61291

-0.49397

-0.37265

-0.24951

-0.12513

-0.00013

-1.90399

-1.84946

-1.79074

-1.72767

-1.66009

-1.58787

-1.51092

-1.42917

-1.34262

-1.25130

-1.15529

-1.05475

-0.94987

-0.84094

-0.72828

-0.61230

-0.49344

-0.37221

-0.24917

-0.12490

0.00000

3.49E-02

3.63E-02

3.74E-02

3.82E-02

3.85E-02

3.85E-02

3.80E-02

3.67E-02

3.42E-02

3.09E-02

2.81E-02

2.57E-02

2.33E-02

2.17E-02

2.11E-02

2.06E-02

2.04E-02

2.02E-02

2.01E-02

2.01E-02

2.02E-02

5.18E-04

6.10E-04

6.83E-04

7.72E-04

8.27E-04

9.08E-04

9.41E-04

1.00E-03

1.05E-03

1.06E-03

1.07E-03

1.02E-03

9.17E-04

7.96E-04

6.95E-04

6.18E-04

5.28E-04

4.38E-04

3.37E-04

2.34E-04

1.29E-04

Максимальный модуль разности

при t = 1.0

3.85e-2

1.07e-3

Узлы пространственной сетки,

Численное решение задачи в узлах сетки в момент времени

Точное решение задачи в узлах сетки,

Модуль разности между точным и приближенным решением в узлах сетки,

Аппроксимация 1 порядка

Аппроксимация 2 порядка

Аппроксимация 1 порядка

Аппроксимация 2 порядка

0.0

0.05

0.10

0.15

0.20

0.25

0.30

0.35

0.40

0.45

0.50

0.55

0.60

0.65

0.70

0.75

0.80

0.85

0.90

0.95

1.00

-2.2268

-2.2006

-2.1738

-2.1452

-2.1138

-2.0799

-2.0437

-2.0042

-1.9613

-1.9153

-1.8653

-1.8114

-1.7534

-1.6908

-1.6239

-1.5521

-1.4754

-1.3937

-1.3070

-1.2148

-1.1165

-2.26282

-2.23921

-2.21343

-2.18524

-2.15448

-2.12094

-2.08446

-2.04476

-2.00170

-1.95503

-1.90455

-1.85008

-1.79140

-1.72836

-1.66082

-1.58858

-1.51162

-1.42987

-1.34326

-1.25185

-1.15578

-2.26287

-2.23923

-2.21338

-2.18513

-2.15431

-2.12071

-2.08414

-2.04439

-2.00125

-1.95452

-1.90399

-1.84946

-1.79074

-1.72767

-1.66009

-1.58787

-1.51092

-1.42917

-1.34262

-1.25130

-1.15529

3.61E-02

3.86E-02

3.96E-02

4.00E-02

4.05E-02

4.08E-02

4.04E-02

4.01E-02

3.99E-02

3.92E-02

3.87E-02

3.81E-02

3.73E-02

3.68E-02

3.62E-02

3.57E-02

3.55E-02

3.55E-02

3.56E-02

3.65E-02

3.88E-02

5.28E-05

2.29E-05

4.86E-05

1.06E-04

1.76E-04

2.29E-04

3.20E-04

3.72E-04

4.48E-04

5.11E-04

5.61E-04

6.27E-04

6.55E-04

6.90E-04

7.25E-04

7.05E-04

7.03E-04

6.96E-04

6.39E-04

5.47E-04

4.87E-04

Максимальный модуль разности

при t = 1.5

4.08e-2

7.25e-4

Соседние файлы в папке 4. Волновое уравнение. Шаблон-крест