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

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

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

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

Модернизация канала формирования импульсов запуска блока Т-17М
Постоянное совершенствование средств воздушного нападения, опыт войн и вооруженных конфликтов последних десятилетий предъявляют повышенные требования к возможностям зенитных комплексов ближнего действия по обнаружению, сопровождению и поражению воздушных целей в различных условиях обстановки. Однако подавляющее большинство зенитны ...

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