【同济线性代数第五版】第五章 相似矩阵及二次型

本文涉及知识点

数学

1 向量的内积、长度及正交性

本书用[x,y]表示向量x和向量y的内积(点乘)。

定理1:若n维向量 a 1 , a 2 , ⋯ a r a_1,a_2,\cdots a_r a1,a2,ar是一组两两正交的非零向量,则 a 1 , a 2 ⋯ a r a_1,a_2\cdots a_r a1,a2ar线性无关。
: 设有 λ 1 , l a m b d a 2 , ⋯   , λ r \lambda1,lambda2,\cdots,\lambda_r λ1,lambda2,,λr使
λ 1 a 1 + λ 2 a 2 + ⋯ λ r a r = 0 \lambda_1a_1+\lambda_2a_2+\cdots \lambda_ra_r=0 λ1a1+λ2a2+λrar=0
a 1 T 左乘上式两端,因当 i ≥ 2 时, a 1 T a i = 0 故得: λ i a i T a i = 0 , 故 λ 1 a 1 T a 1 = 0 , 因为 a i ≠ 0 , 故 a 1 T a 1 = ∣ ∣ a 1 ∣ ∣ 2 ≠ 0 , 故 λ 1 = 0 a_1^T左乘上式两端,因当i \ge 2时,a_1^Ta_i=0故得:\lambda_i a_i^Ta_i=0,故\lambda_1 a_1^Ta_1=0,因为a_i\neq 0,故a_1^Ta_1=||a_1||^2\neq 0,故\lambda_1 =0 a1T左乘上式两端,因当i2时,a1Tai=0故得:λiaiTai=0,λ1a1Ta1=0,因为ai=0,a1Ta1=∣∣a12=0,λ1=0
同理 λ i 全部为 0 \lambda_i 全部为0 λi全部为0故线性无关。

例1 已知3维向量空间 R 3 中两个向量 \R^3中两个向量 R3中两个向量
a 1 = ( 1 1 1 ) , a 2 = ( 1 − 2 1 ) a_1=\begin{pmatrix} 1\\1\\1\\ \end{pmatrix}, a_2=\begin{pmatrix} 1\\-2\\1\\ \end{pmatrix} a1= 111 ,a2= 121
正交,试求一个非零向量a_3,使得 a 1 , a 2 , a 3 a_1,a_2,a_3 a1,a2,a3两两正交。
: 此向量就是 ( a 1 T , a 2 T ) x = 0 (a_1^T,a_2^T)x=0 (a1T,a2T)x=0的任意非零解。
定义3 设n维向量 e 1 , e 2 ⋯ e r e_1,e_2\cdots e_r e1,e2er是向量空间V ( V ⊂ R n ) (V\sub \R^n) (VRn)的一个基,如果 e 1 , e 2 , ⋯ e r e_1,e_2,\cdots e_r e1,e2,er两两正交,且都是单位向量,则称 e 1 , ⋯   , e r e1,\cdots ,e_r e1,,er是V的一个规范正交基。

例2 设 a 1 = ( 1 2 − 1 ) , a 2 = ( − 1 3 1 ) , a 3 = ( 4 − 1 0 ) a_1=\begin{pmatrix}1\\2\\-1\end{pmatrix},a_2=\begin{pmatrix}-1\\3\\1\end{pmatrix},a_3=\begin{pmatrix}4\\-1\\0\end{pmatrix} a1= 121 ,a2= 131 ,a3= 410 试用施密特正交化过程把这组向量规范正交化。
b 1 = a 1 b_1=a_1 b1=a1
b 2 = a 2 − [ a 2 , b 1 ] ∣ ∣ b 1 ∣ ∣ 2 b 1 = ( − 1 3 1 ) − 4 6 ( 1 2 − 1 ) = 5 3 ( − 1 1 1 ) b_2=a_2-\frac{[a_2,b_1]}{||b1||^2}b_1=\begin{pmatrix}-1\\3\\1\end{pmatrix}-\frac 4 6 \begin{pmatrix}1\\2\\-1\end{pmatrix}=\frac 5 3\begin{pmatrix}-1\\1\\1\end{pmatrix} b2=a2∣∣b1∣2[a2,b1]b1= 131 64 121 =35 111
b 3 = a 3 − [ a 3 , b 1 ] ∣ ∣ b 1 ∣ ∣ 2 b 1 − [ a 3 , b 2 ] ∣ ∣ b 2 ∣ ∣ 2 b 2 b_3=a_3-\frac{[a_3,b_1]}{||b1||^2}b_1-\frac{[a_3,b_2]}{||b2||^2}b_2 b3=a3∣∣b1∣2[a3,b1]b1∣∣b2∣2[a3,b2]b2
= ( 4 − 1 0 ) − 1 3 ( 1 2 − 1 ) + 3 5 ( − 1 1 1 ) = \begin{pmatrix}4\\-1\\0\end{pmatrix}-\frac 1 3 \begin{pmatrix}1\\2\\-1\end{pmatrix}+ \frac 3 5 \begin{pmatrix}-1\\1\\1\end{pmatrix} = 410 31 121 +53 111
= 2 ( 1 0 1 ) =2\begin{pmatrix}1\\0\\1\end{pmatrix} =2 101

