矩陣指數- 維基百科，自由的百科全書 - Wikipedia

矩陣指數（matrix exponential）是方塊矩陣的一種矩陣函數，與指數函數類似。

矩陣指數給出了矩陣 ... X的指數，用eX或exp(X)來表示，是由以下冪級數所給出的n×n矩陣：. 矩陣指數 維基百科，自由的百科全書 跳至導覽 跳至搜尋 矩陣指數（matrixexponential）是方塊矩陣的一種矩陣函數，與指數函數類似。

矩陣指數給出了矩陣李代數與對應的李群之間的關係。

設X為n×n的實數或複數矩陣。

X的指數，用eX或exp(X)來表示，是由以下冪級數所給出的n×n矩陣： e X = ∑ k = 0 ∞ 1 k ! X k {\displaystylee^{X}=\sum_{k=0}^{\infty}{1\overk!}X^{k}} 以上的級數總是收斂的，因此X的指數是定義良好的。

注意，如果X是1×1的矩陣，則X的矩陣指數就是由X的元素的指數所組成的1×1矩陣。

目次 1性質 1.1基本性質 1.2線性微分方程 1.3矩陣指數的行列式 1.4指數相加 1.5李乘積公式 1.6貝克爾-坎貝爾-豪斯多夫公式 1.7指數映射 2矩陣指數的計算 2.1可對角化矩陣 2.2冪零矩陣 2.3推廣 3計算 4應用 4.1線性微分方程 4.1.1例子（齊次） 4.1.2非齊次的情況──參數變換 4.1.3例子（非齊次） 5註釋 6參考文獻 7參閱 8外部連結 性質[編輯] 基本性質[編輯] 設X和Y為n×n的複數矩陣，並設a和b為任意的複數。

我們把n×n的單位矩陣記為I，把零矩陣記為0。

我們可以從指數級數的定義直接得到矩陣指數的如下性質[1]： e0=I exp(XT)=(expX)T，其中XT表示X的轉置。

從中可以推出，如果X是對稱矩陣，則eX也是對稱矩陣；如果X是斜對稱矩陣，則eX是正交矩陣。

exp(X*)=(expX)*，其中X*表示X的共軛轉置。

可以推出，如果X是埃爾米特矩陣，則eX也是埃爾米特矩陣；如果X是斜埃爾米特矩陣，則eX是酉矩陣。

如果Y是可逆矩陣，那麼eYXY−1=YeXY−1 接下來是一個關鍵性質： 如果 X Y = Y X {\displaystyleXY=YX} 那麼 e X e Y = e X + Y {\displaystylee^{X}e^{Y}=e^{X+Y}} 由此導出的推論有： eaXebX=e(a+b)X eXe−X=I 線性微分方程[編輯] 矩陣指數的一個重要性，是它可以用來解微分方程。

從(1)可知，以下微分方程 d d t y ( t ) = A y ( t ) , y ( 0 ) = y 0 {\displaystyle{\frac{d}{dt}}y(t)=Ay(t),\quady(0)=y_{0}} 其中A是矩陣，具有解 y ( t ) = e A t y 0   {\displaystyley(t)=e^{At}y_{0}\} 矩陣指數也可以用來解非齊次方程： d d t y ( t ) = A y ( t ) + z ( t ) , y ( 0 ) = y 0 {\displaystyle{\frac{d}{dt}}y(t)=Ay(t)+z(t),\quady(0)=y_{0}} 參見以下的例子。

當A不是常數時，以下形式的微分方程沒有閉式解： d d t y ( t ) = A ( t ) y ( t ) , y ( 0 ) = y 0 {\displaystyle{\frac{d}{dt}}y(t)=A(t)\,y(t),\quady(0)=y_{0}} 但馬格努斯級數可以給出無窮級數形式的解。

矩陣指數的行列式[編輯] 根據雅可比公式，對任意復矩陣，下列跡等式成立：[2] det ( e A ) = e tr ⁡ ( A )   {\displaystyle\det(e^{A})=e^{\operatorname{tr}(A)}~} 除了提供一種額外的計算工具，這個等式還表明矩陣指數總是可逆矩陣。

