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

0 0 1;

-10 -19 -7.4];

%матрица B

B=[0 5 10]';

Y=[0 0 1]';

k=1;

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

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

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

while(k<=3)

%шаг

if(k==1) h=0.1; end;

if(k==2) h=1; end;

if(k==3) h=0.01; end;

%---------------------------%

n=1;

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

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

n=n+1;

while (n<=T/h)

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

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

n=n+1;

end

t=0:h:10;

%k=t/h+1;

i=1;

while(i<=n)

if(k==1) t1=t; x1(i)=X(i,1); Xa1=1-0.9202*exp(-0.6983*t)-0.4636*exp(-3.3508*t).*cos(1.7584*t+4.382)+0.2433*exp(-3.3508*t).*sin(1.7584*t+4.382); end;