定义4 如果n阶矩阵A满足 A T A = E ,即 A − 1 = A T A^TA=E,即A^{-1}=A^T ATA=E,即A1=AT,那么称A为正交矩阵,简称正交阵。
( a 1 T , a 2 T , ⋯ a n T ) ( a 1 , a 2 , ⋯ a n ) = E (a_1^T,a_2^T,\cdots a_n^T)(a_1,a_2,\cdots a_n)=E (a1T,a2T,anT)(a1,a2,an)=E
[ a i T , a j ] = { 1 i = = j a i 是单位向量 0 i ≠ 0 a i 和 a j 正交 [a_i^T,a_j]=\begin{cases}1&i==j& a_i是单位向量\\0&i\neq 0& a_i和a_j正交 \end{cases} [aiT,aj]={10i==ji=0ai是单位向量aiaj正交
一,若A为正交矩阵,则 A − 1 = A T A^{-1}=A^T A1=AT也是正交矩阵,且|A|=1或-1。
∣ A T A ∣ = ∣ E ∣ = 1 , ∣ A T A ∣ = ∣ A T ∣ ∣ A ∣ , 故 ∣ A ∣ 等于 ± 1 |A^TA|=|E|=1,|A^TA|=|A^T||A|,故|A|等于\pm 1 ATA=E=1ATA=AT∣∣A,A等于±1
二,若A和B都是正交阵,则AB也是正交阵。
证明 A B ( A B ) T = A B B T A T = A E A T = A A T = E AB(AB)^T=ABB^TA^T=AEA^T=AA^T=E AB(AB)T=ABBTAT=AEAT=AAT=E
定义5:若P为正交矩阵,则线性变换 y = P x y=Px y=Px称为正交变换。
∣ ∣ y ∣ ∣ = y T y = x T P T P x = x T x = ∣ ∣ x ∣ ∣ ||y||=\sqrt{y^Ty}=\sqrt{x^TP^TPx}=\sqrt{x^Tx}=||x|| ∣∣y∣∣=yTy =xTPTPx =xTx =∣∣x∣∣

2 方阵的特征值与特征向量

定义6 设A是n阶矩阵,如果数 λ 和 n 维非零列向量 x 使关系式 \lambda和n维非零列向量x使关系式 λn维非零列向量x使关系式
A x = λ x Ax=\lambda x Ax=λx成立,那么,这样的数 λ \lambda λ称为矩阵A的特征值,非零向量x称为A的对应与特征值 λ \lambda λ的特征向量。
上式也可写成(A- λ E ) x = 0 \lambda E)x=0 λE)x=0
∣ a 11 − λ a 12 ⋯ a 1 n a 21 a 22 − λ ⋯ a 2 n ⋮ ⋮ ⋮ a 31 a 32 ⋯ a 3 n − λ ∣ \begin{vmatrix} a_{11}-\lambda&a_{12}&\cdots& a_{1n}\\ a_{21}&a_{22}-\lambda&\cdots& a_{2n}\\ \vdots&\vdots&&\vdots\\ a_{31}&a_{32}&\cdots& a_{3n}-\lambda\\ \end{vmatrix} a11λa21a31a12a22λa32a1na2na3nλ
上式是以 λ \lambda λ为未知数的一元n次方程,称为矩阵A的特征方程。其左端是|A-
λ \lambda λ E|是 λ \lambda λ的n次多项式,记作 f ( λ ) f(\lambda) f(λ),称为矩阵A的特征多项式。A的特征值就是特征方程的解。特征方程在复数范围内恒有解,其个数为方程的次数(重根重复记数),因此,n阶矩阵A在复数范围内有n个特征值。
设n阶矩阵A= ( a i j ) 的特征值为 λ 1 , λ 2 ⋯   , λ n (a_ij)的特征值为\lambda_1,\lambda_2\cdots,\lambda_n (aij)的特征值为λ1,λ2,λn
一, λ 1 + λ 2 + ⋯ + λ n = a 11 + a 22 + ⋯ a n n \lambda_1+\lambda_2+\cdots+\lambda_n=a_{11}+a_{22}+\cdots a_{nn} λ1+λ2++λn=a11+a22+ann
二, λ 1 λ 2 ⋯ λ n = ∣ A ∣ \lambda_1 \lambda_2 \cdots \lambda_n=|A| λ1λ2λn=A

