Dalrate.ru

Построение рабочей локальной сети

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

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 of real;,d,c0,c1,c2,c3,c4:real;,n,j:integer;,K1,K2,K3,K4:mas;:array of real;:char;;:=0;:=b0/a5;:=a0/a5;:=a1/a5;:=a2/a5;:=a3/a5;:=a4/a5;j:=1 to 5 do:=0;j:=1 to 5 do begin K1:=0; K2:=0; K3:=0; K4:=0; end; begin n:=0; repeat for j:=1 to 4 do K1:=h*Z; K1:=h*(d-c0*Z-c1*Z-c2*Z-c3*Z-c4*Z); j:=1 to 4 do K2:=h*(Z+K1/2); K2:=h*(d-c0*z-c1*z-c2*Z-c3*Z-c4*Z); for j:=1 to 4 do K3:=h*(Z+K2/2); K3:=h*(d-c0*Z-c1*Z-c2*Z-c3*Z-c4*Z); for j:=1 to 4 do K4:=h*(Z+K3); K4:=h*(d-c0*Z-c1*Z-c2*Z-c3*Z-c4*Z); for j:=1 to 5 do begin Z:=Z+(K1+2*K2+2*K3+K4)/6; end; y:=sqr(z); if k=3 then begin integral:=integral+(h/3)*(y+4*y+y); y:=y; k:=1; end; textcolor(3); writeln('No итерациии - ',n,' Z1 = ', Z: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

Популярное:

Испытания схем увязки САУТ-ЦМ и МПЦ в лаборатории Главной задачей решаемой на железнодорожном транспорте всеми его службами, хозяйствами, техническими и организационными мероприятиями является обеспечение непрерывного, нормального (штатного) протекания основного технологического процесса (ТП) движения поездов. Движение поездов - это комплексный технологический процесс, состоящий из ...