抱歉,您的浏览器无法访问本站
本页面需要浏览器支持(启用)JavaScript
了解详情 >

高斯消元

image-20230103220433806

image-20230103220522818

旋转矩阵:

[cosθsinθsinθcosθ]\begin{bmatrix} \cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{bmatrix}

逆时针的旋转。

转置

tr(ATA)=0A=0tr(A^TA)=0 \Rightarrow A=0

行列式

余子式、代数余子式。

k=1naikAjk=A,i=j;0,ij\sum_{k=1}^n a_{ik}A_{jk}=|A|,i=j;0,i\not=j

λnλIm+AB=λmλIn+BA\lambda^n |\lambda I_m+AB|=\lambda^m|\lambda I_n+BA|

求:

kIαTβ|kI-\alpha^T\beta|

得到:

knkIαTβ=kkβαT=k(kβαT)k^n |kI-\alpha^T\beta|=k|k-\beta\alpha^T|=k(k-\beta\alpha^T)

范德蒙德行列式:

[111a1a2a3a12a22a32]\begin{bmatrix} 1 & 1 & 1\\ a_1 & a_2 & a_3 \\ a_1 ^2 &a_2 ^2 &a_3^2 \end{bmatrix}

det=i<j(ajai)det=\sum_{i<j} (a_j-a_i)

解:类似于拉格朗日插值,

ik(akb)ik(akai)\frac{\prod_{i \not=k} (a_k-b)}{\prod_{i\not=k}(a_k-a_i)}

可逆矩阵

BB 是元素全为 11nn 阶方阵,n2n \ge 2,求

  1. B2017B^{2017}

    不妨设α=(1,1,,1)T\alpha=(1,1,\cdots,1)^TB=ααTB=\alpha\alpha^TB2017=α(αTα)2016αT=n2016BB^{2017}=\alpha(\alpha^T\alpha)^{2016}\alpha^T=n^{2016}B

    这也告诉我们秩一矩阵的 kk 次幂是好求的。

  2. 试求 tt 的取值范围,使得矩阵 EtBE-tB 可逆,并求出这个矩阵的逆。

    一般地,设 f(x)f(x) 是一个多项式,kk 为常数,若 f(A)=Of(A)=O,而且

    f(x)=(xk)g(x)+c(c0)f(x)=(x-k)g(x)+c(c\not=0)

    则,

    (AkE)1=1cg(A)(A-kE)^{-1}=-\frac{1}{c}g(A)

    t=0t=0 时,显然可逆,t0t\not=0 时,

    EtB=t(B1tE)E-tB=-t(B-\frac{1}{t}E),当 t1nt\not=\frac{1}{n} 时,可逆得,到 (B1tE)1=t2nt1(B(n1t)E)(B-\frac{1}{t}E)^{-1}=\frac{t^2}{nt-1}(B-(n-\frac{1}{t})E),于是,结果应该是:

    tnt1(B(n1t)E)=tnt1B+E\frac{-t}{nt-1}(B-(n-\frac{1}{t})E)=-\frac{t}{nt-1}B+E

image-20230104201832480

在右边用列变换。

[OABO]1=[OB1A1O]\begin{bmatrix} O & A\\ B & O \end{bmatrix}^{-1}=\begin{bmatrix} O & B^{-1}\\ A^{-1} & O \end{bmatrix}

可逆包含的可交换的条件

A+B=ABA,BA+B=AB,证明 A,B 可交换

(AI)(BI)=I(BI)(AI)=I(A-I)(B-I)=I \Rightarrow (B-I)(A-I)=I

ABA1=BA1+3EABA^{-1}=BA^{-1}+3E

先转化为两个矩阵相乘为单位阵的情况:

(AE)BA1=3E(A-E)BA^{-1}=3E

那么:

BA1(AE)=B(EA1)=B(EA/A)=3EBA^{-1}(A-E)=B(E-A^{-1})=B(E-A^*/|A|)=3E

则:

B=3(EA/A)1B=3(E-A^*/|A|)^{-1}

伴随矩阵

A=(Aji)AA=AA=AEA^*=(A_{ji})\\A^*A=AA^*=|A|E


A,BA,B 均为 nn 阶矩阵,A=2,B=3|A|=2,|B|=-3,求 2AB1|2A^*B^{-1}|

