资源描述：
《系统辨识实验答案》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验一辩识离散线性系统脉冲响应特性的相关分析法•、递推算法和批量算法流程图图1递推算法流程图图2批量算法流程图二、实验步骤：1、所产生的序列如下图所示:10.5-0.5S0炳序列L序列-1k2、输出如下图所示:输入为冲激函数时的输出输入为何序列时的输出输入为U^#
2、J时的输出3、用M序列作为输入，批量算法所得的辨识结果为用M序列作为输入，递推算法所得的辨识结果为:4、结论：通过比较上面两个图，可得这两种算法的辨识结果基本上是一致的，只是递推算法可以用于在线辨识，就这一点来说，递推算法优于一次完成算法。5、源程序如下：
3、%M序列的产生Xl=l;X2=0;X3=l;X4=0;%移位寄存器输入Xi初态(1010)m=100;%置M序列总长度m值fori=l:m%开始循环Y4=X4;Y3=X3;Y2=X2;Y1=X1;%Yi为移位寄存器各级输出，%在移位之前先将各自的输入传给输出X4=Y3;X3=Y2;X2=Y1;%实现移位寄存器的连接方式X1=xor(Y3,Y4);%实现异或运算ifY4==0%将输出'0'态转化为态U(i)=-1;elseU(i)=Y4;End%转换结束End%60次循环结朿M=U订=m;k=l:l:il;figure
4、(l);subplot(2,2,l)plot(k,U,k,U；rx*)%绘图xlabelCk*)ylabefM序列Jtitle(移位寄存器产生的M序列)%L序列的产生fori=l:m;t=(-DAi；ift==-1s(i)=0;elses(i)=l；endiL(i)=xor(U(i),s(i));ifiL⑴=0L(i)=-1;elseL(i)=l；endendi=l:l:m;subplot(2,2,2);plot(i,L)titleCL序列J%白噪声的产生A=6；X0=l；N=100；X0=l；M=255；fork=
5、l:NX2=A*X0;Xl=mod(X2,M);Vl=Xl/256;V(:,k)=V1-0.5;XO=X1;V0=V1;EndV2=VSubplot(2,2„3)k=l:l:N;plot(k,V)xlabel(，k')ylabel('V')title(X-0.5,+0.5)均匀分布的白噪声〃伙)J%脉冲序列的冲激响应a1=-1;a2=0.5;b1=1;b2=0.5;d=0;y(k)=-a1y(k-l)-a2y(k-2)+b1u(k-l-d)+b2u(k-2-d);Y0=0.2;X0=0.1；d=();Np=15;y(
6、l)二0;y(2)=0;%输入M序列的初始值yt(l)=O;yt⑵二0;%输入冲激函数吋初始值al=-l;a2=0.5;b1=1;b2=0.5;Y(l)=YO+y(1)+V(l);Y(2)二Y0+y(2)+V(2);impulse=[1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0];fork=3:l:2*Npy(k)=b1*U(k-1-d)+b2*U(k-2-d)-a1*y(k-1)-a2*y(k-2);Y(k)=Y0+y(k)+V(k);end
7、fork=3:l:2*Npyt(k)=b1*impulse(k-l-d)+b2*impulse(k-2-d)-a1*yt(k-1)-a2*yt(k-2);endfigure(2);k=l:l:2*Np;plot(k,yt(l:2*Np),k,yt(l:2*Np),'rx')xlabel('k')ylabeledLH1)titleC脉冲响应估计值?figure(3);k=l:l:2*Np;plot(k,Y(l:2*Np),k,Y(l:2*Np),,rx,)xlabel('k')ylabel(输出')title。脉冲响应
8、估计值J实验二离散线性系统参数估计的递推］1、基本递推最小二乘(RLS)算法。⑴程序流程图最小二乘递推算法辨识的程序流程图°eoo-eoo-uoo-coo
