Changeset 931
- Timestamp:
- 05/09/10 22:45:24 (15 years ago)
- Location:
- applications/dual/IterativeLocal/text
- Files:
-
- 2 modified
Legend:
- Unmodified
- Added
- Removed
-
applications/dual/IterativeLocal/text/bpkomplet.lyx
r922 r931 1008 1008 \end_layout 1009 1009 1010 \begin_layout Section 1011 Dynamické programování 1012 \end_layout 1013 1010 1014 \begin_layout Subsection 1011 1015 Formulace problému … … 1217 1221 1218 1222 \begin_layout Subsection 1219 Dynamicképrogramování1223 Přístup dynamického programování 1220 1224 \end_layout 1221 1225 … … 1601 1605 Nicméně dynamické programování je pouze obecný přístup pro iterativní optimaliz 1602 1606 aci při uvažování nejistoty v systému. 1607 \end_layout 1608 1609 \begin_layout Section 1610 Vliv neznalosti na systém 1603 1611 \end_layout 1604 1612 … … 2379 2387 \end_layout 2380 2388 2389 \begin_layout Section 2390 Spojité systémy 2391 \end_layout 2392 2381 2393 \begin_layout Subsection 2382 2394 Deterministické systémy se spojitým časem … … 2831 2843 \end_layout 2832 2844 2833 \begin_layout S ubsection2845 \begin_layout Section 2834 2846 Algoritmy pro duální řízení 2835 2847 \end_layout … … 2976 2988 použita pro srovnání s ústřední metodou této práce iLDP. 2977 2989 Dále pak do kategorie smíšených metod spadá právě i metoda iLDP, která 2978 bude podrobně popsána dále.2990 bude podrobně popsána v následující kapitole. 2979 2991 2980 2992 \end_layout 2981 2993 2994 \begin_layout Chapter 2995 Algoritmy pro návrh řízení 2996 \end_layout 2997 2982 2998 \begin_layout Section 2983 Výběr konkrétníchalgoritmů pro srovnání2999 Výběr algoritmů pro srovnání 2984 3000 \end_layout 2985 3001 … … 3141 3157 \end_layout 3142 3158 3143 \begin_layout Paragraph3159 \begin_layout Subsection 3144 3160 Zobecněné iterativní LQG řízení 3145 3161 \end_layout … … 3992 4008 \emph default 3993 4009 parciálních derivací podle složek vektoru 3994 \begin_inset Formula $ \mathbf{x}$4010 \begin_inset Formula $x$ 3995 4011 \end_inset 3996 4012 … … 4110 4126 4111 4127 \begin_layout Itemize 4112 vyšší přesnost4113 \end_layout4114 4115 \begin_layout Itemize4116 4128 rychlejší dosažení požadované hodnoty 4117 4129 \end_layout … … 4170 4182 \begin_layout Chapter 4171 4183 Systémy pro testování 4184 \begin_inset CommandInset label 4185 LatexCommand label 4186 name "cha:Systémy-pro-testování" 4187 4188 \end_inset 4189 4190 4172 4191 \end_layout 4173 4192 … … 5239 5258 \end_layout 5240 5259 5241 \begin_layout Subsection5242 Metodika testování algoritmů5243 \end_layout5244 5245 \begin_layout Standard5246 (metodika získávání výsledků)5247 \end_layout5248 5249 5260 \begin_layout Section 5250 5261 Synchronní motor s permanentními magnety … … 5670 5681 \end_inset 5671 5682 5672 čtyři složky a kov erianční matice5683 čtyři složky a kovarianční matice 5673 5684 \begin_inset Formula $P_{k}$ 5674 5685 \end_inset … … 6038 6049 6039 6050 \begin_layout Subsubsection 6040 LQG ř izení6051 LQG řízení 6041 6052 \end_layout 6042 6053 … … 6044 6055 Nyní můžeme na matice popisující systém 6045 6056 \begin_inset Formula \begin{eqnarray*} 6046 A_{k}=\frac{d}{dx_{k}}g(x_{k},u_{k}) & = & \left(\begin{array}{ccccc} 6057 A_{k} & = & \frac{d}{dx_{k}}g(x_{k},u_{k})\\ 6058 & = & \left(\begin{array}{ccccc} 6047 6059 a & 0 & b\sin\vartheta_{k} & b\left(\psi_{k}+\overline{\omega}_{k}\right)\cos\vartheta_{k} & b\sin\vartheta_{k}\\ 6048 6060 0 & a & -b\cos\vartheta_{k} & b\left(\psi_{k}+\overline{\omega}_{k}\right)\sin\vartheta_{k} & -b\cos\vartheta_{k}\\ 6049 6061 -e\sin\vartheta_{k} & e\cos\vartheta_{k} & d & -e\left(i_{\alpha,k}\cos\vartheta_{k}+i_{\beta,k}\sin\vartheta_{k}\right) & \left(d-1\right)\\ 6050 6062 0 & 0 & \Delta k & 1 & \Delta k\\ 6051 0 & 0 & 0 & 0 & 1\end{array}\right),\\ 6063 0 & 0 & 0 & 0 & 1\end{array}\right),\end{eqnarray*} 6064 6065 \end_inset 6066 6067 6068 \end_layout 6069 6070 \begin_layout Standard 6071 \begin_inset Formula \begin{eqnarray*} 6052 6072 B_{k}=\frac{d}{du_{k}}g(x_{k},u_{k}) & = & \left(\begin{array}{cc} 6053 6073 c & 0\\ … … 6102 6122 6103 6123 \begin_layout Section 6124 Metodika zpracování a získávání výsledků 6125 \end_layout 6126 6127 \begin_layout Standard 6128 Pro testované systémy viz kapitola 6129 \begin_inset CommandInset ref 6130 LatexCommand ref 6131 reference "cha:Systémy-pro-testování" 6132 6133 \end_inset 6134 6135 , byly jednotlivé algoritmy implementovány jako funkce programu 6136 \emph on 6137 Matlab 6138 \emph default 6139 . 6140 Vstupními hodnotami byla nastavení jednotlivých počátečních podmínek a 6141 parametrů pro výpočet. 6142 Jednotlivé algoritmy pak byly volány ze skriptu programu 6143 \emph on 6144 Matlab 6145 \emph default 6146 se stejným nastavením hodnot. 6147 Návratovými hodnotami jednotlivých funkcí reprezentujících algoritmy pak 6148 byla dosažená hodnota ztráty a posloupnost reprezentující diskrétní trajektorii 6149 systému ve stavovém prostoru. 6150 Hodnota ztráty 6151 \begin_inset Formula $J_{alg}$ 6152 \end_inset 6153 6154 byla pro účely porovnání algoritmů formulována stejně, jako součet kvadrátů 6155 odchylek výstupu od požadované hodnoty, tedy 6156 \begin_inset Formula \[ 6157 J_{alg}=\sum_{k=0}^{N-1}\left(y_{k+1}-r_{k+1}\right)^{2}.\] 6158 6159 \end_inset 6160 6161 Při porovnávání jednotlivých algoritmů nebyly samozřejmě uvažovány výsledky 6162 jednoho běhu výpočtu, ale provedlo se běhů více a následně byly výsledky 6163 zpracovány statisticky. 6164 \end_layout 6165 6166 \begin_layout Subsection 6167 Funkce pro jednoduchý systém 6168 \end_layout 6169 6170 \begin_layout Standard 6171 Všechny konkrétní funkce programu 6172 \emph on 6173 Matlab 6174 \emph default 6175 reprezentující jednotlivé algoritmy mají stejnou strukturu. 6176 Na počátku dostanou jako své vstupní hodnoty parametry: 6177 \end_layout 6178 6179 \begin_layout Labeling 6180 \labelwidthstring 00.00.0000 6181 \begin_inset Formula $y_{0}$ 6182 \end_inset 6183 6184 počáteční hodnota proměnné 6185 \begin_inset Formula $y$ 6186 \end_inset 6187 6188 , tedy hodnota 6189 \begin_inset Formula $y_{k}$ 6190 \end_inset 6191 6192 v čase 6193 \begin_inset Formula $k=0$ 6194 \end_inset 6195 6196 ; 6197 \end_layout 6198 6199 \begin_layout Labeling 6200 \labelwidthstring 00.00.0000 6201 \begin_inset Formula $y_{r}$ 6202 \end_inset 6203 6204 požadovaná hodnota 6205 \begin_inset Formula $y$ 6206 \end_inset 6207 6208 , referenční signál; 6209 \end_layout 6210 6211 \begin_layout Labeling 6212 \labelwidthstring 00.00.0000 6213 \begin_inset Formula $\mathrm{E}b$ 6214 \end_inset 6215 6216 střední hodnota neznámého parametru 6217 \begin_inset Formula $b$ 6218 \end_inset 6219 6220 ; 6221 \end_layout 6222 6223 \begin_layout Labeling 6224 \labelwidthstring 00.00.0000 6225 \begin_inset Formula $P$ 6226 \end_inset 6227 6228 variance neznámého parametru 6229 \begin_inset Formula $b$ 6230 \end_inset 6231 6232 ; 6233 \end_layout 6234 6235 \begin_layout Labeling 6236 \labelwidthstring 00.00.0000 6237 \begin_inset Formula $\sigma^{2}$ 6238 \end_inset 6239 6240 variance šumu; 6241 \end_layout 6242 6243 \begin_layout Labeling 6244 \labelwidthstring 00.00.0000 6245 \begin_inset Formula $K$ 6246 \end_inset 6247 6248 časový horizont, pro který navrhujeme řízení; 6249 \end_layout 6250 6251 \begin_layout Labeling 6252 \labelwidthstring 00.00.0000 6253 \begin_inset Formula $N$ 6254 \end_inset 6255 6256 počet vzorkových trajektorií pro simulaci. 6257 \end_layout 6258 6259 \begin_layout Standard 6260 Následně proběhne výpočet řízení na základě konkrétního algoritmu, zpravidla 6261 podle rovnic a vzorců popsaných v části 6262 \begin_inset CommandInset ref 6263 LatexCommand ref 6264 reference "sec:Jednoduchý-systém-pro-testovani" 6265 6266 \end_inset 6267 6268 . 6269 Dále je provedena simulace běhu systému s použítím řízení získaného v předchozí 6270 m kroku. 6271 V průběhu této simulace je vypočítána hodnota dosažené ztráty a posloupnost 6272 odpovídající trajektorie systému. 6273 Tyto veličiny charakterizující použití daného algoritmu na jednoduchý systém 6274 jsou na závěr vráceny jako návratové hodnoty funkce. 6275 \end_layout 6276 6277 \begin_layout Section 6104 6278 Výsledky algoritmu iLDP 6105 6279 \end_layout … … 6159 6333 \begin_layout Subsection 6160 6334 Porovnání algoritmů 6335 \end_layout 6336 6337 \begin_layout Subsection 6338 Konfrontace s prvotními očekáváními 6161 6339 \end_layout 6162 6340