Show
Ignore:
Timestamp:
05/09/10 22:10:09 (14 years ago)
Author:
zimamiro
Message:
 
Files:
1 modified

Legend:

Unmodified
Added
Removed
  • applications/dual/SIDP/text/ch2.tex

    r919 r930  
    5252 
    5353\section{�zen�yst� s nezn�mi parametry} 
    54 Pokud chceme � syst� jeho� v�z�s�a n�k�nezn�m konstant�parametru $\theta$, m� vyu��znalosti ��robl� s ne�m pozorov�m.  
     54Pokud chceme � syst� jeho� v�z�s�a n�k�nezn�m konstant�parametru $\theta$, m� vyu��znalosti ��robl� s ne�m pozorov�m. Parametr $\theta$ bude reprezentovat stav syst� $x_t$, kter�yn� �e nem�. 
    5555 
    5656V t� � m� v� syst� $y_t$ pops� jako  
    5757\begin{equation} 
    5858\label{poz2} 
    59 y_0=h_0(\theta,v_0),\qquad y_t=h_t( I_{t-1},\theta,u_{t-1},v_t), \qquad t=1,\ldots,N-1, 
     59y_0=h_0(\theta,v_0),\qquad y_t=h_t(I_{t-1},\theta,u_{t-1},v_t), \qquad t=1,\ldots,N-1, 
    6060\end{equation} 
    6161 
     
    6565\end{equation} 
    6666 
    67 Ozna� $\theta_t$ rozlo�en�dhadu pro parametr $\theta$ v �e $t$. P�kl�jme d�, �e o parametru $\theta$ m� n�kou apriorn�nformaci v podob�ustoty pravd�dobnosti $f(\theta_0)$. Aposteriorn�ustotu  $f(\theta_{t+1}|I_t)$ z�� pomoc�ayesova vzorce 
     67Ozna� $T_t$ testovac�tatistiku pro parametr $\theta$ zalo�enou na informac� dostupn� �e $t$. Do $T_t$ zahrneme rovnez ty cleny $I_t$, kter�ystupuj� \eqref{poz2}, abychom mohli ps� 
     68\begin{equation} 
     69\label{poz3} 
     70y_{t+1}=h_t(T_t,\theta,u_t,v_{t+1}). 
     71\end{equation} 
     72 
     73P�kl�jme d�, �e o parametru $\theta$ m� n�kou apriorn�nformaci v podob�ustoty pravd�dobnosti $f(\theta|T_0)$. Aposteriorn�ustotu $f(\theta|T_{t+1})$ z�� pomoc�ayesova vzorce 
    6874\begin{equation} 
    6975\label{bay} 
    70 f(\theta_{t+1}|I_t)=\frac{f(I_t|\theta_{t+1})f(\theta_t)}{\int f(I_t|\theta_{t+1})f(\theta_t)\mathrm{d}\theta_t} 
     76f(\theta|T_{t+1})=\frac{f(T_{t+1}|\theta,T_t)f(\theta|T_t)}{\int f(T_{t+1}|\theta,T_t)f(\theta|T_t)\mathrm{d}\theta} 
    7177\end{equation} 
    72 Rekurzivn�ou�it�zorce \eqref{bay} pro odhad parametru $\theta$ je postup Bayesovsk� u��ref]. 
     78Rekurzivn�ou�it�zorce \eqref{bay} pro odhad parametru $\theta$ je postup Bayesovsk� u��cite{peterka1981bayesian}. 
    7379 
    74 Ozna� $T_t$ doste�u statistiku pro $\theta_t$. Potom dle \eqref{bay} m� ps� 
     80Pro v�estovac�tatistiky v �e m� podle \eqref{bay} ps� 
    7581\begin{equation} 
    7682\label{the} 
    77 T_{t+1}=f_t(I_t,T_t,u_t,y_{t+1}), \qquad  t=1,\ldots,N-1. 
     83T_{t+1}=f_t(T_t,u_t,y_{t+1}), \qquad  t=1,\ldots,N-1. 
    7884\end{equation}  
    79 Rovnici \eqref{the} m� podobn�ako \eqref{nep} pova�ovat za rovnici syst� \eqref{sys} pro stav $(I_t, T_t)$, vstup $u_t$ s �umem $y_{t+1}$.  
     85Rovnici \eqref{the} m� podobn�ako \eqref{nep} pova�ovat za rovnici syst� \eqref{sys} pro stav $T_t$ a vstup $u_t$ s �umem $y_{t+1}$.  
    8086 
    81 Hustota pravd�dobnosti parametru $\theta$ v rovnici pro v�\eqref{poz2} je v �e $t$ ur�a dostate�u statistikou $T_t$. Rovnice \eqref{the}, \eqref{poz2} a \eqref{los2} potom p�avuj�lohu stochastick� �� nep�mi daty. 
     87Hustota pravd�dobnosti pro odhad parametru $\theta$ v rovnici pro v�\eqref{poz3} je v �e $t$ ur�a testovac�tatistikou $T_t$. Rovnice \eqref{the}, \eqref{poz3} a \eqref{los2} potom p�avuj�lohu stochastick� �� nep�mi daty. 
    8288 
    8389\subsection{Kalman�ltr}