《系统建模与仿真》实验1的程序框图.docVIP

《系统建模与仿真》实验1 程序框图及编程思路 产生初始信号 ， 对于本实验而言，模型的阶次。 程序中变量的维数设计：数组，为数组的脚标。 ——标量 ——标量 ——列向量 代号sita_mao ——方阵 ——列向量 代号fai ——列向量 参考程序源代码： 辨识结果： 3 由按照原模型计算 置初值 置初值 计算 N=1 计算 计算 开始 计算 计算 计算 计算 N=N+1 ？ 结束 Y N a1=2; a2=1.3; b0=0.4; b1=0.88; b2=2.2; c=10000; for k=1:1:3 u(k)=1.5*sin(0.2*k); end y(1)=b0*u(1)+ksai(1); y(2)=-a1*y(1)+b0*u(2)+b1*u(1)+ksai(2); y(3)=-a1*y(2)-a2*y(1)+b0*u(3)+b1*u(2)+b2*u(1)+ksai(3); sita_mao(:,1)=[0,0,0,0,0]; z=c^2; P(:,:,1)=diag([z,z,z,z,z],0); for N=1:1:100 u(N+3)=1.5*sin(0.2*(N+3)); y(N+3)=-a1*y(N+2)-a2*y(N+1)+b0*u(N+3)+b1*u(N+2)+b2*u(N+1)+ksai(N+3); fai(:,N+1)=[-y(N+2),-y(N+1),u(N+3),u(N+2),u(N+1)]; W=P(:,:,N)*fai(:,N+1); r=1/(1+W*fai(:,N+1)); K(:,N+1)=W*r; P(:,:,N+1)=P(:,:,N)-K(:,N+1)*W; sita_mao(:,N+1)=sita_mao(:,N)+K(:,N+1)*(y(N+3)-(fai(:,N+1))*sita_mao(:,N)); end function y=ksai(x) y=0.2*rand(1)-0.1;