我如何阴谋粒子(粒子滤波)颜色梯度的方式与他们的重量吗?

2视图(30天)

显示旧的评论

达尼埃莱 2015年11月18日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/255886-how-can-i-plot-particles-of-a-particle-filter-with-their-weights-in-a-color-gradient-way

评论道: bjan912015年11月28日

我写一个粒子滤波,将绘制一个图形,我不知道怎么做,我疯狂的增长。

我绘制的颗粒,观察它的行为与原始状态的行为。基本上我有k粒子向量(轴的名字倒我认为)每个Ns粒子。我还k关联和归一化权重向量(右边的梯度栏)。

所以最后我绘制每个粒子的每一个向量与它的重量与颜色梯度像下面的图片。

任何帮助将不胜感激,并为我的英语感到惋惜。

1评论
显示没有隐藏没有

bjan91 2015年11月28日

嗨,丹尼尔,

“最后我将阴谋每个粒子的每一个向量与它的重量颜色梯度像下面图片里的。”

我只是工作在这个同样的问题(从Arulampalam看似相同的纸),我有一个解决方案,对我来说已经足够了,希望它会让你开始你的需要。下面我将分享的步骤和代码。

收集/保存在一个向量加权粒子每一步
排序向量,粒子是为了从最小最重要
定义一个colormap和linspace阴影的点按重量
情节在你使用散点图。

                               %操作重量和粒子的策划
                              
                               为我= 1:ltvec-1
                              
                               set2 (1: Ns, 1:2)中的= [update1(我:s)“P_w1(我:s) ');%(颗粒、重量)
                              
                               关于我校(1:Ns, 1:2) =(更新2(我:s)“P_w2(我:s) ');
                              
                               一个= sortrows (set2(:,:我)中的2);sorted1(:,我,s) = (: 1);
                              
                               B = sortrows((:,:我),2);sorted2(:,我,s) = B (: 1);
                              
                               % sorted1列当前粒子的步伐,
                              
                               根据重量%排序和命令((至少,.....;大多数));
                              
                               结束
                              
                               …. .
                              
                               …
                              
                               …. .
                              
                               情节(ttvec x1_set (:, j),“b -”);持有在;情节(ttvec x1hat_interp (:, j),“r——”);
                              
                               colormap(灰色)
                              
                               c = linspace(1、10、大小(sorted1 1));
                              
                               为我= 1:ltvec-1
                              
                               的散射(tvec (i + 1) * (Ns, 1), sorted1 (:, i, j), 15日,c);
                              
                               结束