Typesetting math: 85%

【数学】矩阵论学习笔记(一)

符号约定

  • II:单位矩阵

  • AHAHAA共轭转置

  • rankArankAAA

  • cond(A)cond(A)AA条件数

  • Cm×nrCm×nr秩为 rr 的复 m×nm×n 矩阵 集合

第一章 矩阵的相似变换

1.1 矩阵的特征值与特征向量

定义 11(特征值、特征向量)

  • ACn×nACn×n,若 λCλCxCn,x0xCn,x0 使得

    Ax=λxAx=λx

    则称 λλAA特征值xx 称为 AA 的属于 λλ特征向量.

定义 22

ACn×nACn×n

  • λInAλInAAA特征矩阵.

  • det(λInA)=|λInA|det(λInA)=|λInA|AA特征多项式.

  • |λInA|=0|λInA|=0AA特征方程.

📝 Notations

  1. AA 的特征值就是 AA 的特征方程的根;

  2. nn 阶方阵 AA 在复数范围内一定有 nn 个特征值.

定义 33(矩阵的迹)

  • A=(aij)n×nCn×nA=(aij)n×nCn×n,称 a11+a22++anna11+a22++annAA,记为 trAtrA.

定理 11

  • nn 阶方阵 A=(aij)n×nA=(aij)n×n 的特征值为 λ1,λ2,,λnλ1,λ2,,λn,则

    1. λ1+λ2,+λn=a11+a22++ann=trAλ1+λ2,+λn=a11+a22++ann=trA.

    2. nk=1λk=|A|nk=1λk=|A| AA 可逆的充要条件是 AA 的特征值非零.

    3. ATAT 的特征值是 λ1,λ2,,λnλ1,λ2,,λn

    4. AHAH 的特征值是 ¯λ1,¯λ2,,¯λn¯¯¯¯¯λ1,¯¯¯¯¯λ2,,¯¯¯¯¯¯λn

定理 22

  • λiλiACn×nACn×nriri 重特征值(称 riri 为特征值 λiλi 的代数重数),对应 λiλisisi 个线性无关的特征向量(称 sisi 为特征值 λiλi 的几何重数),则 1siri1siri.

定义 44

  • f(λ)f(λ)λλ 的多项式

    f(λ)=asλs+as1λs1++a1λ+a0f(λ)=asλs+as1λs1++a1λ+a0

    对于 ACn×nACn×n,规定

    f(A)=asAs+as1As1++a1A+a0If(A)=asAs+as1As1++a1A+a0I

    f(A)f(A) 为矩阵 AA 的多项式.

定理 33

  • ACn×nACn×nAAnn 个特征值为 λ1,λ2,,λnλ1,λ2,,λn,对应的特征向量分别为 x1,x2,,xnx1,x2,,xn,又设 f(λ)f(λ) 为一多项式,则 f(A)f(A) 的特征值为 f(λ1),f(λ2),,f(λn)f(λ1),f(λ2),,f(λn),对应的特征向量分别仍为 x1,x2,,xnx1,x2,,xn.

  • 特别地f(A)=0f(A)=0 时,AA 的任意特征值 λiλi 都满足 f(λi)=0f(λi)=0.

定理 44

  • λ1,λ2,,λsλ1,λ2,,λs 是方阵 AA 的互不相同的特征值,x1,x2,,xsx1,x2,,xs 是分别与之对应的特征向量,则 x1,x2,,xsx1,x2,,xs 线性无关.

  • 特别地,属于实对称阵的不同特征值的特征向量是正交的.

1.2 矩阵的相似对角化

定义 11

  • A,BCn×nA,BCn×n,若 PCn×nPCn×n 使得 P1AP=BP1AP=B,则称 AABB 相似,记为 ABAB;并称可逆阵 PP 为把 AA 变成 BB相似变换矩阵.

  • AA 与对角阵相似时,称 AA 可对角化.

定理 11(性质)

  • A,B,CCn×nA,B,CCn×n,则

    1. 反身性:AAAA

    2. 对称性:若 ABAB,则 BABA.

    3. 传递性:若 AB,BCAB,BC,则 ACAC.

  • 相似关系为等价关系.

