文章目录
- 十一、方程组解的结构和性质
- 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} ⎩⎨⎧a11x1+a12x2+⋯+a1nxn=0,a21x1+a22x2+⋯+a2nxn=0,⋯⋯am1x1+am2x2+⋯+amnxn=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} ⎩⎨⎧a11x1+a12x2+⋯+a1nxn=b1,a21x1+a22x2+⋯+a2nxn=b2,⋯⋯am1x1+am2x2+⋯+amnxn=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]=101011−11001−10−101
对矩阵作初等行变换
[1001010−1001−20000]\begin{bmatrix} 1 & 0 & 0 & 1 \\ 0 & 1 & 0 & -1 \\ 0 & 0 & 1 & -2 \\ 0 & 0 & 0 & 0 \end{bmatrix}1000010000101−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=k10010+k2−1101=−k2k2k1k2
(其中 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=−k2k2k1k2代入方程组(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=k10010+k2−1101=−k2k2k1k2
(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=l10110+l2−1−101=−l2l1−l2l1l2
由上式可得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(λiE−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λ1y12+λ2y22+⋯+λnyn2,其中 λ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}(λiE−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}(λiE−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=122212221 的正交矩阵 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=31111,γ2=21−110,γ3=61−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/31/31/3−1/21/20−1/6−1/62/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=5000−1000−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=d1y12+d2y22+⋯+dnyn2,对应可逆线性变换 x=Cy\boldsymbol{x} = \boldsymbol{C}\boldsymbol{y}x=Cy(C\boldsymbol{C}C 为可逆矩阵)
- 含平方项的变量优先配方:
若二次型含某个变量(如 x1x_1x1)的平方项,将所有含 x1x_1x1 的项集中,配成完全平方形式,剩余项中重复此操作。
-
不含平方项时构造平方项:
若二次型仅含交叉项(如 x1x2x_1x_2x1x2),令 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+2x1x2+2x22+4x2x3+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+2x1x2+x22)+(x22+4x2x3+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−1102−21y。 -
标准形: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−1102−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=λ1y12+λ2y22+⋯+λnyn2。
- 写出二次型矩阵 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_jxixj 系数的一半)。
-
求 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=λ1y12+λ2y22+⋯+λnyn2
示例:
用正交化法化 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+2x1x2+2x1x3+2x2x3 为标准形。
-
二次型矩阵:A=(211121112)\boldsymbol{A} = \begin{pmatrix} 2 & 1 & 1 \\ 1 & 2 & 1 \\ 1 & 1 & 2 \end{pmatrix}A=211121112。
-
特征值:λ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=31(1,1,1)T,γ2=12(−1,1,0)T\boldsymbol{\gamma}_2 = \frac{1}{\sqrt{2}}(-1, 1, 0)^\text{T}γ2=21(−1,1,0)T,γ3=16(−1,−1,2)T\boldsymbol{\gamma}_3 = \frac{1}{\sqrt{6}}(-1, -1, 2)^\text{T}γ3=61(−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}(1000)与(0100)\begin{pmatrix}0&1\\0&0\end{pmatrix}(0010)等价(秩均为1) | (1101)\begin{pmatrix}1&1\\0&1\end{pmatrix}(1011)与(1011)\begin{pmatrix}1&0\\1&1\end{pmatrix}(1101)相似(特征值均为1) | (100−1)\begin{pmatrix}1&0\\0&-1\end{pmatrix}(100−1)与(200−3)\begin{pmatrix}2&0\\0&-3\end{pmatrix}(200−3)合同(惯性指数均为p=1,q=1p=1,q=1p=1,q=1) |