這點可以如下證明：因為上述等式的右邊恆不等於0，所以左邊det(eA)≠0，從而eA必可逆。

指數相加[編輯] 我們知道，對於任何實數（純量）x和y，指數函數都滿足公式ex+y=exey。

類似的等式對於可交換矩陣也成立：如果矩陣X和Y是可交換的（即XY=YX），則： e X + Y = e X e Y   {\displaystylee^{X+Y}=e^{X}e^{Y}\} 但是，如果它們不是可交換的，則以上的等式不一定成立。

這個命題反過來不成立：eX+Y=eXeY並不一定就意味著X和Y是可交換的。

但是，如果X和Y只含有代數數，而且它們的大小至少為2×2，則反過來也成立[3]。

X和Y不可交換的情況可以用以下方法計算： 李乘積公式[編輯] 即使 X {\displaystyleX} 和 Y {\displaystyleY} 不可交換， e X + Y {\displaystylee^{X+Y}} 可以用李乘積公式來計算[4] e X + Y = lim n → ∞ ( e X / n e Y / n ) n {\displaystylee^{X+Y}=\lim_{n\rightarrow\infty}(e^{X/n}e^{Y/n})^{n}} 貝克爾-坎貝爾-豪斯多夫公式[編輯] 從另一個方向講，如果 X {\displaystyleX} 和 Y {\displaystyleY} 是元素足夠小（但不一定可交換）的矩陣，我們有： e X e Y = e Z {\displaystylee^{X}e^{Y}=e^{Z}} 其中 Z {\displaystyleZ} 可以通過 X {\displaystyleX} 和 Y {\displaystyleY} 的交換子的級數（貝克爾-坎貝爾-豪斯多夫公式）來計算：[5] Z = X + Y + 1 2 [ X , Y ] + 1 12 [ X , [ X , Y ] ] + ⋯ {\displaystyleZ=X+Y+{\frac{1}{2}}[X,Y]+{\frac{1}{12}}[X,[X,Y]]+\cdots} 其中餘項均為與 X {\displaystyleX} 和 Y {\displaystyleY} 相關的迭代交換子。

指數映射[編輯] 注意矩陣的指數總是非奇異方陣。

eX的逆矩陣由e−X給出。

這與複數的指數總是非零的事實類似。

這樣，矩陣指數就給出了一個映射： exp : M n ( C ) → GL ( n , C ) {\displaystyle\exp\colonM_{n}(\mathbb{C})\to{\mbox{GL}}(n,\mathbb{C})} 這是從所有n×n矩陣的空間到一般線性群（所有非奇異方陣所組成的群）的映射。

實際上，這個映射是滿射，就是說每一個非奇異方陣都可以寫成某個矩陣的指數。

矩陣對數就是這個映射的逆映射。

對於任何兩個矩陣X和Y，我們有： ‖ e X + Y − e X ‖ ≤ ‖ Y ‖ e ‖ X ‖ e ‖ Y ‖ {\displaystyle\|e^{X+Y}-e^{X}\|\leq\|Y\|e^{\|X\|}e^{\|Y\|}} 其中|| · ||表示任意的矩陣範數。

從中可以推出，指數映射在Mn(C)的緊子集內是連續和利普希茨連續的。

以下的映射 t ↦ e t X , t ∈ R {\displaystylet\mapstoe^{tX},\qquadt\in\mathbb{R}} 定義了一般線性群中的一條光滑曲線，當t=0時穿過單位元。

實際上，這給出了一般線性群的一個單參數子群，這是由於： e t X e s X = e ( t + s ) X   {\displaystylee^{tX}e^{sX}=e^{(t+s)X}\} 這條曲線在點t的導數（或切向量）由以下等式給出： d d t e t X = X e t X ( 1 ) {\displaystyle{\frac{d}{dt}}e^{tX}=Xe^{tX}\qquad(1)} t=0時的導數就是矩陣X，所以我們可以說，X是這個單參數子群的推廣。

