Changeset 931

Show
Ignore:
Timestamp:
05/09/10 22:45:24 (14 years ago)
Author:
vahalam
Message:
 
Location:
applications/dual/IterativeLocal/text
Files:
2 modified

Legend:

Unmodified
Added
Removed
  • applications/dual/IterativeLocal/text/bpkomplet.lyx

    r922 r931  
    10081008\end_layout 
    10091009 
     1010\begin_layout Section 
     1011Dynamické programování 
     1012\end_layout 
     1013 
    10101014\begin_layout Subsection 
    10111015Formulace problému 
     
    12171221 
    12181222\begin_layout Subsection 
    1219 Dynamické programování 
     1223Přístup dynamického programování 
    12201224\end_layout 
    12211225 
     
    16011605 Nicméně dynamické programování je pouze obecný přístup pro iterativní optimaliz 
    16021606aci při uvažování nejistoty v systému. 
     1607\end_layout 
     1608 
     1609\begin_layout Section 
     1610Vliv neznalosti na systém 
    16031611\end_layout 
    16041612 
     
    23792387\end_layout 
    23802388 
     2389\begin_layout Section 
     2390Spojité systémy 
     2391\end_layout 
     2392 
    23812393\begin_layout Subsection 
    23822394Deterministické systémy se spojitým časem 
     
    28312843\end_layout 
    28322844 
    2833 \begin_layout Subsection 
     2845\begin_layout Section 
    28342846Algoritmy pro duální řízení 
    28352847\end_layout 
     
    29762988 použita pro srovnání s ústřední metodou této práce iLDP. 
    29772989 Dále pak do kategorie smíšených metod spadá právě i metoda iLDP, která 
    2978  bude podrobně popsána le. 
     2990 bude podrobně popsána v následující kapitole. 
    29792991  
    29802992\end_layout 
    29812993 
     2994\begin_layout Chapter 
     2995Algoritmy pro návrh řízení 
     2996\end_layout 
     2997 
    29822998\begin_layout Section 
    2983 Výběr konkrétních algoritmů pro srovnání 
     2999Výběr algoritmů pro srovnání 
    29843000\end_layout 
    29853001 
     
    31413157\end_layout 
    31423158 
    3143 \begin_layout Paragraph 
     3159\begin_layout Subsection 
    31443160Zobecněné iterativní LQG řízení 
    31453161\end_layout 
     
    39924008\emph default 
    39934009 parciálních derivací podle složek vektoru  
    3994 \begin_inset Formula $\mathbf{x}$ 
     4010\begin_inset Formula $x$ 
    39954011\end_inset 
    39964012 
     
    41104126 
    41114127\begin_layout Itemize 
    4112 vyšší přesnost 
    4113 \end_layout 
    4114  
    4115 \begin_layout Itemize 
    41164128rychlejší dosažení požadované hodnoty 
    41174129\end_layout 
     
    41704182\begin_layout Chapter 
    41714183Systémy pro testování 
     4184\begin_inset CommandInset label 
     4185LatexCommand label 
     4186name "cha:Systémy-pro-testování" 
     4187 
     4188\end_inset 
     4189 
     4190 
    41724191\end_layout 
    41734192 
     
    52395258\end_layout 
    52405259 
    5241 \begin_layout Subsection 
    5242 Metodika testování algoritmů 
    5243 \end_layout 
    5244  
    5245 \begin_layout Standard 
    5246 (metodika získávání výsledků) 
    5247 \end_layout 
    5248  
    52495260\begin_layout Section 
    52505261Synchronní motor s permanentními magnety 
     
    56705681\end_inset 
    56715682 
    5672  čtyři složky a koverianční matice  
     5683 čtyři složky a kovarianční matice  
    56735684\begin_inset Formula $P_{k}$ 
    56745685\end_inset 
     
    60386049 
    60396050\begin_layout Subsubsection 
    6040 LQG řizení 
     6051LQG řízení 
    60416052\end_layout 
    60426053 
     
    60446055Nyní můžeme na matice popisující systém 
    60456056\begin_inset Formula \begin{eqnarray*} 
    6046 A_{k}=\frac{d}{dx_{k}}g(x_{k},u_{k}) & = & \left(\begin{array}{ccccc} 
     6057A_{k} & = & \frac{d}{dx_{k}}g(x_{k},u_{k})\\ 
     6058 & = & \left(\begin{array}{ccccc} 
    60476059a & 0 & b\sin\vartheta_{k} & b\left(\psi_{k}+\overline{\omega}_{k}\right)\cos\vartheta_{k} & b\sin\vartheta_{k}\\ 
    604860600 & a & -b\cos\vartheta_{k} & b\left(\psi_{k}+\overline{\omega}_{k}\right)\sin\vartheta_{k} & -b\cos\vartheta_{k}\\ 
    60496061-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)\\ 
    605060620 & 0 & \Delta k & 1 & \Delta k\\ 
    6051 0 & 0 & 0 & 0 & 1\end{array}\right),\\ 
     60630 & 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*} 
    60526072B_{k}=\frac{d}{du_{k}}g(x_{k},u_{k}) & = & \left(\begin{array}{cc} 
    60536073c & 0\\ 
     
    61026122 
    61036123\begin_layout Section 
     6124Metodika zpracování a získávání výsledků 
     6125\end_layout 
     6126 
     6127\begin_layout Standard 
     6128Pro testované systémy viz kapitola 
     6129\begin_inset CommandInset ref 
     6130LatexCommand ref 
     6131reference "cha:Systémy-pro-testování" 
     6132 
     6133\end_inset 
     6134 
     6135, byly jednotlivé algoritmy implementovány jako funkce programu  
     6136\emph on 
     6137Matlab 
     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 
     6144Matlab 
     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 \[ 
     6157J_{alg}=\sum_{k=0}^{N-1}\left(y_{k+1}-r_{k+1}\right)^{2}.\] 
     6158 
     6159\end_inset 
     6160 
     6161Př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 
     6167Funkce pro jednoduchý systém 
     6168\end_layout 
     6169 
     6170\begin_layout Standard 
     6171Všechny konkrétní funkce programu  
     6172\emph on 
     6173Matlab 
     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 
     6260Ná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 
     6263LatexCommand ref 
     6264reference "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í 
     6270m 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 
    61046278Výsledky algoritmu iLDP 
    61056279\end_layout 
     
    61596333\begin_layout Subsection 
    61606334Porovnání algoritmů 
     6335\end_layout 
     6336 
     6337\begin_layout Subsection 
     6338Konfrontace s prvotními očekáváními 
    61616339\end_layout 
    61626340