Changeset 1090 for applications/dual/SIDP/text/ch4.tex
- Timestamp:
- 06/13/10 17:27:15 (14 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
applications/dual/SIDP/text/ch4.tex
r930 r1090 1 V t� kapitole je pops�jednoduch�� na kter�jsou porovn� ��lgoritmy uveden� p�l�apitole. Syst�byl podrobn�koum�v \cite{astrom1986dual}. Pro srovn� uv�me tam���y.1 V t� kapitole je pops�jednoduch��zkouman�ite{astrom1986dual}. Na n�jsou porovn� ��lgoritmy uveden� p�l�apitole. 2 2 3 3 \section{Popis syst�} … … 5 5 \begin{gather} 6 6 \label{simple} 7 y_{t+1}=y_t+\theta_tu_t+v_{t+1} \qquad t=0,\ldots,N-1,\\ 8 v_t\sim N(0,\sigma^2).\\ 9 \theta_t\sim N(\hat{\theta},P_t),\\ 7 y_{t+1}=y_t+\theta u_t+v_{t+1} \qquad t=0,\ldots,N-1,\\ 8 v_{t+1}\sim N(0,\sigma^2), 9 \end{gather} 10 kde rozptyl �umu $\sigma$ je zn� 11 12 O nezn�m parametru $\theta$ m� v �e $t$ informaci v podob�ostate� statistiky $T_t=(\hat{\theta},P_t)$, tvo�st� hodnotou a rozptylem. P�kl�me nekorelovanost $\theta$ s �umem, tedy �e 13 \begin{equation} 10 14 \cov(v_{t+1},\theta)=0. 11 \end{ gather}15 \end{equation} 12 16 13 17 Ztr�vou funkci vol� kvadratickou, tedy 14 18 \begin{equation} 15 g(y_{0:N},u_{0:N-1} ,v_{0:N-1})=\sum_{t=0}^{N-1}y_{t+1}^2.19 g(y_{0:N},u_{0:N-1})=\sum_{t=0}^{N-1}y_{t+1}^2. 16 20 \end{equation} 17 21 … … 24 28 \end{gather} 25 29 26 O��n�tr� je30 Hyperstav syst� $H_t$ tvo�ktor $(y_t,\hat{\theta}_t,P_t)$. O��n�tr� je 27 31 \begin{equation} 28 J_t( y_t,\theta_t)=\min_{u_t \in U_t}\E_{y_{t+1},v_t}\left\{y_{t+1}^2+J_{t+1}(y_{t+1},\theta_{t+1})|y_t,\theta_t,u_t\right\}, \qquad t=0,\ldots,N-1.32 J_t(H_t)=\min_{u_t \in U_t}\E_{y_{t+1},v_t}\left\{y_{t+1}^2+J_{t+1}(H_{t+1})|H_t,u_t\right\}, \qquad t=0,\ldots,N-1. 29 33 \end{equation} 30 34 … … 35 39 \end{gather} 36 40 37 ZDE BY MEL BYT ANGSTROM+...38 39 41 \section{Specifika jednotliv��up� tomto odd� jsou pops� n�er�spekty algoritm�er�udeme srovn�t, p�likaci na syst�\eqref{simple}. 40 42 … … 42 44 O��n�tr� \eqref{CE} prejde v 43 45 \begin{gather} 44 J_t( y_t, \theta_t)=\min_{u_t \in U_t}\left\{\hat{y}_{t+1}^2 +J_{t+1}(y_{t+1},\theta_{t+1})|I_t,\theta_t,u_t\right\}.46 J_t(H_t)=\min_{u_t \in U_t}\left\{\hat{y}_t^2 +J_{t+1}(\hat{H}_{t+1})|I_t,\theta_t,u_t\right\}. 45 47 \end{gather} 46 48 St� hodnota v� je