当前位置:   article > 正文

各种滤波算法的比较(GF、KF、EKF、UKF、PF),内附简单实现代码_kf和ukf

kf和ukf

目录

一、前言

二、滤波算法介绍

1、GF(高斯滤波)

2、KF(卡尔曼滤波)

3、EKF(可扩展卡尔曼滤波)

4、UKF(无迹卡尔曼滤波)

5、PF(粒子滤波)

三、不同滤波算法对比

四、不同滤波算法的表现

五、简单实现


一、前言

为了使基于RSSI/CSI等室内定位的结果更加地稳定,让滤波后的RSSI/CSI值更接近真实值,针对不同场景引入合适的滤波算法是很有必要的。当然文章并非只针对室内定位这个领域,其它领域,特别是信号或者通信领域,了解各种滤波算法也是非常重要的!

二、滤波算法介绍

1、GF(高斯滤波

高斯滤波是一种根据高斯函数的形状来选择权值的线性平滑滤波器,它对抑制服从正态分布的噪声非常有效,从而达到平滑数据的目的。

2、KF(卡尔曼滤波)

卡尔曼滤波算法是一种抑制高斯噪声有效的最优化自回归数据处理算法,对于系统过程的滤波比较好,能够有效滤除信号发生的突变。并且对于变化快速、实时更新的线性系统有着非常好的寻优及滤波效果。

卡尔曼滤波原理的详细解释可参考下面文章:

CSI数据预处理之卡尔曼滤波、高斯滤波、简单平均-CSDN博客

3、EKF(可扩展卡尔曼滤波

扩展卡尔曼滤波器是对KF的改进衍生,它将非线性函数进行泰勒展开,然后省略高阶项,保留展开项的一阶项,以此来实现非线性函数线性化,最后通过卡尔曼滤波算法近似计算系统的状态估计值和方差估计值。

4、UKF(无迹卡尔曼滤波

无迹卡尔曼滤波算法也是对KF的改进衍生,它舍弃了EKF算法的线性化过程,采用无迹变换(Unscented Transform,UT)巧妙地避免了线性化所带来的误差,同时减少了算法的复杂度UKF不需要在每个实例中计算雅可比矩阵),克服了EKF算法的估计精度低、稳定性差的缺陷。

5、PF(粒子滤波

粒子滤波法是指通过寻找一组在状态空间传播的随机样本对概率密度函数进行近似,以样本均值代替积分运算,从而获得状态最小方差分布的过程。和其它非线性滤波方法相同,它也是一种次优的滤波方法

三、不同滤波算法对比

四、不同滤波算法的表现

五、简单实现

本文就是做一个简单实现,其中高斯滤波和卡尔曼滤波比较基础,就不再展示。首先,定义一个简单的非线性状态空间模型,例如一个带有噪声的一维运动模型:

  1. % 非线性状态空间模型
  2. A = 1; % 状态转移矩阵
  3. B = 0; % 控制输入矩阵
  4. H = 1; % 观测矩阵
  5. % 状态转移方程
  6. f = @(x) A*x + B*randn;
  7. % 观测方程
  8. h = @(x) H*x + randn;

可扩展卡尔曼滤波:

  1. % 初始状态和协方差
  2. x0 = 0;
  3. P0 = 1;
  4. % EKF算法
  5. x_EKF = zeros(1, 100); % 存储估计状态
  6. P_EKF = zeros(1, 100); % 存储估计协方差
  7. x_pred = x0;
  8. P_pred = P0;
  9. for k = 1:100
  10. % 预测步骤
  11. x_pred = f(x_pred);
  12. P_pred = A*P_pred*A';
  13. % 更新步骤
  14. K = P_pred*H'/(H*P_pred*H' + 1);
  15. z = h(x_pred);
  16. x_pred = x_pred + K*(z - H*x_pred);
  17. P_pred = (1 - K*H)*P_pred;
  18. x_EKF(k) = x_pred;
  19. P_EKF(k) = P_pred;
  20. end

无迹卡尔曼滤波:

  1. % UKF算法
  2. x_UKF = zeros(1, 100); % 存储估计状态
  3. P_UKF = zeros(1, 100); % 存储估计协方差
  4. % 参数设置
  5. alpha = 1e-3;
  6. beta = 2;
  7. kappa = 0;
  8. % 初始化Sigma点
  9. n = numel(x0);
  10. lambda = alpha^2 * (n + kappa) - n;
  11. weights = 1 / (2 * (n + lambda)) * ones(1, 2*n+1);
  12. weights(1) = lambda / (n + lambda);
  13. X = zeros(n, 2*n+1);
  14. X(:,1) = x0;
  15. P_sqrt = chol(P0, 'lower');
  16. for i = 1:n
  17. X(:,i+1) = x0 + sqrt(n + lambda) * P_sqrt(:,i);
  18. X(:,n+i+1) = x0 - sqrt(n + lambda) * P_sqrt(:,i);
  19. end
  20. for k = 1:100
  21. % 预测步骤
  22. X_pred = zeros(size(X));
  23. for i = 1:2*n+1
  24. X_pred(:,i) = f(X(:,i));
  25. end
  26. x_pred = X_pred * weights';
  27. P_pred = zeros(size(P0));
  28. for i = 1:2*n+1
  29. P_pred = P_pred + weights(i) * (X_pred(:,i) - x_pred) * (X_pred(:,i) - x_pred)';
  30. end
  31. % 更新步骤
  32. Z = zeros(1, 2*n+1);
  33. for i = 1:2*n+1
  34. Z(i) = h(X_pred(:,i));
  35. end
  36. z_pred = Z * weights';
  37. Pzz = zeros(size(H*P_pred*H'));
  38. Pxz = zeros(size(P_pred));
  39. for i = 1:2*n+1
  40. Pzz = Pzz + weights(i) * (Z(i) - z_pred) * (Z(i) - z_pred)';
  41. Pxz = Pxz + weights(i) * (X_pred(:,i) - x_pred) * (Z(i) - z_pred)';
  42. end
  43. K = Pxz / Pzz;
  44. x_pred = x_pred + K * (h(x_pred) - z_pred);
  45. P_pred = P_pred - K * Pzz * K';
  46. x_UKF(k) = x_pred;
  47. P_UKF(k) = P_pred;
  48. end

 粒子滤波:

  1. % PF算法
  2. x_PF = zeros(1, 100); % 存储估计状态
  3. P_PF = zeros(1, 100); % 存储估计协方差
  4. % 参数设置
  5. num_particles = 1000;
  6. % 初始化粒子
  7. particles = randn(1, num_particles);
  8. for k = 1:100
  9. % 预测步骤
  10. particles = f(particles);
  11. % 更新步骤
  12. weights_PF = exp(-(h(particles) - 1).^2 / (2 * 1^2)); % 高斯权重函数
  13. % 归一化权重
  14. weights_PF = weights_PF / sum(weights_PF);
  15. % 重采样
  16. indices = randsample(1:num_particles, num_particles, 'true', weights_PF);
  17. particles = particles(indices);
  18. % 估计状态和协方差
  19. x_PF(k) = mean(particles);
  20. P_PF(k) = var(particles);
  21. end

生成符合高斯分布的状态和观测值:

  1. % 生成符合高斯分布的状态和观测值
  2. true_states = zeros(1, 100);
  3. observations = zeros(1, 100);
  4. true_states(1) = randn(); % 初始状态符合高斯分布
  5. observations(1) = true_states(1) + 0.5*randn(); % 初始观测值加入噪声
  6. for k = 2:100
  7. % 状态转移模型
  8. true_states(k) = 0.9 * true_states(k-1) + 0.1 * randn();
  9. % 观测模型
  10. observations(k) = true_states(k) + 0.5*randn(); % 添加观测噪声
  11. end

三种滤波算法可视化:

  1. % 画图
  2. figure;
  3. subplot(3,1,1);
  4. plot(1:100, true_states, 'LineWidth', 1.2, 'DisplayName', 'True State');
  5. hold on;
  6. plot(1:100, x_EKF, 'LineWidth', 1.2, 'DisplayName', 'EKF');
  7. title('Extended Kalman Filter (EKF)');
  8. xlabel('Time step');
  9. ylabel('State');
  10. legend;
  11. subplot(3,1,2);
  12. plot(1:100, true_states, 'LineWidth', 1.2, 'DisplayName', 'True State');
  13. hold on;
  14. plot(1:100, x_UKF, 'LineWidth', 2, 'DisplayName', 'UKF');
  15. title('Unscented Kalman Filter (UKF)');
  16. xlabel('Time step');
  17. ylabel('State');
  18. legend;
  19. subplot(3,1,3);
  20. plot(1:100, true_states, 'LineWidth', 1.2, 'DisplayName', 'True State');
  21. hold on;
  22. plot(1:100, x_PF, 'LineWidth', 2, 'DisplayName', 'PF');
  23. title('Particle Filter (PF)');
  24. xlabel('Time step');
  25. ylabel('State');
  26. legend;

 效果:

博主的每篇博文都是用心去写的,喜欢的可以多多支持和收藏,创作不易,未经作者允许,请勿转载或者抄袭。因为抄袭风气盛行,故一些细节或者代码没有展示,敬请谅解,如有疑问,可以加入我们的室内定位大家庭!

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/从前慢现在也慢/article/detail/744801
推荐阅读
相关标签
  

闽ICP备14008679号