定理2:设 λ 1 , λ 2 , ⋯   , λ m \lambda_1,\lambda_2,\cdots ,\lambda_m λ1,λ2,,λm是方阵A的m个特征值, p 1 , p 2 , ⋯ p m 依次是与之对应的特征向量,如果 λ 1 , l a m b d a 2 , ⋯   , λ m 各不相等,则 p 1 , p 2 , ⋯ p m p_1,p_2,\cdots p_m依次是与之对应的特征向量,如果\lambda1,lambda_2,\cdots,\lambda_m各不相等,则p_1,p_2,\cdots p_m p1,p2,pm依次是与之对应的特征向量,如果λ1,lambda2,,λm各不相等,则p1,p2,pm线性无关。
m=1 特征向量不为0,故线性无关。
假定m=k-1线性无关,则m=k线性无关。
a 1 p 1 + a 2 p 2 + ⋯ a k p k = 0 式子一 a_1p_1+a_2p_2+\cdots a_kp_k=0 式子一 a1p1+a2p2+akpk=0式子一假定式子一成立。
式子一左乘以A, a 1 A p 1 + a 2 A p 2 + ⋯ a k A p k = 0 即 a 1 λ 1 p 1 + a 2 λ 2 p 2 + ⋯ a k λ k p k = 0 式子二 a_1Ap_1+a_2Ap_2+\cdots a_kAp_k=0即a_1\lambda_1 p_1+a_2\lambda_2p_2+\cdots a_k\lambda_k p_k=0 式子二 a1Ap1+a2Ap2+akApk=0a1λ1p1+a2λ2p2+akλkpk=0式子二
式子二-式子一乘以 λ k \lambda_k λk= a 1 p 1 ( λ 1 − λ k ) + a 2 p 2 ( λ 2 − λ k ) + ⋯ + a k − 1 p k − 1 ( λ k − 1 − λ k ) = 0 式子 3 a_1p_1(\lambda_1-\lambda_k)+a_2p_2(\lambda_2-\lambda_k)+\cdots+a_{k-1}p_{k-1}(\lambda_{k-1}-\lambda_k)=0式子3 a1p1(λ1λk)+a2p2(λ2λk)++ak1pk1(λk1λk)=0式子3
由于前k-1个向量线性无关,故 a i ( λ i − λ k ) = 0 , i < k ,由于特征值各不相等,故 a i = 0 , i < k a_{i}(\lambda_i-\lambda_k)=0,i< k,由于特征值各不相等,故a_i=0,i<k ai(λiλk)=0,i<k,由于特征值各不相等,故ai=0,i<k
故式子一可以简化为: a k p k = 0 a_kp_k=0 akpk=0,即 a k = 0 a_k=0 ak=0,即式子一成立,则 a i 全为 0 a_i全为0 ai全为0,故前k个特征向量也线性无关。

相似矩阵

定义7 设A,B都是n阶矩阵,若有可逆矩阵P,使 P − 1 A P = B P^{-1}AP=B P1AP=B则称B是A的相似矩阵,或说矩阵A与B相似。对A进行运算 P − 1 A P = B P^{-1}AP=B P1AP=B称为对A进行相似变换。可逆矩阵P称为把A变成B的相似变换矩阵。
自己想到的引理:n阶方阵AB的特征多项式相同,则特征值相同。如果 λ \lambda λ是A的特征值,则特征多项式为0,故也是B的特征值。 λ \lambda λ不是A的特征值,则特征多项式不为0,故不是B的特征值。
定理3:若n阶可逆矩阵A与B相似,则A与B的特征多项式相同,从而A与B的特征值亦相同。
证明:对于任意 λ ∣ B − λ E ∣ = ∣ P − 1 A P − P − 1 λ E P ∣ = ∣ P − 1 ( A − λ E ) P ∣ \lambda \begin{vmatrix}B- \lambda E\end{vmatrix}= \begin{vmatrix}P^{-1}AP-P^{-1}\lambda E P\end{vmatrix}= \begin{vmatrix}P^{-1}(A-\lambda E)P\end{vmatrix} λ BλE = P1APP1λEP = P1(AλE)P
= ∣ P − 1 ∣ ∣ A − λ E ∣ ∣ P ∣ = ∣ A − λ E ∣ =|P^{-1}||A-\lambda E||P|=|A-\lambda E| =P1∣∣AλE∣∣P=AλE
推论 若n阶矩阵A与对角阵 Λ = ( λ 1 λ 2 ⋯ λ n ) \Lambda=\begin{pmatrix} \lambda_1\\&\lambda_2\\&&\cdots\\&&&\lambda_n \end{pmatrix} Λ= λ1λ2λn 相似,则 λ 1 , λ 2 ⋯ λ n \lambda_1,\lambda_2 \cdots \lambda_n λ1,λ2λn即A的n个特征值。

