用脉冲响应求传递函数PPT精选文档

1、11.连续系统的传递函数任何一个SISO系统都可以用差分方程来表示。若系统的输入为 函数，则输出为脉冲响应函数g(t)。因为 函数只作用于t=0，而在其他时刻系统的输入为0，所以系统的输出是从t=0开始的脉冲响应函数g(t)。如果采样间隔时间为T。并设系统可以用n阶差分方程表示，则：用脉冲响应来求解传递函数用脉冲响应来求解传递函数0100( )()()0ng ta g tTa g tnT( ) t( ) t等式中a1,a2,an为待定的n个常数。2根据上式，将时间依次延迟T，可以得到：100010001000102000()()( )(2 )(1) )()(3 )(2) )(2 ) ()(1)

2、 )(2)(1) )nnnna g tTa g tnTg ta g tTa g tnTg tTa g tTa g tnTg tTa g tnTa g tnTa g tnTg tnT联立求解上述n个方程，就可以得到差分方程的n个系数a1,a2,an。3等式中s1,s2,sn和c1,c2,.,cn为待求的2n个未知数。对上式求Laplace反变换，得到脉冲响应函数： 任何一个线性定常系统，如果其传递函数G(s)的特征根为s1s2sn，则其传递函数可以表示为：121212()()()12(2 )(2 )(2 )12()()()12()(2 )()nnnst Ts t Tst TnstTs tTstT

3、nst nTs t nTst nTng tTc ec ec eg tTc ec ec eg tnTc ec ec e1212( )nncccG sssssss1212( )ns ts ts tng tc ec ec e4要使上式为成立，应令方括号内的值为0，即：将上面等式带入到下列脉冲响应的差分方程中2121()0, 1,2iiisTsTsTnna ea ea ein211212111121()1()1()0nnns tssTsTs Ts Tnnnns Ts Tnttnsnaea eaeaceceaeaeceeisTex1( )()()0ng ta g tTa g tnT得到：令，则可以得到：

4、110nna xa x5解方程可以得到x的n个解x1,x2,xn。设：1212lnlnln,nnxxxsssTTT1212,ns Ts Ts Tnex exex至此可以得到s1,s2sn，下面求解c1,c2cn。121 1221111 122(0)( )(1) )nnnnnnnngcccg Tc xc xc xg nTc xc xc x6例：有一个三阶系统，脉冲响应数据如下： t 0 1 2 3 4 5 g(t) 0 1 2 2 4 0试求解该系统的线性定常脉冲传递函数：312123( )cccG sssssss7等式中。因而有2. 离散系统的脉冲传递函数( )(),0,1,2.g ig iT

5、i110111()1nnnnbb zb zG za za z112()(0)(1)(2)G zggzgz设系统脉冲传递函数形式为：根据脉冲传递函数的定义可以得到：1120111(0)(1)(2)1nnnnbb zb zggzgza za z8进一步得到:1110111(1)12(0)(1)(1)(1)(2 )(2)(0)( )()niinnnnniniiingbb zb zgga gzg na g nizzzna g nigna gni 90112211-211 0 0 0 0(0) 1 0 0 0(1) 1 0 0(2) ( ) 1nnnnbgbagbaagg nba aaa11(1) (2

6、) ( )(2) (3) (1) ( ) (1(1)(2) () (2 -1)2 )nnggg nggg ngag nag ngnang ngn令上式两边z-i的同次项系数相等，可以得到：10例：设采样间隔时间为0.5s，系统的脉冲响应序列g(k)如下表所示，求系统的脉冲传递函数。t00.050.10.150.200.250.3k0123456g(k)07.5159.4918.5645.9312.8460.145110111()1nnnnbb zb zG za za z11例：有一个三阶系统，脉冲响应数据如下： k 0 1 2 3 4 5 6 g(t) 0 1 4 2 6 2 2试用Hanke

7、l矩阵法求解该系统的脉冲传递函数。12第七章系统阶次的辨识系统的阶次，对传递函数而言，指极点个数；对于状态空间而言，是指最小实现的状态个数；本章讨论单输入单输出系统的阶次辨识问题，主要介绍F检验法和AIC准则这两种基本的阶次辨识方法；阶次辨识和参数估计两者是互相依赖的，参数估计时需要已知阶次，而辨识阶次时又要利用参数估计值，两者密不可分。13如何根据脉冲响应的采样值来判定模型的阶次？7.1. 7.1. 根据根据HankelHankel矩阵判定模型的阶次矩阵判定模型的阶次1112122 ( , ) kkk lkkk lk lk lklggggggH l kggg 已知系统的脉冲响应序列g0,g1