定理 22

  • A,BCn×n,ABA,BCn×n,ABf(λ)f(λ) 是一多项式,则

    1. rank(A)=rank(B)rank(A)=rank(B)|A|=|B||A|=|B|

    2. det(λIA)=det(λIB)det(λIA)=det(λIB),即 AABB 有相同的特征多项式,从而有相同的特征值.

    3. A1B1,ATBT,f(A)f(B)A1B1,ATBT,f(A)f(B).

定理 33

  • ACn×nACn×n,则 AA 可对角化的充分必要条件是 AAnn 个线性无关的特征向量 ξ1,ξ2,,ξnξ1,ξ2,,ξn.

  • 推论:设 ACn×nACn×n,如果 AAnn 个不同的特征值,则 AA 可对角化.

定理 44

  • λ1,λ2,,λsλ1,λ2,,λsnn 阶方阵 AA 的所有互不相同的特征值,其重数分别为 r1,r2,,rsr1,r2,,rs. 则 AA 可对角化的充要条件是 riri 重特征值 λiλiri(i=1,,s)ri(i=1,,s) 个线性无关的特征向量.

求解一阶线性常系数微分方程组

1.3 矩阵的 Jordan 标准形

定义 11

  • 形如

    Ji=(λi1λi1λi)ri×riJi=⎜ ⎜ ⎜ ⎜ ⎜ ⎜λi1λi1λi⎟ ⎟ ⎟ ⎟ ⎟ ⎟ri×ri

    的矩阵称为 ririJordan 块.

  • 由若干个 Jordan 块构成的分块对角阵

    J=(J1J2Js)J=⎜ ⎜ ⎜ ⎜ ⎜J1J2Js⎟ ⎟ ⎟ ⎟ ⎟

    称为 Jordan 矩阵.

定理 11

  • 对于 riri 阶 Jordan 块

    Ji=(λi1λi1λi)ri×riJi=⎜ ⎜ ⎜ ⎜ ⎜ ⎜λi1λi1λi⎟ ⎟ ⎟ ⎟ ⎟ ⎟ri×ri

    Jki=(λkiC1kλk1iC2kλk2iCri1kλkri+1iλkiC1kλk1iCri2kλkri+2iC1kλk1iλki)=(λk11!(λk)12!(λk)1(ri1)!(λk)(ri1)λk11!(λk)1(ri2)!(λk)(ri2)11!(λk)λk)|λ=λiJki=⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜λkiC1kλk1iC2kλk2iCri1kλkri+1iλkiC1kλk1iCri2kλkri+2iC1kλk1iλki⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟=⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜λk11!(λk)12!(λk)1(ri1)!(λk)(ri1)λk11!(λk)1(ri2)!(λk)(ri2)11!(λk)λk⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟∣ ∣ ∣ ∣ ∣ ∣ ∣ ∣ ∣ ∣λ=λi

对于 Jordan 阵

J=(J1J2Js)J=⎜ ⎜ ⎜ ⎜ ⎜J1J2Js⎟ ⎟ ⎟ ⎟ ⎟

Jk=(Jk1Jk2Jks)Jk=⎜ ⎜ ⎜ ⎜ ⎜ ⎜Jk1Jk2Jks⎟ ⎟ ⎟ ⎟ ⎟ ⎟

定理 22

  • ACn×nACn×n,则 AA 与一个 Jordan 矩阵 JJ 相似,即存在可逆矩阵 PCn×nnPCn×nn,使得

    P1AP=JP1AP=J

    这个 Jordan 矩阵 JJ 除 Jordan 块的排列次序外由 AA 唯一确定,并称 JJAAJordan 标准形.

λλ 矩阵的初等变换包括:

  1. 交换两行(列);

  2. 某一行(列)同乘以一个非零常数;

  3. 某一行(列)同乘以多项式 φ(λ)φ(λ) 加到另一行(列).

定理 33

  • ACn×nACn×n,则 λInAλInA 经过初等变换可化为如下 Smith 标准形

    S(λ)=(d1(λ)d2(λ)dn(λ))S(λ)=⎜ ⎜ ⎜ ⎜ ⎜d1(λ)d2(λ)dn(λ)⎟ ⎟ ⎟ ⎟ ⎟

    其中 di(λ)di(λ) 是首一多项式(最高次项系数为 11 的多项式),i=1,2,,ni=1,2,,n,且 di(λ)di+1(λ),i=1,2,,n1di(λ)di+1(λ),i=1,2,,n1. di(λ)di(λ) 称为 AA不变因子.

