《系统辨识M序列》由会员分享，可在线阅读，更多相关《系统辨识M序列（5页珍藏版）》请在装配图网上搜索。
1、1.M序列的产生因为Np=15，所以应该用4个移位寄存器。连接图如下：仿真图如下:134671.U111-14-1f1晦1OO1Dn11D1OL11.二1.1OO1OO1.1OI1O113L111OO1OO1.101O1la1111DO110O11O1町.1OO10111O11110表如下:用程序实现如下:%定义移位寄存器Xi初态（1000）,Yi为移位寄存器各级输出X1=0;X2=0;X3=0;X4=0;m=32;%置M序列总长度fori=1:m%1#Y4=X4;Y3=X3;Y2=X2;X4=Y3;X3=Y2;X2=Y1;X1=xor(Y3,Y4);%异或运算ifY4=0Y1=X1;U(i)
2、=0.03;elseU(i)=Y4*(-0.03);endendM=U;figure(1)%绘M序列图i1=ik=1:1:i1;stairs(M,g-)xlabel(k)ylabel(M序列)title(移位寄存器产生的M序列)2从图形得到数据丫=2.051.801.841.791.080.680.440.801.912.382.472.513.052.691.941.821.822.032.031.030.680.520.861.782.502.502.323.282.822.042.011.671.701.821.040.590.380.811.912.552.282.563.132.70
3、2.01;3相关函数s=0;c=0;%设置中间变量，初始化R=zeros(1,15);r=zeros(1,32);%初始化相关数组fori=1:15%计算相关函数forj=16:32r(j)=U(j-i+1)*Y(j)/30;s=s+r(j)endR(i)=s;s=0;endc=(-1)*R(15);%求补偿量cfori=1:15%求脉冲响应g(k)g(i)=15*(R(i)+c)/(4*0.03A2*16)endfigure(3)plot(R)xlim(1,15)ylabel(Rym)title(互相关函数);互相关函数4脉冲响应曲线holdonfigure(4)plot(g)xlim(1,
4、15)ylabel(g)title(脉冲响应);脉冲响应5阶跃响应曲线holdonp=0;fork=1:15%用积分法求阶跃响应forl=1:kp=g(i)+pendv(k)=p;p=0；endfigure(5)plot(v)title(阶跃响应)xlim(1,15);holdon阶跃响应6得出加热炉的数学模型选择采样间隔T=1，于是得g(t)在t=1:1:15这15个时刻对应的值为：0.2813,0.3177,0.6094,3.896，9.021,12.76，15.97,18.14,19.82,21.21,21.79,21.7421.86,22.1,22.1用Z变换的方法来求解系统的传递函数，设0ZG(z)-(1)zaiZa2由实验曲线阶跃响应曲线可得其Z变换G(z)(采样间隔为T=1)为0.2813/z+0.3177/zA2+0.6094/zA3(2)要使(1)式与(2)式相等，则应取6=0.3201a1=-1.2158a?=-1.2369则可以得到：0.3201sGs2s-1.2158s1.2369
系统辨识M序列
