矩阵方程的解法

整理矩阵方程的基本化简、可逆与不可逆情形、拆 E 思想,以及矩阵方程常见易错点与典型例题。

矩阵方程的解法

矩阵方程是线性代数中的重要题型。所谓矩阵方程,是指方程中含有未知矩阵的方程。例如,在方程

AX=BAX=B

中,若 A,BA,B 为已知矩阵,XX 为未知矩阵,则该方程就是一个矩阵方程。

矩阵方程的核心思想是:能够直接解出未知矩阵时,就直接利用逆矩阵求解;若相关矩阵不可逆,则不能强行消去矩阵因子,而应将矩阵方程转化为线性方程组进行求解。


1. 矩阵方程的基本化简

解矩阵方程时,首先应根据题设条件和矩阵运算规则,将方程化为常见形式,例如

AX=B,XA=B,AXB=C.AX=B,\qquad XA=B,\qquad AXB=C.

在化简过程中,常用以下方法。

1.1 消去可逆因子

CA=CB,CA=CB,

CC 可逆,则两端左乘 C1C^{-1},得

C1CA=C1CB,C^{-1}CA=C^{-1}CB,

EA=EB,EA=EB,

所以

A=B.A=B.

同理,若

AC=BC,AC=BC,

CC 可逆,则两端右乘 C1C^{-1},也可得

A=B.A=B.

需要注意的是,只有可逆矩阵才能被消去。若 CC 不可逆,则不能由 CA=CBCA=CB 推出 A=BA=B

1.2 提取公因式

矩阵乘法满足分配律。因此,若出现相同的左因子或右因子,可以提取公因式。

例如,

CA+CB=C(A+B).CA+CB=C(A+B).

又如,

AC+BC=(A+B)C.AC+BC=(A+B)C.

但必须注意,左公因式只能从左侧提出,右公因式只能从右侧提出。由于矩阵乘法一般不满足交换律,所以不能随意改变矩阵因子的位置。

1.3 移项整理

在矩阵方程中,通常把含有未知矩阵的项移到一边,把已知矩阵项移到另一边。例如

AX+C=BAX+C=B

可化为

AX=BC.AX=B-C.

AA 可逆,则

X=A1(BC).X=A^{-1}(B-C).

1.4 常用矩阵公式

解矩阵方程时,常用以下公式。

AAnn 阶方阵,则

AA=AA=AE.AA^*=A^*A=|A|E.

AA 可逆,则

A=AA1.A^*=|A|A^{-1}.

并且

(A)=An2A,n2.(A^*)^*=|A|^{\,n-2}A,\qquad n\ge 2.

对于矩阵乘积的转置,有

(AB)T=BTAT.(AB)^T=B^TA^T.

对于矩阵乘积的逆,有

(AB)1=B1A1.(AB)^{-1}=B^{-1}A^{-1}.

对于伴随矩阵,有

(AB)=BA.(AB)^*=B^*A^*.

此外,若 AA 可逆,则

(A1)=(A)1,(AT)1=(A1)T,(A)T=(AT).(A^{-1})^*=(A^*)^{-1},\qquad (A^T)^{-1}=(A^{-1})^T,\qquad (A^*)^T=(A^T)^*.

这些公式在含有伴随矩阵、转置矩阵、逆矩阵的方程中十分常用。


2. 矩阵方程的基本求解方法

2.1 可逆情形

AX=B,AX=B,

AA 可逆,则

X=A1B.X=A^{-1}B.

XA=B,XA=B,

AA 可逆,则

X=BA1.X=BA^{-1}.

AXB=C,AXB=C,

A,BA,B 均可逆,则

X=A1CB1.X=A^{-1}CB^{-1}.

这里必须注意左右顺序。例如 AXB=CAXB=C 中,A1A^{-1} 应从左侧乘入,B1B^{-1} 应从右侧乘入,不能写成 B1CA1B^{-1}CA^{-1}

2.2 不可逆情形

AA 不可逆,方程

AX=BAX=B

不能写成