初等变换求矩阵 AA 的 Jordan 标准形 的方法:

  1. 用初等变换化特征矩阵 λInAλInA 为 Smith 标准形,求出 AA 的不变因子 d1(λ),d2(λ),,dn(λ)d1(λ),d2(λ),,dn(λ).

  2. AA 的每个次数大于零的不变因子 di(λ)di(λ) 分解为互不相同的一次因式方幂的乘积,这些一次因式的方幂称为 AA初等因子. 设 AA 的全部初等因子为:

    (λλ1)r1,(λλ2)r2,,(λλs)rs(λλ1)r1,(λλ2)r2,,(λλs)rs

    其中 λ1,λ2,,λsλ1,λ2,,λs 可能有相同的值,且 r1+r2++rs=nr1+r2++rs=n.

  3. 写出每个初等因子 (λλi)ri,i=1,2,,s(λλi)ri,i=1,2,,s 对应的 Jordan 块

    Ji=(λi1λi1λi)ri×ri(i=1,2,,s)Ji=⎜ ⎜ ⎜ ⎜ ⎜ ⎜λi1λi1λi⎟ ⎟ ⎟ ⎟ ⎟ ⎟ri×ri(i=1,2,,s)

    以这些 Jodan 块构成的 Jordan 阵

    J=(J1J2Js)J=⎜ ⎜ ⎜ ⎜ ⎜J1J2Js⎟ ⎟ ⎟ ⎟ ⎟

    即为 AA 的 Jordan 标准形.

定义 22

  • ACm×nACm×nλInAλInA 的所有 kk 阶子式 首一最大公因式 Dk(λ)Dk(λ) 称为 AAkk 阶行列式因子k=1,2,,nk=1,2,,n.

定理 44

  • ACm×nACm×nDk(λ)Dk(λ)AAkk 阶行列式因子,dk(λ)dk(λ)AA 的不变因子,k=1,2,,nk=1,2,,n,则

    d1(λ)=D1(λ)dk(λ)=Dk(λ)/Dk1(λ),k=2,3,,n.d1(λ)=D1(λ)dk(λ)=Dk(λ)/Dk1(λ),k=2,3,,n.

行列式因子求矩阵 AA 的 Jordan 标准形 方法:

  1. λInAλInAnn 个行列式因子 D1(λ),D2(λ),,Dn(λ)D1(λ),D2(λ),,Dn(λ)

  2. d1(λ)=D1(λ),dk(λ)=Dk(λ)/Dk1(λ)d1(λ)=D1(λ),dk(λ)=Dk(λ)/Dk1(λ),求 AA 的不变因子.

  3. AA 的初等因子和 Jordan 标准形.

1.4 Hamilton-Cayley 定理

定理 11(Hamilton-Cayley 定理)

  • ACn×nACn×nφ(λ)=det(λInA)φ(λ)=det(λInA),则 φ(A)=Oφ(A)=O.

📝 Notations

  • 利用定理 11 可以简化矩阵运算.

  • 可逆矩阵逆的多项式表示:设 ACn×nACn×n

    φ(λ)=|λInA|=λn+a1λn1++an1λ+anφ(λ)=|λInA|=λn+a1λn1++an1λ+an

    detA0detA0AA 可逆),则

    A1=1an(An1+a1An2++an1I)A1=1an(An1+a1An2++an1I)

定义 11(零化多项式)

  • ACn×nACn×nf(λ)f(λ) 是多项式,如果 f(A)=Of(A)=O,则称 f(λ)f(λ)AA零化多项式.

定义 22(最小多项式)

  • ACn×nACn×n,在 AA 的零化多项式中,次数最低的首一多项式,称为 AA最小多项式,记为 mA(λ)mA(λ).

定理 22

  • ACn×nACn×nφ(λ)=det(λInA)φ(λ)=det(λInA),又设 Dn1(λ)Dn1(λ)λInAλInAn1n1 阶行列式因子,则

    mA(λ)=φ(λ)Dn1(λ)mA(λ)=φ(λ)Dn1(λ)

