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

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

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

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

Измерение спектральных характеристик волоконных световодов с органическими красителями
Измерение температуры является одной из важнейших и неотъемлемых составляющих многих технологических процессов. Однако в областях техники с воздействием сильных электромагнитных полей, например, в силовой энергетике [1-3] (силовые электрические машины, мощные трансформаторы, усилители и т.д.), в системах с СВЧ-излучением (мощная р ...

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