Changeset 1448 for applications/dual/vahala/DP
- Timestamp:
- 04/10/12 15:30:42 (13 years ago)
- Location:
- applications/dual/vahala/DP
- Files:
-
- 3 modified
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified applications/dual/vahala/DP/prac_verz.lyx ¶
r1447 r1448 621 621 \begin_inset Formula 622 622 \[ 623 \alpha=k\left(a+b\cdot\cos(120^{\circ})+c\cdot\cos(-120^{\circ})\right)=k\left(a-\frac{1}{2}b-\frac{1}{2}c\right) 623 \alpha=k\left(a+b\cdot\cos(120^{\circ})+c\cdot\cos(-120^{\circ})\right)=k\left(a-\frac{1}{2}b-\frac{1}{2}c\right), 624 624 \] 625 625 … … 660 660 \begin_inset Formula 661 661 \[ 662 \beta=k\left(b\cdot\sin(120^{\circ})+c\cdot\sin(-120^{\circ})\right)=k\left(\frac{\sqrt{3}}{2}b-\frac{\sqrt{3}}{2}c\right) 662 \beta=k\left(b\cdot\sin(120^{\circ})+c\cdot\sin(-120^{\circ})\right)=k\left(\frac{\sqrt{3}}{2}b-\frac{\sqrt{3}}{2}c\right). 663 663 \] 664 664 … … 668 668 \begin_inset Formula 669 669 \begin{eqnarray*} 670 \alpha & = & \frac{2}{3}\left(a-\frac{1}{2}b-\frac{1}{2}c\right) \\671 \beta & = & \frac{\sqrt{3}}{3}\left(b-c\right) 670 \alpha & = & \frac{2}{3}\left(a-\frac{1}{2}b-\frac{1}{2}c\right),\\ 671 \beta & = & \frac{\sqrt{3}}{3}\left(b-c\right). 672 672 \end{eqnarray*} 673 673 … … 681 681 \begin_inset Formula 682 682 \begin{eqnarray*} 683 a & = & \alpha+\theta \\684 b & = & \left(-\frac{1}{2}\alpha+\frac{\sqrt{3}}{2}\beta\right)+\theta \\685 c & \text{=} & \left(-\frac{1}{2}\alpha-\frac{\sqrt{3}}{2}\beta\right)+\theta 683 a & = & \alpha+\theta,\\ 684 b & = & \left(-\frac{1}{2}\alpha+\frac{\sqrt{3}}{2}\beta\right)+\theta,\\ 685 c & \text{=} & \left(-\frac{1}{2}\alpha-\frac{\sqrt{3}}{2}\beta\right)+\theta, 686 686 \end{eqnarray*} 687 687 … … 746 746 \alpha\\ 747 747 \beta 748 \end{array}\right) \label{eq:transformace_al-be_na_d-q}748 \end{array}\right).\label{eq:transformace_al-be_na_d-q} 749 749 \end{eqnarray} 750 750 … … 767 767 d\\ 768 768 q 769 \end{array}\right) \label{eq:transformace_d-q_na_al-be}769 \end{array}\right).\label{eq:transformace_d-q_na_al-be} 770 770 \end{eqnarray} 771 771 … … 859 859 \begin_inset Formula 860 860 \begin{equation} 861 u_{s}=R_{s}i_{s}+\frac{d\psi_{s}}{dt} \label{eq:odvoz-statorove-napeti}861 u_{s}=R_{s}i_{s}+\frac{d\psi_{s}}{dt},\label{eq:odvoz-statorove-napeti} 862 862 \end{equation} 863 863 … … 913 913 \begin_inset Formula 914 914 \begin{equation} 915 \psi_{s}=L_{s}i_{s}+\psi_{pm}e^{j\vartheta} \label{eq:odvoz-magneticky-tok}915 \psi_{s}=L_{s}i_{s}+\psi_{pm}e^{j\vartheta},\label{eq:odvoz-magneticky-tok} 916 916 \end{equation} 917 917 … … 947 947 \begin_inset Formula 948 948 \[ 949 u_{s}=R_{s}i_{s}+\frac{d\left(L_{s}i_{s}+\psi_{pm}e^{j\vartheta}\right)}{dt}=R_{s}i_{s}+L_{s}\frac{di_{s}}{dt}+j\psi_{pm}\frac{d\vartheta}{dt}e^{j\vartheta} 949 u_{s}=R_{s}i_{s}+\frac{d\left(L_{s}i_{s}+\psi_{pm}e^{j\vartheta}\right)}{dt}=R_{s}i_{s}+L_{s}\frac{di_{s}}{dt}+j\psi_{pm}\frac{d\vartheta}{dt}e^{j\vartheta}. 950 950 \] 951 951 … … 959 959 \begin_inset Formula 960 960 \begin{equation} 961 \omega=\frac{d\vartheta}{dt} \label{eq:definice-otacek}961 \omega=\frac{d\vartheta}{dt}.\label{eq:definice-otacek} 962 962 \end{equation} 963 963 … … 980 980 \begin_inset Formula 981 981 \begin{eqnarray*} 982 u_{\alpha} & = & R_{s}i_{\alpha}+L_{s}\frac{di_{\alpha}}{dt}-\psi_{pm}\omega\sin\vartheta \\983 u_{\beta} & = & R_{s}i_{\beta}+L_{s}\frac{di_{\beta}}{dt}+\psi_{pm}\omega\cos\vartheta 982 u_{\alpha} & = & R_{s}i_{\alpha}+L_{s}\frac{di_{\alpha}}{dt}-\psi_{pm}\omega\sin\vartheta,\\ 983 u_{\beta} & = & R_{s}i_{\beta}+L_{s}\frac{di_{\beta}}{dt}+\psi_{pm}\omega\cos\vartheta, 984 984 \end{eqnarray*} 985 985 … … 989 989 \begin_inset Formula 990 990 \begin{eqnarray} 991 \frac{di_{\alpha}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\alpha}+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta+\frac{1}{L_{s}}u_{\alpha} \nonumber \\992 \frac{di_{\beta}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\beta}-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta+\frac{1}{L_{s}}u_{\beta} \label{eq:rovnice-proudy-ls}991 \frac{di_{\alpha}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\alpha}+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta+\frac{1}{L_{s}}u_{\alpha},\nonumber \\ 992 \frac{di_{\beta}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\beta}-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta+\frac{1}{L_{s}}u_{\beta}.\label{eq:rovnice-proudy-ls} 993 993 \end{eqnarray} 994 994 … … 1107 1107 \begin_inset Formula 1108 1108 \begin{equation} 1109 \omega=p_{p}\omega_{mech} \label{eq:vztah-el-a-mech-omega}1109 \omega=p_{p}\omega_{mech},\label{eq:vztah-el-a-mech-omega} 1110 1110 \end{equation} 1111 1111 … … 1125 1125 \begin_inset Formula 1126 1126 \begin{equation} 1127 T_{1}+\ldots+T_{n}=\frac{dl}{dt}=\frac{d\left(J\omega_{mech}\right)}{dt}=J\frac{d\omega_{mech}}{dt} \label{eq:rovnice-momenty-preddosaz}1127 T_{1}+\ldots+T_{n}=\frac{dl}{dt}=\frac{d\left(J\omega_{mech}\right)}{dt}=J\frac{d\omega_{mech}}{dt}.\label{eq:rovnice-momenty-preddosaz} 1128 1128 \end{equation} 1129 1129 … … 1184 1184 \begin_inset Formula 1185 1185 \begin{equation} 1186 T_{el}-T_{L}-B\omega_{mech}=J\frac{d\omega_{mech}}{dt} \label{eq:rovnice-momenty-dosazeno}1186 T_{el}-T_{L}-B\omega_{mech}=J\frac{d\omega_{mech}}{dt}.\label{eq:rovnice-momenty-dosazeno} 1187 1187 \end{equation} 1188 1188 … … 1227 1227 \begin_inset Formula 1228 1228 \begin{equation} 1229 P=k_{p}\left(u_{\alpha}i_{\alpha}+u_{\beta}i_{\beta}\right) \label{eq:rovnice-vykon}1229 P=k_{p}\left(u_{\alpha}i_{\alpha}+u_{\beta}i_{\beta}\right),\label{eq:rovnice-vykon} 1230 1230 \end{equation} 1231 1231 … … 1271 1271 \begin_inset Formula 1272 1272 \[ 1273 u_{ind}=\frac{d\psi_{s}}{dt}=\frac{d\left(L_{s}i_{s}+\psi_{pm}e^{j\vartheta}\right)}{dt}=L_{s}\frac{di_{s}}{dt}+j\psi_{pm}\omega e^{j\vartheta} 1273 u_{ind}=\frac{d\psi_{s}}{dt}=\frac{d\left(L_{s}i_{s}+\psi_{pm}e^{j\vartheta}\right)}{dt}=L_{s}\frac{di_{s}}{dt}+j\psi_{pm}\omega e^{j\vartheta}. 1274 1274 \] 1275 1275 … … 1286 1286 \begin_inset Formula 1287 1287 \begin{eqnarray*} 1288 u_{ind,\alpha} & = & -\psi_{pm}\omega\sin\vartheta \\1288 u_{ind,\alpha} & = & -\psi_{pm}\omega\sin\vartheta,\\ 1289 1289 u_{ind,\beta} & = & \psi_{pm}\omega\cos\vartheta 1290 1290 \end{eqnarray*} … … 1302 1302 \begin_inset Formula 1303 1303 \begin{equation} 1304 P=k_{p}\left(-\psi_{pm}i_{\alpha}\omega\sin\vartheta+\psi_{pm}i_{\beta}\omega\cos\vartheta\right) \label{eq:rovnice-vykon-dosazano}1304 P=k_{p}\left(-\psi_{pm}i_{\alpha}\omega\sin\vartheta+\psi_{pm}i_{\beta}\omega\cos\vartheta\right).\label{eq:rovnice-vykon-dosazano} 1305 1305 \end{equation} 1306 1306 … … 1333 1333 \begin_inset Formula 1334 1334 \[ 1335 T_{el}=\frac{P}{\omega_{mech}}=\frac{k_{p}}{\omega_{mech}}\left(-\psi_{pm}i_{\alpha}\omega\sin\vartheta+\psi_{pm}i_{\beta}\omega\cos\vartheta\right) 1335 T_{el}=\frac{P}{\omega_{mech}}=\frac{k_{p}}{\omega_{mech}}\left(-\psi_{pm}i_{\alpha}\omega\sin\vartheta+\psi_{pm}i_{\beta}\omega\cos\vartheta\right), 1336 1336 \] 1337 1337 … … 1348 1348 \begin_inset Formula 1349 1349 \[ 1350 T_{el}=k_{p}p_{p}\left(-\psi_{pm}i_{\alpha}\sin\vartheta+\psi_{pm}i_{\beta}\cos\vartheta\right) 1350 T_{el}=k_{p}p_{p}\left(-\psi_{pm}i_{\alpha}\sin\vartheta+\psi_{pm}i_{\beta}\cos\vartheta\right). 1351 1351 \] 1352 1352 … … 1375 1375 \begin_inset Formula 1376 1376 \[ 1377 k_{p}p_{p}\psi_{pm}\left(-i_{\alpha}\sin\vartheta+i_{\beta}\cos\vartheta\right)-T_{L}-B\omega_{mech}=J\frac{d\omega_{mech}}{dt} 1377 k_{p}p_{p}\psi_{pm}\left(-i_{\alpha}\sin\vartheta+i_{\beta}\cos\vartheta\right)-T_{L}-B\omega_{mech}=J\frac{d\omega_{mech}}{dt}. 1378 1378 \] 1379 1379 … … 1394 1394 \begin_inset Formula 1395 1395 \begin{equation} 1396 \frac{d\omega}{dt}=\frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\left(i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta\right)-\frac{p_{p}}{J}T_{L}-\frac{B}{J}\omega \label{eq:rovnice-pro-omega-ls}1396 \frac{d\omega}{dt}=\frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\left(i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta\right)-\frac{p_{p}}{J}T_{L}-\frac{B}{J}\omega.\label{eq:rovnice-pro-omega-ls} 1397 1397 \end{equation} 1398 1398 … … 1502 1502 \begin_inset Formula 1503 1503 \[ 1504 u_{r}e^{j\vartheta}=R_{s}i_{r}e^{j\vartheta}+\frac{d\psi_{r}}{dt}e^{j\vartheta}+j\psi_{r}\omega e^{j\vartheta} 1504 u_{r}e^{j\vartheta}=R_{s}i_{r}e^{j\vartheta}+\frac{d\psi_{r}}{dt}e^{j\vartheta}+j\psi_{r}\omega e^{j\vartheta}. 1505 1505 \] 1506 1506 … … 1514 1514 \begin_inset Formula 1515 1515 \begin{equation} 1516 u_{r}=R_{s}i_{r}+\frac{d\psi_{r}}{dt}+j\psi_{r}\omega \label{eq:odvoz-ldq-rovnice-napeti}1516 u_{r}=R_{s}i_{r}+\frac{d\psi_{r}}{dt}+j\psi_{r}\omega.\label{eq:odvoz-ldq-rovnice-napeti} 1517 1517 \end{equation} 1518 1518 … … 1527 1527 \begin_inset Formula 1528 1528 \begin{eqnarray*} 1529 \psi_{d} & = & L_{d}i_{d}+\psi_{pm} \\1529 \psi_{d} & = & L_{d}i_{d}+\psi_{pm},\\ 1530 1530 \psi_{q} & = & L_{q}i_{q} 1531 1531 \end{eqnarray*} … … 1558 1558 \begin_inset Formula 1559 1559 \begin{eqnarray*} 1560 u_{d} & = & R_{s}i_{d}+L_{d}\frac{di_{d}}{dt}-L_{q}i_{q}\omega \\1561 u_{q} & = & R_{s}i_{q}+L_{q}\frac{di_{q}}{dt}+\left(L_{d}i_{d}+\psi_{pm}\right)\omega 1560 u_{d} & = & R_{s}i_{d}+L_{d}\frac{di_{d}}{dt}-L_{q}i_{q}\omega,\\ 1561 u_{q} & = & R_{s}i_{q}+L_{q}\frac{di_{q}}{dt}+\left(L_{d}i_{d}+\psi_{pm}\right)\omega. 1562 1562 \end{eqnarray*} 1563 1563 … … 1572 1572 \begin_inset Formula 1573 1573 \begin{eqnarray} 1574 \frac{di_{d}}{dt} & = & -\frac{R_{s}}{L_{d}}i_{d}+\frac{L_{q}}{L_{d}}i_{q}\omega+\frac{1}{L_{d}}u_{d} \nonumber \\1575 \frac{di_{q}}{dt} & = & -\frac{R_{s}}{L_{q}}i_{q}-\frac{L_{d}}{L_{q}}i_{d}\omega-\frac{\psi_{pm}}{L_{q}}\omega+\frac{1}{L_{q}}u_{q} \label{eq:rovnice-proudy-ldq-v-dq}1574 \frac{di_{d}}{dt} & = & -\frac{R_{s}}{L_{d}}i_{d}+\frac{L_{q}}{L_{d}}i_{q}\omega+\frac{1}{L_{d}}u_{d},\nonumber \\ 1575 \frac{di_{q}}{dt} & = & -\frac{R_{s}}{L_{q}}i_{q}-\frac{L_{d}}{L_{q}}i_{d}\omega-\frac{\psi_{pm}}{L_{q}}\omega+\frac{1}{L_{q}}u_{q}.\label{eq:rovnice-proudy-ldq-v-dq} 1576 1576 \end{eqnarray} 1577 1577 … … 1624 1624 \begin_inset Formula 1625 1625 \[ 1626 T_{el}-T_{L}-B\omega_{mech}=J\frac{d\omega_{mech}}{dt} 1626 T_{el}-T_{L}-B\omega_{mech}=J\frac{d\omega_{mech}}{dt}, 1627 1627 \] 1628 1628 … … 1651 1651 \begin_inset Formula 1652 1652 \[ 1653 P=k_{p}\left(u_{\alpha}i_{\alpha}+u_{\beta}i_{\beta}\right)=k_{p}\left(u_{d}i_{d}+u_{q}i_{q}\right) 1653 P=k_{p}\left(u_{\alpha}i_{\alpha}+u_{\beta}i_{\beta}\right)=k_{p}\left(u_{d}i_{d}+u_{q}i_{q}\right). 1654 1654 \] 1655 1655 … … 1660 1660 \begin_inset Formula 1661 1661 \begin{eqnarray*} 1662 u_{ind,d} & = & -L_{q}i_{q}\omega \\1662 u_{ind,d} & = & -L_{q}i_{q}\omega,\\ 1663 1663 u_{ind,q} & = & \left(L_{d}i_{d}+\psi_{pm}\right)\omega 1664 1664 \end{eqnarray*} … … 1669 1669 \begin_inset Formula 1670 1670 \[ 1671 P=k_{p}\left(\left(L_{d}-L_{q}\right)i_{d}i_{q}+\psi_{pm}i_{q}\right)\omega 1671 P=k_{p}\left(\left(L_{d}-L_{q}\right)i_{d}i_{q}+\psi_{pm}i_{q}\right)\omega. 1672 1672 \] 1673 1673 … … 1721 1721 \begin_inset Formula 1722 1722 \begin{equation} 1723 \frac{d\omega}{dt}=\frac{k_{p}p_{p}^{2}}{J}\left(\left(L_{d}-L_{q}\right)i_{d}i_{q}+\psi_{pm}i_{q}\right)-\frac{p_{p}}{J}T_{L}-\frac{B}{J}\omega \label{eq:rovnice-pro-omega-ruzne-ldq}1723 \frac{d\omega}{dt}=\frac{k_{p}p_{p}^{2}}{J}\left(\left(L_{d}-L_{q}\right)i_{d}i_{q}+\psi_{pm}i_{q}\right)-\frac{p_{p}}{J}T_{L}-\frac{B}{J}\omega,\label{eq:rovnice-pro-omega-ruzne-ldq} 1724 1724 \end{equation} 1725 1725 … … 1771 1771 \begin_inset Formula 1772 1772 \begin{eqnarray} 1773 \frac{di_{\alpha}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\alpha}+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta+\frac{1}{L_{s}}u_{\alpha} \nonumber \\1774 \frac{di_{\beta}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\beta}-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta+\frac{1}{L_{s}}u_{\beta} \label{eq:rovnice-pmsm-albe-ls}\\1775 \frac{d\omega}{dt} & = & \frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\left(i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta\right)-\frac{p_{p}}{J}T_{L}-\frac{B}{J}\omega \nonumber \\1776 \frac{d\vartheta}{dt} & = & \omega \nonumber1773 \frac{di_{\alpha}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\alpha}+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta+\frac{1}{L_{s}}u_{\alpha},\nonumber \\ 1774 \frac{di_{\beta}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\beta}-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta+\frac{1}{L_{s}}u_{\beta},\label{eq:rovnice-pmsm-albe-ls}\\ 1775 \frac{d\omega}{dt} & = & \frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\left(i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta\right)-\frac{p_{p}}{J}T_{L}-\frac{B}{J}\omega,\nonumber \\ 1776 \frac{d\vartheta}{dt} & = & \omega.\nonumber 1777 1777 \end{eqnarray} 1778 1778 … … 1815 1815 \begin_inset Formula 1816 1816 \begin{eqnarray} 1817 \frac{di_{d}}{dt} & = & -\frac{R_{s}}{L_{d}}i_{d}+\frac{L_{q}}{L_{d}}i_{q}\omega+\frac{1}{L_{d}}u_{d} \nonumber \\1818 \frac{di_{q}}{dt} & = & -\frac{R_{s}}{L_{q}}i_{q}-\frac{L_{d}}{L_{q}}i_{d}\omega-\frac{\psi_{pm}}{L_{q}}\omega+\frac{1}{L_{q}}u_{q} \label{eq:rovnice-pmsm-dq-ldq}\\1819 \frac{d\omega}{dt} & = & \frac{k_{p}p_{p}^{2}}{J}\left(\left(L_{d}-L_{q}\right)i_{d}i_{q}+\psi_{pm}i_{q}\right)-\frac{p_{p}}{J}T_{L}-\frac{B}{J}\omega \nonumber \\1820 \frac{d\vartheta}{dt} & = & \omega \nonumber1817 \frac{di_{d}}{dt} & = & -\frac{R_{s}}{L_{d}}i_{d}+\frac{L_{q}}{L_{d}}i_{q}\omega+\frac{1}{L_{d}}u_{d},\nonumber \\ 1818 \frac{di_{q}}{dt} & = & -\frac{R_{s}}{L_{q}}i_{q}-\frac{L_{d}}{L_{q}}i_{d}\omega-\frac{\psi_{pm}}{L_{q}}\omega+\frac{1}{L_{q}}u_{q},\label{eq:rovnice-pmsm-dq-ldq}\\ 1819 \frac{d\omega}{dt} & = & \frac{k_{p}p_{p}^{2}}{J}\left(\left(L_{d}-L_{q}\right)i_{d}i_{q}+\psi_{pm}i_{q}\right)-\frac{p_{p}}{J}T_{L}-\frac{B}{J}\omega,\nonumber \\ 1820 \frac{d\vartheta}{dt} & = & \omega.\nonumber 1821 1821 \end{eqnarray} 1822 1822 … … 1954 1954 \end_inset 1955 1955 1956 ) 1957 \begin_inset Formula 1958 \begin{eqnarray*} 1959 \frac{di_{\alpha}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\alpha}+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta+\frac{1}{L_{s}}u_{\alpha}\\ 1960 \frac{di_{\beta}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\beta}-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta+\frac{1}{L_{s}}u_{\beta}\\ 1961 \frac{d\omega}{dt} & = & \frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\left(i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta\right)-\frac{p_{p}}{J}T_{L}-\frac{B}{J}\omega\\ 1962 \frac{d\vartheta}{dt} & = & \omega 1963 \end{eqnarray*} 1964 1965 \end_inset 1966 1967 a užijeme zmiňované Eulerovy metody. 1956 )a užijeme zmiňované Eulerovy metody. 1968 1957 Derivaci tedy nahradíme konečnou diferencí 1969 1958 \begin_inset Formula 1970 1959 \[ 1971 \frac{dx}{dt}\left(t\right)=\frac{x_{t+1}-x_{t}}{\Delta t} 1960 \frac{dx}{dt}\left(t\right)=\frac{x_{t+1}-x_{t}}{\Delta t}, 1972 1961 \] 1973 1962 … … 1982 1971 \begin_inset Formula 1983 1972 \begin{eqnarray*} 1984 i_{\alpha,t+1} & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{\alpha,t}+\frac{\psi_{pm}\Delta t}{L_{s}}\omega_{t}\sin\vartheta_{t}+\frac{\Delta t}{L_{s}}u_{\alpha,t} \\1985 i_{\beta,t+1} & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{\beta,t}-\frac{\psi_{pm}\Delta t}{L_{s}}\omega_{t}\cos\vartheta_{t}+\frac{\Delta t}{L_{s}}u_{\beta,t} \\1986 \omega_{t+1} & \text{=} & \left(1-\frac{B}{J}\Delta t\right)\omega_{t}+\frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\Delta t\left(i_{\beta,t}\cos\vartheta_{t}-i_{\alpha,t}\sin\vartheta_{t}\right)-\frac{p_{p}}{J}T_{L}\Delta t \\1987 \vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t 1973 i_{\alpha,t+1} & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{\alpha,t}+\frac{\psi_{pm}\Delta t}{L_{s}}\omega_{t}\sin\vartheta_{t}+\frac{\Delta t}{L_{s}}u_{\alpha,t},\\ 1974 i_{\beta,t+1} & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{\beta,t}-\frac{\psi_{pm}\Delta t}{L_{s}}\omega_{t}\cos\vartheta_{t}+\frac{\Delta t}{L_{s}}u_{\beta,t},\\ 1975 \omega_{t+1} & \text{=} & \left(1-\frac{B}{J}\Delta t\right)\omega_{t}+\frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\Delta t\left(i_{\beta,t}\cos\vartheta_{t}-i_{\alpha,t}\sin\vartheta_{t}\right)-\frac{p_{p}}{J}T_{L}\Delta t,\\ 1976 \vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t. 1988 1977 \end{eqnarray*} 1989 1978 … … 1993 1982 \begin_inset Formula 1994 1983 \begin{eqnarray} 1995 a & = & 1-\frac{R_{s}}{L_{s}}\Delta t \nonumber \\1996 b & = & \frac{\psi_{pm}}{L_{s}}\Delta t \nonumber \\1997 c & = & \frac{\Delta t}{L_{s}} \label{eq:zjednodus-znaceni-konstant}\\1998 d & = & 1-\frac{B}{J}\Delta t \nonumber \\1984 a & = & 1-\frac{R_{s}}{L_{s}}\Delta t,\nonumber \\ 1985 b & = & \frac{\psi_{pm}}{L_{s}}\Delta t,\nonumber \\ 1986 c & = & \frac{\Delta t}{L_{s}},\label{eq:zjednodus-znaceni-konstant}\\ 1987 d & = & 1-\frac{B}{J}\Delta t,\nonumber \\ 1999 1988 e & = & \frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\Delta t\nonumber 2000 1989 \end{eqnarray} … … 2009 1998 \begin_inset Formula 2010 1999 \begin{eqnarray} 2011 i_{\alpha,t+1} & \text{=} & ai_{\alpha,t}+b\omega_{t}\sin\vartheta_{t}+cu_{\alpha,t} \nonumber \\2012 i_{\beta,t+1} & \text{=} & ai_{\beta,t}-b\omega_{t}\cos\vartheta_{t}+cu_{\beta,t} \label{eq:diskretni-system-albe-ls}\\2013 \omega_{t+1} & \text{=} & d\omega_{t}+e\left(i_{\beta,t}\cos\vartheta_{t}-i_{\alpha,t}\sin\vartheta_{t}\right) \nonumber \\2014 \vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t \nonumber2000 i_{\alpha,t+1} & \text{=} & ai_{\alpha,t}+b\omega_{t}\sin\vartheta_{t}+cu_{\alpha,t},\nonumber \\ 2001 i_{\beta,t+1} & \text{=} & ai_{\beta,t}-b\omega_{t}\cos\vartheta_{t}+cu_{\beta,t},\label{eq:diskretni-system-albe-ls}\\ 2002 \omega_{t+1} & \text{=} & d\omega_{t}+e\left(i_{\beta,t}\cos\vartheta_{t}-i_{\alpha,t}\sin\vartheta_{t}\right),\nonumber \\ 2003 \vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t.\nonumber 2015 2004 \end{eqnarray} 2016 2005 … … 2036 2025 \end_inset 2037 2026 2038 ) 2039 \begin_inset Formula 2040 \begin{eqnarray*} 2041 \frac{di_{\alpha}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\alpha}+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta+\frac{1}{L_{s}}u_{\alpha}\\ 2042 \frac{di_{\beta}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\beta}-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta+\frac{1}{L_{s}}u_{\beta}\\ 2043 \frac{d\omega}{dt} & = & \frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\left(i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta\right)-\frac{p_{p}}{J}T_{L}-\frac{B}{J}\omega\\ 2044 \frac{d\vartheta}{dt} & = & \omega 2045 \end{eqnarray*} 2046 2047 \end_inset 2048 2049 a pomocí přavodního vztahu ( 2027 )a pomocí přavodního vztahu ( 2050 2028 \begin_inset CommandInset ref 2051 2029 LatexCommand ref … … 2065 2043 \begin_inset Formula 2066 2044 \begin{eqnarray} 2067 \frac{d}{dt}\left(i_{d}\cos\vartheta-i_{q}\sin\vartheta\right) & = & -\frac{R_{s}}{L_{s}}\left(i_{d}\cos\vartheta-i_{q}\sin\vartheta\right)+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta+\frac{1}{L_{s}}\left(u_{d}\cos\vartheta-u_{q}\sin\vartheta\right) \nonumber \\2068 \frac{d}{dt}\left(i_{d}\sin\vartheta+i_{q}\cos\vartheta\right) & = & -\frac{R_{s}}{L_{s}}\left(i_{d}\sin\vartheta+i_{q}\cos\vartheta\right)-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta+\frac{1}{L_{s}}\left(u_{d}\sin\vartheta+u_{q}\cos\vartheta\right) \label{eq:uprava-rovnic-z-albe-na-dq-pro-ls}2045 \frac{d}{dt}\left(i_{d}\cos\vartheta-i_{q}\sin\vartheta\right) & = & -\frac{R_{s}}{L_{s}}\left(i_{d}\cos\vartheta-i_{q}\sin\vartheta\right)+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta+\frac{1}{L_{s}}\left(u_{d}\cos\vartheta-u_{q}\sin\vartheta\right),\nonumber \\ 2046 \frac{d}{dt}\left(i_{d}\sin\vartheta+i_{q}\cos\vartheta\right) & = & -\frac{R_{s}}{L_{s}}\left(i_{d}\sin\vartheta+i_{q}\cos\vartheta\right)-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta+\frac{1}{L_{s}}\left(u_{d}\sin\vartheta+u_{q}\cos\vartheta\right).\label{eq:uprava-rovnic-z-albe-na-dq-pro-ls} 2069 2047 \end{eqnarray} 2070 2048 … … 2081 2059 \begin_inset Formula 2082 2060 \begin{eqnarray*} 2083 \frac{d}{dt}\left(i_{d}\cos\vartheta-i_{q}\sin\vartheta\right) & = & \frac{di_{d}}{dt}\cos\vartheta-\frac{di_{q}}{dt}\sin\vartheta-i_{d}\omega\sin\vartheta-i_{q}\omega\cos\vartheta \\2061 \frac{d}{dt}\left(i_{d}\cos\vartheta-i_{q}\sin\vartheta\right) & = & \frac{di_{d}}{dt}\cos\vartheta-\frac{di_{q}}{dt}\sin\vartheta-i_{d}\omega\sin\vartheta-i_{q}\omega\cos\vartheta,\\ 2084 2062 \frac{d}{dt}\left(i_{d}\sin\vartheta+i_{q}\cos\vartheta\right) & = & \frac{di_{d}}{dt}\sin\vartheta+\frac{di_{q}}{dt}\cos\vartheta+i_{d}\omega\cos\vartheta-i_{q}\omega\sin\vartheta 2085 2063 \end{eqnarray*} … … 2152 2130 \begin_inset Formula 2153 2131 \begin{eqnarray*} 2154 \frac{di_{d}}{dt} & \text{=} & -\frac{R_{s}}{L_{s}}i_{d}-i_{q}\omega+\frac{u_{d}}{L_{s}} \\2155 \frac{di_{q}}{dt} & \text{=} & i_{d}\omega-\frac{R_{s}}{L_{s}}i_{q}-\frac{\psi_{pm}}{L_{s}}\omega+\frac{u_{q}}{L_{s}} 2132 \frac{di_{d}}{dt} & \text{=} & -\frac{R_{s}}{L_{s}}i_{d}-i_{q}\omega+\frac{u_{d}}{L_{s}},\\ 2133 \frac{di_{q}}{dt} & \text{=} & i_{d}\omega-\frac{R_{s}}{L_{s}}i_{q}-\frac{\psi_{pm}}{L_{s}}\omega+\frac{u_{q}}{L_{s}}. 2156 2134 \end{eqnarray*} 2157 2135 … … 2173 2151 \begin_inset Formula 2174 2152 \[ 2175 \frac{d\omega}{dt}\text{=}\frac{k_{p}p_{p}^{2}\psi_{pm}}{J}i_{q}-\frac{B}{J}\omega-\frac{p_{p}}{J}T_{L} 2153 \frac{d\omega}{dt}\text{=}\frac{k_{p}p_{p}^{2}\psi_{pm}}{J}i_{q}-\frac{B}{J}\omega-\frac{p_{p}}{J}T_{L}. 2176 2154 \] 2177 2155 … … 2181 2159 \begin_inset Formula 2182 2160 \[ 2183 \frac{d\vartheta}{dt}\text{=}\omega 2161 \frac{d\vartheta}{dt}\text{=}\omega. 2184 2162 \] 2185 2163 … … 2198 2176 \begin_inset Formula 2199 2177 \begin{eqnarray*} 2200 i_{d,t+1} & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{d,t}+\Delta t\cdot i_{q,t}\omega_{t}+\frac{\Delta t}{L_{s}}u_{d,t} \\2201 i_{q,t+1} & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{q,t}-\Delta t\cdot i_{d,t}\omega_{t}-\frac{\psi_{pm}\Delta t}{L_{s}}\omega_{t}+\frac{\Delta t}{L_{s}}u_{q,t} \\2202 \omega_{t+1} & \text{=} & \left(1-\frac{B}{J}\Delta t\right)\omega_{t}+\Delta t\frac{k_{p}p_{p}^{2}\psi_{pm}}{J}i_{q,t}-\frac{p_{p}}{J}T_{L}\Delta t \\2203 \vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t 2178 i_{d,t+1} & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{d,t}+\Delta t\cdot i_{q,t}\omega_{t}+\frac{\Delta t}{L_{s}}u_{d,t},\\ 2179 i_{q,t+1} & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{q,t}-\Delta t\cdot i_{d,t}\omega_{t}-\frac{\psi_{pm}\Delta t}{L_{s}}\omega_{t}+\frac{\Delta t}{L_{s}}u_{q,t},\\ 2180 \omega_{t+1} & \text{=} & \left(1-\frac{B}{J}\Delta t\right)\omega_{t}+\Delta t\frac{k_{p}p_{p}^{2}\psi_{pm}}{J}i_{q,t}-\frac{p_{p}}{J}T_{L}\Delta t,\\ 2181 \vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t. 2204 2182 \end{eqnarray*} 2205 2183 … … 2220 2198 \begin_inset Formula 2221 2199 \begin{eqnarray} 2222 i_{d,t+1} & \text{=} & ai_{d,t}+\Delta t\cdot i_{q,t}\omega_{t}+cu_{d,t} \nonumber \\2223 i_{q,t+1} & \text{=} & ai_{q,t}-\Delta t\cdot i_{d,t}\omega_{t}-b\omega_{t}+cu_{q,t} \label{eq:diskretni-system-dq-ls}\\2224 \omega_{t+1} & \text{=} & d\omega_{t}+ei_{q,t} \nonumber \\2225 \vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t \nonumber2200 i_{d,t+1} & \text{=} & ai_{d,t}+\Delta t\cdot i_{q,t}\omega_{t}+cu_{d,t},\nonumber \\ 2201 i_{q,t+1} & \text{=} & ai_{q,t}-\Delta t\cdot i_{d,t}\omega_{t}-b\omega_{t}+cu_{q,t},\label{eq:diskretni-system-dq-ls}\\ 2202 \omega_{t+1} & \text{=} & d\omega_{t}+ei_{q,t},\nonumber \\ 2203 \vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t.\nonumber 2226 2204 \end{eqnarray} 2227 2205 … … 2247 2225 \end_inset 2248 2226 2249 ) 2250 \begin_inset Formula 2251 \begin{eqnarray*} 2252 \frac{di_{d}}{dt} & = & -\frac{R_{s}}{L_{d}}i_{d}+\frac{L_{q}}{L_{d}}i_{q}\omega+\frac{1}{L_{d}}u_{d}\\ 2253 \frac{di_{q}}{dt} & = & -\frac{R_{s}}{L_{q}}i_{q}-\frac{L_{d}}{L_{q}}i_{d}\omega-\frac{\psi_{pm}}{L_{q}}\omega+\frac{1}{L_{q}}u_{q}\\ 2254 \frac{d\omega}{dt} & = & \frac{k_{p}p_{p}^{2}}{J}\left(\left(L_{d}-L_{q}\right)i_{d}i_{q}+\psi_{pm}i_{q}\right)-\frac{p_{p}}{J}T_{L}-\frac{B}{J}\omega\\ 2255 \frac{d\vartheta}{dt} & = & \omega 2256 \end{eqnarray*} 2257 2258 \end_inset 2259 2260 a diskretizaci provedeme opět stejným způsobem pomocí Eulerovy metody. 2227 )a diskretizaci provedeme opět stejným způsobem pomocí Eulerovy metody. 2261 2228 Popis PMSM pomocí diferenčních rovnic v souřadné soustavě 2262 2229 \begin_inset Formula $d-q$ … … 2274 2241 \begin_inset Formula 2275 2242 \begin{eqnarray} 2276 i_{d,t+1} & = & \left(1-\frac{R_{s}}{L_{d}}\Delta t\right)i_{d,t}+\frac{L_{q}\Delta t}{L_{d}}i_{q,t}\omega_{t}+\frac{\Delta t}{L_{d}}u_{d,t} \nonumber \\2277 i_{q,t+1} & = & \left(1-\frac{R_{s}}{L_{q}}\Delta t\right)i_{q,t}-\frac{L_{d}\Delta t}{L_{q}}i_{d,t}\omega_{t}-\frac{\psi_{pm}\Delta t}{L_{q}}\omega_{t}+\frac{\Delta t}{L_{q}}u_{q,t} \label{eq:diskretni-system-dq-ldq}\\2278 \omega_{t+1} & = & \left(1-\frac{B}{J}\Delta t\right)\omega_{t}+\frac{k_{p}p_{p}^{2}\Delta t}{J}\left(\left(L_{d}-L_{q}\right)i_{d,t}i_{q,t}+\psi_{pm}i_{q,t}\right) \nonumber \\2279 \vartheta_{t+1} & = & \vartheta_{t}+\omega_{t}\Delta t \nonumber2243 i_{d,t+1} & = & \left(1-\frac{R_{s}}{L_{d}}\Delta t\right)i_{d,t}+\frac{L_{q}\Delta t}{L_{d}}i_{q,t}\omega_{t}+\frac{\Delta t}{L_{d}}u_{d,t},\nonumber \\ 2244 i_{q,t+1} & = & \left(1-\frac{R_{s}}{L_{q}}\Delta t\right)i_{q,t}-\frac{L_{d}\Delta t}{L_{q}}i_{d,t}\omega_{t}-\frac{\psi_{pm}\Delta t}{L_{q}}\omega_{t}+\frac{\Delta t}{L_{q}}u_{q,t},\label{eq:diskretni-system-dq-ldq}\\ 2245 \omega_{t+1} & = & \left(1-\frac{B}{J}\Delta t\right)\omega_{t}+\frac{k_{p}p_{p}^{2}\Delta t}{J}\left(\left(L_{d}-L_{q}\right)i_{d,t}i_{q,t}+\psi_{pm}i_{q,t}\right),\nonumber \\ 2246 \vartheta_{t+1} & = & \vartheta_{t}+\omega_{t}\Delta t.\nonumber 2280 2247 \end{eqnarray} 2281 2248 … … 3259 3226 \begin_layout Subsection 3260 3227 Přiblížení metody vysokofrekvenční injektáží 3228 \begin_inset CommandInset label 3229 LatexCommand label 3230 name "sub:Priblizeni-metody-vysokofrekvenc" 3231 3232 \end_inset 3233 3234 3261 3235 \end_layout 3262 3236 … … 3296 3270 \begin_inset Formula 3297 3271 \[ 3298 u_{d}^{inj}=A_{inj}\cos\left(\omega_{inj}t\right) 3272 u_{d}^{inj}=A_{inj}\cos\left(\omega_{inj}t\right), 3299 3273 \] 3300 3274 … … 3372 3346 \begin_inset Formula 3373 3347 \begin{eqnarray} 3374 \frac{di_{d}}{dt} & = & \frac{1}{L_{d}}u_{d} \nonumber \\3375 \frac{di_{q}}{dt} & = & \frac{1}{L_{q}}u_{q} \label{eq:inj-hf-model}3348 \frac{di_{d}}{dt} & = & \frac{1}{L_{d}}u_{d},\nonumber \\ 3349 \frac{di_{q}}{dt} & = & \frac{1}{L_{q}}u_{q}.\label{eq:inj-hf-model} 3376 3350 \end{eqnarray} 3377 3351 … … 3395 3369 3396 3370 . 3397 \end_layout 3398 3399 \begin_layout Standard 3400 Průběh injektáže je pak následující: 3401 \end_layout 3402 3403 \begin_layout Enumerate 3404 injektování vysokofrekvenčního signálu do estimované osy 3371 Průběh injektáže je pak následující: 3372 \end_layout 3373 3374 \begin_layout Standard 3375 Nejdříve je injektován vysokofrekvenční signál do estimované osy 3405 3376 \begin_inset Formula $d$ 3406 3377 \end_inset … … 3413 3384 \begin_inset Formula 3414 3385 \begin{eqnarray*} 3415 \tilde{u}_{\hat{d}} & = & u_{\hat{d}}+A_{inj}\cos\left(\omega_{inj}t\right) \\3416 \tilde{u}_{\hat{q}} & = & u_{\hat{q}} 3386 \tilde{u}_{\hat{d}} & = & u_{\hat{d}}+A_{inj}\cos\left(\omega_{inj}t\right),\\ 3387 \tilde{u}_{\hat{q}} & = & u_{\hat{q}}, 3417 3388 \end{eqnarray*} 3418 3389 … … 3427 3398 \end_inset 3428 3399 3429 řídící zásah s injektáží 3430 \end_layout 3431 3432 \begin_layout Enumerate 3433 provedeme transformaci z estimovaného rotorového 3400 řídící zásah s injektáží. 3401 Následně provedeme transformaci z estimovaného rotorového 3434 3402 \begin_inset Formula $d-q$ 3435 3403 \end_inset … … 3453 3421 \begin_inset Formula 3454 3422 \begin{eqnarray*} 3455 \tilde{u}_{\alpha} & = & u_{\alpha}+A_{inj}\cos\left(\omega_{inj}t\right)\cos\hat{\vartheta} \\3456 \tilde{u}_{\beta} & = & u_{\beta}+A_{inj}\cos\left(\omega_{inj}t\right)\sin\hat{\vartheta} 3423 \tilde{u}_{\alpha} & = & u_{\alpha}+A_{inj}\cos\left(\omega_{inj}t\right)\cos\hat{\vartheta},\\ 3424 \tilde{u}_{\beta} & = & u_{\beta}+A_{inj}\cos\left(\omega_{inj}t\right)\sin\hat{\vartheta}, 3457 3425 \end{eqnarray*} 3458 3426 … … 3468 3436 \end_inset 3469 3437 3470 3471 \end_layout 3472 3473 \begin_layout Enumerate 3474 řídící zásahy 3438 . 3439 3440 \end_layout 3441 3442 \begin_layout Standard 3443 Řídící zásahy 3475 3444 \begin_inset Formula $\tilde{u}_{\alpha\beta}$ 3476 3445 \end_inset … … 3503 3472 \begin_inset Formula 3504 3473 \begin{eqnarray*} 3505 \tilde{u}_{d} & = & u_{d}+A_{inj}\cos\left(\omega_{inj}t\right)\cos\hat{\vartheta}\cos\vartheta+A_{inj}\cos\left(\omega_{inj}t\right)\sin\hat{\vartheta}\sin\vartheta \\3506 \tilde{u}_{q} & = & u_{q}-A_{inj}\cos\left(\omega_{inj}t\right)\cos\hat{\vartheta}\sin\vartheta+A_{inj}\cos\left(\omega_{inj}t\right)\sin\hat{\vartheta}\cos\vartheta 3474 \tilde{u}_{d} & = & u_{d}+A_{inj}\cos\left(\omega_{inj}t\right)\cos\hat{\vartheta}\cos\vartheta+A_{inj}\cos\left(\omega_{inj}t\right)\sin\hat{\vartheta}\sin\vartheta,\\ 3475 \tilde{u}_{q} & = & u_{q}-A_{inj}\cos\left(\omega_{inj}t\right)\cos\hat{\vartheta}\sin\vartheta+A_{inj}\cos\left(\omega_{inj}t\right)\sin\hat{\vartheta}\cos\vartheta, 3507 3476 \end{eqnarray*} 3508 3477 … … 3521 3490 \end_inset 3522 3491 3523 a nikoliv v estimované 3524 \end_layout 3525 3526 \begin_layout Enumerate 3527 řídící zásahy 3492 a nikoliv v estimované. 3493 Transformované řízení 3528 3494 \begin_inset Formula $\tilde{u}_{dq}$ 3529 3495 \end_inset … … 3548 3514 \begin_inset Formula 3549 3515 \begin{eqnarray*} 3550 \tilde{i}_{d} & = & i_{d}+\frac{A_{inj}}{L_{d}\omega_{inj}}\sin\left(\omega_{inj}t\right)\cos\theta \\3551 \tilde{i}_{q} & = & i_{q}-\frac{A_{inj}}{L_{q}\omega_{inj}}\sin\left(\omega_{inj}t\right)\sin\theta 3516 \tilde{i}_{d} & = & i_{d}+\frac{A_{inj}}{L_{d}\omega_{inj}}\sin\left(\omega_{inj}t\right)\cos\theta,\\ 3517 \tilde{i}_{q} & = & i_{q}-\frac{A_{inj}}{L_{q}\omega_{inj}}\sin\left(\omega_{inj}t\right)\sin\theta, 3552 3518 \end{eqnarray*} 3553 3519 … … 3566 3532 \end_inset 3567 3533 3568 a případné integrační konstanty 3569 \end_layout 3570 3571 \begin_layout Enumerate3572 návrh systému předpokládá měření proudů ve statorových souřadnicích a tedy3534 a případné integrační konstanty. 3535 \end_layout 3536 3537 \begin_layout Standard 3538 Návrh systému předpokládá měření proudů ve statorových souřadnicích a tedy 3573 3539 je nutné provést transformaci ( 3574 3540 \begin_inset CommandInset ref … … 3585 3551 \begin_inset Formula 3586 3552 \begin{eqnarray*} 3587 \tilde{i}_{\alpha} & = & i_{\alpha}+\frac{A_{inj}}{\omega_{inj}}\sin\left(\omega_{inj}t\right)\left(\frac{\cos\theta\cos\vartheta}{L_{d}}+\frac{\sin\theta\sin\vartheta}{L_{q}}\right) \\3588 \tilde{i}_{\beta} & = & i_{\beta}+\frac{A_{inj}}{\omega_{inj}}\sin\left(\omega_{inj}t\right)\left(\frac{\cos\theta\sin\vartheta}{L_{d}}-\frac{\sin\theta\cos\vartheta}{L_{q}}\right) 3553 \tilde{i}_{\alpha} & = & i_{\alpha}+\frac{A_{inj}}{\omega_{inj}}\sin\left(\omega_{inj}t\right)\left(\frac{\cos\theta\cos\vartheta}{L_{d}}+\frac{\sin\theta\sin\vartheta}{L_{q}}\right),\\ 3554 \tilde{i}_{\beta} & = & i_{\beta}+\frac{A_{inj}}{\omega_{inj}}\sin\left(\omega_{inj}t\right)\left(\frac{\cos\theta\sin\vartheta}{L_{d}}-\frac{\sin\theta\cos\vartheta}{L_{q}}\right), 3589 3555 \end{eqnarray*} 3590 3556 … … 3599 3565 \end_inset 3600 3566 3601 3602 \end_layout 3603 3604 \begin_layout Enumerate 3605 dále je ještě třeba převést proudy pomocí transformace ( 3567 . 3568 Dále je ještě třeba převést proudy pomocí transformace ( 3606 3569 \begin_inset CommandInset ref 3607 3570 LatexCommand ref … … 3617 3580 \begin_inset Formula 3618 3581 \begin{eqnarray*} 3619 \tilde{i}_{\hat{d}} & = & i_{\hat{d}}+\frac{A_{inj}}{\omega_{inj}}\sin\left(\omega_{inj}t\right)\left(\frac{\cos^{2}\theta}{L_{d}}+\frac{\sin^{2}\theta}{L_{q}}\right) \\3620 \tilde{i}_{\hat{q}} & = & i_{\hat{q}}+\frac{A_{inj}}{\omega_{inj}}\sin\left(\omega_{inj}t\right)\left(\frac{\sin\theta\cos\theta}{L_{d}}-\frac{\sin\theta\cos\theta}{L_{q}}\right) 3582 \tilde{i}_{\hat{d}} & = & i_{\hat{d}}+\frac{A_{inj}}{\omega_{inj}}\sin\left(\omega_{inj}t\right)\left(\frac{\cos^{2}\theta}{L_{d}}+\frac{\sin^{2}\theta}{L_{q}}\right),\\ 3583 \tilde{i}_{\hat{q}} & = & i_{\hat{q}}+\frac{A_{inj}}{\omega_{inj}}\sin\left(\omega_{inj}t\right)\left(\frac{\sin\theta\cos\theta}{L_{d}}-\frac{\sin\theta\cos\theta}{L_{q}}\right). 3621 3584 \end{eqnarray*} 3622 3585 … … 3626 3589 \end_layout 3627 3590 3628 \begin_layout Enumerate3629 následuje izolování modulovaného vysokofrekvenčního signálu na frekvenci3591 \begin_layout Standard 3592 Následuje izolování modulovaného vysokofrekvenčního signálu na frekvenci 3630 3593 3631 3594 \begin_inset Formula $\omega_{inj}$ … … 3643 3606 \begin_inset Formula 3644 3607 \[ 3645 i_{q}^{inj}=\frac{A_{inj}}{\omega_{inj}}\sin\left(\omega_{inj}t\right)\sin\theta\cos\theta\left(\frac{1}{L_{d}}-\frac{1}{L_{q}}\right) 3608 i_{q}^{inj}=\frac{A_{inj}}{\omega_{inj}}\sin\left(\omega_{inj}t\right)\sin\theta\cos\theta\left(\frac{1}{L_{d}}-\frac{1}{L_{q}}\right), 3646 3609 \] 3647 3610 … … 3655 3618 \begin_inset Formula 3656 3619 \begin{equation} 3657 \frac{A_{inj}\left(L_{q}-L_{d}\right)}{2\omega_{inj}L_{d}L_{q}}\sin2\theta \label{eq:inj-modul-signal}3620 \frac{A_{inj}\left(L_{q}-L_{d}\right)}{2\omega_{inj}L_{d}L_{q}}\sin2\theta.\label{eq:inj-modul-signal} 3658 3621 \end{equation} 3659 3622 3660 3623 \end_inset 3661 3624 3662 tento výsledek lze nalézt například v (3625 Tento výsledek lze nalézt například v ( 3663 3626 \series bold 3664 3627 citace - ale bohužel všede jsem to našel se 4 místo 2 … … 3901 3864 \begin_inset Formula 3902 3865 \begin{eqnarray*} 3903 u_{\alpha} & = & V\cos(2\pi ft) \\3904 u_{\beta} & = & V\sin(2\pi ft) 3866 u_{\alpha} & = & V\cos(2\pi ft),\\ 3867 u_{\beta} & = & V\sin(2\pi ft). 3905 3868 \end{eqnarray*} 3906 3869 … … 4026 3989 \begin_inset Formula 4027 3990 \[ 4028 x_{t}=\mathrm{PI}\left(e_{t},K_{p},K_{i}\right)=K_{p}e_{t}+K_{i}\intop_{0}^{t}e_{\tau}d\tau 3991 x_{t}=\mathrm{PI}\left(e_{t},K_{p},K_{i}\right)=K_{p}e_{t}+K_{i}\intop_{0}^{t}e_{\tau}d\tau. 4029 3992 \] 4030 3993 … … 4034 3997 \begin_inset Formula 4035 3998 \[ 4036 x_{t}=\mathrm{PI}\left(e_{t},K_{p},K_{i}\right)=K_{p}e_{t}+K_{i}\sum_{k=0}^{t}e_{k} 3999 x_{t}=\mathrm{PI}\left(e_{t},K_{p},K_{i}\right)=K_{p}e_{t}+K_{i}\sum_{k=0}^{t}e_{k}. 4037 4000 \] 4038 4001 … … 4128 4091 \begin_inset Formula 4129 4092 \[ 4130 \omega_{t+1}\text{=}d\omega_{t}+ei_{q,t} 4093 \omega_{t+1}\text{=}d\omega_{t}+ei_{q,t}, 4131 4094 \] 4132 4095 … … 4148 4111 \begin_inset Formula 4149 4112 \[ 4150 \overline{\omega}-d\omega=ei_{q} 4113 \overline{\omega}-d\omega=ei_{q}. 4151 4114 \] 4152 4115 … … 4160 4123 \begin_inset Formula 4161 4124 \[ 4162 \overline{i_{q}}=\mathrm{PI}(\overline{\omega}-\omega,K_{p,i},K_{i,i}) 4125 \overline{i_{q}}=\mathrm{PI}(\overline{\omega}-\omega,K_{p,i},K_{i,i}). 4163 4126 \] 4164 4127 … … 4203 4166 \begin_inset Formula 4204 4167 \begin{eqnarray*} 4205 i_{d,t+1} & = & ai_{d,t}+\Delta t\cdot i_{q,t}\omega_{t}+cu_{d,t} \\4206 i_{q,t+1} & = & ai_{q,t}-\Delta t\cdot i_{d,t}\omega_{t}-b\omega_{t}+cu_{q,t} 4168 i_{d,t+1} & = & ai_{d,t}+\Delta t\cdot i_{q,t}\omega_{t}+cu_{d,t},\\ 4169 i_{q,t+1} & = & ai_{q,t}-\Delta t\cdot i_{d,t}\omega_{t}-b\omega_{t}+cu_{q,t}, 4207 4170 \end{eqnarray*} 4208 4171 … … 4229 4192 \begin_inset Formula 4230 4193 \begin{eqnarray*} 4231 -ai_{d} & = & cu_{d} \\4232 \overline{i_{q}}-ai_{q} & = & cu_{q} 4194 -ai_{d} & = & cu_{d},\\ 4195 \overline{i_{q}}-ai_{q} & = & cu_{q}. 4233 4196 \end{eqnarray*} 4234 4197 … … 4563 4526 \begin_layout Subsection 4564 4527 Úloha duálního řízení 4528 \begin_inset CommandInset label 4529 LatexCommand label 4530 name "sub:uloha-dualniho-rizeni" 4531 4532 \end_inset 4533 4534 4565 4535 \end_layout 4566 4536 … … 4586 4556 \begin_inset Formula 4587 4557 \begin{eqnarray*} 4588 x_{t+1} & = & f_{t}\left(x_{t},p_{t},u_{t},\xi_{t}\right),\quad t=0,1,\ldots,T-1 \\4589 p{}_{t+1} & = & \upsilon_{t}\left(p_{t},\varepsilon_{t}\right) \\4590 y_{t} & = & h_{t}\left(x_{t},\eta_{t}\right) 4558 x_{t+1} & = & f_{t}\left(x_{t},p_{t},u_{t},\xi_{t}\right),\quad t=0,1,\ldots,T-1,\\ 4559 p{}_{t+1} & = & \upsilon_{t}\left(p_{t},\varepsilon_{t}\right),\\ 4560 y_{t} & = & h_{t}\left(x_{t},\eta_{t}\right), 4591 4561 \end{eqnarray*} 4592 4562 … … 4677 4647 \begin_inset Formula 4678 4648 \begin{equation} 4679 J=\mathrm{\mathbf{E}}\left\{ \sum_{t=0}^{T-1}g_{t+1}\left(x_{t+1},u_{t}\right)\right\} \label{eq:dclossfunc}4649 J=\mathrm{\mathbf{E}}\left\{ \sum_{t=0}^{T-1}g_{t+1}\left(x_{t+1},u_{t}\right)\right\} ,\label{eq:dclossfunc} 4680 4650 \end{equation} 4681 4651 … … 4752 4722 \begin_inset Formula 4753 4723 \begin{eqnarray*} 4754 J_{T-1}\left(I_{T-1}\right) & = & \min_{u_{T-1}\in U_{T-1}}\mathrm{\mathbf{E}}\left\{ g_{T}\left(x_{T},u_{T-1}\right)\mid I_{T-1}\right\}\\4755 J_{t}\left(I_{t}\right) & = & \min_{u_{t}\in U_{t}}\mathrm{\mathbf{E}}\left\{ g_{t+1}\left(x_{t+1},u_{t}\right)+J_{t+1}\left(I_{t+1}\right)\mid I_{t}\right\} 4724 V_{T-1}\left(I_{T-1}\right) & = & \min_{u_{T-1}\in U_{T-1}}\mathrm{\mathbf{E}}\left\{ g_{T}\left(x_{T},u_{T-1}\right)\mid I_{T-1}\right\} ,\\ 4725 V_{t}\left(I_{t}\right) & = & \min_{u_{t}\in U_{t}}\mathrm{\mathbf{E}}\left\{ g_{t+1}\left(x_{t+1},u_{t}\right)+V_{t+1}\left(I_{t+1}\right)\mid I_{t}\right\} , 4756 4726 \end{eqnarray*} 4757 4727 … … 4760 4730 pro 4761 4731 \begin_inset Formula $t=T-2,T-3,\ldots,0$ 4732 \end_inset 4733 4734 . 4735 Funkce 4736 \begin_inset Formula $V$ 4737 \end_inset 4738 4739 vystupující v předchozích rovnicích je nazývána jako 4740 \emph on 4741 Bellmanova 4742 \emph default 4743 funkce 4744 \begin_inset CommandInset citation 4745 LatexCommand cite 4746 key "utiaBDM2005" 4747 4762 4748 \end_inset 4763 4749 … … 4777 4763 Hlavními komplikacemi jsou jednak výpočet střední hodnoty a minimalizace, 4778 4764 ale hlavně problémy spojené s funkcí 4779 \begin_inset Formula $ J$4780 \end_inset 4781 4782 . 4783 Funkce4784 \begin_inset Formula $ J$4765 \begin_inset Formula $V$ 4766 \end_inset 4767 4768 . 4769 Bellmanova funkce 4770 \begin_inset Formula $V$ 4785 4771 \end_inset 4786 4772 … … 4790 4776 Tuto funkci je navíc třeba uchovávat mezi jednotlivými časovými kroky v 4791 4777 její plné reprezentaci jako funkce, ne pouze její hodnotu. 4792 ( 4793 \series bold 4794 možná citace 4795 \series default 4796 ) 4778 4779 \begin_inset CommandInset citation 4780 LatexCommand cite 4781 key "utiaBDM2005" 4782 4783 \end_inset 4784 4785 4797 4786 \end_layout 4798 4787 … … 4910 4899 \begin_inset Formula 4911 4900 \[ 4912 J_{t}\left(I_{t},\rho_{t}\right)=\mathrm{\mathbf{E}}_{\rho_{t}}\left\{ \sum_{i=t}^{T-1}g_{t+1}\left(x_{i+1},u_{i}\right)\mid I_{k}\right\} 4901 J_{t}\left(I_{t},\rho_{t}\right)=\mathrm{\mathbf{E}}_{\rho_{t}}\left\{ \sum_{i=t}^{T-1}g_{t+1}\left(x_{i+1},u_{i}\right)\mid I_{k}\right\} . 4913 4902 \] 4914 4903 … … 4927 4916 \begin_inset Formula 4928 4917 \[ 4929 \rho_{t}=\mathrm{p}\left(x_{t+i},p_{t+i}\mid I_{t+i}\right) 4918 \rho_{t}=\mathrm{p}\left(x_{t+i},p_{t+i}\mid I_{t+i}\right), 4930 4919 \] 4931 4920 … … 4961 4950 \begin_inset Formula 4962 4951 \[ 4963 \rho_{t}=\left\{ \mathrm{p}\left(x_{t+i},p_{t+i}\mid I_{t+i}\right)=\mathrm{p}\left(x_{t+i},p_{t+i}\mid I_{0}\right),i=0,\ldots,T-t-1\right\} 4952 \rho_{t}=\left\{ \mathrm{p}\left(x_{t+i},p_{t+i}\mid I_{t+i}\right)=\mathrm{p}\left(x_{t+i},p_{t+i}\mid I_{0}\right),i=0,\ldots,T-t-1\right\} . 4964 4953 \] 4965 4954 … … 5005 4994 \begin_inset Formula 5006 4995 \[ 5007 \rho_{t}=\left\{ \mathrm{p}\left(x_{t+i},p_{t+i}\mid I_{t+i}\right)=\mathrm{p}\left(x_{t+i},p_{t+i}\mid I_{t}\right),i=0,\ldots,T-t-1\right\} 4996 \rho_{t}=\left\{ \mathrm{p}\left(x_{t+i},p_{t+i}\mid I_{t+i}\right)=\mathrm{p}\left(x_{t+i},p_{t+i}\mid I_{t}\right),i=0,\ldots,T-t-1\right\} . 5008 4997 \] 5009 4998 … … 5023 5012 \begin{align*} 5024 5013 \rho_{t} & =\left\{ \mathrm{p}\left(x_{t+i},p_{t+i}\mid I_{t+i}\right)\right.\\ 5025 = & \left.\delta\left(x_{t+i}-\hat{x}_{t+i}\right)\delta\left(p_{t+i}-\hat{p}_{t+i}\right),i=0,\ldots,T-t-1\right\} 5014 = & \left.\delta\left(x_{t+i}-\hat{x}_{t+i}\right)\delta\left(p_{t+i}-\hat{p}_{t+i}\right),i=0,\ldots,T-t-1\right\} , 5026 5015 \end{align*} 5027 5016 … … 5081 5070 \begin{align*} 5082 5071 \rho_{t} & =\left\{ \mathrm{p}\left(z_{1,t+i},z_{2,t+i}\mid I_{t+i}\right)\right.\\ 5083 = & \left.\delta\left(z_{1,t+i}-\hat{z}_{1,t+i}\right)\mathrm{p}\left(z_{2,t+i}\mid I_{t}\right),i=0,\ldots,T-t-1\right\} 5072 = & \left.\delta\left(z_{1,t+i}-\hat{z}_{1,t+i}\right)\mathrm{p}\left(z_{2,t+i}\mid I_{t}\right),i=0,\ldots,T-t-1\right\} , 5084 5073 \end{align*} 5085 5074 … … 5180 5169 5181 5170 . 5182 Nyní definuj eme vektor5171 Nyní definujme vektor 5183 5172 \emph on 5184 5173 hyperstavu … … 5373 5362 \begin_inset Formula 5374 5363 \[ 5375 P=\mathtt{E}\left\{ \left[g(x)-\theta\right]\left[g(x)-\theta\right]^{T}\right\} \geq J^{-1} 5364 P=\mathtt{E}\left\{ \left[g(x)-\theta\right]\left[g(x)-\theta\right]^{T}\right\} \geq J^{-1}, 5376 5365 \] 5377 5366 … … 5389 5378 \begin_inset Formula 5390 5379 \[ 5391 J_{ij}=\mathtt{E}\left[-\frac{\partial^{2}\log p_{x,\theta}(X,\Theta)}{\partial\Theta_{i}\partial\Theta_{j}}\right] 5380 J_{ij}=\mathtt{E}\left[-\frac{\partial^{2}\log p_{x,\theta}(X,\Theta)}{\partial\Theta_{i}\partial\Theta_{j}}\right], 5392 5381 \] 5393 5382 … … 5431 5420 \begin_inset Formula 5432 5421 \begin{eqnarray} 5433 x_{t+1} & = & f_{t}(x_{t},w_{t}) \nonumber \\5434 z_{t} & = & h_{t}(x_{t},v_{t}) \label{eq:PCRB-system}5422 x_{t+1} & = & f_{t}(x_{t},w_{t}),\nonumber \\ 5423 z_{t} & = & h_{t}(x_{t},v_{t}),\label{eq:PCRB-system} 5435 5424 \end{eqnarray} 5436 5425 … … 5484 5473 \begin_inset Formula 5485 5474 \[ 5486 J_{t+1}=D_{t}^{22}-D_{t}^{21}\left(J_{t}+D_{t}^{11}\right)^{-1}D_{t}^{12} 5475 J_{t+1}=D_{t}^{22}-D_{t}^{21}\left(J_{t}+D_{t}^{11}\right)^{-1}D_{t}^{12}, 5487 5476 \] 5488 5477 … … 5499 5488 \begin_inset Formula 5500 5489 \begin{eqnarray} 5501 D_{t}^{11} & = & \mathtt{E}\left\{ -\triangle_{x_{t}}^{x_{t}}\log p(x_{t+1}\mid x_{t})\right\} \nonumber \\5502 D_{t}^{12} & = & \mathtt{E}\left\{ -\triangle_{x_{t}}^{x_{t+1}}\log p(x_{t+1}\mid x_{t})\right\} \label{eq:PCRB-rovnice-pro-D}\\5503 D_{t}^{21} & = & \mathtt{E}\left\{ -\triangle_{x_{t+1}}^{x_{t}}\log p(x_{t+1}\mid x_{t})\right\} =\left(D_{t}^{12}\right)^{T} \nonumber \\5504 D_{t}^{22} & = & \mathtt{E}\left\{ -\triangle_{x_{t+1}}^{x_{t+1}}\log p(x_{t+1}\mid x_{t})\right\} +\mathtt{E}\left\{ -\triangle_{x_{t+1}}^{x_{t+1}}\log p(z_{t+1}\mid x_{t+1})\right\} \nonumber5490 D_{t}^{11} & = & \mathtt{E}\left\{ -\triangle_{x_{t}}^{x_{t}}\log p(x_{t+1}\mid x_{t})\right\} ,\nonumber \\ 5491 D_{t}^{12} & = & \mathtt{E}\left\{ -\triangle_{x_{t}}^{x_{t+1}}\log p(x_{t+1}\mid x_{t})\right\} ,\label{eq:PCRB-rovnice-pro-D}\\ 5492 D_{t}^{21} & = & \mathtt{E}\left\{ -\triangle_{x_{t+1}}^{x_{t}}\log p(x_{t+1}\mid x_{t})\right\} =\left(D_{t}^{12}\right)^{T},\nonumber \\ 5493 D_{t}^{22} & = & \mathtt{E}\left\{ -\triangle_{x_{t+1}}^{x_{t+1}}\log p(x_{t+1}\mid x_{t})\right\} +\mathtt{E}\left\{ -\triangle_{x_{t+1}}^{x_{t+1}}\log p(z_{t+1}\mid x_{t+1})\right\} .\nonumber 5505 5494 \end{eqnarray} 5506 5495 … … 5528 5517 \begin_inset Formula 5529 5518 \begin{eqnarray} 5530 x_{t+1} & = & f_{t}(x_{t})+w_{t} \nonumber \\5519 x_{t+1} & = & f_{t}(x_{t})+w_{t},\nonumber \\ 5531 5520 z_{t} & = & h_{t}(x_{t})+v_{t}\label{eq:PCRB-system-adsum} 5532 5521 \end{eqnarray} … … 5563 5552 \begin_inset Formula 5564 5553 \begin{eqnarray} 5565 D_{t}^{11} & = & \mathtt{E}\left\{ \left[\nabla_{x_{t}}f_{t}^{T}(x_{t})\right]Q_{t}^{-1}\left[\nabla_{x_{t}}f_{t}^{T}(x_{t})\right]^{T}\right\} \nonumber \\5566 D_{t}^{12} & = & -\mathtt{E}\left\{ \nabla_{x_{t}}f_{t}^{T}(x_{t})\right\} Q_{t}^{-1} \label{eq:PCRB-rovnice-pro-D-gauss}\\5567 D_{t}^{22} & = & Q_{t}^{-1}+\mathtt{E}\left\{ \left[\nabla_{x_{t+1}}h_{t+1}^{T}(x_{t+1})\right]R_{t+1}^{-1}\left[\nabla_{x_{t+1}}h_{t+1}^{T}(x_{t+1})\right]^{T}\right\} \nonumber5554 D_{t}^{11} & = & \mathtt{E}\left\{ \left[\nabla_{x_{t}}f_{t}^{T}(x_{t})\right]Q_{t}^{-1}\left[\nabla_{x_{t}}f_{t}^{T}(x_{t})\right]^{T}\right\} ,\nonumber \\ 5555 D_{t}^{12} & = & -\mathtt{E}\left\{ \nabla_{x_{t}}f_{t}^{T}(x_{t})\right\} Q_{t}^{-1},\label{eq:PCRB-rovnice-pro-D-gauss}\\ 5556 D_{t}^{22} & = & Q_{t}^{-1}+\mathtt{E}\left\{ \left[\nabla_{x_{t+1}}h_{t+1}^{T}(x_{t+1})\right]R_{t+1}^{-1}\left[\nabla_{x_{t+1}}h_{t+1}^{T}(x_{t+1})\right]^{T}\right\} .\nonumber 5568 5557 \end{eqnarray} 5569 5558 … … 5602 5591 5603 5592 \begin_layout Section 5593 \begin_inset CommandInset label 5594 LatexCommand label 5595 name "sec:LQG-obecne" 5596 5597 \end_inset 5598 5604 5599 Lineárně kvadraticky Gaussovské řízení 5605 5600 \end_layout … … 5688 5683 \begin_inset Formula 5689 5684 \begin{eqnarray*} 5690 x_{t} & = & f\left(x_{t-1},u_{t-1}\right)+w_{t-1} \\5691 y_{t} & = & h\left(x_{t}\right)+v_{t} 5685 x_{t} & = & f\left(x_{t-1},u_{t-1}\right)+w_{t-1},\\ 5686 y_{t} & = & h\left(x_{t}\right)+v_{t}, 5692 5687 \end{eqnarray*} 5693 5688 … … 5793 5788 \begin_inset Formula 5794 5789 \begin{eqnarray} 5795 \overline{\hat{x}}_{t} & = & f\left(\hat{x}_{t-1},u_{t-1},0\right) \nonumber \\5796 \overline{P}_{t} & = & A_{t}P_{t-1}A_{t}^{T}+Q_{t-1} \label{eq:EKF-rovnice-time-upd}5790 \overline{\hat{x}}_{t} & = & f\left(\hat{x}_{t-1},u_{t-1},0\right),\nonumber \\ 5791 \overline{P}_{t} & = & A_{t}P_{t-1}A_{t}^{T}+Q_{t-1}.\label{eq:EKF-rovnice-time-upd} 5797 5792 \end{eqnarray} 5798 5793 … … 5815 5810 \begin_inset Formula 5816 5811 \begin{eqnarray} 5817 K_{t} & = & \overline{P}_{t}C_{t}^{T}\left(C_{t}\overline{P}_{t}C_{t}^{T}+R_{t}\right)^{-1} \nonumber \\5818 \hat{x}_{t} & = & \overline{\hat{x}}_{t}+K_{t}\left(y_{t}-h\left(\overline{\hat{x}}_{t},0\right)\right) \label{eq:EKF-rovnice-data-upd}\\5819 P_{t} & = & \left(I-K_{t}C_{t}\right)\overline{P}_{t} \nonumber5812 K_{t} & = & \overline{P}_{t}C_{t}^{T}\left(C_{t}\overline{P}_{t}C_{t}^{T}+R_{t}\right)^{-1},\nonumber \\ 5813 \hat{x}_{t} & = & \overline{\hat{x}}_{t}+K_{t}\left(y_{t}-h\left(\overline{\hat{x}}_{t},0\right)\right),\label{eq:EKF-rovnice-data-upd}\\ 5814 P_{t} & = & \left(I-K_{t}C_{t}\right)\overline{P}_{t},\nonumber 5820 5815 \end{eqnarray} 5821 5816 5822 5817 \end_inset 5823 5818 5824 Pro úplnost je ještě třeba dodat počáteční apriorní odhady 5819 kde 5820 \begin_inset Formula $I$ 5821 \end_inset 5822 5823 značí jednotkovou matici vhodného rozměru. 5824 Pro úplnost je ještě třeba dodat počáteční apriorní odhady 5825 5825 \begin_inset Formula $\hat{x}_{0}$ 5826 5826 \end_inset … … 5864 5864 \begin_inset Formula 5865 5865 \begin{equation} 5866 x_{t+1}=A_{t}x_{t}+B_{t}u_{t}+w_{t},\quad t=0,1,\ldots,T-1 \label{eq:lq-obecny-lin-system}5866 x_{t+1}=A_{t}x_{t}+B_{t}u_{t}+w_{t},\quad t=0,1,\ldots,T-1,\label{eq:lq-obecny-lin-system} 5867 5867 \end{equation} 5868 5868 … … 5903 5903 \begin_inset Formula 5904 5904 \begin{equation} 5905 \mathbf{E}\left\{ x_{T}^{T}Q_{T}x_{T}+\sum_{t=0}^{T-1}\left(x_{t}^{T}Q_{t}x_{t}+u_{t}^{T}R_{t}u_{t}\right)\right\} \label{eq:lq-adit-kv-ztrata}5905 \mathbf{E}\left\{ x_{T}^{T}Q_{T}x_{T}+\sum_{t=0}^{T-1}\left(x_{t}^{T}Q_{t}x_{t}+u_{t}^{T}R_{t}u_{t}\right)\right\} ,\label{eq:lq-adit-kv-ztrata} 5906 5906 \end{equation} 5907 5907 … … 5962 5962 \begin_inset Formula 5963 5963 \[ 5964 \mu_{t}(I_{t})=L_{t}\mathrm{\mathbf{E}}\left\{ x_{t}\mid I_{t}\right\} 5964 \mu_{t}(I_{t})=L_{t}\mathrm{\mathbf{E}}\left\{ x_{t}\mid I_{t}\right\} , 5965 5965 \] 5966 5966 … … 5974 5974 \begin_inset Formula 5975 5975 \begin{equation} 5976 L_{t}=-\left(R_{t}+B_{t}^{T}K_{t+1}B_{t}\right)^{-1}B_{t}^{T}K_{t+1}A_{t} \label{eq:riccati-lqg-matice-L}5976 L_{t}=-\left(R_{t}+B_{t}^{T}K_{t+1}B_{t}\right)^{-1}B_{t}^{T}K_{t+1}A_{t},\label{eq:riccati-lqg-matice-L} 5977 5977 \end{equation} 5978 5978 … … 5986 5986 \begin_inset Formula 5987 5987 \begin{eqnarray} 5988 K_{T} & = & Q_{T} \label{eq:riccati-lqg-matice-K}\\5989 K_{t} & = & A_{t}^{T}\left(K_{t+1}-K_{t+1}B_{t}\left(R_{t}+B_{t}^{T}K_{t+1}B_{t}\right)^{-1}B_{t}^{T}K_{t+1}\right)A_{t}+Q_{t} \nonumber5988 K_{T} & = & Q_{T},\label{eq:riccati-lqg-matice-K}\\ 5989 K_{t} & = & A_{t}^{T}\left(K_{t+1}-K_{t+1}B_{t}\left(R_{t}+B_{t}^{T}K_{t+1}B_{t}\right)^{-1}B_{t}^{T}K_{t+1}\right)A_{t}+Q_{t}.\nonumber 5990 5990 \end{eqnarray} 5991 5991 … … 6041 6041 \begin_inset Formula 6042 6042 \begin{equation} 6043 x_{t+1}^{T}Q_{t}x_{t+1}+u_{t}^{T}R_{t}u_{t}=x_{t+1}^{T}\sqrt{Q_{t}}^{T}\sqrt{Q_{t}}x_{t+1}+u_{t}^{T}\sqrt{R_{t}}^{T}\sqrt{R_{t}}u_{t} \label{eq:lq-kv-ztrata-prepis-odmoc}6043 x_{t+1}^{T}Q_{t}x_{t+1}+u_{t}^{T}R_{t}u_{t}=x_{t+1}^{T}\sqrt{Q_{t}}^{T}\sqrt{Q_{t}}x_{t+1}+u_{t}^{T}\sqrt{R_{t}}^{T}\sqrt{R_{t}}u_{t},\label{eq:lq-kv-ztrata-prepis-odmoc} 6044 6044 \end{equation} 6045 6045 … … 6067 6067 \begin_inset Formula 6068 6068 \begin{equation} 6069 x_{t+1}^{T}\sqrt{Q_{t}}^{T}\sqrt{Q_{t}}x_{t+1}+u_{t}^{T}\sqrt{R_{t}}^{T}\sqrt{R_{t}}u_{t}+x_{t+1}^{T}\sqrt{\Sigma_{t}}^{T}\sqrt{\Sigma_{t}}x_{t+1} \label{eq:lq-odm-ztrata-se-sigma}6069 x_{t+1}^{T}\sqrt{Q_{t}}^{T}\sqrt{Q_{t}}x_{t+1}+u_{t}^{T}\sqrt{R_{t}}^{T}\sqrt{R_{t}}u_{t}+x_{t+1}^{T}\sqrt{\Sigma_{t}}^{T}\sqrt{\Sigma_{t}}x_{t+1},\label{eq:lq-odm-ztrata-se-sigma} 6070 6070 \end{equation} 6071 6071 … … 6100 6100 u_{t}\\ 6101 6101 x_{t} 6102 \end{array}\right) \label{eq:lq-matic-ztrata-pro-qr}6102 \end{array}\right).\label{eq:lq-matic-ztrata-pro-qr} 6103 6103 \end{equation} 6104 6104 … … 6154 6154 R_{uu} & R_{ux}\\ 6155 6155 0 & R_{xx} 6156 \end{array}\right] 6156 \end{array}\right]. 6157 6157 \] 6158 6158 … … 6175 6175 R_{xx}x_{t} 6176 6176 \end{array}\right)\\ 6177 & = & \left(R_{uu}u_{t}+R_{ux}x_{t}\right)^{T}\left(R_{uu}u_{t}+R_{ux}x_{t}\right)+x_{t}^{T}R_{xx}^{T}R_{xx}x_{t} 6177 & = & \left(R_{uu}u_{t}+R_{ux}x_{t}\right)^{T}\left(R_{uu}u_{t}+R_{ux}x_{t}\right)+x_{t}^{T}R_{xx}^{T}R_{xx}x_{t}, 6178 6178 \end{eqnarray*} 6179 6179 … … 6195 6195 \begin_inset Formula 6196 6196 \[ 6197 u_{t}=-R_{uu}^{-1}R_{ux}x_{t} 6197 u_{t}=-R_{uu}^{-1}R_{ux}x_{t}. 6198 6198 \] 6199 6199 … … 6232 6232 \begin_layout Section 6233 6233 Úloha řízení PMSM 6234 \begin_inset CommandInset label 6235 LatexCommand label 6236 name "sec:uloha-rizeni-PMSM" 6237 6238 \end_inset 6239 6240 6234 6241 \end_layout 6235 6242 … … 6497 6504 \begin_inset Formula 6498 6505 \begin{eqnarray*} 6499 x_{t} & = & \left(i_{\alpha,t},i_{\beta,t},\omega_{t},\vartheta_{t}\right)^{T} \\6500 y_{t} & = & \left(y_{\alpha,t},y_{\beta,t}\right)^{T} 6506 x_{t} & = & \left(i_{\alpha,t},i_{\beta,t},\omega_{t},\vartheta_{t}\right)^{T},\\ 6507 y_{t} & = & \left(y_{\alpha,t},y_{\beta,t}\right)^{T}, 6501 6508 \end{eqnarray*} 6502 6509 … … 6513 6520 \begin_inset Formula 6514 6521 \begin{eqnarray*} 6515 x_{t+1} & = & f\left(x_{t},u_{t}\right)+w_{t} \\6516 y_{t} & = & h\left(x_{t}\right)+v_{t} 6522 x_{t+1} & = & f\left(x_{t},u_{t}\right)+w_{t},\\ 6523 y_{t} & = & h\left(x_{t}\right)+v_{t}, 6517 6524 \end{eqnarray*} 6518 6525 … … 6595 6602 -e\sin\hat{\vartheta}_{t-1} & e\cos\hat{\vartheta}_{t-1} & d & -e\left(\hat{i}_{\beta,t-1}\sin\hat{\vartheta}_{t-1}+\hat{i}_{\alpha,t-1}\cos\hat{\vartheta}_{t-1}\right)\\ 6596 6603 0 & 0 & \Delta t & 1 6597 \end{array}\right] \nonumber \\6604 \end{array}\right],\nonumber \\ 6598 6605 C_{t} & = & C=\left[\begin{array}{cccc} 6599 6606 1 & 0 & 0 & 0\\ 6600 6607 0 & 1 & 0 & 0 6601 \end{array}\right] \label{eq:matice-ekf-plnymodel-ls}6608 \end{array}\right].\label{eq:matice-ekf-plnymodel-ls} 6602 6609 \end{eqnarray} 6603 6610 … … 6661 6668 \begin_inset Formula 6662 6669 \begin{eqnarray*} 6663 x_{t} & = & \left(\omega_{t},\vartheta_{t}\right)^{T} \\6664 y_{t} & = & \left(i_{\alpha,t},i_{\beta,t}\right)^{T} 6670 x_{t} & = & \left(\omega_{t},\vartheta_{t}\right)^{T},\\ 6671 y_{t} & = & \left(i_{\alpha,t},i_{\beta,t}\right)^{T}. 6665 6672 \end{eqnarray*} 6666 6673 … … 6681 6688 d & -e\left(\hat{i}_{\beta,t-1}\sin\hat{\vartheta}_{t-1}+\hat{i}_{\alpha,t-1}\cos\hat{\vartheta}_{t-1}\right)\\ 6682 6689 \Delta t & 1 6683 \end{array}\right] \nonumber \\6690 \end{array}\right],\nonumber \\ 6684 6691 C_{t} & = & \left[\begin{array}{cc} 6685 6692 b\sin\hat{\vartheta}_{t-1} & b\hat{\omega}_{t-1}\cos\hat{\vartheta}_{t-1}\\ 6686 6693 -b\cos\hat{\vartheta}_{t-1} & b\hat{\omega}_{t-1}\sin\hat{\vartheta}_{t-1} 6687 \end{array}\right] \label{eq:matice-ekf-redukovanymodel-ls}6694 \end{array}\right].\label{eq:matice-ekf-redukovanymodel-ls} 6688 6695 \end{eqnarray} 6689 6696 … … 6713 6720 Q_{1}\\ 6714 6721 & Q_{2} 6715 \end{array}\right] 6722 \end{array}\right]. 6716 6723 \] 6717 6724 … … 6722 6729 \begin_inset Formula 6723 6730 \begin{eqnarray*} 6724 Q_{red} & = & Q_{2} \\6725 R_{red} & = & R+Q_{1} 6731 Q_{red} & = & Q_{2},\\ 6732 R_{red} & = & R+Q_{1}. 6726 6733 \end{eqnarray*} 6727 6734 … … 7060 7067 7061 7068 \begin_layout Plain Layout 7062 a) oříznutí pevnou mezí 7063 \begin_inset Formula $\frac{\pi^{2}}{3}$ 7064 \end_inset 7065 7066 (čárkovaně) 7069 a) pevná mez 7067 7070 \end_layout 7068 7071 … … 7073 7076 7074 7077 \begin_layout Plain Layout 7075 b) oříznut í pomocí oříznutého normálníhorozdělení7078 b) oříznuté normální rozdělení 7076 7079 \end_layout 7077 7080 … … 7090 7093 7091 7094 \begin_layout Plain Layout 7092 Hodnoty PCRB polohy7095 Srovnání metod omezování hodnoty PCRB polohy 7093 7096 \begin_inset Formula $\vartheta$ 7094 7097 \end_inset 7095 7098 7096 v závislosti na amplitudě injektovaného konstantního signálu (viz legenda). 7099 : První možností je oříznutí pevnou mezí 7100 \begin_inset Formula $\frac{\pi^{2}}{3}$ 7101 \end_inset 7102 7103 (znázorněna čárkovaně), druhou pak užití oříznutého normálního rozdělení. 7097 7104 7098 7105 \begin_inset CommandInset label … … 7142 7149 7143 7150 \begin_layout Standard 7144 Oříznuté normální rozdělení pro skalární váhodnou veličinu7151 Oříznuté normální rozdělení pro skalární náhodnou veličinu 7145 7152 \begin_inset Formula $x$ 7146 7153 \end_inset … … 7158 7165 \begin_inset Formula 7159 7166 \begin{eqnarray*} 7160 \hat{x} & = & \mu-\sqrt{r}\varphi(\mu,r) \\7161 \hat{x^{2}} & = & r+\mu\hat{x}-\sqrt{r}\kappa(\mu,r) 7167 \hat{x} & = & \mu-\sqrt{r}\varphi(\mu,r),\\ 7168 \hat{x^{2}} & = & r+\mu\hat{x}-\sqrt{r}\kappa(\mu,r), 7162 7169 \end{eqnarray*} 7163 7170 … … 7167 7174 \begin_inset Formula 7168 7175 \begin{eqnarray*} 7169 \varphi(\mu,r) & = & \frac{\sqrt{2}\left(\exp(-\beta^{2})-\exp(-\alpha^{2})\right)}{\sqrt{\pi}\left(\mathrm{erf}(\beta)-\mathrm{erf}(\alpha)\right)} \\7176 \varphi(\mu,r) & = & \frac{\sqrt{2}\left(\exp(-\beta^{2})-\exp(-\alpha^{2})\right)}{\sqrt{\pi}\left(\mathrm{erf}(\beta)-\mathrm{erf}(\alpha)\right)},\\ 7170 7177 \kappa(\mu,r) & = & \frac{\sqrt{2}\left(b\exp(-\beta^{2})-a\exp(-\alpha^{2})\right)}{\sqrt{\pi}\left(\mathrm{erf}(\beta)-\mathrm{erf}(\alpha)\right)} 7171 7178 \end{eqnarray*} … … 7176 7183 \begin_inset Formula 7177 7184 \begin{eqnarray*} 7178 \alpha & = & \frac{a-\mu}{\sqrt{2r}} \\7179 \beta & = & \frac{b-\mu}{\sqrt{2r}} 7185 \alpha & = & \frac{a-\mu}{\sqrt{2r}},\\ 7186 \beta & = & \frac{b-\mu}{\sqrt{2r}}. 7180 7187 \end{eqnarray*} 7181 7188 … … 7232 7239 \begin_inset Formula 7233 7240 \[ 7234 \kappa=\frac{2\sqrt{2}\pi\exp\left(-\frac{\pi^{2}}{2r}\right)}{2\sqrt{\pi}\mathrm{erf}\left(\frac{\pi}{\sqrt{2r}}\right)} 7241 \kappa=\frac{2\sqrt{2}\pi\exp\left(-\frac{\pi^{2}}{2r}\right)}{2\sqrt{\pi}\mathrm{erf}\left(\frac{\pi}{\sqrt{2r}}\right)}. 7235 7242 \] 7236 7243 … … 7244 7251 \begin_inset Formula 7245 7252 \[ 7246 \mathrm{Var}(x)=r-\sqrt{2\pi r}\frac{\exp\left(-\frac{\pi^{2}}{2r}\right)}{\mathrm{erf}\left(\frac{\pi}{\sqrt{2r}}\right)} 7253 \mathrm{Var}(x)=r-\sqrt{2\pi r}\frac{\exp\left(-\frac{\pi^{2}}{2r}\right)}{\mathrm{erf}\left(\frac{\pi}{\sqrt{2r}}\right)}. 7247 7254 \] 7248 7255 … … 7286 7293 \begin_inset Formula 7287 7294 \[ 7288 f\left(x\right)\cong f\left(x_{0}\right)+\frac{\partial f}{\partial x}\left(x_{0}\right)\left(x-x_{0}\right) 7295 f\left(x\right)\cong f\left(x_{0}\right)+\frac{\partial f}{\partial x}\left(x_{0}\right)\left(x-x_{0}\right), 7289 7296 \] 7290 7297 … … 7313 7320 \begin_inset Formula 7314 7321 \[ 7315 f\left(x\right)\cong Ax+\left(f\left(x_{0}\right)-Ax_{0}\right)=Ax+\gamma 7322 f\left(x\right)\cong Ax+\left(f\left(x_{0}\right)-Ax_{0}\right)=Ax+\gamma, 7316 7323 \] 7317 7324 … … 7356 7363 x\\ 7357 7364 1 7358 \end{array}\right) 7365 \end{array}\right), 7359 7366 \] 7360 7367 … … 7367 7374 A & \gamma\\ 7368 7375 0 & 1 7369 \end{array}\right] 7376 \end{array}\right], 7370 7377 \] 7371 7378 … … 7428 7435 a zavést substituci 7429 7436 \begin_inset Formula 7430 \ [7431 \psi_{t}=\omega_{t}-\overline{\omega}_{t} 7432 \ ]7437 \begin{equation} 7438 \psi_{t}=\omega_{t}-\overline{\omega}_{t}\label{eq:substituce-psi} 7439 \end{equation} 7433 7440 7434 7441 \end_inset … … 7483 7490 0 & 0 & 0 & 0 & 0\\ 7484 7491 0 & 0 & 0 & 0 & 0 7485 \end{array}\right] \label{eq:matice-Q-lq}7492 \end{array}\right],\label{eq:matice-Q-lq} 7486 7493 \end{equation} 7487 7494 … … 7518 7525 \begin_inset Formula 7519 7526 \begin{equation} 7520 \left|u_{k,t}\right|\leq U_{max} \label{eq:omezeni}7527 \left|u_{k,t}\right|\leq U_{max},\label{eq:omezeni} 7521 7528 \end{equation} 7522 7529 7523 7530 \end_inset 7524 7531 7525 jakoomezení na každou složku7532 tedy omezení na každou složku 7526 7533 \begin_inset Formula $k$ 7527 7534 \end_inset … … 7639 7646 \begin_inset Formula 7640 7647 \[ 7641 \left(u_{t}-u_{t-1}\right)^{T}S_{t}\left(u_{t}-u_{t-1}\right)=\left(u_{t}-u_{t-1}\right)^{T}\sqrt{S_{t}}^{T}\sqrt{S_{t}}\left(u_{t}-u_{t-1}\right) 7648 \left(u_{t}-u_{t-1}\right)^{T}S_{t}\left(u_{t}-u_{t-1}\right)=\left(u_{t}-u_{t-1}\right)^{T}\sqrt{S_{t}}^{T}\sqrt{S_{t}}\left(u_{t}-u_{t-1}\right). 7642 7649 \] 7643 7650 … … 7662 7669 u_{t-1}\\ 7663 7670 x_{t} 7664 \end{array}\right) 7671 \end{array}\right), 7665 7672 \] 7666 7673 … … 7679 7686 \sqrt{S_{t}} & -\sqrt{S_{t}} & 0\\ 7680 7687 \left(\sqrt{\Sigma_{t}}\right)_{3,\ldots,7}B_{t} & \left(\sqrt{\Sigma_{t}}\right)_{1,2} & \left(\sqrt{\Sigma_{t}}\right)_{3,\ldots,7}A_{t} 7681 \end{array}\right] 7688 \end{array}\right], 7682 7689 \] 7683 7690 … … 7716 7723 \begin_inset Formula 7717 7724 \[ 7718 \sum_{t=0}^{T-1}\left(\left(x_{t}-\overline{x}_{t}\right)^{T}Q\left(x_{t}-\overline{x}_{t}\right)+u_{t}^{T}Ru_{t}+\left(u_{t}-u_{t-1}\right)^{T}S\left(u_{t}-u_{t-1}\right)\right) 7725 \sum_{t=0}^{T-1}\left(\left(x_{t}-\overline{x}_{t}\right)^{T}Q\left(x_{t}-\overline{x}_{t}\right)+u_{t}^{T}Ru_{t}+\left(u_{t}-u_{t-1}\right)^{T}S\left(u_{t}-u_{t-1}\right)\right), 7719 7726 \] 7720 7727 7721 7728 \end_inset 7722 7729 7723 7730 kde řízení 7724 7731 \begin_inset Formula $u_{t}$ 7725 7732 \end_inset … … 7786 7793 \begin_inset Formula 7787 7794 \[ 7788 R^{dq}=diag(r_{d},r_{q}) 7795 R^{dq}=diag(r_{d},r_{q}). 7789 7796 \] 7790 7797 … … 7808 7815 \cos\vartheta_{t} & \sin\vartheta_{t}\\ 7809 7816 -\sin\vartheta_{t} & \cos\vartheta_{t} 7810 \end{array}\right] 7817 \end{array}\right], 7811 7818 \] 7812 7819 … … 7871 7878 e\vartheta_{0}\left(i_{\beta,0}\sin\vartheta_{0}+i_{\alpha,0}\cos\vartheta_{0}\right)\\ 7872 7879 0 7873 \end{array}\right) 7880 \end{array}\right), 7874 7881 \end{eqnarray*} 7875 7882 … … 7910 7917 0 & 0 & \Delta t & 1 & 0\\ 7911 7918 0 & 0 & 0 & 0 & 1 7912 \end{array}\right] 7919 \end{array}\right], 7913 7920 \] 7914 7921 … … 7918 7925 \begin_inset Formula 7919 7926 \[ 7920 \zeta=\left(i_{\beta,0}\sin\vartheta_{0}+i_{\alpha,0}\cos\vartheta_{0}\right) 7927 \zeta=\left(i_{\beta,0}\sin\vartheta_{0}+i_{\alpha,0}\cos\vartheta_{0}\right). 7921 7928 \] 7922 7929 … … 7948 7955 0 & 0\\ 7949 7956 0 & 0 7950 \end{array}\right] 7957 \end{array}\right], 7951 7958 \] 7952 7959 … … 7976 7983 . 7977 7984 V důsledku této substituce 7978 \begin_inset Formula $\psi_{t}=\omega_{t}-\overline{\omega}_{t}$ 7985 \begin_inset CommandInset ref 7986 LatexCommand ref 7987 reference "eq:substituce-psi" 7988 7979 7989 \end_inset 7980 7990 … … 7989 7999 \begin_inset Formula 7990 8000 \begin{eqnarray} 7991 i_{\alpha,t+1} & = & ai_{\alpha,t}+b\left(\psi_{t}+\overline{\omega}_{t}\right)\sin\vartheta_{t}+cu_{\alpha,t} \nonumber \\7992 i_{\beta,t+1} & = & ai_{\beta,t}-b\left(\psi_{t}+\overline{\omega}_{t}\right)\cos\vartheta_{t}+cu_{\beta,t} \label{eq:subst-system-s-psi}\\7993 \psi_{t+1} & = & d\psi_{t}+e\left(i_{b,t}\cos\vartheta_{t}-i_{\alpha,t}\sin\vartheta_{t}\right) \nonumber \\7994 \vartheta_{t+1} & = & \vartheta_{t}+\Delta t\left(\psi_{t}+\overline{\omega}_{t}\right) \nonumber8001 i_{\alpha,t+1} & = & ai_{\alpha,t}+b\left(\psi_{t}+\overline{\omega}_{t}\right)\sin\vartheta_{t}+cu_{\alpha,t},\nonumber \\ 8002 i_{\beta,t+1} & = & ai_{\beta,t}-b\left(\psi_{t}+\overline{\omega}_{t}\right)\cos\vartheta_{t}+cu_{\beta,t},\label{eq:subst-system-s-psi}\\ 8003 \psi_{t+1} & = & d\psi_{t}+e\left(i_{b,t}\cos\vartheta_{t}-i_{\alpha,t}\sin\vartheta_{t}\right),\nonumber \\ 8004 \vartheta_{t+1} & = & \vartheta_{t}+\Delta t\left(\psi_{t}+\overline{\omega}_{t}\right),\nonumber 7995 8005 \end{eqnarray} 7996 8006 … … 8022 8032 -e\sin\vartheta_{t} & e\cos\vartheta_{t} & d & -e\left(i_{\beta,t}\sin\vartheta_{t}+i_{\alpha,t}\cos\vartheta_{t}\right)\\ 8023 8033 0 & 0 & \Delta t & 1 8024 \end{array}\right] 8034 \end{array}\right], 8025 8035 \] 8026 8036 … … 8055 8065 \end_inset 8056 8066 8057 , která do něj vstupuje jako časově proměnný parametr .8067 , která do něj vstupuje jako časově proměnný parametr: 8058 8068 \begin_inset Formula 8059 8069 \begin{eqnarray*} … … 8063 8073 e\vartheta_{0}\left(i_{\beta,0}\sin\vartheta_{0}+i_{\alpha,0}\cos\vartheta_{0}\right)\\ 8064 8074 \Delta t\overline{\omega}_{t} 8065 \end{array}\right) 8075 \end{array}\right). 8066 8076 \end{eqnarray*} 8067 8077 … … 8084 8094 0 & 0 & \Delta t & 1 & \Delta t\overline{\omega}_{t}\\ 8085 8095 0 & 0 & 0 & 0 & 1 8086 \end{array}\right] 8096 \end{array}\right], 8087 8097 \] 8088 8098 … … 8092 8102 \begin_inset Formula 8093 8103 \[ 8094 \zeta=\left(i_{\beta,0}\sin\vartheta_{0}+i_{\alpha,0}\cos\vartheta_{0}\right) 8104 \zeta=\left(i_{\beta,0}\sin\vartheta_{0}+i_{\alpha,0}\cos\vartheta_{0}\right). 8095 8105 \] 8096 8106 … … 8237 8247 e\left(i_{b,t}\cos\vartheta_{t}-i_{\alpha,t}\sin\vartheta_{t}\right)\\ 8238 8248 0 8239 \end{array}\right) 8249 \end{array}\right). 8240 8250 \] 8241 8251 … … 8267 8277 -e\sin\vartheta_{t} & e\cos\vartheta_{t}\\ 8268 8278 0 & 0 8269 \end{array}\right] 8279 \end{array}\right]. 8270 8280 \] 8271 8281 … … 8311 8321 ai_{\alpha,t}+b\omega_{t}\sin\vartheta_{t}+cu_{\alpha,t}\\ 8312 8322 ai_{\beta,t}-b\omega_{t}\cos\vartheta_{t}+cu_{\beta,t} 8313 \end{array}\right) 8323 \end{array}\right), 8314 8324 \] 8315 8325 … … 8353 8363 \end_inset 8354 8364 8355 souřadnice. 8356 Vyjdeme ze soustavy rovnic ( 8365 souřadnice, vyjdeme však ze soustavy rovnic ( 8357 8366 \begin_inset CommandInset ref 8358 8367 LatexCommand ref … … 8361 8370 \end_inset 8362 8371 8363 ) 8364 \begin_inset Formula 8365 \begin{eqnarray*} 8366 i_{d,t+1} & \text{=} & ai_{d,t}+\Delta t\cdot i_{q,t}\omega_{t}+cu_{d,t}\\ 8367 i_{q,t+1} & \text{=} & ai_{q,t}-\Delta t\cdot i_{d,t}\omega_{t}-b\omega_{t}+cu_{q,t}\\ 8368 \omega_{t+1} & \text{=} & d\omega_{t}+ei_{q,t}\\ 8369 \vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t 8370 \end{eqnarray*} 8371 8372 \end_inset 8373 8374 8375 \end_layout 8376 8377 \begin_layout Standard 8378 Tento tvar rovnic je z hlediska linearizece daleko příznivější, protože 8372 ). 8373 Tento tvar rovnic je z hlediska linearizece daleko příznivější, protože 8379 8374 jedinými nelineárními členy jsou 8380 8375 \begin_inset Formula $\pm\Delta t\cdot i_{q,d}\omega$ … … 8453 8448 0 & 0 & \Delta t & 1 & \Delta t\overline{\omega}\\ 8454 8449 0 & 0 & 0 & 0 & 1 8455 \end{array}\right] \\8450 \end{array}\right],\\ 8456 8451 B & = & \left[\begin{array}{cc} 8457 8452 c & 0\\ … … 8460 8455 0 & 0\\ 8461 8456 0 & 0 8462 \end{array}\right] 8457 \end{array}\right]. 8463 8458 \end{eqnarray*} 8464 8459 … … 8472 8467 \end_layout 8473 8468 8469 \begin_layout Standard 8470 Jak již bylo uvedeno v části 8471 \begin_inset CommandInset ref 8472 LatexCommand ref 8473 reference "sub:Prehled-metod-dualniho-rizeni" 8474 8475 \end_inset 8476 8477 , hlavní myšlenka využití hyperstavu spočívá v aplikaci EKF v jistém smyslu 8478 dvakrát. 8479 To umožňuje získat kromě odhadu samotného stavu i odhad jeho kovarianční 8480 matice. 8481 Proč je právě znalost kovarianční matice pro konkrétní uvažovaný systém 8482 PMSM výhodná bude nejdříve ukázáno pomocí Bellmanovy funkce. 8483 Pak již bude následovat odvození samotného algoritmu založeného na hyperstavu. 8484 \end_layout 8485 8474 8486 \begin_layout Subsection 8475 Bellmanova funkce 8487 Bellmanova funkce pro PMSM 8476 8488 \begin_inset CommandInset label 8477 8489 LatexCommand label … … 8484 8496 8485 8497 \begin_layout Standard 8486 Jak bylo již uvedeno v odstavci8487 \begin_inset CommandInset ref8488 LatexCommand ref8489 reference "sub:Prehled-metod-dualniho-rizeni"8490 8491 \end_inset8492 8493 , využití hyperstavu umožňuje řídícímu algoritmu pracovat kromě odhadu stavu8494 i s odhadem jeho kovariance.8495 Proč se toto může jevit obzvláště výhodným ilustrují následující vztahy8496 týkající se Bellmanovi funkce.8497 Z tohoto důvodu je následující odstavec zařazen ještě před samotný popis8498 hyperstavu.8499 \end_layout8500 8501 \begin_layout Standard8502 8498 Cílem řídícího algoritmu je minimalizovat ztrátovou funkci uvažovanou ve 8503 8499 tvaru ( … … 8510 8506 ). 8511 8507 Klasickým postupem pro nalezení optimálního řešení této úlohy je užítí 8512 Bellmanovy funkce a algoritmu dynamického programování: 8508 Bellmanovy funkce a algoritmu dynamického programování jak bylo popsáno 8509 v odstavci 8510 \begin_inset CommandInset ref 8511 LatexCommand ref 8512 reference "sub:uloha-dualniho-rizeni" 8513 8514 \end_inset 8515 8516 . 8517 8513 8518 \end_layout 8514 8519 … … 8529 8534 \begin_inset Formula 8530 8535 \begin{equation} 8531 V_{t-1}\left(x_{t-1},u_{t-1}\right)=\min_{u_{t-1}\in U_{t-1}}\mathrm{E}\left\{ x_{t}^{T}Q_{t}x_{t}+u_{t}^{T}R_{t}u_{t}+V_{t}\left(x_{t},u_{t}\right)\mid\mathcal{I}_{t}\right\} \label{eq:bellVrek}8536 V_{t-1}\left(x_{t-1},u_{t-1}\right)=\min_{u_{t-1}\in U_{t-1}}\mathrm{E}\left\{ x_{t}^{T}Q_{t}x_{t}+u_{t}^{T}R_{t}u_{t}+V_{t}\left(x_{t},u_{t}\right)\mid\mathcal{I}_{t}\right\} ,\label{eq:bellVrek} 8532 8537 \end{equation} 8533 8538 … … 8546 8551 \end_inset 8547 8552 8548 , kde střední hodnota je podmíněna 8553 , kde střední hodnota je podmíněna informačním vektorem 8549 8554 \begin_inset Formula $\mathcal{I}_{t}$ 8550 8555 \end_inset 8551 8556 8552 , kter éreprezentuje současně dostupnou informaci o systému zahrnující všechna8557 , který reprezentuje současně dostupnou informaci o systému zahrnující všechna 8553 8558 měření a řídící vstupy do času 8554 8559 \begin_inset Formula $t$ … … 8562 8567 \begin_inset Formula 8563 8568 \[ 8564 x_{t}^{T}Q_{t}x_{t}+u_{t}^{T}R_{t}u_{t} 8569 x_{t}^{T}Q_{t}x_{t}+u_{t}^{T}R_{t}u_{t}, 8565 8570 \] 8566 8571 … … 8578 8583 \end_inset 8579 8584 8580 ) přejde na 8585 ) a zavedení substituce 8586 \begin_inset CommandInset ref 8587 LatexCommand ref 8588 reference "eq:substituce-psi" 8589 8590 \end_inset 8591 8592 přejde na 8581 8593 \begin_inset Formula 8582 8594 \[ 8583 q\left(x_{t}^{(3)}-\overline{\omega}_{t}\right) +u_{t}^{T}R_{t}u_{t}8595 q\left(x_{t}^{(3)}-\overline{\omega}_{t}\right)^{2}+u_{t}^{T}R_{t}u_{t}, 8584 8596 \] 8585 8597 … … 8598 8610 \begin{alignat}{1} 8599 8611 V_{t-1} & \left(x_{t-1},u_{t-1}\right)\text{=}\min_{u_{t-1}}\mathrm{E}\left\{ x_{t}^{T}Q_{t}x_{t}+u_{t}^{T}R_{t}u_{t}+V_{t}\left(x_{t},u_{t}\right)\mid\mathcal{I}_{t}\right\} \nonumber \\ 8600 = & \min_{u_{t-1}}\left(\mathrm{E}\left\{ q\left(x_{t}^{(3)}-\overline{\omega}_{t}\right)\right\} +\mathrm{E}\left\{ u_{t}^{T}R_{t}u_{t}+V_{t}\left(x_{t},u_{t}\right)\mid\mathcal{I}_{t}\right\} \right)\nonumber \\ 8601 \text{=} & \min_{u_{t-1}}\left(q\mathrm{E}\left\{ \left(x_{t}^{(3)}\right)^{2}+\overline{\omega}_{t}^{2}+2x_{t}^{(3)}\overline{\omega}_{t}\right\} +\mathrm{E}\left\{ u_{t}^{T}R_{t}u_{t}+V_{t}\left(x_{t},u_{t}\right)\mid\mathcal{I}_{t}\right\} \right)\label{eq:bellamn-s-Pome}\\ 8612 = & \min_{u_{t-1}}\left(\mathrm{E}\left\{ q\left(x_{t}^{(3)}-\overline{\omega}_{t}\right)^{2}\right\} +\mathrm{E}\left\{ u_{t}^{T}R_{t}u_{t}+V_{t}\left(x_{t},u_{t}\right)\mid\mathcal{I}_{t}\right\} \right)\label{eq:bellamn-s-Pome}\\ 8602 8613 = & \min_{u_{t-1}}\left(q\left(\mathrm{E}\left\{ \left(x_{t}^{(3)}\right)^{2}\right\} +\mathrm{E}\left\{ \overline{\omega}_{t}^{2}\right\} +\mathrm{E}\left\{ 2x_{t}^{(3)}\overline{\omega}_{t}\right\} \right)+\mathrm{E}\left\{ u_{t}^{T}R_{t}u_{t}+V_{t}\left(x_{t},u_{t}\right)\mid\mathcal{I}_{t}\right\} \right)\nonumber \\ 8603 8614 = & \min_{u_{t-1}}\left(q\left(\left(\hat{x}_{t}^{(3)}\right)^{2}+\mathrm{Var}\left(x_{t}^{(3)}\right)+\overline{\omega}_{t}^{2}+2\hat{x}_{t}^{(3)}\overline{\omega}_{t}\right)+\mathrm{E}\left\{ u_{t}^{T}R_{t}u_{t}+V_{t}\left(x_{t},u_{t}\right)\mid\mathcal{I}_{t}\right\} \right)\nonumber \\ 8604 = & \min_{u_{t-1}}\left(q\left(\hat{x}_{t}^{(3)}-\overline{\omega}_{t}\right)+q\mathrm{Var}\left(x_{t}^{(3)}\right)+\mathrm{E}\left\{ u_{t}^{T}R_{t}u_{t}+V_{t}\left(x_{t},u_{t}\right)\mid\mathcal{I}_{t}\right\} \right) \nonumber8615 = & \min_{u_{t-1}}\left(q\left(\hat{x}_{t}^{(3)}-\overline{\omega}_{t}\right)+q\mathrm{Var}\left(x_{t}^{(3)}\right)+\mathrm{E}\left\{ u_{t}^{T}R_{t}u_{t}+V_{t}\left(x_{t},u_{t}\right)\mid\mathcal{I}_{t}\right\} \right),\nonumber 8605 8616 \end{alignat} 8606 8617 … … 8615 8626 \end_inset 8616 8627 8617 a dále jsme využili toho, že 8628 a dále jsme využili vztahu 8629 \begin_inset Formula $\mathrm{Var}\left(x\right)=\mathrm{E}\left\{ x^{2}\right\} -\left(\mathrm{E}\left\{ x\right\} \right)^{2}$ 8630 \end_inset 8631 8632 a toho, že 8618 8633 \begin_inset Formula $\overline{\omega}_{t}$ 8619 8634 \end_inset 8620 8635 8621 je daný parametr a tedy je pro výpočet střední hodnoty konstantou a vztahu 8622 8623 \begin_inset Formula $\mathrm{Var}\left(x\right)=\mathrm{E}\left\{ x^{2}\right\} -\left(\mathrm{E}\left\{ x\right\} \right)^{2}$ 8624 \end_inset 8625 8626 . 8627 Tedy ve výpočtu Bellmanovy funkce 8636 je známý parametr a tedy je pro výpočet střední hodnoty konstantou. 8637 Následně můžeme ve výpočtu Bellmanovy funkce 8628 8638 \begin_inset Formula $V$ 8629 8639 \end_inset … … 8636 8646 \end_inset 8637 8647 8638 ) můžemenáhodnou veličinu8648 ) náhodnou veličinu 8639 8649 \begin_inset Formula $x_{t}$ 8640 8650 \end_inset … … 8648 8658 \end_inset 8649 8659 8650 , tj. 8651 varianci otáček stroje. 8660 , to jest varianci otáček stroje. 8652 8661 \end_layout 8653 8662 … … 8657 8666 8658 8667 \begin_layout Standard 8659 Následující postup s hyperstavem vychází z 8660 \begin_inset CommandInset citation 8661 LatexCommand cite 8662 key "Kim2006" 8663 8664 \end_inset 8665 8666 . 8667 V tomto článku je však narozdíl od následujícího postupu používán spojitý 8668 čas. 8669 \end_layout 8670 8671 \begin_layout Standard 8672 Jedná se o analogii s LQG v předchozí části, s tím rozdílem, že použijem 8673 EKF algoritmus v jistém smyslu jakoby dvakrát. 8668 Následující algoritmus využívající hyperstav je praktickou realizací postupu 8669 nastíněného v ostavci 8670 \begin_inset CommandInset ref 8671 LatexCommand ref 8672 reference "sub:Prehled-metod-dualniho-rizeni" 8673 8674 \end_inset 8675 8676 . 8677 Jedná se o analogii s LQG popsaným v předchozích částech, s tím rozdílem, 8678 že místo stavu je aplikován právě na hyperstav. 8674 8679 Protože tímto přístupem již značně narůstá dimenzionalita úlohy je z výpočetníc 8675 8680 h důvodů výhodnější užití redukovaného modelu, i přes komplikace, které … … 8691 8696 \end_inset 8692 8697 8693 i odhad jeho variance v podobě matice8698 i odhad jeho kovariance v podobě matice 8694 8699 \begin_inset Formula 8695 8700 \[ … … 8716 8721 \end_inset 8717 8722 8718 ) představují předpis pro výpočet 8719 \begin_inset Formula $P$ 8720 \end_inset 8721 8722 : 8723 \end_layout 8724 8725 \begin_layout Standard 8726 \begin_inset Formula 8727 \begin{eqnarray} 8728 \overline{P} & = & APA^{T}+V\nonumber \\ 8729 S & = & C\overline{P}C^{T}+W\nonumber \\ 8730 K & = & \overline{P}C^{T}S^{-1}\label{eq:ekf-pro-hyperstav}\\ 8731 P^{+} & \text{=} & \left(I-KC\right)\overline{P}\nonumber 8732 \end{eqnarray} 8733 8734 \end_inset 8735 8736 kde jsou z důvodu jednoduššího zápisu vynechány časové indexy 8723 ) představují předpis její výpočet. 8724 Z důvodu jednoduššího zápisu budou vynechány časové indexy 8737 8725 \begin_inset Formula $t$ 8738 8726 \end_inset … … 8747 8735 8748 8736 . 8737 Použité rovnice pro EKF jsou tedy ve tvaru 8738 \end_layout 8739 8740 \begin_layout Standard 8741 \begin_inset Formula 8742 \begin{eqnarray} 8743 \overline{P} & = & APA^{T}+V,\nonumber \\ 8744 S & = & C\overline{P}C^{T}+W,\nonumber \\ 8745 K & = & \overline{P}C^{T}S^{-1},\label{eq:ekf-pro-hyperstav}\\ 8746 P^{+} & \text{=} & \left(I-KC\right)\overline{P}.\nonumber 8747 \end{eqnarray} 8748 8749 \end_inset 8750 8751 8749 8752 \end_layout 8750 8753 … … 8765 8768 \begin_inset Formula 8766 8769 \[ 8767 \xi_{t}=\left(\omega_{t},\vartheta_{t},P_{\omega},P_{\omega\vartheta},P_{\vartheta}\right)^{T} 8770 \xi_{t}=\left(\omega_{t},\vartheta_{t},P_{\omega},P_{\omega\vartheta},P_{\vartheta}\right)^{T}. 8768 8771 \] 8769 8772 8770 8773 \end_inset 8771 8774 8772 Na hyperstav již aplikujeme algoritmus pro LQG, jak byl popsán v předchozí 8773 části. 8774 Problém však představuje nalezení matice derivací 8775 \begin_inset Formula $A$ 8776 \end_inset 8777 8778 , protože je třeba derivovat maticové rovnice pro výpočet EKF ( 8775 Na hyperstav následně aplikujeme algoritmus LQG viz část 8776 \begin_inset CommandInset ref 8777 LatexCommand ref 8778 reference "sec:LQG-obecne" 8779 8780 \end_inset 8781 8782 , podobně jako pro EKF a LQ regulátor pro redukovaný model jak bylo popsáno 8783 v předchozích částech 8784 \begin_inset CommandInset ref 8785 LatexCommand ref 8786 reference "sub:EKF-Redukovany-model" 8787 8788 \end_inset 8789 8790 a 8791 \begin_inset CommandInset ref 8792 LatexCommand ref 8793 reference "sub:LQ-řízení-pro-red-model" 8794 8795 \end_inset 8796 8797 . 8798 Větší problém však představuje nalezení matice derivací 8799 \begin_inset Formula $A_{t}$ 8800 \end_inset 8801 8802 , protože je třeba navíc derivovat maticové rovnice pro výpočet EKF ( 8779 8803 \begin_inset CommandInset ref 8780 8804 LatexCommand ref … … 8796 8820 \end_inset 8797 8821 8798 ) dle jednotlivých složek vektoru hyperstavu 8799 \begin_inset Formula $\xi$ 8822 ) dle jednotlivých složek 8823 \begin_inset Formula $\xi_{i}$ 8824 \end_inset 8825 8826 vektoru hyperstavu 8827 \begin_inset Formula $\xi_{t}$ 8800 8828 \end_inset 8801 8829 … … 8803 8831 \begin_inset Formula 8804 8832 \begin{eqnarray} 8805 \frac{\partial\overline{P}}{\partial\xi_{i}} & = & \frac{\partial A}{\partial\xi_{i}}PA^{T}+A\frac{\partial P}{\partial\xi_{i}}A^{T}+AP\frac{\partial A^{T}}{\partial\xi_{i}} \nonumber \\8806 \frac{\partial S}{\partial\xi_{i}} & = & \frac{\partial C}{\partial\xi_{i}}\overline{P}C^{T}+C\frac{\partial\overline{P}}{\partial\xi_{i}}C^{T}+C\overline{P}\frac{\partial C^{T}}{\partial\xi_{i}} \nonumber \\8807 \frac{\partial K}{\partial\xi_{i}} & = & \frac{\partial\overline{P}}{\partial\xi_{i}}C^{T}S^{-1}+\overline{P}\frac{\partial C^{T}}{\partial\xi_{i}}S^{-1}-\overline{P}C^{T}S^{-1}\frac{\partial S}{\partial\xi_{i}}S^{-1} \label{eq:ekf-stav-derivace}\\8808 \frac{\partial P^{+}}{\partial\xi_{i}} & \text{=} & \frac{\partial\overline{P}}{\partial\xi_{i}}-\frac{\partial K}{\partial\xi_{i}}C\overline{P}-K\frac{\partial C}{\partial\xi_{i}}\overline{P}-KC\frac{\partial\overline{P}}{\partial\xi_{i}} \nonumber8833 \frac{\partial\overline{P}}{\partial\xi_{i}} & = & \frac{\partial A}{\partial\xi_{i}}PA^{T}+A\frac{\partial P}{\partial\xi_{i}}A^{T}+AP\frac{\partial A^{T}}{\partial\xi_{i}},\nonumber \\ 8834 \frac{\partial S}{\partial\xi_{i}} & = & \frac{\partial C}{\partial\xi_{i}}\overline{P}C^{T}+C\frac{\partial\overline{P}}{\partial\xi_{i}}C^{T}+C\overline{P}\frac{\partial C^{T}}{\partial\xi_{i}},\nonumber \\ 8835 \frac{\partial K}{\partial\xi_{i}} & = & \frac{\partial\overline{P}}{\partial\xi_{i}}C^{T}S^{-1}+\overline{P}\frac{\partial C^{T}}{\partial\xi_{i}}S^{-1}-\overline{P}C^{T}S^{-1}\frac{\partial S}{\partial\xi_{i}}S^{-1},\label{eq:ekf-stav-derivace}\\ 8836 \frac{\partial P^{+}}{\partial\xi_{i}} & \text{=} & \frac{\partial\overline{P}}{\partial\xi_{i}}-\frac{\partial K}{\partial\xi_{i}}C\overline{P}-K\frac{\partial C}{\partial\xi_{i}}\overline{P}-KC\frac{\partial\overline{P}}{\partial\xi_{i}},\nonumber 8809 8837 \end{eqnarray} 8810 8838 … … 8846 8874 A_{1} & A_{2} & 0 & 0 & 0\\ 8847 8875 \left(\frac{\partial P^{+}}{\partial\omega}\right)^{sl} & \left(\frac{\partial P^{+}}{\partial\vartheta}\right)^{sl} & \left(\frac{\partial P^{+}}{\partial P_{\omega}}\right)^{sl} & \left(\frac{\partial P^{+}}{\partial P_{\omega\vartheta}}\right)^{sl} & \left(\frac{\partial P^{+}}{\partial P_{\vartheta}}\right)^{sl} 8848 \end{array}\right] 8876 \end{array}\right], 8849 8877 \] 8850 8878 … … 8859 8887 \end_inset 8860 8888 8861 -tý sloupec matice 8862 \begin_inset Formula $A$ 8863 \end_inset 8864 8865 , zápis 8889 -tý sloupec původní matice 8890 \begin_inset Formula $A_{t}$ 8891 \end_inset 8892 8893 viz ( 8894 \begin_inset CommandInset ref 8895 LatexCommand ref 8896 reference "eq:matice-ekf-redukovanymodel-ls" 8897 8898 \end_inset 8899 8900 ), 8866 8901 \begin_inset Formula $0$ 8867 8902 \end_inset 8868 8903 8869 je sloupec nul vhodné délky a parciální derivace8904 označuje sloupec nul vhodné délky a parciální derivace 8870 8905 \begin_inset Formula $P^{+}$ 8871 8906 \end_inset … … 8879 8914 sl 8880 8915 \emph default 8881 8916 , tedy 8882 8917 \begin_inset Formula $\left(\frac{\partial P^{+}}{\partial\xi_{i}}\right)^{sl}$ 8883 8918 \end_inset 8884 8919 8885 je myšlena v tom smyslu, že po vypočtení příslušné derivace8920 , je myšlena v tom smyslu, že po vypočtení příslušné derivace 8886 8921 \begin_inset Formula $\frac{\partial P^{+}}{\partial\xi_{i}}$ 8887 8922 \end_inset … … 8895 8930 8896 8931 ) jsou z této matice vybrány 3 z jejích 4 prvků tvořící horní nebo dolní 8897 trojúhelník a zapísány ve smyslu tvorby vektoru hyperstavu do sloupce: 8932 trojúhelník a zapísány ve smyslu tvorby vektoru hyperstavu do sloupce. 8933 Tedy matice 8898 8934 \begin_inset Formula 8899 8935 \[ … … 8906 8942 \end_inset 8907 8943 8908 8944 je zapsána jako 8909 8945 \begin_inset Formula 8910 8946 \[ 8911 8947 \left(\frac{\partial P^{+}}{\partial\xi_{i}}\right)^{sl}=\left(\begin{array}{ccc} 8912 \frac{\partial P_{\omega}^{+}}{\partial\xi_{i}} & \frac{\partial P_{\omega\vartheta}^{+}}{\partial\xi_{i}} & \frac{\partial P_{\vartheta}^{+}}{\partial\xi_{i}}\end{array}\right)^{T} 8948 \frac{\partial P_{\omega}^{+}}{\partial\xi_{i}} & \frac{\partial P_{\omega\vartheta}^{+}}{\partial\xi_{i}} & \frac{\partial P_{\vartheta}^{+}}{\partial\xi_{i}}\end{array}\right)^{T}. 8913 8949 \] 8914 8950 … … 8933 8969 8934 8970 pouze doplněnou nulami na vhodný rozměr. 8935 Pro lineárně kvadratick é řízení platí opět totéž, co pro jednoduché (tj.8936 bez hyperstavu) a matici8971 Pro lineárně kvadratický regulátor platí opět totéž, co pro jednoduchý 8972 (tedy bez hyperstavu) a matici 8937 8973 \begin_inset Formula $A_{hyp}$ 8938 8974 \end_inset 8939 8975 8940 je tedy třeba rozšířit zahrnutím konstantních členů, dále je třeba ošetřit 8941 substitucí řízení na nenulové požadované otáčky 8976 je třeba rozšířit zahrnutím konstantních členů, dále je třeba ošetřit substituc 8977 í řízení ( 8978 \begin_inset CommandInset ref 8979 LatexCommand ref 8980 reference "eq:substituce-psi" 8981 8982 \end_inset 8983 8984 ) na nenulové požadované otáčky 8942 8985 \begin_inset Formula $\overline{\omega}$ 8943 8986 \end_inset 8944 8987 8945 8988 . 8989 Kvůli větší složitosti tohoto výpočtu bylo rozšíření matice 8990 \begin_inset Formula $A_{hyp}$ 8991 \end_inset 8992 8993 prováděno pomocí symbolický výpočtů v programu Matlab, výsledek je pak 8994 uveden v příloze ( 8995 \series bold 8996 zatím teda není 8997 \series default 8998 ). 8999 Protože uvažujeme redukovaný model je třeba dále užít zřetězení dvou LQ 9000 regulátorů, podobně jako v případě bez hyperstavu v odstavci 9001 \begin_inset CommandInset ref 9002 LatexCommand ref 9003 reference "sub:LQ-řízení-pro-red-model" 9004 9005 \end_inset 9006 9007 . 8946 9008 8947 9009 \end_layout 8948 9010 8949 9011 \begin_layout Standard 8950 Protože uvažujeme redukovaný model je třeba užít zřetězení dvou LQ regulátorů. 8951 Výhodou využití hyperstavu ale je, že máme k dispozici i odhady variancí 8952 9012 Výhodou využití hyperstavu je, že máme k dispozici i odhady kovarianční 9013 matice 8953 9014 \begin_inset Formula $P$ 8954 9015 \end_inset … … 8980 9041 8981 9042 \begin_layout Standard 8982 Tedy pro stav9043 Pro stav 8983 9044 \begin_inset Formula 8984 9045 \[ … … 8988 9049 \end_inset 8989 9050 8990 vypočteme z rovnic pro EKF kovarianční matici 9051 vypočteme z rovnic ( 9052 \begin_inset CommandInset ref 9053 LatexCommand ref 9054 reference "eq:ekf-pro-hyperstav" 9055 9056 \end_inset 9057 9058 ) EKF kovarianční matici 8991 9059 \begin_inset Formula 8992 9060 \[ … … 9016 9084 \begin_inset Formula 9017 9085 \[ 9018 \xi_{t}=\left(i_{\alpha,t},i_{\beta,t},\omega_{t},\vartheta_{t},P_{5},P_{6},P_{7},P_{8},P_{9},P_{10},P_{11},P_{12},P_{13},P_{14}\right)^{T} 9086 \xi_{t}=\left(i_{\alpha,t},i_{\beta,t},\omega_{t},\vartheta_{t},P_{5},P_{6},P_{7},P_{8},P_{9},P_{10},P_{11},P_{12},P_{13},P_{14}\right)^{T}. 9019 9087 \] 9020 9088 … … 9029 9097 \end_inset 9030 9098 9031 , jsou formálně shodné s rovnicemi pro redukovaný model, pouze rozměry vystupují 9032 cích matic jsou větší. 9033 A matice 9099 , jsou formálně shodné s rovnicemi ( 9100 \begin_inset CommandInset ref 9101 LatexCommand ref 9102 reference "eq:ekf-stav-derivace" 9103 9104 \end_inset 9105 9106 ) pro redukovaný model, pouze rozměry vystupujících matic jsou větší. 9107 Matice 9034 9108 \begin_inset Formula $A_{hyp}$ 9035 9109 \end_inset 9036 9110 9037 je vetvaru9111 je pak v blokovém tvaru 9038 9112 \begin_inset Formula 9039 9113 \[ … … 9042 9116 \\ 9043 9117 \left(\frac{\partial P^{+}}{\partial\xi_{i}}\right)_{i\in\left\{ 1\ldots14\right\} }^{sl} 9044 \end{array}\right] 9118 \end{array}\right], 9045 9119 \] 9046 9120 9047 9121 \end_inset 9048 9122 9049 9123 kde 9124 \begin_inset Formula $A$ 9125 \end_inset 9126 9127 značí původní matici 9128 \begin_inset Formula $A_{t}$ 9129 \end_inset 9130 9131 ( 9132 \begin_inset CommandInset ref 9133 LatexCommand ref 9134 reference "eq:matice-ekf-plnymodel-ls" 9135 9136 \end_inset 9137 9138 ) pro EKF standartního plného stavu, 9139 \begin_inset Formula $0$ 9140 \end_inset 9141 9142 je nulová matice vhodného rozměru a zápis 9143 \begin_inset Formula $\left(\frac{\partial P^{+}}{\partial\xi_{i}}\right)^{sl}$ 9144 \end_inset 9145 9146 má stejný význam jako v předchozím odstavci, jen vybírá trojúhelník z větší 9147 matice; tyto sloupcové vektory jsou pak po řadě, dle indexu 9148 \begin_inset Formula $i$ 9149 \end_inset 9150 9151 , zapsány do blokové matice 9152 \begin_inset Formula $A_{hyp}$ 9153 \end_inset 9154 9155 . 9156 Pro LQ regulátor je matici opět třeba rozšířit o konstantní člen, který 9157 byl počítán symbolicky a je možné jej nalézt v příloze ( 9158 \series bold 9159 zatím teda ne 9160 \series default 9161 ). 9162 Výhodou je opět dostupnost odhadu kovarianční matice 9163 \begin_inset Formula $P$ 9164 \end_inset 9165 9166 a oproti redukovanému modelu se zde navíc nevyskytují komplikace při návrhu 9167 řízení jako pro redukovaný model. 9050 9168 \end_layout 9051 9169 9052 9170 \begin_layout Section 9053 Bikriteriální metoda( 9054 \series medium 9055 nutno přepracovat a doplnit další verze 9056 \series default 9057 ) 9171 Bikriteriální metoda 9058 9172 \end_layout 9059 9173 … … 9063 9177 9064 9178 \begin_layout Standard 9065 Posledním z implementovaných algoritmů je následující 9066 \emph on 9067 jednoduchý duální návrh 9068 \emph default 9069 . 9070 Hlavní myšlenka je založena na 9179 Dalším z implementovaných algoritmů je jednoduchý návrh založený na 9071 9180 \emph on 9072 9181 bikriteriální metodě … … 9075 9184 \begin_inset CommandInset ref 9076 9185 LatexCommand ref 9077 reference "sub: Vybrané-algoritmy-proDC"9186 reference "sub:Prehled-metod-dualniho-rizeni" 9078 9187 9079 9188 \end_inset … … 9092 9201 , které se pod tímto pojmem obvykle rozumí není v případě zde uvažovaného 9093 9202 systému snadné. 9094 Proto místo něj využijeme opět LQ řízení v 9095 \begin_inset Formula $d-q$ 9096 \end_inset 9097 9098 souřadnicích. 9203 Proto místo něj využijeme některé standartní řízení, například vektorové 9204 založené na PI nebo LQ regulátorech. 9099 9205 Toto není z hlediska bikriteriální metody korektní, zde uvažovaný postup 9100 je ale myšlen jako jednoduchý duálnínávrh a je pouze jejím jistým přiblížením.9206 je ale myšlen jako jednoduchý návrh a je pouze jejím jistým přiblížením. 9101 9207 \end_layout 9102 9208 … … 9113 9219 9114 9220 ve tvaru 9115 \begin_inset Formula $\left\langle \tilde{u}_{d}-\varepsilon ,\tilde{u}_{d}+\varepsilon\right\rangle \times\left\langle \tilde{u}_{q}-\varepsilon,\tilde{u}_{q}+\varepsilon\right\rangle $9221 \begin_inset Formula $\left\langle \tilde{u}_{d}-\varepsilon_{d},\tilde{u}_{d}+\varepsilon_{d}\right\rangle \times\left\langle \tilde{u}_{q}-\varepsilon_{q},\tilde{u}_{q}+\varepsilon_{q}\right\rangle $ 9116 9222 \end_inset 9117 9223 … … 9122 9228 9123 9229 \begin_layout Standard 9124 Jak již bylo uvedeno v kapitole 9125 \begin_inset CommandInset ref 9126 LatexCommand ref 9127 reference "sec:Estimace-stavových-veličin" 9128 9129 \end_inset 9130 9131 , čím jsou vyšší otáčky, tím získáváme lepší odhad stavových veličin 9132 \begin_inset Formula $\omega$ 9133 \end_inset 9134 9135 a 9136 \begin_inset Formula $\vartheta$ 9137 \end_inset 9138 9139 , protože na otáčkách přímo úměrně závisí velikost zpětné elektromotorické 9140 síly. 9141 Na tomto základě můžeme uvažovat, že optimální buzení pro PMSM je takové, 9142 které se snaží maximalizovat otáčky 9230 Velikost zpětné elektromotorické síly, na základě které jsou odhadovány 9231 mechanické veličiny, je přímo úměrná otáčkám a tedy čím jsou vyšší otáčky, 9232 tím získáváme lepší odhad těchto veličin. 9233 Můžeme tedy uvažovat, že optimální buzení pro PMSM je takové, které se 9234 snaží maximalizovat otáčky 9143 9235 \begin_inset Formula $\omega$ 9144 9236 \end_inset … … 9152 9244 \begin_inset CommandInset ref 9153 9245 LatexCommand ref 9154 reference "eq:rovnice_jedn_dq" 9155 9156 \end_inset 9157 9158 ) 9159 \begin_inset Formula 9160 \begin{eqnarray*} 9161 i_{d,t+1} & = & ai_{d,t}+\Delta t\cdot i_{q,t}\omega_{t}+cu_{d,t},\\ 9162 i_{q,t+1} & = & ai_{q,t}-\Delta t\cdot i_{d,t}\omega_{t}-b\omega_{t}+cu_{q,t},\\ 9163 \omega_{t+1} & \text{=} & d\omega_{t}+ei_{q,t},\\ 9164 \vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t, 9165 \end{eqnarray*} 9166 9167 \end_inset 9168 9169 kde do třetí rovnice dosadíme z prvních dvou 9246 reference "eq:diskretni-system-dq-ls" 9247 9248 \end_inset 9249 9250 ), kde do třetí rovnice dosadíme z prvních dvou 9170 9251 \begin_inset Formula 9171 9252 \begin{equation} … … 9175 9256 \end_inset 9176 9257 9177 Dosazovat by šlo samozřejmě dále, ale již teď je vidět, jak je vhodné volit9178 9258 Dosazovat by šlo samozřejmě dále, ale již teď je zřejmé, jak volit řídící 9259 zásahy 9179 9260 \begin_inset Formula $u_{d}$ 9180 9261 \end_inset … … 9184 9265 \end_inset 9185 9266 9186 . 9187 Chceme maximalizovat 9267 : Chceme maximalizovat 9188 9268 \begin_inset Formula $\left|\omega\right|$ 9189 9269 \end_inset … … 9220 9300 Výsledné řízení je tedy 9221 9301 \begin_inset Formula 9302 \begin{eqnarray} 9303 u_{d} & = & \tilde{u}_{d}-\varepsilon_{d}\,\mathrm{sign}\,\omega,\nonumber \\ 9304 u_{q} & = & \tilde{u}_{q}+\varepsilon_{q}\,\mathrm{sign}\,\omega.\label{eq:zakladni-bikriterialní} 9305 \end{eqnarray} 9306 9307 \end_inset 9308 9309 9310 \end_layout 9311 9312 \begin_layout Subsection 9313 Časový posun hodnot 9314 \end_layout 9315 9316 \begin_layout Standard 9317 Předchozí návrh byl poměrně přímočarým vyvozením rovnic ( 9318 \begin_inset CommandInset ref 9319 LatexCommand ref 9320 reference "eq:zakladni-bikriterialní" 9321 9322 \end_inset 9323 9324 ) pro řízení z rovnice pro otáčky ( 9325 \begin_inset CommandInset ref 9326 LatexCommand ref 9327 reference "eq:rovnice_pro_opt_buz" 9328 9329 \end_inset 9330 9331 ) a nebyly při něm respektovány správné časy jednotlivých veličin. 9332 Chceme-li respektovat správný časový posun uvažovaných veličin, je třeba 9333 rovnice ( 9334 \begin_inset CommandInset ref 9335 LatexCommand ref 9336 reference "eq:zakladni-bikriterialní" 9337 9338 \end_inset 9339 9340 ) přepsat vzhledem k rovnice ( 9341 \begin_inset CommandInset ref 9342 LatexCommand ref 9343 reference "eq:rovnice_pro_opt_buz" 9344 9345 \end_inset 9346 9347 ) do tvaru 9348 \begin_inset Formula 9222 9349 \begin{eqnarray*} 9223 u_{d } & = & \tilde{u}_{d}-\varepsilon\,\mathrm{sign}\,\omega,\\9224 u_{q } & = & \tilde{u}_{q}+\varepsilon\,\mathrm{sign}\,\omega.9350 u_{d,t} & = & \tilde{u}_{d,t}-\varepsilon_{d}\,\mathrm{sign}\,\omega_{t+3},\\ 9351 u_{q,t} & = & \tilde{u}_{q,t}+\varepsilon_{q}\,\mathrm{sign}\,\omega_{t+2}. 9225 9352 \end{eqnarray*} 9226 9353 9227 9354 \end_inset 9228 9355 9229 9230 \end_layout 9231 9232 \begin_layout Standard 9233 Tento postup je relativně jednoduchou modifikací předchozího LQ algoritmu, 9234 ale jak ukazují simulace, může přinést značnou výhodu při určování počátečního 9235 natočení rotoru. 9356 Zde vystupující veličiny 9357 \begin_inset Formula $\omega_{t+2}$ 9358 \end_inset 9359 9360 a 9361 \begin_inset Formula $\omega_{t+3}$ 9362 \end_inset 9363 9364 lze získat, podobně jako rovnici ( 9365 \begin_inset CommandInset ref 9366 LatexCommand ref 9367 reference "eq:rovnice_pro_opt_buz" 9368 9369 \end_inset 9370 9371 ), dosazováním do třetí rovnice z ostatních v soustavě ( 9372 \begin_inset CommandInset ref 9373 LatexCommand ref 9374 reference "eq:diskretni-system-dq-ls" 9375 9376 \end_inset 9377 9378 ). 9379 Jejich tvar je pak následující 9380 \begin_inset Formula 9381 \begin{eqnarray*} 9382 \omega_{t+2} & \text{=} & \left(d^{2}-be\right)\omega_{t}+\left(a+d\right)ei_{q,t}-e\Delta t\cdot i_{d,t}\omega_{t}+ceu_{q,t},\\ 9383 \omega_{t+3} & \text{=} & d\omega_{t+2}+ceu_{q,t+1}-be\left(a+d\right)\omega_{t}+e\left(a^{2}-be\right)i_{q,t}+aceu_{q,t}-\\ 9384 & - & e\Delta t\left(a+ad\right)i_{d,t}\omega_{t}-ae^{2}\Delta t\cdot i_{d,t}i_{q,t}-ed\left(\Delta t\right)^{2}i_{q,t}\omega_{t}^{2}-e^{2}\left(\Delta t\right)^{2}i_{q,t}^{2}\omega_{t}-\\ 9385 & - & cde\Delta t\cdot\omega_{t}u_{d,t}-ce^{2}\Delta t\cdot i_{q,t}u_{d,t}. 9386 \end{eqnarray*} 9387 9388 \end_inset 9389 9390 Výpočet řídících zásahů pak probíhá v následujícím pořadí: 9391 \end_layout 9392 9393 \begin_layout Standard 9394 Nejdříve je vypočtena pomocná veličina 9395 \begin_inset Formula 9396 \[ 9397 \phi=\left(d^{2}-be\right)\omega_{t}+\left(a+d\right)ei_{q,t}-e\Delta ti_{d,t}\omega_{t} 9398 \] 9399 9400 \end_inset 9401 9402 a následně řídící zásah v ose 9403 \begin_inset Formula $q$ 9404 \end_inset 9405 9406 9407 \begin_inset Formula 9408 \[ 9409 u_{q,t}=\tilde{u}_{q,t}+\varepsilon_{q}\,\mathrm{sign}\,\phi. 9410 \] 9411 9412 \end_inset 9413 9414 To umožní výpočet druhé pomocné veličiny 9415 \begin_inset Formula 9416 \begin{eqnarray*} 9417 \zeta & = & d\phi+aceu_{q,t}-e\left(ab+bd+a\Delta t\left(1+d\right)i_{d,t}\right)\omega_{t}+\\ 9418 & + & e\left(a^{2}-be-ae\Delta t\cdot i_{d,t}\right)i_{q,t}-e\left(\Delta t\right)^{2}\left(d\omega_{t}+ei_{q,t}\right)i_{q,t}\omega_{t}. 9419 \end{eqnarray*} 9420 9421 \end_inset 9422 9423 Na závěr je vypočten řídící zásah v ose 9424 \begin_inset Formula $d$ 9425 \end_inset 9426 9427 9428 \begin_inset Formula 9429 \[ 9430 u_{d,t}=\tilde{u}_{d,t}-\varepsilon_{d}\,\mathrm{sign}\,\zeta. 9431 \] 9432 9433 \end_inset 9434 9435 9436 \end_layout 9437 9438 \begin_layout Standard 9439 Ještě je třeba zmínit, že bylo zanedbáno budoucí řízení 9440 \begin_inset Formula $u_{q,t+1}$ 9441 \end_inset 9442 9443 , které není známe a předchozím výpočtu je předpokládáno, že 9444 \begin_inset Formula $u_{q,t+1}=0$ 9445 \end_inset 9446 9447 . 9448 Lepších výsledků by mohlo být dosaženo užitím předpokladu 9449 \begin_inset Formula $u_{q,t+1}=u_{q,t}$ 9450 \end_inset 9451 9452 , pak by se rovnice pro výpočet výsledného řízení změnily na 9453 \begin_inset Formula $u_{d,t}=\tilde{u}_{d,t}-\varepsilon_{d}\,\mathrm{sign}\,\zeta'$ 9454 \end_inset 9455 9456 , kde 9457 \begin_inset Formula $\zeta'=\zeta+ceu_{q,t+1}$ 9458 \end_inset 9459 9460 . 9236 9461 \end_layout 9237 9462 9238 9463 \begin_layout Subsection 9239 9464 Další verze bikriteriální metody 9240 \end_layout9241 9242 \begin_layout Itemize9243 posun hodnot pro správné časy9244 9465 \end_layout 9245 9466 … … 9280 9501 Injektáž ( 9281 9502 \series medium 9282 nutno přepracovat 9503 poupravit dle skutečné implementace 9283 9504 \series default 9284 9505 ) … … 9286 9507 9287 9508 \begin_layout Standard 9288 V této části bude popsán jednoduchý návrh řízení využívajícího injektáží. 9289 Jedná se o velmi základní návrh. 9290 \end_layout 9291 9292 \begin_layout Standard 9293 Základní myšlenka je následující: Pomocí techniky injektáží se nepodařilo 9294 získat dostatečně kvalitní odhad úhlu natočení 9509 V této části bude popsán návrh regulátoru využívajícího vysokofrekvenční 9510 injektáže. 9511 Základní myšlenka je následující: Protože se pomocí techniky injektáží 9512 se nepodařilo získat dostatečně kvalitní odhad úhlu natočení 9295 9513 \begin_inset Formula $\vartheta$ 9296 9514 \end_inset 9297 9515 9298 , aby byl použit přímo pro řízení. 9299 Je tedy užíváno současně i EKF, kdy odhad 9516 , aby byl použit přímo pro řízení, je užíván současně i EKF, kdy odhad 9300 9517 \begin_inset Formula $\vartheta$ 9301 9518 \end_inset 9302 9519 9303 9520 z injektáže slouží jako další zdroj informace pro EKF. 9304 Kompletní odhad stavu pro řízení pak poskytuje EKF. 9305 Jako řízení je využíváno LQ řízení v 9306 \begin_inset Formula $d-q$ 9307 \end_inset 9308 9309 souřadné soustavě. 9521 Kompletní odhad stavu pro návrh řízení pak poskytuje EKF. 9522 Jako řídící zásahy jsou generovány vektorovým řízením s LQ regulátorem. 9310 9523 \end_layout 9311 9524 … … 9352 9565 \begin_layout Standard 9353 9566 Větší část zde používaných algoritmů (LQ, EKF) již byla popsána výše v textu, 9354 proto zde uvedeme pouze případné změny. 9355 Mění se matice 9567 proto zde uvedeme pouze případné změny, ty jsou v maticích 9356 9568 \begin_inset Formula $C$ 9357 9569 \end_inset … … 9364 9576 \begin_inset Formula 9365 9577 \begin{eqnarray*} 9366 C& = & \left[\begin{array}{ccccc}9578 \overline{C} & = & \left[\begin{array}{ccccc} 9367 9579 1 & 0 & 0 & 0 & 0\\ 9368 9580 0 & 1 & 0 & 0 & 0\\ 9369 9581 0 & 0 & 0 & 1 & 0 9370 9582 \end{array}\right],\\ 9371 R & = & \left[\begin{array}{ccc} 9372 r & 0 & 0\\ 9373 0 & r & 0\\ 9374 0 & 0 & r_{\vartheta} 9375 \end{array}\right]. 9583 \overline{R} & = & \left[\begin{array}{cc} 9584 R & 0\\ 9585 0 & r_{\vartheta} 9586 \end{array}\right], 9376 9587 \end{eqnarray*} 9377 9588 9378 9589 \end_inset 9379 9590 9380 9591 kde 9592 \begin_inset Formula $\overline{C}$ 9593 \end_inset 9594 9595 a 9596 \begin_inset Formula $\overline{R}$ 9597 \end_inset 9598 9599 představují upravené matice a 9600 \begin_inset Formula $R$ 9601 \end_inset 9602 9603 je matice původní, 9604 \begin_inset Formula $0$ 9605 \end_inset 9606 9607 pak představuje nulový blok vhodné velikosti. 9381 9608 \end_layout 9382 9609 … … 9389 9616 \begin_inset CommandInset ref 9390 9617 LatexCommand ref 9391 reference "s ec:Injektáže"9618 reference "sub:Priblizeni-metody-vysokofrekvenc" 9392 9619 9393 9620 \end_inset … … 9397 9624 \end_inset 9398 9625 9399 . 9400 Z tohoto důvodu je třeba upravit i samotný simulátor a založit jej na rovnicích 9401 ( 9402 \begin_inset CommandInset ref 9403 LatexCommand ref 9404 reference "eq:rovnice_ruzneL_proi" 9405 9406 \end_inset 9407 9408 ) a ( 9409 \begin_inset CommandInset ref 9410 LatexCommand ref 9411 reference "eq:rovnice_ruzneL_omega" 9412 9413 \end_inset 9414 9415 ), které uvažují různé indukčnosti. 9416 Pro jednodušší zpracování byly zvoleny indukčnosti 9417 \begin_inset Formula 9418 \begin{eqnarray*} 9419 L_{d} & = & 1.5L_{s},\\ 9420 L_{q} & = & L_{s}. 9421 \end{eqnarray*} 9422 9423 \end_inset 9424 9425 Tato volba samozřejmě neodpovídá SMPMSM, kde je rozdíl indukčností v osách 9426 9427 \begin_inset Formula $d$ 9428 \end_inset 9429 9430 a 9431 \begin_inset Formula $q$ 9432 \end_inset 9433 9434 velmi malý. 9435 Zde užité hodnoty jsou voleny pro usnadnění návrhu. 9626 , platnost tohoto předpokladu byla stanovena v odstavci 9627 \begin_inset CommandInset ref 9628 LatexCommand ref 9629 reference "sec:uloha-rizeni-PMSM" 9630 9631 \end_inset 9632 9633 . 9436 9634 \end_layout 9437 9635 … … 9439 9637 Vysokofrekvenční signál užitý pro injektáž byl zvolen jako kosinový signál 9440 9638 o amplitudě 9441 \begin_inset Formula $ 10V$9639 \begin_inset Formula $A_{inj}$ 9442 9640 \end_inset 9443 9641 9444 9642 a frekvenci 9445 \begin_inset Formula $ 500Hz$9446 \end_inset 9447 9448 . 9449 Volba velikosti amplitudy je opět komplikovanou záležitostí.9643 \begin_inset Formula $\omega_{inj}$ 9644 \end_inset 9645 9646 . 9647 Volba velikosti amplitudy je problematickou záležitostí. 9450 9648 Obecně platí, že větší amplituda umožní snadnější zpracování signálu, především 9451 9649 z důvodu většího odstupu signálu od šumu. 9452 9650 Naopak ale větší amplituda způsobuje i větší rušení v samotném PMSM. 9453 Obvykle je v injektážních technikách užívána amplituda menší, zde zvolená9454 hodnota je vyšší aby opět usnadnila zpracování.9651 Obvykle je v injektážních technikách užívána amplituda menší, řádově v 9652 jednokách voltů. 9455 9653 Dalším problémem může být, že zde předkládaný návrh amplitudu nijak neomezuje 9456 9654 s rostoucími otáčkami, stále je tedy injektován signál o stejné amplitudě. … … 9461 9659 \begin_layout Standard 9462 9660 Asi největší komplikací tohoto přístupu, ale i injektáží obecně je vhodný 9463 návrh low-pass filtru. 9464 Používá se k získání amplitudově modulované informace o poloze rotoru. 9465 Návrh filtrů je obecně netriviální záležitostí a může mít značný dopad 9466 na kvalitu výsledného odhadu 9661 postup demodulace. 9662 Používá se k k tomu obecně digitálních filtrů, jejichž návrh je netriviální 9663 záležitostí a může mít značný dopad na kvalitu výsledného odhadu 9467 9664 \begin_inset Formula $\vartheta$ 9468 9665 \end_inset … … 9487 9684 \begin_inset Formula 9488 9685 \[ 9489 \frac{ V_{hf}}{\omega_{hf}}\frac{L_{q}-L_{d}}{2L_{d}L_{q}}\sin2\theta,9686 \frac{A_{inj}}{\omega_{inj}}\frac{L_{q}-L_{d}}{2L_{d}L_{q}}\sin2\theta, 9490 9687 \] 9491 9688 … … 9493 9690 9494 9691 kde 9495 \begin_inset Formula $ V_{hf}$9692 \begin_inset Formula $A_{inj}$ 9496 9693 \end_inset 9497 9694 9498 9695 představuje amplitudu a 9499 \begin_inset Formula $\omega_{ hf}$9696 \begin_inset Formula $\omega_{inj}$ 9500 9697 \end_inset 9501 9698 … … 9505 9702 9506 9703 je chyba odhadu 9507 \begin_inset Formula $\theta=\vartheta _{sys}-\vartheta_{est}$9704 \begin_inset Formula $\theta=\vartheta-\hat{\vartheta}$ 9508 9705 \end_inset 9509 9706 … … 9520 9717 9521 9718 je však náročný a nedává příliš dobré výsledky z důvodu omezení na jeho 9522 definiční obor, proto je v ýhodné využítaproximaci9719 definiční obor, proto je využita aproximaci 9523 9720 \begin_inset Formula $\sin x\approx x$ 9524 9721 \end_inset … … 9690 9887 9691 9888 \begin_layout Plain Layout 9692 0,19889 1e-1 9693 9890 \end_layout 9694 9891 … … 9699 9896 9700 9897 \begin_layout Plain Layout 9701 0,01 9898 1e-2 9702 9899 \end_layout 9703 9900 … … 9708 9905 9709 9906 \begin_layout Plain Layout 9710 0,001 9907 1e-3 9711 9908 \end_layout 9712 9909 … … 9717 9914 9718 9915 \begin_layout Plain Layout 9719 0,0001 9916 1e-4 9720 9917 \end_layout 9721 9918 … … 9726 9923 9727 9924 \begin_layout Plain Layout 9728 0,00001 9925 1e-5 9729 9926 \end_layout 9730 9927 … … 9735 9932 9736 9933 \begin_layout Plain Layout 9737 0,000001 9934 1e-6 9738 9935 \end_layout 9739 9936 … … 9755 9952 9756 9953 \begin_layout Plain Layout 9757 1,361 1e49954 1,361e4 9758 9955 \end_layout 9759 9956 … … 9764 9961 9765 9962 \begin_layout Plain Layout 9766 2,79 75e39963 2,798e3 9767 9964 \end_layout 9768 9965 … … 9773 9970 9774 9971 \begin_layout Plain Layout 9775 2,713 2e19972 2,713e1 9776 9973 \end_layout 9777 9974 … … 9782 9979 9783 9980 \begin_layout Plain Layout 9784 7,715 2e-29981 7,715e-2 9785 9982 \end_layout 9786 9983 … … 9791 9988 9792 9989 \begin_layout Plain Layout 9793 1,809 0e-19990 1,809e-1 9794 9991 \end_layout 9795 9992 … … 9800 9997 9801 9998 \begin_layout Plain Layout 9802 2,011 2e-19999 2,011e-1 9803 10000 \end_layout 9804 10001 … … 9809 10006 9810 10007 \begin_layout Plain Layout 9811 2,143 4e-110008 2,143e-1 9812 10009 \end_layout 9813 10010 … … 9820 10017 9821 10018 \begin_layout Plain Layout 9822 0,110019 1e-1 9823 10020 \end_layout 9824 10021 … … 9829 10026 9830 10027 \begin_layout Plain Layout 9831 1,36 08e410028 1,361e4 9832 10029 \end_layout 9833 10030 … … 9838 10035 9839 10036 \begin_layout Plain Layout 9840 2,782 1e310037 2,782e3 9841 10038 \end_layout 9842 10039 … … 9847 10044 9848 10045 \begin_layout Plain Layout 9849 2,692 3e110046 2,692e1 9850 10047 \end_layout 9851 10048 … … 9856 10053 9857 10054 \begin_layout Plain Layout 9858 6,10 79e-210055 6,108e-2 9859 10056 \end_layout 9860 10057 … … 9865 10062 9866 10063 \begin_layout Plain Layout 9867 6,52 45e-210064 6,525e-2 9868 10065 \end_layout 9869 10066 … … 9874 10071 9875 10072 \begin_layout Plain Layout 9876 6,82 55e-210073 6,826e-2 9877 10074 \end_layout 9878 10075 … … 9883 10080 9884 10081 \begin_layout Plain Layout 9885 7,2 499e-210082 7,250e-2 9886 10083 \end_layout 9887 10084 … … 9894 10091 9895 10092 \begin_layout Plain Layout 9896 0,01 10093 1e-2 9897 10094 \end_layout 9898 10095 … … 9903 10100 9904 10101 \begin_layout Plain Layout 9905 1,36 09e410102 1,361e4 9906 10103 \end_layout 9907 10104 … … 9912 10109 9913 10110 \begin_layout Plain Layout 9914 2,77 17e310111 2,772e3 9915 10112 \end_layout 9916 10113 … … 9921 10118 9922 10119 \begin_layout Plain Layout 9923 2,6 895e110120 2,690e1 9924 10121 \end_layout 9925 10122 … … 9930 10127 9931 10128 \begin_layout Plain Layout 9932 8,353 1e-210129 8,353e-2 9933 10130 \end_layout 9934 10131 … … 9939 10136 9940 10137 \begin_layout Plain Layout 9941 3,64 47e-210138 3,645e-2 9942 10139 \end_layout 9943 10140 … … 9948 10145 9949 10146 \begin_layout Plain Layout 9950 3,038 0e-210147 3,038e-2 9951 10148 \end_layout 9952 10149 … … 9957 10154 9958 10155 \begin_layout Plain Layout 9959 3,653 2e-210156 3,653e-2 9960 10157 \end_layout 9961 10158 … … 9968 10165 9969 10166 \begin_layout Plain Layout 9970 0,001 10167 1e-3 9971 10168 \end_layout 9972 10169 … … 9977 10174 9978 10175 \begin_layout Plain Layout 9979 1,361 2e410176 1,361e4 9980 10177 \end_layout 9981 10178 … … 9986 10183 9987 10184 \begin_layout Plain Layout 9988 2,76 79e310185 2,768e3 9989 10186 \end_layout 9990 10187 … … 9995 10192 9996 10193 \begin_layout Plain Layout 9997 2,673 4e110194 2,673e1 9998 10195 \end_layout 9999 10196 … … 10004 10201 10005 10202 \begin_layout Plain Layout 10006 9,154 4e-210203 9,154e-2 10007 10204 \end_layout 10008 10205 … … 10013 10210 10014 10211 \begin_layout Plain Layout 10015 3,427 3e-210212 3,427e-2 10016 10213 \end_layout 10017 10214 … … 10022 10219 10023 10220 \begin_layout Plain Layout 10024 3,17 28e-210221 3,173e-2 10025 10222 \end_layout 10026 10223 … … 10031 10228 10032 10229 \begin_layout Plain Layout 10033 3,363 3e-210230 3,363e-2 10034 10231 \end_layout 10035 10232 … … 10042 10239 10043 10240 \begin_layout Plain Layout 10044 0,0001 10241 1e-4 10045 10242 \end_layout 10046 10243 … … 10051 10248 10052 10249 \begin_layout Plain Layout 10053 1,36 16e410250 1,362e4 10054 10251 \end_layout 10055 10252 … … 10060 10257 10061 10258 \begin_layout Plain Layout 10062 2,77 16e310259 2,772e3 10063 10260 \end_layout 10064 10261 … … 10069 10266 10070 10267 \begin_layout Plain Layout 10071 2,67 87e110268 2,679e1 10072 10269 \end_layout 10073 10270 … … 10078 10275 10079 10276 \begin_layout Plain Layout 10080 8,93 57e-210277 8,936e-2 10081 10278 \end_layout 10082 10279 … … 10087 10284 10088 10285 \begin_layout Plain Layout 10089 3,391 1e-210286 3,391e-2 10090 10287 \end_layout 10091 10288 … … 10096 10293 10097 10294 \begin_layout Plain Layout 10098 3,92 16e-210295 3,922e-2 10099 10296 \end_layout 10100 10297 … … 10105 10302 10106 10303 \begin_layout Plain Layout 10107 3,06 16e-210304 3,062e-2 10108 10305 \end_layout 10109 10306 … … 10116 10313 10117 10314 \begin_layout Plain Layout 10118 0,00001 10315 1e-5 10119 10316 \end_layout 10120 10317 … … 10125 10322 10126 10323 \begin_layout Plain Layout 10127 1,36 07e410324 1,361e4 10128 10325 \end_layout 10129 10326 … … 10134 10331 10135 10332 \begin_layout Plain Layout 10136 2,769 4e310333 2,769e3 10137 10334 \end_layout 10138 10335 … … 10143 10340 10144 10341 \begin_layout Plain Layout 10145 2,67 39e110342 2,674e1 10146 10343 \end_layout 10147 10344 … … 10152 10349 10153 10350 \begin_layout Plain Layout 10154 8,48 37e-210351 8,484e-2 10155 10352 \end_layout 10156 10353 … … 10161 10358 10162 10359 \begin_layout Plain Layout 10163 3,484 0e-210360 3,484e-2 10164 10361 \end_layout 10165 10362 … … 10170 10367 10171 10368 \begin_layout Plain Layout 10172 3,924 0e-210369 3,924e-2 10173 10370 \end_layout 10174 10371 … … 10179 10376 10180 10377 \begin_layout Plain Layout 10181 3,765 1e-210378 3,765e-2 10182 10379 \end_layout 10183 10380 … … 10190 10387 10191 10388 \begin_layout Plain Layout 10192 0,000001 10389 1e-6 10193 10390 \end_layout 10194 10391 … … 10199 10396 10200 10397 \begin_layout Plain Layout 10201 1,36 17e410398 1,362e4 10202 10399 \end_layout 10203 10400 … … 10208 10405 10209 10406 \begin_layout Plain Layout 10210 2,77 06e310407 2,771e3 10211 10408 \end_layout 10212 10409 … … 10217 10414 10218 10415 \begin_layout Plain Layout 10219 2,678 3e110416 2,678e1 10220 10417 \end_layout 10221 10418 … … 10226 10423 10227 10424 \begin_layout Plain Layout 10228 8,868 1e-210425 8,868e-2 10229 10426 \end_layout 10230 10427 … … 10235 10432 10236 10433 \begin_layout Plain Layout 10237 3,821 0e-210434 3,821e-2 10238 10435 \end_layout 10239 10436 … … 10244 10441 10245 10442 \begin_layout Plain Layout 10246 3,53 46e-210443 3,535e-2 10247 10444 \end_layout 10248 10445 … … 10253 10450 10254 10451 \begin_layout Plain Layout 10255 3,57 47e-210452 3,575e-2 10256 10453 \end_layout 10257 10454 … … 10325 10522 10326 10523 \begin_layout Plain Layout 10327 0,110524 1e-1 10328 10525 \end_layout 10329 10526 … … 10334 10531 10335 10532 \begin_layout Plain Layout 10336 0,01 10533 1e-2 10337 10534 \end_layout 10338 10535 … … 10343 10540 10344 10541 \begin_layout Plain Layout 10345 0,001 10542 1e-3 10346 10543 \end_layout 10347 10544 … … 10352 10549 10353 10550 \begin_layout Plain Layout 10354 0,0001 10551 1e-4 10355 10552 \end_layout 10356 10553 … … 10361 10558 10362 10559 \begin_layout Plain Layout 10363 0,00001 10560 1e-5 10364 10561 \end_layout 10365 10562 … … 10370 10567 10371 10568 \begin_layout Plain Layout 10372 0,000001 10569 1e-6 10373 10570 \end_layout 10374 10571 … … 10390 10587 10391 10588 \begin_layout Plain Layout 10392 4,444 4e110589 4,444e1 10393 10590 \end_layout 10394 10591 … … 10399 10596 10400 10597 \begin_layout Plain Layout 10401 3,629 0e110598 3,629e1 10402 10599 \end_layout 10403 10600 … … 10408 10605 10409 10606 \begin_layout Plain Layout 10410 1,79 48e110607 1,795e1 10411 10608 \end_layout 10412 10609 … … 10417 10614 10418 10615 \begin_layout Plain Layout 10419 1,82 28e110616 1,823e1 10420 10617 \end_layout 10421 10618 … … 10426 10623 10427 10624 \begin_layout Plain Layout 10428 1,72 69e110625 1,727e1 10429 10626 \end_layout 10430 10627 … … 10435 10632 10436 10633 \begin_layout Plain Layout 10437 1,7 295e110634 1,730e1 10438 10635 \end_layout 10439 10636 … … 10444 10641 10445 10642 \begin_layout Plain Layout 10446 1,14 18e010643 1,142e0 10447 10644 \end_layout 10448 10645 … … 10455 10652 10456 10653 \begin_layout Plain Layout 10457 0,110654 1e-1 10458 10655 \end_layout 10459 10656 … … 10464 10661 10465 10662 \begin_layout Plain Layout 10466 4,444 4e110663 4,444e1 10467 10664 \end_layout 10468 10665 … … 10473 10670 10474 10671 \begin_layout Plain Layout 10475 2,890 1e110672 2,890e1 10476 10673 \end_layout 10477 10674 … … 10482 10679 10483 10680 \begin_layout Plain Layout 10484 1,383 1e110681 1,383e1 10485 10682 \end_layout 10486 10683 … … 10491 10688 10492 10689 \begin_layout Plain Layout 10493 1,58 45e010690 1,585e0 10494 10691 \end_layout 10495 10692 … … 10500 10697 10501 10698 \begin_layout Plain Layout 10502 1,8 896e110699 1,890e1 10503 10700 \end_layout 10504 10701 … … 10509 10706 10510 10707 \begin_layout Plain Layout 10511 2,333 3e010708 2,333e0 10512 10709 \end_layout 10513 10710 … … 10518 10715 10519 10716 \begin_layout Plain Layout 10520 6,725 0e-210717 6,725e-2 10521 10718 \end_layout 10522 10719 … … 10529 10726 10530 10727 \begin_layout Plain Layout 10531 0,01 10728 1e-2 10532 10729 \end_layout 10533 10730 … … 10538 10735 10539 10736 \begin_layout Plain Layout 10540 4,44 45e110737 4,445e1 10541 10738 \end_layout 10542 10739 … … 10547 10744 10548 10745 \begin_layout Plain Layout 10549 1,141 1e110746 1,141e1 10550 10747 \end_layout 10551 10748 … … 10556 10753 10557 10754 \begin_layout Plain Layout 10558 1,83 47e010755 1,835e0 10559 10756 \end_layout 10560 10757 … … 10565 10762 10566 10763 \begin_layout Plain Layout 10567 1,181 4e010764 1,181e0 10568 10765 \end_layout 10569 10766 … … 10574 10771 10575 10772 \begin_layout Plain Layout 10576 5,37 88e010773 5,379e0 10577 10774 \end_layout 10578 10775 … … 10583 10780 10584 10781 \begin_layout Plain Layout 10585 5,89 28e010782 5,893e0 10586 10783 \end_layout 10587 10784 … … 10592 10789 10593 10790 \begin_layout Plain Layout 10594 5,474 4e-310791 5,474e-3 10595 10792 \end_layout 10596 10793 … … 10603 10800 10604 10801 \begin_layout Plain Layout 10605 0,001 10802 1e-3 10606 10803 \end_layout 10607 10804 … … 10612 10809 10613 10810 \begin_layout Plain Layout 10614 4,44 45e110811 4,445e1 10615 10812 \end_layout 10616 10813 … … 10621 10818 10622 10819 \begin_layout Plain Layout 10623 1,826 3e110820 1,826e1 10624 10821 \end_layout 10625 10822 … … 10630 10827 10631 10828 \begin_layout Plain Layout 10632 1,99 87e010829 1,999e0 10633 10830 \end_layout 10634 10831 … … 10639 10836 10640 10837 \begin_layout Plain Layout 10641 1,179 2e010838 1,179e0 10642 10839 \end_layout 10643 10840 … … 10648 10845 10649 10846 \begin_layout Plain Layout 10650 1,260 0e010847 1,260e0 10651 10848 \end_layout 10652 10849 … … 10657 10854 10658 10855 \begin_layout Plain Layout 10659 1,12 46e010856 1,125e0 10660 10857 \end_layout 10661 10858 … … 10666 10863 10667 10864 \begin_layout Plain Layout 10668 5,79 68e-310865 5,797e-3 10669 10866 \end_layout 10670 10867 … … 10677 10874 10678 10875 \begin_layout Plain Layout 10679 0,0001 10876 1e-4 10680 10877 \end_layout 10681 10878 … … 10686 10883 10687 10884 \begin_layout Plain Layout 10688 4,444 4e110885 4,444e1 10689 10886 \end_layout 10690 10887 … … 10695 10892 10696 10893 \begin_layout Plain Layout 10697 1,897 0e110894 1,897e1 10698 10895 \end_layout 10699 10896 … … 10704 10901 10705 10902 \begin_layout Plain Layout 10706 7,45 75e010903 7,458e0 10707 10904 \end_layout 10708 10905 … … 10713 10910 10714 10911 \begin_layout Plain Layout 10715 6,752 2e010912 6,752e0 10716 10913 \end_layout 10717 10914 … … 10722 10919 10723 10920 \begin_layout Plain Layout 10724 5,838 4e010921 5,838e0 10725 10922 \end_layout 10726 10923 … … 10731 10928 10732 10929 \begin_layout Plain Layout 10733 5,782 2e010930 5,782e0 10734 10931 \end_layout 10735 10932 … … 10740 10937 10741 10938 \begin_layout Plain Layout 10742 4,013 1e010939 4,013e0 10743 10940 \end_layout 10744 10941 … … 10751 10948 10752 10949 \begin_layout Plain Layout 10753 0,00001 10950 1e-5 10754 10951 \end_layout 10755 10952 … … 10760 10957 10761 10958 \begin_layout Plain Layout 10762 4,444 4e110959 4,444e1 10763 10960 \end_layout 10764 10961 … … 10769 10966 10770 10967 \begin_layout Plain Layout 10771 1,16 25e110968 1,163e1 10772 10969 \end_layout 10773 10970 … … 10778 10975 10779 10976 \begin_layout Plain Layout 10780 1,76 69e010977 1,767e0 10781 10978 \end_layout 10782 10979 … … 10787 10984 10788 10985 \begin_layout Plain Layout 10789 1,1 195e010986 1,120e0 10790 10987 \end_layout 10791 10988 … … 10796 10993 10797 10994 \begin_layout Plain Layout 10798 1,184 4e010995 1,184e0 10799 10996 \end_layout 10800 10997 … … 10805 11002 10806 11003 \begin_layout Plain Layout 10807 1,044 3e011004 1,044e0 10808 11005 \end_layout 10809 11006 … … 10814 11011 10815 11012 \begin_layout Plain Layout 10816 1,700 2e-211013 1,700e-2 10817 11014 \end_layout 10818 11015 … … 10825 11022 10826 11023 \begin_layout Plain Layout 10827 0,000001 11024 1e-6 10828 11025 \end_layout 10829 11026 … … 10834 11031 10835 11032 \begin_layout Plain Layout 10836 4,444 4e111033 4,444e1 10837 11034 \end_layout 10838 11035 … … 10843 11040 10844 11041 \begin_layout Plain Layout 10845 1,07 88e111042 1,079e1 10846 11043 \end_layout 10847 11044 … … 10852 11049 10853 11050 \begin_layout Plain Layout 10854 1,97 59e011051 1,976e0 10855 11052 \end_layout 10856 11053 … … 10861 11058 10862 11059 \begin_layout Plain Layout 10863 1,212 2e011060 1,212e0 10864 11061 \end_layout 10865 11062 … … 10870 11067 10871 11068 \begin_layout Plain Layout 10872 1,712 0e011069 1,712e0 10873 11070 \end_layout 10874 11071 … … 10879 11076 10880 11077 \begin_layout Plain Layout 10881 1,16 76e011078 1,168e0 10882 11079 \end_layout 10883 11080 … … 10888 11085 10889 11086 \begin_layout Plain Layout 10890 1,740 3e-211087 1,740e-2 10891 11088 \end_layout 10892 11089 … … 10987 11184 10988 11185 \begin_layout Plain Layout 10989 3,115 1e-211186 3,115e-2 10990 11187 \end_layout 10991 11188 … … 10996 11193 10997 11194 \begin_layout Plain Layout 10998 2,834 2e-211195 2,834e-2 10999 11196 \end_layout 11000 11197 … … 11005 11202 11006 11203 \begin_layout Plain Layout 11007 3,23 26e-211204 3,233e-2 11008 11205 \end_layout 11009 11206 … … 11025 11222 11026 11223 \begin_layout Plain Layout 11027 3,03 35e-211224 3,034e-2 11028 11225 \end_layout 11029 11226 … … 11034 11231 11035 11232 \begin_layout Plain Layout 11036 3,288 4e-211233 3,288e-2 11037 11234 \end_layout 11038 11235 … … 11043 11240 11044 11241 \begin_layout Plain Layout 11045 2,731 2e-211242 2,731e-2 11046 11243 \end_layout 11047 11244 … … 11063 11260 11064 11261 \begin_layout Plain Layout 11065 2,84 57e-211262 2,846e-2 11066 11263 \end_layout 11067 11264 … … 11072 11269 11073 11270 \begin_layout Plain Layout 11074 2,94 18e-211271 2,942e-2 11075 11272 \end_layout 11076 11273 … … 11081 11278 11082 11279 \begin_layout Plain Layout 11083 3,518 3e-211280 3,518e-2 11084 11281 \end_layout 11085 11282 … … 11101 11298 11102 11299 \begin_layout Plain Layout 11103 3,099 4e-211300 3,099e-2 11104 11301 \end_layout 11105 11302 … … 11110 11307 11111 11308 \begin_layout Plain Layout 11112 2,537 0e-211309 2,537e-2 11113 11310 \end_layout 11114 11311 … … 11119 11316 11120 11317 \begin_layout Plain Layout 11121 2,579 3e-211318 2,579e-2 11122 11319 \end_layout 11123 11320 … … 11139 11336 11140 11337 \begin_layout Plain Layout 11141 3,03 75e-211338 3,038e-2 11142 11339 \end_layout 11143 11340 … … 11148 11345 11149 11346 \begin_layout Plain Layout 11150 2,917 0e-211347 2,917e-2 11151 11348 \end_layout 11152 11349 … … 11157 11354 11158 11355 \begin_layout Plain Layout 11159 2,73 15e-211356 2,732e-2 11160 11357 \end_layout 11161 11358 … … 11177 11374 11178 11375 \begin_layout Plain Layout 11179 4,202 4e-211376 4,202e-2 11180 11377 \end_layout 11181 11378 … … 11186 11383 11187 11384 \begin_layout Plain Layout 11188 4,385 1e-211385 4,385e-2 11189 11386 \end_layout 11190 11387 … … 11195 11392 11196 11393 \begin_layout Plain Layout 11197 5,556 1e-211394 5,556e-2 11198 11395 \end_layout 11199 11396 … … 11215 11412 11216 11413 \begin_layout Plain Layout 11217 3,44 87e-211414 3,449e-2 11218 11415 \end_layout 11219 11416 … … 11224 11421 11225 11422 \begin_layout Plain Layout 11226 4,442 3e-211423 4,442e-2 11227 11424 \end_layout 11228 11425 … … 11233 11430 11234 11431 \begin_layout Plain Layout 11235 3,709 3e-211432 3,709e-2 11236 11433 \end_layout 11237 11434 … … 11313 11510 11314 11511 \begin_layout Plain Layout 11315 9,83 46e-111512 9,835e-1 11316 11513 \end_layout 11317 11514 … … 11322 11519 11323 11520 \begin_layout Plain Layout 11324 2,225 3e011521 2,225e0 11325 11522 \end_layout 11326 11523 … … 11331 11528 11332 11529 \begin_layout Plain Layout 11333 2,8 597e-211530 2,860e-2 11334 11531 \end_layout 11335 11532 … … 11351 11548 11352 11549 \begin_layout Plain Layout 11353 3,332 2e-211550 3,332e-2 11354 11551 \end_layout 11355 11552 … … 11360 11557 11361 11558 \begin_layout Plain Layout 11362 3,826 3e-211559 3,826e-2 11363 11560 \end_layout 11364 11561 … … 11369 11566 11370 11567 \begin_layout Plain Layout 11371 4,14 18e-211568 4,142e-2 11372 11569 \end_layout 11373 11570 … … 11501 11698 11502 11699 \begin_layout Plain Layout 11503 3,531 3e-311700 3,531e-3 11504 11701 \end_layout 11505 11702 … … 11510 11707 11511 11708 \begin_layout Plain Layout 11512 4,03 39e-311709 4,034e-3 11513 11710 \end_layout 11514 11711 … … 11519 11716 11520 11717 \begin_layout Plain Layout 11521 3,695 0e-311718 3,695e-3 11522 11719 \end_layout 11523 11720 … … 11539 11736 11540 11737 \begin_layout Plain Layout 11541 4,69 29e-311738 4,693e-3 11542 11739 \end_layout 11543 11740 … … 11548 11745 11549 11746 \begin_layout Plain Layout 11550 4,499 0e-311747 4,499e-3 11551 11748 \end_layout 11552 11749 … … 11557 11754 11558 11755 \begin_layout Plain Layout 11559 4,699 3e-311756 4,699e-3 11560 11757 \end_layout 11561 11758 … … 11577 11774 11578 11775 \begin_layout Plain Layout 11579 1,345 3e-211776 1,345e-2 11580 11777 \end_layout 11581 11778 … … 11586 11783 11587 11784 \begin_layout Plain Layout 11588 1,250 3e-211785 1,250e-2 11589 11786 \end_layout 11590 11787 … … 11595 11792 11596 11793 \begin_layout Plain Layout 11597 1,29 49e-211794 1,295e-2 11598 11795 \end_layout 11599 11796 … … 11615 11812 11616 11813 \begin_layout Plain Layout 11617 7,04 18e-311814 7,042e-3 11618 11815 \end_layout 11619 11816 … … 11624 11821 11625 11822 \begin_layout Plain Layout 11626 5,06 45e-311823 5,065e-3 11627 11824 \end_layout 11628 11825 … … 11633 11830 11634 11831 \begin_layout Plain Layout 11635 7,09 85e-311832 7,099e-3 11636 11833 \end_layout 11637 11834 … … 11653 11850 11654 11851 \begin_layout Plain Layout 11655 6,667 0e211852 6,667e2 11656 11853 \end_layout 11657 11854 … … 11662 11859 11663 11860 \begin_layout Plain Layout 11664 7,26 89e211861 7,269e2 11665 11862 \end_layout 11666 11863 … … 11671 11868 11672 11869 \begin_layout Plain Layout 11673 6,653 4e211870 6,653e2 11674 11871 \end_layout 11675 11872 … … 11691 11888 11692 11889 \begin_layout Plain Layout 11693 8,620 2e111890 8,620e1 11694 11891 \end_layout 11695 11892 … … 11700 11897 11701 11898 \begin_layout Plain Layout 11702 8,582 2e111899 8,582e1 11703 11900 \end_layout 11704 11901 … … 11709 11906 11710 11907 \begin_layout Plain Layout 11711 9,314 1e111908 9,314e1 11712 11909 \end_layout 11713 11910 … … 11814 12011 \begin{eqnarray*} 11815 12012 \frac{di_{\alpha}}{dt} & = & \cos\vartheta\left(\frac{u_{d}}{L_{d}}-\frac{R_{s}i_{d}}{L_{d}}-i_{\beta}\omega\cos\vartheta+i_{\alpha}\omega\sin\vartheta+\frac{L_{q}\omega i_{q}}{L_{d}}\right)\\ 11816 & - & \sin\vartheta\left(\frac{u_{q}}{L_{q}}-\frac{R_{s}i_{q}}{L_{q}}+i_{\alpha}\omega\cos\vartheta+i_{\beta}\omega\sin\vartheta-\frac{\omega\psi_{pm}}{L_{q}}-\frac{L_{d}\omega i_{d}}{L_{q}}\right) 12013 & - & \sin\vartheta\left(\frac{u_{q}}{L_{q}}-\frac{R_{s}i_{q}}{L_{q}}+i_{\alpha}\omega\cos\vartheta+i_{\beta}\omega\sin\vartheta-\frac{\omega\psi_{pm}}{L_{q}}-\frac{L_{d}\omega i_{d}}{L_{q}}\right), 11817 12014 \end{eqnarray*} 11818 12015 … … 11826 12023 \begin{eqnarray*} 11827 12024 \frac{di_{\beta}}{dt} & = & \cos\vartheta\left(\frac{u_{q}}{L_{q}}-\frac{R_{s}i_{q}}{L_{q}}+i_{\alpha}\omega\cos\vartheta+i_{\beta}\omega\sin\vartheta-\frac{\omega\psi_{pm}}{L_{q}}+\frac{L_{d}\omega i_{d}}{L_{q}}\right)\\ 11828 & + & \sin\vartheta\left(\frac{u_{d}}{L_{d}}-\frac{R_{s}i_{d}}{L_{d}}-i_{\beta}\omega\cos\vartheta+i_{\alpha}\omega\sin\vartheta+\frac{L_{q}\omega i_{q}}{L_{d}}\right) 12025 & + & \sin\vartheta\left(\frac{u_{d}}{L_{d}}-\frac{R_{s}i_{d}}{L_{d}}-i_{\beta}\omega\cos\vartheta+i_{\alpha}\omega\sin\vartheta+\frac{L_{q}\omega i_{q}}{L_{d}}\right), 11829 12026 \end{eqnarray*} 11830 12027 … … 11837 12034 \begin_inset Formula 11838 12035 \begin{eqnarray*} 11839 \frac{d\omega}{dt} & = & \frac{k_{p}p_{p}^{2}\left(\psi_{pm}i_{q}+\left(L_{d}-L_{q}\right)i_{d}i_{q}\right)}{J}-\frac{B\omega}{J} 12036 \frac{d\omega}{dt} & = & \frac{k_{p}p_{p}^{2}\left(\psi_{pm}i_{q}+\left(L_{d}-L_{q}\right)i_{d}i_{q}\right)}{J}-\frac{B\omega}{J}, 11840 12037 \end{eqnarray*} 11841 12038 … … 11845 12042 \begin_inset Formula 11846 12043 \begin{eqnarray*} 11847 i_{d} & = & i_{\alpha}\cos\vartheta+i_{\beta}\sin\vartheta \\11848 i_{q} & = & i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta \\11849 u_{d} & = & u_{\alpha}\cos\vartheta+u_{\beta}\sin\vartheta \\11850 u_{q} & = & u_{\beta}\cos\vartheta-u_{\alpha}\sin\vartheta 12044 i_{d} & = & i_{\alpha}\cos\vartheta+i_{\beta}\sin\vartheta,\\ 12045 i_{q} & = & i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta,\\ 12046 u_{d} & = & u_{\alpha}\cos\vartheta+u_{\beta}\sin\vartheta,\\ 12047 u_{q} & = & u_{\beta}\cos\vartheta-u_{\alpha}\sin\vartheta. 11851 12048 \end{eqnarray*} 11852 12049 … … 12024 12221 \frac{d\omega_{t+1}}{di_{\alpha,t}} & \frac{d\omega_{t+1}}{di_{\beta,t}} & \frac{d\omega_{t+1}}{d\omega_{t}} & \frac{d\omega_{t+1}}{d\vartheta_{t}}\\ 12025 12222 \frac{d\vartheta_{t+1}}{di_{\alpha,t}} & \frac{d\vartheta_{t+1}}{di_{\beta,t}} & \frac{d\vartheta_{t+1}}{d\omega_{t}} & \frac{d\vartheta_{t+1}}{d\vartheta_{t}} 12026 \end{array}\right] 12223 \end{array}\right]. 12027 12224 \] 12028 12225 … … 12044 12241 v důsledku linearizace a řízení na nenulové požadované otáčky, tedy substituce 12045 12242 12046 \begin_inset Formula $\psi=\omega-\overline{\omega}$ 12243 \begin_inset CommandInset ref 12244 LatexCommand ref 12245 reference "eq:substituce-psi" 12246 12047 12247 \end_inset 12048 12248 … … 12054 12254 A_{t} & \gamma\\ 12055 12255 0 & 1 12056 \end{array}\right] 12256 \end{array}\right], 12057 12257 \] 12058 12258 … … 12087 12287 & - & 2L_{d}R_{s}i_{\alpha}\vartheta\sin2\vartheta+2L_{q}R_{s}i_{\alpha}\vartheta\sin2\vartheta+2L_{d}^{2}i_{\alpha}\psi\vartheta\cos2\vartheta-2L_{q}^{2}i_{\alpha}\psi\vartheta\cos2\vartheta+\\ 12088 12288 & + & 2L_{d}^{2}i_{\alpha}\overline{\omega}\vartheta\cos2\vartheta-2L_{q}^{2}i_{\alpha}\overline{\omega}\vartheta\cos2\vartheta+2L_{d}^{2}i_{\beta}\psi\vartheta\sin2\vartheta-2L_{q}^{2}i_{\beta}\psi\vartheta\sin2\vartheta+\\ 12089 & + & \left.2L_{d}^{2}i_{\beta}\overline{\omega}\vartheta\sin2\vartheta-2L_{q}^{2}i_{\beta}\overline{\omega}\vartheta\sin2\vartheta+2L_{d}\psi_{pm}\psi\vartheta\cos\vartheta+2L_{d}\psi_{pm}\overline{\omega}\vartheta\cos\vartheta\right) 12289 & + & \left.2L_{d}^{2}i_{\beta}\overline{\omega}\vartheta\sin2\vartheta-2L_{q}^{2}i_{\beta}\overline{\omega}\vartheta\sin2\vartheta+2L_{d}\psi_{pm}\psi\vartheta\cos\vartheta+2L_{d}\psi_{pm}\overline{\omega}\vartheta\cos\vartheta\right), 12090 12290 \end{eqnarray*} 12091 12291 … … 12102 12302 & - & 2L_{d}R_{s}i_{\beta}\vartheta\sin2\vartheta+2L_{q}R_{s}i_{\beta}\vartheta\sin2\vartheta+2L_{d}^{2}i_{\beta}\psi\vartheta\cos2\vartheta-2L_{q}^{2}i_{\beta}\psi\vartheta\cos2\vartheta+\\ 12103 12303 & + & 2L_{d}^{2}i_{\beta}\overline{\omega}\vartheta\cos2\vartheta-2L_{q}^{2}i_{\beta}\overline{\omega}\vartheta\cos2\vartheta-2L_{d}^{2}i_{\alpha}\psi\vartheta\sin2\vartheta+2L_{q}^{2}i_{\alpha}\psi\vartheta\sin2\vartheta-\\ 12104 & - & \left.2L_{d}^{2}i_{\alpha}\overline{\omega}\vartheta\sin2\vartheta+2L_{q}^{2}i_{\alpha}\overline{\omega}\vartheta\sin2\vartheta-2L_{d}\psi_{pm}\psi\vartheta\sin\vartheta-2L_{d}\psi_{pm}\overline{\omega}\vartheta\sin\vartheta\right) 12304 & - & \left.2L_{d}^{2}i_{\alpha}\overline{\omega}\vartheta\sin2\vartheta+2L_{q}^{2}i_{\alpha}\overline{\omega}\vartheta\sin2\vartheta-2L_{d}\psi_{pm}\psi\vartheta\sin\vartheta-2L_{d}\psi_{pm}\overline{\omega}\vartheta\sin\vartheta\right), 12105 12305 \end{eqnarray*} 12106 12306 … … 12114 12314 & - & 2L_{d}i_{\alpha}i_{\beta}k_{p}p_{p}^{2}\cos2\vartheta+2L_{q}i_{\alpha}i_{\beta}k_{p}p_{p}^{2}\cos2\vartheta+2L_{d}i_{\alpha}^{2}k_{p}p_{p}^{2}\vartheta\cos2\vartheta-\\ 12115 12315 & - & 2L_{d}i_{\beta}^{2}k_{p}p_{p}^{2}\vartheta\cos2\vartheta-2L_{q}i_{\alpha}^{2}k_{p}p_{p}^{2}\vartheta\cos2\vartheta+2L_{q}i_{\beta}^{2}k_{p}p_{p}^{2}\vartheta\cos2\vartheta+\\ 12116 & + & 4L_{d}i_{\alpha}i_{\beta}k_{p}p_{p}^{2}\vartheta\sin2\vartheta-4L_{q}i_{\alpha}i_{\beta}k_{p}p_{p}^{2}\vartheta\sin2\vartheta12316 & + & \left.4L_{d}i_{\alpha}i_{\beta}k_{p}p_{p}^{2}\vartheta\sin2\vartheta-4L_{q}i_{\alpha}i_{\beta}k_{p}p_{p}^{2}\vartheta\sin2\vartheta\right), 12117 12317 \end{eqnarray*} 12118 12318 … … 12122 12322 \begin_inset Formula 12123 12323 \[ 12124 \gamma_{4}=\Delta t\overline{\omega} 12324 \gamma_{4}=\Delta t\overline{\omega}. 12125 12325 \] 12126 12326 … … 12135 12335 \end_inset 12136 12336 12137 , ta je nyní závislá na čase a má zápis ve tvaru12337 , ta je nyní závislá na čase a má zápis 12138 12338 \begin_inset Formula 12139 12339 \[ … … 12144 12344 0 & 0\\ 12145 12345 0 & 0 12146 \end{array}\right] 12346 \end{array}\right]. 12147 12347 \] 12148 12348 -
TabularUnified applications/dual/vahala/DP/vyz_texty.bib ¶
r1446 r1448 72 72 note = {{UNC-Chapel Hill}}, 73 73 year = {2006}, 74 howpublished = {[cit. 2012-04-08] Dostupn �: http://www.cs.unc.edu},74 howpublished = {[cit. 2012-04-08] Dostupn\'{e} z: http://www.cs.unc.edu}, 75 75 } 76 76 … … 80 80 note= {{Institute for Systems and Robotics, Lisboa Technical University}}, 81 81 year = {2004}, 82 howpublished = {[cit. 2012-04-08] Dostupn �: http://users.isr.ist.utl.pt},82 howpublished = {[cit. 2012-04-08] Dostupn\'{e} z: http://users.isr.ist.utl.pt}, 83 83 } 84 85 @MISC{utiaBDM2005, 86 author = {I. Nagy and L. Pavelkov\'{a} and E. Suzdaleva and aj.}, 87 title = {Bayesian decision making: Theory and examples}, 88 note= {{\'{U}TIA AV\v{C}R, Prague}}, 89 year = {2005}, 90 howpublished = {[cit. 2012-04-10] Dostupn\'{e} z: http://www.utia.cas.cz/cs/AS}, 91 }