8、,gN，定义Hankel矩阵H(l,k)为：我们根据Hankel矩阵的秩来判定系统模型阶次。14定理：若Hankel矩阵的维数l大于系统的阶次n，则Hankel矩阵的秩等于系统的阶次n。当Hankel矩阵维数l=n+1时，对于所有的k，Hankel矩阵的行列式为零。当我们对于每个k值以及不同的维数l值，计算Hankel的行列式，就可以判定模型的阶次n。15( , )(1, )lH l kDH lk实际上，由于噪声存在，当维数l=n+1时，这些行列式的值并不恒等于零，但会突然变小。我们必须引入某个准则，以确定显著性水平。有一种方法是对于每一个不同的维数l值，计算Hankel矩阵的行列式的平均值。

9、然后对于不同的l值，比较行列式比值Dl。Dl值为最大时的维数l值，就是系统模型的阶次。16以自相关系数 作为Hankel矩阵的元素，再按新的Hankel矩阵来确定矩阵的秩。同样，由于噪声的影响，所得的行列式也不恒等于零。( )g 01( ), =0,1,2.1NgkkkRg gN另一种方法是根据脉冲响应序列，求出它们的自相关序列的估计值，以及自相关系数值。( )( )(0)gggRR 171.0,0.80,0.65,0.54,0.46,0.39,0.35,0.31,0.28,0.26,0.24,0.23,0.22,0.21,0.20,0.19,0.19,0.18,0.18,0.18,0.17,

10、0.17,0.17,0.16,0.16,0.15,0.15,0.15,0.15,0.14,0.14,0.14,0.13,0.13,0.13,0.13,0.12,0.12,0.12,0.12,0.12,0.11,0.11,0.11,0.11,0.10,0.10,0.10试判定该模型的阶次。例：已知系统的脉冲响应序列g(k)为18第一种方法，求得各Hankel矩阵行列式的平均值，以及行列式比分别为：矩阵H(2,k)行列式的平均值为0.00087872矩阵H(3,k)行列式的平均值为-0.00029311矩阵H(4,k)行列式的平均值为-3.21410-7矩阵H(5,k)行列式的平均值为-5.709

11、10-9D2=2.998, D3=913.1, D4=64.2因此，可以确定系统的阶数为3。19第二种方法，求出脉冲响应序列的相关系数为：以为元素构造Hankel矩阵并计算Hankel矩阵的行列式，得到：detH(2,0)=0.014937detH(3,0)=-1.28210-5detH(4,0)=-5.810-8由行列式的值可知，系统模型的阶次可以定为3阶，也可以定义为2阶。因为detH(3,0)已经很小了。012345671, 0.88052126, 0.79025506,0.72231277, 0.67060564, 0.62999127,0.60107303, 0.57697552,

12、.i20用最小二乘法求出参数的估值，则目标函数为：1. 阶和目标函数7.2 根据残差特性判定模型的阶次11TJN 11() ( )() ( )( )A zy kB zu kk考虑系统模型为：1111() ( )() ( )() ( )A zy kB zu kC zk如果系统模型为：则目标函数为：2111TJN 21当n=1,2,时，J1(n)和J2(n)随着n的增加而减小。如果n0为正确的阶次，则n=n0-1时，J(n)出现最后一次陡峭的下降，n再增大，则J(n)保持不变或者只有微小的变化。如果阶次已给定，估计参数，则要求J1和J2最小值。如果阶次未知，则估计参数个数就未知，也要求J1和J2取

13、极小值。那么，当阶次递增时，J1和J2的变化规律如何呢？对于不同阶次，目标函数为：11( )() ()TJ nYYN2111( )() ()TJnYYN22假设检验与参数估计区别 参数估计和假设检验都是统计推断的两个组成部分，都是利用样本对总体进行某种推断，但推断的角度不同。 参数估计是在总体参数未知的情况下用样本统计量估计总体参数。 假设检验是先对总体参数提出一个假设，然后利用样本信息去检验这个假设是否成立，如果成立，就接受这个假设，否则就放弃。23在实际工作中，前人对某些问题得到初步的结论。这些结论可能正确、可能错误。若视这些结论为假设，问 题 在 于 我 们 是 否 应 该 接 受 这