注意 kA=kn1A,A=An1|kA|=k^{n-1}|A|,|A^*|=|A|^{n-1}


A=ATA^*=A^T

An1=A|A|^{n-1}=|A|

推出:若 A=0|A|=0,可能成为一组解,但若 AA 非零,注意到 r(A)=nr(A)=n 不满足 A=0|A|=0,而当 r(A)=n1r(A)=n-1 时,r(A)=1r(A^*)=1,只有在 n=2n=2 的情况成立,当 r(A)n2r(A)\le n-2 时,r(A)=0r(A^*)=0,更不成立。

nn 为奇数时,A=1|A|=1

nn 为偶数时,A=1|A|=-1A=1|A|=1

AE=AA=ATA|A|E=A^* A=A^TA

A=1|A|=1 时,AA 正交。

分块矩阵

AOdB=AB\begin{vmatrix} A & O\\ d & B \end{vmatrix} = |A||B|

证明:归纳证明。

AB=AB|AB|=|A||B|

构造:

D=A0EBD=\begin{vmatrix}A&0\\-E&B\end{vmatrix}

一方面 D=AB|D|=|A||B|,另一方面,

D=OABEB=(1)nABOBE=ABD=\begin{vmatrix}O&AB\\-E &B\end{vmatrix}=(-1)^n\begin{vmatrix}AB&O\\B &-E\end{vmatrix}=|AB|

于是成立。

[OABO]1=[OB1A1O]\begin{bmatrix}O&A\\B & O\end{bmatrix}^{-1}=\begin{bmatrix}O&B^{-1}\\A^{-1} & O\end{bmatrix}

矩阵的秩

矩阵的秩的定义

r+1r+1 阶子式全为 00rr 阶子式存在一个不为 00

矩阵的秩和方程有解的关系

r(A)=r(A,β)r(A)=r(A,\beta)

平时应该是 r(A)r(A,β)r(A)\le r(A,\beta)

想要证明 r(A)=r(B)r(A)=r(B),需要证明 Ax=0Bx=0Ax=0 \Rightarrow Bx=0,而且 Bx=0Ax=0Bx=0 \Rightarrow Ax=0

例如:

Ax=0ATAx=0ATAx=0xTATAx=0Ax=0Ax=0 \Rightarrow A^TAx=0\\ A^TAx=0 \Rightarrow x^TA^TAx=0 \Rightarrow Ax=0

A,BA,Bnn 阶方阵,且 r(A)=r(BA)r(A)=r(BA),证明 r(A2)=r(BA2)r(A^2)=r(BA^2)

A2x=0BA2x=0BA2x=0(BA)Ax=0A^2x=0 \Rightarrow BA^2 x=0\\ BA^2x=0 \Rightarrow (BA)Ax=0

BAx=0Ax=0BAx=0 \Rightarrow Ax=0,推出 Ax=0Ax=0,于是 A2x=0A^2x=0


Ax=0Ax=0Bx=0Bx=0 同解的充要条件是存在可逆矩阵 PP,满足 B=PAB=PA。(对应的 rref 相等)

矩阵的秩的大小关系

r(AB)min(r(A),r(B))r(AB)\le \min(r(A),r(B))


Am×n,Bn×sA_{m\times n},B_{n\times s},则:

r(A)+r(B)r(AB)+n(AB)r(A)+r(B) \le r(AB)+n(A的列数,B的行数)

AB=OAB=O 时,r(A)+r(B)nr(A)+r(B)\le n

image-20230106172719328


r(A+B)r(A,B)r(A,O;O,B)=r(A)+r(B)r(A+B)\le r(A,B)\le r(A,O;O,B)=r(A)+r(B)

r(A,O;C,B)r(A)+r(B)r(A,O;C,B)\ge r(A)+r(B)


r(A)=1,r(A)=n1r(A^*)=1, 当 r(A)=n-1 时。首先至少 r(A)=1r(A^*)=1。其次,r(AA)=r(AE)=0r(A^*A)=r(|A|E)=0,当 AA=OA^*A=O 时,r(A)+r(A)nr(A)+r(A^*)\le n,推出 r(A)1r(A^*)\le 1,所以 r(A)=1r(A^*)=1