更加一般地： d d t e X ( t ) = ∫ 0 1 e ( 1 − α ) X ( t ) d X ( t ) d t e α X ( t ) d α {\displaystyle{\frac{d}{dt}}e^{X(t)}=\int_{0}^{1}e^{(1-\alpha)X(t)}{\frac{dX(t)}{dt}}e^{\alphaX(t)}\,d\alpha} 矩陣指數的計算[編輯] 尋找計算矩陣指數的可靠和準確的方法是困難的，目前在數學和數值分析領域中仍然是一個正在研究的話題。

有些方法列舉如下。

可對角化矩陣[編輯] 如果矩陣是對角的： A = [ a 1 0 … 0 0 a 2 … 0 ⋮ ⋮ ⋱ ⋮ 0 0 … a n ] {\displaystyleA={\begin{bmatrix}a_{1}&0&\ldots&0\\0&a_{2}&\ldots&0\\\vdots&\vdots&\ddots&\vdots\\0&0&\ldots&a_{n}\end{bmatrix}}} 則把主對角線上的所有元素取指數，就是原矩陣的指數： e A = [ e a 1 0 … 0 0 e a 2 … 0 ⋮ ⋮ ⋱ ⋮ 0 0 … e a n ] {\displaystylee^{A}={\begin{bmatrix}e^{a_{1}}&0&\ldots&0\\0&e^{a_{2}}&\ldots&0\\\vdots&\vdots&\ddots&\vdots\\0&0&\ldots&e^{a_{n}}\end{bmatrix}}} 這也允許了我們計算可對角化矩陣的指數。

如果 A = U D U − 1 {\displaystyleA=UDU^{-1}} ，且D是對角矩陣，則 e A = U e D U − 1 {\displaystylee^{A}=Ue^{D}U^{-1}} 。

用西爾維斯特公式，也可以得到相同的結果。

冪零矩陣[編輯] 如果對於某個整數q，有Nq=0，則矩陣N稱為冪零矩陣。

在這種情況下，矩陣指數eN可以直接從級數展開式來計算，這是因為級數在有限個項後就終止了： e N = I + N + 1 2 N 2 + 1 6 N 3 + ⋯ + 1 ( q − 1 ) ! N q − 1 {\displaystylee^{N}=I+N+{\frac{1}{2}}N^{2}+{\frac{1}{6}}N^{3}+\cdots+{\frac{1}{(q-1)!}}N^{q-1}} 推廣[編輯] 當矩陣X的最小多項式可以分解為一次多項式的積時，它就可以表示為以下的和： X = A + N   {\displaystyleX=A+N\} 其中： A是可對角化矩陣； N是冪零矩陣； A與N是可交換的（也就是說，AN=NA）。

這稱為Dunford分解。

這就是說，我們可以通過化為前兩種情況，來計算X的指數： e X = e A + N = e A e N   {\displaystylee^{X}=e^{A+N}=e^{A}e^{N}\} 注意為了讓最後一步成立，A和N必須是可交換的。

另外一個密切相關的方法，是利用X的若爾當標準型。

假設X=PJP −1，其中J是X的若爾當標準型。

那麼： e X = P e J P − 1   {\displaystylee^{X}=Pe^{J}P^{-1}\} 另外，由於 J = J a 1 ( λ 1 ) ⊕ J a 2 ( λ 2 ) ⊕ ⋯ ⊕ J a n ( λ n ) {\displaystyleJ=J_{a_{1}}(\lambda_{1})\oplusJ_{a_{2}}(\lambda_{2})\oplus\cdots\oplusJ_{a_{n}}(\lambda_{n})} e J = exp ⁡ ( J a 1 ( λ 1 ) ⊕ J a 2 ( λ 2 ) ⊕ ⋯ ⊕ J a n ( λ n ) ) = exp ⁡ ( J a 1 ( λ 1 ) ) ⊕ exp ⁡ ( J a 2 ( λ 2 ) ) ⊕ ⋯ ⊕ exp ⁡ ( J a k ( λ k ) ) {\displaystyle{\begin{aligned}e^{J}&{}=\exp{\big(}J_{a_{1}}(\lambda_{1})\oplusJ_{a_{2}}(\lambda_{2})\oplus\cdots\oplusJ_{a_{n}}(\lambda_{n}){\big)}\\&{}=\exp{\big(}J_{a_{1}}(\lambda_{1}){\big)}\oplus\exp{\big(}J_{a_{2}}(\lambda_{2}){\big)}\oplus\cdots\oplus\exp{\big(}J_{a_{k}}(\lambda_{k}){\big)}\end{aligned}}} 因此，我們只需要知道怎樣計算若爾當塊的矩陣指數。

