程序及注释.doc

第3章程序及注释 例3.3 考虑仿真对象 (3.41) 其中，是服从正态分布的白噪声N。输入信号采用4阶M序列，幅度为1。选择如下形式的辨识模型 (3.42) 设输入信号的取值是从k =1到k =16的M序列，则待辨识参数为=。其中，被辨识参数、观测矩阵z L、H L的表达式为 , ， (3.43) 程序框图如图3.2所示。Matlab6.0仿真程序如下： %二阶系统的最小二乘一次完成算法辨识程序，在光盘中的文件名：FLch3LSeg1.m u=[-1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,1,1]; %系统辨识的输入信号为一个周期的M序列 z=zeros(1,16); %定义输出观测值的长度 for k=3:16 z(k)=1.5*z(k-1)-0.7*z(k-2)+u(k-1)+0.5*u(k-2); %用理想输出值作为观测值 end subplot(3,1,1) %画三行一列图形窗口中的第一个图形 stem(u) %画输入信号u的径线图形 subplot(3,1,2) %画三行一列图形窗口中的第二个图形 i=1:1:16; %横坐标范围是1到16，步长为1 plot(i,z) %图形的横坐标是采样时刻i, 纵坐标是输出观测值z, 图形格式为连续曲线 subplot(3,1,3) %