drawpar.asv

基于SIR粒子滤波器的故障诊断方法

function xnext = drawpar(xprev, stdw, k)
M = length(xprev);
b0=25;
for i=1:M
    %if i<301
     %   b=b0;
    %else
     %   b=b0*0.5;
    %end
xstate = xprev(i)/2 + b0*xprev(i)/(1+xprev(i)^2) +8*cos(1.2*(k-1));
xnext(i,1) = stdw*randn + xstate;
end