2.卡尔曼滤波简介

  卡尔曼滤波是对随机信号作估计的算法之一。与最小二乘、维纳滤波等诸多估计算法相比,卡尔曼滤波具有显著的优点:采用状态空间法在时域内设计滤波器,用状态方程描述任何复杂多维信号的动力学特性,避开了在频域内对信号功率谱作分解带来的麻烦,滤波器设计简单易行;采用递推算法,实时量测信息经提炼被浓缩在估计值中,而不必存储时间过程中的量测量。所以,卡尔曼滤波能适用于白噪声激励的任何平稳或非平稳随机向量过程的估计,所得估计在线性估计中精度最佳。正由于其独特的优点,卡尔曼滤波在 20 世纪 60 年代初一经提出,立即受到工程界,特别是空间技术和航空界的高度重视。阿波罗登月计划中的导航系统设计和C-5A飞机的多模式导航系统的设计是卡尔曼滤波早期应用中最为成功的实例。随着计算机技术的发展,目前卡尔曼滤波的应用几乎涉及通讯、导航、遥感、地震测量、石油勘探、经济和社会学研究等所有领域。

卡尔曼滤波要解决的问题

  所谓滤波就是从混合在一起的诸多信号中提取出所需要的信号。

  信号是传递和运载信息的时间或空间函数。有一类信号的变化规律是既定的,如调幅广播中的载波信号、阶跃信号、脉宽固定的矩形脉冲信号等,它们都具有确定的频谱,这类信号称为确定性信号。另一类信号没有既定的变化规律,在相同的初始条件和环境条件下,信号的每次实现都不一样,如陀螺漂移、海浪、作水平飞行的飞机飞越山区时无线电高度表的输出信号、惯导系统的导航输出误差、GPS的SA误差等,它们没有确定的频谱,这类信号称为随机信号。

  由于确定性信号具有确定的频谱,所以可根据各信号频带的不同,设置具有相应频率特性的滤波器,如低通、高通、带通、带阻滤波器,使有用信号无衰减地通过,使干扰信号受到抑制。这类滤波器可用物理方法实现,此即模拟滤波器,也可用计算机通过算法实现,此即数字滤波器。对确定性信号的滤波处理也称常规滤波。

  随机信号没有确定的频谱,无法用常规滤波提取或抑制信号,但随机信号具有确定的功率谱,所以可根据有用信号和干扰信号的功率谱设计滤波器。维纳滤波是解决此类问题的方法之一。但设计维纳滤波器须作功率谱分解,只有当被处理信号为平稳的,干扰信号和有用信号均为一维,且功率谱为有理分式时,维纳滤波器的传递函数才可用伯特——香农设计法较容易地求解出,否则设计维纳滤波器存在着诸多困难。维纳滤波除设计思想与常规滤波不同外,对信号作抑制和选通这一点是相似的。

  卡尔曼滤波从与被提取信号有关的最测量中通过算法估计出所需信号。其中被估计信号是由白噪声激励引起的随机响应,激励源与响应之间的传递结构(系统方程)已知,量测量与被估计量之间的函数关系(量测方程)也已知。估计过程中利用了如下信息:系统方程、量测方程、白噪声激励的统计特性、量测误差的统计特性。由于所用信息都是时域内的量,所以卡尔曼滤波器是在时域内设计的,且适用于多维情况。这就完全避免了维纳滤波器在频域内设计遇到的限制和障碍,适用范围远比维纳滤波器广。

  从以上简述中可看出卡尔曼滤波有如下特点:

  (1)卡尔曼滤波处理的对象是随机信号;

  (2)被处理信号无有用和干扰之分,滤波的目的是要估计出所有被处理信号;

  (3)系统的白噪声激励和量测噪声并不是需要滤除的对象,它们的统计特性正是估计过程中需要利用的信息。

  所以确切地说,卡尔曼滤波应称作最优估计理论,此处称谓的滤波与常规滤波具有完全不同的概念和含意。

  就实现形式而言,卡尔曼滤波器实质上是一套由数字计算机实现的递推算法,每个递推周期中包含对被估计量的时间更新和量测更新两个过程。时间更新由上一步的量测更新结果和设计卡尔曼滤波器时的先验信息确定,量测更新则在时间更新的基础上根据实时获得的量测值确定。因此,量测量可看做卡尔曼滤波器的输入,估计值可看做输出。输入与输出之间由时间更新和量测更新算法联系,这与数字信号处理概念是类似的,所以有些书上称卡尔曼滤波为广义数字信号处理似有一定道理。

