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