Курсовая работа: Математическое моделирование в задачах расчета и проектирования систем автоматического управления

h=0.01;

%начальные условия

X(1,1:4)=[0 0 0 0];

%матрица А (X'=AX+BY)

A=[0 1 0 0;

0 0 1 0;

0 0 0 1;

-a0 -a1 -a2 -a3];

%матрица B

B=[b3 b2 b1 b0]';

Y=[0 0 0 1]';

n=1;

k=1;

while(k<=10)

F(n,1:4)=(A*(X(n,1:4))'+B.*Y)';

X(n+1,1:4)=(X(n,1:4)'+h/10*(F(n,1:4))')';% МетодЭйлера

n=n+1;

k=k+1;

end

X(2,1:4)=X(n,1:4);

n=2;

I=[1 0 0 0; 0 1 0 0; 0 0 1 0; 0 0 0 1];

while (n<=t/h)

F(n,1:4)=(A*(X(n,1:4))'+B.*Y)';

X(n+1,1:4)=(((I-5*h/12*A)^-1)*(X(n,1:4)'+h/12*(5*B.*Y+8*(F(n,1:4))'-(F(n-1,1:4))')))';

n=n+1;

end

i=1;

while(i<=n)