最小二乘参数辨识方法与应用程序清单.docVIP

最小二乘参数辨识方法及应用程序清单 一、3.2.2节利用最小二乘法求取模型参数的程序 程序1： 热敏电阻例子 clear all close all clc T=[20.5 26 32.7 40 51 61 73 80 88 95.7]; R=[765 790 826 850 873 910 942 980 1010 1032]; [m,n]=size(T); figure plot(T,R,b+) t=0; z=0; tz=0; tt=0; for i=1:n t=t+T(i); tt=tt+T(i)*T(i); z=z+R(i); tz=tz+T(i)*R(i); end a=(tt*z-t*tz)/(n*tt-t*t) b=(n*tz-t*z)/(n*tt-t*t) %%%%%%%%%%%%%%%%%%%%%% A=polyfit(T,R,1); z=polyval(A,T); figure plot(T,z); figure plot(T,R,b+) hold on plot(T,z,r); hold off 程序2 t=[0.1:0.1:1]; N=10; y=a+b*t; v=rands(1,10); z=y+v; sum_t1=0;sum_t2=0;sum_tz=0; sum_N=0; sum_z=0; for i=1:1:N sum_t1=sum_t1+t(i); sum_t2=sum_t2+t(i)^2; sum_tz=sum_tz+t(i)*z(i); sum_N=sum_N+z(i); sum_z=sum_z+z(i); end ap_num=sum_t2*sum_N-sum_t1*sum_tz; ap_den=N*sum_t2-sum_t1^2; ap=ap_num/ap_den bp_num=N*sum_tz-sum_t1*sum_z; bp_den=N*sum_t2-sum_t1^2; bp=bp_num/bp_den 二、3.3 加权最小二乘算法 程序3：最小二乘参数辨识程序 clear all%清理工作间变量 close all%关闭所有图形 clc%清屏 z(1)=440,z(2)=430,z(3)=420,z(4)=380,z(5)=370,z(6)=360,z(7)=320,z(8)=310,z(9)=300,z(10)=260,z(11)=250,z(12)=240,z(13)=220,z(14)=210,z(15)=170,z(16)=160; u(1)=3,u(2)=2.7,u(3)=2.4,u(4)=2.1,u(5)=2.0,u(6)=1.9,u(7)=1.6,u(8)=1.54,u(9)=1.48,u(10)=1.2,u(11)=1.14,u(12)=1.08,u(13)=0.95,u(14)=0.9,u(15)=0.7,u(16)=0.6; HL=[-z(1) u(1);-z(2) u(2);-z(3) u(3);-z(4) u(4);-z(5) u(5);-z(6) u(6); -z(7) u(7); -z(8) u(8);-z(9) u(9); -z(10) u(10);-z(11) u(11);-z(12) u(12);-z(13) u(13);-z(14) u(14)] %给样本矩阵HL赋值 ZL=[z(3);z(4);z(5);z(6);z(7);z(8);z(9);z(10);z(11);z(12);z(13);z(14);z(15);z(16)]% 给样本矩阵zL赋值 %calculating parameters%计算参数 c1=HL*HL; c2=inv(c1); c3=HL*ZL; c=c2*c3 %计算并显示 %DISPLAY PARAMETERS a2=c(1), b2=c(2) 三、3.4 递推最小二乘算法 程序4：递推最小二乘参数辨识程序 clear all%清理工作间变量 close all%关闭所有图形 clc%清屏 %%%%%%%%%%%%%%%%%%%%%%%%%% 产生白噪声的程序%%%%%%%%%%%%%%%%%%%%%%%%%%%% A=6; x0=1;M=255; for k=1:10000 x2=A*x0; x1=mod (x2,M); v1=x1/256; v(:,k)=(v1-0.5)*2; x0=x1; v0=v1; end num=v; k1=k; % num=zeros(1,1000); % randn(seed,100);