文章目錄
- 十一、方程組解的結構和性質
- 1、齊次線性方程組
- 2、非齊次線性方程組
- 十二、Ax=0的基礎解系
- 十三、兩個方程組的公共解
- 十四、同解方程
- 十五、求特征值、特征向量
- 十六、判斷A能否相似對角化
- 十七、若A可以相似對角化,求P(Q)
- 十八、二次型化標準型
- 1、拉格朗日配方法
- 2、正交化法
- 十九、二次型正定
- 二十、等價、相似、合同
十一、方程組解的結構和性質
1、齊次線性方程組
方程組
{a11x1+a12x2+?+a1nxn=0,a21x1+a22x2+?+a2nxn=0,??am1x1+am2x2+?+amnxn=0\begin{cases} a_{11}x_1 + a_{12}x_2 + \cdots + a_{1n}x_n = 0, \\ a_{21}x_1 + a_{22}x_2 + \cdots + a_{2n}x_n = 0, \\ \quad \quad \quad \quad \cdots \cdots \\ a_{m1}x_1 + a_{m2}x_2 + \cdots + a_{mn}x_n = 0 \end{cases} ????a11?x1?+a12?x2?+?+a1n?xn?=0,a21?x1?+a22?x2?+?+a2n?xn?=0,??am1?x1?+am2?x2?+?+amn?xn?=0? (Ⅰ)
稱為m個方程,n個未知量的齊次線性方程組
(1)有解的條件
①當r(A) = n(α1,α2,???,αnα_1,α_2,···,α_nα1?,α2?,???,αn?線性無關)時,方程組(Ⅰ)有唯一零解
②當r(A) = r < n(α1,α2,???,αnα_1,α_2,···,α_nα1?,α2?,???,αn?線性相關)時,方程組(Ⅰ)有非零解(無窮多解),且有n-r個線性無關解
(2)求解方法
①將系數矩陣A作初等行變換化為行階梯形矩陣B,求出r(A)
②按列找出一個秩為r的子矩陣,剩余列位置的未知數設為自由變量 n - r(A)個自由變量
③ 按基礎解系定義求出 ξ1,ξ2,…,ξn?r\boldsymbol{\xi}_1, \boldsymbol{\xi}_2, \dots, \boldsymbol{\xi}_{n-r}ξ1?,ξ2?,…,ξn?r? ,并寫出通解。
2、非齊次線性方程組
{a11x1+a12x2+?+a1nxn=b1,a21x1+a22x2+?+a2nxn=b2,??am1x1+am2x2+?+amnxn=bm\begin{cases} a_{11}x_1 + a_{12}x_2 + \cdots + a_{1n}x_n = b_1, \\ a_{21}x_1 + a_{22}x_2 + \cdots + a_{2n}x_n = b_2, \\ \quad \quad \quad \quad \cdots \cdots \\ a_{m1}x_1 + a_{m2}x_2 + \cdots + a_{mn}x_n = b_m \end{cases} ????a11?x1?+a12?x2?+?+a1n?xn?=b1?,a21?x1?+a22?x2?+?+a2n?xn?=b2?,??am1?x1?+am2?x2?+?+amn?xn?=bm?? (Ⅱ)
稱為m個方程,n個未知量的非齊次線性方程組
(1)有解的條件
①若r(A)≠r([A,b])(b不能由α1,α2,???,αnα_1,α_2,···,α_nα1?,α2?,???,αn?線性表示),則方程組(Ⅱ)無解
②若r(A)=r([A,b]) = n(即α1,α2,???,αnα_1,α_2,···,α_nα1?,α2?,???,αn?線性無關,α1,α2,???,αn,bα_1,α_2,···,α_n,bα1?,α2?,???,αn?,b線性相關),則方程組(Ⅱ)有唯一解
③若r(A)=r([A,b]) < n,則方程組(Ⅱ)有無窮多解
(2)求解方法
① 寫出 Ax=bA\boldsymbol{x} = \boldsymbol{b}Ax=b 的導出方程組 Ax=0A\boldsymbol{x} = \boldsymbol{0}Ax=0 ,并求 Ax=0A\boldsymbol{x} = \boldsymbol{0}Ax=0 的通解:
x=k1ξ1+k2ξ2+?+kn?rξn?r\boldsymbol{x} = k_1\boldsymbol{\xi}_1 + k_2\boldsymbol{\xi}_2 + \dots + k_{n-r}\boldsymbol{\xi}_{n-r} x=k1?ξ1?+k2?ξ2?+?+kn?r?ξn?r?
② 求 Ax=bA\boldsymbol{x} = \boldsymbol{b}Ax=b 的一個特解 η\boldsymbol{\eta}η
③非齊次線性方程組 Ax=bA\boldsymbol{x} = \boldsymbol{b}Ax=b 的通解為:
x=k1ξ1+k2ξ2+?+kn?rξn?r+η\boldsymbol{x} = k_1\boldsymbol{\xi}_1 + k_2\boldsymbol{\xi}_2 + \dots + k_{n-r}\boldsymbol{\xi}_{n-r} + \boldsymbol{\eta} x=k1?ξ1?+k2?ξ2?+?+kn?r?ξn?r?+η
其中 k1,k2,…,kn?rk_1, k_2, \dots, k_{n-r}k1?,k2?,…,kn?r?為任意常數
十二、Ax=0的基礎解系
設 ξ1,ξ2,…,ξn?r\boldsymbol{\xi}_1, \boldsymbol{\xi}_2, \dots, \boldsymbol{\xi}_{n-r}ξ1?,ξ2?,…,ξn?r? 滿足以下充要條件,則稱 ξ1,ξ2,…,ξn?r\boldsymbol{\xi}_1, \boldsymbol{\xi}_2, \dots, \boldsymbol{\xi}_{n-r}ξ1?,ξ2?,…,ξn?r? 為齊次線性方程組 Ax=0A\boldsymbol{x} = \boldsymbol{0}Ax=0 的基礎解系:
-
是方程組的解:
ξ1,ξ2,…,ξn?r\boldsymbol{\xi}_1, \boldsymbol{\xi}_2, \dots, \boldsymbol{\xi}_{n-r}ξ1?,ξ2?,…,ξn?r? 均滿足 Ax=0A\boldsymbol{x} = \boldsymbol{0}Ax=0(即屬于解空間);
-
線性無關:
ξ1,ξ2,…,ξn?r\boldsymbol{\xi}_1, \boldsymbol{\xi}_2, \dots, \boldsymbol{\xi}_{n-r}ξ1?,ξ2?,…,ξn?r? 作為向量組線性無關(是解空間的一組“基”的候選);
-
能表示所有解:
方程組 Ax=0A\boldsymbol{x} = \boldsymbol{0}Ax=0的任一解都可由 ξ1,ξ2,…,ξn?r\boldsymbol{\xi}_1, \boldsymbol{\xi}_2, \dots, \boldsymbol{\xi}_{n-r}ξ1?,ξ2?,…,ξn?r? 線性表示(即它們構成解空間的一組基)。
十三、兩個方程組的公共解
已知線性方程組:
{(I){x1+x2=0,x2?x4=0(II){x1?x2+x3=0,x2?x3+x4=0\begin{cases} \text{(I)} & \begin{cases} x_1 + x_2 = 0, \\ x_2 - x_4 = 0 \end{cases} \\[1em] \text{(II)} & \begin{cases} x_1 - x_2 + x_3 = 0, \\ x_2 - x_3 + x_4 = 0 \end{cases} \end{cases} ????(I)(II)?{x1?+x2?=0,x2??x4?=0?{x1??x2?+x3?=0,x2??x3?+x4?=0??
(1) 求方程組 (I)、(II) 的基礎解系
(2) 求方程組 (I)、(II) 的公共解
【解】
(1)
方程組 (I)的基礎解析為
ξ1=(0010),ξ2=(?1101)\boldsymbol{\xi}_1 = \begin{pmatrix} 0 \\ 0 \\ 1 \\ 0 \end{pmatrix}, \quad \boldsymbol{\xi}_2 = \begin{pmatrix} -1 \\ 1 \\ 0 \\ 1 \end{pmatrix} ξ1?=?0010??,ξ2?=??1101??
方程組(II)的基礎解析為
η1=(0110),η2=(?1?101)\boldsymbol{\eta}_1 = \begin{pmatrix} 0 \\ 1 \\ 1 \\ 0 \end{pmatrix}, \quad \boldsymbol{\eta}_2 = \begin{pmatrix} -1 \\ -1 \\ 0 \\ 1 \end{pmatrix} η1?=?0110??,η2?=??1?101??
(2)
方法一:聯立方程
聯立后的系數矩陣為:
[AB]=[1100010?11?11001?11]\begin{bmatrix} \boldsymbol{A} \\ \boldsymbol{B} \end{bmatrix} = \begin{bmatrix} 1 & 1 & 0 & 0 \\ 0 & 1 & 0 & -1 \\ 1 & -1 & 1 & 0 \\ 0 & 1 & -1 & 1 \end{bmatrix}[AB?]=?1010?11?11?001?1?0?101??
對矩陣作初等行變換
[1001010?1001?20000]\begin{bmatrix} 1 & 0 & 0 & 1 \\ 0 & 1 & 0 & -1 \\ 0 & 0 & 1 & -2 \\ 0 & 0 & 0 & 0 \end{bmatrix}?1000?0100?0010?1?1?20??
因此方程組 (I)、(II)的公共解為
x=k[?1121],k∈R\boldsymbol{x} = k \begin{bmatrix} -1 \\ 1 \\ 2 \\ 1 \end{bmatrix}, \quad k \in \mathbb{R}x=k??1121??,k∈R
方法二:通解代入
在方程組 (I) 的通解中,篩選出同時滿足方程組 (II) 的解,即為 (I)、(II) 的公共解(也可在 (II) 的通解中篩選滿足 (I) 的解 )
已知方程組 (I) 的基礎解系為 ξ1,ξ2\boldsymbol{\xi}_1, \boldsymbol{\xi}_2ξ1?,ξ2? ,因此其通解為:
x=k1ξ1+k2ξ2=k1[0010]+k2[?1101]=[?k2k2k1k2]\boldsymbol{x} = k_1\boldsymbol{\xi}_1 + k_2\boldsymbol{\xi}_2 = k_1 \begin{bmatrix} 0 \\ 0 \\ 1 \\ 0 \end{bmatrix} + k_2 \begin{bmatrix} -1 \\ 1 \\ 0 \\ 1 \end{bmatrix} = \begin{bmatrix} -k_2 \\ k_2 \\ k_1 \\ k_2 \end{bmatrix}x=k1?ξ1?+k2?ξ2?=k1??0010??+k2???1101??=??k2?k2?k1?k2???
(其中 k1,k2∈Rk_1, k_2 \in \mathbb{R}k1?,k2?∈R 為任意常數 )
將通解 x=[?k2k2k1k2]代入方程組(II)\boldsymbol{x} = \begin{bmatrix} -k_2 \\ k_2 \\ k_1 \\ k_2 \end{bmatrix} 代入方程組 (II)x=??k2?k2?k1?k2???代入方程組(II):
{x1?x2+x3=0,x2?x3+x4=0\begin{cases}x_1 - x_2 + x_3 = 0, \\ x_2 - x_3 + x_4 = 0 \end{cases}{x1??x2?+x3?=0,x2??x3?+x4?=0?
代入第1個方程:
(?k2)?k2+k1=0?k1?2k2=0?k1=2k2(-k_2) - k_2 + k_1 = 0 \implies k_1 - 2k_2 = 0 \implies k_1 = 2k_2(?k2?)?k2?+k1?=0?k1??2k2?=0?k1?=2k2?
代入第2個方程:
k2?k1+k2=0?2k2?k1=0k_2 - k_1 + k_2 = 0 \implies 2k_2 - k_1 = 0k2??k1?+k2?=0?2k2??k1?=0
結合 $k_1 = 2k_2 $,得方程組(I)、(II) 的公共解為
x=k2[?1121],k2∈R\boldsymbol{x} = k_2 \begin{bmatrix} -1 \\ 1 \\ 2 \\ 1 \end{bmatrix}, \quad k_2 \in \mathbb{R}x=k2???1121??,k2?∈R
方法三:通解相等
(I)的通解:k1ξ1+k2ξ2=k1[0010]+k2[?1101]=[?k2k2k1k2](I) 的通解:k_1\boldsymbol{\xi}_1 + k_2\boldsymbol{\xi}_2 = k_1 \begin{bmatrix} 0 \\ 0 \\ 1 \\ 0 \end{bmatrix} + k_2 \begin{bmatrix} -1 \\ 1 \\ 0 \\ 1 \end{bmatrix} = \begin{bmatrix} -k_2 \\ k_2 \\ k_1 \\ k_2 \end{bmatrix}(I)的通解:k1?ξ1?+k2?ξ2?=k1??0010??+k2???1101??=??k2?k2?k1?k2???
(II)的通解:l1η1+l2η2=l1[0110]+l2[?1?101]=[?l2l1?l2l1l2](II) 的通解:l_1\boldsymbol{\eta}_1 + l_2\boldsymbol{\eta}_2 = l_1 \begin{bmatrix} 0 \\ 1 \\ 1 \\ 0 \end{bmatrix} + l_2 \begin{bmatrix} -1 \\ -1 \\ 0 \\ 1 \end{bmatrix} = \begin{bmatrix} -l_2 \\ l_1 - l_2 \\ l_1 \\ l_2 \end{bmatrix}(II)的通解:l1?η1?+l2?η2?=l1??0110??+l2???1?101??=??l2?l1??l2?l1?l2???
由上式可得k2=l2,k2=l1?l2,k1=l1k_2 = l_2, k_2 = l_1 - l_2, k_1 = l_1k2?=l2?,k2?=l1??l2?,k1?=l1?
故
k1=2k2k_1 = 2k_2k1?=2k2? 或 l1=2l2l_1 = 2l_2l1?=2l2?
因此公共解為
x=k2[?1121],k2∈R\boldsymbol{x} = k_2 \begin{bmatrix} -1 \\ 1 \\ 2 \\ 1 \end{bmatrix}, \quad k_2 \in \mathbb{R}x=k2???1121??,k2?∈R
或
x=l2[?1121],l2∈R\boldsymbol{x} = l_2 \begin{bmatrix} -1 \\ 1 \\ 2 \\ 1 \end{bmatrix}, \quad l_2 \in \mathbb{R}x=l2???1121??,l2?∈R
十四、同解方程
Ax=0A\boldsymbol{x} = \boldsymbol{0}Ax=0 與 Bx=0B\boldsymbol{x} = \boldsymbol{0}Bx=0 同解
<=> 基礎解系為等價向量組
<=>A、BA、BA、B行向量組為等價向量組
<=> Ax=0A\boldsymbol{x} = \boldsymbol{0}Ax=0的解均為Bx=0B\boldsymbol{x} = \boldsymbol{0}Bx=0的解且 r(A)=r(B)r(\boldsymbol{A}) = r(\boldsymbol{B})r(A)=r(B)
<=>r(A)=r(B)=r(AB)r(\boldsymbol{A}) = r(\boldsymbol{B}) = r\begin{pmatrix} \boldsymbol{A} \\ \boldsymbol{B} \end{pmatrix}r(A)=r(B)=r(AB?)
十五、求特征值、特征向量
方法一:|λE-A|=0,求λ,回代(λiE?A)x=0(λ_iE-A)x=0(λi?E?A)x=0求α
方法二:常用結論
-
行列式與跡(對 n 階矩陣 A,λ1,λ2,…,λn\boldsymbol{A} , \lambda_1, \lambda_2, \dots, \lambda_nA,λ1?,λ2?,…,λn?為特征值 )
- ∣A∣=λ1λ2?λn|\boldsymbol{A}| = \lambda_1 \lambda_2 \cdots \lambda_n∣A∣=λ1?λ2??λn?(行列式等于特征值之積 )
- tr(A)=λ1+λ2+?+λn\text{tr}(\boldsymbol{A}) = \lambda_1 + \lambda_2 + \cdots + \lambda_ntr(A)=λ1?+λ2?+?+λn?(跡等于特征值之和 )
-
多項式矩陣的特征值(若 Aα=λα\boldsymbol{A}\boldsymbol{\alpha} = \lambda\boldsymbol{\alpha}Aα=λα,則 )
對多項式 f(x)f(x)f(x),有:
f(A)α=f(λ)αf(\boldsymbol{A})\boldsymbol{\alpha} = f(\lambda)\boldsymbol{\alpha}f(A)α=f(λ)α
具體應用:
- Akα=λkα\boldsymbol{A}^k\boldsymbol{\alpha} = \lambda^k\boldsymbol{\alpha}Akα=λkα(kkk 次冪 )
- (A+kE)α=(λ+k)α(\boldsymbol{A} + k\boldsymbol{E})\boldsymbol{\alpha} = (\lambda + k)\boldsymbol{\alpha}(A+kE)α=(λ+k)α(加數量矩陣 )
- 若 A\boldsymbol{A}A 可逆,則 A?1α=1λα\boldsymbol{A}^{-1}\boldsymbol{\alpha} = \frac{1}{\lambda}\boldsymbol{\alpha}A?1α=λ1?α,A?α=∣A∣λα\boldsymbol{A}^*\boldsymbol{\alpha} = \frac{|\boldsymbol{A}|}{\lambda}\boldsymbol{\alpha}A?α=λ∣A∣?α(伴隨矩陣 )
- 相似變換:P?1APα=λα\boldsymbol{P}^{-1}\boldsymbol{A}\boldsymbol{P}\boldsymbol{\alpha} = \lambda\boldsymbol{\alpha}P?1APα=λα(相似矩陣特征值相同,特征向量變換為 P?1α\boldsymbol{P}^{-1}\boldsymbol{\alpha}P?1α )
-
特殊特征值
- 若 A\boldsymbol{A}A 為對合矩陣(A2=E\boldsymbol{A}^2 = \boldsymbol{E}A2=E ),則 λ=±1\lambda = \pm 1λ=±1
- 若 A\boldsymbol{A}A 行和為 $ a $,則 λ=a\lambda = aλ=a 是一個特征值,對應特征向量 α=(11?1)\boldsymbol{\alpha} = \begin{pmatrix} 1 \\ 1 \\ \vdots \\ 1 \end{pmatrix}α=?11?1??(所有分量為1 )
-
特征值的重數
若 AB=λB\boldsymbol{A}\boldsymbol{B} = \lambda\boldsymbol{B}AB=λB 且 B≠0\boldsymbol{B} \neq \boldsymbol{0}B=0,則 λ\lambdaλ 是 A\boldsymbol{A}A 的特征值,且 B\boldsymbol{B}B 的非零列是對應特征向量;若 B\boldsymbol{B}B 有 nnn 個線性無關列滿足,則 λ\lambdaλ 至少是 nnn 重特征值
-
二次型與特征值
- 二次型 f=xTAxf = \boldsymbol{x}^\text{T}\boldsymbol{A}\boldsymbol{x}f=xTAx 經正交變換 x=Qy\boldsymbol{x} = \boldsymbol{Q}\boldsymbol{y}x=Qy 化為標準型 λ1y12+λ2y22+?+λnyn2\lambda_1 y_1^2 + \lambda_2 y_2^2 + \cdots + \lambda_n y_n^2λ1?y12?+λ2?y22?+?+λn?yn2?,其中 λi\lambda_iλi? 是 A\boldsymbol{A}A 的特征值
-
相似矩陣的特征值
若 A~B\boldsymbol{A} \sim \boldsymbol{B}A~B(相似 ),則 A\boldsymbol{A}A 與 B\boldsymbol{B}B 特征值完全相同(包括重數 ),但特征向量不同(滿足 Aα=λα?B(P?1α)=λ(P?1α)\boldsymbol{A}\boldsymbol{\alpha} = \lambda\boldsymbol{\alpha} \iff \boldsymbol{B}(\boldsymbol{P}^{-1}\boldsymbol{\alpha}) = \lambda(\boldsymbol{P}^{-1}\boldsymbol{\alpha})Aα=λα?B(P?1α)=λ(P?1α),P\boldsymbol{P}P 為相似變換矩陣 )
十六、判斷A能否相似對角化
方法一:基于特征值和特征向量的個數判斷(適用于一般矩陣)
- 判斷條件:n階矩陣A\boldsymbol{A}A可相似對角化的充分必要條件是A\boldsymbol{A}A有n個線性無關的特征向量。
- 具體步驟
- 計算特征值:根據特征方程∣λE?A∣=0\vert\lambda\boldsymbol{E} - \boldsymbol{A}\vert = 0∣λE?A∣=0 ,求出矩陣A\boldsymbol{A}A的所有特征值λ1,λ2,?,λs\lambda_1,\lambda_2,\cdots,\lambda_sλ1?,λ2?,?,λs?,以及它們對應的代數重數n1,n2,?,nsn_1,n_2,\cdots,n_sn1?,n2?,?,ns?(特征值λi\lambda_iλi?的代數重數是指它在特征方程的根中出現的重數,且n1+n2+?+ns=nn_1 + n_2+\cdots + n_s = nn1?+n2?+?+ns?=n)。
- 計算特征向量并判斷線性無關性:對于每個特征值λi\lambda_iλi?,求解齊次線性方程組(λiE?A)x=0(\lambda_i\boldsymbol{E} - \boldsymbol{A})\boldsymbol{x} = \boldsymbol{0}(λi?E?A)x=0,得到其基礎解系,基礎解系中的向量就是屬于λi\lambda_iλi?的線性無關的特征向量,設其個數為mim_imi?,mim_imi?也被稱為特征值λi\lambda_iλi?的幾何重數, 即屬于λi\lambda_iλi?的線性無關特征向量的個數)。若對于每一個特征值λi\lambda_iλi?,都有其代數重數nin_ini?等于幾何重數mim_imi?,即ni=min_i = m_ini?=mi?,i=1,2,?,si = 1,2,\cdots,si=1,2,?,s,則矩陣A\boldsymbol{A}A有n個線性無關的特征向量,A\boldsymbol{A}A可以相似對角化;若存在某個特征值,其代數重數不等于幾何重數,則A\boldsymbol{A}A不能相似對角化。
方法二:判斷矩陣是否為實對稱矩陣(適用于實矩陣)
- 判斷條件:實對稱矩陣一定可以相似對角化,并且可以正交相似對角化(即存在正交矩陣Q\boldsymbol{Q}Q,使得Q?1AQ=QTAQ\boldsymbol{Q}^{-1}\boldsymbol{A}\boldsymbol{Q}=\boldsymbol{Q}^T\boldsymbol{A}\boldsymbol{Q}Q?1AQ=QTAQ為對角矩陣)。
- 具體步驟:只需判斷矩陣A\boldsymbol{A}A是否滿足AT=A\boldsymbol{A}^T = \boldsymbol{A}AT=A,若滿足,則A\boldsymbol{A}A可相似對角化。
十七、若A可以相似對角化,求P(Q)
若矩陣 A\boldsymbol{A}A 可相似對角化,求可逆矩陣 P\boldsymbol{P}P(或正交矩陣 Q\boldsymbol{Q}Q)的步驟
一、求可逆矩陣 P\boldsymbol{P}P 使 P?1AP=Λ\boldsymbol{P}^{-1}\boldsymbol{A}\boldsymbol{P} = \boldsymbol{\Lambda}P?1AP=Λ(Λ\boldsymbol{\Lambda}Λ 為對角矩陣)
- 求特征值 解特征方程 ∣λE?A∣=0|\lambda\boldsymbol{E} - \boldsymbol{A}| = 0∣λE?A∣=0,得所有特征值 λ1,λ2,…,λn\lambda_1, \lambda_2, \dots, \lambda_nλ1?,λ2?,…,λn?(含重數)。
- 求特征向量 對每個特征值 λi\lambda_iλi?,解齊次方程組 (λiE?A)x=0(\lambda_i\boldsymbol{E} - \boldsymbol{A})\boldsymbol{x} = \boldsymbol{0}(λi?E?A)x=0,得基礎解系 ξi1,ξi2,…,ξiki\boldsymbol{\xi}_{i1}, \boldsymbol{\xi}_{i2}, \dots, \boldsymbol{\xi}_{ik_i}ξi1?,ξi2?,…,ξiki??(kik_iki? 為幾何重數,且 ∑ki=n\sum k_i = n∑ki?=n)。
- 構造矩陣 P\boldsymbol{P}P 與對角矩陣 Λ\boldsymbol{\Lambda}Λ
-
將所有線性無關的特征向量按列排列,組成可逆矩陣: P=(ξ11,ξ12,…,ξ1k1,ξ21,…,ξnkn)\boldsymbol{P} = (\boldsymbol{\xi}_{11}, \boldsymbol{\xi}_{12}, \dots, \boldsymbol{\xi}_{1k_1}, \boldsymbol{\xi}_{21}, \dots, \boldsymbol{\xi}_{nk_n})P=(ξ11?,ξ12?,…,ξ1k1??,ξ21?,…,ξnkn??)
-
對角矩陣 Λ\boldsymbol{\Lambda}Λ 的對角線元素為對應特征值,順序與 P\boldsymbol{P}P 的列向量一致: Λ=(λ1λ2?λn)\boldsymbol{\Lambda} = \begin{pmatrix} \lambda_1 & & & \\ & \lambda_2 & & \\ & & \ddots & \\ & & & \lambda_n \end{pmatrix}Λ=?λ1??λ2????λn???
二、求正交矩陣 Q\boldsymbol{Q}Q 使 Q?1AQ=QTAQ=Λ\boldsymbol{Q}^{-1}\boldsymbol{A}\boldsymbol{Q} = \boldsymbol{Q}^\text{T}\boldsymbol{A}\boldsymbol{Q} = \boldsymbol{\Lambda}Q?1AQ=QTAQ=Λ(適用于實對稱矩陣)
- 完成“求可逆矩陣 P\boldsymbol{P}P”的步驟1-2 得特征值 λ1,…,λn\lambda_1, \dots, \lambda_nλ1?,…,λn? 和對應特征向量 ξi1,…,ξiki\boldsymbol{\xi}_{i1}, \dots, \boldsymbol{\xi}_{ik_i}ξi1?,…,ξiki??。
- 正交化 對同一特征值 λi\lambda_iλi? 的線性無關特征向量 ξi1,…,ξiki\boldsymbol{\xi}_{i1}, \dots, \boldsymbol{\xi}_{ik_i}ξi1?,…,ξiki??,用施密特正交化法化為正交向量組: βi1=ξi1,βij=ξij?∑m=1j?1(ξij,βim)(βim,βim)βim(j=2,…,ki)\boldsymbol{\beta}_{i1} = \boldsymbol{\xi}_{i1}, \quad \boldsymbol{\beta}_{ij} = \boldsymbol{\xi}_{ij} - \sum_{m=1}^{j-1} \frac{(\boldsymbol{\xi}_{ij}, \boldsymbol{\beta}_{im})}{(\boldsymbol{\beta}_{im}, \boldsymbol{\beta}_{im})}\boldsymbol{\beta}_{im} \quad (j=2, \dots, k_i)βi1?=ξi1?,βij?=ξij??m=1∑j?1?(βim?,βim?)(ξij?,βim?)?βim?(j=2,…,ki?)
- 單位化 將正交向量組 βi1,…,βiki\boldsymbol{\beta}_{i1}, \dots, \boldsymbol{\beta}_{ik_i}βi1?,…,βiki?? 單位化: γij=βij∥βij∥(∥β∥=(β,β)為向量模長)\boldsymbol{\gamma}_{ij} = \frac{\boldsymbol{\beta}_{ij}}{\|\boldsymbol{\beta}_{ij}\|} \quad (\|\boldsymbol{\beta}\| = \sqrt{(\boldsymbol{\beta}, \boldsymbol{\beta})} \text{ 為向量模長})γij?=∥βij?∥βij??(∥β∥=(β,β)??為向量模長)
- 構造正交矩陣 Q\boldsymbol{Q}Q 將所有單位正交特征向量按列排列,組成正交矩陣: Q=(γ11,…,γ1k1,γ21,…,γnkn)\boldsymbol{Q} = (\boldsymbol{\gamma}_{11}, \dots, \boldsymbol{\gamma}_{1k_1}, \boldsymbol{\gamma}_{21}, \dots, \boldsymbol{\gamma}_{nk_n})Q=(γ11?,…,γ1k1??,γ21?,…,γnkn??)
【例】
求實對稱矩陣 A=(122212221)\boldsymbol{A} = \begin{pmatrix} 1 & 2 & 2 \\ 2 & 1 & 2 \\ 2 & 2 & 1 \end{pmatrix}A=?122?212?221?? 的正交矩陣 Q\boldsymbol{Q}Q
-
特征值:λ1=5\lambda_1 = 5λ1?=5,λ2=λ3=?1\lambda_2 = \lambda_3 = -1λ2?=λ3?=?1(代數重數均等于幾何重數)。
-
特征向量:
-
λ1=5\lambda_1 = 5λ1?=5 對應 ξ1=(111)\boldsymbol{\xi}_1 = \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix}ξ1?=?111??
-
λ2=?1\lambda_2 = -1λ2?=?1 對應 ξ2=(?110),ξ3=(?101)\boldsymbol{\xi}_2 = \begin{pmatrix} -1 \\ 1 \\ 0 \end{pmatrix}, \boldsymbol{\xi}_3 = \begin{pmatrix} -1 \\ 0 \\ 1 \end{pmatrix}ξ2?=??110??,ξ3?=??101??
-
-
正交化:
- β1=ξ1=(111)\boldsymbol{\beta}_1 = \boldsymbol{\xi}_1 = \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix}β1?=ξ1?=?111??
- β2=ξ2=(?110)\boldsymbol{\beta}_2 = \boldsymbol{\xi}_2 = \begin{pmatrix} -1 \\ 1 \\ 0 \end{pmatrix}β2?=ξ2?=??110??
- β3=ξ3?(ξ3,β2)(β2,β2)β2=(?1/2?1/21)\boldsymbol{\beta}_3 = \boldsymbol{\xi}_3 - \frac{(\boldsymbol{\xi}_3, \boldsymbol{\beta}_2)}{(\boldsymbol{\beta}_2, \boldsymbol{\beta}_2)}\boldsymbol{\beta}_2 = \begin{pmatrix} -1/2 \\ -1/2 \\ 1 \end{pmatrix}β3?=ξ3??(β2?,β2?)(ξ3?,β2?)?β2?=??1/2?1/21??
-
單位化: γ1=13(111),γ2=12(?110),γ3=16(?1?12)\boldsymbol{\gamma}_1 = \frac{1}{\sqrt{3}}\begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix}, \quad \boldsymbol{\gamma}_2 = \frac{1}{\sqrt{2}}\begin{pmatrix} -1 \\ 1 \\ 0 \end{pmatrix}, \quad \boldsymbol{\gamma}_3 = \frac{1}{\sqrt{6}}\begin{pmatrix} -1 \\ -1 \\ 2 \end{pmatrix}γ1?=3?1??111??,γ2?=2?1???110??,γ3?=6?1???1?12??
-
正交矩陣: Q=(1/3?1/2?1/61/31/2?1/61/302/6)\boldsymbol{Q} = \begin{pmatrix} 1/\sqrt{3} & -1/\sqrt{2} & -1/\sqrt{6} \\ 1/\sqrt{3} & 1/\sqrt{2} & -1/\sqrt{6} \\ 1/\sqrt{3} & 0 & 2/\sqrt{6} \end{pmatrix}Q=?1/3?1/3?1/3???1/2?1/2?0??1/6??1/6?2/6??? 滿足 QTAQ=(5000?1000?1)\boldsymbol{Q}^\text{T}\boldsymbol{A}\boldsymbol{Q} = \begin{pmatrix} 5 & 0 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & -1 \end{pmatrix}QTAQ=?500?0?10?00?1??。
- P\boldsymbol{P}P 是可逆矩陣,由線性無關特征向量組成,適用于所有可對角化矩陣;
- Q\boldsymbol{Q}Q 是正交矩陣(Q?1=QT\boldsymbol{Q}^{-1} = \boldsymbol{Q}^\text{T}Q?1=QT),由單位正交特征向量組成,僅適用于實對稱矩陣(必可對角化且可正交對角化)。
十八、二次型化標準型
1、拉格朗日配方法
通過代數配方將二次型 f=xTAxf = \boldsymbol{x}^\text{T}\boldsymbol{A}\boldsymbol{x}f=xTAx 轉化為只含平方項的標準形 f=d1y12+d2y22+?+dnyn2f = d_1y_1^2 + d_2y_2^2 + \cdots + d_ny_n^2f=d1?y12?+d2?y22?+?+dn?yn2?,對應可逆線性變換 x=Cy\boldsymbol{x} = \boldsymbol{C}\boldsymbol{y}x=Cy(C\boldsymbol{C}C 為可逆矩陣)
- 含平方項的變量優先配方:
若二次型含某個變量(如 x1x_1x1?)的平方項,將所有含 x1x_1x1? 的項集中,配成完全平方形式,剩余項中重復此操作。
-
不含平方項時構造平方項:
若二次型僅含交叉項(如 x1x2x_1x_2x1?x2?),令 x1=y1+y2x_1 = y_1 + y_2x1?=y1?+y2?,x2=y1?y2x_2 = y_1 - y_2x2?=y1??y2?,xi=yi(i≥3)x_i = y_i \ (i \geq 3)xi?=yi??(i≥3),引入平方項后再配方。
-
寫出標準形和變換矩陣:
配方后得到標準形,根據變量替換關系寫出可逆矩陣 C\boldsymbol{C}C,滿足 f=yT(CTAC)yf = \boldsymbol{y}^\text{T}(\boldsymbol{C}^\text{T}\boldsymbol{A}\boldsymbol{C})\boldsymbol{y}f=yT(CTAC)y 為標準形。
示例:
化二次型 f=x12+2x1x2+2x22+4x2x3+4x32f = x_1^2 + 2x_1x_2 + 2x_2^2 + 4x_2x_3 + 4x_3^2f=x12?+2x1?x2?+2x22?+4x2?x3?+4x32? 為標準形。
-
配方過程:
f=(x12+2x1x2+x22)+(x22+4x2x3+4x32)=(x1+x2)2+(x2+2x3)2\begin{align*} f &= (x_1^2 + 2x_1x_2 + x_2^2) + (x_2^2 + 4x_2x_3 + 4x_3^2) \\ &= (x_1 + x_2)^2 + (x_2 + 2x_3)^2 \end{align*} f?=(x12?+2x1?x2?+x22?)+(x22?+4x2?x3?+4x32?)=(x1?+x2?)2+(x2?+2x3?)2? -
變量替換:
令 y1=x1+x2y_1 = x_1 + x_2y1?=x1?+x2?,y2=x2+2x3y_2 = x_2 + 2x_3y2?=x2?+2x3?,y3=x3y_3 = x_3y3?=x3?,則 x=(1?1201?2001)y\boldsymbol{x} = \begin{pmatrix} 1 & -1 & 2 \\ 0 & 1 & -2 \\ 0 & 0 & 1 \end{pmatrix}\boldsymbol{y}x=?100??110?2?21??y。 -
標準形:f=y12+y22f = y_1^2 + y_2^2f=y12?+y22?,變換矩陣 C=(1?1201?2001)\boldsymbol{C} = \begin{pmatrix} 1 & -1 & 2 \\ 0 & 1 & -2 \\ 0 & 0 & 1 \end{pmatrix}C=?100??110?2?21??(可逆)。
2、正交化法
通過正交變換 x=Qy\boldsymbol{x} = \boldsymbol{Q}\boldsymbol{y}x=Qy(Q\boldsymbol{Q}Q 為正交矩陣)將二次型化為標準形,標準形的系數為矩陣 A\boldsymbol{A}A 的特征值,即 f=λ1y12+λ2y22+?+λnyn2f = \lambda_1y_1^2 + \lambda_2y_2^2 + \cdots + \lambda_ny_n^2f=λ1?y12?+λ2?y22?+?+λn?yn2?。
- 寫出二次型矩陣 A\boldsymbol{A}A:
二次型 f=xTAxf = \boldsymbol{x}^\text{T}\boldsymbol{A}\boldsymbol{x}f=xTAx 中,A\boldsymbol{A}A 為實對稱矩陣(aiia_{ii}aii? 是 xi2x_i^2xi2? 的系數,aij=ajia_{ij} = a_{ji}aij?=aji? 是 xixjx_ix_jxi?xj? 系數的一半)。
-
求 A\boldsymbol{A}A 的特征值和特征向量:
解特征方程 ∣λE?A∣=0|\lambda\boldsymbol{E} - \boldsymbol{A}| = 0∣λE?A∣=0 得特征值 λ1,…,λn\lambda_1, \dots, \lambda_nλ1?,…,λn?,對應特征向量 ξ1,…,ξn\boldsymbol{\xi}_1, \dots, \boldsymbol{\xi}_nξ1?,…,ξn?。 -
特征向量正交化與單位化:
對同一特征值的線性無關特征向量用施密特正交化,再將所有特征向量單位化,得單位正交向量組 γ1,…,γn\boldsymbol{\gamma}_1, \dots, \boldsymbol{\gamma}_nγ1?,…,γn?。
-
構造正交矩陣 Q\boldsymbol{Q}Q 和標準形:
Q=(γ1,…,γn)\boldsymbol{Q} = (\boldsymbol{\gamma}_1, \dots, \boldsymbol{\gamma}_n)Q=(γ1?,…,γn?),則正交變換 x=Qy\boldsymbol{x} = \boldsymbol{Q}\boldsymbol{y}x=Qy 化二次型為標準形:
f=λ1y12+λ2y22+?+λnyn2f = \lambda_1y_1^2 + \lambda_2y_2^2 + \cdots + \lambda_ny_n^2 f=λ1?y12?+λ2?y22?+?+λn?yn2?
示例:
用正交化法化 f=2x12+2x22+2x32+2x1x2+2x1x3+2x2x3f = 2x_1^2 + 2x_2^2 + 2x_3^2 + 2x_1x_2 + 2x_1x_3 + 2x_2x_3f=2x12?+2x22?+2x32?+2x1?x2?+2x1?x3?+2x2?x3? 為標準形。
-
二次型矩陣:A=(211121112)\boldsymbol{A} = \begin{pmatrix} 2 & 1 & 1 \\ 1 & 2 & 1 \\ 1 & 1 & 2 \end{pmatrix}A=?211?121?112??。
-
特征值:λ1=4\lambda_1 = 4λ1?=4,λ2=λ3=1\lambda_2 = \lambda_3 = 1λ2?=λ3?=1。
-
單位正交特征向量:
γ1=13(1,1,1)T\boldsymbol{\gamma}_1 = \frac{1}{\sqrt{3}}(1, 1, 1)^\text{T}γ1?=3?1?(1,1,1)T,γ2=12(?1,1,0)T\boldsymbol{\gamma}_2 = \frac{1}{\sqrt{2}}(-1, 1, 0)^\text{T}γ2?=2?1?(?1,1,0)T,γ3=16(?1,?1,2)T\boldsymbol{\gamma}_3 = \frac{1}{\sqrt{6}}(-1, -1, 2)^\text{T}γ3?=6?1?(?1,?1,2)T。
- 標準形:f=4y12+y22+y32f = 4y_1^2 + y_2^2 + y_3^2f=4y12?+y22?+y32?,正交矩陣 Q=(γ1,γ2,γ3)\boldsymbol{Q} = (\boldsymbol{\gamma}_1, \boldsymbol{\gamma}_2, \boldsymbol{\gamma}_3)Q=(γ1?,γ2?,γ3?)。
十九、二次型正定
若n元二次型f=xTAxf=x^TAxf=xTAx正定 <=> 對任意x≠0,有f=xTAxf=x^TAxf=xTAx>0
<=> f的正慣性指數p = n
<=> 存在可逆矩陣D,使 A = DT^TTD
<=>A合同與E
<=>A的特征值λi>0(i=1,2,???,n)λ_i >0 (i = 1, 2,···,n)λi?>0(i=1,2,???,n)
<=>A的全部順序主子式均大于0
二十、等價、相似、合同
關系 | 等價(矩陣A\boldsymbol{A}A與B\boldsymbol{B}B等價) | 相似(矩陣A\boldsymbol{A}A與B\boldsymbol{B}B相似) | 合同(矩陣A\boldsymbol{A}A與B\boldsymbol{B}B合同) |
---|---|---|---|
定義 | 存在可逆矩陣P,Q\boldsymbol{P},\boldsymbol{Q}P,Q,使B=PAQ\boldsymbol{B}=\boldsymbol{P}\boldsymbol{A}\boldsymbol{Q}B=PAQ | 存在可逆矩陣P\boldsymbol{P}P,使B=P?1AP\boldsymbol{B}=\boldsymbol{P}^{-1}\boldsymbol{A}\boldsymbol{P}B=P?1AP | 存在可逆矩陣C\boldsymbol{C}C,使B=CTAC\boldsymbol{B}=\boldsymbol{C}^\text{T}\boldsymbol{A}\boldsymbol{C}B=CTAC |
核心本質 | 矩陣經初等變換可互化(體現秩的一致性) | 線性變換在不同基下的矩陣表示(保持特征值等核心屬性) | 二次型經可逆線性變換的等價性(保持正定性等慣性性質) |
充要條件 | 同型且秩相等:r(A)=r(B)r(\boldsymbol{A})=r(\boldsymbol{B})r(A)=r(B) | ① 特征值完全相同(含重數); ② 存在可逆矩陣P,Q\boldsymbol{P},\boldsymbol{Q}P,Q使A=PBQ\boldsymbol{A}=\boldsymbol{P}\boldsymbol{B}\boldsymbol{Q}A=PBQ且P?1=Q\boldsymbol{P}^{-1}=\boldsymbol{Q}P?1=Q(特殊等價) | ① 均為實對稱矩陣且慣性指數相同(正、負慣性指數分別相等); ② 存在可逆矩陣P,Q\boldsymbol{P},\boldsymbol{Q}P,Q使A=PBQ\boldsymbol{A}=\boldsymbol{P}\boldsymbol{B}\boldsymbol{Q}A=PBQ且PT=Q\boldsymbol{P}^\text{T}=\boldsymbol{Q}PT=Q(特殊等價) |
包含關系 | 等價是最寬泛的關系: 相似?\subset?等價,合同?\subset?等價(實對稱矩陣中相似?\subset?合同) | 相似矩陣必等價,但等價矩陣不一定相似; 實對稱矩陣相似必合同,但合同不一定相似 | 合同矩陣必等價,但等價矩陣不一定合同; 實對稱矩陣合同?\nRightarrow?相似(特征值可不同) |
不變量 | 矩陣的秩r(A)r(\boldsymbol{A})r(A) | 特征值、行列式、跡、秩、可逆性 | 慣性指數(正慣性指數ppp、負慣性指數qqq)、秩、對稱性(若原矩陣對稱) |
適用場景 | 矩陣秩的比較、方程組同解性等 | 特征值與特征向量、矩陣對角化、線性變換等 | 二次型化簡、正定性判定、曲面分類等 |
示例 | (1000)\begin{pmatrix}1&0\\0&0\end{pmatrix}(10?00?)與(0100)\begin{pmatrix}0&1\\0&0\end{pmatrix}(00?10?)等價(秩均為1) | (1101)\begin{pmatrix}1&1\\0&1\end{pmatrix}(10?11?)與(1011)\begin{pmatrix}1&0\\1&1\end{pmatrix}(11?01?)相似(特征值均為1) | (100?1)\begin{pmatrix}1&0\\0&-1\end{pmatrix}(10?0?1?)與(200?3)\begin{pmatrix}2&0\\0&-3\end{pmatrix}(20?0?3?)合同(慣性指數均為p=1,q=1p=1,q=1p=1,q=1) |