X=A1B.X=A^{-1}B.

此时应把 XXBB 按列分块。设

X=(ξ1,ξ2,,ξn),B=(β1,β2,,βn).X=(\xi_1,\xi_2,\cdots,\xi_n),\qquad B=(\beta_1,\beta_2,\cdots,\beta_n).

于是

AX=A(ξ1,ξ2,,ξn)=(Aξ1,Aξ2,,Aξn).AX=A(\xi_1,\xi_2,\cdots,\xi_n)=(A\xi_1,A\xi_2,\cdots,A\xi_n).

AX=BAX=B,得

(Aξ1,Aξ2,,Aξn)=(β1,β2,,βn).(A\xi_1,A\xi_2,\cdots,A\xi_n)=(\beta_1,\beta_2,\cdots,\beta_n).

因此原矩阵方程等价于

Aξi=βi,i=1,2,,n.A\xi_i=\beta_i,\qquad i=1,2,\cdots,n.

也就是说,一个矩阵方程被转化成了若干个同系数矩阵的线性方程组。

若方程不能化为上述几种标准形式,则可设未知矩阵

X=(xij),X=(x_{ij}),

再代入原方程,转化为关于元素 xijx_{ij} 的线性方程组。这就是待定元素法。


解答题例 1

A=(111111111),A= \begin{pmatrix} 1&1&-1\\ -1&1&1\\ 1&-1&1 \end{pmatrix},

且满足

AB(12A)=8A1B+16E.A^*B\left(\frac12 A^*\right)^*=8A^{-1}B+16E.

求矩阵 BB

解答

先计算

A=4.|A|=4.

由于 AA 为三阶矩阵,所以

(A)=AA.(A^*)^*=|A|A.

又因为对于三阶矩阵有

(kA)=k2A,(kA)^*=k^2A^*,

于是

(12A)=(12)2(A)=14AA=A.\left(\frac12 A^*\right)^*=\left(\frac12\right)^2(A^*)^*=\frac14|A|A=A.

又由

A=AA1=4A1,A^*=|A|A^{-1}=4A^{-1},

原方程左边可化为

AB(12A)=4A1BA.A^*B\left(\frac12 A^*\right)^*=4A^{-1}BA.

所以原方程化为

4A1BA=8A1B+16E.4A^{-1}BA=8A^{-1}B+16E.

两边同除以 44,得

A1BA=2A1B+4E.A^{-1}BA=2A^{-1}B+4E.

移项整理,

A1B(A2E)=4E.A^{-1}B(A-2E)=4E.

两边左乘 AA,得

B(A2E)=4A.B(A-2E)=4A.

接下来计算

(A2E)1=12(110011101).(A-2E)^{-1} = -\frac12 \begin{pmatrix} 1&1&0\\ 0&1&1\\ 1&0&1 \end{pmatrix}.

因此

B=4A(A2E)1=(040004400).B=4A(A-2E)^{-1} = \begin{pmatrix} 0&-4&0\\ 0&0&-4\\ -4&0&0 \end{pmatrix}.

本题的关键在于正确使用伴随矩阵公式,尤其是

(A)=An2A.(A^*)^*=|A|^{\,n-2}A.

n=3n=3 时,该式化为

(A)=AA.(A^*)^*=|A|A.
解答题例 2

已知 aa 是常数,且矩阵

A=(12a13027a)A= \begin{pmatrix} 1&2&a\\ 1&3&0\\ 2&7&-a \end{pmatrix}

可经初等列变换化为矩阵

B=(1a2011111).B= \begin{pmatrix} 1&a&2\\ 0&1&1\\ -1&1&1 \end{pmatrix}.

求:

(1) a;(1)\ a; (2) 满足 AP=B 的可逆矩阵 P.(2)\ 满足\ AP=B\ 的可逆矩阵\ P.
解答

由于初等列变换不改变矩阵的秩,所以

r(A)=r(B).r(A)=r(B).

先对 AA 作初等行变换,可得

