plotting RMSE

Hi,

i am calculating root mean square error (RMSE) of true state(x) and estimated state (xhat)using the equation below:

err=sqrt(sum(x(:,1)-xhat(:,1))^2/numel(x))

this will gave me one answer at the last of simulation time. Now, if i want to plot the RMSE over simulation time, how should i coding it?

Answer 1

[bym]

Assuming you are calculating x and xhat in a loop with an index variable(e.g. k), then

err(k)=sqrt(sum(x(1:k,1)-xhat(1:k,1))^2/k);
t(k) = dt*k;  %simulation time, where dt is time for one increment of loop
%outside of loop:

plot(err,t)