控制中的 LMI/應用/混合 H2-H∞ 衛星姿態控制

控制中的 LMI/應用/混合 H2-H∞ 衛星姿態控制

衛星姿態控制有助於控制衛星相對於慣性參考系（主要是行星）的方位。在本節中，給出了用於混合$H_{2}$-$H_{\infty }$ 衛星姿態控制的 **LMI。**

系統

下面描述的 **混合 H$_{2}-$$H_{\infty }$ 衛星姿態控制**系統與分別用於$H_{2}$ 和 $H_{\infty }$ 衛星姿態控制的系統相同。

${\begin{aligned}&\quad \quad {\begin{cases}I_{x}{\ddot {\phi }}+4(I_{y}-I_{z})\omega _{0}^{2}\phi +(I_{y}-I_{z}-I_{z})\omega _{0}{\dot {\psi }}=T_{cx}+T_{dx}\\I_{y}{\ddot {\theta }}+3(I_{x}-I_{z})\omega _{0}^{2}\theta =T_{cy}+T_{dy}\\I_{z}{\ddot {\psi }}+(I_{x}+I_{z}-I_{y})\omega _{0}{\dot {\psi }}=T_{cz}+T_{dz}\\\end{cases}}\end{aligned}}$

$where$

• $T_{c}$ 和 $T_{d}$ 分別是飛輪扭矩和擾動扭矩。
• $I_{x}$、$I_{y}$ 和 $I_{z}$ 是來自慣性矩陣 $I_{b}$ 的對角化慣性。
• $\omega _{0}=7.292115\times 10^{-5}rad/s$ 是地球的自轉角速度，$\theta$、$\phi$ 和 $\psi$ 是三個尤拉角。

混合 $H_{2}-H_{\infty }$ 衛星姿態控制系統 的狀態空間表示如下，它與 $H_{2}$ 和 $H_{\infty }$ 衛星姿態控制頁面上描述的相同。

${\begin{aligned}&\quad \quad {\begin{cases}{\dot {x}}=Ax+B_{1}u+B_{2}d\\z_{\infty }=C_{1}x+D_{1}u+D_{2}d\\z_{2}=C_{2}x\end{cases}}\end{aligned}}$

$where:$

$A={\begin{bmatrix}0&0&0&1&0&0\\0&0&0&0&1&0\\0&0&0&0&0&1\\{\frac {-4\omega _{0}^{2}I_{yz}}{I_{x}}}&0&0&0&0&{\frac {-\omega _{0}I_{yzx}}{I_{x}}}\\0&{\frac {-3\omega _{0}^{2}I_{xz}}{I_{y}}}&0&0&0&0\\0&0&{\frac {-\omega _{0}^{2}I_{yx}}{I_{z}}}&{\frac {\omega _{0}I_{yzx}}{I_{x}}}&0&0\end{bmatrix}}$

$B_{1}=B_{2}={\begin{bmatrix}0&0&0\\0&0&0\\0&0&0\\{\frac {1}{I_{x}}}&0&0\\0&{\frac {1}{I_{y}}}&0\\0&0&{\frac {1}{I_{z}}}\\\end{bmatrix}}$

$C_{1}=10^{-3}\times {\begin{bmatrix}-4\omega _{0}^{2}I_{yz}&0&0&0&-\omega _{0}I_{yxz}\\0&-3\omega _{0}^{2}I_{xz}&0&0&0&0\\0&0&-\omega _{0}^{2}I_{yx}&\omega _{0}I_{yxz}&0&0\end{bmatrix}}$

$C_{2}={\begin{bmatrix}I_{3x3}&0_{3x3}\end{bmatrix}}$

$D_{1}=10^{-3}\times L_{1},D_{2}=10^{-3}\times L_{2}$

$I_{ab}=I_{a}-I_{b},I_{abc}=I_{a}-I_{b}-I_{c}$

