卡尔曼滤波.docx
- 文档编号:25827750
- 上传时间:2023-06-15
- 格式:DOCX
- 页数:21
- 大小:81.25KB
卡尔曼滤波.docx
《卡尔曼滤波.docx》由会员分享,可在线阅读,更多相关《卡尔曼滤波.docx(21页珍藏版)》请在冰豆网上搜索。
卡尔曼滤波
卡尔曼滤波
卡尔曼滤波(Kalmanfiltering)一种利用线性系统状态方程,通过系统输入输出观测数据,对系统状态进行最优估计的算法。
由于观测数据中包括系统中的噪声和干扰的影响,所以最优估计也可看作是滤波过程。
斯坦利施密特(StanleySchmidt)首次实现了卡尔曼滤波器。
卡尔曼在NASA埃姆斯研
究中心访问时,发现他的方法对于解决阿波罗计划的轨道预测很有用,后来阿波罗飞船的导
航电脑使用了这种滤波器。
关于这种滤波器的论文由Swerling(1958),Kalman(I960)与
KalmanandBucy(1961)发表。
数据滤波是去除噪声还原真实数据的一种数据处理技术,Kalman滤波在测量方差已知
的情况下能够从一系列存在测量噪声的数据中,估计动态系统的状态•由于,它便于计算机
编程实现,并能够对现场采集的数据进行实时的更新和处理,Kalman滤波是目前应用最为
广泛的滤波方法,在通信,导航,制导与控制等多领域得到了较好的应用•
中文名
卡尔曼滤波器,Kalman滤波,卡曼滤波
外文名
KALMANFILTER
表达式
X(k)=AX(k-1)+BU(k)+W(k)
提岀者
斯坦利施密特
提岀时间
1958
应用学科
天文,宇航,气象
适用领域范围
雷达跟踪去噪声
适用领域范围
控制、制导、导航、通讯等现代工程
斯坦利施密特(StanleySchmidt)首次实现了卡尔曼滤波器。
卡尔曼在NASA埃姆斯研
究中心访问时,发现他的方法对于解决阿波罗计划的轨道预测很有用,后来阿波罗飞船的导—
航电脑使用了这种滤波器。
关于这种滤波器的论文由Swerling(1958),Kalman(1960)与
KalmanandBucy(1961)发表。
2定义
传统的滤波方法,只能是在有用信号与噪声具有不同频带的条件下才能实现.20世纪
40年代,N.维纳和A.H.柯尔莫哥罗夫把信号和噪声的统计性质引进了滤波理论,在假设信号和噪声都是平稳过程的条件下,利用最优化方法对信号真值进行估计,达到滤波目的,
从而在概念上与传统的滤波方法联系起来,被称为维纳滤波。
这种方法要求信号和噪声都必
须是以平稳过程为条件。
60年代初,卡尔曼(R.E.Kalman)和布塞(R.S.Bucy)发表了
一篇重要的论文《线性滤波和预测理论的新成果》,提出了一种新的线性滤波和预测理由论,被称之为卡尔曼滤波。
特点是在线性状态空间表示的基础上对有噪声的输入和观测信号进行处理,求取系统状态或真实信号。
这种理论是在时间域上来表述的,基本的概念是:
在线性系统的状态空间表示基础上,
从输出和输入观测数据求系统状态的最优估计。
这里所说的系统状态,是总结系统所有过去
的输入和扰动对系统的作用的最小参数的集合,知道了系统的状态就能够与未来的输入与系
统的扰动一起确定系统的整个行为。
卡尔曼滤波不要求信号和噪声都是平稳过程的假设条件。
对于每个时刻的系统扰动和观测误差(即噪声),只要对它们的统计性质作某些适当的假定,通过对含有噪声的观测信号进行处理,就能在平均的意义上,求得误差为最小的真实信号的估计值。
因此,自从卡尔
曼滤波理论问世以来,在通信系统、电力系统、航空航天、环境污染控制、工业控制、雷达信号处理等许多部门都得到了应用,取得了许多成功应用的成果。
例如在图像处理方面,应
用卡尔曼滤波对由于某些噪声影响而造成模糊的图像进行复原。
在对噪声作了某些统计性质
的假定后,就可以用卡尔曼的算法以递推的方式从模糊图像中得到均方差最小的真实图像,使模糊的图像得到复原。
3性质
1卡尔曼滤波是一个算法,它适用于线性、离散和有限维系统。
每一个有外部变量的自回归移动平均系统(ARMAX)或可用有理传递函数表示的系统都可以转换成用状态空间表示的系统,从而能用卡尔曼滤波进行计算。
2任何一组观测数据都无助于消除x(t)的确定性。
增益K(t)也同样地与观测数据无关。
3当观测数据和状态联合服从高斯分布时用卡尔曼递归公式计算得到的是高斯随机变量的条件均值和条件方差,从而卡尔曼滤波公式给出了计算状态的条件概率密度的更新过程线性最小方差估计,也就是最小方差估计。
卡尔曼滤波已经有很多不同的实现,卡尔曼最初提出的形式一般称为简单卡尔曼滤波
器。
除此以外,还有施密特扩展滤波器、信息滤波器以及很多Bierman,Thornton开发的平方根滤波器的变种。
最常见的卡尔曼滤波器是锁相环,它在收音机、计算机和几乎任何视频
或通讯设备中广泛存在。
5实例
■
卡尔曼滤波的一个典型实例是从一组有限的,对物体位置的,包含噪声的观察序列中
预测出物体的坐标位置及速度。
在很多工程应用(雷达、计算机视觉)中都可以找到它的身影。
同时,卡尔曼滤波也是控制理论以及控制系统工程中的一个重要话题。
比如,在雷达中,人们感兴趣的是跟踪目标,但目标的位置、速度、加速度的测量值往往在任何时候都有噪声。
卡尔曼滤波利用目标的动态信息,设法去掉噪声的影响,得到一
个关于目标位置的好的估计。
这个估计可以是对当前目标位置的估计(滤波),也可以是对于
将来位置的估计(预测),也可以是对过去位置的估计(插值或平滑)。
扩展卡尔曼滤波(EXTENDKALMANFILTER,EKF)
扩展卡尔曼滤波器
是由kalmanfilter考虑时间非线性的动态系统,常应用于目标跟踪系统。
7状态估计
状态估计
状态估计是卡尔曼滤波的重要组成部分。
一般来说,根据观测数据对随机量进行定量推断就是估计问题,特别是对动态行为的状态估计,它能实现实时运行状态的估计和预测功
能。
比如对飞行器状态估计。
状态估计对于了解和控制一个系统具有重要意义,所应用的方
法属于统计学中的估计理论。
最常用的是最小二乘估计,线性最小方差估计、最小方差估计、递推最小二乘估计等。
其他如风险准则的贝叶斯估计、最大似然估计、随机逼近等方法也都有应用。
状态量
受噪声干扰的状态量是个随机量,不可能测得精确值,但可对它进行一系列观测,并
依据一组观测值,按某种统计观点对它进行估计。
使估计值尽可能准确地接近真实值,这就是最优估计。
真实值与估计值之差称为估计误差。
若估计值的数学期望与真实值相等,这种估计称为无偏估计。
卡尔曼提出的递推最优估计理论,采用状态空间描述法,在算法采用递推形式,卡尔曼滤波能处理多维和非平稳的随机过程。
理论
卡尔曼滤波理论的提出,克服了威纳滤波理论的局限性使其在工程上得到了广泛的应用,尤其在控制、制导、导航、通讯等现代工程方面。
8MATLABS序
MATLAB
N=200;
w
(1)=0;
w=randn(1,N)
x
(1)=0;
a=1;
fork=2:
N;
x(k)=a*x(k_1)+w(k_1);
end
V=randn(1,N);
q1=std(V);
Rvv=q1.A2;
q2=std(x);
Rxx=q2.A2;
q3=std(w);
Rww=q3.A2;
c=0.2;
Y=c*x+V;
p
(1)=0;
s
(1)=0;
fort=2:
N;
p1(t)=a.A2*p(t-1)+Rww;
b(t)=c*p1(t)/(c.A2*p1(t)+Rvv);
s(t)=a*s(t-1)+b(t)*(Y(t)-a*c*s(t-1));
p(t)=p1(t)-c*b(t)*p1(t);
end
t=1:
N;
plot(t,s,'r',t,Y,‘g',t,x,'b');
function[x,V,VV,loglik]=kalman_filter(y,A,C,Q,R,init_x,init_V,varargin)%Kalmanfilter.
%[x,V,VV,loglik]=kalman_filter(y,A,C,Q,R,init_x,init_V,...)
%
%INPUTS:
%y(:
t)-theobservationattimet
%A-thesystemmatrix
%C-theobservationmatrix
%Q-thesystemcovarianee
%R-theobservationcovarianee%init_x-theinitialstate(column)vector
%init_V-theinitialstatecovarianee
%
%OPTIONALINPUTS(string/valuepairs[defaultinbrackets])
%'model'-model(t)=mmeansuseparamsfrommodelmattimet[ones(1,T)]
%Inthiscase,alltheabovematricestakeanadditionalfinaldimension,
%i.e.,A(:
:
m),C(:
:
m),Q(:
:
m),R(:
:
m).
%However,init_xandinit_Vareindependentofmodel
(1).
%'u'-u(:
t)thecontrolsignalattimet[[]]
%'B'-B(:
,:
m)theinputregressionmatrixformodelm
%
%OUTPUTS(whereXisthehiddenstatebeingestimated)
%x(:
t)=E[X(:
t)|y(:
1:
t)]
%V(:
,:
t)=Cov[X(:
t)|y(:
1:
t)]
%VV(:
,:
t)=Cov[X(:
t),X(:
t-1)|y(:
1:
t)]t>=2
%loglik=sum{t=1}ATlogP(y(:
t))
%
%Ifaninputsignalisspecified,wealsoconditiononit:
%e.g.,x(:
t)=E[X(:
t)|y(:
1:
t),u(:
1:
t)]
%Ifamodelsequeneeisspecified,wealsoconditiononit:
%e.g.,x(:
t)=E[X(:
t)|y(:
1:
t),u(:
1:
t),m(1:
t)]
[osT]=size(y);
ss=size(A,1);%sizeofstatespace
%setdefaultparams
model=ones(1,T);
u=[];
B=[];
ndx=[];
args=varargin;
nargs=length(args);
fori=1:
2:
nargs
switchargs
case'model',model=args{i+1};
case'u',u=args{i+1};
case'B',B=args{i+1};
case'ndx',ndx=args{i+1};
otherwise,error(['unrecognizedargument'args])
endend
x=zeros(ss,T);
V=zeros(ss,ss,T);
VV=zeros(ss,ss,T);
loglik=0;
fort=1:
T
m=model(t);
ift==1
%prevx=init_x(:
m);
%prevV=init_V(:
,:
m);
prevx=init_x;
prevV=init_V;
initial=1;
else
prevx=x(:
t-1);
prevV=V(:
,:
t-1);
initial=0;
end
ifisempty(u)
[x(:
t),V(:
,:
t),LL,VV(:
,:
t)]=...
kalman_update(A(:
:
m),C(:
,:
m),Q(:
,:
m),R(:
,:
m),y(:
t),prevx,prevV,'initial',initial);else
ifisempty(ndx)
[x(:
t),V(:
,:
t),LL,VV(:
,:
t)]=...
kalman_update(A(:
:
m),C(:
,:
m),Q(:
,:
m),R(:
,:
m),y(:
t),prevx,prevV,...
'initial',initial,'u',u(:
t),'B',B(:
,:
m));
else
i=ndx;
%copyoverallelements;onlysomewillgetupdated
x(:
t)=prevx;
prevP=inv(prevV);
prevPsmall=prevP(i,i);
prevVsmall=inv(prevPsmall);
[x(i,t),smallV,LL,VV(i,i,t)]=...
kalman_update(A(i,i,m),C(:
i,m),Q(i,i,m),R(:
,:
m),y(:
t),prevx(i),prevVsmall,...
'initial',initial,'u',u(:
t),'B',B(i,:
m));
smallP=inv(smallV);
prevP(i,i)=smallP;
V(:
:
t)=inv(prevP);
end
end
loglik=loglik+LL;
end
另一个例子:
⑴
Z=(1:
100);%观测值
noise=randn(1,100);%方差为1的高斯噪声
Z=Z+noise;
X=[0;0];%状态
P=[10;01];%状态协方差矩阵F=[11;01];%状态转移矩阵
Q=[0.0001,0;00.0001];%状态转移协方差矩阵H=[10];%观测矩阵
R=1;%观测噪声方差
figure;
holdon;
fori=1:
100
X_=F*X;
P_=F*P*F'+Q;
K=P_*H'/(H*P_*H'+R);
X=X_+K*(Z(i)-H*X_);
P=(eye
(2)-K*H)*P_;
plot(X
(1),X
(2));%画点,横轴表示位置,纵轴表示速度
end
9通俗解释
简单来说,卡尔曼滤波器是一个“optimalrecursivedataprocessingalgorithm(最优化
自回归数据处理算法)”。
对于解决很大部分的问题,他是最优,效率最高甚至是最有用的。
他的广泛应用已经超过30年,包括机器人导航,控制,传感器数据融合甚至在军事方面的
雷达系统以及导弹追踪等等。
近来更被应用于计算机图像处理,例如头脸识别,图像分割,图像边缘检测等等。
卡尔曼滤波器的介绍:
为了可以更加容易的理解卡尔曼滤波器,这里会应用形象的描述方法来讲解,而不是像大多数参考书那样罗列一大堆的数学公式和数学符号。
但是,他的5条公式是其核心内
容。
结合现代的计算机,其实卡尔曼的程序相当的简单,只要你理解了他的那5条公式。
在介绍他的5条公式之前,先让我们来根据下面的例子一步一步的探索。
假设我们要研究的对象是一个房间的温度。
根据你的经验判断,这个房间的温度是恒
定的,也就是下一分钟的温度等于现在这一分钟的温度(假设我们用一分钟来做时间单位)。
假设你对你的经验不是100%的相信,可能会有上下偏差几度。
我们把这些偏差看成是高斯
白噪声(WhiteGaussianNoise),也就是这些偏差跟前后时间是没有关系的而且符合高斯分布(GaussianDistribution)。
另外,我们在房间里放一个温度计,但是这个温度计也不准确的,测量值会比实际值偏差。
我们也把这些偏差看成是高斯白噪声。
好了,现在对于某一分钟我们有两个有关于该房间的温度值:
你根据经验的预测值(系
统的预测值)和温度计的值(测量值)。
下面我们要用这两个值结合他们各自的噪声来估算
出房间的实际温度值。
假如我们要估算k时刻的实际温度值。
首先你要根据k-1时刻的温度值,来预测k时
刻的温度。
因为你相信温度是恒定的,所以你会得到k时刻的温度预测值是跟k-1时刻一样的,假设是23度,同时该值的高斯噪声的偏差是5度(5是这样得到的:
如果k-1时刻估
算出的最优温度值的偏差是3,你对自己预测的不确定度是4度,他们平方相加再开方,就
是5)。
然后,你从温度计那里得到了k时刻的温度值,假设是25度,同时该值的偏差是4
度。
由于我们用于估算k时刻的实际温度有两个温度值,分别是23度和25度。
究竟实际温度是多少呢?
相信自己还是相信温度计呢?
究竟相信谁多一点,我们可以用他们的协方差
(covarianee)来判断。
因为Kg=5A2/(5A2+4A2),所以Kg=0.61,我们可以估算出k时刻的实际温度值是:
23+0.61*(25-23)=24.22度。
可以看出,因为温度计的协方差(covarianee)比
较小(比较相信温度计),所以估算出的最优温度值偏向温度计的值。
现在我们已经得到k时刻的最优温度值了,下一步就是要进入k+1时刻,进行新的最
优估算。
到现在为止,好像还没看到什么自回归的东西出现。
对了,在进入k+1时刻之前,
我们还要算出k时刻那个最优值(24.22度)的偏差。
算法如下:
((1-勺)*5人2)人0.5=3.12。
这里的5就是上面的k时刻你预测的那个23度温度值的偏差,得出的3.12就是进入k+1时刻以后k时刻估算出的最优温度值的偏差(对应于上面的3)。
就是这样,卡尔曼滤波器就不断的把(协方差(covarianee)递归,从而估算出最优的温
度值。
他运行的很快,而且它只保留了上一时刻的协方差(covarianee)。
上面的Kg,就是卡
尔曼增益(KalmanGain)。
他可以随不同的时刻而改变他自己的值,是不是很神奇!
下面就要言归正传,讨论真正工程系统上的卡尔曼。
卡尔曼滤波器算法:
在这一部分,我们就来描述源于DrKalman的卡尔曼滤波器。
下面的描述,会涉及一
些基本的概念知识,包括概率(Probability),随机变量(RandomVariable),高斯或正态
分配(GaussianDistribution)还有State-spaceModel等等。
但对于卡尔曼滤波器的详细证明,这里不能一一描述。
首先,我们先要引入一个离散控制过程的系统。
该系统可用一个线性随机微分方程
(LinearStochasticDifferenceequation)来描述:
X(k)=AX(k-1)+BU(k)+W(k)
再加上系统的测量值:
Z(k)=HX(k)+V(k)
上两式子中,X(k)是k时刻的系统状态,U(k)是k时刻对系统的控制量。
A和B是系
统参数,对于多模型系统,他们为矩阵。
Z(k)是k时刻的测量值,H是测量系统的参数,对
于多测量系统,H为矩阵。
W(k)和V(k)分别表示过程和测量的噪声。
他们被假设成高斯白噪声(WhiteGaussianNoise),他们的协方差(covarianee)分别是Q,R(这里我们假设他们不随系统状态变化而变化)。
对于满足上面的条件(线性随机微分系统,过程和测量都是高斯白噪声),卡尔曼滤波器是最优的信息处理器。
下面我们结合他们的协方差来估算系统的最优化输出(类似上一节那
个温度的例子)。
首先我们要利用系统的过程模型,来预测下一状态的系统。
假设现在的系统状态是k,
根据系统的模型,可以基于系统的上一状态而预测出现在状态:
X(k|k-1)=AX(k-1|k-1)+BU(k).
(1)
式⑴中,X(k|k-1)是利用上一状态预测的结果,X(k-1|k-1)是上一状态最优的结果,U(k)
为现在状态的控制量,如果没有控制量,它可以为0。
到现在为止,我们的系统结果已经更新了,可是,对应于X(k|k-1)的协方差还没更新。
我们用P表示协方差(covarianee):
P(k|k-1)=AP(k-1|k-1)A'+Q
(2)
式⑵中,P(k|k-1)是X(k|k-1)对应的协方差,P(k-1|k-1)是X(k-1|k-1)对应的协方差,A
表示A的转置矩阵,Q是系统过程的协方差。
式子1,2就是卡尔曼滤波器5个公式当中的
前两个,也就是对系统的预测。
现在我们有了现在状态的预测结果,然后我们再收集现在状态的测量值。
结合预测值和测量值,我们可以得到现在状态(k)的最优化估算值X(k|k):
X(k|k)=X(k|k-1)+Kg(k)(Z(k)-HX(k|k-1))(3)
其中Kg为卡尔曼增益(KalmanGain):
Kg(k)=P(k|k-1)H'/(HP(k|k1)H'+R)(4)
到现在为止,我们已经得到了k状态下最优的估算值X(k|k)。
但是为了要令卡尔曼滤
波器不断的运行下去直到系统过程结束,我们还要更新k状态下X(k|k)的协方差:
P(k|k)=(I-Kg(k)H)P(k|k-1)(5)
其中I为1的矩阵,对于单模型单测量,1=1。
当系统进入k+1状态时,P(k|k)就是式子⑵的P(k-1|k-1)。
这样,算法就可以自回归的运算下去。
卡尔曼滤波器的原理基本描述了,式子1,2,3,4和5就是他的5个基本公式。
根
据这5个公式,可以很容易用计算机编程实现。
在上面的例子中,过程误差和测量误差设定为4是为了讨论的方便。
实际中,温度的
变化速度以及温度计的测量误差都没有这么大。
假设如下一个系统:
•房间内连续两个时刻温度差值的标准差为0.02度
•温度计的测量值误差的标准差为0.5度
•房间温度的真实值为24度
•对温度的初始估计值为23.5度,误差的方差为1
MatLab仿真的代码如下:
%KalmanfilterexampleoftemperaturemeasurementinMatlab
%ThisMcodeismodifiedfromXuchenYao'smatlabon2013/4/18
%房间当前温度真实值为24度,认为下一时刻与当前时刻温度相同,误差为0.02度
(即认为连续的两个时刻最多变化0.02度)。
%温度计的测量误差为0.5度。
%开始时,房间温度的估计为23.5度,误差为1度。
%Kalmanfilter
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 卡尔 滤波