二、维纳滤波
维纳滤波与卡尔曼滤波都是最小平方偏差滤波,但它们依据的已知条件、计算方法及适用范围等都有所不同。
[维纳滤波准则] 假设希望得到序列z(t)(t只取某些整数值),而实际得到的是序列x(t)。于是设计一个脉冲响应为h(t)的线性定常系统,使得它的输入量为x(t),而输出量
与z(t)的偏差在某种准则下尽可能地小。
所谓维纳滤波问题,就是适当地取h(t),使得
当x(t)是确定序列时,
当x(t)是随机序列时,
[单路最小平方滤波] 假设滤波因子h(t)是由s+1个等跨距的系数组成的序列,如
s+1称为滤波因子的长度,则最小平方滤波因子h(t)满足
或用褶积表示为
其中
称为输入x(t)的自相关函数(第十六章§3),
称为z(t)与x(t)的互相关函数。
滤波输出与希望输出的符合程度可用标准化均方误差:
来度量,其中
显然。当时,滤波作用最好;当时,滤波作用最差。
当滤波因子的长度是无限的时候,可以得到类似的结果。用,,分别表示h(t),的频率谱*,则根据两个函数的褶积相当于其谱的乘积,可得到其频率特征
它表明用互相关函数与自相关函数的频率谱可以决定滤波因子的频率谱。
[多路最小平方滤波] 多路最小平方滤波方法是利用多路讯号的重复性,从而提供更多的有益讯息。
设为n路输入,为m路希望输出,写成矩阵分别为
为多路滤波因子,写成矩阵为
-------------------------------
离散时间序列的频率谱定义为
其中
当离散时间序列只有有限项时,可把它补上许多0,成为无穷序列而用上述定义。
是m路滤波输出,写成矩阵为
若每路输入的长度为k+1,所有的滤波因子的长度为s+1,则每路的滤波输出的长度为k+s+1。
所谓多路最小平方滤波,就是适当选取多路滤波因子(矩阵)H(t)使总均方差
当是确定序列时,
当是随机序列时,
多路最小平方滤波因子满足下列方程
(1)
其中
是多路最小平方滤波因子,是矩阵,
为输入的自相关函数方阵,称为多路自相关矩阵,的转置矩阵,
为希望输出和输入的互相关函数的矩阵。
方程(1)是以滤波因子 为未知数的mn(s+1)个线性联立方程。它的解即所求的多路最小平方滤波因子。
滤波输出与希望输出之间的符合程度,可用标准化均方误差
来度量,显然,。当时,多路最小平方滤波输出与希望输出符合得最好;当时,符合的程度最差。
当滤波因子的长度是无限时,最小平方滤波因子满足下列方程
用频率谱表示,则有
=
上述方程的解就是所求的多路最小平方滤波因子。