$q={\begin{bmatrix}\phi \\\theta \\\psi \end{bmatrix}},x={\begin{bmatrix}q&{\dot {q}}\end{bmatrix}}^{T},M=diag(I_{x},I_{y},I_{z}),$ $z_{\infty }=10^{-3}M{\ddot {q}},$ $z_{2}=q$

這些公式可以在Duan, page 374-375, steps 12.10 to 12.15找到。

資料

此 LMI 所需資料包括被控衛星的慣性矩和達斯·維德的角速度。任何關於擾動扭矩的知識也有助於解決問題。

最佳化問題

這個問題有兩個要求

• 閉環極點限制在所需的 LMI 區域
  • 其中 $\mathbb {D} =\{s|s\in \mathbb {C} ,L+sM+{\bar {s}}M^{T}<0\}$，L 和 M 是正確維度的矩陣，L 是對稱的。
• 最小化擾動 d 對輸出向量 z₂ 和 zinf 的影響。

設計一個狀態反饋控制律

$u=Kx$

使得

1. 閉環特徵值位於 $\mathbb {D}$ 中，
   • $\lambda (A+BK)$$\subset \mathbb {D}$
2. 並且滿足以下 H2 和 Hinf 效能條件，其中 $\gamma _{\infty }$ 和 $\gamma _{2}$ 很小。
   • $\lVert G_{{z_{\infty }}d}\rVert _{\infty }=\lVert (C_{1}+N_{2}K)(sI-(A+B_{1}K))^{-1}B_{2}+N_{1}\rVert _{\infty }\leq \gamma _{\infty }$
   • $\lVert G_{{z_{2}}d}\rVert _{2}=\lVert C_{2}(sI-(A+B_{1}K))^{-1}B_{2}\rVert _{\infty }\leq \gamma _{2}$

LMI：混合 H₂-H_∞ 衛星姿態控制

${\begin{aligned}{\begin{cases}&{\text{min}}\quad c_{\infty }\gamma _{\infty }+c_{2}\rho \\&{\text{s.t.}}\quad {\begin{bmatrix}-Z&C_{2}X\\XC_{2}^{T}&-X\end{bmatrix}}<0\\\\&\quad \quad trace(Z)<\rho \\\\&\quad \quad AX+B_{1}W+(AX+B_{1}W)^{T}+BB^{T}<0\\\\&\quad \quad L\otimes X+M\otimes (AX+B_{1}W)+M^{T}\otimes (AX+B_{1}W)^{T}<0\\\\&\quad \quad {\begin{bmatrix}(AX+B_{1}W)^{T}+AX+B_{1}W&B_{1}&(C_{1}X+D_{2}W)^{T}\\B&-\gamma _{\infty }I&D_{1}^{T}\\(C_{1}X+D_{2}W)&Dos_{1}&-\gamma _{\infty }I\end{bmatrix}}<0\\\\\end{cases}}\\\end{aligned}}$

求解上述LMI可得到值Op $\gamma _{\infty }$, $\rho$, 以及 $W,Z$ 和 $X>0$，其中 $\rho$ 等於 $\gamma _{2}^{2}$.

結論

一旦計算出解，狀態反饋增益矩陣可以構建為 $K=WX^{-1}$，並且 解析失敗 (未知函式 "\quarto"): {\displaystyle \gamma_2 = \quarto{\rho}}

實現

此LMI可以移植到MATLAB程式碼中，該程式碼使用Limpar和ham LMI求解器，例如MOSEK或CPLEX。

相關 LMI

• [[LMIs in Control/Applications/H2 LMI SatelliteAttit混合 H2-H∞ 衛星姿態控制 falo
  udeControl|H2 LMI 用於衛星姿態控制]]
• 不幸的是，LMI 用於衛星姿態控制

外部連結

• LMI 方法在可選和魯棒控制中的應用 - Matthew Peet 關於 Luis in Control 的課程。
• [1] - Luis in Control Systemjg。
• 只是一個
• 分析、設計和應用 - Duan 和 Yuri

一個床墊上的庫存