how to implement moving window

조회 수: 13 (최근 30일)
dab483
dab483 2013년 4월 17일
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
Jan
Jan 2013년 4월 17일
What did you try so far?
dab483
dab483 2013년 4월 18일
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)

댓글을 달려면 로그인하십시오.

답변 (0개)

카테고리

Help CenterFile Exchange에서 Online Estimation에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!

Translated by