但是，每一個若爾當塊都具有形式 J a ( λ ) = λ I + N   {\displaystyleJ_{a}(\lambda)=\lambdaI+N\} 其中N是冪零矩陣。

則這個區塊的矩陣指數由下式給出： e λ I + N = e λ e N   {\displaystylee^{\lambdaI+N}=e^{\lambda}e^{N}\} 計算[編輯] 假設我們想要計算以下矩陣的指數。

B = [ 21 17 6 − 5 − 1 − 6 4 4 16 ] {\displaystyleB={\begin{bmatrix}21&17&6\\-5&-1&-6\\4&4&16\end{bmatrix}}} 它的若爾當型為： J = P − 1 B P = [ 4 0 0 0 16 1 0 0 16 ] {\displaystyleJ=P^{-1}BP={\begin{bmatrix}4&0&0\\0&16&1\\0&0&16\end{bmatrix}}} 其中矩陣P由下式給出： P = [ − 1 4 2 5 4 1 4 − 2 − 1 4 0 4 0 ] {\displaystyleP={\begin{bmatrix}-{\frac{1}{4}}&2&{\frac{5}{4}}\\{\frac{1}{4}}&-2&-{\frac{1}{4}}\\0&4&0\end{bmatrix}}} 我們首先來計算exp(J)。

我們有： J = J 1 ( 4 ) ⊕ J 2 ( 16 ) {\displaystyleJ=J_{1}(4)\oplusJ_{2}(16)} 1×1矩陣的指數僅僅是該矩陣的元素的指數，因此exp(J1(4))=[e4]。

J 2 ( 16 ) {\displaystyleJ_{2}(16)} 的指數可以用以上提到的公式exp(λ I {\displaystyleI} +N)=eλexp(N)來算出： exp ⁡ ( [ 16 1 0 16 ] ) = e 16 exp ⁡ ( [ 0 1 0 0 ] ) = e 16 ( [ 1 0 0 1 ] + [ 0 1 0 0 ] + 1 2 ! [ 0 0 0 0 ] + ⋯ ) = [ e 16 e 16 0 e 16 ] {\displaystyle\exp\left({\begin{bmatrix}16&1\\0&16\end{bmatrix}}\right)=e^{16}\exp\left({\begin{bmatrix}0&1\\0&0\end{bmatrix}}\right)=e^{16}\left({\begin{bmatrix}1&0\\0&1\end{bmatrix}}+{\begin{bmatrix}0&1\\0&0\end{bmatrix}}+{1\over2!}{\begin{bmatrix}0&0\\0&0\end{bmatrix}}+\cdots\right)={\begin{bmatrix}e^{16}&e^{16}\\0&e^{16}\end{bmatrix}}} 因此，原矩陣B的指數為： exp ⁡ ( B ) = P exp ⁡ ( J ) P − 1 = P [ e 4 0 0 0 e 16 e 16 0 0 e 16 ] P − 1 = 1 4 [ 13 e 16 − e 4 13 e 16 − 5 e 4 2 e 16 − 2 e 4 − 9 e 16 + e 4 − 9 e 16 + 5 e 4 − 2 e 16 + 2 e 4 16 e 16 16 e 16 4 e 16 ] {\displaystyle\exp(B)=P\exp(J)P^{-1}=P{\begin{bmatrix}e^{4}&0&0\\0&e^{16}&e^{16}\\0&0&e^{16}\end{bmatrix}}P^{-1}={1\over4}{\begin{bmatrix}13e^{16}-e^{4}&13e^{16}-5e^{4}&2e^{16}-2e^{4}\\-9e^{16}+e^{4}&-9e^{16}+5e^{4}&-2e^{16}+2e^{4}\\16e^{16}&16e^{16}&4e^{16}\end{bmatrix}}} 應用[編輯] 線性微分方程[編輯] 矩陣指數在解線性微分方程時十分有用。

