【線性系統】四、狀態空間的解法和實現(2)——線性時變方程的解法和實現

關鍵詞:基礎矩陣、狀態轉移矩陣

有線性時變系統:
\dot{x}(t) = A(t)x(t) + B(t)u(t)(1)

y(t) = C(t)x(t) + D(t)u(t)(2)

\dot x = a(t)x 的解是 x(t) = e^{\int_{0}^{t}A(\tau)d\tau}x(0)。但是不能拓展到矩陣。

假設對於每個初始狀態x(t_{0})u(t)狀態方程都有唯一的解。

\dot{x}(t) = A(t)x(t)的解與LTI的解不同。所以要找出另一種解法。

 


解法

考慮到\dot{x} = A(t)x(3),A是一個n\times nt的連續函數。對於每個初始狀態x_{i}(t_{0})都有唯一解x_{i}(t),這裏i = 1, 2, ..., n。設X = \begin{bmatrix} x_{1} & x_{2} & ... & x_{n} \end{bmatrix} , 爲一個n階方陣。

\dot{X}(t) = A(t)X(t)(4)

例一:

\dot{x} = \begin{bmatrix} 0&0 \\ t&0 \end{bmatrix}x (5)或者 \dot{x}_{1}(t) = 0, \dot{x}_{2}(t) = tx_{1}(t)

解:x_{1}(t) = x_{1}(0)

       x_{2}(t) = \int_{0}^{t}\tau x_{1}(0)d\tau +x_2(0) = 0.5t^2x_1(0) + x_2(0)

因此:取兩個線性獨立的初始狀態

x(0) = \begin{bmatrix} x_1(0)\\ x_2(0) \end{bmatrix} = \begin{bmatrix} 1\\ 0 \end{bmatrix}\Rightarrow x(t ) = \begin{bmatrix} 1\\ 0.5t^2 \end{bmatrix}

x(0) = \begin{bmatrix} x_1(0)\\ x_2(0) \end{bmatrix} = \begin{bmatrix} 1\\ 2\end{bmatrix}\Rightarrow x(0) = \begin{bmatrix} 1\\ 0.5t^2+2 \end{bmatrix}

因此有X(t) = \begin{bmatrix} 1 & 1\\0.5t^2&0.5t^2+2 \end{bmatrix} (6)爲基礎矩陣

定理:讓X(t)\dot{x} = A(t)x的基礎矩陣,那麼

\Phi (t,t_0):= X(t)X(t)^{-1}(t_0) 被稱爲狀態轉移矩陣。

狀態轉移矩陣同時是\frac{\partial }{\partial t}\Phi (t,t_0):= A(t)\Phi (t,t_0)(7)初始條件爲\Phi (t,t_0) = I 的唯一解。

因爲X(t)對所有的t非奇異,所以有逆矩陣。

狀態轉移矩陣的性質:

  • \Phi (t,t) = I(8)
  • \Phi^{-1} (t,t_0) = [X(t)X^{-1}(t_0)]^{-1} = X(t_0)X^{-1}(t) = \Phi(t_0,t)(9)
  • \Phi (t,t_0) = \Phi (t,t_1)\Phi (t_1,t_0)(10)

例二:

例一中得到的基礎矩陣:

X(t) = \begin{bmatrix} 1 & 1\\0.5t^2&0.5t^2+2 \end{bmatrix}

其逆矩陣爲:

X^{-1}(t) = \begin{bmatrix} 0.25t^2+1 & -0.5\\-0.25t^2&0.5\end{bmatrix}

因此,它的狀態轉移矩陣爲:\Phi (t,t_0):= \begin{bmatrix} 1 &1 \\ 0.5t^2 & 0.5t^2+2 \end{bmatrix} \begin{bmatrix} 0.25t_0^2+1 &-0.5 \\ -0.25t_0^2 & 0.5 \end{bmatrix}=\begin{bmatrix} 1 & 0\\ 0.5(t^2-t_0^2)& 1 \end{bmatrix}

顯然,轉移矩陣滿足以上(8)(9)(10)三條性質。


 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章