定理 33

  • ACn×nACn×n,则 AA 的最小多项式 mA(λ)mA(λ) 能整除 AA 的任一零化多项式,且最小多项式是唯一的.

定理 44

  • ACn×nACn×nλ1,λ2,,λsλ1,λ2,,λsAA 的所有互不相同的特征值,则

    mA(λ)=(λλ1)m1(λλ2)m2mA(λ)=(λλ1)m1(λλ2)m2

    其中,mimiAA 的 Jordan 标准形 JJ 中含 λiλi 的 Jordan 块的最高阶数

1.5 向量的内积

定义 11(向量的内积)

  • x=(x1,x2,,xn)TCnx=(x1,x2,,xn)TCny=(y1,y2,,yn)TCny=(y1,y2,,yn)TCn,令

    (x,y)=ni=1xi¯yi=yHx(x,y)=ni=1xi¯¯¯¯yi=yHx

    (x,y)(x,y)xxyy内积.(这里,yH=(¯y1,¯y2,,¯yn)C1×nyH=(¯¯¯¯¯y1,¯¯¯¯¯y2,,¯¯¯¯¯yn)C1×n

定理 11

  • x,y,zCnx,y,zCnλCλC,则

    1. 酉对称性

      (x,y)=¯(y,x)(x,y)=¯¯¯¯¯¯¯¯¯¯¯¯¯(y,x)

    2. 单线性性

      (x+y,z)=(x,z)+(y,z)(λx,y)=λ(x,y),(x,λy)=¯λ(x,y)(x+y,z)=(x,z)+(y,z)(λx,y)=λ(x,y),(x,λy)=¯¯¯λ(x,y)

    3. 正定性

      (x,x)0(x,x)0

      当且仅当 x=0x=0 时,有 (x,x)=0(x,x)=0.

定义 22

  • x=(x1,x2,,xn)TCnx=(x1,x2,,xn)TCn,令

    x=(x,x)=ni=1|xi|2=xHxx=(x,x)= ni=1|xi|2=xHx

    xx 为向量 xx长度22 - 范数.

定理 22

  • x,yCnx,yCnλCλC,则

    1. 非负性:当 x0x0 时,x>0x>0;当 x=0x=0 时,x=0x=0

    2. 齐次性

    λx=|λ|xλx=|λ|x

    1. 三角不等式

    x+yx+yx+yx+y

定理 33(Cauchy-Schwarz 不等式)

  • 对任意 x,yCnx,yCn,有

    |(x,y)|xy|(x,y)|xy

    (x,y)(y,x)(x,x)(y,y)(x,y)(y,x)(x,x)(y,y)

定义 33

  • x,yCnx,yCn,若 (x,y)=0(x,y)=0,则称 xxyy 正交,记作 xyxy.

定理 44

  • CnCn 中两两正交的非零向量组线性无关.

定义 44

  • x1,x2,,xnx1,x2,,xnCnCn 的一组基,若

    (xi,xj)=δij={0,ij,1,i=j.(xi,xj)=δij={0,ij,1,i=j.

    则称 x1,x2,,xnx1,x2,,xnCnCn 空间中的一组 标准正交基.

Schmidt 正交化

  • x1,x2,,xnx1,x2,,xnCnCn 空间这的一组基,令

    y1=x1y2=x2(x2,y1)(y1,y1)y1y3=x3(x3,y1)(y1,y1)y1(x3,y2)(y2,y2)y2yn=xn(xn,y1)(y1,y1)y1(xn,y2)(y2,y2)y2(xn,yn1)(yn1,yn1)yn1y1=x1y2=x2(x2,y1)(y1,y1)y1y3=x3(x3,y1)(y1,y1)y1(x3,y2)(y2,y2)y2yn=xn(xn,y1)(y1,y1)y1(xn,y2)(y2,y2)y2(xn,yn1)(yn1,yn1)yn1

    此时,y1,y2,,yny1,y2,,yn 两两正交(再进行单位化,即可得到一组标准正交基).

定义 55(酉矩阵)

  • ACn×nACn×n,若 AA 满足

    AHA=AAH=IAHA=AAH=I

    则称 AA酉矩阵.

📝 Notations

  • ACn×nACn×n,则

    AHA=AAH=IAHA=IAAH=IA1=AHAHA=AAH=IAHA=IAAH=IA1=AH

  • ARn×nARn×n 时,酉矩阵就是正交矩阵.

定理 55

  • ACn×nACn×n 是酉矩阵,则

    1. |detA|=1|detA|=1

    2. AT,AH,A1AT,AH,A1 仍为酉矩阵;

    3. BCn×nBCn×n 是酉矩阵,则 ABAB 也是酉矩阵.

定理 66

  • A=(x1,x2,,xn)Cn×nA=(x1,x2,,xn)Cn×n,则 AA 是酉矩阵的充要条件是 x1,x2,,xnx1,x2,,xnCnCn 中的一组标准正交基.

1.6 矩阵的酉相似

定义 11(酉相似)

  • A,BCn×nA,BCn×n,若存在酉矩阵 UU 使得

    U1AU=UHAU=BU1AU=UHAU=B

    则称 AABB 酉相似.

定理 11(Schur 分解定理)

  • ACn×nACn×n,则 AA 可酉相似于 上三角矩阵 TT,即存在 nn 阶酉矩阵 UU,使得

    U1AU=UHAU=T=(λ1λ2λn)U1AU=UHAU=T=⎜ ⎜ ⎜ ⎜ ⎜ ⎜λ1λ2λn⎟ ⎟ ⎟ ⎟ ⎟ ⎟

    其中 λ1,λ2,,λnλ1,λ2,,λnAA 的特征值.

定义 22(正规矩阵)

  • ACn×nACn×n,若 AA 满足

    AHA=AAHAHA=AAH

    则称 AA正规矩阵.

📝 Notations(正规矩阵的例子)

  • 实对称阵:ARn×n,AT=AARn×n,AT=A

  • 实反对称阵:ARn×n,AT=AARn×n,AT=A

  • 实正交矩阵:ARn×n,ATA=AAT=IARn×n,ATA=AAT=I

  • Hermite 矩阵:ACn×n,AH=AACn×n,AH=A

  • 反 Hermite 矩阵:ACn×n,AH=AACn×n,AH=A

  • 酉矩阵:ACn×n,AHA=AAH=IACn×n,AHA=AAH=I.

定理 22

  • ACn×nACn×n,则 AA 可以酉相似对角化的充要条件是

    AHA=AAHAHA=AAH

    AA 为正规矩阵.

定理 33

  • ACn×nACn×n 是正规矩阵,则

    1. AA 是 Hermite 矩阵 AA 的特征值全是实数;

    2. AA 是反 Hermite 矩阵 AA 的特征值是 00 或纯虚数;

    3. AA 是酉矩阵 AA 的特征值的模是 11

    4. λλAA 的特征值,xx 是对应 λλ 的特征向量,则 ¯λ¯¯¯λAHAH 的特征值,对应 ¯λ¯¯¯λ 的特征向量仍为 xx.

酉相似对角化方法:

  1. 求出 AA 的全部特征值,设 λ1,λ2,,λsλ1,λ2,,λsAA 的互不相同的特征值,其重数分别为 r1,r2,,rsr1,r2,,rs,且 r1+r2++rs=nr1+r2++rs=n.

  2. 对于特征值 λi,i=1,2,,sλi,i=1,2,,s,求出对应的 riri 个线性无关的特征向量 pi1,pi2,,piri,i=1,2,,spi1,pi2,,piri,i=1,2,,s.

  3. 用 Schmidt 正交化方法将 pi1,pi2,,piripi1,pi2,,piri 正交化再单位化得 ui1,ui2,,uiri,i=1,2,,sui1,ui2,,uiri,i=1,2,,s,则酉矩阵

    U=(u11,,u1r1,u21,,u2r2,,us1,,usrs,)U1AU=UHAU=Λ=(λ1Ir1λ2Ir2λsIrs)U=(u11,,u1r1,u21,,u2r2,,us1,,usrs,)U1AU=UHAU=Λ=⎜ ⎜ ⎜ ⎜ ⎜λ1Ir1λ2Ir2λsIrs⎟ ⎟ ⎟ ⎟ ⎟

定义 33

  • ACn×nACn×n 是一个 Hermite 矩阵,如果

    xHAx>0,xCn,x0xHAx>0,xCn,x0

    则称 AA 是一个 正定 的 Hermite 矩阵.

    如果

    xHAx0,xCnxHAx0,xCn

    则称 AA 是一个 非负定半正定)的 Hermite 矩阵.

定理 44(正定 Hermite 矩阵的性质)

  • ACn×nACn×n 是 Hermite 矩阵,则下述条件等价:

    1. AA 是正定的 Hermite 矩阵;

    2. AA 的特征值全为正数;

    3. 存在可逆矩阵 PP,使得 A=PHPA=PHP

    4. AA 的顺序主子式全为正数.

定理 55(非负定 Hermite 矩阵的性质)

  • ACn×nACn×n 是 Hermite 矩阵,则下述条件等价:

    1. AA 是非负定的 Hermite 矩阵;

    2. AA 的特征值全为非负数;

    3. 存在矩阵 PP,使得 A=PHPA=PHP.

定理 66

  • ACn×nACn×n,则

    1. AHAAHAAAHAAH 的特征值全为非负实数;

    2. AHAAHAAAHAAH 的非零特征值相同;

    3. rank(AHA)=rank(AAH)=rank(A)rank(AHA)=rank(AAH)=rank(A).

第二章 范数理论

2.1 向量的范数

定义 11(向量范数)

  • CnCn 上一个泛函,满足

    1. 正定性

      xCn,x0xCn,x0

      x=0x=0 当且仅当 x=0x=0.

    2. 齐次性

      λC,xCn,λx=|λ|xλC,xCn,λx=|λ|x

    3. 三角不等式

      x,yCn,x+yx+yx,yCn,x+yx+y

      则称 CnCn 上的一个 向量范数.

定理 11

  • 对任意 x,yCnx,yCn,有

    1. x=xx=x

    2. |xy|xyxyxy

定义 22

  • x=(x1,x2,,xn)TCnx=(x1,x2,,xn)TCn,定义

    x1=ni=1|xi|x2=(ni=1|xi|2)12xp=(ni=1|xi|p)1p,1p<+x=max1in|xi|x1=ni=1|xi|x2=(ni=1|xi|2)12xp=(ni=1|xi|p)1p,1p<+x=max1in|xi|

    可以验证 1,2,p,1,2,p, 均是 CnCn 上的向量范数,分别称为 11 - 范数22 - 范数pp - 范数 - 范数.

Hölder 不等式

  • 对任意 xi,yiC,i=1,2,,nxi,yiC,i=1,2,,n,有

    ni=1|xiyi|(ni=1|xi|p)1p(ni=1|yi|q)1qni=1|xiyi|(ni=1|xi|p)1p(ni=1|yi|q)1q

    其中 1p,q+1p,q+1p+1q=11p+1q=1.

Minkowski 不等式pp - 范数的三角不等式)

  • 对任意 xi,yiC,i=1,2,,nxi,yiC,i=1,2,,n,有

    (ni=1|xi+yi|p)1p(ni=1|xi|p)1p+(ni=1|yi|p)1p(ni=1|xi+yi|p)1p(ni=1|xi|p)1p+(ni=1|yi|p)1p

    其中 1p<+1p<+.

定理 22

  • xCn,limpxp=xxCn,limpxp=x.

📝 Notations

  • 如果 ACn×nACn×n 是 Hermite 正定矩阵,则 xA=xHAx,xCnxA=xHAx,xCn 也是 CnCn 上的向量范数.

定义 33

  • v1v1v2v2CnCn 上两个向量范数,如果存在常数 c1,c2>0c1,c2>0,使得

    xCn,c1xv1xv2c2xv1xCn,c1xv1xv2c2xv1

    则称向量范数 v1v1v2v2 等价.

向量范数的等价是 等价关系

  • 自反性:所有范数与自己等价.

  • 对称性:若向量范数 v1v1v2v2 等价,则向量范数 v2v2v1v1 等价.

  • 传递性:若向量范数 v1v1v2v2 等价,且向量范数 v2v2v3v3 等价,则向量范数 v1v1v3v3 等价.

定理 33

  • CnCn 上所有向量范数等价.

定义 44

  • x(k)=(x(k)1,x(k)2,,x(k)n)TCnx(k)=(x(k)1,x(k)2,,x(k)n)TCnk=1,2,k=1,2, 是一个向量序列,如果

    limkx(k)i=xi,i=1,2,,nlimkx(k)i=xi,i=1,2,,n

    则称向量序列 {x(k)}{x(k)} 收敛于 x=(x1,x2,,xn)Tx=(x1,x2,,xn)T,记为

    limkx(k)=xlimkx(k)=x

    否则,称 {x(k)}{x(k)} 发散.

定理 44

  • limkx(k)=xlimkx(k)x=0limkx(k)=xlimkx(k)x=0

    其中,CnCn 上任一向量范数.

2.2 矩阵范数

定义 11(矩阵范数)

  • Cn×nCn×n 上的一个泛函,满足

    1. 正定性

      ACn×n,A0ACn×n,A0

      其中 A=0A=0 当且仅当 A=OA=O.

    2. 齐次性

      λC,ACn×n,λA=|λ|AλC,ACn×n,λA=|λ|A

    3. 三角不等式

      A,BCn×n,A+BA+BA,BCn×n,A+BA+B

    4. 乘积不等式相容性):

      A,BCn×n,ABABA,BCn×n,ABAB

    则称 Cn×nCn×n 上的一个 矩阵范数.