卡尔曼滤波理论的发展和工程应用

  随机信号没有既定的变化规律,对它们的估计也不可能完全准确,所谓最优估计也仅仅是指在某一准则下的最优。根据不同的最优准则,可获得随机信号的不同最优估计。使贝叶斯风险达到最小的估计为贝叶斯估计;使关于条件概率密度的似然函数达到极大的估计为极大似然估计;使验后概率密度达到极大的估计为极大验后估计;使估计误差的均方误差达到最小的估计为最小方差估计,若估计具有线性形式,则估计为线性最小方差估计,卡尔曼滤波即属此类估计

  卡尔曼滤波理论的创立是科学技术和社会需要发展到一定程度的必然结果。早在1795年,高斯(Karl Gauss)为测定行星运动轨道而提出了最小二乘估计法。20世纪40年代,为了解决火力控制系统精确跟踪问题,维纳(N. Weaner)于1942年提出了维纳滤波理论。维纳根据有用信号和干扰信号的功率谱确定出线性滤波器的频率特性,首次将数理统计理论与线性系统理论有机地联系在一起,形成了对随机信号作平滑、估计或预测的最优估计新理论。比维纳稍早,前苏联科学家戈尔莫克洛夫(A. N. Kolmogorov)于1941年也曾提出过类似理论。维纳给出了由功率谱求解维纳滤波器频率特性闭合解的一般方法,包括对功率谱的上、下平面分解及傅里叶变换和反变换,运算繁杂,解析求解十分困难。 1950年,伯特和香农给出了在功率谱为有理谱这一特殊条件下,由功率谱直接求取维纳滤波器传递函数的设计方法口, 这一方法简单易行,具有一定的工程实用价值。维纳滤波的最大缺点是适用范围极其有限,它 要求被处理信号必须是平稳的,且是一维的。人们试图将维纳滤波推广到非平稳和多维的情况,都因无法突破计算上的困难而难以推广和应用。

  采用频域设计法是造成维纳滤波器设计困难的根本原因。因此人们逐渐转向寻求在时域内直接设计最优滤波器的新方法。其中卡尔曼的研究最具有代表性,他提出的递推最优估计理论也因此而被称为卡尔曼滤波。由于采用了状态空间法描述系统,算法采用递推形式,所以卡尔曼滤波能处理多维和非平稳的随机过程。

  卡尔曼滤波理论一经提出,立即受到了工程界的重视,而工程应用中遇到的实际问题又使卡尔曼滤波的研究更深入更完善。1959年起美国太空署即NASA(National Aeronautics and Space Administration)开始研究载人太空船登月方案,当时提出了两个主要问题:(1)中途导航和制导改,(2)液体燃料助推器大挠度条件下的自动驾驶问题。因这两项研究的工作量都很庞大,无力同时进行,所以选择前者作为重点,即宇宙飞船的测轨间题。导航问题中主要解决对太空船运动状态的估计。量测信息来自三个子系统:飞船装备的惯性测量装置和天文观测仪,地面测轨系统,测轨数据经数据链传送至太空船。估计方法曾试图采用递推加权最小二乘和维纳滤波,均因精度满足不了要求和计算过于繁杂而不得不放弃。1960年秋,卡尔曼访问了NASA,提出了卡尔曼滤波算法,立即引起重视并投入研究。由于最初提出的卡尔曼滤波仅适用于线性系统,而实际系统是非线性系统,滤波初值应如何取才合理,这些都迫使卡尔曼做进一步的思考,广义卡尔曼滤波就是在此情况下提出来的。阿波罗计划中的导航系统后由麻省理工学院研制完成。卡尔曼滤波早期应用中的另一成功实例为C-5A飞机的多模式导航系统。

  卡尔曼滤波比维纳滤波的应用范围广,设计方法也简单易行得多,但它必须在计算机上执行,而 20 世纪 60 年代初,无论是速度、字长、容址等方面,计算机还处于低水平阶段。为了适应当时的技术水平,避免由于字长不够产生的舍入误差引起卡尔曼滤波的计算发散,Bierman,Carlson和Schmidt等人提出了平方根滤波算法和\(UDU^T\)分解滤波算法,以确保卡尔曼滤波增益回路中的滤波方差阵始终正定。

  卡尔曼最初提出的滤波基本理论只适用于线性系统,并且要求量测也必须是线性的。在之后的十多年时间内,Bucy,Sunahara等人致力于研究卡尔曼滤波理论在非线性系统和非线性最测情况下的推广,拓宽了卡尔曼滤波理论的适用范围。

  卡尔曼滤波最成功的工程应用是设计运载体的高精度组合导航系统。20 世纪 80 年代起,可供运载体装备的导航系统越来越多,非相似导航子系统的增加使量测信息增多,这对提高组合导航系统的精度十分有利。但是,如果采用集中式卡尔曼滤波器实现组合,则存在两个致命问题:(1)滤波器计算量以状态维数的三次方剧增,无法满足导航的实时性要求;(2)导航子系统的增加使故障率也随之增加,只要有一个子系统发生故障又没有及时检测出并隔离掉,则整个导航系统都会被污染。为了解决这一矛盾,1979年至1985年间,Speyer,Bierman和Kerr等人先后提出了分散滤波思想。并行计算技术的成熟为分散滤波的发展创造了有利条件。1988年起,Carlson提出了联邦滤波理论(Federated Filtering),旨在为容错组合导航系统提供设计理论。Carlson在装备运载体的诸多非相似导航子系统中选择导航信息全面、输出速率高、可靠性绝对保证的子系统作为公共参考系统,与其余子系统两两组合,形成若干个子滤波器。各子滤波器并行运行,获得建立在子滤波器局部量测基础上的局部最优估计,这些局部最优估计在第二级滤波器即主滤波器内按融合算法合成,获得建立在所有量测基础上的全局估计。全局估计再按信息守恒原则反馈给各子滤波器。实际设计的联邦滤波器是全局次优的,但是对于自主性要求特别高的重要运载体来说,导航系统的可靠性比精度更为重要。采用联邦滤波结构设计组合导航系统,虽然相对最优损失了少许精度,但换来的却是组合导航系统的高容错能力。目前美国空军已将联邦滤波器确定为新一代导航系统的通用滤波器。

Leave a Reply

Your email address will not be published. Required fields are marked *