r(A4B4)=2r(A^*_4B_4)=2,首先,推出 r(A4)2,r(B4)2r(A_4^*)\ge 2,r(B_4)\ge 2,其次,r(A)r(A^*) 只可能取 4,1,04,1,0,于是 AA 可逆。推出 r(B)=2r(B)=2


r(ATA)=r(A)r(A^TA)=r(A)

(通过解相同证明)

xTATAx=0(Ax,Ax)=0Ax=0x^TA^TAx=0 \Rightarrow (Ax,Ax)=0\Rightarrow Ax=0


初等变换不改变矩阵的秩。

r(A)=r(AB)=r(BA)Br(A)=r(AB)=r(BA) 当 B 可逆


An×mBm×n=En(n<m)A_{n\times m}B_{m\times n}=E_n (n<m),证明 r(Bm×n)=nr(B_{m\times n})=n

n=r(En)=r(An×mBm×n)r(Bm×n)=nn=r(E_n)=r(A_{n\times m}B_{m\times n}) \le r(B_{m\times n})=n(必须)

image-20230103231644260

ACm×nA \in \mathbb C^{m \times n}r(A)=m<nr(A)=m<nImI_mmm 阶单位矩阵,则下列说法正确的是:

  • 若矩阵 BB 满足 BA=0BA=0,则 B=0B=0,利用 r(A)+r(B)r(BA)+mr(A)+r(B)\le r(BA)+m,其中 r(A)=mr(A)=m,则 r(B)=0r(B)=0。而且必有解。

ACm×nA\in \mathbb C^{m\times n}AX=bAX=b 为一非齐次线性方程组(b0b\not=0),则:

  • m<nm<n,则方程组 AX=bAX=b 有非零解。r(A)r(A) 不一定等于 r(A,b)r(A,b),例如,取 A=0A=0
  • r(A)=mr(A)=m,则方程组 AX=0AX=0 有非零解。不一定,比如取 A=[1,0],b=[0;1]A=[1,0],b=[0;1],则无解。
  • r(A)=nr(A)=n,则 AX=0AX=0 只有零解。r(A)+r(X)r(AX)+nr(A)+r(X)\le r(AX)+n,得到 r(X)=0r(X)=0

nn 阶矩阵 AAA(A+kE)=0r(A)+r(A+kE)=nA(A+kE)=0 \Leftrightarrow r(A)+r(A+kE)=n

证明:

r(A)+r(A+kE)r(0)+nr(A)+r(A+kE)\le r(0)+n

r(A)+r(A+kE)r(A+kEA)=nr(A)+r(A+kE)\ge r(A+kE-A)=n

线性空间与线性变换

线性空间

  • α+β=β+α\alpha+\beta=\beta+\alpha
  • α+(β+γ)=(α+β)+γ\alpha+(\beta+\gamma)=(\alpha+\beta)+\gamma
  • 零元:0V,0+α=α\exists 0 \in \mathbb V,0+\alpha=\alpha
  • 负元:α,β,s.t.α+β=0\forall \alpha,\exists \beta ,s.t. \alpha+\beta=0
  • 单位元:1F,1α=α1\in \mathbb F,1\alpha=\alpha
  • (kl)α=k(lα)(kl)\alpha=k(l\alpha)
  • (k+l)α=lα+kα(k+l)\alpha=l\alpha+k\alpha
  • k(α+β)=kα+kβk(\alpha+\beta)=k\alpha+k\beta

子空间

子空间需要对加法、数乘封闭即可。

W,UW,UVV 的子空间,则 WUW\cap U 也是 VV 的子空间。 如果互不包含,则 WVW\cup V 不是 VV 的子空间。

和子空间:

W+U={w+uwW,uU}W+U=\{w+u|w\in W,u\in U\}

生成子空间

L(α1,α2,α3,,αm)L(\alpha_1,\alpha_2,\alpha_3,\cdots,\alpha_m)

线性相关,线性组合。

α1,α2,,αn\alpha_1,\alpha_2,\cdots,\alpha_n

线性无关,而且 VV 中任意向量都可以用 α1,α2,,αn\alpha_1,\alpha_2,\cdots,\alpha_n 线性表示,则称 VVnn 维线性空间。