定义 22

  • A=(aij)Cn×nA=(aij)Cn×n,定义

    • m1m1 - 范数

      Am1=ni=1nj=1|aij|Am1=ni=1nj=1|aij|

    • FF - 范数

      AF=(ni=1nj=1|aij|2)12=tr(AAH)AF=(ni=1nj=1|aij|2)12=tr(AAH)

    • mm - 范数

      Am=nmax1i,jn|aij|Am=nmax1i,jn|aij|

定义 33

  • mmCn×nCn×n 上的矩阵范数,vvCnCn 上的向量范数,如果对任意 ACn×nACn×nxCnxCn

    AxvAmxvAxvAmxv

    则称矩阵范数 mm 与向量范数 vv 相容.

定理 11

  • 矩阵范数 m1,Fm1,F 分别与 1,21,2 相容;

  • 矩阵范数 mm 与向量范数 1,2,1,2, 相容.

定义 44(由矩阵范数诱导的向量范数)

  • mmCn×nCn×n 上的一个矩阵范数,取 aCnaCn,且 a0a0. 定义

    xv=xaHm,xCnxv=xaHm,xCn

    可以证明,它是 CnCn 上的向量范数,称为由矩阵范数 mm诱导的 向量范数.

定理 22

  • Cn×nCn×n 是任意一矩阵范数 mm 与它所诱导的向量范数 vv 相容.

