实验2相关分析法辨识脉冲响应

1、1实验题目:相关分析法辨识系统脉冲响应实验2 相关分析法辨识脉冲响应实验报告哈尔滨工业大学航天学院控制科学与工程系专业：自动化班级：姓名：籍洋日期： 2013年10月10日2.实验目的通过仿真实验掌握利用相关分析法辨识脉冲响应的原理和方法。3实验主要原理g( R)描y(i)一个单入单出线性定常系统的动态特性可用它的脉冲响应函数 述。线性系纺则 y(t) 0 g( )x(t)d上式两端同乘x(t1 T lim y(t)x(t t t 0 则)，)dt进而取时间均值，有1 T0 g( )lTm t 0 x(t )x(t )dtdo g()Rx()dRxy()这就是著名的维纳霍夫积分方程。1实验题目

2、:相关分析法辨识系统脉冲响应1实验题目:相关分析法辨识系统脉冲响应Rx( ) k ( ),Rx(则根据维纳霍夫积分方程可得Rxy()0 g( )Rx()d kg()如果输入是 白噪声，这时x(t)的自相关函数为)k ()1实验题目:相关分析法辨识系统脉冲响应1实验题目:相关分析法辨识系统脉冲响应g()Rxy()kx(t)、y(t)的值，并计算它们的互相关函数，即可求得脉或者这样，只要记录 冲响应函数g( T)而在系统有正常输入的情形下，辨识脉冲响应的原理图如下图所示。线性系统g(r)1实验题目:相关分析法辨识系统脉冲响应1实验题目:相关分析法辨识系统脉冲响应白枳分器延迟TMH4.实验对象或参数

3、下图为本实验的原理框图。系统的传递函数为G(s),其中K 120, Ti 8.3Sec,T2 6.2Sec ； u(k)和z(k)分别为系统的输入和输出变量;v(k)为测量白噪声，服从正态分布，均值为零，方差为：，记作v(k)N(0, 2)； go(k)为系统的脉冲响应理论值，g(k)为系统脉冲响应估计值，(k)为系统脉冲响应估计误差。系统的输入采用M序列(采用实验1中的M序列即可)，输出受到白噪，利用相关分析法计算声v(k)的污染。根据过程的输入和输出数据u(k), z(k)出系统的脉冲响应值?g(k)，并与系统的脉冲响应理论值go(k)比较，得到系统脉冲响应估计误差值k t/T? g(k)

4、(k)，当k 时，应该有(k)1、模拟过程传递函数G(s)，获得过程的输入和输出数据u(k), z(k)(采样时间取1秒)。(1)惯性环节其中，T为惯性环节的时间常数，K为惯性环节的静态放大倍数。若采样时间 记作To，则惯性环节的输出可写成：y(k) e To/Ty(k 1) TK(1 e To/T)u(k 1)TKT(eT/T 1) To) u(kLu(1)To传递函数G(s)仿真(串联)K 11G(s)T1T2 s 1/T1 s 1/T2K令K1=，则G (s)的表达框图为：T1T2u(k)K1x(k)1y(k)s 1/T1s 1/T22、互相关函数的计算1 (r 1)NpRMz(k) M

5、 u(i k)z(i)rN p i Np 1其中，r为周期数，i Np 1表示计算互相关函数所用的数据是从第二个周期 开始的，目的是等过程仿真数据进入平稳状态。(可分别令r =1、3,对比仿真 结果)3、c的补偿补偿量c应取Rmz(Np 1)，不能取Rmz(Np)。因为RMz(k)是周期函数,则有 Rmz(Np)= Rmz(0)，故不能取Rmz(Np)。4、计算脉冲响应估计值?脉冲响应估计值 g（k）脉冲响应估计误差g_Z 寫2 t RMz（k） c5 程序框图绘图结束6 .程序代码fun cti onsigma=resp on se(r)x=0,1,0,1,1,0;%初始化Np=2A6-1;

6、%游列长度a=1;%振幅t=1;for i=1:Np* (叶1)y(i)=x(6);temp=xor(x(5),x(6);for j=5:-1:1x(j+1)=x(j);endx(1)=temp;endfor i=1:Np* (叶1)if (y(i)=0)u(i)=a;elseu(i)=-a;endendK=120;T1=8.3;T2=6.2;T0=1;K1=K/T1/T2;x(1)=0;y(1)=0;for k=2:Np*(r+1)x(k)=exp(-T0/T1)*x(k-1)+T1*K1*(1-exp(-T0/T1)*u(k-1)+T1*K1*(T1*(exp( -T0/T1)-1)+T0

7、)*(u(k)-u(k-1)/T0;y(k)=exp(-T0/T2)*y(k-1)+T2*(1-exp(-T0/T2)*x(k-1)+T1* (T 2*(exp(-T0/T2 )-1)+T0)*(x(k)-x(k-1)/T0;%未经白噪声污染的输岀endv=white noise(1,le ngth(y);%产生白噪声z=y+v;%系统实际输岀for k=1:Npsum=0;for i=Np+1:(r+1)*Npsum=u(i-k)*z(i)+sum;endRmz(k)=1/(r*Np)*sum;endc=-Rmz(Np-1);%补偿量 c% 计算脉冲响应估计值for k=1:NpG(k)=N

8、p/(Np+1)*aA2*t)*(Rmz(k)+c);gO(k)=K/(T1-T2)*(exp(-k*t/T1)-exp(-k*t/T2);end%计算脉冲响应估计误差SUM1=0;SUM2=0;for k=1:Npe(k)=g0(k)-G(k);SUM1=e(k)A2+SUM1;SUM2=g0(k)A2+SUM2;endsigma=sqrt(SUM1/SUM2);step=0:Np-1;plot(step,Rmz(63),Rmz(1:62);hold on;plot(step,G(63),G(1:62),r)；plot(step,g0(63),g0(1:62),g)；legend(互相关函数

9、,脉冲响应估计值,脉冲响应理论值)end产生白噪声的函数：lenfunctionsig=whitenoise(sigma,len)%白噪声产生函数，sigma 为均方差,为白噪声序列数据长度a=65539;M=2147483647;b=100;x(1)=12345;r(1)=x(1)/M;%第一部分为产生0-1的均匀分布随机数for i=1:12*lenx(i+1)=mod(a*x(i)+b,M);r(i+1)=x(i+1)/M;end%n = 12;for i = 1 : lensig(i)=0;for j = 1 : nsig(i) = sig(i) + r(n * (i - 1) + j

10、);%第二部分产生正态分布，方差为sigma的随机序列endendsig = (sig - 12*0.5)*sigma;end7 实验结果及分析 另r = 1,白噪声均方差 sigma=0.5,运行命令： sigma=response(1)运行结果如下：得 sigma = 0.0416543210-1-2010 2030405060互相关函数脉冲响应估计值脉冲响应理论值另 r = 3,白噪声均方差 sigma=0.5,运行命令：sigma=response(3,0.5) 运行结果如下：得 sigma = 0.0373图像如下：76543210-1互相关函数脉冲响应估计值脉冲响应理论值-2 t010 203040506070比较发现r=1和r=3时产生的曲线基本相似，但是脉冲响应误差在r=3时更小。另r = 1,白噪声均方差 sigma=1,运行命令：sigma= response(1,1)运行结果如下：得 sigma = 0.0594互相关函数脉冲响应估计值脉冲响应理论值76543210-1-2-3 010203040506070比较第一个与第三个结果，发现脉冲响应的估计误差是随着输入白噪声标准差的增大而 增大的