實(shí)驗(yàn)四離散相似法數(shù)字仿真

1、精選優(yōu)質(zhì)文檔-傾情為你奉上精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)專心-專注-專業(yè)精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)實(shí)驗(yàn)四離散相似法數(shù)字仿真【實(shí)驗(yàn)?zāi)康摹空莆针x散相似法仿真方法【實(shí)驗(yàn)內(nèi)容】用離散相似法仿真程序（參考課本sp4-4m）重現(xiàn)下圖輸出Y1的數(shù)據(jù)和曲線，并與四階龍格一庫塔法比較精度?！緦?shí)驗(yàn)的程序】P=1 0.01 1 0;0 0.085 1 0.17;1 0.01 1 0;0 0.051 1 0.15; 1 0.0067 70 0;1 0.15 0.21 0;0 1 130 0;1 0.01 0.1 0;1 0.01 0.0044 0;WIJ=1 0 1;2 1 1;2 9 -

2、1;3 2 1;4 3 1;4 8 -1;5 4 1;6 5 1;6 7 -0.212;7 6 1;8 6 1;9 7 1;n=9;Y0=1;Yt0=0 0 0 0 0 0 0 0 0;h=0.001;T=0;T0=0;Tf=5;Nout=7;A=P(:,1);B=P(:,2);C=P(:,3);D=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;for i=

3、1:n if(A(i)=0);%求積分環(huán)節(jié)或積分比例環(huán)節(jié)各系數(shù) FI(i)=1; FIM(i)=h*C(i)/B(i); FIJ(i)=h*h*C(i)/B(i)/2; FIC(i)=1;FID(i)=0; %積分環(huán)節(jié)系數(shù)求取完畢 if(D(i)=0);%若為積分比例，修正d FID(i)=D(i)/B(i); else %求慣性環(huán)節(jié)或慣性環(huán)節(jié)比例系數(shù) end; else FI(i)=exp(-h*A(i)/B(i); FIM(i)=(1-FI(i)*C(i)/A(i); FIJ(i)=h*C(i)/A(i)-FIM(i)*B(i)/A(i); FIC(i)=1;FID(i)=0; %慣性環(huán)節(jié)

4、或慣性環(huán)節(jié)比例系數(shù)求取完畢 if(D(i)=0); %若為慣性比例環(huán)節(jié)修正 c d FIM(i)=(1-FI(i)*D(i)/A(i); FIJ(i)=h*D(i)/A(i)-FIM(i)*B(i)/A(i); FIC(i)=C(i)/D(i)-A(i)/B(i); FID(i)=D(i)/B(i); else end; end;end;%求輸入作用Y=zeros(n,1);X=Y;y=0;Uk=zeros(n,1);Ub=Uk;%置初值t=T0:h:Tf;N=length(t);%建立時(shí)間序列向量 n維for k=1:N-1 Ub=Uk;%保存當(dāng)前一次輸入值 Uk=W*Y+W0*Y0;%保存

5、當(dāng)前uk Udot=(Uk-Ub)/h;%求當(dāng)前uk導(dǎo)數(shù) Uf=2*Uk-Ub;%求下一步輸入uf X=FI.*X+FIM.*Uk+FIJ.*Udot;%點(diǎn)運(yùn)算求取狀態(tài)和輸出 Y=FIC.*X+FID.*Uf; y=y,Y(Nout);end;figure(1)plot(t,y,g)hold以下運(yùn)用龍格-庫塔法實(shí)現(xiàn)：Q=B-D*w;%求出各陣Qn=inv(Q);R=C*w-A;V1=C*w0;Ab=Qn*R;b1=Qn*V1;Y=yt0;y=Y(nout);t=T0;N=round(Tf-T0)/h);%龍哥庫塔法for i=1:N 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; y=y,Y(nout); T=T,t+h; t=t+h;end;figu