本文已参与「新人创作礼」活动,一起开启掘金创作之路。
矩阵的初等变换
一、初等变换
1. 初等变换的定义
下面三种变换称为矩阵的初等行变换
-
对调两行(对调i,j两行,记作ri↔rj
-
以数(k=0)乘某一行中的所有元素(第i行乘k,记作ri×k)
-
把某一行所有元素的k倍加到另一行对应的元素上去(第j行的k倍加到第i行,记作ri+krj)
把定义中的“行”换成“列”,即得矩阵的初等列变换的定义(所用记号是把“r“换成”c“)。矩阵的初等行变换与初等列变换,统称初等变换
2. 行最简形矩阵的定义
非零行的第一个非零元为1,且这些非零元所在的列的其他元素都为0
例如:⎝⎛1000010000102300⎠⎞
二、矩阵等价
1. 矩阵等价的定义
如果矩阵A经过有限次初等行变换变成矩阵B,就称矩阵A与矩阵B行等价,记作A∼rB;如果矩阵A经过有限次初等列变换变成矩阵B,就称矩阵A与矩阵B列等价,记作A∼cB;如果A经过有限次初等变换变成矩阵B,就称矩阵A与矩阵B等价,记作A∼B。(A,B是同型矩阵)
2. 矩阵等价的性质
矩阵之间的等价关系具有以下性质
-
反身性:A∼A
-
对称性:若A∼B,则B∼A
-
传递性:若A∼B,B∼C,则A∼C
3. 矩阵等价的定理
设A与B为m×n矩阵,那么:
-
A∼rB的充分必要条件是存在m阶可逆矩阵P,使PA=B
-
A∼cB的充分必要条件是存在n阶可逆矩阵Q,使AQ=B
-
A∼B的充分必要条件是存在m阶可逆矩阵及n阶可逆矩阵Q,使PAQ=B
三、初等矩阵
1. 初等矩阵的定义
由单位阵E经过一次初等变换得到的矩阵称为初等矩阵
2. 三种初等矩阵
-
把单位阵中的第i,j两行(或第i,j两列)对调,得初等矩阵
-
以数k=0乘单位阵的第i行(或第i列),得初等矩阵
-
以k乘E的第j行加到第i行上或以k乘E的第i列加到第j列上,得初等矩阵
3. 初等矩阵的性质
-
设A是一个m×n矩阵,对A施行一次初等行变换,相当于在A的左边乘以相应的m阶初等矩阵;对A施行一次初等列变换,相当于在A的右边乘以相应的n阶初等矩阵
-
方阵A可逆的充分必要条件是,存在有限个初等矩阵P1,P2,⋯,Pl,使A=P1P2⋯Pl
例1:设A=⎝⎛03−2−2031−20⎠⎞,证明A可逆,并求A−1
思路:构建矩阵(A∣E),将A的部分通过初等行变换变成E,同时E变成A−1
(A∣E)=⎝⎛03−2−2031−201000100 0 1⎠⎞E的部分变成上三角矩阵→⎝⎛3000−20−211019104006⎠⎞将非主对角线上的数变成0⎝⎛3000−2000118−899−4412−66⎠⎞主对角线上的值变为1⎝⎛100010001649324436⎠⎞=(E∣A−1)
故A−1=⎝⎛649324436⎠⎞
矩阵的秩
一、k阶子式的定义
在m×n矩阵A中,任取k行与k列(k≤m,k≤n),位于这些行列交叉处的k2个元素,不改变它们在A中所处的位置次序而得的k阶行列式,称为矩阵A的k阶子式
二、矩阵的秩的定义
设在矩阵A中有一个不等于0的r阶子式D,且所有r+1阶子式(如果存在的话)全等于0,那么D称为矩阵A的最高阶非零子式,数r称为矩阵A的秩,记作R(A)(也可写作r(A))。并规定零矩阵的秩等于0
三、矩阵秩的性质
-
0≤R(Am×n)≤min{m,n}
-
R(AT)=R(A)
-
R(AB)≤min{R(A),R(B)}
-
若P,Q可逆,则R(PAQ)=R(A)
证明:
r(PAQ)≤r(AQ)
r(AQ)=r(EAQ)=r(P−1PAQ)≤r(PAQ)
∴r(PAQ)=r(AQ)
r(AQ)≤r(A)
r(A)=r(AE)=r(AQQ−1)≤r(AQ)
∴r(A)=r(AQ)=r(PAQ)
-
max{R(A),R(B)}≤R(A,B)≤R(A)+R(B),特别的,当B=b为非零向量时,有R(A)≤R(A,b)≤R(A)+1
-
R(A+B)≤R(A)+R(B)
-
若A∼B,则R(A)=R(B)
-
若Am×nBn×l=O,则R(A)+R(B)≤n
求矩阵的秩,可以用定义求,建议将矩阵化为阶梯形矩阵,阶梯形矩阵的秩为行列式的行值
例1:设A=⎝⎛12−2−1−2−44−628−20−103−6⎠⎞,b=⎝⎛1234⎠⎞,求矩阵A及矩阵B=(A,b)的秩
B=[A∣b]=⎝⎛12−2−1−2−44−628−20−103−61234⎠⎞→⎝⎛1000−200−82422−121−71055⎠⎞→⎝⎛1000−2−8002224−1−7121550⎠⎞→⎝⎛1000−2−8002220−1−7101551⎠⎞
∴r(A)=3,r(B)=4
线性方程组的解
一、线性方程组的定义
设有n个未知数m个方程的线性方程组⎩⎨⎧a11x1+a12x2+⋯+a1nxn=b1a21x1+a22x2+⋯+a2nxn=b2⋯am1x1+am2x2+⋯+amnxn=bm,可以写成以向量x为未知元的向量方程Ax=b,则n元线性方程Ax=b
-
无解的充分必要条件是R(A)<R(A,b)
-
有唯一解得充分必要条件是R(A)=R(A,b)=n
-
有无限多解得充分必要条件是R(A)=R(A,b)<n
n是未知数的个数,也是系数矩阵A的列数
二、步骤
-
对[A∣b]进行初等行变换,化为行最简形
-
令自由未知量分别取1,其余取0,得到非自由未知量的值。得到基础解系(n−r(A)),即齐次方程通解k1ξ1+⋯+kn−rξn−r
-
令自由未知量取0,得到非自由未知量,即特解η,故非齐次方程通解y+k1η1+⋯+kn−rηn−r
首先要先判断有没有解!
齐次方程组没有第三步。做了也是O矩阵
例1:求其次线性方程组⎩⎨⎧x1+2x2+2x3+x4=02x1+x2−2x3−2x4x1−x2−4x3−3x4
A=⎝⎛12121−12−2−41−2−3⎠⎞→⎝⎛100010−220−35340⎠⎞
n−r=4−2=2
(此处根据令自由未知量分别取1,其余取0,得到非自由未知量的值。得到基础解系(n−r(A)),即齐次方程通解k1ξ1+⋯+kn−rξn−r一步一步演示)
先设ξ1=⎝⎛10⎠⎞,ξ2=⎝⎛ 01⎠⎞
(注意,无论是非齐次还是齐次线性方程组,此处都不需要使用增广矩阵,只使用系数矩阵A)
将ξ1=⎝⎛10⎠⎞分别代入矩阵
第一行得x1的系数为2
第二行得x2的系数为−2
则ξ1=⎝⎛2−210⎠⎞
同理ξ2=⎝⎛−35−3401⎠⎞
(齐次线性方程组没有最后一步)
故通解为x=⎝⎛x1x2x3x4⎠⎞=k1ξ1+k2ξ2=k1⎝⎛2−210⎠⎞+k2⎝⎛35−3401⎠⎞(k1,k2为任意常数)
例2:设有线性方程组⎩⎨⎧(1+λ)x1+x2+x3=0x1+(1+λ)x2+x3=3x1+x2+(1+λ)x3=λ,讨论方程组解的情况,并求无限解时的通解
法1:
(A∣b)=⎝⎛1001λ01+x−λ−λ(3+λ)λ3−λ(1−λ)(3+λ)⎠⎞
-
当λ=0且λ=−3时,r(A)=3=r(A∣b)=3,有唯一解
-
当λ=0时,r(A)=1=r(A∣b)=2,无解
-
当λ=−3时,r(A)=r(A∣b)=2<3,无限多解
当λ=−3时
(A∣b)=⎝⎛100010−1−10−1−20⎠⎞
n−r=3−2=1
∴ξ=⎝⎛111⎠⎞,y=⎝⎛−1−20⎠⎞
∴通解为y+kξ(k为任意常数)
法2:
观察到系数矩阵A为方阵
有唯一解即r(A)=r(A∣b)=3⇔∣A∣=0
有无限多解r(A)=r(A∣b)<3⇔∣A∣=0
∣A∣=∣∣1+λ1111+λ1111+λ∣∣=(3+λ)λ2=0
解得λ=−3或λ=0
(此处只能得到唯一解得解集,其余的需要代入讨论)
(A∣b)=⎝⎛111111111030⎠⎞=⎝⎛100100100020⎠⎞
r(A)=r(A∣b)无解
(A∣b)=⎝⎛−2111−2111−20−3−3⎠⎞→⎝⎛100010−1−10−1−20⎠⎞
n−r=3−2=1
ξ=⎝⎛111⎠⎞,y=⎝⎛−1−20⎠⎞
∴通解为y+kξ(k为任意常数)