root/applications/dual/SIDP/text/ch1.tex @ 1090

Revision 1090, 3.8 kB (checked in by zimamiro, 14 years ago)
Line 
1DEFINICNI OBORY
2\section{Formulace � stochastick� ��
3�t�m pojmem v teorii ��e \emph{syst�. Syst�je �t sv�, kterou chceme poznat ��. Budeme-li p�kl�t diskr��ovahu �u, stav syst� v �ov�okam�iku $t$ pod��� horizontu d�y $N$ popisuje syst�rovnic
4\begin{equation}
5\label{sys}
6x_{t+1}=f_k(x_t,u_t,w_t), \qquad t=0,1,\ldots,N-1,
7\end{equation}
8kde $x_t$ je stav syst� v �e $t$, $u_t$ je vstup v �e $t$ a $w_t$ n�dn�eli�a reprezentuj� p�nost �umu. V t� kapitole budeme p�kl�t, �e m� stav syst� pozorovat. P�em ne�ho pozorov� se zab�sleduj� kapitola.
9
10V � ��� v�dy p�sanou ztr�vou (resp. �vou) funkci
11\begin{equation}
12g(x_{1:N},u_{0:N-1}).
13\end{equation}
14
15Ozna� $U(x_t)$ mno�inu p�tn�d�ch z�h� syst�ve stavu $x_t$. P�tnou ��trategi�\pi=\mu_{0:N-1}$ budeme rozum�posloupnost zobrazen�\begin{equation}
16\label{con}
17\mu_t(x_t)=u_t \, \qquad t=0,1,\ldots,N-1,
18\end{equation}
19kde $u_t \in U(x_t)$ je p�tn�c��h.
20
21Pro danou ��trategii ozna� o��nou ztr� jako
22\begin{equation}
23\label{los}
24J_\pi(x_0)=\E_{w_{0:N-1}}\left\{g(x_{1:N},\mu_{0:N-1}(x_{0:N-1}))\right\}.
25\end{equation}
26
27�ohou je potom naj�takovou $\pi^*$, pro kterou plat�\begin{equation}
28J_{\pi^*}(x_0)=\min_{\pi \in \Pi}J_\pi(x_0),
29\end{equation}
30kde $\Pi$ zna�mno�inu v�ech p�tn�d�ch strategi�
31
32Celkov�e tedy jedn� optimaliza� � nal� takovou posloupnost funkc�eqref{con}, kter�inimalizuje o��nou ztr�vu \eqref{los} za podm�k \eqref{sys}.
33
34\section{Pou�it�ynamick� programov� p��en�lohy stochastick� �� aditivn�tr�u}
35�ohu stochastick� ��ak, jak byla definov� v p�oz��i, nelze obecn�e�it. Je tedy pot�� n�k bl� specifikovat. V tomto sm� je mo�n�mezit se na n�k�i��var ztr�v�unkce \eqref{los}. Jako vhodn�e ukazuje uva�ovat tzv. aditivn�var ztr�v�unkce, tedy �e existuj�unkce $g_t$ takov��e m� ps�
36\begin{equation}
37\label{adi}
38g(x_{1:N},u_{0:N-1})=\sum_{t=1}^{N-1}g_t(x_{t+1},u_t).
39\end{equation}
40
41O��nou ztr� \eqref{los} potom m� p�t do tvaru
42\begin{equation}
43J_\pi(x_0)=\E_{w_{0:N-1}}\left\{\sum_{t=0}^{N-1}g_t(x_{t+1},\mu_t(x_t))\right\}.
44\end{equation}
45
46Takto specifikovan�loha se d�e�it pou�it�dynamick� programov� \cite{bellman1957dynamic}. Dynamick�rogramov� je p�p k ��ptimaliza�ch � na kter�e m� d�t jako na posloupnost rozhodnut�pro kter�lat�zv. princip optimality.  Ten � �e optim��osloupnost rozhodnut��u vlastnost, �e pro libovoln�te� stav a rozhudnut�us��chna n�eduj� rozhodnut�ptim��zhledem k v��zhodnut�rvn�. D� �e pro ztr� tvaru \eqref{adi} plat�rincip optimality je snadn�e ho nal� nap�d v \cite{bertsekas1995dynamic}.
47
48P��en�lohy stochastick� �� aditivn�tr�u je tedy mo�n�ostupovat, jak je u ���moc�ynamick� programov� zvykem. Minim��odnotu st� ztr� od okam�iku $t$ do $N$ v z�slosti na $x_t$ ozna�e $J_t(x_t)$. M� pro ni ps�
49\begin{gather}
50J_N(x_N)=0\\
51J_t(x_t)=\min_{u_t \in U(x_t)}\E_{w_t}\left\{g_k(x_{t+1},u_t)+J_{t+1}(x_{t+1})\right\} \qquad t=0,\ldots,N-1.
52\end{gather}
53
54P��en�udeme postupovat od konce �� horizontu a postupn�ledat $J_t(x_t)$. Pro v� $x_{t+1}$ se pou�ije rovnice \eqref{sys}.
55Libovolnou ��trategii $\pi=\{\mu_0,\ldots,\mu_{N-1}\}$, kter�pl� syst�rovnic
56\begin{equation}
57\label{impl}
58J_t(x_t)=\E_{w_t}\left\{g_k(x_t,\mu_t(x_t),w_t)+J_{t+1}(f_t(x_t,\mu_t(x_t),w_t))\right\} \qquad t=0,\ldots,N-1
59\end{equation}
60pak nazveme optim��osloupnost�ozhodnut�
Note: See TracBrowser for help on using the browser.