Курсовая работа: Проектирование системы автоматического управления

K2=h*(F+K1/3);

K3=h*(F+K2/6+K1/6);

K4=h*(F+K1/8+3/8*K2);

K5=h*(F+K1/2-3/2*K3+2*K4);

Xt=Xt+(1./6)*(K1+4*K4+K5);

X(2,1)=Xt(1);

X(2,2)=Xt(2);

X(2,3)=Xt(3);

X(2,4)=Xt(4);

t(2)=t(1)+h;

F=A*Xt+B;

i=2;

%Неявный метод второго порядка

while t(i)<5

X1(1)=X(i-1,1);

X1(2)=X(i-1,2);

X1(3)=X(i-1,3);

X1(4)=X(i-1,4);

Xt=Xt+(h./12)*(5*B+8*(A*Xt+B)-(A*X1'+B));

Xt=((eye(4)-(5./12)*h*A)^-1)*Xt;

X(i+1,1)=Xt(1);

X(i+1,2)=Xt(2);

X(i+1,3)=Xt(3);

X(i+1,4)=Xt(4);

t(i+1)=t(i)+h;

i=i+1;

end

for j=1:i

V(j)=X(j,1);