• 状态估计理论的目的是
    • 对目标过去的运动状态进行平滑
    • 对目标现在的运动状态滤波
    • 对目标未来的运动状态进行预测
  • 这些运动状态包括目标位置、速度和加速度等
  • 多传感器跟踪系统中广泛应用的状态估计技术,是多传感器信息系统的最基本要素,也是形成多目标自适应跟踪滤波的前提和基础

    2.1 估计问题的构成

  • 用 n 维向量表示被估计参数

【数据融合】状态估计理论 - 图1

  • 用 m 维向量表示量测值

【数据融合】状态估计理论 - 图2

  • 通常假设量测值 Z 与估计量 X 满足关系式

【数据融合】状态估计理论 - 图3

  • 这里 j 是指 j 时刻, v(j) 是 l 维量测噪声并满足一定概率分布,通常是满足 Gauss 白噪声
  • 经过 k 个时刻的量测

【数据融合】状态估计理论 - 图4

  • 对 X(k) 进行估计,其估计值为【数据融合】状态估计理论 - 图5。当被估计量 X 不随时间变化时,则称对参数 X 的估计为静态估计
  • 当参数 X 随时间变化时,一般认为 X 满足某一动态方程

【数据融合】状态估计理论 - 图6

  • 其中 U(k) 是 p 维输入向量,W(k) 是 q 维过程噪声,也满足一定的概率分布,通常为Gauss白噪声
  • 更一般的是连续型方程

【数据融合】状态估计理论 - 图7

  • U(t) 是输入向量,W(t) 是过程噪声
  • 要解决的问题是
  • 已知 k 次量测值【数据融合】状态估计理论 - 图8,对 k 时刻的状态向量 X(k) 进行最优估计

    2.2 状态估计问题

  • 状态与系统相联系

  • 状态估计是指对动态随机系统的状态估计
  • 设有动态系统,它满足一定的数学模型

【数据融合】状态估计理论 - 图9

  • 其有关随机向量满足一定的统计性质
  • 所指系统的状态估计问题就是根据选定的估计准则和已获得的量测信息

【数据融合】状态估计理论 - 图10
【数据融合】状态估计理论 - 图11

  • 对系统进行估计,其中状态方程确定了被估计量的随机状态过程
  • 状态向量与量测量之间在时间上有不同的对应关系
  • 【数据融合】状态估计理论 - 图12为已知 j 和 j 以前时刻的量测值,对 k 时刻状态 X(k) 作出某种估计

    • 当k=j时,称为滤波问题,称【数据融合】状态估计理论 - 图13为X(k)的最优滤波估计量
    • 当k>j时,称为预测问题,称【数据融合】状态估计理论 - 图14为X(k)的最优预测估计量
    • 当k

      2.3 离散线性系统的最优估计

      2.3.1 线性关系

  • 动态系统是离散的和线性的,即式【数据融合】状态估计理论 - 图15中的 f 是线性函数

  • 状态方程满足

【数据融合】状态估计理论 - 图16

  • 【数据融合】状态估计理论 - 图17是n维向量,表示k时刻的状态向量
  • 【数据融合】状态估计理论 - 图18是k时刻 n*n 阶的状态转移矩阵
  • 【数据融合】状态估计理论 - 图19为(已知) p 维输入或控制信号(向量)
  • 【数据融合】状态估计理论 - 图20为 n*p 阶输入矩阵
  • 【数据融合】状态估计理论 - 图21为 n*n 维过程噪声分布矩阵
  • 【数据融合】状态估计理论 - 图22为 n 维过程噪声,满足Gauss白噪声,并且满足

【数据融合】状态估计理论 - 图23

  • 其中【数据融合】状态估计理论 - 图24为Dirichlet函数,即满足
  • 【数据融合】状态估计理论 - 图25
    • 量测方程满足

【数据融合】状态估计理论 - 图26

  • 【数据融合】状态估计理论 - 图27是 m 维向量,表示k时刻的量测向量
  • 【数据融合】状态估计理论 - 图28是 m*n 阶量测矩阵
  • 【数据融合】状态估计理论 - 图29是 m 维量测噪声,并满足Gauss白噪声,有

【数据融合】状态估计理论 - 图30

  • 初始状态描述
    • 假定初始状态【数据融合】状态估计理论 - 图31是高斯的,具有均值【数据融合】状态估计理论 - 图32和协方差【数据融合】状态估计理论 - 图33,且

【数据融合】状态估计理论 - 图34

  • 其中【数据融合】状态估计理论 - 图35

    2.3.2 Kalman推导

  • 【数据融合】状态估计理论 - 图36是已知 k 时刻及 k 以前时刻量测值,对【数据融合】状态估计理论 - 图37的最小均方误差估计
  • (称【数据融合】状态估计理论 - 图38【数据融合】状态估计理论 - 图39的滤波)