Λ 是 n 阶对角矩阵 , A 是 n 阶矩阵 , 假定存在可逆矩阵 P ,使得 P − 1 A P = Λ \Lambda是n阶对角矩阵,A是n阶矩阵,假定存在可逆矩阵P,使得P^{-1}AP=\Lambda Λn阶对角矩阵,An阶矩阵,假定存在可逆矩阵P,使得P1AP=Λ。式子左乘P,得 A P = P Λ AP=P\Lambda AP=PΛ。将前式的第二个矩阵按列拆开,分别处理。 A P i = P Λ i AP_i=P\Lambda_i APi=PΛi右式的(j,0)元= a ( j , i ) λ i 即 P i λ i a(j,i)\lambda_i即P_i\lambda_i a(j,i)λiPiλi。即 P i P_i Pi是特征向量, λ i \lambda_i λi是特征值。 P i 之间一定线性无关,否则没有逆矩阵 P_i之间一定线性无关,否则没有逆矩阵 Pi之间一定线性无关,否则没有逆矩阵
定理4 n阶矩阵A与对角线阵相似(即A能对角化)的充分必要条件是A有n个线性无关的特征向量。
推论 如果n阶矩阵A的n个特征值互不相等,则A与对角阵相似。
否命题不一定成立。特征值相等,也可能对角化。
特征值不同,则特征向量一定不同,用反证发证明。令 特征值 λ 1 ≠ λ 2 , 对应的特征向量都是 p , λ 1 p ≠ λ 2 , λ 1 p = A p , λ 2 p = A p 特征值\lambda_1 \neq \lambda_2,对应的特征向量都是p,\lambda_1 p \neq \lambda_2,\lambda_1p=Ap,\lambda_2p=Ap 特征值λ1=λ2,对应的特征向量都是p,λ1p=λ2,λ1p=Ap,λ2p=Ap矛盾。

4 对称矩阵的对角化

定理5:对称矩阵的特征值为实数。
定理6:设 λ 1 , λ 2 \lambda_1,\lambda_2 λ1,λ2是对称矩阵A的两个特征值, p 1 , p 2 p_1,p_2 p1,p2是对应的特征向量。若 λ 1 ≠ λ 2 , \lambda_1\neq \lambda_2, λ1=λ2, p 1 与 p 2 p_1与p_2 p1p2正交。
: λ 1 p 1 T = ( λ 1 p 1 ) T = ( A p 1 ) T = p 1 T A T = p 1 T A \lambda_1 p_1^T=(\lambda_1 p_1)^T=(Ap_1)^T=p_1^TA^T=p_1^TA λ1p1T=(λ1p1)T=(Ap1)T=p1TAT=p1TA
λ 1 p 1 T p 2 = p 1 T A p 2 = p 1 T ( λ 2 p 2 ) = λ 2 p 1 T p 2 \lambda_1p_1^Tp_2=p_1^TAp_2=p_1^T(\lambda_2p_2)=\lambda_2p_1^Tp_2 λ1p1Tp2=p1TAp2=p1T(λ2p2)=λ2p1Tp2
即: ( λ 1 − λ 2 ) p 1 T p 2 = 0 (\lambda_1-\lambda_2)p_1^Tp_2=0 (λ1λ2)p1Tp2=0
λ 1 ≠ λ 2 , 故 p 1 T p 2 = 0 即 p 1 与 p 2 正交 \lambda_1\neq \lambda_2,故p_1^Tp_2=0即p_1与p_2正交 λ1=λ2,p1Tp2=0p1p2正交
定理7 设A为n阶对称阵,则必有正交阵P,使得 P − 1 A P = P T A P = Λ , 其中 Λ 是以 A 的 n 个特征值为对角元的对角阵 P^{-1}AP=P^TAP=\Lambda,其中\Lambda是以A的n个特征值为对角元的对角阵 P1AP=PTAP=Λ,其中Λ是以An个特征值为对角元的对角阵

推论 设A为n阶段对称阵, λ 是 A 的特征方程的重根,则矩阵 A − λ E 的秩 = n − k , 从而对应的特征值 λ 恰有 k 个线性无关的特征向量 \lambda是A的特征方程的重根,则矩阵A-\lambda E的秩=n-k,从而对应的特征值\lambda 恰有k个线性无关的特征向量 λA的特征方程的重根,则矩阵AλE的秩=nk,从而对应的特征值λ恰有k个线性无关的特征向量

