系统辨识-最小二乘法MATLAB仿真.doc

《系统辨识-最小二乘法MATLAB仿真.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、《系统辨识》基于MATLAB的最小二乘法（一阶）的仿真clcclear%①白噪声的生成过程如下：e=randn(1,500);e=e/std(e);e=e-mean(e);A=0;%白噪声的均值为0B=sqrt(0.1);%白噪声的方差为0.1e=A+B*e;%绘制白噪声图k=1:500;subplot(4,1,1)%画四行一列图形窗口中的第一个图形plot(k,e,'r');xlabel('k'),ylabel('e');title('(0,1)均匀分布的随机序列')%②生成M序列的过程如下：X1=1;X2=0;X3=1;X4=0;%移位寄存器输入Xi初始状态（0101），Y

2、i寄存器的各级输出m=500;%M序列的总长度fori=1:mY4=X4;Y3=X3;Y2=X2;Y1=X1;X4=Y3;X3=Y2;X2=Y1;X1=xor(Y3,Y4);%异或运算ifY4==0U(i)=-1;elseU(i)=Y4;endendM=U;u=U;%绘制M序列图¼i1=ik=1:1:i1;subplot(4,1,2)%画四行一列图形窗口中的第二个图形plot(k,U,k,U,'rx')stem(M)xlabel('k')ylabel('M序列')title('移位寄存器产生的M序列')%③参数估计的过程如下：%绘制参数估计的相关图形z=zeros(1,500)

3、;%定义输出观测值的长度fork=2:500z(k)=0.9*z(k-1)+u(k-1)+e(k);%用理想输出值作为观测值endsubplot(4,1,3)%画四行一列图形窗口中的第三个图形i=1:1:500;%横坐标的范围从1到500，步长为1plot(i,z)%图形的横坐标是采样时刻i，纵坐标是输出观测值Z，图形格式为连续曲线subplot(4,1,4)%画四行一列图形窗口中的第四个图形stem(z),gridon%画出输出观测值Z的经线图形，并显示坐标网络u,z%显示输入信号和输出观测信号%给矩阵HL和ZL赋初值HL(:,1)=-z(1,1:499);HL(:,2)=u

4、(1,1:499);ZL(:,1)=z(1,2:500);%计算参数c1=HL'*HL;c2=inv(c1);c3=HL'*ZL;c=c2*c3a1=c(1),b1=c(2)%从中分离出a1¡b1仿真截图仿真结果之一：c=-9.018a1=b1=1.0108-9.0181.0108