A(103a01a000).A\sim \begin{pmatrix} 1&0&3a\\ 0&1&-a\\ 0&0&0 \end{pmatrix}.

因此

r(A)=2.r(A)=2.

再对 BB 作初等行变换,可得

B(102a011002a).B\sim \begin{pmatrix} 1&0&2-a\\ 0&1&1\\ 0&0&2-a \end{pmatrix}.

要使

r(B)=2,r(B)=2,

必须有

2a=0.2-a=0.

所以

a=2.a=2.

a=2a=2 时,

A=(122130272),B=(122011111).A= \begin{pmatrix} 1&2&2\\ 1&3&0\\ 2&7&-2 \end{pmatrix}, \qquad B= \begin{pmatrix} 1&2&2\\ 0&1&1\\ -1&1&1 \end{pmatrix}.

因为 AA 不可逆,不能直接写 P=A1BP=A^{-1}B。此时应将 PP 按列分块。设

P=(ξ1,ξ2,ξ3),B=(β1,β2,β3).P=(\xi_1,\xi_2,\xi_3),\qquad B=(\beta_1,\beta_2,\beta_3).

AP=BAP=B

可得

Aξi=βi,i=1,2,3.A\xi_i=\beta_i,\qquad i=1,2,3.

对增广矩阵 (AB)(A\mid B) 作初等行变换,有

(AB)(106344012111000000).(A\mid B) \sim \left( \begin{array}{ccc|ccc} 1&0&6&3&4&4\\ 0&1&-2&-1&-1&-1\\ 0&0&0&0&0&0 \end{array} \right).

于是矩阵方程 AP=BAP=B 等价于

(106012000)P=(344111000).\begin{pmatrix} 1&0&6\\ 0&1&-2\\ 0&0&0 \end{pmatrix}P= \begin{pmatrix} 3&4&4\\ -1&-1&-1\\ 0&0&0 \end{pmatrix}.

A=(106012000).A'= \begin{pmatrix} 1&0&6\\ 0&1&-2\\ 0&0&0 \end{pmatrix}.

则齐次方程组 Ax=0A'x=0 的一个基础解系为

(6,2,1)T.(-6,2,1)^T.

于是三个非齐次方程组的通解分别为

ξ1=k1(6,2,1)T+(3,1,0)T,\xi_1=k_1(-6,2,1)^T+(3,-1,0)^T, ξ2=k2(6,2,1)T+(4,1,0)T,\xi_2=k_2(-6,2,1)^T+(4,-1,0)^T, ξ3=k3(6,2,1)T+(4,1,0)T.\xi_3=k_3(-6,2,1)^T+(4,-1,0)^T.

因此

P=(6k1+36k2+46k3+42k112k212k31k1k2k3),P= \begin{pmatrix} -6k_1+3&-6k_2+4&-6k_3+4\\ 2k_1-1&2k_2-1&2k_3-1\\ k_1&k_2&k_3 \end{pmatrix},

其中 (k1,k2,k3)(k_1,k_2,k_3) 为任意常数。

题目要求 PP 可逆,所以还需满足

P0.|P|\neq 0.

计算得

P=k3k2.|P|=k_3-k_2.

因此

k3k2.k_3\neq k_2.

综上,所求可逆矩阵为

P=(6k1+36k2+46k3+42k112k212k31k1k2k3),k3k2.P= \begin{pmatrix} -6k_1+3&-6k_2+4&-6k_3+4\\ 2k_1-1&2k_2-1&2k_3-1\\ k_1&k_2&k_3 \end{pmatrix}, \qquad k_3\neq k_2.

本题体现了矩阵方程的重要思想:当 AA 不可逆时,不能用 A1A^{-1} 求解,而应按列分块,把矩阵方程转化为线性方程组。

解答题例 3

A=(112237202),B=(101211211).A= \begin{pmatrix} 1&1&2\\ 2&3&7\\ 2&0&-2 \end{pmatrix}, \qquad B= \begin{pmatrix} 1&0&-1\\ 2&1&1\\ 2&1&1 \end{pmatrix}.

