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

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

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

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

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

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

при t = 3.0

.e-

.e-

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

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

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

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

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

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

при t = 4.0

.e-

.e-

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

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

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

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

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

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

при t = 5.0

.e-

.e-

4. Полулогарифмический график зависимости максимальной по­греш­ности от времени:

Код с 1 и 2 порядком аппроксимации начальных и граничных условий, шаблон-крест

uses crt;

Function f(x,t:real):real;

var tangent:real;

begin

tangent:=sin(x-t)/cos(x-t);

f:=-1+(1+2*tangent*tangent)/4/cos(x-t) ;

end;

Function phi1(x:real):real;

begin

phi1:=x*x+1/2/cos(x);

end;

Function dphi1(x:real):real;

var cosx:real;

begin

cosx:=cos(x);

dphi1:=2+(2-cosx*cosx)/2/cosx/cosx/cosx;

end;

Function phi2(x:real):real;

begin

phi2:=-sin(x)/2/cos(x)/cos(x);

end;

Function psi1(t:real):real;

begin

psi1:=(1+sin(t)/cos(t))/2/cos(t);

end;

Function psi2(t:real):real;

begin

psi2:=1+1/2/cos(1-t);

end;

Function U0(x,t:real):real;

begin

U0:=x*x+1/2/cos(x-t);

end;

var

j,k,nx,nt:integer;

hx,ht,xmax,tmax,aa,sr,m:real;

x,t:array[0..500] of real;

u1,u2: array [0..500,0..500] of real;

Output:text;

begin

clrscr;

Assign(output,'out4.txt');

Rewrite(output);

aa:=1/2;

hx:=0.05;

ht:=0.05;

xmax:=1;

Write('T max = '); Readln(tmax);

nx:=round(xmax/hx);

nt:=round(tmax/ht);

x[0]:=0;

t[0]:=0;

For j:=1 to nx do

x[j]:=x[j-1]+hx;

For k:=1 to nt do

t[k]:=t[k-1]+ht;

{1 poryadok approks.}

For j:=0 to nx do

begin

U1[j,0]:=phi1(x[j]);

U1[j,1]:=U1[j,0]+ht*phi2(x[j]);

end;

U1[nx,1]:=psi2(ht);

For k:=1 to nt-1 do

begin

For j:=1 to nx-1 do

U1[j,k+1]:=2*U1[j,k]-U1[j,k-1]+aa*ht*ht/(hx*hx)*(U1[j+1,k]-2*U1[j,k]+U1[j-1,k])+f(x[j],t[k])*ht*ht;

U1[0,k+1]:=(U1[1,k+1]+hx*psi1(t[k+1])) / (1+hx);

U1[nx,k+1]:=psi2(t[k+1]);

end;

{2 poryadok approks.}

For j:=0 to nx do

begin

U2[j,0]:=phi1(x[j]);

U2[j,1]:=phi1(x[j])+ht*phi2(x[j])+ht*ht/2*(aa*dphi1(x[j])+f(x[j],0));

end;

U2[nx,1]:=psi2(ht);

For k:=1 to nt-1 do

begin

For j:=1 to nx-1 do

U2[j,k+1]:=2*U2[j,k]-U2[j,k-1]+aa*ht*ht/(hx*hx)*(U2[j+1,k]-2*U2[j,k]+U2[j-1,k])+f(x[j],t[k])*ht*ht;

U2[0,k+1]:=2*aa*ht*ht/(hx*hx)*(U2[1,k]-(1+hx)*U2[0,k]+hx*psi1(t[k]))+2*U2[0,k]-U2[0,k-1]+ht*ht*f(0,t[k]);

U2[nx,k+1]:=psi2(t[k+1]);

end;

k:=nt;

Writeln(output,'T= ',t[k]);

Writeln(output,'Tochnoe reshenie');

For j:=0 to nx do

Writeln(output,U0(x[j],t[k]));

Writeln(output);

Writeln(output,'Shislennoe 1-go poryadka');

For j:=0 to nx do

Writeln(output,U1[j,k]);

Writeln(output);

Writeln(output,'Raznost 1-go poryadka s tochnim');

For j:=0 to nx do

Writeln(output,abs(u1[j,k]-U0(x[j],t[k])));

Writeln(output);

Writeln(output,'Shislennoe 2-go poryadka');

For j:=0 to nx do

Writeln(output,U2[j,k]);

Writeln(output);

Writeln(output,'Raznost 2-go poryadka s tochnim');

For j:=0 to nx do

Writeln(output,abs(u2[j,k]-U0(x[j],t[k])));

Writeln(output);

Writeln(output,'-------');

Writeln(output);

{1}

Writeln(output,'maximaln. modul raznosti dlya 1-go poryadka');

k:=nt;

m:=0;

For j:=0 to nx do

begin

sr:=abs(u1[j,k]-U0(x[j],t[k]));

If m<sr

then m:=sr;

end;

Writeln(output,m);

{2}

Writeln(output,'maximaln. modul raznosti dlya 2-go poryadka');

k:=nt;

m:=0;

For j:=0 to nx do

begin

sr:=abs(u2[j,k]-U0(x[j],t[k])) ;

If m<sr

then m:=sr;

end;

Writeln(output,m);

Close(output);

Readln;

end.

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