how to implement moving window

13 Ansichten (letzte 30 Tage)
dab483
dab483 am 17 Apr. 2013
Hi, I have a question regarding moving window implementation. I have this equation. how to implement/coding it in Matlab where m = window size; k=epoch ?
C = 1/m * (summation(i=1 to m) of V(k-i)*V(k-i)')
where V is the residual between observation and measurement size 3x1
Any advise? Thank you.
  2 Kommentare
Jan
Jan am 17 Apr. 2013
What did you try so far?
dab483
dab483 am 18 Apr. 2013
Hi, Actually i am implementing extended kalman filter. i am extending it to estimate QekfA(:,:,t) where have moving window as shown below in AAA. i have no idea how i should continue this..Need help to continue this.
for t=2:T
mu_ekfPred(:,t) = feval('ffunasp',mu_ekf(:,t-1),t);
Ax_ekf = ....
PPred_ekf(:,:,t) = QekfA(:,:,t-1) + Ax_ekf*P_ekf(:,:,t-1)*Ax_ekf';
Cy_ekf = [0 0 1 0];
yPred_ekf(:,t) = Cy_ekf*mu_ekfPred(:,t);
M_ekf = Rekf + Cy_ekf*PPred_ekf(:,:,t)*Cy_ekf';
K_ekf = PPred_ekf(:,:,t)*Cy_ekf'*(M_ekf^(-1));
mu_ekf(:,t) = mu_ekfPred(:,t) + K_ekf*(y(:,t)-yPred_ekf(:,t));
P_ekf(:,:,t) =(eye(4)-K_ekf*Cy_ekf)*PPred_ekf(:,:,t)......
QekfA(:,:,t)=QekfA(:,:,t-1)*sqrt(AAA)
end
where AAA = 1/m * (summation(i=1 to m) of V(t-i)*V(t-i)') - Rekf
where V(t-i) = residual = y(:,t)-yPred_ekf(:,t)

Melden Sie sich an, um zu kommentieren.

Antworten (0)

Kategorien

Mehr zu Online Estimation finden Sie in Help Center und File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by