【数据融合】状态估计理论 - 图40

  • 相应的状态滤波值【数据融合】状态估计理论 - 图41的协方差矩阵为

【数据融合】状态估计理论 - 图42

  • 下面导出 k+1 时刻滤波值【数据融合】状态估计理论 - 图43和协方差矩阵 【数据融合】状态估计理论 - 图44
    • 【数据融合】状态估计理论 - 图45基础上增加了新的量测值【数据融合】状态估计理论 - 图46
    • 【数据融合】状态估计理论 - 图47是已知 k 时刻及 k 以前时刻量测值下,对 【数据融合】状态估计理论 - 图48 的最小均方误差估计
    • 状态预测值

【数据融合】状态估计理论 - 图49

  • 由 k+1 时刻状态预测值【数据融合】状态估计理论 - 图50得到 k+1 时刻的状态预测误差

【数据融合】状态估计理论 - 图51

  • k+1 时刻状态预测协方差【数据融合】状态估计理论 - 图52

【数据融合】状态估计理论 - 图53

  1. - 现需证明上式中第二、三项为0
  2. - 几个假设:
  3. - 1、测量噪声![](https://cdn.nlark.com/yuque/__latex/1f180697af0129a82a5435c4d58500e3.svg#card=math&code=W%28k%29&id=fhbKY)与状态值是独立的
  4. - 2、过程噪声![](https://cdn.nlark.com/yuque/__latex/99d7ae883ad9a5e00a7a9bbe38ddc521.svg#card=math&code=V%28k%29&id=S7QIm)与测量值是独立的
  5. - 3、过程噪声![](https://cdn.nlark.com/yuque/__latex/99d7ae883ad9a5e00a7a9bbe38ddc521.svg#card=math&code=V%28k%29&id=lCxsy)与测量噪声![](https://cdn.nlark.com/yuque/__latex/1f180697af0129a82a5435c4d58500e3.svg#card=math&code=W%28k%29&id=qTIUb)是独立的
  6. - 4k 时刻测量噪声![](https://cdn.nlark.com/yuque/__latex/1f180697af0129a82a5435c4d58500e3.svg#card=math&code=W%28k%29&id=uMKWW)与测量值![](https://cdn.nlark.com/yuque/__latex/23a05a27fca48c3bf21e1ab49abf173a.svg#card=math&code=Z%28k%29&id=zSwNY)是相关的
  7. - 5k 时刻的过程噪声![](https://cdn.nlark.com/yuque/__latex/99d7ae883ad9a5e00a7a9bbe38ddc521.svg#card=math&code=V%28k%29&id=bsjzt)与 k 时刻的状态值是独立的
  8. - 下一时刻噪声与当前量测值是独立的,估计量是个常数![](https://cdn.nlark.com/yuque/__latex/304fc1635c3e254b8fc689bfc1e67aa8.svg#card=math&code=%5Chat%20X%28k%2Fk%29&id=qLPQ1)
  9. - 第二项:

【数据融合】状态估计理论 - 图54

  1. - 第三项:

【数据融合】状态估计理论 - 图55

  • 从而有

【数据融合】状态估计理论 - 图56

  • 【数据融合】状态估计理论 - 图57是已知 k 及 k 以前时刻量测值下,对量测【数据融合】状态估计理论 - 图58作最小均方误差估计
  • 量测预测值

【数据融合】状态估计理论 - 图59


【数据融合】状态估计理论 - 图60


【数据融合】状态估计理论 - 图61

  1. - 上式第二项

【数据融合】状态估计理论 - 图62

  1. - 同理第三项也为0
  • 状态预测值与量测预测值协方差:

【数据融合】状态估计理论 - 图63

  • 量测预测值与状态预测值协方差:

【数据融合】状态估计理论 - 图64

  • 已知 k+1 时刻量测值【数据融合】状态估计理论 - 图65,导出相应 k+1 时刻 X 的滤波值【数据融合】状态估计理论 - 图66及滤波协方差【数据融合】状态估计理论 - 图67
  • 根据递推的最小均方误差估计公式

【数据融合】状态估计理论 - 图68

  1. - 需要注意
  2. - X 进行估计
  3. - Z 是新的量测值, ![](https://cdn.nlark.com/yuque/__latex/593d9b96ceba775ac243cb137dd63e4b.svg#card=math&code=Z_0&id=rxQjN)是以往量测值的集合
  4. - ![](https://cdn.nlark.com/yuque/__latex/58c3d0d0ebe063b8eb10b3dcf9c10b0c.svg#card=math&code=%5Cbar%7BX%7D&id=MFmm5)表示在![](https://cdn.nlark.com/yuque/__latex/593d9b96ceba775ac243cb137dd63e4b.svg#card=math&code=Z_0&id=uxmMq)条件下,对 X 的一个估计
  5. - ![](https://cdn.nlark.com/yuque/__latex/e09dfee7a4be42f8f9b465fda6debc80.svg#card=math&code=%5Cbar%7BZ%7D&id=JAFnV)是![](https://cdn.nlark.com/yuque/__latex/593d9b96ceba775ac243cb137dd63e4b.svg#card=math&code=Z_0&id=yyKBR)条件下,对 Z 的一个估计
  6. - ![](https://cdn.nlark.com/yuque/__latex/fd2b1eb3a3aacb801a8c6d0b7ec448b5.svg#card=math&code=%5Chat%7BX%7D&id=h5Vlc)是![](https://cdn.nlark.com/yuque/__latex/593d9b96ceba775ac243cb137dd63e4b.svg#card=math&code=Z_0&id=Jp2st)及 Z 条件下,对 X 的一个估计
  7. - ![](https://cdn.nlark.com/yuque/__latex/eb524a38ab1083099444bdbda2040829.svg#card=math&code=P_%7BXZ%7D&id=MbsfZ)是状态 X 与测量 Z 的协方差,而且是在![](https://cdn.nlark.com/yuque/__latex/593d9b96ceba775ac243cb137dd63e4b.svg#card=math&code=Z_0&id=r48HP)条件下

【数据融合】状态估计理论 - 图69

  1. - 这样,X 对应![](https://cdn.nlark.com/yuque/__latex/b548c1d047b1d568dc9fe7c5868b9320.svg#card=math&code=X%28k%2B1%29&id=xqdnp),![](https://cdn.nlark.com/yuque/__latex/593d9b96ceba775ac243cb137dd63e4b.svg#card=math&code=Z_0&id=P3utk)对应![](https://cdn.nlark.com/yuque/__latex/bc3ee99511099ae277a40bd113a4e8f6.svg#card=math&code=Z%5Ek&id=OUEMB),Z 对应![](https://cdn.nlark.com/yuque/__latex/ccb77b673da9e04e573d757cb7147bb3.svg#card=math&code=Z%28k%2B1%29&id=DCU6s),从而![](https://cdn.nlark.com/yuque/__latex/58c3d0d0ebe063b8eb10b3dcf9c10b0c.svg#card=math&code=%5Cbar%7BX%7D&id=P1gQd)对应![](https://cdn.nlark.com/yuque/__latex/c70502debadc608bb5c76155b9456079.svg#card=math&code=%5Chat%7BX%7D%28k%2B1%2Fk%29&id=puaCK),![](https://cdn.nlark.com/yuque/__latex/e09dfee7a4be42f8f9b465fda6debc80.svg#card=math&code=%5Cbar%7BZ%7D&id=zUGjS)对应![](https://cdn.nlark.com/yuque/__latex/461727cb1063272ea008ddc386cca156.svg#card=math&code=%5Chat%7BZ%7D%28k%2B1%2Fk%29&id=I8uOl),![](https://cdn.nlark.com/yuque/__latex/fd2b1eb3a3aacb801a8c6d0b7ec448b5.svg#card=math&code=%5Chat%7BX%7D&id=P0HvL)对应![](https://cdn.nlark.com/yuque/__latex/6e2e3ba9444269b7d1d20f31a7c89ad2.svg#card=math&code=%5Chat%7BX%7D%28k%2B1%2Fk%2B1%29&id=tM0GQ)
  2. - 则有

【数据融合】状态估计理论 - 图70

  1. - 利用

【数据融合】状态估计理论 - 图71

  1. - 可得

【数据融合】状态估计理论 - 图72

  • 定义增益

【数据融合】状态估计理论 - 图73

  • 定义新息

【数据融合】状态估计理论 - 图74

  • 增益的两种形式

【数据融合】状态估计理论 - 图75

  1. - 证明两种形式等价

【数据融合】状态估计理论 - 图76

  • 推导

【数据融合】状态估计理论 - 图77

  1. - 利用

【数据融合】状态估计理论 - 图78

  1. - 上式左乘![](https://cdn.nlark.com/yuque/__latex/a5b34f73a7d0cca60f80677e1ca68a56.svg#card=math&code=%7BP%5E%7B%20-%201%7D%7D%28k%20%2B%201%2Fk%20%2B%201%29&id=MA308),右乘![](https://cdn.nlark.com/yuque/__latex/22da61cc7e94c3d6950ddbb2da1c94b3.svg#card=math&code=%7BP%5E%7B%20-%201%7D%7D%28k%20%2B%201%2Fk%29&id=Rq0rY)得:

【数据融合】状态估计理论 - 图79