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

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

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

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

Отчёт по лабораторной работе № 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.98597

-1.97216

-1.95242

-1.92730

-1.89744

-1.86244

-1.82199

-1.77694

-1.72792

-1.67426

-1.61604

-1.55419

-1.48814

-1.41814

-1.34478

-1.26792

-1.18712

-1.10319

-1.01643

-0.92586

-0.83229

-1.98052

-1.96396

-1.94249

-1.91616

-1.88504

-1.84920

-1.80874

-1.76376

-1.71436

-1.66067

-1.60284

-1.54099

-1.47528

-1.40589

-1.33298

-1.25673

-1.17734

-1.09500

-1.00991

-0.92226

-0.83229

-1.97998

-1.96340

-1.94192

-1.91557

-1.88444

-1.84860

-1.80814

-1.76316

-1.71378

-1.66011

-1.60229

-1.54046

-1.47479

-1.40543

-1.33255

-1.25635

-1.17700

-1.09472

-1.00969

-0.92215

-0.83229

5.99E-03

8.76E-03

1.05E-02

1.17E-02

1.30E-02

1.38E-02

1.38E-02

1.38E-02

1.41E-02

1.42E-02

1.37E-02

1.37E-02

1.33E-02

1.27E-02

1.22E-02

1.16E-02

1.01E-02

8.48E-03

6.73E-03

3.71E-03

4.44E-16

5.36E-04

5.59E-04

5.72E-04

5.83E-04

5.95E-04

5.94E-04

5.95E-04

5.92E-04

5.80E-04

5.67E-04

5.52E-04

5.26E-04

4.95E-04

4.69E-04

4.31E-04

3.82E-04

3.33E-04

2.84E-04

2.16E-04

1.19E-04

4.44E-16

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

при t = 3.0

1.42e-2

5.95e-4

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

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

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

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

Аппроксимация 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.2721

-1.3460

-1.4169

-1.4835

-1.5463

-1.6054

-1.6600

-1.7106

-1.7582

-1.8024

-1.8418

-1.8765

-1.9067

-1.9335

-1.9563

-1.9737

-1.9857

-1.9924

-1.9936

-1.9891

-1.9800

-1.30685

-1.38088

-1.45145

-1.51839

-1.58153

-1.64071

-1.69578

-1.74660

-1.79307

-1.83507

-1.87250

-1.90525

-1.93323

-1.95639

-1.97466

-1.98801

-1.99638

-1.99977

-1.99816

-1.99156

-1.97998

-1.30729

-1.38130

-1.45186

-1.51880

-1.58194

-1.64112

-1.69620

-1.74704

-1.79352

-1.83551

-1.87291

-1.90564

-1.93360

-1.95672

-1.97496

-1.98826

-1.99659

-1.99993

-1.99827

-1.99162

-1.97998

3.52E-02

3.53E-02

3.50E-02

3.53E-02

3.56E-02

3.58E-02

3.62E-02

3.65E-02

3.53E-02

3.31E-02

3.11E-02

2.92E-02

2.69E-02

2.32E-02

1.87E-02

1.45E-02

1.09E-02

7.50E-03

4.69E-03

2.50E-03

0.00E+00

4.37E-04

4.27E-04

4.16E-04

4.06E-04

4.02E-04

4.06E-04

4.20E-04

4.39E-04

4.47E-04

4.40E-04

4.16E-04

3.88E-04

3.68E-04

3.38E-04

2.97E-04

2.53E-04

2.08E-04

1.63E-04

1.10E-04

5.26E-05

0.00E+00

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

при t = 4.0

3.65e-2

4.47e-4

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

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

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

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

Аппроксимация 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

0.602

0.508

0.413

0.316

0.217

0.118

0.019

-0.081

-0.181

-0.282

-0.384

-0.485

-0.586

-0.684

-0.781

-0.876

-0.967

-1.056

-1.143

-1.227

-1.307

0.5687

0.4721

0.3744

0.2757

0.1764

0.0766

-0.0235

-0.1235

-0.2232

-0.3223

-0.4206

-0.5179

-0.6139

-0.7083

-0.8010

-0.8917

-0.9801

-1.0661

-1.1494

-1.2299

-1.3073

0.5673

0.4708

0.3730

0.2744

0.1750

0.0752

-0.0248

-0.1247

-0.2243

-0.3234

-0.4216

-0.5188

-0.6147

-0.7090

-0.8016

-0.8922

-0.9805

-1.0664

-1.1496

-1.2300

-1.3073

3.50E-02

3.75E-02

3.98E-02

4.16E-02

4.24E-02

4.31E-02

4.40E-02

4.41E-02

4.29E-02

4.09E-02

3.77E-02

3.34E-02

2.89E-02

2.45E-02

2.02E-02

1.66E-02

1.39E-02

1.07E-02

6.83E-03

3.04E-03

0.00E+00

1.37E-03

1.38E-03

1.38E-03

1.38E-03

1.38E-03

1.35E-03

1.30E-03

1.24E-03

1.15E-03

1.06E-03

9.67E-04

8.81E-04

7.90E-04

7.03E-04

6.13E-04

5.15E-04

4.19E-04

3.18E-04

2.12E-04

1.08E-04

0.00E+00

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

при t = 5.0

4.41e-2

1.38e-3

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