Info
This question is closed. Reopen it to edit or answer.
How to use Wald Sequence detection method to test the filtering residual?
1 view (last 30 days)
Show older comments
Hi,
I have this EKF Equation.
PMinus(:,:,i) = A(:,:,i)*P(:,:,i)*A(:,:,i)'+ Q(:,:,i);
xhat(:,:,i+1) = Fs(xhat(:,:,i),u(:,:,i),wp(:,:,i),dt);
K(:,:,i) = PMinus(:,:,i)*H'*(H*PMinus(:,:,i)*H'+R)^-1;
xhat(:,:,i+1) = xhat(:,:,i+1)+K(:,:,i)*(y(:,:,i+1)-H*xhat(:,:,i+1));
P(:,:,i+1) = (eye(4)-K(:,:,i)*H)*P(:,:,i);
Now I need to use Wald Sequence detection method to test the filtering residual, r=y-H*xhat. to observe if the residual contain any fault information where its distribution N(m,sg^2) The general simplified decision law of wald sequence is:
Ly(n-1)=Ly(n) + m/sg^2*[r(n)-0.5m]
Any one can help me to code it to use with my above equation? Thanks.
0 Comments
Answers (0)
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!