求所有满足

PA=BPA=B

的可逆矩阵 PP

解答

本题与上一题的区别在于,未知矩阵 PP 位于左侧。由于 AA 不可逆,不能写成

P=BA1.P=BA^{-1}.

此时可对方程两端同时取转置:

(PA)T=BT.(PA)^T=B^T.

由转置公式

(PA)T=ATPT,(PA)^T=A^TP^T,

ATPT=BT.A^TP^T=B^T.

这就把原方程转化成了

ATX=BT,A^TX=B^T,

其中 X=PTX=P^T

注意:

AT=(122130272),BT=(122011111).A^T= \begin{pmatrix} 1&2&2\\ 1&3&0\\ 2&7&-2 \end{pmatrix}, \qquad B^T= \begin{pmatrix} 1&2&2\\ 0&1&1\\ -1&1&1 \end{pmatrix}.

这与例 2 中 a=2a=2 时的方程完全相同。因此由上一题可得

PT=(6k1+36k2+46k3+42k112k212k31k1k2k3).P^T= \begin{pmatrix} -6k_1+3&-6k_2+4&-6k_3+4\\ 2k_1-1&2k_2-1&2k_3-1\\ k_1&k_2&k_3 \end{pmatrix}.

于是

P=(6k1+32k11k16k2+42k21k26k3+42k31k3).P= \begin{pmatrix} -6k_1+3&2k_1-1&k_1\\ -6k_2+4&2k_2-1&k_2\\ -6k_3+4&2k_3-1&k_3 \end{pmatrix}.

又因为

P=PT,|P|=|P^T|,

而上一题中

PT=k3k2,|P^T|=k_3-k_2,

所以

P=k3k2.|P|=k_3-k_2.

因此 PP 可逆的条件为

k3k2.k_3\neq k_2.

综上,所求可逆矩阵为

P=(6k1+32k11k16k2+42k21k26k3+42k31k3),k3k2.P= \begin{pmatrix} -6k_1+3&2k_1-1&k_1\\ -6k_2+4&2k_2-1&k_2\\ -6k_3+4&2k_3-1&k_3 \end{pmatrix}, \qquad k_3\neq k_2.

本题的关键是:当未知矩阵在左侧时,可以通过两端取转置,把方程化为未知矩阵在右侧的形式。


3. “拆 EE”思想

在矩阵方程中,还有一个非常重要的技巧,叫作“拆 EE”。

所谓“拆 EE”,就是根据题设条件,把单位矩阵 EE 表示成某些矩阵的乘积,从而判断这些矩阵可逆,并进一步求解矩阵方程。

例如,若

A3=O,A^3=O,

EA3=E.E-A^3=E.

EA3=(EA)(E+A+A2).E-A^3=(E-A)(E+A+A^2).

因此

(EA)(E+A+A2)=E.(E-A)(E+A+A^2)=E.

这说明 EAE-A 可逆,且

(EA)1=E+A+A2.(E-A)^{-1}=E+A+A^2.

同理,由于 A3=OA^3=O 可推出 A4=OA^4=O,所以

(EA2)(E+A2)=EA4=E.(E-A^2)(E+A^2)=E-A^4=E.

因此 EA2E-A^2 可逆,且

(EA2)1=E+A2.(E-A^2)^{-1}=E+A^2.

这类方法在含有 A2,A3A^2,A^3 以及 Am=OA^m=O 的矩阵方程中非常常见。

解答题例 4

设矩阵

A=(a101a101a),A= \begin{pmatrix} a&1&0\\ 1&a&-1\\ 0&1&a \end{pmatrix},

A3=O.A^3=O.

求:

(1) a 的值;(1)\ a\ 的值; (2) 若矩阵 X 满足XXA2AX+AXA2=E,(2)\ 若矩阵\ X\ 满足 X-XA^2-AX+AXA^2=E,

其中 EE 为三阶单位矩阵,求 XX

解答

A3=OA^3=O 可知

A3=0.|A^3|=0.

