本文涉及知识点
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,a2⋯ar线性无关。
证: 设有
λ
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左乘上式两端,因当i≥2时,a1Tai=0故得:λiaiTai=0,故λ1a1Ta1=0,因为ai=0,故a1Ta1=∣∣a1∣∣2=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=
1−21
正交,试求一个非零向量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,e2⋯er是向量空间V
(
V
⊂
R
n
)
(V\sub \R^n)
(V⊂Rn)的一个基,如果
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=
12−1
,a2=
−131
,a3=
4−10
试用施密特正交化过程把这组向量规范正交化。
解:
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
12−1
=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}
=
4−10
−31
12−1
+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,即A−1=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是单位向量ai和aj正交
一,若A为正交矩阵,则
A
−
1
=
A
T
A^{-1}=A^T
A−1=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∣=1,∣ATA∣=∣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−λa21⋮a31a12a22−λ⋮a32⋯⋯⋯a1na2n⋮a3n−λ
上式是以
λ
\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=0即a1λ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)+⋯+ak−1pk−1(λk−1−λ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
P−1AP=B则称B是A的相似矩阵,或说矩阵A与B相似。对A进行运算
P
−
1
A
P
=
B
P^{-1}AP=B
P−1AP=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
=
P−1AP−P−1λEP
=
P−1(A−λE)P
=
∣
P
−
1
∣
∣
A
−
λ
E
∣
∣
P
∣
=
∣
A
−
λ
E
∣
=|P^{-1}||A-\lambda E||P|=|A-\lambda E|
=∣P−1∣∣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阶对角矩阵,A是n阶矩阵,假定存在可逆矩阵P,使得P−1AP=Λ。式子左乘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)λi即Piλ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
p1与p2正交。
证:
λ
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=0即p1与p2正交。
定理7 设A为n阶对称阵,则必有正交阵P,使得
P
−
1
A
P
=
P
T
A
P
=
Λ
,
其中
Λ
是以
A
的
n
个特征值为对角元的对角阵
P^{-1}AP=P^TAP=\Lambda,其中\Lambda是以A的n个特征值为对角元的对角阵
P−1AP=PTAP=Λ,其中Λ是以A的n个特征值为对角元的对角阵。
推论 设A为n阶段对称阵, λ 是 A 的特征方程的重根,则矩阵 A − λ E 的秩 = n − k , 从而对应的特征值 λ 恰有 k 个线性无关的特征向量 \lambda是A的特征方程的重根,则矩阵A-\lambda E的秩=n-k,从而对应的特征值\lambda 恰有k个线性无关的特征向量 λ是A的特征方程的重根,则矩阵A−λE的秩=n−k,从而对应的特征值λ恰有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中列向量的排列次序相对应。
P−1AP=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+⋯+2an−1,nxn−1xn称为二次型。式子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=1∑naijxixj式子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,k2⋯kn只有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=1∑naijxixy,总有正交变换
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=
0−11−101110
此矩阵即例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=
−31−3131−21210616162
使得
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=21z1,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)2−x22−x32−2x2x3+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=y1−y2+y3x2=y2−2y3x3=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=
100−1101−21
例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+2x1x3−6x2x3称规范型,并求所用的变换矩阵。
解:在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=y1−y2,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=2y12−2y22+2y1y3+2y2y3−6y1y3+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
2y12−2y22−4y1y3+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(y1−y3)2−2(y2−2y3)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(y1−y3),z2=2(y2−2y3),z3=6y3,即⎩
⎨
⎧y1=21z1+61z3,y2=21z2+62z3y3=61z3就把f化成规范行
f
=
z
1
2
−
z
2
2
+
z
3
2
f=z_1^2-z_2^2+z_3^2
f=z12−z22+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=
1101−10001
21000210616261
=
2121021−21063−6161
7 正定二次型
定理9(惯性定理) 设有二次型
f
=
x
T
A
x
f=x^TAx
f=xTAx,它的秩为r,有两个可逆变换
x
=
C
y
及
x
=
P
z
x=Cy及x=Pz
x=Cy及x=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+⋯+yp2−yp+12−⋯yr2
定义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⋯
a11⋮an1⋯⋮⋯a1n⋮ann
>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
a11⋮ar1⋯⋯a1r⋮arr
>0(r=1,2,3⋯n)

扩展阅读
| 算法为骨,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++**实现。

1271

被折叠的 条评论
为什么被折叠?



