![雷达数据处理及应用(第四版)](https://wfqqreader-1252317822.image.myqcloud.com/cover/900/47379900/b_47379900.jpg)
3.2.2 滤波模型
在所有的线性形式的滤波器中,线性均方估计滤波器是最优的[17-20]。线性均方误差准则下的滤波器包括:维纳滤波器和卡尔曼滤波器,稳态条件下二者是一致的,但卡尔曼滤波器适用于有限观测间隔的非平稳问题,它是适合于计算机计算的递推算法。
2.5节给出静态(非时变)情况下随机向量x的最小均方误差估计为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_1.jpg?sign=1740019185-IghEi1DL85xIg7jhlefGApcj0fmU3x6s-0-3991e7d3f5bcd2e824469e4f81a012b7)
其对应的条件误差协方差矩阵为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_2.jpg?sign=1740019185-lyGesaD6Abz6G0u7q1mD83ohCQw8vBbB-0-45a5eff424c1f0d80430c61c1f3c3a20)
类似地,动态(时变)情况下的最小均方误差估计可定义为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_3.jpg?sign=1740019185-vToB5BO4mEqy9qtOHNP2cYh6Xmd6ECkT-0-ff511e0dea00f0e7f181de25d15b3c15)
式中
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_4.jpg?sign=1740019185-8WVXDHKlh1qC5isCXaJEi0j8QfrGGGo8-0-a858a516cb3cccfb9e6f0f6c8faed53c)
与式(3.41)相伴的状态误差协方差矩阵定义为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_5.jpg?sign=1740019185-ltrFQSpQWUTwNwil4yfFGFPH8XSeppr5-0-5740cb126ed90a427e5a58b8920892a2)
把以Zk为条件的期望算子应用到式(3.31)中,得到状态的一步预测为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_6.jpg?sign=1740019185-RlVdAun8Im306j7lg2ITBxhfp3K1NDvp-0-7115539bd69b9bfcb608d833aa83c59a)
预测值的误差为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_7.jpg?sign=1740019185-NvIij94OvZ9McW5x4pd51UMb1dh8Mqxn-0-e48329956a0568c3caa6fe691447be3a)
一步预测协方差为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_8.jpg?sign=1740019185-c4T2yFfsoh74p86QfSlqUrvbNdODYu0D-0-92993323146ecc402018cf6af4d2b20e)
注意:一步预测协方差P(k+1|k)为对称阵,它可用来衡量预测的不确定性,P(k+1|k)越小则预测越精确。
通过对式(3.32)取在k+1时刻、以Zk为条件的期望值,可以类似地得到量测的预测是
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_9.jpg?sign=1740019185-BxXKPsusDgJfznv35FgzJ8NIDhEiZeCJ-0-258c04956dd3947190b7862b70947642)
进而可求得量测的预测值和量测值之间的差值为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_10.jpg?sign=1740019185-v2aEl4tlqIbvQzgAZwWOpWYTJ5fWrFUY-0-b41aaabde5d601f42e7bb289f4922534)
量测的预测协方差(或新息协方差)为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_11.jpg?sign=1740019185-mRfvyswjKrxDRq6yfrZyUUeW4JvyxZlU-0-7ef4fa215ee1eac5787049936334f62c)
注意:新息协方差S(k+1)也为对称阵,它用来衡量新息的不确定性,新息协方差越小,则说明量测值越精确。
状态和量测之间的协方差为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_12.jpg?sign=1740019185-GQfzJ3Um6Dd6fqgkndyfbLfK9L475HdX-0-80a2b3777479bbd56fd7cc9bfcb4bd54)
增益为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_68_13.jpg?sign=1740019185-j07yP94rLvnFIXBxJPd4ghethglq6XwI-0-6c96dbc3f116ebd0ce72ef1d3c54040d)
增益的大小反映了最新观测信息对状态估计量的贡献大小。
进而,可求得k+1时刻的状态更新方程为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_69_1.jpg?sign=1740019185-jd7QDOtqv0NruZsT5yxi6Ubj8LsgtbUO-0-7cfba327d88e9f72ebff8339cd766ebd)
式中,v(k+1)为新息或量测残差,即
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_69_2.jpg?sign=1740019185-k9N682OzYra9dQEE0C5Ro25Udztw4x8h-0-eff4b82f4ef01e0de6c1b7a89624482e)
式(3.52)说明k+1时刻的估计等于该时刻的状态预测值
再加上一个修正项,而这个修正项与增益K(k+1)和新息有关。
协方差更新方程为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_69_5.jpg?sign=1740019185-LSXOFwkXVG8BsRQ9R1Z47CkVgt6gi1Fh-0-d31bc6bcfa63e115a18994cdbde81941)
式中,I为与协方差同维的单位阵。式(3.57)可保证协方差矩阵P的对称性和正定性。
滤波器增益的另一种表示形式为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_69_6.jpg?sign=1740019185-GY8nTdFkZisYcHlG1pSko3XFYlSuco6P-0-98dd8504efe5fd03315520a159d0aa69)
卡尔曼滤波除了系统噪声和量测噪声为高斯白噪声且已知其二阶矩之外,不需任何其他条件,因而完全适用于非平稳、多维的随机序列的估计问题。图3.3给出了卡尔曼滤波所包含的方程及滤波流程,而卡尔曼滤波算法单次循环流程图如图3.4所示,其余的依次类推。
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_69_7.jpg?sign=1740019185-hTKFSI9ufl2cFeETgMK3YM412YJApcKv-0-07fa033bf6f26e2ec2150a2781d3bc4f)
图3.3 卡尔曼滤波算法框图
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_70_1.jpg?sign=1740019185-hoyb48sq5k6Xn1AuwgrAv8mzGvKgOBRd-0-c7915deea5ede0fe0453037a88b31ee9)
图3.4 卡尔曼滤波算法单次循环流程图