系统辨识M序列.doc

《系统辨识M序列.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、1.M序列的产生因为Np=15，所以应该用4个移位寄存器。连接图如下:仿真图如下：表如下：用程序实现如下：X1=0;X2=0;X3=0;X4=0;%定义移位寄存器Xi初态（1000），Yi为移位寄存器各级输出m=32;%置M序列总长度fori=1:m%1#Y4=X4;Y3=X3;Y2=X2;Y1=X1;X4=Y3;X3=Y2;X2=Y1;X1=xor(Y3,Y4);%异或运算ifY4==0U(i)=0.03;elseU(i)=Y4*(-0.03);endendM=U;figure(1)%绘M序列图i1=i

2、k=1:1:i1;stairs(M,'g-')xlabel('k')ylabel('M序列')title('移位寄存器产生的M序列')2．从图形得到数据Y=[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.03^2*16)endfigure(3)plot(R)xlim([1,15])ylabel('Rym')tit

4、le('互相关函数');4.脉冲响应曲线holdonfigure(4)plot(g)xlim([1,15])ylabel('g')title('脉冲响应');5.阶跃响应曲线holdonp=0;fork=1:15%用积分法求阶跃响应forl=1:kp=g(l)+pendv(k)=p;p=0;endfigure(5)plot(v)title('阶跃响应')xlim([1,15]);holdon6.得出加热炉的数学模型选择采样间隔T=1，于是得g(t)在t=[1：1：15]这15个时刻对应的值为：0.2813

5、,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变换的方法来求解系统的传递函数，设(1)由实验曲线阶跃响应曲线可得其Z变换G（z）（采样间隔为T=1）为0.2813/z+0.3177/z^2+0.6094/z^3(2)要使（1）式与（2）式相等，则应取b1=0.3201a1=-1.2158a2=-1.2369则可以得到：