一,求出A的全部互不相等的特征值 λ 1 , ⋯   , λ s \lambda_1,\cdots ,\lambda_s λ1,,λs,它的重数分别为 k 1 , k 2 , ⋯ k s ( k 1 + ⋯ + k + s = n ) k_1,k_2,\cdots k_s(k_1+\cdots +k+s=n) k1,k2,ks(k1++k+s=n)
二,对每个 k i k_i ki重特征值 λ i \lambda_i λi,求方程 ( A − λ i E ) x = 0 (A-\lambda_iE)x=0 (AλiE)x=0的基础解系,得 k i 个线性无关的特征向量。在把它们正交化、单位化,得 k i 个两两正交的特征向量。共得到 n 个两两正交的单位特征向量 k_i个线性无关的特征向量。在把它们正交化、单位化,得k_i个两两正交的特征向量。共得到n个两两正交的单位特征向量 ki个线性无关的特征向量。在把它们正交化、单位化,得ki个两两正交的特征向量。共得到n个两两正交的单位特征向量
三,把这n个两两正交的单位特征向量构成的正交阵P,便有 P − 1 A P = P T A P = Λ . 注意 Λ 中对角元的排列顺序与 P 中列向量的排列次序相对应。 P^{-1}AP=P^TAP=\Lambda.注意\Lambda中对角元的排列顺序与P中列向量的排列次序相对应。 P1AP=PTAP=Λ.注意Λ中对角元的排列顺序与P中列向量的排列次序相对应。

5 二次型及其标准形