dim(V)\dim(V) 基,坐标。

A={AxxV}kerA={xVAx=0}\Im \mathcal A=\{\mathcal Ax|x \in V\} \\ \ker \mathcal A=\{x\in V| \mathcal Ax =0\}

核和像相等的条件:A2=0\mathcal A^2=0

A(x,y)=A(y,0)A=kerA={(x,0)xR}\mathcal A(x,y)=\mathcal A(y,0)\\ \Im A=\ker A=\{(x,0)|x\in \R \}

dim+dimker=n\dim \Im +\dim \ker=n

V={Xn×mAX=0}V=\{X_{n\times m} | AX=0\} 的维数:

XX 的每一列都有 nr(A)n-r(A) 种选择,因此是 (nr(A))m(n-r(A))^m

dimV1+dimV2=dimV1V2+dimV1+V2\dim V_1+\dim V_2=\dim V_1 \cap V_2 + \dim V_1+V_2

内积,实数域上的线性空间,欧氏空间。

坐标

A=(α1,α2,α3)xAxA=(\alpha_1,\alpha_2,\alpha_3)\\x 是列向量,Ax

过渡矩阵

[η]=[ε]C[\eta]=[\varepsilon]C

CC 是可逆矩阵。

[η][\eta] 下坐标是 YY[ε][\varepsilon] 下坐标是 XX,则:

CY=XCY=X

内积的性质

  1. 对称性
  2. 线性性
  3. 正定性 (α,α)0(\alpha,\alpha)\ge0(α,α)=0(\alpha,\alpha)=0 充要条件是 α=0\alpha=0

如果定义 (α,β)=(Aα)T(Aβ)(\alpha,\beta)=(A\alpha)^T(A\beta),其中 AA 非奇异,可以满足正定性、线性性的要求,而对称性:

(α,β)=αTATAβ=(αTATAβ)T=βTATAα=(β,α)(\alpha,\beta)=\alpha^T A^T A\beta=(\alpha^T A^T A\beta)^T= \beta^T A^TA\alpha=(\beta,\alpha)

柯西施瓦茨不等式

(α,β)αβ|(\alpha,\beta)| \le |\alpha| \cdot |\beta|

等号成立当且仅当 α,β\alpha,\beta 线性相关。

可以进行角度的定义:

cosθ=(α,β)αβ\cos \theta =\frac{(\alpha,\beta)}{|\alpha||\beta|}

可以进行长度的定义:

α=(α,α)|\alpha|=\sqrt{(\alpha,\alpha)}

施密特正交化

αs\alpha_sβt\beta_t 上的投影:

(αs,βt)(βt,βt)βt\frac{(\alpha_s,\beta_t)}{(\beta_t,\beta_t)}\beta_t

原因:因为 (αs,βt)(\alpha_s,\beta_t) 代表 αs\alpha_sβt\beta_t 上的投影的长度乘以 βt\beta_t 的长度。

正交矩阵的定义

AAT=ATA=EAA^T=A^TA=E

注意两个都要写到。


AA 为实反对称矩阵,α\alphann 维实列向量,且 Aα=βA\alpha=\beta,证明 α\alphaβ\beta 正交。

αTAα+(αTAα)T=0αTAα=αTβ=0\alpha^TA\alpha+(\alpha^TA\alpha)^T=0 \Rightarrow \alpha^TA\alpha=\alpha^T\beta=0

线性变换

σ(A+B)=σ(A)+σ(B)σ(kA)=kσ(A)\sigma(A+B)=\sigma(A)+\sigma(B)\\\sigma(kA)=k\sigma(A)

记为:

β=σ(α),β=σα\beta=\sigma(\alpha),\beta=\sigma\alpha

线性变换在给定基下的矩阵

A(ε1,ε2,,εn)=(Aε1,Aε2,,Aεn)=(ε1,ε2,,εn)A\mathcal A(\varepsilon_1,\varepsilon_2,\cdots,\varepsilon_n)=(\mathcal A \varepsilon_1,\mathcal A \varepsilon_2,\cdots,\mathcal A \varepsilon_n)=(\varepsilon_1,\varepsilon_2,\cdots,\varepsilon_n)A

