rsmp.m

来自「matlab环境下粒子滤波器demo」· M 代码 · 共 17 行

function Inew=rsmp(pik,N)
%picks N indeces from pik randomly. The probability of a certain element is proportional to the size of the corresponding element in pik.
%Author PSL 200411
akkupik=cumsum(pik);

Inew=zeros(N,1);
y=rand(1)/N;
j=1;
for i=1:N        
    while(akkupik(j)<y)%flush all values that are below the current level of integrated propability
        j=j+1;        
    end
    Inew(i)=j;
    y=y+1/N;
end