又因为

A3=A3,|A^3|=|A|^3,

所以

A=0.|A|=0.

直接计算

A=a101a101a=a3.|A|= \begin{vmatrix} a&1&0\\ 1&a&-1\\ 0&1&a \end{vmatrix} =a^3.

因此

a3=0,a^3=0,

所以

a=0.a=0.

a=0a=0 时,

A=(010101010).A= \begin{pmatrix} 0&1&0\\ 1&0&-1\\ 0&1&0 \end{pmatrix}.

接下来解矩阵方程

XXA2AX+AXA2=E.X-XA^2-AX+AXA^2=E.

先按左右因子分组:

XXA2=X(EA2),X-XA^2=X(E-A^2), AX+AXA2=AX(EA2).-AX+AXA^2=-AX(E-A^2).

所以原方程化为

(EA)X(EA2)=E.(E-A)X(E-A^2)=E.

由于 A3=OA^3=O,所以

(EA)1=E+A+A2,(E-A)^{-1}=E+A+A^2,

又因为

(EA2)1=E+A2.(E-A^2)^{-1}=E+A^2.

因此

X=(EA)1(EA2)1=(E+A+A2)(E+A2).X=(E-A)^{-1}(E-A^2)^{-1} =(E+A+A^2)(E+A^2).

计算

A2=(101000101).A^2= \begin{pmatrix} 1&0&-1\\ 0&0&0\\ 1&0&-1 \end{pmatrix}.

于是

E+A+A2=(211111110),E+A+A^2= \begin{pmatrix} 2&1&-1\\ 1&1&-1\\ 1&1&0 \end{pmatrix}, E+A2=(201010100).E+A^2= \begin{pmatrix} 2&0&-1\\ 0&1&0\\ 1&0&0 \end{pmatrix}.

相乘得

X=(312111211).X= \begin{pmatrix} 3&1&-2\\ 1&1&-1\\ 2&1&-1 \end{pmatrix}.

所以

X=(312111211).X= \begin{pmatrix} 3&1&-2\\ 1&1&-1\\ 2&1&-1 \end{pmatrix}.

4. 易错点

在例 4 中,很多同学会犯一个典型错误:对方程

XXA2AX+AXA2=EX-XA^2-AX+AXA^2=E

两端同时右乘 AA,得到

XAXA3AXA+AXA3=A.XA-XA^3-AXA+AXA^3=A.

由于 A3=OA^3=O,于是化简为

XAAXA=A.XA-AXA=A.

然后再去求 XX

这种做法是错误的。原因在于:由 A3=OA^3=O 可知 AA 不可逆。若在方程两端同时乘以一个不可逆矩阵,则新方程不一定与原方程等价。

也就是说,原方程的解一定满足右乘 AA 后得到的新方程,但新方程的解未必满足原方程。这样可能会扩大解集,产生多余解。

因此,在矩阵方程中必须牢记:

只有两端乘以可逆矩阵,才不会改变方程的解集。


5. 总结

矩阵方程的解题原则可以概括为以下几句话。

若方程可以化为

AX=B,XA=B,AXB=CAX=B,\qquad XA=B,\qquad AXB=C

且相关矩阵可逆,则直接利用逆矩阵求解。

若矩阵不可逆,则不能强行消去矩阵因子,应将未知矩阵按列分块,把矩阵方程转化为若干个线性方程组。

若未知矩阵在左侧,例如

XA=BXA=B

PA=B,PA=B,

可以考虑两端取转置,把问题转化为未知矩阵在右侧的形式。

若题目给出

Am=OA^m=O

这类条件,应优先想到“拆 EE”思想,即把单位矩阵写成矩阵乘积,从而判断某些矩阵可逆。

矩阵方程的难点不在计算本身,而在于判断每一步变形是否等价。尤其要注意:不可逆矩阵不能随意消去,也不能通过乘以不可逆矩阵来获得等价方程。

Back to archive

Discussion

Comments

Post

Share questions, corrections, or extra notes about this post.