前面曾提到，以下形式的微分方程 y ′ = C y {\displaystyle\mathbf{y}'=C\mathbf{y}} 具有解eCty(0)。

如果我們考慮以下向量 y ( t ) = [ y 1 ( t ) ⋮ y n ( t ) ] {\displaystyle\mathbf{y}(t)={\begin{bmatrix}y_{1}(t)\\\vdots\\y_{n}(t)\end{bmatrix}}} 我們就可以把線性微分方程表示為： y ′ ( t ) = A y ( t ) + b ( t ) {\displaystyle\mathbf{y}'(t)=A\mathbf{y}(t)+\mathbf{b}(t)} 如果我們作一個猜想，把兩邊乘以一個積分因子e−At，便得到： e − A t y ′ − e − A t A y = e − A t b {\displaystylee^{-At}\mathbf{y}'-e^{-At}A\mathbf{y}=e^{-At}\mathbf{b}} d d t ( e − A t y ) = e − A t b {\displaystyle{\frac{d}{dt}}(e^{-At}\mathbf{y})=e^{-At}\mathbf{b}} 如果我們可以計算eAt，那麼就得到了微分方程的解。

例子（齊次）[編輯] 假設我們有以下的微分方程組： x ′ = 2 x − y + z y ′ = 3 y − 1 z z ′ = 2 x + y + 3 z {\displaystyle{\begin{matrix}x'&=&2x&-y&+z\\y'&=&&3y&-1z\\z'&=&2x&+y&+3z\end{matrix}}} 相關的矩陣為： M = [ 2 − 1 1 0 3 − 1 2 1 3 ] {\displaystyleM={\begin{bmatrix}2&-1&1\\0&3&-1\\2&1&3\end{bmatrix}}} 在以上的例子中，我們計算了矩陣指數 e t M = [ 2 e t − 2 t e 2 t − 2 t e 2 t 0 − 2 e t + 2 ( t + 1 ) e 2 t 2 ( t + 1 ) e 2 t 0 2 t e 2 t 2 t e 2 t 2 e t ] {\displaystylee^{tM}={\begin{bmatrix}2e^{t}-2te^{2t}&-2te^{2t}&0\\-2e^{t}+2(t+1)e^{2t}&2(t+1)e^{2t}&0\\2te^{2t}&2te^{2t}&2e^{t}\end{bmatrix}}} 因此微分方程組的通解為： [ x y z ] = C 1 [ 2 e t − 2 t e 2 t − 2 e t + 2 ( t + 1 ) e 2 t 2 t e 2 t ] + C 2 [ − 2 t e 2 t 2 ( t + 1 ) e 2 t 2 t e 2 t ] + C 3 [ 0 0 2 e t ] {\displaystyle{\begin{bmatrix}x\\y\\z\end{bmatrix}}=C_{1}{\begin{bmatrix}2e^{t}-2te^{2t}\\-2e^{t}+2(t+1)e^{2t}\\2te^{2t}\end{bmatrix}}+C_{2}{\begin{bmatrix}-2te^{2t}\\2(t+1)e^{2t}\\2te^{2t}\end{bmatrix}}+C_{3}{\begin{bmatrix}0\\0\\2e^{t}\end{bmatrix}}} 也就是說， x = C 1 ( 2 e t − 2 t e 2 t ) + C 2 ( − 2 t e 2 t ) y = C 1 ( − 2 e t + 2 ( t + 1 ) e 2 t ) + C 2 ( 2 ( t + 1 ) e 2 t ) z = ( C 1 + C 2 ) ( 2 t e 2 t ) + 2 C 3 e t {\displaystyle{\begin{aligned}x&=C_{1}(2e^{t}-2te^{2t})+C_{2}(-2te^{2t})\\y&=C_{1}(-2e^{t}+2(t+1)e^{2t})+C_{2}(2(t+1)e^{2t})\\z&=(C_{1}+C_{2})(2te^{2t})+2C_{3}e^{t}\end{aligned}}} 非齊次的情況──參數變換[編輯] 對於非齊次的情況，我們可以用積分因子的方法（類似於參數變換的方法）。

我們找到形為yp(t)=exp(tA)z(t)一個特解： y p ′ = ( e t A ) ′ z ( t ) + e t A z ′ ( t ) {\displaystyle\mathbf{y}_{p}'=(e^{tA})'\mathbf{z}(t)+e^{tA}\mathbf{z}'(t)} = A e t A z ( t ) + e t A z ′ ( t ) {\displaystyle=Ae^{tA}\mathbf{z}(t)+e^{tA}\mathbf{z}'(t)} = A y p ( t ) + e t A z ′ ( t ) {\displaystyle=A\mathbf{y}_{p}(t)+e^{tA}\mathbf{z}'(t)} 為了讓yp為方程的解，必須有： e t A z ′ ( t ) = b ( t ) {\displaystylee^{tA}\mathbf{z}'(t)=\mathbf{b}(t)} z ′ ( t ) = ( e t A ) − 1 b ( t ) {\displaystyle\mathbf{z}'(t)=(e^{tA})^{-1}\mathbf{b}(t)} z ( t ) = ∫ 0 t e − u A b ( u ) d u + c {\displaystyle\mathbf{z}(t)=\int_{0}^{t}e^{-uA}\mathbf{b}(u)\,du+\mathbf{c}} 因此， y p = e t A ∫ 0 t e − u A b ( u ) d u + e t A c = ∫ 0 t e ( t − u ) A b ( u ) d u + e t A c {\displaystyle{\begin{aligned}\mathbf{y}_{p}&{}=e^{tA}\int_{0}^{t}e^{-uA}\mathbf{b}(u)\,du+e^{tA}\mathbf{c}\\&{}=\int_{0}^{t}e^{(t-u)A}\mathbf{b}(u)\,du+e^{tA}\mathbf{c}\end{aligned}}} 其中c由問題的初始條件決定。

例子（非齊次）[編輯] 假設我們有以下的微分方程組： x ′ = 2 x − y + z + e 2 t y ′ = 3 y − 1 z z ′ = 2 x + y + 3 z + e 2 t {\displaystyle{\begin{matrix}x'&=&2x&-y&+z&+e^{2t}\\y'&=&&3y&-1z&\\z'&=&2x&+y&+3z&+e^{2t}\end{matrix}}} 那麼我們有 M = [ 2 − 1 1 0 3 − 1 2 1 3 ] {\displaystyleM={\begin{bmatrix}2&-1&1\\0&3&-1\\2&1&3\end{bmatrix}}} 以及 b = e 2 t [ 1 0 1 ] {\displaystyle\mathbf{b}=e^{2t}{\begin{bmatrix}1\\0\\1\end{bmatrix}}} 用前面的方法，我們可以得出齊次微分方程的解。

由於齊次方程的通解與非齊次方程的特解的和就是非齊次方程的通解，因此我們只需要找到一個特解（用參數變換法）。

我們有： y p = e t ∫ 0 t e ( − u ) M [ e 2 u 0 e 2 u ] d u + e t M c {\displaystyle\mathbf{y}_{p}=e^{t}\int_{0}^{t}e^{(-u)M}{\begin{bmatrix}e^{2u}\\0\\e^{2u}\end{bmatrix}}\,du+e^{tM}\mathbf{c}} y p = e t ∫ 0 t [ 2 e u − 2 u e 2 u − 2 u e 2 u 0 − 2 e u + 2 ( u + 1 ) e 2 u 2 ( u + 1 ) e 2 u 0 2 u e 2 u 2 u e 2 u 2 e u ] [ e 2 u 0 e 2 u ] d u + e t M c {\displaystyle\mathbf{y}_{p}=e^{t}\int_{0}^{t}{\begin{bmatrix}2e^{u}-2ue^{2u}&-2ue^{2u}&0\\\\-2e^{u}+2(u+1)e^{2u}&2(u+1)e^{2u}&0\\\\2ue^{2u}&2ue^{2u}&2e^{u}\end{bmatrix}}{\begin{bmatrix}e^{2u}\\0\\e^{2u}\end{bmatrix}}\,du+e^{tM}\mathbf{c}} y p = e t ∫ 0 t [ e 2 u ( 2 e u − 2 u e 2 u ) e 2 u ( − 2 e u + 2 ( 1 + u ) e 2 u ) 2 e 3 u + 2 u e 4 u ] + e t M c {\displaystyle\mathbf{y}_{p}=e^{t}\int_{0}^{t}{\begin{bmatrix}e^{2u}(2e^{u}-2ue^{2u})\\\\e^{2u}(-2e^{u}+2(1+u)e^{2u})\\\\2e^{3u}+2ue^{4u}\end{bmatrix}}+e^{tM}\mathbf{c}} y p = e t [ − 1 24 e 3 t ( 3 e t ( 4 t − 1 ) − 16 ) 1 24 e 3 t ( 3 e t ( 4 t + 4 ) − 16 ) 1 24 e 3 t ( 3 e t ( 4 t − 1 ) − 16 ) ] + [ 2 e t − 2 t e 2 t − 2 t e 2 t 0 − 2 e t + 2 ( t + 1 ) e 2 t 2 ( t + 1 ) e 2 t 0 2 t e 2 t 2 t e 2 t 2 e t ] [ c 1 c 2 c 3 ] {\displaystyle\mathbf{y}_{p}=e^{t}{\begin{bmatrix}-{1\over24}e^{3t}(3e^{t}(4t-1)-16)\\\\{1\over24}e^{3t}(3e^{t}(4t+4)-16)\\\\{1\over24}e^{3t}(3e^{t}(4t-1)-16)\end{bmatrix}}+{\begin{bmatrix}2e^{t}-2te^{2t}&-2te^{2t}&0\\\\-2e^{t}+2(t+1)e^{2t}&2(t+1)e^{2t}&0\\\\2te^{2t}&2te^{2t}&2e^{t}\end{bmatrix}}{\begin{bmatrix}c_{1}\\c_{2}\\c_{3}\end{bmatrix}}} 進一步簡化，就可以得到原方程的特解。

註釋[編輯] ^Hall2015harvnberror:notarget:CITEREFHall2015(help)Proposition2.3 ^Hall2015harvnberror:notarget:CITEREFHall2015(help)Theorem2.12 ^（Horn&Johnson1991，pp.435–437） ^Hall2015harvnberror:notarget:CITEREFHall2015(help)Theorem2.11 ^Hall2015harvnberror:notarget:CITEREFHall2015(help)Chapter5 參考文獻[編輯] Horn,RogerA.;Johnson,CharlesR.,TopicsinMatrixAnalysis,CambridgeUniversityPress,1991,ISBN 978-0-521-46713-1 . Moler,Cleve;VanLoan,CharlesF.,NineteenDubiousWaystoComputetheExponentialofaMatrix,Twenty-FiveYearsLater(PDF),SIAMReview,2003,45(1):3–49[2008-08-14],ISSN 1095-7200,（原始內容存檔(PDF)於2008-12-08） . 參閱[編輯] 矩陣對數 指數函數 指數映射 向量流 高登─湯普森不等式 位相型分布 外部連結[編輯] 埃里克·韋斯坦因.MatrixExponential.MathWorld.  矩陣指數的教程 取自「https://zh.wikipedia.org/w/index.php?title=矩阵指数&oldid=67860958」 分類：​矩陣論李群指數隱藏分類：​含有哈佛參考文獻格式系列模板連結指向錯誤的頁面 導覽選單 個人工具 沒有登入討論貢獻建立帳號登入 命名空間 條目討論 臺灣正體 不转换简体繁體大陆简体香港繁體澳門繁體大马简体新加坡简体臺灣正體 查看 閱讀編輯檢視歷史 更多 搜尋 導航 首頁分類索引特色內容新聞動態近期變更隨機條目資助維基百科 說明 說明維基社群方針與指引互助客棧知識問答字詞轉換IRC即時聊天聯絡我們關於維基百科 工具 連結至此的頁面相關變更上傳檔案特殊頁面靜態連結頁面資訊引用此頁面維基數據項目 列印/匯出 下載為PDF可列印版 其他語言 CatalàDeutschEnglishEspañolFrançaisItaliano日本語한국어PolskiPortuguêsРусскийSlovenščinaSvenskaУкраїнськаاردو 編輯連結