[R-br] Dúvida em erro na função Filtro de Kalman

[Victor Eduardo]

Ivan, no for é  *for(t in 1:n)*, sem querer coloquei 2 ali. No caso o erro
é nas dimensões da matriz? Só que o ar_A está sendo multiplicado, e ele vai
ficar na mesma dimensão de D. No caso o array ar_A tem n+1 e o os outros
tem n, só que seria a previsão.



Abraços!




Em 14 de maio de 2012 09:39, Ivan Bezerra Allaman
<ivanalaman em yahoo.com.br>escreveu:

> O erro é claro!
>
> Verifique a linha que eu grifei em vermelho no seu código. As dimensão do
> objeto "ar_A e ar_D" estão diferentes, impossibilitando a subtração.
>
> .
> .
> .
> for(t in 1:2){
> ar_U[,,t]<-ar_Y[,,t]-ar_Z[,,t,drop=FALSE]%*%
> ar_A[,,t,drop=FALSE]-ar_D[,,t,drop=FALSE]
>
> ar_F[,,t]<-ar_Z[,,t,drop=FALSE]%*%ar_P[,,t,drop=FALSE]%*%t(ar_Z[,,t,drop=FALSE])+ar_H[,,t,drop=FALSE]
>
> ar_AA[,,t]<-ar_A[,,t,drop=FALSE]+ar_P[,,t,drop=FALSE]%*%ar_Z[,,t,drop=FALSE]%*%ginv(ar_F[,,t,drop=FALSE])%*%ar_U[,,t,drop=FALSE]
>
> ar_PP[,,t]<-ar_P[,,t,drop=FALSE]-ar_P[,,t,drop=FALSE]%*%t(ar_Z[,,t,drop=FALSE])%*%ginv(ar_F[,,t,drop=FALSE])%*%ar_Z[,,t,drop=FALSE]%*%ar_P[,,t]
>
> ar_A[,,t+1]<-ar_T[,,t,drop=FALSE]%*%ar_AA[,,t,drop=FALSE]+ar_C[,,t,drop=FALSE]
>
> ar_P[,,t+1]<-ar_T[,,t,drop=FALSE]%*%ar_PP[,,t,drop=FALSE]%*%t(ar_PP[,,t,drop=FALSE])+ar_R[,,t,drop=FALSE]%*%ar_Q[,,t,drop=FALSE]%*%t(ar_R[,,t,drop=FALSE])
>  }
> return(list(ar_A,ar_P))
> }
>
> Filtro_Kalman(ar_Z,ar_T,ar_C,ar_D,ar_R,ar_H,ar_Q,a1,p1)
>
> (S,f,P)
> Allaman
>
> *
> *
> \begin{signature}
> <<>>=
> Prof. Dr. Ivan Bezerra Allaman
> Universidade Estadual de Santa Cruz
> Departamento de Ciências Exatas e Tecnológicas
> Ilhéus/BA - Brasil
> Fone: +55 73 3680-5596
> E-mail: ivanalaman em yahoo.com.br/ivanalaman em gmail.com
> @
> \end{signature}
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120514/a8bf2e9c/attachment.html>