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

f31.m

A=diag(P(:, 1));B=diag(P(:, 2));

C=diag(P(:, 3));D=diag(P(:, 4));

m=length(WIJ(: ,1));

W0=zeros(n,1);W=zeros(n,n);

for k=1:m

if (WIJ(k,2)==0);W0(WIJ(k,1))=WIJ(k,3);

else W(WIJ(k,1),WIJ(k,2))=WIJ(k,3);

end;

end;

Q=B-D*W;Qn=inv(Q);

R=C*W-A; V1=C*W0;

Ab=Qn*R; b1=Qn*V1;

Y=Yt0;y=0;t=T0;

N=round((Tf-T0)/(h*L1));

for i=1:N;

for j=1:L1;

K1=Ab*Y+b1*y0;

K2=Ab*(Y+h*K1/2)+b1*y0;

K3=Ab*(Y+h*K2/2)+b1*y0;

K4=Ab*(Y+h*K3)+b1*y0;

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

end ;

y=[y;Y(nout)];

t=[t;t(i)+h*L1];

end

P=[0 0.07 1 0.14;

1 0.012 1 0;

0 0.05 1 0.15;

10 1 1 0;

1 0.01 0 0.0008];

WIJ=[1,0,1;1,4,-1;2,1,1;3,2,1;3,5,-1;4,3,1;5,4,1 ]; n=5;

y0=1;

Yt0=zeros(n,1);

h=0.01;

L1=5;

T0=0;

Tf=10;

nout=4;

f31

plot(t,y)

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库自控系统数字仿真实验代码(4)在线全文阅读。