资源描述：
《最小二乘参数辨识方法及应用程序清单》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第3章最小二乘参数辨识方法及应用程序清单一、3.2.2节利用最小二乘法求取模型参数的程序程序1：热敏电阻例子clearallcloseallclcT=[20.52632.740516173808895.7];R=[76579082685087391094298010101032];[m,n]=size(T);figureplot(T,R,'b+')t=0;z=0;tz=0;tt=0;fori=1:nt=t+T(i);tt=tt+T(i)*T(i);z=z+R(i);tz=tz+T(i)*R(i);enda=(tt*z-t*tz)/(n*tt-t*t)b=(n*tz-t*z)/
2、(n*tt-t*t)%%%%%%%%%%%%%%%%%%%%%%A=polyfit(T,R,1);z=polyval(A,T);figureplot(T,z);figureplot(T,R,'b+')holdonplot(T,z,'r');holdoff程序2t=[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;fori=1:1:Nsum_t1=sum_t1+t(i);sum_t2=sum_t2+t(i)^2;sum_tz=sum_tz+t(i)*z(
3、i);sum_N=sum_N+z(i);sum_z=sum_z+z(i);endap_num=sum_t2*sum_N-sum_t1*sum_tz;ap_den=N*sum_t2-sum_t1^2;ap=ap_num/ap_denbp_num=N*sum_tz-sum_t1*sum_z;bp_den=N*sum_t2-sum_t1^2;bp=bp_num/bp_den二、3.3加权最小二乘算法程序3：最小二乘参数辨识程序clearall%清理工作间变量closeall%关闭所有图形clc%清屏z(1)=440,z(2)=430,z(3)=420,z(4)=380,z(5)=3
4、70,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
5、)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赋值%calculatingparameters%计算参数c1=HL'*HL;c2=inv(c1);c3=HL'*ZL;c=c2*c
6、3%计算并显示%DISPLAYPARAMETERSa2=c(1),b2=c(2)三、3.4递推最小二乘算法程序4：递推最小二乘参数辨识程序clearall%清理工作间变量closeall%关闭所有图形clc%清屏%%%%%%%%%%%%%%%%%%%%%%%%%%产生白噪声的程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%A=6;x0=1;M=255;fork=1:10000x2=A*x0;x1=mod(x2,M);v1=x1/256;v(:,k)=(v1-0.5)*2;x0=x1;v0=v1;endnum=v;k1=k;%num=zeros(1,1000);%r
7、andn('seed',100);%num=randn(1,1000);%白噪声程序结束%%%%%%%%%%%%%M序列产生程序%%%L=15;%M序列的周期y1=1;y2=1;y3=1;y4=0;%四个移位积存器的输出初始值fori=1:L;%开始循环，长度为Lx1=xor(y3,y4);%第一个移位积存器的输入是第3个与第4个移位积存器的输出的“或”x2=y1;%第二个移位积存器的输入是第3个移位积存器的输出x3=y2;%第三个移位积存器的输入是第2个移位积存器的输出x4=y3;%第四个移位积存器的