14、些 假 设 呢 ？例：我们对某产品进行了一些工艺改造，或研制了新的产品。要比较原产品和新产品在某一项指标上的差异，这样我们面临选择是否接受假设“新产品的某一项指标优于老产品”。我们必须作一些试验，也就是抽样。根据得到的样本观察值来作出决定。nxxx,21 假设检验问题就是根据样本的信息检验关于总体的某个假设是否正确。假设检验的思想24假设检验的方法先介绍一条所谓实际推断原理(小概率原理)。通过大量实践，人们对小概率事件(即在一次试验中发生的概率很小的事情)总结出一条原理：小概率事件在一次试验中几乎不会发生小概率事件在一次试验中几乎不会发生并称此为实际推断原理，其为判断假设的根据。在假设检验时，

15、若一次试验中小概率事件发生了，就认为是不合理的。小概率事件在一次试验中发生的概率记为，一般取1 . 001. 005. 0在假设检验中，称为显著水平、检验水平。25 即先对所关心的问题提出原假设H0 ,然后运用样本信息看在H0成立的条件下会不会发生矛盾。最后对H0成功与否作出判断：若小概率事件发生了，则否定H0；若不发生，则接受H0。概率反证法的逻辑是：如果小概率事件在一次试验中居然发生，我们就以很大的把握否定原假设。假设检验使用的方法是概率论的反证法反证法：26某日开工后为检验包装机是否正常。随机地抽取它所包装的糖9袋，称得净重为(公斤)：(=0.05)例：某车间用一台包装机包装葡萄糖。包得

16、的袋装糖重是一个随机变量X，且。当机器正常时，其均值为=0.5公斤，标准差=0.015公斤。0.497;0.506;0.518;0.524;0.498;0.511;0.520;0.515;0.512问机器是否正常?),(2NX0:0.5H1:0.5H那么，如何判断原假设H0是否成立呢？即看在=0.5的条件下会不会产生不合理现象。解解：先提出假设27即在= 0.5的条件下)015. 0,5 . 0(2NX则21910.015() (0.5,)99XXXN921,XXX为X的一个样本，它们是随机变量，00 1XZNn（ ， ）于是有统计量问Z大到什么程度可以否定H0呢？这就要确定一个否定H0的值/

17、2kz为此令：/2P zz/20.50.015 30.05PXz即2802xZzn即02xZzn022.21.96xZzn.0常，认为包装机工作不正于是拒绝H511. 0)512. 0497. 0(91x/20.0251.96zz0 x)(x22z2z拒绝H0接受H0接受域接受域29假设检验的步骤由实际问题提出原假设 和备择假设 ；确定适当的检验统计量，并在原假设为真的条件下确定该统计量的分布；根据问题的要求规定显著性水平 (一般题目中会给定)，从而得到拒绝域；由样本观测值计算检验统计量的值，看是否属 于拒绝域，从而对原假设 作出判断。0H0H1H30 Astrom(1968)提出的F检验法，

18、引入一个假设检验，将模型阶次的判定问题归结为当阶次从n1增加到n2时，J2(n2)较J2(n1)下降是否显著的问题。 Astrom证明，当N足够大时，若n2n1=n成立，则：2. 确定阶的F检验法122222()(2)NJnNn122122212()()(22 )NJnJnnn122112()(2)NJnNn且J2(n2)和J2(n1)-J2(n2)是相互独立的随机变量。31存在显著性水平 ，当 时，n2n1=n成立。否则，n比n1和n2大。这里取 ，则引入统计量t也就是说，当n2n1=n时，统计量t服从自由度分别为2n2-2n1与N-2n2的F分布。212222221221()()2(222

19、2,)2)JnJnNnFnnnnNtnJntt12,1nn nn 222( )(1)22 ( ,1)(2,12)22)JnJnNnt n nFNnnJ32当时，则成立， 即为对象的阶令，阶次逐渐增加，用F检验判断阶次增加，减少。令 ( ,1)t n nt2( )Jn(2,22)tFNn nn1,2,.n nn例如，取置信度，在100,200,400和时，查表可以得到F分布的值：F(2,100)=3.09; F(2,200)=3.04; F(2,400)=3.02; F(2,)=3.00;0.05因此，当 =100时，若统计量，则接受假设，即认为 下降也不显著，判定系统阶次为 ，否则否定假设，继续增加阶次并考察统计量22Nn ( ,1) 3.09t nn 2 ( )J n n ( ,1)t n n n33其中 表示n阶模型未知参数的个数， 表示参数的极大似然估计值， 为似然函数，反映拟合精度。称该准则为AIC准则。Akaika证明了使为极小的阶，即系统的阶。 Akaika(1972)提出一种具有客观标准的阶次判别方法，所采用阶次判定准则为：2. 确定阶的AIC准则n( , )Ln()2ln( , )2AIC nLnn ()AIC n n34这里对AIC准则作一个定性的解释。设系统的阶次为n，当阶次