意义是,只要找到一组 (ε1,ε2,εn)(\varepsilon_1,\varepsilon_2\cdots,\varepsilon_n) 与和其对应变换过后的基,就可以知道 AA 的全貌。

基变换

[α1,α2,α3]γ=[β1,β2,β2]γ[\alpha_1,\alpha_2,\alpha_3] \gamma=[\beta_1,\beta_2,\beta_2] \gamma'

希望 γ=Aγ\gamma=A\gamma',因此:

A=[α1,α2,α3]1[β1,β2,β2]A=[\alpha_1,\alpha_2,\alpha_3]^{-1}[\beta_1,\beta_2,\beta_2]

也可以这样快速记忆,要从 α\alpha 转到 β\beta,需要先消除 α\alpha 的信息,再转到 β\beta 的信息。

β\beta 转换到 α\alpha,矩阵就是 A1A^{-1}

线性变换在不同基矩阵下的关系

[ε],[η][\varepsilon],[\eta] 是基,已知 A([ε])\mathcal A([\varepsilon]) 对应矩阵 AA,和 [ε]C=[η][\varepsilon]C=[\eta]A([η])\mathcal A([\eta]) 对应矩阵 BB

A([η])=A([ε]C)=A([ε])C=[ε]AC=[η]C1AC\mathcal A([\eta])=\mathcal A([\varepsilon]C)=\mathcal A([\varepsilon])C=[\varepsilon]AC=[\eta]C^{-1}AC

也可以考虑 BBη\etaη\eta 的矩阵,知道 εη\varepsilon \to \eta 的变换,就需要 (εη)1×εε×(εη)(\varepsilon \to \eta)^{-1} \times \varepsilon \to \varepsilon \times (\varepsilon \to \eta)


在线性空间 P[x]n+1P[x]_{n+1} 中,求向量:f

f(x)=a0+a1x+a2x2++anxnf(x)=a_0+a_1x+a_2x^2+\cdots+a_nx^n

在基:

1,(xa),(xa)2,,(xa)n1,(x-a),(x-a)^2,\cdots,(x-a)^n

下的坐标。

