
Добавил:
Studfiles2
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:
uses crt;
const
k1=960;
t1=0.006;
k2=3;
t2=0.03;
k3=2;
t3=0.2;
k4=0.001;
TM=4;
dt=0.1;
hp=1;
var a,b,time,g,eps,Iy,ud,amega,L:real;
i,IT:integer;
x:array [1..5] of real;
begin
clrscr;
L:=0;
IY:=0;
ud:=0;
amega:=0;
time:=0;
g:=1;
x[1]:=0;
x[2]:=0;
x[3]:=0;
x[4]:=0;
IT:=trunc(hp/dt);
repeat
for i:=1 to IT do begin
time:=time+dt;
eps:=g-L;
x[1]:=x[1]+(K1*eps-Iy)*dt;
Iy:=x[1]/t1;
x[2]:=x[2]+(k2*Iy-ud)*dt;
ud:=x[2]/t2;
x[3]:=x[3]+(k3*ud-amega)*dt;
amega:=x[3]/t3;
x[4]:=x[4]+(k4*amega-L)*dt;
L:=x[4]
end;
writeln('',eps:3:3,' ',l:3:3)
until time>TM;
readln
end.
const
k1=960;
t1=0.006;
k2=3;
t2=0.03;
k3=2;
t3=0.2;
k4=0.001;
TM=4;
dt=0.1;
hp=1;
var a,b,time,g,eps,Iy,ud,amega,L:real;
i,IT:integer;
x:array [1..5] of real;
begin
clrscr;
L:=0;
IY:=0;
ud:=0;
amega:=0;
time:=0;
g:=1;
x[1]:=0;
x[2]:=0;
x[3]:=0;
x[4]:=0;
IT:=trunc(hp/dt);
repeat
for i:=1 to IT do begin
time:=time+dt;
eps:=g-L;
x[1]:=x[1]+(K1*eps-Iy)*dt;
Iy:=x[1]/t1;
x[2]:=x[2]+(k2*Iy-ud)*dt;
ud:=x[2]/t2;
x[3]:=x[3]+(k3*ud-amega)*dt;
amega:=x[3]/t3;
x[4]:=x[4]+(k4*amega-L)*dt;
L:=x[4]
end;
writeln('',eps:3:3,' ',l:3:3)
until time>TM;
readln
end.
Соседние файлы в папке СПО