定义 55(由向量范数诱导的矩阵范数)

  • vvCnCn 上一个向量范数,定义

    Am=maxxv=1Axv=maxx0Axvxv,ACn×nAm=maxxv=1Axv=maxx0Axvxv,ACn×n

    可以证明(下面的定理 33),它是 Cn×nCn×n 上的矩阵范数,称为由向量范数 vv 诱导的 矩阵范数,简称 从属范数.

定理 33

  • vvCnCn 上任一向量范数,则它所诱导的 mm 为一个矩阵范数且与其相容.

定理 44

  • A=(aij)n×nCn×nA=(aij)n×nCn×n,记由向量 1122 范数导出的矩阵范数分别为 A1A1A2A2AA,则有

    A1=maxjni=1|aij|A2=λ1A=maxinj=1|aij|A1=maxjni=1|aij|A2=λ1A=maxinj=1|aij|

    其中,λ1λ1AHAAHA 的最大特征值.

    分别称 A1A1A2A2AA 为矩阵的 11 - 范数22 - 范数 - 范数.

定理 55

  • ACn×nACn×n,则

    1. AHF=AFAHF=AFAH2=A2AH2=A2

    2. 酉不变性:对任意酉矩阵 U,VCn×nU,VCn×n,有

      UAF=AVF=UAVFUA2=AV2=UAV2UAF=AVF=UAVFUA2=AV2=UAV2

    3. AA 是正规矩阵,且 λ1,λ2,,λnλ1,λ2,,λnAAnn 个特征值,则

      A2=maxk|λk|A2=maxk|λk|

