自控系统数字仿真实验代码(3)

X=X+h*(K1+2*K2+2*K3+K4)/6;

y=[y,C*X];

t=[t,t(i)+h];

end

k=1 ;

V=1;

a=[0.25 1 0 0];

b=[2*k k];

X0=[0 0 0 ];

T0=0;Tf=10;

n=3 ;

h=0.25;

R=1 ;

f23

plot(t,y)

2.4

A=[-21,19,-20;19,-21,20;40,-40,-40];

T=1 ; h=0.01 h=0.01 ;

x=[1;0;-1];

y= [0 ; 0 ; 0 ; ];

t=0 ;

N=T/h ;

for i=1:N

k1=A*x;

k2=A*(x+h*k1/2);

k3=A*(x+h*k2/2);

k4=A*(x+h*k3);

x=x+h*(k1+2*k2+2*k3+k4)/6 ;

y=[y,x];

t=[t,t(i)+h];

end

plot(t,y) h=0.04