小知识,大挑战!本文正在参与“程序员必备小知识”创作活动。
前言
Hello!小伙伴!
非常感谢您阅读海轰的文章,倘若文中有错误的地方,欢迎您指出~
自我介绍 ଘ(੭ˊᵕˋ)੭
昵称:海轰
标签:程序猿|C++选手|学生
简介:因C语言结识编程,随后转入计算机专业,有幸拿过一些国奖、省奖...已保研。目前正在学习C++/Linux/Python
学习经验:扎实基础 + 多做笔记 + 多敲代码 + 多思考 + 学好英语!
机器学习小白阶段
文章仅作为自己的学习笔记 用于知识体系建立以及复习
知其然 知其所以然!
5.5 二次型及其标准形
定义8:二次型
含有n个变量x1,x2,...,xn的二次齐次函数
f(x1,x2,...,xn)=a11x12+a22x22+....+annxn2+2a12x1x2+2a13x1x3+....+2an−1,nxn−1xn
称为二次型
f中的每一项的次数都是2称为二次齐次函数 比如x12,x1x2....
取aji=aij,则有
2aijxixj=aijxixj+ajixjxi
所以
f=a11x12+a12x1x2+...+a1nx1xn+a21x2x1+a22x22+...+a2nx2xn+...+an1xnx1+an2xnx2+....+annxn2=∑i,j=1naijxixj
对于二次型,寻求可逆的线性变换
⎩⎨⎧x1=c11y1+c12y2+....+c1nynx2=c21y1+c22y2+....+c2nyn....xn=cn1y1+cn2y2+....+cnnyn
也就是将上式代入f,替换xi,使得二次型只含有平方项,得到
f=k1y12+k2y22+....+knyn2
这种只含有平方项的二次型,称为二次型的标准形(或法式)
如果标准型的系数k1,k2,...,kn只在1,−1,0三个数中取值,即
f=y12+...+yp2−yp+12−...−yr2
则称上式为二次型的规范形
f=a11x12+a12x1x2+...+a1nx1xn+a21x2x1+a22x22+...+a2nx2xn+...+an1xnx1+an2xnx2+....+annxn2=∑i,j=1naijxixj
可以变形为
f=x1(a11x1+a12x2+...+a1nxn)+x2(a21x1+a22x2+...+a2nxn)+....+xn(an1x1+an2x2+...+annxn)=(x1,x2,...,xn)⎣⎡a11x1+a12x2+...+a1nxna21x1+a22x2+...+a2nxn...an1x1+an2x2+...+annxn⎦⎤=(x1,x2,...,xn)⎣⎡a11a21...an1a12a22...an2.........a1na2n...ann⎦⎤⎣⎡x1x2...xn⎦⎤
记
A=⎣⎡a11a21...an1a12a22...an2.........a1na2n...ann⎦⎤,x=⎣⎡x1x2...xn⎦⎤
则二次型可记作
f=xTAx
其中A为对称阵
例如,二次型f=x2−3z2−4xy+yz用矩阵记号写出来,就是
f=(x,y,z)⎣⎡1−20−1021021−3⎦⎤⎣⎡xyz⎦⎤
对称阵A的求法
注意xx、yy、zz...系数就是多项式中相对应的系数,而xy,xz,yz,...的系数就是多项式中相对应的系数的一半
任给一个二次型,就可以惟一确定一个对称阵;反之认给一个对称阵,也可以惟一地确定一个二次型。说明二次阵与对称阵之间存在一一对应的关系
把对称阵A叫做二次型f的矩阵,也把f叫做对称阵A的二次型,对称阵A的秩就叫做二次型f的秩
在式子中
⎩⎨⎧x1=c11y1+c12y2+....+c1nynx2=c21y1+c22y2+....+c2nyn....xn=cn1y1+cn2y2+....+cnnyn
记C=(cij),则上式可以变为
x=Cy
代入f=xTAx,得到
f=xTAx=(Cy)TA(Cy)=yT(CTAC)y
定义9:合同
设A和B是n阶矩阵,若存在可逆矩阵C,使得B=CTAC,则称矩阵A与B合同
若A是对称阵
BT=(CTAC)T=CTATC=CTAC=B
得到B也是对称阵
在B=CTAC中,因为C、CT可逆 所以
R(B)=R(A)
由此可知,经可逆变换x=Cy后,二次型f的矩阵由A变为与A合同的矩阵CTAC,且二次型的秩不变
如果要使二次型f经可逆变换x=Cy变成标准形,也就是
f=yT(CTAC)y=k1y12+k2y22+...+knyn2=(y1,y2,...,yn)⎣⎡k1k2.kn⎦⎤⎣⎡y1y2.yn⎦⎤
也就是要使得CTAC成为对角阵
所以就是寻找一个可逆矩阵C 使得CTAC变为对角阵,这个过程就称为对称阵A合同对角化
定理8
由定理7得,任一对称阵A,总存在正交阵P,使得P−1AP=PTAP=Λ
将此结论运用到二次型得到
任一二次型f=∑i,j=1naijxixj(aij=aji),总有正交变换x=Py,使f化为标准
形
f=λ1y12+λ2y22+....+λnyn2
其中λ1,λ2,...,λn是f的矩阵A=(aij)的特征值
推论
任给n元二次型f(x)=xTAx(AT=A),总有可逆变换x=Cz,使得f(Cz)为规范形
举例
例14
求一个正交变换x=Py,把二次型
f=−2x1x2+2x1x3+2x2x3
化为标准形
思路:
- 需要将二次型变为标准形,其实就是需要找到一个可逆矩阵C,使得CTAC变成一个对角阵Λ
- 又因为正交阵P,可以使得P−1AP=PTAP=Λ
- 所以这道题实质是求一个正交阵P
解答:
二次型的矩阵为
A=⎣⎡0−11−101110⎦⎤
再求一个正交阵P,使得P−1AP=PTAP=Λ
由
∣A−λE∣=A=⎣⎡−λ−11−1−λ111−λ⎦⎤=−(λ−12)(λ+1)
解得A的特征值为
λ1=−2,λ2=λ3=1
对应λ1=−2,解方程(A+2E)x=0
A+2E=⎣⎡2−11−121112⎦⎤∼⎣⎡100010110⎦⎤
得基础解系ζ1=⎣⎡−1−11⎦⎤
对ζ1进行单位化,得
p1=31⎣⎡−1−11⎦⎤
对应λ2=λ3=1,解方程(A−E)x=0
∣A−E∣=⎣⎡−1−11−1−1111−1⎦⎤∼⎣⎡100100−100⎦⎤
得基础解系
ζ2=⎣⎡−110⎦⎤,ζ3=⎣⎡101⎦⎤
将ζ2、ζ2正交化
令η2=ζ2
η3=ζ3−∣∣ζ2∣∣[η2,ζ3]η2=⎣⎡101⎦⎤+21⎣⎡−110⎦⎤=21⎣⎡112⎦⎤
再将η2、ζ3单位化,得
p2=21⎣⎡−110⎦⎤,p3=61⎣⎡112⎦⎤
将p1,p2,p3构成正交矩阵
P=(p1,p2,p3)=⎣⎡−31−3131−21210616162⎦⎤
有
P−1AP=PTAP=Λ=⎣⎡−200010001⎦⎤
综上,有正交变换x=Py
⎣⎡x1x2x3⎦⎤=⎣⎡−31−3131−21210616162⎦⎤⎣⎡y1y2y3⎦⎤
结语
说明:
- 参考于 课本《线性代数》第五版 同济大学数学系编
- 配合书中概念讲解 结合了自己的一些理解及思考
文章仅作为学习笔记,记录从0到1的一个过程
希望对您有所帮助,如有错误欢迎小伙伴指正~
我是 海轰ଘ(੭ˊᵕˋ)੭
如果您觉得写得可以的话,请点个赞吧
谢谢支持 ❤️
