Предыдущая страница реферата | 12  13  14  15  16  17  18  19  20  21  22 | Следующая страница реферата



procedure Right_Part;

begin

if s>smin then A2:=(R2t+(R2pt-R2t)*(S-Smin)/(1-Smin))/L21

else A2:=R2t/L21;

U1a:=Uss*cos(W1*t);

U1b:=Uss*sin(W1*t);

f[1]:=U1a-A1*x[1]+A1*K2*x[3];

f[2]:=U1b-A1*x[2]+A1*K2*x[4];

f[3]:=A2*(K1*x[1]-x[3])-x[5]*x[4];

f[4]:=A2*(K1*x[2]-x[4])+x[5]*x[3];

M:=3/2*p*L0/(L1*L11)*(x[2]*x[3]-x[1]*x[4]);

f[5]:=p/J*(M-Ms);

S:=(W1-x[5])/W1;

end;

{ Модифицированный метод Рунге-Кутта 4-го порядка }

procedure Runge;

begin

for k:=1 to 5 do h[k]:=x[k];

Right_Part;

for k:=1 to 5 do

begin

f1[k]:=f[k];

x[k]:=h[k]+f1[k]*dt/2;

end;

Right_Part;

for k:=1 to 5 do

begin

f2[k]:=f[k];

x[k]:=h[k]+f2[k]*dt/2;

end;

Right_Part;

for k:=1 to 5 do

begin

f3[k]:=f[k];

x[k]:=h[k]+f3[k]*dt;

end;

Right_Part;

for k:=1 to 5 do

begin

f4[k]:=f[k];

x[k]:=h[k]+(f1[k]+2*f2[k]+2*f3[k]+f4[k])/6*dt;

end;

end;

{ Инициализация графики }

procedure Init_Graph;

var GraphDriver,GraphMode:integer;

i:integer;

color:word;

Рекомендуем скачать другие рефераты по теме: мировая торговля, банк курсовых работ бесплатно.

Предыдущая страница реферата | 12  13  14  15  16  17  18  19  20  21  22 | Следующая страница реферата

Поделитесь этой записью или добавьте в закладки