《限定记忆最小二乘法》由会员分享，可在线阅读，更多相关《限定记忆最小二乘法（5页珍藏版）》请在金锄头文库上搜索。
1、%限定记忆最小二乘的递推算法 %Z(k+2)=1.5*Z(k+1)-0.7*Z(k)+u(k+1)+0.5*u(k)+v(k) %= clear clc%=4 0 0 产生 M 序列作为输入=x=0 1 0 1 1 0 1 1 1; %initial value n=403; %n 为脉冲数目M=; %存放M序列for i=1:ntemp=xor(x(4),x(9);M(i)=x(9);for j=9:-1:2x(j)=x(j-1);endx(1)=temp;end%=产生均值为 0，方差为1的高斯白噪声= v=randn(1,402);%=产生观测序列=z=zeros(402,1);z(1)
2、=-1;z(2)=0;for i=3:402 z(i)=1.5*z(i-1)-0.7*z(i-2)+M(i-1)+0.5*M(i-2)+v(i); end% 递推求解P_a=100*eye(4); %估计方差Theta_a=3;3;3;3;Pstore=zeros(4,381); % 存放方差中间过程Theta_Store=zeros(4,381); %参数的估计值，存放中间过程估值L=20; %记忆长度for i=3:L-1 h=-z(i-1);-z(i-2);M(i-1);M(i-2);K=P_a*h*inv(h*P_a*h+1); Theta_a=Theta_a+K*(z(i)-h*Th
3、eta_a);P_a=(eye(4)-K*h)*P_a;endfor k=0:380 hL=-z(k+L-1);-z(k+L-2);M(k+L-1);M(k+L-2); K_b=P_a*hL*inv(1+hL*P_a*hL);Theta_b=Theta_a+K_b*(z(k+L)-hL*Theta_a); P_b=(eye(4)-K_b*hL)*P_a;hk=-z(k+L);-z(k+L-1);M(k+L);M(k+L-1); K_a=P_b*hk*inv(1+hk*P_b*hk);Theta_a=Theta_b-K_a*(z(k+L+1)-hk*Theta_b); P_a=(eye(4)+K
4、_a*hk)*P_b;Pstore(:,k+1)=P_a(1,1),P_a(2,2),P_a(3,3),P_a(4,4); Theta_Store(:,k+1)=Theta_a;end %=输出结果及作图= disp(参数al a2 bl b2的估计值为：)Theta_ai=1:381;figure(1)plot(i,Theta_Store(1,:),i,Theta_Store(2,:),i,Theta_Store(3,:),i,Theta_S tore(4,:)title(待估参数过渡过程)figure(2)plot(i,Pstore(1,:),i,Pstore(2,:),i,Pstore(
5、3,:),i,Pstore(4,:)title( 估计方差过渡过程)运行结果如下：参数 a1 a2 b1 b2 的估计值为：Theta_a =-1.54260.72971.00690.3365彳 MATLAB 7.11.0 (R2010b)File Edit Debug Parallel Desktop Window Help:D tS | # 町 O| 掷国蓟 I I P:Program FilesMATLABR2010bbinShortcuts How to Add 创 Whats NewCurrent Folder? x卫；bin Q l|Name 勺 yiwangyinNi.m 巧 y
6、asuochuli2.m 巧 yasuochulil.m 勺 yasuochuli.m *1 xindianxinhaol.m *1 xindianxinhao.m 勺 xiaozaochuli.mxiaoguoduibi.m 巧 xiandingjiyi.m m worker.bat 巧 tuxiang.m 巧 tubianxinhao.m 勺 tubiandian.m _| tubiandian.asv *1 qiyidian.m _, qiyidian.asv D ProductRoots HI mw_mpiexec.bat 一 msvc_modules_installer.pm ,me
7、xutils.pm _| mexsetup.pm 鬲 mexext.bat .,mex.pl 1(強1 may DetailsWorkspace* 口？ x肓i到甩 11国叨“ ”JHddJHJHddJHddJdJHddValue0.0080;0.0( * -0.0068;0.00.0035;-0.0 三 201x403 doi 4x4 doubl 4x4 doubl 4x381 doi 4x381 doi -1.5426;0.7 Comma nd History-?image(IJ :”colormap (map): nbcol=size(map, 1)colorbar :；% 2011/6/30 19:08 %；% 2011/6/30 20:00 %i % 2011/7/1 9:29 % 2011/7/1 16:22 % 2011/7/2 9:46 % 2011/7/2 14:58 % 2011/7/2 19:47 %仿真图如下所示：图 1 待估参数过渡过程图 2 估计方差过渡过程
