最优估计与卡尔曼滤波基本原理

本文详细介绍了估计理论的核心概念——最优估计,包括均方误差的计算和比较。接着,深入探讨了卡尔曼滤波的原理,阐述了卡尔曼滤波方程及其在状态估计中的应用。同时,提到了滤波初值的选择策略以及连续时间随机系统的离散化方法。通过对不同随机过程的离散化示例,展示了如何将连续时间模型转换为离散时间模型,以便于实际应用中进行滤波处理。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

一、最优估计

1、估计理论:

估计理论是概率论与数理统计的一个分支,是根据受扰动的观测数据来提取系统某些参数或状态的一种数学方法。

定义估计量\breve{X}的误差为:

 对于随机变量,常用均方误差(Mean Square Error,MSE)来衡量估计量与被估计量之间的误差大小,类似的,对于随机向量定义均方误差阵为:

 特别地,如果估计\breve{X}是无偏的,则均方误差阵就等于其误差的方差阵,即有MSE[\breve{X}]=var[\breve{X}]成立。一般也称标量函数E[\breve{X^{T}\breve{X}}]为均方误差,显然有:

 对两个均方误差阵进行大小比较,通常是指比较均方误差阵之迹的大小,它等价于比较均方误差的大小。

二、卡尔曼滤波方程

1、kalman方程:

<1>随机系统状态空间模型:

其中:

X_{k}是n维的状态向量;

Z_{k}是m维的量测向量;

W_{k-1}l维的系统噪声向量;

V_{k}m维的量测噪声向量;

\Phi _{k/k-1}n维的状态一步转移矩阵;

\Gamma _{k/k-1}n\times l维的系统噪声分配矩阵;

H_{k}m\times n维的量测矩阵;

 <2>kalman滤波的五个公式:

Ⅰ状态一步预测:

Ⅱ 状态一步预测均方误差阵:

Ⅲ 滤波增益:


Ⅳ状态估计:

Ⅴ状态估计均方误差阵:

 注:

Ⅰ:量测一步预测均方误差阵P_{ZZ k/k-1}

Ⅱ:状态一步预测与量测一步预测之间的协均方误差阵P_{XZ k/k-1}

 Ⅲ:滤波增益计算公式,等价:

Ⅳ: K时刻,状态估计均方误差阵:

2、滤波流程框图和滤波初值的选择:

 

量测信息Z_{k}是Kalman滤波的最主要输入,但对于时变系统而言,系统结构参数\Phi _{k/k-1}\Gamma _{k/k-1}H_{k}及噪声Q_{k-1}R_{k}中的全部或部分是时变的,也是为滤波算法的输入,需实时更新。除了状态估计\breve{X}外,状态估计均方误差阵P_{k}也是Kalman滤波输出的重要组成部分,P_{k}评价状态估计的质量发挥着非常重要的作用。

tips:

一般将初始均方误差阵P_{0}设置为对角矩阵,各对角线元素的平方根粗略反映了状态分量初值的不确定度。

(1)对于可观测性较强的状态分量,对应的状态初值和均方误差阵设置偏差容许适当大些, 它们随着滤波更新将会快速收敛,如果均方误差阵设置太小,则会使收敛速度变慢。

(2)对于可观测性较弱的状态, 对应的状态初值和均方误差阵应该设置尽量准确,如果均方误差阵设置过大, 容易引起状态估计的剧烈波动, 反之,如果均方误差阵设置过小,同样会使状态收敛速度变慢,这两种情况下均方误差阵都不宜用于评估相应状态估计的精度。

(3)对于不可观测的状态分量, 其状态估计及其均方误差阵不会随滤波更新而变化, 即不会有滤波效果。 

 3、连续时间随机系统的离散化

<1>系统方程

连续时间线性随机系统(线性随机微分方程):

 离散化如下1:

其中:

注意到Q_{k}=q(t_{k-1})T_{s},它与等效噪声\eta _{k-1}W_{k}X_{k})方差的量纲单位是一致的;

 离散化如下2:

其中:

--%%--%%--%%--%%--%%--%%--%%--%%--%%--%%--%%--%%--%%--%%--%%--%%--%%--

几种典型的连续时间随机过程及其离散化结果:

(1)一阶马尔可夫过程

连续时间:

离散化:

其中:

(2)二阶马尔可夫过程:

连续时间:

离散化:

其中:

 (3)随机游走:

连续时间:

离散化:

其中:

(4)随机常值:

连续时间:

离散化:

 <2>量测方程:

在实际应用中, 大多数系统的量测方程是以离散形式直接给出的, 无需再进行离散化处理。

值得一提的是, 如果在一定量测频率范围内量测噪声的方差大小基本不变,则在量测设备允许的情况下选用较高的量测频率对提高滤波估计精度是有益的;

另一方面,如果系统状态变化比较平缓, 为了减小量测更新频率和计算量,则可将相继多次量测作平均处理, 并相应减少量测噪声大小, 利用平均量测进行滤波量测更新与进行多次量测更新是基本等效的。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值