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