drawpar.m

基于粒子滤波的故障检测

function xnext = drawpar(xprev, stdw, k)
M = length(xprev);

for i=1:M
   xstate = xprev(i)/2 + 25*xprev(i)/(1+xprev(i)^2) +8*cos(1.2*(k-1));
xnext(i,1) = stdw*randn+xstate;
end