定义8:含有n个变量 x 1 , x 2 , ⋯ x n x_1,x_2,\cdots x_n x1,x2,xn的二次齐次函数
f ( x 1 , x 2 , ⋯   , x n ) = a 11 x 1 2 + a 22 x 2 2 + ⋯ + a n n x n 2 f(x_1,x_2,\cdots ,x_n)=a_{11}x_1^2+a_{22}x_2^2+\cdots + a_{nn}x_n^2 f(x1,x2,,xn)=a11x12+a22x22++annxn2
+ 2 a 12 x 1 x 2 + 2 a 13 x 1 x 3 + ⋯ + 2 a n − 1 , n x n − 1 x n +2a_{12}x_1x_2+2a_{13}x_1x_3+\cdots+2a_{n-1,n}x_{n-1}x_n +2a12x1x2+2a13x1x3++2an1,nxn1xn称为二次型。式子5
a j i = a i j , 则 2 a i j x i x j = a i j x i x j + a j i x i x j a_{ji}=a_{ij},则2a_{ij}x_ix_j=a_{ij}x_ix_j+a_{ji}x_ix_j aji=aij,2aijxixj=aijxixj+ajixixj则上式可写成:
= ∑ i , j = 1 n a i j x i x j =\sum\limits_{i,j=1}^na_{ij}x_ix_j =i,j=1naijxixj式子6
对于二次型,我们讨论的主要问题是:寻求可逆的线性变换
{ x 1 = c 11 y 1 + c 12 y 2 + ⋯ + c 1 n y n , x 2 = c 21 y 1 + c 22 y 2 + ⋯ + c 2 n y n , ⋯ ⋯ ⋯ x n = c n 1 y 1 + c n 2 y 2 + ⋯ + c n n y n \begin{cases} x_1=c_{11}y_1+c_{12}y_2+\cdots+ c_{1n}y_n,\\ x_2=c_{21}y_1+c_{22}y_2+\cdots + c_{2n}y_n,\\ \cdots\cdots\cdots\\ x_n=c_{n1}y_1+c_{n2}y_2+\cdots +c_{nn}y_n \end{cases} x1=c11y1+c12y2++c1nyn,x2=c21y1+c22y2++c2nyn,⋯⋯⋯xn=cn1y1+cn2y2++cnnyn式子7
使得二次型只含平方项,也就是式子7带入式子5,能使
f = k 1 y 1 2 + k 2 y 2 2 + ⋯ + k n y n 2 f=k_1y_1^2+k_2y_2^2+\cdots +k_ny_n^2 f=k1y12+k2y22++knyn2
这种只含平方项的二次型,称为二次型的标准型(或法式)。
如果标准型的系数 k 1 , k 2 ⋯ k n k_1,k_2\cdots k_n k1,k2kn只有1,-1,0三个数中取值,则称二次型的规范型。

我们把对称阵A叫做二次型f的矩阵,也把f叫做对称阵A的二次型。对称阵A的秩就叫做二次型f的秩。

定义9 设A和B的n阶矩阵,若有可逆矩阵C,使得B= C T A C C^TAC CTAC,则称矩阵A与B合同。

若A为对称阵,则 B = C T A C B=C^TAC B=CTAC也为对称阵,且R(B)=R(A)。
B T = ( C T A C ) T = C T A C = B B^T=(C^TAC)^T=C^TAC=B BT=(CTAC)T=CTAC=B
因为C可逆,故 C T C^T CT可逆,故AB等价,故秩相等。

定理8 任何二次型 f = ∑ i , j = 1 n a i j x i x y f=\sum\limits_{i,j=1}^na_{ij}x_ix_y f=i,j=1naijxixy,总有正交变换 x = P y x=Py x=Py,使f化为标准形 f = λ 1 y 1 2 + λ 2 y 2 2 + ⋯ λ n y n 2 f=\lambda_1y_1^2+\lambda_2y_2^2+\cdots \lambda_ny_n^2 f=λ1y12+λ2y22+λnyn2其中 λ 1 , l a b m d a 2 , ⋯   , λ n \lambda_1,labmda_2,\cdots,\lambda_n λ1,labmda2,,λn是f的矩阵 A = ( a i j ) 的特征值 A=(a_ij)的特征值 A=(aij)的特征值
推论 任给n元二次型 f ( x ) = x T A x ( A T ) = A f(x)=x^TAx(A^T)=A f(x)=xTAx(AT)=A,总有可逆变换x=Cz,使得f(Cz)为规范型。

例14:求一个正交变换 x = P y x=Py x=Py,把二次型 f = − 2 x 1 x 2 + 2 x 1 x 3 + 2 x 2 x 3 f=-2x_1x_2+2x_1x_3+2x_2x_3 f=2x1x2+2x1x3+2x2x3为化标准形。
二次型的矩阵为: A = ( 0 − 1 1 − 1 0 1 1 1 0 ) A=\begin{pmatrix}0&-1&1\\ -1&0&1\\ 1&1&0\\ \end{pmatrix} A= 011101110 此矩阵即例12的矩阵,有正交阵:
P = ( − 1 3 − 1 2 1 6 − 1 3 1 2 1 6 1 3 0 2 6 ) P=\begin{pmatrix}-\frac 1 {\sqrt 3}&-\frac 1 {\sqrt 2}&\frac 1 {\sqrt 6}\\ -\frac 1 {\sqrt 3}&\frac 1 {\sqrt 2}&\frac 1 {\sqrt 6}\\ \frac 1 {\sqrt 3}&0&\frac 2 {\sqrt 6}\\ \end{pmatrix} P= 3 13 13 12 12 106 16 16 2 使得 P T A P = Λ = ( − 2 0 0 0 1 0 0 0 1 ) P^TAP=\Lambda=\begin{pmatrix}-2&0&0\\ 0&1&0\\ 0&0&1\\ \end{pmatrix} PTAP=Λ= 200010001
把二次型f化成标准形 f = − 2 y 1 2 + y 2 2 + y 3 2 f=-2y_1^2+y_2^2+y_3^2 f=2y12+y22+y32
如果要把二次型f化成规范型,只需令
{ y 1 = 1 2 z 1 , y 2 = z 2 , y 3 = z 3 , \begin{cases} y1 = \frac 1 {\sqrt 2}z_1,\\ y_2=z2,\\ y3 = z3,\\ \end{cases} y1=2 1z1,y2=z2,y3=z3,
f= − z 1 2 + z 2 2 + z 3 2 -z_1^2+z_2^2+z_3^2 z12+z22+z32

6 用配方发化二次型称标准型

例15 化二次型 f = x 1 2 + 2 x 2 2 + 5 x 3 2 + 2 x 1 x 2 + 2 x 1 x 3 + 6 x 2 x 3 f=x_1^2+2x_2^2+5x_3^2+2x_1x_2+2x_1x_3+6x_2x_3 f=x12+2x22+5x32+2x1x2+2x1x3+6x2x3称标准形,并求所用的变换矩阵。
( x 1 + x 2 + x 3 ) 2 − x 2 2 − x 3 2 − 2 x 2 x 3 + 2 x 2 2 + 5 x 3 2 + 6 x 2 x 3 (x_1+x_2+x_3)^2-x_2^2-x_3^2-2x_2x_3+2x_2^2+5x_3^2+6x_2x_3 (x1+x2+x3)2x22x322x2x3+2x22+5x32+6x2x3
= ( x 1 + x 2 + x 3 ) 2 + x 2 2 + 4 x 2 x 3 + 4 x 3 2 (x_1+x_2+x_3)^2+x_2^2+4x_2x_3+4x_3^2 (x1+x2+x3)2+x22+4x2x3+4x32
= ( x 1 + x 2 + x 3 ) 2 + ( x 2 + 2 x 3 2 ) 2 (x_1+x_2+x_3)^2+(x_2+2x_3^2)^2 (x1+x2+x3)2+(x2+2x32)2

{ y 1 = x 1 + x 2 + x 3 , y 2 = x 2 + 2 x 3 , y 3 = x 3 → { x 1 = y 1 − y 2 + y 3 x 2 = y 2 − 2 y 3 x 3 = y 3 \begin{cases} y_1=x_1+&x_2+&x_3,\\ y_2= &x_2+2&x_3,\\ y_3=&&x_3\\ \end{cases}\to\begin{cases} x_1=y_1-y_2+y_3\\ x_2=y_2-2y_3\\ x_3=y_3\\ \end{cases} y1=x1+y2=y3=x2+x2+2x3,x3,x3 x1=y1y2+y3x2=y22y3x3=y3
f = y 1 2 + y 2 2 f=y_1^2+y_2^2 f=y12+y22,所用变换矩阵为:
C = ( 1 − 1 1 0 1 − 2 0 0 1 ) C=\begin{pmatrix}1&-1&1\\ 0 & 1&-2\\ 0 & 0 &1\\ \end{pmatrix} C= 100110121
例16 化二次型 f = 2 x 1 x 2 + 2 x 1 x 3 − 6 x 2 x 3 f=2x_1x_2+2x_1x_3-6x_2x_3 f=2x1x2+2x1x36x2x3称规范型,并求所用的变换矩阵。

:在f中不包含平方项,由于含有 x 1 x 2 x_1x_2 x1x2乘积项,故令
{ x 1 = y 1 + y 2 , x 2 = y 1 − y 2 , x 3 = y 3 , \begin{cases} x_1=y_1+y_2,\\ x_2=y_1-y_2,\\ x_3=y_3, \end{cases} x1=y1+y2,x2=y1y2,x3=y3,
带入可得: f = 2 y 1 2 − 2 y 2 2 + 2 y 1 y 3 + 2 y 2 y 3 − 6 y 1 y 3 + 6 y 2 y 3 f=2y_1^2-2y_2^2+2y_1y_3+2y_2y_3-6y_1y_3+6y_2y_3 f=2y122y22+2y1y3+2y2y36y1y3+6y2y3
= 2 y 1 2 − 2 y 2 2 − 4 y 1 y 3 + 8 y 2 y 3 2y_1^2-2y_2^2-4y_1y_3+8y_2y_3 2y122y224y1y3+8y2y3
= 2 ( y 1 − y 3 ) 2 − 2 ( y 2 − 2 y 3 ) 2 + 6 y 3 2 2(y_1-y_3)^2-2(y_2-2y_3)^2+6y_3^2 2(y1y3)22(y22y3)2+6y32

令 { z 1 = 2 ( y 1 − y 3 ) , z 2 = 2 ( y 2 − 2 y 3 ) , z 3 = 6 y 3 , 即 { y 1 = 1 2 z 1 + 1 6 z 3 , y 2 = 1 2 z 2 + 2 6 z 3 y 3 = 1 6 z 3 令\begin{cases}z_1=\sqrt 2(y_1-y_3), \\ z_2 =\sqrt 2(y_2-2y_3),\\ z_3=\sqrt 6 y_3,\\ \end{cases}即 \begin{cases} y_1=\frac 1 {\sqrt 2}z_1 + \frac 1 {\sqrt 6}z_3,\\ y_2=\frac 1 {\sqrt 2}z_2+\frac 2 {\sqrt 6}z_3\\ y_3 = \frac 1 {\sqrt 6}z_3 \end{cases}\\ z1=2 (y1y3),z2=2 (y22y3),z3=6 y3, y1=2 1z1+6 1z3,y2=2 1z2+6 2z3y3=6 1z3就把f化成规范行 f = z 1 2 − z 2 2 + z 3 2 f=z_1^2-z_2^2+z_3^2 f=z12z22+z32
在这里插入图片描述
C = ( 1 1 0 1 − 1 0 0 0 1 ) ( 1 2 0 1 6 0 1 2 2 6 0 0 1 6 ) = ( 1 2 1 2 3 6 1 2 − 1 2 − 1 6 0 0 1 6 ) C=\begin{pmatrix} 1&1&0\\ 1&-1&0\\ 0&0&1\\ \end{pmatrix} \begin{pmatrix} \frac 1 {\sqrt 2}&0&\frac 1 {\sqrt 6}\\ 0&\frac 1 {\sqrt 2}&\frac 2 {\sqrt 6}\\ 0&0&\frac 1 {\sqrt 6}\\ \end{pmatrix}= \begin{pmatrix} \frac 1 {\sqrt 2}&\frac 1 {\sqrt 2}&\frac 3 {\sqrt 6}\\ \frac 1 {\sqrt 2}&-\frac 1 {\sqrt 2}&-\frac 1 {\sqrt 6}\\ 0&0&\frac 1 {\sqrt 6}\\ \end{pmatrix} C= 110110001 2 10002 106 16 26 1 = 2 12 102 12 106 36 16 1

7 正定二次型

定理9(惯性定理) 设有二次型 f = x T A x f=x^TAx f=xTAx,它的秩为r,有两个可逆变换
x = C y 及 x = P z x=Cy及x=Pz x=Cyx=Pz
使 f = k 1 y 1 2 + k 2 y 2 2 + ⋯ + k r y r 2 , ( k i ≠ 0 ) f=k_1y_1^2+k_2y_2^2+\cdots +k_ry_r^2,(k_i\neq 0) f=k1y12+k2y22++kryr2,(ki=0)
使 f = λ 1 z 1 2 + λ 2 z 2 2 + ⋯ + λ r z r 2 , ( λ i ≠ 0 ) f=\lambda_1z_1^2+\lambda_2z_2^2+\cdots +\lambda_rz_r^2,(\lambda_i \neq 0) f=λ1z12+λ2z22++λrzr2,(λi=0)
k 1 , ⋯   , k r 中正数的个数与 λ 1 , λ 2 , ⋯ λ r k_1,\cdots ,k_r中正数的个数与\lambda_1,\lambda_2,\cdots \lambda_r k1,,kr中正数的个数与λ1,λ2,λr中正数的个数相等。
二次型的标准形中正系数的个数称为二次型的正惯性指数,负系数的个数称为负惯性指数。若二次型f的正惯性指数为p,秩为r,则f的规范型便可确定为
f = y 1 2 + ⋯ + y p 2 − y p + 1 2 − ⋯ y r 2 f=y_1^2+\cdots+y_p^2-y_{p+1}^2-\cdots y_r^2 f=y12++yp2yp+12yr2

定义10:设有二次型f(x)= x T A x ,如果对任何 x ≠ 0 ,都有 f ( x ) > 0 ,则称 f 正定二次型,并称对称阵 A 是正定的;如果对任何 x ≠ 0 ,都有 f ( x ) < 0 ,则称 f 为负定二次型,并称对称阵 A 是负定的 x^TAx,如果对任何x \neq 0,都有f(x)>0,则称f正定二次型,并称对称阵A是正定的;如果对任何x \neq 0,都有f(x)<0,则称f为负定二次型,并称对称阵A是负定的 xTAx,如果对任何x=0,都有f(x)>0,则称f正定二次型,并称对称阵A是正定的;如果对任何x=0,都有f(x)<0,则称f为负定二次型,并称对称阵A是负定的
定理10 n元二次型f= x T A x x^TAx xTAx为正定的重复必要条件是:它的标准型的n个系统全为正,即它的规范型的n个系数全为1,亦即它的正惯性指数等于n。
推论 对称阵A为正定的重复必要条件是:A的特征值全为正。

定理11(赫尔维茨定理) 对称阵A为正定的充分必要条件是:A的各阶主子式都为正,即: a 11 > 0 , ∣ a 11 a 12 a 21 a 22 ∣ > 0 ⋯ ∣ a 11 ⋯ a 1 n ⋮ ⋮ ⋮ a n 1 ⋯ a n n ∣ > 0 对称阵负定的充分必要条件是:奇数阶主子式为负,而偶数阶主子式为正, a_{11}>0,\begin{vmatrix}a_{11}&a_{12}\\ a_{21}&a_{22}\end{vmatrix}>0 \cdots \begin{vmatrix}a_{11}&\cdots& a_{1n}\\ \vdots&\vdots&\vdots\\ a_{n1}&\cdots&a_{nn}\end{vmatrix}>0 对称阵负定的充分必要条件是:奇数阶主子式为负,而偶数阶主子式为正, a11>0, a11a21a12a22 >0 a11an1a1nann >0对称阵负定的充分必要条件是:奇数阶主子式为负,而偶数阶主子式为正,
( − 1 ) r ∣ a 11 ⋯ a 1 r ⋮ ⋮ a r 1 ⋯ a r r ∣ > 0 ( r = 1 , 2 , 3 ⋯ n ) (-1)^r\begin{vmatrix} a_{11}&\cdots&a_{1r}\\ \vdots&&\vdots\\ a_{r1}&\cdots&a_{rr}\\ \end{vmatrix}>0(r=1,2,3\cdots n) (1)r a11ar1a1rarr >0(r=1,2,3n)

扩展阅读

算法为骨,CAD为魂
亲士工具箱:支持中望CAD2024、AutoCad2013及以上,多年承接CAD项目的精华
工作中遇到的问题,可以按类别查阅鄙人的算法文章,请点击《算法与数据汇总》。
学习算法:按章节学习《喜缺全书算法册》,大量的题目和测试用例,打包下载。重视操作
活到老,学到老。明朝中后期,大约50%的进士能当上堂官(副部及更高);能当上堂官的举人只有十余人。
子墨子言之:事无终始,无务多业。也就是我们常说的专业的人做专业的事。

视频课程

先学简单的课程,请移步CSDN学院,听白银讲师(也就是鄙人)的讲解。
https://edu.csdn.net/course/detail/38771
如何你想快速形成战斗了,为老板分忧,请学习C#入职培训、C++入职培训等课程
https://edu.csdn.net/lecturer/6176

测试环境

操作系统:win7 开发环境: VS2019 C++17
或者 操作系统:win10 开发环境: VS2022 C++17
如无特殊说明,本算法用**C++**实现。

评论 69
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

软件架构师何志丹

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值