第七章-频响函数的估计

7. 频响函数的估计（相干分析）7.1. SISO系统的频响函数及其估计对于SISO系统，其频响函数的估计有很多计算方法，主要的有三种估计式。在没有噪声污染的情况下，它们的估计是等价的。但是实际上，由于不可避免的存在噪声，三种估计有所差异。本节讨论在主要的三种噪声污染下，三种传统估计式与真值之间的误差。7.1.1. 随机激励下的频响函数考虑一个SISO时不变线性系统，其频率响应函数为。设随机输入和响应信号分别为和，其傅立叶变换分别为和，则有上式两端乘以，取时间平均及集合平均，并注意与平均无关，则即如果不为零，则可得系统的频响函数的第一种计算式同样，如果在系统输入/出频谱式两端乘以，取时间平均和集合平均，得如果不为零，则可得系统的频响函数的第二种计算式将系统输入/出频谱式两端取共轭，得乘以原输入/出频谱式，并去时间平均和集合平均，得可得系统的频响函数的幅值计算式7.1.2. 频响函数的估计方法考虑一个SISO时不变线性系统，其频率响应函数为。设系统的实际输入和响应信号分别为和，其傅立叶变换分别为和，它们的测量信号分别为和，其傅立叶变换分别为和。(1) 输出端噪声的影响若只有输出端受到噪声信号的污染，并设它与系统的和无关。则有A、 第一估计式根据第一估计式的定义，有由于噪声与激励，亦即无关，故根据大数定律，平均次数足够多时为零，则结果表明：只有输出端响应受到噪声污染时，通过平均，根据第一估计式得到的频响函数估计是实际频响函数的无偏估计。B、 第二估计式根据第二估计式的定义，有由于噪声与响应，以及激励亦即无关，故根据大数定律，平均次数足够多时、和都为零，则式中，可见，只有输出端响应受到噪声污染时，通过平均，根据第二估计式得到的频响函数估计是实际频响函数的有偏估计，是过估计。C、 第三估计式根据第二估计式的定义，有由于噪声与响应无关，故根据大数定律，平均次数足够多时和都为零，则式中，可见，只有输出端响应受到噪声污染时，通过平均，根据第三估计式得到的频响函数幅值的估计是实际频响函数幅值的有偏估计，也是过估计。因此，只有输出响应受到噪声污染时，频响函数的三种估计式有如下关系：(2) 输入端噪声的影响若只有输入端受到噪声信号的污染，并设它与系统的和无关。则有A、 第一估计式根据第一估计式的定义，有结果表明：只有输入端激励受到噪声污染时，通过平均，根据第一估计式得到的频响函数估计是实际频响函数的欠估计。B、 第二估计式根据第二估计式的定义，有可见，只有输入端激励受到噪声污染时，通过平均，根据第二估计式得到的频响函数估计是实际频响函数的无偏估计。C、 第三估计式根据第二估计式的定义，有可见，只有输入端激励受到噪声污染时，通过平均，根据第三估计式得到的频响函数幅值的估计是实际频响函数幅值的有偏估计，也是欠估计。所以，只有输入激励受到噪声污染时，频响函数的三种估计式有如下关系：(3) 输入/输出端复合噪声的影响若系统的输入/输出端分别受到噪声信号/的污染，并设它们与系统的和无关。则有A、 第一估计式根据第一估计式的定义，可以推得可见，在系统输入激励和输出端响应都受到噪声污染时，通过平均，根据第一估计式得到的频响函数估计是实际频响函数的欠估计，并且与响应信号中的噪声无关。B、 第二估计式同理，根据第二估计式的定义，可以推导出上式表明，在系统输入激励和输出端响应都受到噪声污染时，通过平均，根据第二估计式得到的频响函数估计是实际频响函数的有偏估计，是过估计，并且与系统的输入噪声无关。C、 第三估计式根据第二估计式的定义，可以推得进一步分析上式，可以得出下述结论：在系统输入激励和输出端响应都受到噪声污染时，通过平均，根据第三估计式得到的频响函数幅值的估计是较和都更为接近实际频响函数幅值的估计。因此，在系统输入激励和输出端响应都受到噪声污染的情况下，频响函数的三种估计式有如下关系：综上所述，根据上述所有分析过程，可以进一步获得以下几点结论：（1） 频响函数的第一估计式可以抑制输出响应噪声，第二估计式则可以抑制系统的输入激励噪声；（2） 在各种噪声情况下，第一、第二估计式给出了实际频响函数的范围：；（3） 在受到噪声的情况下，第三估计式虽然与实际频响函数最为接近，但是由于它只能给出频响幅值的估计，无法获得相频特性信息，故在实际中很少使用；7.1.3. 常相干函数根据上面的第二点综合结论，可知：将第一、二无估计式代入上式，可知：上式成为相干函数，记为根据互谱不等式，显然有，如果测试信号不受噪声污染，等于，则等于1。如果测试信号完全被噪声淹没，将趋于零。所以，相干函数反映了测试信号受噪声污染的程度。相干函数数值越大，说明噪声污染越小。在无噪声的情况下，从这个意义上而言，对于两个任意信号，相干函数反映了两个信号的相干关系。如果两个信号的相干函数等于1，则说明这两个信号之间必定存在类似理想无噪声的线性系统输入/出之间的相干关系。值得指出的是，线性系统这一点非常重要。如果系统具有非线性特性，其输入输出之间的关系不能由传递函数完全表达，相干函数将不等于1，会降低。通常，在实际情况下，当相干函数大于零但小于1时，有下面四个主要原因的一个或多个成立：(1) 测量中有外界噪声；(2) 两个信号之间的系统不是线性的；(3) 系统有测量输入之外的其它激励源；(4) 谱估计中有偏度误差，主要包括采样混叠、频率泄漏、分辨率不够，以及计算误差等等。当系统只有输出响应噪声时，而且与系统的输入/输出都不相关时，频响函数的第一估计式是系统频响函数的无偏估计，则系统的输出谱可以表示为通常将相干函数与输出自功率谱的乘积称为相干输出谱。上式表明，相干输出谱可以解释为输出功率谱中，由线性系统形成的部分。输出噪声功率谱可以通过下式计算它可以解释为线性运算时没有计及的由于上述四个原因形成的那部分输出。上式可以改写成因此，相干函数的下降可以理解为在各个频率点的输出信噪比的下降。7.2. MIMO系统的频响函数及其估计7.2.1. MIMO系统的频响函数考虑一个MIMO多自由度的时不变线性系统，其频率响应函数为。在m个输入平稳随机激励的作用下，其n个输出响应也是平稳的随机过程。设随机输入和响应信号分别为和，其傅立叶变换分别为和，则有展开写成矩阵形式，即上式两端右乘，取时间平均及集合平均，并注意与平均无关，则即式中互谱矩阵为自谱矩阵为如果不为零，则可得系统的频响函数的第一种计算式同样，如果在系统输入/出频谱式两端右乘，取时间平均和集合平均，得从上式可以得知，当输入/输出数目不相等时，不是方阵，与SISO系统，不能利用上式估计频响函数矩阵。将系统输入/出频谱式两端取转置共轭，得以上式右乘原输入/出频谱式，并取时间平均和集合平均，得从上式可以看出，对于一般情况下，上式也很难用于MIMO系统的频响函数矩阵的估计计算。所以，对于MIMO系统，通常采用第一种估计方法去计算频响函数矩阵。采用这种估计方法的缺点是，在某些频率点处，有可能出现输入自谱矩阵条件数很低时，这时，这种估计将会引起很大的估计偏差；当输入自谱矩阵奇异时，估计失败。所以，在实际应用中，通常采用估计算法。一般认为，估计与上述估计在理论上等价。7.2.2. 频响函数的MISO估计方法对于MISO的情况的系统频响函数估计相对简单，在理论上可以进行分析，有助于理解频响函数和相干函数在MIMO的应用。在MISO的估计模型中，由于输入情况很复杂，通常假定输入没有噪声，所有噪声污染都在输出端引入一个输出噪声一并考虑。对于这种情况下，系统的输入/输出频谱特性之间有，其中，由上式，噪声项及其复共轭可以写成于是上面两式相乘，得上式取时间平均和集合平均，得这就为估计频响函数确定了噪声功率谱函数。现在定义最优系统为这样的系统：所取得的频响函数，应能使噪声功率谱函数取得最小。令或者可以得到于是如果进一步假定输入之间两两不相关，则上式可以简化为这就意味着整个系统只是简单的SISO系统的线性组合。7.2.3. 重相干函数根据SISO系统的相干函数的定义，可以知道对于MISO和MIMO系统，通常的相干函数，即常相干函数不能用来判别所估计的MISO系统的频响函数的好坏。根据MISO系统的特点，我们可以将所有输入的线性叠加定义为一个理想线性输出，模仿只有输出噪声的SISO系统的相干分析来分析。定义理想的线性输出与实际测量之间的常相干函数为MISO系统的重相干函数，记作记号的含义是，在所有输入都作用下的理想输出与实际输出之比。理想输出功率谱的计算为与SISO系统一样，这个理想输出谱通常也称为重相干输出谱。对于输入之间相互独立的情况，重相干函数是所有输入输出之间的常相干函数的和。重相干函数提供了一种方法来判别MISO系统的频响函数的估计的准则。但应注意，利用重相干函数，不能判断单一输入与输出之间的传递特性估计的好坏。为了弥补这个不足，引入条件谱估计的概念，并进而定义偏相干函数。7.2.4. 条件谱函数对于前面定义的MISO系统，定义条件功率谱为式中，条件频谱的计算为各自的频谱函数，按SISO的输出噪声无偏估计式估计信号之间的频率特性，再减去该输入的影响。计算式为：上式表明：所谓条件功率谱，是按照SISO的线性系统理论，在系统输入和输出功率谱中，扣除了某一或某些输入引起的分量。这种扣除，也可以推广到以多个输入为条件的条件谱。7.2.5. 偏相干函数偏相干函数的是根据扣除某一或者某些输入为条件，求解的条件功率谱，并按照SISO的常相干函数来定义的：根据偏相干函数的上述定义，我们不难理解，偏相干函数主要用于判定某一输入或某些输入对输出的线性贡献，条件是根据先验知识，先扣除其他的输入的线性贡献。从而，据此可以判别所估计的传递函数是否有偏。