f(x)f(x)x=ax=a 点处展开,可以得到 f(x)=f(a)+f(a)1(xa)++f(n)(a)n!(xa)nf(x)=f(a)+\frac{f'(a)}{1}(x-a)+\cdots+\frac{f^{(n)}(a)}{n!}(x-a)^n

利用泰勒公式展开的唯一性。

x2,(x+1)2,(x+2)2,(x+3)2x^2,(x+1)^2,(x+2)^2,(x+3)^2

线性相关。

HouseHolder矩阵

A=E+kαTαA=E+k\alpha^T\alpha

第一种特征向量是 αx=0\alpha x=0 的基础解系中的向量。

例如 α=(a1,a2,a3)\alpha=(a_1,a_2,a_3),得出 α1=(a2,a1,0)T,α2=(a3,0,a1)T\alpha_1=(-a_2,a_1,0)^T,\alpha_2=(-a_3,0,a_1)^T,一共 n1n-1 个。对应的特征值是 11

第二种特征向量就是 αT\alpha^T,对应特征值是 1+k(α,α)1+k(\alpha,\alpha)。如果是正交矩阵,需要特征值 =1=-1,那么 k=2(α,α)k=\frac{-2}{(\alpha,\alpha)},如果需要正定,则 k>2(α,α)k>\frac{-2}{(\alpha,\alpha)}

矩阵的特征值

(λ,α)A,(λ1,α)A1(λ,α)A,(g(λ),α)g(A)A,AT(\lambda,\alpha) A,(\lambda^{-1},\alpha)A^{-1}\\ (\lambda,\alpha) A,(g(\lambda),\alpha)g(A)\\ A,A^T 相同

nn 阶方阵的属于不同特征值的特征向量线性无关。

λi=tr(A)λi=A\sum \lambda_i=tr(A)\\ \prod \lambda_i=|A|

AA^* 的特征值为 λ1A\lambda^{-1}|A|。若 AA 可逆,则 AA 对应于 λ\lambda 的特征向量也为 AA^* 对应于 A/λ|A|/\lambda 的特征向量。

B=P1APB=P^{-1}AP,对应的特征向量为 P1αP^{-1}\alpha

AA 的每行元素都是 aa,取特征向量为 (1,1,,1)T(1,1,\cdots,1)^T ,则得出一个特征值为 aa


证明:

i=1nλi2=i=1nk=1naikaki\sum_{i=1}^n \lambda_i^2 = \sum_{i=1}^n\sum_{k=1}^n a_{ik}a_{ki}

注意到 i=1nλi2=tr(A2)\sum_{i=1}^n \lambda_i^2=tr(A^2),而 tr(A2)tr(A^2) 展开就是右边的形式。


秩1矩阵的特征值为 n1n-1 个 0 和 一个 tr(A)tr(A)。对应的向量为 αTx=0\alpha^Tx=0 的解空间和 αT\alpha^T

代入 A=α,B=αTA=\alpha,B=\alpha^T,就可以得到秩一矩阵的性质。

判断矩阵是否相似于一个对角阵,主要是看是否有 nn 个线性无关的特征向量。


经典题型:

  1. 已知 AA 的所有特征值,求 f(A)|f(A)| 的行列式:AA^*A1A^{-1} 代替,

f(λi)\prod f(\lambda_i)

  1. 实对称矩阵 AA,已知所有特征值和部分特征对,求 AA。不同特征值对应的特征向量正交。
  2. AA^* 的特征值,长的是 λ1λi1λi+1,λn\lambda_1\cdots\lambda_{i-1}\lambda_{i+1},\cdots\lambda_n

AA 的特征对 (λ,α)(\lambda,\alpha),则 f(A)f(A) 的特征对 (f(λ),α)(f(\lambda),\alpha)

AA 实反对称矩阵,证明:

  1. A2A^2 半负定矩阵。

    xTA2x=xTATAx=(Ax)T(Ax)0x^TA^2x=-x^TA^TAx=-(Ax)^T(Ax)\le0

  2. EA2>1|E-A^2|>1,相似的,半负定矩阵的所有特征值 0\le0,令 f(x)=1xf(x)=1-x,则 EA2E-A^2 的所有特征值为 f(λi)1f(\lambda_i)\ge1,而且特征值不全为 00,推出 EA2>1|E-A^2|>1,这道题主要用到了特征值和行列式的关系。


限定 kk 重特征值:r(λEA)=nkr(\lambda E-A)=n-k。几何重数+原矩阵的秩=nn。例如 (AE)B=0,r(B)=r(A-E)B=0,r(B)=r,得出 r(AE)=nrr(A-E)=n-r,则 11 对应的几何重数应该是 rr


A2=A,Adiag(1,1,1,,1,,0)A^2=A,A \sim diag(1,1,1,\cdots,1,\cdots,0)

A2=0,Adiag(0,,0)A^2=0,A\sim diag(0,\cdots,0)

可对角化等价于最小多项式没有重根。

矩阵的相似对角化

正交相似

AQ=QΛ,QTQ=EAQ=Q\Lambda,Q^TQ=E

QTAQ=ΛQ^TAQ=\Lambda

二次型

合同

非奇异 PPs.t.s.t. PTAP=BP^TAP=B

惯性系数相同。

  • 自反性
  • 对称性
  • 传递性
  • 合同矩阵必有相同的秩
  • AA 对称推出 BB 对称
  • f(x)=xTAxf(x)=x^TAx 经过可逆变换 x=Pyx=Py 化为二次型 g(x)=yTByg(x)=y^TBy 充要条件是 PTAP=BP^TAP=B

正定的性质

如果 AA 正定,则:

  • 正定的定义,xTAx>0x^TAx>0x0x\not=0

  • A=QTΛQA=Q^T\Lambda Q,其中 QQ 是正交阵,如果令 Λ\Lambda 分解为 Λ1/2\Lambda^{1/2},则可以得到 A=PTPA=P^TP,但此时 PP 不满足正交的性质。

  • A=PTPA=P^TP,其中 PP 非奇异。从这个性质可以根据定义推出 AA 正定。

  • AA对称阵A=ATA=A^T,由上一条性质可知。

  • xx 取正交替换 x=Qyx=Qy,则还是满足正定性的要求,因为 xTAx=yTQTAQyx^TAx=y^TQ^TAQy,则 B=QTAQB=Q^TAQ 依然是正定矩阵。

如果 AA 正定,则下属命题等价。

  • AA 为正定矩阵。
  • AA 的特征值全部大于 0。
  • AA 合同于单位矩阵 EE
  • 存在非奇异矩阵 MM,使得 A=MTMA=M^TM

判断正定性的方法:主子式法,定义法。

已知 A,BA,B 均为 nn 阶正定矩阵,则 kA+lBkA+lB 为正定矩阵,其中 k,l>0k,l >0

正定矩阵的对角线元素全部 >0>0,如果 0\le 0,取 (0,,1,,0)(0,\cdots,1,\cdots,0),就会出现反例。

AA 正定推出 A1A^{-1} 正定,因为 A=QTQA=Q^TQA1=(QT)T(QT)A^{-1}=(Q^{-T})^T(Q^{-T}),也是正定。反之也能成立。

AA 正定也能推出 AkA^k 正定,因为 A,AA,A 可交换。

A,BA,B 正定,可以推出 BABBAB 正定。

nn 阶方阵 AA 是反对称矩阵当且仅当对于任意的 nn 维向量 α\alpha,有 αTAα=0\alpha^T A\alpha=0

定义:标准型,规范标准型

规范标准型必须全是 11

一些题目

如果 A,BA,B 正定,证明 ABAB 正定与 AB=BAAB=BA 等价。

正向:

AB=(AB)T=BTAT=BAAB=(AB)^T=B^TA^T=BA

反向:

由于 BB 正定,所以 BB 对称,可以分解为 PΛPTP\Lambda P^T,其中 PP 是正交阵。

于是 B1/2B^{1/2} 可以定义为 PΛ1/2PTP\Lambda^{1/2} P^T,易得 B1/2B^{1/2} 是正交矩阵,是对称的。ABABB1/2AB1/2B^{1/2}AB^{1/2} 相似,具有相同的特征值,只需证明 B1/2AB1/2B^{1/2}AB^{1/2 } 正定。

xTB1/2AB1/2x=(B1/2x)TA(B1/2x)x^TB^{1/2}AB^{1/2}x=(B^{1/2}x)^TA(B^{1/2}x),根据 AA 正定可得 B1/2AB1/2B^{1/2}AB^{1/2} 正定。

于是命题得证。


另外,如果 A,BA,B 正定,则 ABAB 的特征值都大于 00

QABQ1=(PQT)TPQTQABQ^{-1}=(PQ^T)^TPQ^T

可以知道 ABAB 相似于一个正定矩阵,所以特征值全部大于 00。如果要 ABAB 正定,还需要 ABAB 对称,也就是 (AB)T=BA=AB(AB)^T=BA=AB,也就是可交换。


如果 AA 是正交阵而且正定,则证明 A=EA=E

是正定,推出 A=PTPA=P^TPAT=AA^T=A,是正交,推出 AAT=EAA^T=E,则 A2=EA^2=E。可以看成 g(λ)=λ21g(\lambda)=\lambda^2-1g(λ)g(\lambda)A2EA^2-E 的特征值,必为 00,则有两个解 λ=1,1\lambda=1,-1。然而,AA 正定,则 λ=1\lambda=-1 需要舍去,推出 AA 特征值全为 11,则 A=B1EB=EA=B^{-1}EB=E


Bm×nB_{m\times n} 实矩阵,证明 BTBB^TB 正定 \Leftrightarrow r(B)=nr(B)=n

BTBB^TB 正定则 xTBTBx=Bx2>0x^TB^TBx=|Bx|^2 >0(对于 x0,Bx0\forall x\not=0,Bx\not=0),则 r(B)=nr(B)=n

反过来也好证。


A,BA,Bnn 阶正定矩阵,方程 f(x)=xAB=0f(x)=|xA-B|=0 的根均为 11,证明:

A=BA=B

A=MTMA=M^TM

xAB=xMTMB=MTxEMTBM1M=AxEMTBM1|xA-B|=|xM^TM-B|=|M^T||xE-M^{-T}BM^{-1}||M|=|A||xE-M^{-T}BM^{-1}|

xEMTBM1|xE-M^{-T}BM^{-1}| 的根均为 11

MTBM1M^{-T}BM^{-1} 是正定的,因为定义。对其进行分解,为 PTΛPP^{T}\Lambda P,其中 PP 正交。

xEMTBM1=xEΛ|xE-M^{-T}BM^{-1}|=|xE-\Lambda|

推出:Λ=E\Lambda=E。则 PTMTBM1P=Λ=EP^{T}M^{-T}BM^{-1}P=\Lambda=E,推出 B=MTM=AB=M^TM=A


假设 α,β\alpha,\beta 为正交的三维单位列向量,二次型 f(x)=xT(2ααT+ββT)xf(x)=x^T(2\alpha\alpha^T+\beta\beta^T)x 经过正交替换可以化为标准型?

正交:经常需要乘一下。

r(A)r(2ααT)+r(ββT)=2r(A)\le r(2\alpha\alpha^T)+r(\beta\beta^T)=2

f(x)=xTAxf(x)=x^TAx,得到 Aα=2αA\alpha=2\alphaAβ=βA\beta=\beta。而且由秩的性质,得到 λ3=0\lambda_3=0


补充矩阵的秩一分解:

A=(α1,,αn)diag(λ1,,λn)(α1,,αn)T=λ1α1α1T+A=(\alpha_1,\cdots,\alpha_n) diag(\lambda_1,\cdots,\lambda_n) (\alpha_1,\cdots,\alpha_n)^T\\=\lambda_1\alpha_1\alpha_1^T+\cdots

其中 α1,,αn\alpha_1,\cdots,\alpha_n 是标准正交基。

于是可以很容易看出 AA 的三个特征值为 2,1,02,1,0

计算 AnA^n,我们也可以直接将 λ1,\lambda_1,\cdots 进行 mm 次方。


二次型 f(x)=xTααTxf(x)=x^T\alpha\alpha^Tx 化为标准型:

ααT\alpha\alpha^T 的特征值:0(n1)αTα(0)0(n-1个),\alpha^T\alpha(\ge0)


正定矩阵可以定义“开根号”

AA 正定,证明存在正定矩阵 BB,使得 A=B2A=B^2

A=P1ΛPA=P^{-1}\Lambda P

其中 PP 是正交阵,Λ\Lambda 每一个对角线元素都是正的,则令 B=P1Λ1/2PB=P^{-1}\Lambda^{1/2} P


AAmnm*n 阶实矩阵,证明 r(A)=nr(A)=n 的充要条件是 ATAA^TA 正定。

充分条件:$r(A)=n \Rightarrow $ 若存在 xx 使得 Ax=0Ax=0,则 r(A)+r(x)r(Ax)+nr(A)+r(x) \le r(Ax)+n,推出 r(x)=0,x=0r(x)=0,x=0,因此,xTATAx=(Ax,Ax)0x^TA^TAx=(Ax,Ax) \ge 0,等号当且仅当 x=0x=0 取到。

必要条件:xTATAx0x^TA^TAx \ge 0 等号当且仅当 x=0x=0 取到,推出 r(A)=nr=nr(A)=n-r=n


An,BnA_n,B_n 实对称,AA 正定,证明存在可逆矩阵 PP,使得 PTAP,PTBPP^TAP,P^TBP 均为对角阵。


AAmm 阶正定矩阵,BBm×nm\times n 实矩阵,证明 BTABB^TAB 为正定矩阵的充分必要条件是矩阵 BB 的秩数 r(B)=nr(B)=n

xTBTABx=(Bx)TA(Bx)x^TB^TABx=(Bx)^TA(Bx)

需要 x=0Bx=0x=0 \Leftrightarrow Bx=0。因此 r(B)=nr(B)=n。因此,不能随便说 BTBB^TB 正定,还要求 r(B)=nr(B)=n


证明:可逆实矩阵 AA 可以表示为一个正定矩阵与一个正交矩阵的积,且表示唯一。

ATA=B12,B1A^TA=B_1^2,B_1 正定

Q=AB11,QTQ=(B1)TATAB11=B1TB1=(B1B11)T=EQ=AB_1^{-1},Q^TQ=(B_1)^{-T} A^T A B_1^{-1}=B_1^{-T}B_1 = (B_1 B_1^{-1})^T=E

于是,A=QB1A=QB_1

评论