定理 66

  • Cn×nCn×n 上所有矩阵范数等价.

长方阵的范数

对方阵范数的定义稍作修改可以推广到 m×nm×n 矩阵的情形:

  • 矩阵范数的 相容性:对任意 ACm×nACm×nBCn×lBCn×l

    ABABABAB

  • 矩阵范数 与向量范数的相容性:对任意 ACm×nACm×nxCnxCn

    AxAxAxAx

  • 从属范数

    A=maxxv=1Axv=maxx0Axvxv,ACm×nA=maxxv=1Axv=maxx0Axvxv,ACm×n

    其中 AxvCm 上的范数,xvCn 上的范数,它们应取为同类的向量范数.

常用矩阵范数

ACm×n

  • m1 范数

    Am1=mi=1nj=1|aij|

  • F 范数

    AF=mi=1nj=1|aij|2=tr(AHA)

  • M 范数(最大范数)

    AM=max{m,n}maxi,j|aij|

  • G 范数(几何平均范数)

    AG=nmmaxi,j|aij|

  • 1 范数(列和范数)

    A1=maxjmi=1|aij|

  • 2 范数(谱范数)

    A2=AHA的最大特征值

  • 范数(行和范数)

    A=maxinj=1|aij|

  • 矩阵范数的性质:

    1. F 范数和 2 范数是 酉不变 的.

    2. m1 范数与向量 1 范数相容;

      F 范数、G 范数与向量 2 范数相容;

      M 范数与向量 1,2, 范数相容.

    3. 矩阵的 1,2, 范数分别由向量 1,2, 范数导出,从而相容.

    4. Cm×n 上所有矩阵范数 等价.

