Разделы сайта

Текст программы

program kurs2;crt,dos;Kv=7.2; b0=7.2; a0=7.2; a1=1; a2=0.275; a3=0.017107; a4=0.000761; a5=0.0000211; h=0.05;mas=array [1 5,0 100] of real;,d,c0,c1,c2,c3,c4:real;,n,j:integer;,K1,K2,K3,K4:mas;:array [1 3] of real;:char;;:=0;:=b0/a5;:=a0/a5;:=a1/a5;:=a2/a5;:=a3/a5;:=a4/a5;j:=1 to 5 do[j,0]:=0;j:=1 to 5 do begin K1[j,0]:=0; K2[j,0]:=0; K3[j,0]:=0; K4[j,0]:=0; end; begin n:=0; repeat for j:=1 to 4 do K1[j,n]:=h*Z[j+1,n]; K1[5,n]:=h*(d-c0*Z[1,n]-c1*Z[2,n]-c2*Z[3,n]-c3*Z[4,n]-c4*Z[5,n]); j:=1 to 4 do K2[j,n]:=h*(Z[j+1,n]+K1[j+1,n]/2); K2[5,n]:=h*(d-c0*z[1,n]-c1*z[2,n]-c2*Z[3,n]-c3*Z[4,n]-c4*Z[5,n]); for j:=1 to 4 do K3[j,n]:=h*(Z[j+1,n]+K2[j+1,n]/2); K3[5,n]:=h*(d-c0*Z[1,n]-c1*Z[2,n]-c2*Z[3,n]-c3*Z[4,n]-c4*Z[5,n]); for j:=1 to 4 do K4[j,n]:=h*(Z[j+1,n]+K3[j+1,n]); K4[5,n]:=h*(d-c0*Z[1,n]-c1*Z[2,n]-c2*Z[3,n]-c3*Z[4,n]-c4*Z[5,n]); for j:=1 to 5 do begin Z[j,n+1]:=Z[j,n]+(K1[j,n]+2*K2[j,n]+2*K3[j,n]+K4[j,n])/6; end; y[k]:=sqr(z[2,n]); if k=3 then begin integral:=integral+(h/3)*(y[1]+4*y[2]+y[3]); y[1]:=y[3]; k:=1; end; textcolor(3); writeln('No итерациии - ',n,' Z1 = ', Z[1,n]:3:8,' T= ',((n*h)-0.05):3:2,' c.', ' Интеграл = ',integral:3:8); normvideo; c:=readkey; k:=k+1; n:=n+1; until c='e' end;

end.

Описание переменных

В программе использованы такие переменные:

Массив Z - решения диф. уравнений.

Массив К1, К2, К3, К4 - коэффиенты для решения диф. уравнений по методу Рунге-Кутта.

С0,с1,с2, с3,d - откорректированные коэффициенты для решения диф. уравнения 5-го порядка.

У - массив для нахождения интегральной оценки., n, k - счетчики столбцов и строк массивов.- переменная, в которую на каждом шаге заносится значение интегральной оценки.

Описание работы программы

На первом шаге программа обнуляет все массивы, вычисляет коэффициенты. Далее программа в цикле производит решение уравнений по методу Рунге-Кутта, при этом коэффициенты К на каждом шаге находятся сразу для всех уравнений.

Далее программа находит интегральную оценку и выводит результаты на экран. Итеррации будут продолжаться до тех пор пока мы не нажмём клавишу ‘e’.

Шаг на каждой итеррации увеличивается на 0,05. Интеграл считается на каждом 3-м шаке.

6.3 Результаты работы программы

Таблица 8.1

T

H(t), (kv/2)

интеграл

H(t), (kv/4)

интеграл

0.00

0,315

0,00

0,00

0.00000

0.15

0,973

0,240

0,158

0,0601

0.3

1,45

3,818

0,506

1,042

0.45

1,52

5,116

0,85

1,608

0.6

1,26

5,522

1,09

2,241

0.75

0,925

5,755

1,21

2,327

0.9

0,734

6,735

1,22

2,344

1.05

0,757

6,969

1,16

2,356

1.2

0,917

7,027

1,08

2,410

1.35

1,09

7,127

1,01

2,434

1.5

1,16

7,375

0,967

2,456

1.65

1,11

7,410

0,95

2,458

1.8

1,02

7,427

0,955

2,460

1.95

0,943

7,462

0,97

2,462

2.1

0,925

7,518

0,988

2,463

2.25

0,957

7,523

1

2,464

2.4

1,00

7,530

1,01

2,464

2.55

1,04

7,540

1,01

2,464

2.7

1,04

7,551

1,01

2,464

2.85

1,02

7,552

1

2,464

3.00

0,992

7,555

1

2,464

3,15

0,98

7,558

0,999

2,464

3,3

0,984

7,560

0,998

2,464

3,45

0,996

7,560

0,998

2,464

3,15

1,01

7,561

0,999

2,464

3,75

1,01

7,562

1

2,464

3,9

1,01

7,562

1

2,464

4,05

1,00

7,563

1

2,464

4,2

0,995

7,563

1

2,464

4,35

1,00

7,563

1

2,464

Перейти на страницу: 1 2

Самое читаемое:

Кодовый замок
Кодовые замки являются эффективным средством предотвращения доступа посторонних лиц к охраняемым помещениям. К их достоинствам можно отнести простоту в обращении, надёжность, возможность обеспечить высокую степень защиты, относительную лёгкость смены кода (по сравнению со сменой обычного механического замка). Также немаловажным ...

www.techstages.ru : Все права защищены! 2022