卡尔曼滤波器综述
卡尔曼滤波器综述
瞿伟军
G10074
1、卡尔曼滤波的起源
1960年，匈牙利数学家卡尔曼发表了一篇关于离散数据线性滤波递推算法的论文，这意味着卡尔曼滤波的诞生。斯坦利.施密特(Stanley Schmidt)首次实现了卡尔曼滤波器，卡尔曼在NASA埃姆斯研究中心访问时，发现他的方法对于解决阿波罗计划的轨道预测很有用，后来阿波罗飞船的导航电脑使用了这种滤波器。关于这种滤波器的论文由Swerling (1958)、Kalman (1960)与 Kalman and Bucy (1961)发表。
2、卡尔曼滤波的发展
卡尔曼滤波是一种有着相当广泛应用的滤波方法，但它既需要假定系统是线性的，又需要认为系统中的各个噪声与状态变量均呈高斯分布，而这两条并不总是确切的假设限制了卡尔曼滤波器在现实生活中的应用。扩展卡尔曼滤波器（EKF）极大地拓宽了卡尔曼滤波的适用范围。EKF的基本思路是，假定卡尔曼滤滤对当前系统状态估计值非常接近于其真实值，于是将非线性函数在当前状态估计值处进行台劳展开并实现线性化。另一种非线性卡尔曼滤波叫线性化卡尔曼滤波。它与EKF的主要区别是前者将非线函数在滤波器对当前系统状态的最优估计值处线性化，而后者因为预先知道非线性系统的实际运行状态大致按照所要求、希望的轨迹变化，所以这些非线性化函数在实际状态处的值可以表达为在希望的轨迹处的台劳展开式，从而完成线性化。
不敏卡尔曼滤波器（UKF）是针对非线性系统的一种改进型卡尔曼滤波器。UKF处理非线性系统的基本思路在于不敏变换，而不敏变换从根本上讲是一种描述高斯随机变量在非线性化变换后的概率分布情况的方法。不敏卡尔曼滤波认为，与其将一个非线性化变换线性化、近似化，还不如将高斯随机变量经非线性变换后的概率分布情况用高斯分布来近似那样简单，因而不敏卡尔曼滤波算法没
有非线性化这一步骤。在每一定位历元，不敏卡尔曼滤波器按照一套公式产生一系列样点，每一样点均配有一个相应的权重，而这些带权的样点被用来完整地描述系统状态向量估计值的分布情况，它们替代了原先卡尔曼滤波器中的状态向量估计值及协方差。不敏卡尔曼滤器让这些样点一一经历非线性状态方程与测量方程，然后再将这些经非线性变换后的样点按照它们的权重而综合出对当前时刻的系统状态向量估计值。
多态自适应（MMA）卡尔曼滤波器是一种受到广泛关注的滤波器，它由好多个并联、同时运行的卡尔曼滤波器组成。在这组卡尔曼滤波器中，每一个滤波器对未知的滤波参数分别做出相互不同的假设，然后各自按照自己的模型假设进行滤波计算，而多态自适应滤波器最后将它们对系统状态的各个估计值进行加权，并以此作为最优估计值输出。
3、卡尔曼滤波器概述
简单来说，卡尔曼滤波器是一个“optimal recursive data processing algorithm（最优化自回归数据处理算法）”。对于解决很大部分的问题，他是最优，效率最高甚至是最有用的。他的广泛应用已经超过30年，包括机器人导航，控制，传感器数据融合甚至在军事方面的雷达系统以及导弹追踪等等。近年来更被应用于计算机图像处理，例如头脸识别，图像分割，图像边缘检测等等。
卡尔曼滤波器由一系列递归数学公式描述，它们提供了一种高效可计算的方法来估计过程的状态，并使估计均方误差最小。卡尔曼滤波器应用广泛且功能强大，它可以估计信号的过去和当前状态，甚至能估计将来的状态，即使并不知道模型的确切性质。
假设我们要研究的对象是一个房间的温度。根据你的经验判断，这个房间的温度是恒定的，也就是下一分钟的温度等于现在这一分钟的温度（假设我们用一分钟来做时间单位）。假设你对你的经验不是100%的相信，可能会有上下偏差几度。我们把这些偏差看成是高斯白噪声（White Gaussian Noise），也就是这些偏差跟前后时间是没有关系的而且符合高斯分配（Gaussian Distribution）。另外，我们在房间里放一个温度计，但是这个温度计也不准确的，测量值会比实际值偏差。我们也把这些偏差看成是高斯白噪声。
现在对于某一分钟我们有两个有关于该房间的温度值：你根据经验的预测值（系统的预测值）和温度计的值（测量值）。下面我们要用这两个值结合他们各自的噪声来估算出房间的实际温度值。
假如我们要估算k 时刻的是实际温度值。首先你要根据1k -时刻的温度值，来预测k 时刻的温度。因为你相信温度是恒定的，所以你会得到k 时刻的温度预测值是跟1k -时刻一样的，假设是23度，同时该值的高斯噪声的偏差是5度（5是这样得到的：如果1k -时刻估算出的最优温度值的偏差是3，你对自己预测的不确定度是4度，他们平方相加再开方，就是5）。然后，你从温度计那里得到了
k 时刻的温度值，假设是25度，同时该值的偏差是4度。
由于我们用于估算k 时刻的实际温度有两个温度值，分别是23度和25度。究竟实际温度是多少呢？相信自己还是相信温度计呢？究竟相信谁多一点，我们
可以用他们的协方差来判断。因为
252(5242)Kg ∧∧∧∧
=+，所以0.78Kg =，我们可以估算出k 时刻的实际温度值是：230.78(2523)24.56+*-=度。可以看出，因为温度计的协方差比较小（比较相信温度计），所以估算出的最优温度值偏向温度计的值。
现在我们已经得到k 时刻的最优温度值了，下一步就是要进入1k -时刻，进行新的最优估算。到现在为止，好像还没看到什么自回归的东西出现。对了，在进入1k -时刻之前，我们还要算出k 时刻那个最优值（24.56度）的偏差。算法
如下：
((1)52)0.5 2.35Kg ∧∧
-*=。这里的5就是上面的k 时刻你预测的那个23度温度值的偏差，得出的2.35就是进入1k +时刻以后k 时刻估算出的最优温度值的偏差（对应于上面的3）。
就是这样，卡尔曼滤波器就不断的把协方差递归，从而估算出最优的温度值。他运行的很快，而且它只保留了上一时刻的协方差。上面的Kg ，就是卡尔曼增益（Kalman Gain ）。他可以随不同的时刻而改变他自己的值。
4、卡尔曼滤波器算法
首先，我们先要引入一个离散控制过程的系统。该系统可用一个线性随机微分方程（Linear Stochastic Difference equation ）来描述：
()(1)()()X k AX k BU k W K =-++ (3.1)
再加上系统的测量值：
()()()Z k HX k V k =+ (3.2)
上两式子中，()X k 是k 时刻的系统状态，
()U k 是k 时刻对系统的控制量。A 和B 是系统参数，对于多模型系统，他们为矩阵。()Z k 是k 时刻的测量值，H 是测量系统的参数，对于多测量系统，H 为矩阵。()W K 和()V k 分别表示过程和测量的噪声。他们被假设成高斯白噪声(White Gaussian Noise)，他们的协方差分别是,Q R （这里我们假设他们不随系统状态变化而变化）。
对于满足上面的条件(线性随机微分系统，过程和测量都是高斯白噪声)，卡尔曼滤波器是最优的信息处理器。下面我们来用他们结合他们的协方差s 来估算系统的最优化输出。
首先我们要利用系统的过程模型，来预测下一状态的系统。假设现在的系统状态是k ，根据系统的模型，可以基于系统的上一状态而预测出现在状态：
(1)(11)()
X k k AX k k BU k -=--+ (3.3)
式(3.3)中，
(1)
X k k -是利用上一状态预测的结果，
(11)
X k k --是上一状
态最优的结果，()U k 为现在状态的控制量，如果没有控制量，它可以为0。
到现在为止，我们的系统结果已经更新了，可是，对应于(1)
X k k -的协方
差(协方差)还没更新。我们用P 表示协方差：
'(1)(11)P k k AP k k A Q
-=--+
(3.4)
式(3.4)中，
(1)
P k k -是
(1)
X k k -对应的协方差，
(11)
P k k --是
(11)
X k k --对应的协方差，'
A 表示A 的转置矩阵，Q 是系统过程的协方差。式
子(3.3),(3.4)就是卡尔曼滤波器5个公式当中的前两个，也就是对系统的预测。
现在我们有了现在状态的预测结果，然后我们再收集现在状态的测量值。结合预测值和测量值，我们可以得到现在状态()k 的最优化估算值
()
X k k ：
()(1)()(()(1))
X k k X k k Kg k Z k HX k k =-+--
(3.5)
其中Kg 为卡尔曼增益(Kalman Gain)： ''()(1)((1))
Kg k P k k H HP k k H R =--+
(3.6)
到现在为止，我们已经得到了k 状态下最优的估算值()
X k k 。但是为了要令
卡尔曼滤波
器不断的运行下去直到系统过程结束，我们还要更新k 状态下()
X k k 的协方
差：
()(())(1)
P k k I Kg k H P k k =-- (3.7)
其中I 为1的矩阵，对于单模型单测量，1I =。当系统进入1k +状态时，
()
P k k 就是式子(3.4)的
(11)
P k k --。这样，算法就可以自回归的运算下去。
式子（3.3），（3.4），（3.5），（3.6）和（3.7）就是卡尔曼滤波器的基本原理了。根据这5个公式，可以很容易的实现计算机的程序。
5、卡尔曼滤波器的原理
在现代，随机最优控制和随机信号处理技术中，信号和噪声往往是多维非平稳随机过程。卡尔曼滤波理论采用时域上的递推算法在数字计算机上进行数据滤波处理。
对于离散域线性系统：
()(1)(()())x k Ax k B u k w k =-++ （3.8）
()()()v y k Cx k v k =+ （3.9）
式中，()w k 为过程噪声信号， ()v k 为测量噪声信号。 离散卡尔曼滤波器递推算法为：
()()()T
n T
P k C M k CP k C R
=+ （3.10） ()(1)T T P k AP k A BQB =-+ （3.11） ()(())()n n P k I M k C P k =- （3.12） ()(1)()(()(1))n v x k Ax k M k y k CAx k =-+-- （3.13） ()()e y k Cx k = （3.14）
误差的协方差为：
cov()()T err k CP k C = （3.15）
卡尔曼滤波器结构如图所示
卡尔曼滤波器结构
6、卡尔曼滤波器的应用实例
6.1、基于VDLL 的GPS 信号跟踪算法
图2 VDLL 的基本结构
系统的状态量为用户状态( 位置、速度、用户钟差、用户钟漂移) ，观测量为所有卫星通道的码环鉴相器输出组成的矢量D(X)
X =(x y z b x ? y ? z ? b ?
)T
D(X)=(D(Δτ1) D(Δτ2) … D(ΔτN ))T （3-1）
其中x ，y ，z ，b ，x ?
，y ?
，z ?
，b ?
分别为在ECEF( Earth Centered Earth Fixed, 地心坐标系) 下的用户接收机三维位置和速度， b , b ?
分别为接收机钟差和钟差变化率。预检测积分时间T 选择为1ms ，即滤波器迭代的k 时刻与k+1时刻间隔1ms 。假设用户为匀加速运动，可写出状态方程如下
1k 1k 44
44k w X I 0I I X --+???
???=T （3-2） 其中w k-1 为状态扰动噪声。 导航滤波器把k- 1 时刻的状态和对应的卫星位置信息反馈给本地信号生成模块，产生超前、即时与滞后码序列，分别与k
时刻的1ms 卫星信号的同相支路与正交支路做相关运算，即图2 中的“积分存储”模块。为了去除对信号幅度变化的敏感性和对载波环路锁定的依赖性，码环鉴相器选择归一化的超前码与滞后码相关能量差( 第i 颗卫星)
2
222
2
222,)(QL IL QE IE QL IL QE IE D k i +++--+=?τ （3-3） 该码环鉴相器在没有噪声的理想情况下，当超前码与滞后码间距为 1 chip( 码片) 时，在±0.5 chip 输入误差范围内，鉴相器输出等于真实的跟踪误差。在实际应用中，超前滞后码的移位通常以采样点为单位，1 个chip 的长度不一定为整数个采样点，但鉴相器仍然与输入误差成线性关系，所以式可写为
k i k i i k i v g D ,,,)(+?=?ττ （3-4）
其中g i 为待求的鉴相器特性在线性区的斜率， v i, k 为观测噪声。鉴相器输入误差k i ,τ?是k 时刻( 接收机时刻) 第i 颗卫星1ms 长度的采样信号与本地即时码之间的相位偏移，而本地即时码的初始相位值等于k- 1时刻对卫星信号的相位估计。所以Δτi, k 实质上是k 时刻卫星信号的相位与k - 1时刻的相位估计之差，等于k 时刻信号传输延迟的真实值与k- 1 时刻信号传输延迟的估计值。 k 时刻接收到第i 颗卫星信号的理论测量延迟为
k k i k k i k k i k k i b z z y y x x c
+-+-+-=
2,2,2,,)()()(1
τ （3-5） 其中( x k , y k , z k ) 为k 时刻用户的三维位置坐标,( x i, k , y i, k , z i, k ) 为k 时刻第i 颗卫星的位置坐标,b k 为k 时刻的接收机钟差。K-1 时刻信号传输延迟的估计值可由k-1 时刻的用户状态估计X ^
( k- 1 |k- 1) 和卫星位置计算得到
1|121,1|121,1|121,1|11,?)?()?()?(1?------------+-+-+-=
k k k i k k k i k k k i k k k i b x z x y x x c
τ （3-6） 由( 4) ) ( 6) 式，可写出第i 颗卫星的观测方程如下
k i k i k k k i k k k i k k k k k k i k k i k k i k i k v x z x y x x c
b b z z y y x x c
g X D ,21,1|121,1|121,1|11|12,2,2,})?()?()?(1
?)()()(1{
)(+-+-+---+-+-+-=----------- （3-7）
其他卫星的观测方程可以类似地写出。
6.2 卡尔曼滤波器的仿真及分析 6.2.1 系统描述
验证卡尔曼滤波器的滤波性能 对象为二阶传递函数：
2
133
()25P G s s s
=
+ （3.16） 取采样时间为1ms ，采样Z 变换将对象离散化，并描述离散状态方程的形式： (1)()(()())x k Ax k B u k w k +=++ （3.17）
()()y k Cx k = （3.18）
带有测量噪声的被控对象输出为：
()()()v y k Cx k v k =+ （3.19） 式中，
[][]1.0000000,0.00098760.0000659,,1,0,00.0000000,0.97530990.1313512A B C D ????====????
????
6.2.2 仿真方法一：采用M 语言进行仿真
控制干扰信号()w k 和测量噪声信号()v k 幅值均为0.10的白噪声信号，输入信号幅值为1.0、频率为1.5Hz 的正图信号。采用卡尔曼滤波器实现信号的滤波，取1,1Q R ==。仿真时间为3s ，误差协方差的变化，原始信号及带有噪声的原始信号和原始信号及滤波后的信号和分别如图1、图2和图3所示。
图1误差协方差
的变化
-3
C o v a r i a n c e o f e s t i m a t i o n e r r o r
00.51
1.52
2.53
time(s)
y -i d e a l s i g n a l ; y v -s i g n a l w i t h n o i s e
图2原始信号及带有噪声的原始信号
00.51
1.52
2.53
-0.2
0.2
0.4
0.6
0.8
1
1.2
time(s)
y -i d e a l s i g n a l ; y e -f i l t e r e d s i g n a l
图3原始信号及滤波后的信号
仿真结果表明，该滤波器对控制干扰和测量噪声具有良好的滤波作用。
6.2.3 仿真方法二 ：采用Simulink 进行仿真
卡尔曼算法由M函数实现。控制干扰信号()
w k和测量噪声信号()
v k幅值均为0.10的白噪声信号，输入信号幅值为1.0、频率为0.5Hz的正弦信号采用卡尔曼滤波器实现信号的滤波，取1,1
==。仿真时间为10s，仿真结果如下图
Q R
4和图5所示
y（黄色线）
图4 原始信号y（红色线）及滤波后的信号
e
y
图5原始信号y及带有噪声的原始信号
v
仿真程序如下图所示
基于Kalman 滤波器的Simulink 仿真
6.3 基于卡尔曼滤波器的 PID 控制算法及仿真 6.3.1 基于卡尔曼滤波器 PID 控制的原理
结合卡尔曼滤波器与经典PID 控制得到一种基于卡尔曼滤波器的PID 控制方法。其控制系统结构如图：
基于卡尔曼滤波器的PID 控制系统结构
6.3.2 仿真程序及分析
被控对象为二阶传递函数[8]：
2133
()25P G s s s
=
+ (5.1)
离散化结果与
3.4的仿真实例相同。采样时间为1ms 。
控制干扰信号()w k 和测量噪声信号()v k 幅值均为0.002的白噪声信号，输
入信号为一阶跃信号。采用卡尔曼滤波器实现信号的滤波，取1,1Q R ==。仿真时间为1s ，分两种情况进行仿真：1M =时为未加滤波，2M =时为加滤波。在PID 控制器中，取8.0,0.80,0.20p i d k k k ===。用MATLAB 进行仿真[9]，加入滤波前后PID 阶跃响应如图6和图7所示。
图6未加入滤波的PID 控制阶跃响应（1M =）
0.1
0.2
0.3
0.4
0.50.6
0.7
0.8
0.9
1
00.20.4
0.6
0.8
1
1.2
1.4
time(s)
r i n ,y o u t
图7 加入滤波后PID 控制阶跃响应（2M ）
仿真结果表明，采用滤波器使控制效果明显改善。
00.10.20.30.4
0.50.60.70.80.91
0.2
0.4
0.6
0.8
1
1.2
1.4
time(s)
r i n ,y o u t
卡尔曼滤波器及其简matlab仿真
卡尔曼滤波器及其简matlab仿真
卡尔曼滤波器及其简matlab仿真 一、卡尔曼滤波的起源 谈到信号的分析与处理，就离不开滤波两个字。通常，信号的频谱处于有限的频率范围内，而噪声的频谱则散布在很广的频率范围内，为了消除噪声，可以进行频域滤波。但在许多应用场合，需要直接进行时域滤波，从带噪声的信号中提取有用信号。虽然这样的过程其实也算是对信号的滤波，但其所依据的理论，即针对随机信号的估计理论，是自成体系的。人们对于随机信号干扰下的有用信号不能“确知”，只能“估计”。为了“估计”，要事先确定某种准则以评定估计的好坏程度。 1960年卡尔曼发表了用递归方法解决离散数据线性滤波问题的论文A New Approach to Linear Filtering and Prediction Problems （线性滤波与预测问题的新方法），在这篇文章里一种克服了维纳滤波缺点的新方法被提出来，这就是我们今天称之为卡尔曼滤波的方法。卡尔曼滤波应用广泛且功能强大，它可以估计信号的过去和当前状态甚至能估计将来的状态即使并不知道模型的确切性质。 其基本思想是以最小均方误差为最佳估计准则，采用信号与噪声的状态空间模型利用前一时刻的估计值和当前时刻的观测值来更新对状态变量的估计，求出当前时刻的估计值。算法根据建立的系统方程和观测方程对需要处理的信号做出满足最小均方误差的估计。 对于解决很大部分的问题，它是最优，效率最高甚至是最有用的。它的广泛应用已经超过30年，包括机器人导航，控制，传感器数据融合甚至在军事方面的雷达系统以及导弹追踪等等。近年来更被应用于计算机图像处理，例如头脸识别，图像分割，图像边缘检测等等。 卡尔曼滤波不要求保存过去的测量数据，当新的数据到来时，根据新的数据和前一时刻的储值的估计，借助于系统本身的状态转移方程，按照一套递推公式，即可算出新的估值。卡尔曼递推算法大大减少了滤波装置的存储量和计算量，并且突破了平稳随机过程的限制，使卡尔曼滤波器适用于对时变信号的实时处理。
基于FPGA的卡尔曼滤波器的设计
基于FPGA的卡尔曼滤波器的设计 时间：2010-04-12 12:52:33 来源：电子科技作者：米月琴，黄军荣西安电子科技大学摘要：针对电路设计中经常碰到数据的噪声干扰现象，提出了一种Kalman滤波的FPGA实现方法。该方法采用了TI公司的高精度模数转换器ADSl25l以及Altera公司的EPlCl2，首先用卡尔曼滤波算法 设计了一个滤波器，然后将该滤波器分解成简单的加、减、乘、除运算。通过基于FPGA平台的硬件与 软件的合理设计，成功地实现了数据噪声的滤除设计，并通过实践仿真计算，验证了所实现滤波的有效性。 关键词：卡尔曼；FPGA；最小方差估计 卡尔曼滤波是一个“Optimal Recursive Data Processing Algorithm(最优化自回归数据处 理算法)”，对于解决很大部分的问题，是最优化的，效率最高甚至是最有用的。传统的卡尔曼滤波是 在DSP上实现的。但是DSP成本相对较高，而且指令是串行执行的，不能满足有些要求较高的场合。而FPGA由于其硬件结构决定了它的并行处理方式，无论在速度还是实时性都更胜一筹。文中以基于FPGA 器件和A／D转换器的数据采集系统为硬件平台，进行了卡尔曼滤波算法设计，详述了基于FPGA的卡尔 曼滤波器的设计实现。 1 卡尔曼滤波算法 工程中，为了了解工程对象(滤波中称为系统)的各个物理量(滤波中称为状态)的确切数值，或为了 达到对工程对象进行控制的目的，必须利用测量手段对系统的各个状态进行测量。但是，量测值可能仅 是系统的部分状态或是部分状态的线性组合，且量测值中有随机误差(常称为量测噪声)。最优估计就是 针对上述问题的一种解决方法。它能将仅与部分状态有关的测量进行处理，得出从统计意义上讲误差最 小的更多状态的估值。误差最小的标准常称为估计准则，根据不同的估计准则和估计计算方法，有各种 不同的最优估计，卡尔曼滤波是一种递推线性最小方差估计的最优估计。 系统的状态方程可设定为 式（3）为系统噪声。设设备的量测噪声为Vk，系统得量测方程为
三阶卡尔曼滤波数字锁频环设计及性能分析
三阶卡尔曼滤波数字锁频环设计及性能分析 作者：李金海， 巴晓辉， 陈杰， LI Jin-hai， BA Xiao-hui， Chen Jie 作者单位：中国科学院微电子研究所,北京,朝阳区,100029 刊名： 电子科技大学学报 英文刊名：JOURNAL OF UNIVERSITY OF ELECTRONIC SCIENCE AND TECHNOLOGY OF CHINA 年，卷(期)：2008，37(5) 被引用次数：1次 参考文献(10条) 1.HINEDI S.STATMAN J I High-dynamic GPS tracking final report[JPL Publication 88-35] 1988 2.AGUIRRE S.HINEDI S Two novel automatic frequency tracking loops 1989(05) 3.HINEDI S An extended Kalmaa filter based automatic frequency control loop[TDA Progress Report 42-95] 1988 4.VILNROTTER V A.HINEDI S.KUMAR R Frequency estimation techniques for high dynamic trajectories 1989(04) 5.BAR-SHALOM Y.LI X R.KIRUBARAJAN T Estimation with applications to tracking and navigation:Theory algorithms and soRware 2001 6.张厥盛.郑继禹.万心平锁相技术 2005 7.JURY E I Theory and application of the z-transform method 1964 8.邓自立.郭一新现代时间序列分析及其应用--建模、滤波、去卷、预报和控制 1988 9.GREWAL M S.ANDREWS A P Kalman filtering:theory and practice using matlab 2001 10.ARNOLD W https://www.wendangku.net/doc/6c2113153.html,UB A J Generalized eigenproblem algorithms and software for algebraic Riccati equations 1984(12) 相似文献(9条) 1.学位论文孙峰高动态多星座接收机捕获和跟踪技术的研究与实现2009 全球导航卫星系统(GNSS)是用于定位用户接收机地理位置的一种卫星系统。目前，GNSS包括现已投入运行的三个卫星定位系统：全球定位系统(GPS)、全球导航卫星系统( GLONASS)、北斗一代系统(BD)。鉴于多星座导航定位系统的建立，多星座接收机将大大提高卫星导航定位的可靠性、精度和实时性。 高动态接收机的捕获和跟踪技术一直是研究的热点和难点。许多学者针对高动态的特殊应用做出了一些卓有成效的研究，提出了多种设计方案，重点为伪码的快速捕获和多普勒频率的跟踪。伪码的快速捕获的主要方法为：基于FFT和匹配滤波的并行捕获方法以及串并结合的滑动相关捕获方法。这些捕获方法在捕获性能和复杂性上各有优劣。本文采用了串并结合的滑动相关捕获方法，这种算法的捕获性能较好，硬件实现简单。 载波多普勒频率跟踪的主流方案是采用锁频环(FLL)+锁相环(PLL)的环路跟踪结构。使用FLL来跟踪频率的快速变化，当频率引导到PLL可处理的范围时，通过PLL来跟踪相位的变化，精确的锁定载波频率。本文采用二阶锁频环辅助三阶锁相环的环路结构，可很好的跟踪接收机的动态。 本文的主要内容为： 1.完成了多星座卫星信号接收机的硬件设计，为系统的实现搭建了硬件平台。 2.在分析了GPS、GLONASS、BD的伪码特性的基础上，采用串并结合的时域相关捕获的方法，缩短了伪码的捕获时间。 3.研究并设计了DLL码跟踪环路。经过测试验证了设计的DLL环路的正确性。 4.载波跟踪环采用三阶二象限反正切Costas环和二阶四相锁环相结合的方法，有效的消除了高动态的影响。 本论文设计的捕获和跟踪的方法最终在高动态多星座接收机上得到了实现，测试结果表明本文的设计满足系统指标要求。 2.学位论文郑宏磊GPS在干扰环境下的可用性研究2006 全球定位系统(GPS)能在全球范围内提供精确的位置、速度和时间信息，在军事和民用领域发挥着极其重要的作用。随着GPS的广泛应用，它易受到干扰的弱点也随之暴露出来，针对GPS进行的抗干扰技术也日益成为研究的热点。本文阐述了全球定位系统的工作原理，系统组成以及信号格式，在此基础上着重分析了GPS受干扰特性，为以后的工作奠定了基础。 本文将理论分析和实验相结合，结合商用GPS接收机的实际测量结果，对GPS信号受干扰前后的特性进行了分析。针对射频RF等干扰源以及多路径 ，本文介绍了抗干扰的总体设计方案，分析了几种可行的抗干扰措施，重点对环路滤波和自适应调零天线进行了研究设计。 论文在环路滤波器设计方面采用了由锁频环(FLL)辅助的锁相环(PLL)滤波器，在自适应调零天线方面设计空间-时间自适应阵列以代替空间自适应阵列，并采用功率最小预处理算法。最后通过实验仿真得到了较为理想的结果，可在一定程度上保证GPS在干扰环境下的可用性。 3.期刊论文李国栋.崔晓伟.尹旭明.冯振明.LI Guodong.CUI Xiaowei.YIN Xuming.FENG Zhenming GPS接收机中锁 频环频率误锁的检测-清华大学学报（自然科学版）2007,47(1) 为了解决全球定位系统(GPS)接收机中的锁频环在载波同步过程中可能出现的频率误锁问题,在分析了锁频环在噪声环境下的工作原理及产生频率误锁原因的基础上,基于有无发生频率误锁时同一信息符号对应的多个预检测积分值的变化规律,提出了一种用于频率误锁检测和快速纠正的算法.仿真结果表明:该方法能够在锁频环完成工作之后及时判决是否有误锁发生,误锁时可在1～2个导航比特时间内把载波频率调整到正确频率上.该方法实现简单,可
卡尔曼滤波器及其简matlab仿真.
卡尔曼滤波器及其简matlab仿真 一、卡尔曼滤波的起源 谈到信号的分析与处理，就离不开滤波两个字。通常，信号的频谱处于有限的频率范围内，而噪声的频谱则散布在很广的频率范围内，为了消除噪声，可以进行频域滤波。但在许多应用场合，需要直接进行时域滤波，从带噪声的信号中提取有用信号。虽然这样的过程其实也算是对信号的滤波，但其所依据的理论，即针对随机信号的估计理论，是自成体系的。人们对于随机信号干扰下的有用信号不能“确知”，只能“估计”。为了“估计”，要事先确定某种准则以评定估计的好坏程度。 1960年卡尔曼发表了用递归方法解决离散数据线性滤波问题的论文A New Approach to Linear Filtering and Prediction Problems（线性滤波与预测问题的新方法），在这篇文章里一种克服了维纳滤波缺点的新方法被提出来，这就是我们今天称之为卡尔曼滤波的方法。卡尔曼滤波应用广泛且功能强大，它可以估计信号的过去和当前状态甚至能估计将来的状态即使并不知道模型的确切性质。 其基本思想是以最小均方误差为最佳估计准则，采用信号与噪声的状态空间模型利用前一时刻的估计值和当前时刻的观测值来更新对状态变量的估计，求出当前时刻的估计值。算法根据建立的系统方程和观测方程对需要处理的信号做出满足最小均方误差的估计。 对于解决很大部分的问题，它是最优，效率最高甚至是最有用的。它的广泛应用已经超过30年，包括机器人导航，控制，传感器数据融合甚至在军事方面的雷达系统以及导弹追踪等等。近年来更被应用于计算机图像处理，例如头脸识别，图像分割，图像边缘检测等等。 卡尔曼滤波不要求保存过去的测量数据，当新的数据到来时，根据新的数据和前一时刻的储值的估计，借助于系统本身的状态转移方程，按照一套递推公式，即可算出新的估值。卡尔曼递推算法大大减少了滤波装置的存储量和计算量，并且突破了平稳随机过程的限制，使卡尔曼滤波器适用于对时变信号的实时处理。 二、卡尔曼滤波的原理
几种卡尔曼滤波算法理论
自适应卡尔曼滤波 卡尔曼滤波发散的原因 如果卡尔曼滤波是稳定的，随着滤波的推进，卡尔曼滤波估计的精度应该越来越高，滤波误差方差阵也应趋于稳定值或有界值。但在实际应用中，随着量测值数目的增加，由于估计误差的均值和估计误差协方差可能越来越大，使滤波逐渐失去准确估计的作用，这种现象称为卡尔曼滤波发散。 引起滤波器发散的主要原因有两点： （1）描述系统动力学特性的数学模型和噪声估计模型不准确，不能直接真实地反映物理过程，使得模型与获得的量测值不匹配而导致滤波发散。这种由于模型建立过于粗糙或失真所引起的发散称为滤波发散。 （2）由于卡尔曼滤波是递推过程，随着滤波步数的增加，舍入误差将逐渐积累。如果计算机字长不够长，这种积累误差很有可能使估计误差方差阵失去非负定性甚至失去对称性，使滤波增益矩阵逐渐失去合适的加权作用而导致发散。这种由于计算舍入误差所引起的发散称为计算发散。 针对上述卡尔曼滤波发散的原因，目前已经出现了几种有效抑制滤波发散的方法，常用的有衰减记忆滤波、限定记忆滤波、扩充状态滤波、有限下界滤波、平方根滤波、和自适应滤波等。这些方法本质上都是以牺牲滤波器的最优性为代价来抑制滤波发散，也就是说，多数都是次优滤波方法。 自适应滤波 在很多实际系统中，系统过程噪声方差矩阵Q和量测误差方差阵R事先是不知道的，有时甚至连状态转移矩阵 或量测矩阵H也不能确切建立。如果所建立的模型与实际模型不符可能回引起滤波发散。自适应滤波就是这样一种具有抑制滤波发散作用的滤波方法。在滤波过程中，自适应滤波一方面利用量测值修正预测值，同时也对未知的或不确切的系统模型参数和噪声统计参数进行估计修正。自适应滤波的方法很多，包括贝叶斯法、极大似然法、相关法与协方差匹配法，其中最基本也是最重要的是相关法，而相关法可分为输出相关法和新息相关法。 在这里只讨论系统模型参数已知，而噪声统计参数Q和R未知情况下的自适应滤波。由于Q和R等参数最终是通过增益矩阵K影响滤波值的，因此进行自适应滤波时，也可以不去估计Q和R等参数而直接根据量测数据调整K就可以了。
维纳最速下降法滤波器卡尔曼滤波器设计及Matlab仿真
信息融合大作业 ——维纳最速下降法滤波器，卡尔曼滤波器设计及Matlab仿真 1.滤波问题浅谈 估计器或滤波器这一术语通常用来称呼一个系统，设计这样的系统是为了从含有噪声的数据中提取人们感兴趣的，接近规定质量的信息。由于这样一个宽目标，估计理论应用于诸如通信、雷达、声纳、导航、地震学、生物医学工程、 金融工程等众多不同的领域。例如，考虑一个数字通信系统，其基本形式由发
射机、信道和接收机连接组成。发射机的作用是把数字源(例如计算机)产生的0、1符号序列组成的消息信号变换成为适合于信道上传送的波形。而由于符号间干扰和噪声的存在，信道输出端收到的信号是含有噪声的或失真的发送信号。接收机的作用是，操作接收信号并把原消息信号的一个可靠估值传递给系统输出端的某个用户。随着通信系统复杂度的提高，对原消息信号的还原成为通信系统中最为重要的环节，而噪声是接收端需要排除的最主要的干扰，人们也设计出了针对各种不同条件应用的滤波器，其中最速下降算法是一种古老的最优化技术，而卡尔曼滤波器随着应用条件的精简成为了普适性的高效滤波器。2．维纳最速下降算法滤波器 2.1 最速下降算法的基本思想 考虑一个代价函数，它是某个未知向量的连续可微分函数。函数 将的元素映射为实数。这里，我们要寻找一个最优解。使它满足如下条件 (2.1) 这也是无约束最优化的数学表示。 特别适合于自适应滤波的一类无约束最优化算法基于局部迭代下降的算法： 从某一初始猜想出发，产生一系列权向量，使得代价函数在算法的每一次迭代都是下降的，即 其中是权向量的过去值，而是其更新值。 我们希望算法最终收敛到最优值。迭代下降的一种简单形式是最速下降法，该方法是沿最速下降方向连续调整权向量。为方便起见，我们将梯度向量表示为
卡尔曼滤波器介绍 --- 最容易理解
10.6 卡尔曼滤波器简介 本节讨论如何从带噪声的测量数据把有用信号提取出来的问题。通常，信号的频谱处于有限的频率范围内，而噪声的频谱则散布在很广的频率范围内。如前所述，为了消除噪声，可以把 FIR滤波器或IIR滤波器设计成合适的频带滤波器，进行频域滤波。但在许多应用场合，需要进行时域滤波，从带噪声的信号中提取有用信号。虽然这样的过程其实也算是对信号的滤波，但所依据的理论,即针对随机信号的估计理论,是自成体系的。人们对随机信号干扰下的有用信号不能“确知”，只能“估计”。为了“估计”，要事先确定某种准则以评定估计的好坏程度。最小均方误差是一种常用的比较简单的经典准则。典型的线性估计器是离散时间维纳滤波器与卡尔曼滤波器。 对于平稳时间序列的最小均方误差估计的第一个明确解是维纳在1942年2月首先给出的。当时美国的一个战争研究团体发表了一个秘密文件，其中就包括维纳关于滤波问题的研究工作。这项研究是用于防空火力控制系统的。维纳滤波器是基于最小均方误差准则的估计器。为了寻求维纳滤波器的冲激响应，需要求解著名的维纳-霍夫方程。这种滤波理论所追求的是使均方误差最小的系统最佳冲激响应的明确表达式。这与卡尔曼滤波（Kalman filtering）是很不相同的。卡尔曼滤波所追求的则是使均方误差最小的递推算法。 在维纳进行滤波理论研究并导出维纳-霍夫方程的十年以前，在1931年，维纳和霍夫在数学上就已经得到了这个方程的解。 对于维纳-霍夫方程的研究，20世纪五十年代涌现了大量文章，特别是将维纳滤波推广到非平稳过程的文章甚多，但实用结果却很少。这时正处于卡尔曼滤波问世的前夜。 维纳滤波的困难问题，首先在上世纪五十年代中期确定卫星轨道的问题上遇到了。1958年斯韦尔林（Swerling）首先提出了处理这个问题的递推算法，并且立刻被承认和应用。1960年卡尔曼进行了比斯韦尔林更有意义的工作。他严格地把状态变量的概念引入到最小均方误差估计中来，建立了卡尔曼滤波理论。空间时代的到来推动了这种滤波理论的发展。 维纳滤波与卡尔曼滤波所研究的都是基于最小均方误差准则的估计问题。 维纳滤波理论的不足之处是明显的。在运用的过程中，它必须把用到的全部数据存储起来，而且每一时刻都要通过对这些数据的运算才能得到所需要的各种量的估值。按照这种滤波方法设置的专用计算机的存储量与计算量必然很大，很难进行实时处理。虽经许多科技工作者的努力，在解决非平稳过程的滤波问题时，给出能用的方法为数甚少。到五十年代中期，随着空间技术的发展，这种方法越来越不能满足实际应用的需要，面临了新的挑战。尽管如此，维纳滤波理论在滤波理论中的开拓工作是不容置疑的，维纳在方法论上的创见，仍然影响着后人。 五十年代中期，空间技术飞速发展，要求对卫星轨道进行精确的测量。为此，人们将滤波问题以微分方程表示，提出了一系列适应空间技术应用的精练算法。1960年
扩展卡尔曼滤波器(EKF)进行信号处理及信号参数估计
% 扩展卡尔曼滤波器估计单相电压幅值、相位、频率参数（含直流）function test2_EKF close all; clc; tic; %计时 %模型：y=A0+A1*cos(omega*t+phy1) %离散化：y(k)=A0(k)+A1(k)*cos(omega(k)*k*Ts+phy1(k)) %状态变量：x1(k)=A0(k),x2(k)=omega(k),x3(k)=A1(k)*cos(omega(k)*k*Ts+phy1(k) ),x4(k)=A1(k)*sin(omega(k)*k*Ts+phy1(k)) %下一时刻状态变量为(假设状态不突变)：A0(k+1)=A0(k),A1(k+1)=A1(k),omega(k+1)=omega(k),phy1(k+1)=phy1 (k); %则对应状态为：x1(k+1)=x1(k),x2(k+1)=x2(k),x3(k+1)=x3(k)*cos(x2(k)*Ts)- x4(k)*sin(x(2)*Ts),x4(k+1)=x3(k)*sin(x2(k)*Ts)+x4(k)*cos(x(2)*Ts); %状态空间描述：X(k+1)=f(X(k))+W(k);y(k)=H*X(k)+v(k) %f(X(k))=[x1(k);x2(k);x3(k)*cos(x2(k)*Ts)- x4(k)*sin(x(2)*Ts);x3(k)*sin(x2(k)*Ts)+x4(k)*cos(x(2)*Ts)] %偏导(只求了三个)：f`(X(k))=[1,0,0;0,1,0;0,-x3(k)*Ts*sin(x2(k)*Ts)-x4(k)*Ts*cos(x2(k)*Ts),cos(x2(k)*Ts);0,x3(k)*Ts*cos(x2(k)*Ts)- x4(k)*Ts*sin(x2(k)*Ts),sin(x2(k)*Ts)]
直流电机运行状态的卡尔曼滤波估计器设计.doc
二 〇 一 五 年 六 月 题 目：直流电机运行状态的卡尔曼滤波估计器设计 学生姓名：张傲 学 院：电力学院 系 别：电力系 专 业：风能与动力工程 班 级：风能11-1 指导教师：董朝轶 教授
摘要 卡尔曼滤波是一个迭代自回归算法，对于连续运动状态用中的大部分问题它都能够给出最优的预测。它已经广泛应用了近半个世纪，例如数据的融合，机械的导航乃至军用雷达的导航等等。卡尔曼滤波一般用于动态数据的处理，是从混沌的信号中提取有用信号消除误差的参数估计法。卡尔曼滤波是依据上一个估计数值和当下的检测数据运用递推估计算出当前的估计值。通过状态方程运用递推的方法进行估计,可以建立物体运动的模型。本文采用的工程设计对运行状态下的直流电机进行参数的计算和校验。而且直流电机的调节性能非常好只需要加上电阻调压就可以了，而且启动曲线非常好，启动的转矩大适合高精度的控制。而交流电机调速需要变频，控制相对复杂一些，而对于设计无论是哪种电机都不影响结果，所以本实验采用直流电机。简单来说卡尔曼滤波就是对被观测量进行一个物理的建模，目的是用‘道理’来约束观测结果，减少噪声的影响。因此卡尔曼滤波是根据一个事物的当前状态预测它的下一个状态的过程。 此设计主要是通过对直流电机的数学模型利用MATLAB来设计卡尔曼滤波估计，进行仿真编程建模，进而对系统进行评估，并且分析估计误差。 关键词：卡尔曼滤波器；直流电机；MATLAB
Abstract Kalman filter is an iterative autoregression algorithm for continuous motion of most of the problems with it are able to give the best prediction. And it has been widely used for nearly half a century, such as the integration of data, as well as military machinery of navigation radar navigation, and so on. Kalman filter is generally used to process dynamic data, extract useful signal parameter estimation method to eliminate errors from the chaotic signal. Kalman filter is based on an estimate on the value and the current detection data is calculated using recursive estimation current estimates. By using recursive state equation method to estimate the movement of objects can be modeled. The paper describes the engineering design of the DC motor running state parameter calculation and verification. The DC motor performance and adjust very well simply by adding resistance regulator on it, and start curve is very good, start torque for precision control. The required frequency AC motor speed control is relatively complicated, and for the design of either the motor does not affect the outcome.In order to facilitate learning, so wo use the DC motor. Simply the Kalman filter is to be observables conduct a physical modeling; the purpose is to use 'sense' to restrict the observations to reduce the influence of noise. Therefore, the Kalman filter is based on the current state of things predict its next state of the process. This design is mainly through the DC motor mathematical model using MATLAB to design the Kalman filter estimation, simulation modeling program, and then to evaluate the system and analyze the estimation error. Keywords：Kalman filter; DC;MATLAB
卡尔曼滤波的原理及应用自己总结
卡尔曼滤波的原理以及应用 滤波，实质上就是信号处理与变换的过程。目的是去除或减弱不想要成分，增强所需成分。卡尔曼滤波的这种去除与增强过程是基于状态量的估计值和实际值之间的均方误差最小准则来实现的，基于这种准则，使得状态量的估计值越来越接近实际想要的值。而状态量和信号量之间有转换的关系，所以估计出状态量，等价于估计出信号量。所以不同于维纳滤波等滤波方式，卡尔曼滤波是把状态空间理论引入到对物理系统的数学建模过程中来，用递归方法解决离散数据线性滤波的问题，它不需要知道全部过去的数据，而是用前一个估计值和最近一个观察数据来估计信号的当前值，从而它具有运用计算机计算方便，而且可用于平稳和不平稳的随机过程(信号)，非时变和时变的系统的优越性。 卡尔曼滤波属于一种软件滤波方法，概括来说其基本思想是：以最小均方误差为最佳估计准则，采用信号与噪声的状态空间模型，利用前一时刻的估计值和当前时刻的观测值来更新对状态变量的估计，求出当前时刻的估计值，算法根据建立的系统方程和观测方程对需要处理的信号做出满足最小均方误差的估计。其所得到的解是以估计值的形式给出的。 卡尔曼滤波过程简单来说主要包括两个步骤：状态变量的预估以及状态变量的校正。预估过程是不考虑过程噪声和量测噪声，只是基于系统本身性质并依靠前一时刻的估计值以及系统控制输入的一种估计；校正过程是用量测值与预估量测值之间的误差乘以一个与过程
噪声和量测噪声相关的增益因子来对预估值进行校正的，其中增益因子的确定与状态量的均方误差有关，用到了使均方误差最小的准则。而这一过程中体现出来的递归思想即是：对于当前时刻的状态量估计值以及均方误差预估值实时进行更新，以便用于下一时刻的估计，使得系统在停止运行之前能够源源不断地进行下去。 下面对于其数学建模过程进行详细说明。 1.状态量的预估 （1）由前一时刻的估计值和送给系统的可控制输入来预估计当前时刻状态量。 X(k|k-1)=A X(k-1|k-1)+B U(k) 其中，X(k-1|k-1)表示前一时刻的估计值，U(k)表示系统的控制输入，X(k|k-1)表示由前一时刻估计出来的状态量的预估计值，A表示由k-1时刻过渡到k时刻的状态转移矩阵，B表示控制输入量与状态量之间的一种转换因子，这两个都是由系统性质来决定的。 （2）由前一时刻的均方误差阵来预估计当前时刻的均方误差阵。 P(k|k-1)=A P(k-1|k-1)A’+Q 其中，P(k-1|k-1)是前一时刻的均方误差估计值，A’代表矩阵A 的转置，Q代表过程噪声的均方误差矩阵。该表达式具体推导过程如下： P(k|k-1)=E{[Xs(k|k)-X(k|k-1)][Xs(k|k)-X(k|k-1)]’}------ 其中Xs(k|k)=A Xs(k-1|k-1)+B U(k)+W(k-1)表示当前时刻的实际值，Xs(k-1|k-1)表示前一时刻的实际值，可以看出与当前时刻的预估计值
卡尔曼滤波的基本原理及应用
卡尔曼滤波的基本原理及应用卡尔曼滤波在信号处理与系统控制领域应用广泛，目前，正越来越广泛地应用于计算机应用的各个领域。为了更好地理解卡尔曼滤波的原理与进行滤波算法的设计工作，主要从两方面对卡尔曼滤波进行阐述：基本卡尔曼滤波系统模型、滤波模型的建立以及非线性卡尔曼滤波的线性化。最后，对卡尔曼滤波的应用做了简单介绍。 卡尔曼滤波属于一种软件滤波方法，其基本思想是：以最小均方误差为最佳估计准则，采用信号与噪声的状态空间模型，利用前一时刻的估计值和当前时刻的观测值来更新对状态变量的估计，求出当前时刻的估计值，算法根据建立的系统方程和观测方程对需要处理的信号做出满足最小均方误差的估计。 最初的卡尔曼滤波算法被称为基本卡尔曼滤波算法，适用于解决随机线性离散系统的状态或参数估计问题。卡尔曼滤波器包括两个主要过程：预估与校正。预估过程主要是利用时间更新方程建立对当前状态的先验估计，及时向前推算当前状态变量和误差协方差估计的值，以便为下一个时间状态构造先验估计值；校正过程负责反馈，利用测量更新方程在预估过程的先验估计值及当前测量变量的基础上建立起对当前状态的改进的后验估计。这样的一个过程，我们称之为预估－校正过程，对应的这种估计算法称为预估－校正算法。以下给出离散卡尔曼滤波的时间更新方程和状态更新方程。 时间更新方程： 状态更新方程： 在上面式中，各量说明如下： A：作用在X k-1上的n×n 状态变换矩阵 B：作用在控制向量U k-1上的n×1 输入控制矩阵 H：m×n 观测模型矩阵，它把真实状态空间映射成观测空间 P k－：为n×n 先验估计误差协方差矩阵 P k：为n×n 后验估计误差协方差矩阵 Q：n×n 过程噪声协方差矩阵 R：m×m 过程噪声协方差矩阵 I：n×n 阶单位矩阵K k：n×m 阶矩阵，称为卡尔曼增益或混合因数 随着卡尔曼滤波理论的发展，一些实用卡尔曼滤波技术被提出来，如自适应滤波，次优滤波以及滤波发散抑制技术等逐渐得到广泛应用。其它的滤波理论也迅速发展，如线性离散系统的分解滤波（信息平方根滤波，序列平方根滤波，UD 分解滤波），鲁棒滤波（H∞波）。 非线性样条自适应滤波：这是一类新的非线性自适应滤波器，它由一个线性组合器后跟挠性无记忆功能的。涉及的自适应处理的非线性函数是基于可在学习
卡尔曼滤波器综述
卡尔曼滤波器综述 瞿伟军 G10074 1、卡尔曼滤波的起源 1960年，匈牙利数学家卡尔曼发表了一篇关于离散数据线性滤波递推算法的论文，这意味着卡尔曼滤波的诞生。斯坦利.施密特(Stanley Schmidt)首次实现了卡尔曼滤波器，卡尔曼在NASA埃姆斯研究中心访问时，发现他的方法对于解决阿波罗计划的轨道预测很有用，后来阿波罗飞船的导航电脑使用了这种滤波器。关于这种滤波器的论文由Swerling (1958)、Kalman (1960)与 Kalman and Bucy (1961)发表。 2、卡尔曼滤波的发展 卡尔曼滤波是一种有着相当广泛应用的滤波方法，但它既需要假定系统是线性的，又需要认为系统中的各个噪声与状态变量均呈高斯分布，而这两条并不总是确切的假设限制了卡尔曼滤波器在现实生活中的应用。扩展卡尔曼滤波器（EKF）极大地拓宽了卡尔曼滤波的适用范围。EKF的基本思路是，假定卡尔曼滤滤对当前系统状态估计值非常接近于其真实值，于是将非线性函数在当前状态估计值处进行台劳展开并实现线性化。另一种非线性卡尔曼滤波叫线性化卡尔曼滤波。它与EKF的主要区别是前者将非线函数在滤波器对当前系统状态的最优估计值处线性化，而后者因为预先知道非线性系统的实际运行状态大致按照所要求、希望的轨迹变化，所以这些非线性化函数在实际状态处的值可以表达为在希望的轨迹处的台劳展开式，从而完成线性化。 不敏卡尔曼滤波器（UKF）是针对非线性系统的一种改进型卡尔曼滤波器。UKF处理非线性系统的基本思路在于不敏变换，而不敏变换从根本上讲是一种描述高斯随机变量在非线性化变换后的概率分布情况的方法。不敏卡尔曼滤波认为，与其将一个非线性化变换线性化、近似化，还不如将高斯随机变量经非线性变换后的概率分布情况用高斯分布来近似那样简单，因而不敏卡尔曼滤波算法没
卡尔曼滤波器
卡尔曼滤波器 来这里几个月，发现有些问题很多人都很感兴趣。所以在这里希望能尽自己能力跟大家讨论一些力所能及的算法。现在先讨论一下卡尔曼滤波器，如果时间和能力允许，我还希望能够写写其他的算法，例如遗传算法，傅立叶变换，数字滤波，神经网络，图像处理等等。 因为这里不能写复杂的数学公式，所以也只能形象的描述。希望如果哪位是这方面的专家，欢迎讨论更正。 卡尔曼滤波器– Kalman Filter 1．什么是卡尔曼滤波器 （What is the Kalman Filter?） 在学习卡尔曼滤波器之前，首先看看为什么叫“卡尔曼”。跟其他著名的理论（例如傅立叶变换，泰勒级数等等）一样，卡尔曼也是一个人的名字，而跟他们不同的是，他是个现代人！ 卡尔曼全名Rudolf Emil Kalman，匈牙利数学家，1930年出生于匈牙利首都布达佩斯。1953，1954年于麻省理工学院分别获得电机工程学士及硕士学位。1957年于哥伦比亚大学获得博士学位。我们现在要学习的卡尔曼滤波器，正是源于他的博士论文和1960年发表的论文《A New Approach to Linear Filtering and Prediction Problems》（线性滤波与预测问题的新方法）。如果对这编论文有兴趣，可以到这里的地址下载：https://www.wendangku.net/doc/6c2113153.html,/~welch/media/pdf/Kalman1960.pdf。 简单来说，卡尔曼滤波器是一个“optimal recursive data processing algorithm（最优化自回归数据处理算法）”。对于解决很大部分的问题，他是最优，效率最高甚至是最有用的。他的广泛应用已经超过30年，包括机器人导航，控制，传感器数据融合甚至在军事方面的雷达系统以及导弹追踪等等。近年来更被应用于计算机图像处理，例如头脸识别，图像分割，图像边缘检测等等。 2．卡尔曼滤波器的介绍 （Introduction to the Kalman Filter） 为了可以更加容易的理解卡尔曼滤波器，这里会应用形象的描述方法来讲解，而不是像大多数参考书那样罗列一大堆的数学公式和数学符号。但是，他的5条公式是其核心内容。结合现代的计算机，其实卡尔曼的程序相当的简单，只要你理解了他的那5条公式。 在介绍他的5条公式之前，先让我们来根据下面的例子一步一步的探索。
卡尔曼滤波简介和实例讲解.
卡尔曼，美国数学家和电气工程师。1930年5月 19日生于匈牙利首都布达佩斯。1953年在美国麻省理工学院毕业获理学士学位,1954年获理学硕士学位,1957年在哥伦比亚大学获科学博士学位。1957～1958年在国际商业机器公司(IBM)研究大系统计算机控制的数学问题。1958～1964年在巴尔的摩高级研究院研究控制和数学问题。1964～1971年到斯坦福大学任教授。1971年任佛罗里达大学数学系统理论研究中心主任，并兼任苏黎世的瑞士联邦高等工业学校教授。1960年卡尔曼因提出著名的卡尔曼滤波器而闻名于世。卡尔曼滤波器在随机序列估计、空间技术、工程系统辨识和经济系统建模等方面有许多重要应用。1960年卡尔曼还提出能控性的概念。能控性是控制系统的研究和实现的基本概念，在最优控制理论、稳定性理论和网络理论中起着重要作用。卡尔曼还利用对偶原理导出能观测性概念，并在数学上证明了卡尔曼滤波理论与最优控制理论对偶。为此获电气与电子工程师学会(IEEE)的最高奖──荣誉奖章。卡尔曼著有《数学系统概论》(1968)等书。 什么是卡尔曼滤波 最佳线性滤波理论起源于40年代美国科学家Wiener和前苏联科学家Ｋолмогоров等人的研究工作，后人统称为维纳滤波理论。从理论上说，维纳滤波的最大缺点是必须用到无限过去的数据，不适用于实时处理。为了克服这一缺点，60年代Kalman把状态空间模型引入滤波理论，并导出了一套递推估计算法，后人称之为卡尔曼
滤波理论。卡尔曼滤波是以最小均方误差为估计的最佳准则，来寻求一套递推估计的算法，其基本思想是：采用信号与噪声的状态空间模型，利用前一时刻地估计值和现时刻的观测值来更新对状态变量的估计，求出现时刻的估计值。它适合于实时处理和计算机运算。 卡尔曼滤波的实质是由量测值重构系统的状态向量。它以“预测—实测—修正”的顺序递推，根据系统的量测值来消除随机干扰，再现系统的状态，或根据系统的量测值从被污染的系统中恢复系统的本来面目。 释文：卡尔曼滤波器是一种由卡尔曼（Kalman）提出的用于时变线性系统的递归滤波器。这个系统可用包含正交状态变量的微分方程模型来描述，这种滤波器是将过去的测量估计误差合并到新的测量误差中来估计将来的误差。 卡尔曼滤波的应用 斯坦利.施密特(Stanley Schmidt)首次实现了卡尔曼滤波器.卡尔曼在NASA埃姆斯研究中心访问时,发现他的方法对于解决阿波罗计划的轨道预测很有用,后来阿波罗飞船的导航电脑使用了这种滤波器. 关于这种滤波器的论文由Swerling (1958), Kalman (1960)与 Kalman and Bucy (1961)发表.
关于卡尔曼滤波器的一个简单介绍
关于卡尔曼滤波器的一个简单介绍 在实际生产过程中，我们经常需要使用某种仪器测量某一物理参数。比如，用一台仪器测量大气中一氧化碳的浓度。由于测量误差永远存在，使得这样两个问题非常突出，首先，是否存在某个计算方法，能够从含有误差的测量结果中获得比较接近真实值的结果；其次，如何证明这样的结果是最优的，也就是说有没有这样一种数学方法，经过其进行处理后得到的测量结果是最接近真实值的。 为了消除测量误差，人们首先想到的方法是取平均值。通过对于同一物理量的多次测量，抵消可能存在的测量误差，从而得到真实值。但是，通过生产实践，人们很快发现这样的测量方法并不是最优的。如果参与平均值计算的数据量太小，就达不到抵消测量误差的目的，如果参与计算的数据量太多，不仅完全消除了可能的物理量变化，而且实现起来非常麻烦。 1960年，匈牙利数学家卡尔曼（Rudolf Emil Kalman，1930年出生于匈牙利首都布达佩斯。1953，1954年于麻省理工学院分别获得电机工程学士及硕士学位，1957年于哥伦比亚大学获得博士学位）在他的博士论文和1960年发表的论文《A New Approach to Linear Filtering and Prediction Problems》（线性滤波与预测问题的新方法）中提出一种计算方法，后来被称为卡尔曼滤波器。这种滤波器比较简单，但是对于绝大多数类似本文开头提出的问题，可以证明，它得出的结果是最优的。卡尔曼滤波器的最大优势是它非常简单，非常容易实现。故而在各领域广泛应用超过三十年。 下面简单介绍卡尔曼滤波器，为了便于阅读，尽量少使用数学公式。下面的用词不是非常严格的。 假如，我们需要用一个温度计测量一个房间的温度。每个1秒钟从温度计上读取一个数值。假设这个房间的真实温度是25度。 由于温度计本身有测量误差，我们得到的测量值是围绕25度上下波动的一组数值。“围绕”这一特性是符合我们的直观感受的。也就是说虽然温度计上的读数在不停的变化，但是我们知道温度计上出现26度或24度的可能性要比出现30度或20度的可能性大。这种特性在数学上被称为高斯分布。当有很多因素影响测量结果，而每种干扰因素都不起主导作用时，测量的结果呈现高斯分布。对于温度计来讲，制造工艺，空气扰动，我们读数时候的误差都可能影响测量结果。所以，温度计的读取结果是符合高斯分布的。高斯分布有两个参数，一个叫做期望，在这里可以等同于平均值，一个叫做方差。方差反映了数据的集中程度。比如一个温度计测量的结果分布在24到26度之间，另一个温度计测量结果在20到30度之间，显然，前一个温度计更好一些。它的数据方差就比较小。这里假设我们温度计的方差是3。于是，温度计的测量结果就是期望是25，方差是3的一组数据。 另外，我们可以根据经验猜测一下房间的温度（就好像把自己当成一个人体温度计），比如，我们猜测房间的温度是23度，当然，我们自己猜测的结果是有方差的，我们假设这个方差是4，于是根据人的感觉，测量结果是期望是23，方差是4的一组数据。 然后，由经验我们知道房间温度变化不大，基本可以认为，这一秒的温度和下一秒的温度是相等的。 现在，在第1秒，我们有了两个值，一个是温度计上显示的数据，25。一个是我们自己 根据感觉得出来的数据，23。那么相信谁呢？这里有一个数学公式 2 2 22 3 34 g k= + ，可以算出 kg＝0.6，于是25×0.6+23×（1－0.6）＝24.2，这个24.2就是我们在第1秒得到的测量结果。也就是说第1秒的温度是24.2度。 现在，我们进入了第2秒，这时，由于我们已经知道了第1秒的温度是24.2度，所以，
matlab对卡尔曼滤波的仿真实现
MATLAB 对卡尔曼滤波器的仿真实现 刘丹，朱毅，刘冰 武汉理工大学信息工程学院，武汉（430070） E-mail ：liudan_ina@https://www.wendangku.net/doc/6c2113153.html, 摘 要：本文以卡尔曼滤波器原理为理论基础，用MATLAB 进行卡尔曼滤波器仿真、对比卡尔曼滤波器的预测效果，对影响滤波其效果的各方面原因进行讨论和比较，按照理论模型进行仿真编程，清晰地表述了编程过程。 关键词：数字信号处理；卡尔曼滤波器；MATLAB ；仿真过程 中图分类号: TN912.3 1. 引言 随着信息时代和数字世界的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。数字信号处理已在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理中,数字滤波占有极其重要的地位,目前对数字滤波器的设计有多种方法,其中著名的MATLAB 软件包在多个研究领域都有着广泛的应用,它的频谱分析[1]和滤波器的分析设计功能很强,从而使数字信号处理变得十分简单、直观。本文分析了数字滤波器的设计方法,举出了基于MATLAB 软件的信号处理工具在数字滤波器设计中的应用。 2. 卡尔曼滤波基本原理 卡尔曼滤波过程实际上是获取维纳解的递推运算过程[2]。从维纳解导出的卡尔曼滤波器实际上是卡尔曼滤波过程结束后达到稳态的情况，这时Kalman Filtering 的结果与Wiener Solution 是相同的[3]。具体推导如下： )()1|1(?)|(?n Gy n n x f n n x +??= )|(?)()(n n x n x n e ?= 已知由此求c a cG a f F G n e E n ,)1(( ..min )]([)(2?=??→?==ε 由 f G f G ,0??????????=??εε ⑴ )]1|1(?)()[()1|1(?)|(????+??=n n x ac n y n G n n x a n n x 可以是时变的，非平稳的随机信号 ⑵ Q n a n P +?=)1()(2 ε均为正数。 ⑶ ) () ()(2n P C R n CP n G += ⑷ )()](1[)()(n P n CG n G C P n ??== ε )(n G 是个随时间变化的量，每次输入输出，)(n G 就调整一次，并逐渐逼近Kalman Filter 的增益G ，而)1()(?