[发明专利]一种基于直方图估计粒子滤波算法的非线性动态系统信号处理方法

说明书

技术领域

本发明涉及基于粒子滤波算法的信号处理方法，要求保护的技术方案属于信号处理、人工智能和计算机视觉领域。

背景技术

动态系统的状态估计问题涉及很多领域，尤其是信号处理、人工智能和计算机视觉领域。传统的卡尔曼滤波只适用于线性高斯系统，而扩展卡尔曼滤波也只能应对系统的弱非线性。因此，适用于非线性、非高斯系统的粒子滤波备受关注。

粒子滤波是一种基于蒙特卡洛模拟和递推贝叶斯估计的滤波方法。它采用粒子描述状态空间，使用一组带权重的粒子近似表示系统的后验概率密度，并通过模型方程和观测信息实现递推的估计过程。常见的粒子滤波算法包括SIR粒子滤波、辅助粒子滤波(APF)、正则粒子滤波(RPF)、高斯粒子滤波(GPF)和无迹粒子滤波(UPF)。

粒子滤波算法中的两个关键技术是建议分布的选取和重采样算法。常用的建议分布虽然容易获取，但易造成粒子权重的急剧退化，使得有效粒子数大大减少，尤其是当似然分布位于先验分布的尾部时。同时，虽然重采样算法解决了粒子权重的退化问题，但常见的重采样算法只是对粒子的简单复制和剔除，从而又导致了粒子匮乏问题(粒子多样性的丧失)。

为了克服现有基于粒子滤波算法的动态系统信号处理方法中的由于粒子权重易急剧退化和粒子匮乏，不能有效代表后验密度，因而无法有效处理非线性和非高斯信号的不足，本发明提出一种能更好地估计系统的后验概率密度、粒子多样性好、能有效处理非线性和非高斯信号的基于直方图估计粒子滤波算法的非线性动态系统信号处理方法。

为了解决上述技术问题提出的技术方案为：

一种基于直方图估计粒子滤波算法的非线性动态系统信号处理方法，设非线性动态系统的状态空间模型为：

x_k＝f(x_k-1)+v_k-1

z_k＝h(x_k)+n_k

其中，x_k和z_k分别表示系统在k时刻的状态和观测值，f(x_k-1)和h(x_k)分别表示系统的状态转移方程和观测方程，x_k-1表示系统在k时刻的状态，x_k-1和n_k分别表示系统噪声和观测噪声；

所述非线性动态系统信号处理方法包括以下步骤：

第一步，根据k-1时刻的N个粒子i＝1，2，...，N，通过状态转移方程得到k时刻的N个预测粒子i＝1，2，...，N；

第二步，确定预测粒子所占据的状态空间，并将其均分为箱，箱的个数记为M；

第三步，根据划分的箱估计系统的先验概率密度，各箱对应的密度d_i为：

di=NiN×h,i=1,2,...,M]]>

其中，Nⁱ表示落入对应箱的预测粒子数，h表示一个箱所占据的超空间的体积。

第四步，估计各箱对应的后验概率密度wⁱ：

wⁱ＝dⁱ×p(z_k|xⁱ)，i＝1，2，...，M

其中，xⁱ表示第i个箱的中心；

第五步，对各箱的后验概率密度进行归一化：