2.3 范数应用举例

定义 1(矩阵的谱半径)

  • ACn×nλ1,λ2,,λnA 的特征值,则称

    ρ(A)=max1in|λi|

    为矩阵 A谱半径.

定理 1

  • ACn×n,则

    1. ρ(AH)=ρ(AT)=ρ(A)

    2. ρ(Ak)=[ρ(A)]k

    3. A 是正规矩阵时,ρ(A)=A2

定理 2

  • ρ(A)A,其中 AA 的任一矩阵范数.

定理 3

  • ACn×n,则 ε>0,必存在 Cn×n 上的矩阵范数 m,使得

    Amρ(A)+ε

引理

  • PCn×n,若对 Cn×n 上的某一矩阵范数 ,有 P<1,则 IP 可逆.

定理 4

  • ACn×nnδACn×n. 若对 Cn×n 上的某一矩阵范数 ,有 A1δA<1 ,则

    1. A+δA 可逆;

    2. (A+δA)1A11A1δA

    3. A1(A+δA)1A1A1δA1A1δA

定义 2(条件数)

  • ACn×nnCn×n 上的某一矩阵范数,称

    cond(A)=AA1

    为矩阵 A条件数(condition number).

推论

  • ACn×nnδACn×n,若对 Cn×n 上的某一矩阵范数 ,有 A1δA<1 ,则

    A1(A+δA)1A1cond(A)δAA1cond(A)δAA

定理 5

  • ACn×nnδACn×nb,δbCn,若对 Cn×n 上的某一矩阵范数 ,有 A1δA<1,则非齐次线性方程组

    Ax=b

    (A+δA)(x+δx)=b+δb

    的解满足

    δxvxvcond(A)1cond(A)δAA(δAA+δbvbv)

    其中 vCn 上与矩阵范数 相容的向量范数.

📝 Notations

  • 如果矩阵 A 的条件数很大,则称 A 对于求逆矩阵或求解线性方程组是 病态的(坏条件的),否则,称为 良态的(好条件的).

参考资料:

  1. 矩阵论 (武汉理工大学)

  2. 矩阵论简明教程(第三版). 徐仲,张凯院,陆全,冷国伟 编著

posted @   铃园轻梦  阅读(546)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· DeepSeek 开源周回顾「GitHub 热点速览」
· 物流快递公司核心技术能力-地址解析分单基础技术分享
· .NET 10首个预览版发布:重大改进与新特性概览!
· AI与.NET技术实操系列(二):开始使用ML.NET
· .NET10 - 预览版1新功能体验(一)
点击右上角即可分享
微信分享提示