Administrator
发布于 2026-09-09 / 0 阅读
0
0

粒子滤波

粒子滤波所要做的就是估计感兴趣的信号,与传统的估计方法不同,它主要针对非线性系统或者乘性噪音这类传统估计无法解决的问题,其核心方法是蒙特卡洛法,它基于随机信号的概率建模。因此我们先回顾概率统计部分的知识。

随机变量

x 表示随机变量,用 \Omega 表示概率空间

离散随机变量

离散随机变量从有限可数的集合 \Omega = \{x_i \mid i \in I\subset \mathbb{Z}\} 中取值,它的统计特性由概率分布给出。

我们用 p_i 来表示取到第几个 x_i ,即:

p_i = P(X = x_i)

每个离散元素概率取值和等于1,即:

\sum_{i \in I} p_i = 1

数学期望

E[X] = \sum_{i \in I} x_i p_i

方差

\sigma^2 = E[(X - E[X])^2] = E[X^2] - (E[X])^2

连续随机变量

我们有一个连续的概率密度函数,在某点 x0 的概率密度值定义为:

p(x_0) = \frac{P(x_0 < X < x_0 + dx)}{dx}

其统计特性由概率密度函数 p(x) 来刻画,具体表现为:

  • 概率密度函数值一定为正
p(x) ≥ 0
  • 区间概率等于区间积分
P[X ∈ [a,b]] = ∫_a^b p(x) dx
  • 概率密度函数的总积分为 1
∫_{−∞}^{+∞} p(x) dx = 1
  • 分布函数 F(x) 定义为
F(x) = P[X ≤ x] = ∫_{−∞}^x p(u) du

数学期望

E[X] = \int x \, p(x) \, dx

方差

\sigma^2 = E[(X - E[X])^2] = E[X^2] - (E[X])^2

对于高斯概率密度函数分布,约有95%的变量取值落在 E[X] 的两侧 2\sigma 的范围内,即 [E[X]-2\sigma, E[X]+2\sigma] 区间之中。

现在我们分别探讨完了连续和离散随机变量的公式,我们想做的是,用一种相同的形式来同时表示离散和连续随机变量的分布律

假设我们有一个离散随机变量 X,其可能取值为 x_1, x_2, ..., x_n,概率简化为 p_i = P[X = x_i]。则它的离散的概率分布可以表示为:

p(x) = ∑_{i∈I} p_i δ(x - x_i)

我们下面给出证明

补充知识,狄拉克 δ 函数的特性是:

∫ f(x) δ(x - x_i) dx = f(x_i)

特别地,当我们对某一区间 J 积分时,有:

∫_J δ(x - x_i) dx= \begin{cases} 1 & \text{if } x_i \in J \\ 0 & \text{if } x_i \notin J \end{cases}

这一性质表明,狄拉克函数可以挑选出 x_i

推导过程:

对于离散随机变量 X,如果 \{x_1, x_2, ..., x_n\}X 的所有可能取值,那么

P(X = x_i) = p_i

对于任意集合(或区间) JX 落入其中的概率为 在这个集合中 x_i 概率的累加

P(X \in J) = \sum_{x_i \in J} p_i

在连续型分布中,概率密度函数 p(x) 满足:

P[X ∈ J] = ∫_J p(x) dx

我们想用连续型分布公式来类比离散情况的形式,因此引入狄拉克函数

P[X ∈ J] = ∫_J (∑_{i=1}^n p_i δ(x - x_i)) dx

![image-20241209160158852](/Users/zehua/Library/Application Support/typora-user-images/image-20241209160158852.png)

未完成

对于离散变量,我们对其区间取一部分,然后计算 x 落在这个区间 J 的概率为:

P[x \in J] = \sum_{x_i \in J} P(x_i) =P(X =x_1)+P(X =x_2) + ...= \sum P_i

然后我们用连续方式来写出区间概率公式:

P[X ∈ J] = ∫_J p(x) dx

将上面的离散分布形式代入

P[X ∈ J] = ∫_J (∑ p_i δ(x - x_i)) dx

由于概率是具体的数值,所以可以提到积分外,交换顺序:

P[X ∈ J] = ∑ p_i ∫_J δ(x - x_i) dx

∫_J δ(x - x_i) dx 的值取决于 x_i 是否在集合 J 中:

∫_J δ(x - x_i) dx= \begin{cases} 1 & \text{if } x_i \in J \\ 0 & \text{if } x_i \notin J \end{cases}

因此:

P[X ∈ J] = ∑_{x_i ∈ J} p_iδ(x - x_i)

这与离散随机变量的基本定义完全一致。当 X 的取值在集合 J 内时,其概率就是落在该集合内的点的所有概率质量之和。


评论