syntax.garden

27 正算子

27 正算子

正算子

这讲讲什么:在自伴算子中挑出「特征值全非负」的那一类,证明它恰好就是能开平方的那一类,并证明这个平方根唯一
前置:22 讲(内积与范数)、25 讲(自伴与正规算子)、26 讲(谱定理)
一句话:正算子是算子世界的非负数——它由 ⟨Tv,v⟩≥0\langle Tv,v\rangle\ge 0 定义,却等价于「有唯一的非负平方根」。


一、为什么需要它

26 讲的谱定理告诉我们:自伴算子一定能用规范正交基对角化,而且对角线上全是实数。

这句话已经把问题解决了一半。剩下的那一半是:实数还分正负。于是立刻冒出一个自然的子类:

哪些自伴算子的特征值全都落在非负半轴上?

这个子类就是正算子。它不是数学家为了分类而分类——有四条独立的理由说明这个类必须单独拎出来。

理由一:它正好是能开方的那一类。 在实数里,「非负」与「有平方根」是同一件事:−1-1 没有实平方根,22 有。在算子里,一般算子未必有平方根,有了也未必唯一。而正算子不但一定有平方根,还只有一个正的。开方是这套理论的基本操作(极分解、SVD、矩阵函数都要用),而它成立的边界恰好就是正算子。

理由二:T∗TT^*T 永远落在这一类里。 对任意线性映射 T:V→WT:V\to W(不要求方阵、不要求自伴),合成 T∗TT^*T 一定是 VV 上的正算子。这不是巧合,它是「长度平方」的算子版本:

∥Tv∥2=⟨Tv,Tv⟩=⟨T∗Tv,v⟩\|Tv\|^2=\langle Tv,Tv\rangle=\langle T^*Tv,v\rangle

你熟的那条长度公式 ∥v∥2=⟨v,v⟩\|v\|^2=\langle v,v\rangle 只是 T=IT=I 的特例。29 讲的奇异值分解,整套东西都建立在「T∗TT^*T 正,所以能取 T∗T\sqrt{T^*T}」这一句话上。

理由三:它是二次型的符号判据。 化学与工程里大量出现二次型 xTAxx^{\mathsf T}Ax(能量、方差、应变能)。判断「这个二次型会不会取负值」,等价于判断 AA 是不是正算子。这不是抽象问题——优化里判断极小值、数值上判断矩阵能否 Cholesky 分解、统计里判断协方差矩阵是否合法,都是同一件事。

理由四:记号统一。 没有这个概念,你会写出一堆平行的小结论:「特征值非负的算子有这个性质」「能写成 R∗RR^*R 的算子有那个性质」「投影算子有第三个性质」——而它们其实是同一类对象的三种描述。有了「正算子」,这些结论合成一条定理(下一节的六条等价)。

一句话概括动机

谱定理把算子变成对角矩阵;正算子是其中「对角线全是非负数」的那一半。
与非负数一样,它们可以比较大小、可以开方、可以由「平方」造出来。


二、定义与最小例子

定义

设 VV 是内积空间,T∈L(V)T\in\mathcal{L}(V)。称 TT 是正算子,如果

  1. TT 是自伴的,即 T∗=TT^*=T;
  2. 对每一个 v∈Vv\in V,都有 ⟨Tv,v⟩≥0\langle Tv,v\rangle\ge 0。

写成一条:

  T 正  ⟺  T=T∗ 且 ∀v∈V: ⟨Tv,v⟩≥0  \boxed{\;T\ \text{正}\iff T=T^*\ \text{且}\ \forall v\in V:\ \langle Tv,v\rangle\ge 0\;}

为什么「自伴」这一条不能省(在实空间里)

只有当 ⟨Tv,v⟩\langle Tv,v\rangle 是实数时,「≥0\ge 0」才有意义。
自伴保证它取实值;反过来,若已知 ⟨Tv,v⟩\langle Tv,v\rangle 对一切 vv 都是实数,
那么在复空间上可以推出 TT 自伴。所以复情形下第一条可以划掉,实情形下不行——
后面第五节会给出一个「⟨Tv,v⟩≡0\langle Tv,v\rangle\equiv 0 却不自伴」的反例。

最小例子:00 与 II

  • 零算子:Tv=0Tv=0,于是 ⟨Tv,v⟩=0≥0\langle Tv,v\rangle=0\ge 0,自伴。正。
  • 恒等算子:⟨Iv,v⟩=⟨v,v⟩=∥v∥2≥0\langle Iv,v\rangle=\langle v,v\rangle=\|v\|^2\ge 0,自伴。正。
  • 数乘:cIcI 正   ⟺  c≥0\iff c\ge 0。(−1)I(-1)I 不正(取 v≠0v\ne 0 得 −∥v∥2<0-\|v\|^2<0)。

这三个例子说明:正算子是「非负标量」在算子层的推广,而且推广得很忠实。

例 1:一个必须动手算的 2×22\times2 矩阵

设 T∈L(F2)T\in\mathcal{L}(\mathbb{F}^2) 关于标准基的矩阵是

A=(2−1−11)A=\begin{pmatrix}2&-1\\ -1&1\end{pmatrix}

它对称,所以自伴。对 v=(w,z)v=(w,z),

T(w,z)=(2w−z, −w+z)T(w,z)=(2w-z,\ -w+z)

于是

⟨T(w,z),(w,z)⟩=(2w−z)w‾+(−w+z)z‾=2∣w∣2−zw‾−wz‾+∣z∣2=2∣w∣2−2Re⁡(wz‾)+∣z∣2\begin{aligned}\langle T(w,z),(w,z)\rangle&=(2w-z)\overline{w}+(-w+z)\overline{z}\\ &=2|w|^2-z\overline{w}-w\overline{z}+|z|^2\\ &=2|w|^2-2\operatorname{Re}(w\overline{z})+|z|^2\end{aligned}

这个式子看着不好判断符号。配方:

2∣w∣2−2Re⁡(wz‾)+∣z∣2=(∣w∣2−2Re⁡(wz‾)+∣z∣2)+∣w∣2=∣w−z∣2+∣w∣22|w|^2-2\operatorname{Re}(w\overline{z})+|z|^2=\big(|w|^2-2\operatorname{Re}(w\overline{z})+|z|^2\big)+|w|^2=|w-z|^2+|w|^2

两个模平方之和,必然 ≥0\ge 0。所以 TT 是正的。

这个例子教了什么

不要试图对每个 vv 硬验。 真正管用的是:算出 ⟨Tv,v⟩\langle Tv,v\rangle 的表达式,
然后把它写成平方和。这也是「正定二次型配方法」的算子版本。
顺便:det⁡A=2⋅1−(−1)(−1)=1>0\det A=2\cdot1-(-1)(-1)=1>0,tr⁡A=3>0\operatorname{tr}A=3>0,特征值 3±52\frac{3\pm\sqrt5}{2} 确实都为正。

例 2:正交投影是正的

设 UU 是 VV 的子空间,PUP_U 是正交投影。它自伴(24 讲已证)。再看内积:

⟨PUv,v⟩=⟨PUv, PUv+(v−PUv)⟩=⟨PUv,PUv⟩+⟨PUv,v−PUv⟩\langle P_Uv,v\rangle=\langle P_Uv,\ P_Uv+(v-P_Uv)\rangle=\langle P_Uv,P_Uv\rangle+\langle P_Uv,v-P_Uv\rangle

第二项为零,因为 PUv∈UP_Uv\in U 而 v−PUv∈U⊥v-P_Uv\in U^\perp。所以

⟨PUv,v⟩=∥PUv∥2≥0\langle P_Uv,v\rangle=\|P_Uv\|^2\ge 0

投影算子全是正的,而且这个证明一行都不涉及坐标——只用了「正交分解」和「自伴」。

反过来成立吗

不正。正算子不必是投影:diag⁡(1,2)\operatorname{diag}(1,2) 是正的,但它把 (1,1)(1,1) 拉长成 (1,2)(1,2),不是投影。
准确的说法是:投影是特征值只取 00 和 11 的正算子。正算子允许特征值取任何非负数。

例 3:T2+bT+cIT^2+bT+cI 型

设 TT 自伴,b,c∈Rb,c\in\mathbb{R} 且 b2<4cb^2<4c(即多项式 x2+bx+cx^2+bx+c 没有实根、恒正)。那么

S=T2+bT+cIS=T^2+bT+cI

是正算子。理由:配方

S=(T+b2I)2+(c−b24)IS=\Big(T+\tfrac b2 I\Big)^2+\Big(c-\tfrac{b^2}{4}\Big)I

第一项是自伴算子的平方,形如 R∗RR^*R(R=T+b2IR=T+\tfrac b2I 自伴),因而是正的;第二项是正数乘 II,也是正的;正算子相加仍正(第六节习题)。所以 SS 正。

一个非例

B=(100−1)B=\begin{pmatrix}1&0\\ 0&-1\end{pmatrix}

自伴,但取 v=(0,1)v=(0,1) 得 ⟨Bv,v⟩=−1<0\langle Bv,v\rangle=-1<0。不正。它的特征值一正一负——这正是下一节等价刻画要抓的东西。


三、几个必须自己算一遍的结论

3.1 六条等价刻画(本讲的主定理)

对 T∈L(V)T\in\mathcal{L}(V),下列六条等价:

编号 条件
(a) TT 是正算子
(b) TT 自伴,且所有特征值非负
(c) 关于 VV 的某个规范正交基,TT 的矩阵是对角矩阵且对角元全为非负数
(d) TT 有正平方根(即存在正算子 RR 使 R2=TR^2=T)
(e) TT 有自伴平方根
(f) 存在 R∈L(V)R\in\mathcal{L}(V) 使 T=R∗RT=R^*R

证明方式是绕一圈:

(a)⇒(b)⇒(c)⇒(d)⇒(e)⇒(f)⇒(a)(a)\Rightarrow(b)\Rightarrow(c)\Rightarrow(d)\Rightarrow(e)\Rightarrow(f)\Rightarrow(a)

(a)⇒(b)(a)\Rightarrow(b)。 正算子按定义自伴。设 λ\lambda 是特征值,v≠0v\ne 0 是对应特征向量:

0≤⟨Tv,v⟩=⟨λv,v⟩=λ⟨v,v⟩=λ∥v∥20\le\langle Tv,v\rangle=\langle \lambda v,v\rangle=\lambda\langle v,v\rangle=\lambda\|v\|^2

因 ∥v∥2>0\|v\|^2>0,得 λ≥0\lambda\ge 0。

这一行是整个证明的枢纽

它说明 ⟨Tv,v⟩≥0\langle Tv,v\rangle\ge 0 在特征向量上取值就足以约束特征值。
换句话说:定义里「对所有 vv」这个量词看起来很重,实际起作用的信息集中在特征方向上。

(b)⇒(c)(b)\Rightarrow(c)。 由 26 讲的谱定理:自伴算子有一组由特征向量组成的规范正交基 e1,…,ene_1,\dots,e_n。设 Tek=λkekTe_k=\lambda_k e_k,每个 λk≥0\lambda_k\ge 0。于是 TT 关于这组基的矩阵是 diag⁡(λ1,…,λn)\operatorname{diag}(\lambda_1,\dots,\lambda_n)。

(c)⇒(d)(c)\Rightarrow(d)。 关键工具是线性映射引理(08 讲:在基上任意指定像,就唯一确定一个线性映射)。在基 e1,…,ene_1,\dots,e_n 上指定

Rek=λk ek,k=1,…,nRe_k=\sqrt{\lambda_k}\,e_k,\qquad k=1,\dots,n

λk\sqrt{\lambda_k} 有定义,因为 λk≥0\lambda_k\ge 0。这个 RR 是正的,而且不需要回头引用任何一条,直接验即可:它关于这组规范正交基的矩阵是实对角矩阵,故自伴;又对 v=∑akekv=\sum a_ke_k 有

⟨Rv,v⟩=⟨∑akλk ek, ∑akek⟩=∑kλk ∣ak∣2≥0\langle Rv,v\rangle=\Big\langle\sum a_k\sqrt{\lambda_k}\,e_k,\ \sum a_ke_k\Big\rangle=\sum_k\sqrt{\lambda_k}\,|a_k|^2\ge 0

(交叉项全部消失,这正是「正交基」的用处。)所以 RR 正。并且

R2ek=λk⋅λk ek=λkek=TekR^2e_k=\sqrt{\lambda_k}\cdot\sqrt{\lambda_k}\,e_k=\lambda_k e_k=Te_k

对每个基向量成立,所以 R2=TR^2=T。RR 就是 TT 的正平方根。

(d)⇒(e)(d)\Rightarrow(e)。 正平方根本身是正算子,正算子自伴(定义)。所以它是自伴平方根。

(e)⇒(f)(e)\Rightarrow(f)。 若 T=R2T=R^2 且 R∗=RR^*=R,则 T=R2=R∗RT=R^2=R^*R。取的就是这个 RR。

(f)⇒(a)(f)\Rightarrow(a)。 设 T=R∗RT=R^*R。先算自伴性:

T∗=(R∗R)∗=R∗(R∗)∗=R∗R=TT^*=(R^*R)^*=R^*(R^*)^*=R^*R=T

再算内积,这一步只需一次伴随的定义 ⟨R∗u,w⟩=⟨u,Rw⟩\langle R^*u,w\rangle=\langle u,Rw\rangle:

⟨Tv,v⟩=⟨R∗Rv,v⟩=⟨Rv,Rv⟩=∥Rv∥2≥0\langle Tv,v\rangle=\langle R^*Rv,v\rangle=\langle Rv,Rv\rangle=\|Rv\|^2\ge 0

■\blacksquare

这六条怎么记

记住「非负数」的三句话,一一对应:
λ≥0\lambda\ge 0 ↔\leftrightarrow (b)(c) 特征值/对角元非负;
λ\lambda 有非负平方根 ↔\leftrightarrow (d)(e);
λ=zˉz\lambda=\bar z z ↔\leftrightarrow (f) T=R∗RT=R^*R。
数系里的「非负」有几种等价说法,算子层就有几条等价刻画。

3.2 正平方根唯一(存在性的另一半)

3.1 只说了「存在」。而 3.1 的 (d) 里那个 RR 是照着一组预先选好的特征基造出来的——如果换一组基,会不会造出另一个正平方根?答案是:不会。

定理:VV 上每个正算子都有唯一的正平方根。

证明。 设 TT 正,RR 是 TT 的任一正平方根。策略是:证明 RR 在 TT 的每个特征向量上的作用被完全钉死;然后利用「TT 的特征向量构成一组基」,得出 RR 被钉死。

第一步:RR 也是正的,所以它也有一条特征基。 RR 正 ⇒\Rightarrow RR 自伴,由谱定理,VV 有由 RR 的特征向量组成的规范正交基 e1,…,ene_1,\dots,e_n。记 Rek=μkekRe_k=\mu_k e_k。因为 RR 正,μk≥0\mu_k\ge 0。令

λk=μk2≥0\lambda_k=\mu_k^2\ge 0

则 R2ek=μk2ek=λkekR^2e_k=\mu_k^2e_k=\lambda_k e_k,即 Tek=λkekTe_k=\lambda_k e_k。顺带得到:RR 的特征基同时也是 TT 的特征基(这是 T=R2T=R^2 的直接后果)。

第二步:设 vv 是 TT 的任一特征向量,Tv=λvTv=\lambda v,其中 λ≥0\lambda\ge 0(3.1 (a)⇒\Rightarrow(b))。把 vv 按 RR 的特征基展开:

v=a1e1+⋯+anenv=a_1e_1+\cdots+a_ne_n

于是

Rv=a1μ1e1+⋯+anμnen=a1λ1 e1+⋯+anλn en(3.1)Rv=a_1\mu_1e_1+\cdots+a_n\mu_ne_n=a_1\sqrt{\lambda_1}\,e_1+\cdots+a_n\sqrt{\lambda_n}\,e_n \tag{3.1}

另一方面,vv 是 TT 的特征向量,而 T=R2T=R^2:

λv=Tv=R2v=a1λ1e1+⋯+anλnen(3.2)\lambda v=Tv=R^2v=a_1\lambda_1e_1+\cdots+a_n\lambda_ne_n \tag{3.2}

同时,由 v=∑akekv=\sum a_ke_k 又有

λv=a1λe1+⋯+anλen(3.3)\lambda v=a_1\lambda e_1+\cdots+a_n\lambda e_n \tag{3.3}

比较 (3.2) 与 (3.3),并利用 e1,…,ene_1,\dots,e_n 线性无关:

ak(λ−λk)=0,k=1,…,na_k(\lambda-\lambda_k)=0,\qquad k=1,\dots,n

第三步:读出结论。 上式说明:只要 λk≠λ\lambda_k\ne\lambda,系数 aka_k 就必须是 00。所以

v=∑{k: λk=λ}akekv=\sum_{\{k:\ \lambda_k=\lambda\}}a_ke_k

代回 (3.1),这些下标上 λk=λ\sqrt{\lambda_k}=\sqrt{\lambda}:

Rv=∑{k: λk=λ}akλ ek=λ vRv=\sum_{\{k:\ \lambda_k=\lambda\}}a_k\sqrt{\lambda}\,e_k=\sqrt{\lambda}\,v

第四步:拼起来。 我们证明了:对 TT 的每个特征向量 vv(特征值 λ\lambda),任何正平方根 RR 都必须满足 Rv=λ vRv=\sqrt\lambda\,v。而由谱定理,TT 的特征向量张成整个 VV(含一组基)。两个线性映射若在一组基上取值相同,则处处相同。故 RR 唯一。■\blacksquare

这个证明的骨架,值得单独记

「在特征向量上确定 ⇒\Rightarrow 在基上确定 ⇒\Rightarrow 处处确定」。
这是一个反复出现的套路:要证某个算子唯一,别看它在任意 vv 上怎样,
先看它在一组基上被逼成什么。3.1 的 (c)⇒\Rightarrow(d) 是反方向——在基上指定像来造算子;
这里是在基上读出像来定算子。同一枚硬币的两面。

「唯一」的唯一性,只对「正的」那一个成立

一个正算子可以有无穷多个平方根,正平方根只是其中一个。
最简单的例子:V=R2V=\mathbb{R}^2 上的恒等算子 II。取

R=(100−1)R=\begin{pmatrix}1&0\\ 0&-1\end{pmatrix}

则 R2=IR^2=I,所以 RR 是 II 的平方根,但 RR 不正(⟨R(0,1),(0,1)⟩=−1\langle R(0,1),(0,1)\rangle=-1)。
更有甚者,Rθ=(cos⁡θ−sin⁡θsin⁡θcos⁡θ)R_\theta=\begin{pmatrix}\cos\theta&-\sin\theta\\ \sin\theta&\cos\theta\end{pmatrix}(旋转 θ\theta)都满足 Rθ2=IR_\theta^2=I,当 θ=π\theta=\pi 时是 −I-I。
但正平方根只有一个:II 自己(因为 II 的特征值全是 11,1=1\sqrt1=1)。

3.3 记号 T\sqrt{T}

有了唯一性,下面这个记号才站得住:

对正算子 TT,T\sqrt{T} 表示 TT 的那个正平方根。

若无唯一性,写 T\sqrt{T} 就是含糊的——读者不知道你说的是哪一个。数学里「先证唯一,再定记号」是标准动作,这里是它的一个干净样本。

3.4 T∗TT^*T 永远是正的

这条重要到值得单列。设 T∈L(V,W)T\in\mathcal{L}(V,W)(允许 V≠WV\ne W、允许不可逆、允许不是方阵),T∗∈L(W,V)T^*\in\mathcal{L}(W,V)。那么 T∗T∈L(V)T^*T\in\mathcal{L}(V),且

(T∗T)∗=T∗T,⟨T∗Tv,v⟩=⟨Tv,Tv⟩=∥Tv∥2≥0(T^*T)^*=T^*T,\qquad \langle T^*Tv,v\rangle=\langle Tv,Tv\rangle=\|Tv\|^2\ge 0

所以 T∗TT^*T 总是正算子。这条不需要任何假设。

它还给出一个常用推论:T∗TT^*T 的零空间就是 TT 的零空间:

v∈null⁡(T∗T)  ⟺  ∥Tv∥2=0  ⟺  Tv=0v\in\operatorname{null}(T^*T)\iff\|Tv\|^2=0\iff Tv=0

(一般的 TT 不满足 null⁡T∗=null⁡T\operatorname{null}T^*=\operatorname{null}T,但 T∗TT^*T 与 TT 的零空间永远相同。)

为什么 SVD 从这里长出来

因为 T∗TT^*T 正,所以能取 T∗T\sqrt{T^*T},记作 ∣T∣|T|。
∣T∣|T| 的特征值就是 TT 的奇异值,∣T∣|T| 的特征向量就是右奇异向量。
29 讲会把这条线索走完。现在只要记住:SVD 的地基是「T∗TT^*T 正」。

3.5 一个「不涉及平方根、但靠平方根证」的结论

定理:设 TT 是 VV 上的正算子,v∈Vv\in V 且 ⟨Tv,v⟩=0\langle Tv,v\rangle=0,则 Tv=0Tv=0。

证明只有两行,但它用到了 3.1 的 (d):

0=⟨Tv,v⟩=⟨TTv,v⟩=⟨Tv,Tv⟩=∥Tv∥2\begin{aligned}0&=\langle Tv,v\rangle=\langle \sqrt{T}\sqrt{T}v,v\rangle=\langle \sqrt{T}v,\sqrt{T}v\rangle=\big\|\sqrt{T}v\big\|^2\end{aligned}

第二个等号把 TT 拆成 T⋅T\sqrt T\cdot\sqrt T,第三个等号把 T\sqrt T 从内积左边搬到右边——这正是「T\sqrt T 自伴」的作用。于是 Tv=0\sqrt Tv=0,进而

Tv=T(Tv)=T(0)=0Tv=\sqrt T\big(\sqrt Tv\big)=\sqrt T(0)=0

■\blacksquare

这条结论**只对正算子**成立

取 T=diag⁡(1,−1)T=\operatorname{diag}(1,-1),v=(1,1)v=(1,1),则 ⟨Tv,v⟩=1−1=0\langle Tv,v\rangle=1-1=0,但 Tv=(1,−1)≠0Tv=(1,-1)\ne 0。
差别在哪一行?第二行:T\sqrt T 不存在(TT 不正),拆不出那个平方。
证明的每一步都要问「这一步凭什么」,这里的凭据就是正性。

3.6 手算:把 T\sqrt{T} 真的算出来

例 A:对角矩阵。

S=(1002),S=(1002)S=\begin{pmatrix}1&0\\0&2\end{pmatrix},\qquad \sqrt S=\begin{pmatrix}1&0\\0&\sqrt2\end{pmatrix}

因为 SS 已经是对角的,特征值 1,21,2 非负,逐项开方即可。这是 3.1 (c)⇒\Rightarrow(d) 的最直白形态。

例 B:秩一矩阵。

T=(1111),⟨T(x,y),(x,y)⟩=x2+2xy+y2=(x+y)2≥0T=\begin{pmatrix}1&1\\1&1\end{pmatrix},\qquad \langle T(x,y),(x,y)\rangle=x^2+2xy+y^2=(x+y)^2\ge 0

自伴、二次型是平方,故正。特征值:tr⁡=2\operatorname{tr}=2,det⁡=0\det=0,故为 22 与 00;特征向量是 (12,12)(\tfrac1{\sqrt2},\tfrac1{\sqrt2})(对应 22)与 (12,−12)(\tfrac1{\sqrt2},-\tfrac1{\sqrt2})(对应 00)。逐项开方得特征值 2\sqrt2 与 00,在同一组特征向量上:

T=(12121212)\sqrt T=\begin{pmatrix}\tfrac1{\sqrt2}&\tfrac1{\sqrt2}\\[2pt]\tfrac1{\sqrt2}&\tfrac1{\sqrt2}\end{pmatrix}

验算:这个矩阵每行每列和分别是 2\sqrt2、2\sqrt2,平方后恰为 (1111)\begin{pmatrix}1&1\\1&1\end{pmatrix}。✓

例 C:不试凑、直接解方程。

A=(5445)A=\begin{pmatrix}5&4\\4&5\end{pmatrix}

它对称,特征值 99 与 11(tr⁡=10\operatorname{tr}=10,det⁡=25−16=9\det=25-16=9),都非负,故正。设正平方根也是对称的(为什么可以这样假设?因为 A\sqrt A 与 AA 共享特征基,见 3.2 第一步),记

A=(abbd)\sqrt A=\begin{pmatrix}a&b\\b&d\end{pmatrix}

平方必须等于 AA:

(a2+b2b(a+d)b(a+d)b2+d2)=(5445)\begin{pmatrix}a^2+b^2&b(a+d)\\ b(a+d)&b^2+d^2\end{pmatrix}=\begin{pmatrix}5&4\\4&5\end{pmatrix}

于是

{a2+b2=5b(a+d)=4b2+d2=5\begin{cases}a^2+b^2=5\\ b(a+d)=4\\ b^2+d^2=5\end{cases}

由第一个与第三个得 a2=d2a^2=d^2,即 a=da=d 或 a=−da=-d。

  • 取 a=−da=-d:则 b(a+d)=0≠4b(a+d)=0\ne 4,矛盾。
  • 取 a=da=d:2ab=42ab=4,即 ab=2ab=2;配合 a2+b2=5a^2+b^2=5 得 (a+b)2=9(a+b)^2=9,(a−b)2=1(a-b)^2=1。

两组解:(a,b)=(2,1)(a,b)=(2,1) 或 (−2,−1)(-2,-1)。前者给出

A=(2112)\sqrt A=\begin{pmatrix}2&1\\1&2\end{pmatrix}

验算 (2112)2=(5445)\begin{pmatrix}2&1\\1&2\end{pmatrix}^2=\begin{pmatrix}5&4\\4&5\end{pmatrix} ✓。后者给出 −A-\sqrt A,它是平方根但不正(特征值 −3,−1-3,-1)。

例 C 的两个收获

① 解方程时出现两组解是正常的——平方根从来不唯一,唯一的是「正」的那一个。
② 「要正」这个条件在这里等价于「特征值取 9=3\sqrt9=3 与 1=1\sqrt1=1 而非 −3,−1-3,-1」,
也就是取非负的那一支。逐项开方时永远取非负根,这就是唯一性的来源。


四、几何直观与记忆钩子

视角一:二次型是「读数的表盘」

⟨Tv,v⟩\langle Tv,v\rangle 可以读成:把 vv 送进 TT,再和原来的 vv 做内积——相当于用 vv 去「量」TvTv 在 vv 方向上的分量。

  • 若这个读数在任何方向都不为负,说明 TT 在所有方向上都不「反向拉扯」空间,只可能压缩、拉伸、或压扁。
  • 允许取 00:说明某些方向被压成零(不可逆)。这正是「半正定」的几何含义——碗底是平的,不是一个尖点。

视角二:正算子 = 非负伸缩的正交组合

由谱定理,正算子总可以写成

T=λ1P1+λ2P2+⋯+λmPm,λk≥0T=\lambda_1P_1+\lambda_2P_2+\cdots+\lambda_mP_m,\qquad \lambda_k\ge 0

其中 PkP_k 是到各特征空间的正交投影(互不相同、两两正交、和为 II)。于是

T=λ1P1+⋯+λmPm\sqrt T=\sqrt{\lambda_1}P_1+\cdots+\sqrt{\lambda_m}P_m

开方就是「沿各个正交方向的伸缩倍数逐项开方,方向一个都不动」。 这是本讲最值得带走的一张图:

方向不动,只把每个方向的倍数开方\text{方向不动,只把每个方向的倍数开方}

视角三:与复数的完整对照

建立这张对照表,本讲的定义就都「有根」了:

复数 c∈Cc\in\mathbb{C} 算子 T∈L(V)T\in\mathcal{L}(V)
c≥0c\ge 0(非负实数) TT 正:T=T∗T=T^* 且 ⟨Tv,v⟩≥0\langle Tv,v\rangle\ge 0
cc 是实数 TT 自伴(25 讲)
cc 有非负平方根 TT 有正平方根
cc 有实平方根 TT 有自伴平方根
c=zˉz=∣z∣2c=\bar zz=|z|^2 T=R∗RT=R^*R
c\sqrt{c} 唯一(取非负根) T\sqrt T 唯一
∣z∣2=zˉz|z|^2=\bar zz ∥Tv∥2=⟨T∗Tv,v⟩\|Tv\|^2=\langle T^*Tv,v\rangle
为什么这张表这么整齐

因为「算子」在谱定理下就是「一组数(特征值)+ 一组正交方向」。
一旦方向已经正交、互不干扰,所有逐点的运算都退化成对每个特征值的标量运算。
所以实数里关于「非负」的每一条性质,只要不涉及大小顺序,都能一字不改地搬过来。
反过来,凡是涉及「谁大谁小」的(比如 a+b≤a+b\sqrt{a+b}\le\sqrt a+\sqrt b),搬运时就要小心——
算子的序不是全序。

记忆钩子

正算子之于自伴算子,如同非负数之于实数,如同半正定矩阵之于对称矩阵。
三句话是同一句话,只是词汇表不同。


五、常见误区

误区 1:以为「元素全为正的对称矩阵」就是正算子

反例:A=(1221)A=\begin{pmatrix}1&2\\2&1\end{pmatrix}。它对称、四个元素全为正,但特征值是 33 与 −1-1。
取 v=(1,−1)v=(1,-1):Av=(1−2, 2−1)=(−1,1)Av=(1-2,\ 2-1)=(-1,1),于是 ⟨Av,v⟩=(−1)⋅1+1⋅(−1)=−2<0\langle Av,v\rangle=(-1)\cdot1+1\cdot(-1)=-2<0。
判据是特征值(或二次型的符号),不是元素的符号。
正确做法只有三条:查特征值全 ≥0\ge0;或配方成平方和;或找出 RR 使 A=R∗RA=R^*R。

误区 2:在实空间里丢掉「自伴」这一条

取 T=(01−10)T=\begin{pmatrix}0&1\\-1&0\end{pmatrix}(逆时针转 90°90° 再……实际是旋转)。对任意 (x,y)(x,y):

⟨T(x,y),(x,y)⟩=yx+(−x)y=0≥0\langle T(x,y),(x,y)\rangle=yx+(-x)y=0\ge 0

条件 2 完全成立,但 TT=(0−110)≠TT^{\mathsf T}=\begin{pmatrix}0&-1\\1&0\end{pmatrix}\ne T,TT 不自伴。
若按「⟨Tv,v⟩≥0\langle Tv,v\rangle\ge0 就算正」的错定义,这个旋转就成了正算子,而它的特征值是 ±i\pm i——
连实特征值都没有,与 3.1(b) 直接冲突。
复空间才允许省掉自伴(因为那时「⟨Tv,v⟩\langle Tv,v\rangle 恒为实数」已蕴含自伴)。

误区 3:以为「有平方根」就能推出「是正的」

反例来自平方根本身的构造:T(z1,z2,z3)=(z3,0,0)T(z_1,z_2,z_3)=(z_3,0,0),取 R(z1,z2,z3)=(z2,z3,0)R(z_1,z_2,z_3)=(z_2,z_3,0),
则 R2=TR^2=T,TT 确实有平方根。但 TT 不正:取 v=(1,0,1)v=(1,0,1),
先取 v=(1,0,1)v=(1,0,1),得 ⟨Tv,v⟩=⟨(1,0,0),(1,0,1)⟩=1\langle Tv,v\rangle=\langle(1,0,0),(1,0,1)\rangle=1;再取 v=(1,0,i)v=(1,0,i),得

⟨Tv,v⟩=⟨(1,0,0),(1,0,i)⟩=i⋅1‾=i\langle Tv,v\rangle=\langle(1,0,0),(1,0,i)\rangle=i\cdot\overline{1}=i

取到了非实数,连自伴都不成立。
取到了非实数,自伴都不成立。
正确的方向是 (d)(e):有「正的」或「自伴的」平方根才等价于正。
只是「有平方根」,谁都可以有。

误区 4:以为「平方根唯一」

唯一的是正平方根。恒等算子 II 在 dim⁡V>1\dim V>1 时有无穷多个平方根
(一切满足 R2=IR^2=I 的算子,包括无数组旋转与反射,见 3.2 的警告框)。
所以看到 T\sqrt{T} 这个记号,心里要补一句「TT 是正的,且取的是正的那一支」。

误区 5:把「⟨Tv,v⟩=0\langle Tv,v\rangle=0 于某个 vv」当成「Tv=0Tv=0」

只对正算子成立(3.5)。对一般自伴算子不成立:diag⁡(1,−1)\operatorname{diag}(1,-1) 与 v=(1,1)v=(1,1) 就是反例。
这类「在特征基下正负相消」的构造,是考试题最爱埋的坑。

误区 6:以为正算子的序是「全序」

S,TS,T 都正,不代表 S−TS-T 或 T−ST-S 中有一个正——两者可能都「不定」。
类比:矩阵没有「谁大谁小」,只有「S−TS-T 非负」这种偏序关系。
所以写 S≤T\sqrt S\le\sqrt T 之前,先确认你指的是「T−S\sqrt T-\sqrt S 正」,
并且不要把实数的运算律(如 S+T=S+T\sqrt{S+T}=\sqrt S+\sqrt T)照搬过来——那是错的。


六、拓展(可跳过)

① T\sqrt T 其实是 TT 的多项式。 对正算子 TT,存在实系数多项式 pp 使 T=p(T)\sqrt T=p(T)(pp 可以取为在 TT 的谱上插值   ⋅  \sqrt{\;\cdot\;} 的那个多项式)。这条说明 T\sqrt T 没有「跳出 TT 生成的代数」——它对数值计算意义很大:不必求特征分解,用多项式(或级数)近似即可。

② 正算子之和、正整数次幂仍正。 若 S,TS,T 正,则 S+TS+T 正(⟨(S+T)v,v⟩=⟨Sv,v⟩+⟨Tv,v⟩≥0\langle(S+T)v,v\rangle=\langle Sv,v\rangle+\langle Tv,v\rangle\ge0,自伴显然)。若 TT 正,TkT^k 正(kk 为正整数):TkT^k 自伴,且 Tk=(Tk/2)∗Tk/2T^k=\big(T^{k/2}\big)^*T^{k/2}(kk 偶)或 =T(k−1)/2 T T(k−1)/2=T^{(k-1)/2}\,T\,T^{(k-1)/2}(kk 奇),形如 R∗RR^*R。

③ STST 何时正:可交换是关键。 若 S,TS,T 都正,则

ST 正  ⟺  ST=TSST\ \text{正}\iff ST=TS

理由:STST 自伴   ⟺  (ST)∗=ST  ⟺  T∗S∗=ST  ⟺  TS=ST\iff (ST)^*=ST\iff T^*S^*=ST\iff TS=ST。有了可交换,S,TS,T 能同时用一组规范正交基对角化(26 讲的推论),STST 的对角元就是 λkμk≥0\lambda_k\mu_k\ge0。不可交换时 STST 连自伴都不是。

④ 一个具体判据:T−αIT-\alpha I。 设 TT 自伴,α∈R\alpha\in\mathbb{R}。则

T−αI 正  ⟺  α≤T 的每个特征值T-\alpha I\ \text{正}\iff\alpha\le T\ \text{的每个特征值}

因为 T−αIT-\alpha I 的特征值是 λk−α\lambda_k-\alpha。由此得到:「αI≤T\alpha I\le T」这种写法在算子层是有确切含义的,它就是「T−αIT-\alpha I 正」。

⑤ 一个让人意外的正算子:−d2dx2-\dfrac{d^2}{dx^2}。 在「二阶可导、满足某组边值条件」的函数空间上,取内积 ⟨f,g⟩=∫fgˉ\langle f,g\rangle=\int f\bar g,分部积分两次给出

⟨−d2fdx2,f⟩=∫∣f′∣2≥0\Big\langle -\tfrac{d^2f}{dx^2},f\Big\rangle=\int|f'|^2\ge 0

所以负的二阶导数(配上合适的边界条件)是正算子。这一条把「正算子」与物理里的能量、与 Sturm–Liouville 理论接上了:能量算子的正性,就是系统稳定性的代数表达。

⑥ 工程里的名字。 数值线代里正算子对应「对称正定 / 半正定矩阵」。Cholesky 分解 A=LL∗A=LL^* 就是 A=R∗RA=R^*R 的计算版;判断它是 AA 是否正;协方差矩阵必须正(否则方差为负),这是「协方差矩阵不合法」报警的来源。


七、自测(先做,再看答案)

1. 证明:若 TT 是正算子,则 TT 的每个特征值都是非负实数;并举出一个自伴但不正的 2×22\times2 例子。

2. 设 T∈L(F3)T\in\mathcal{L}(\mathbb{F}^3) 定义为 T(z1,z2,z3)=(z3,0,0)T(z_1,z_2,z_3)=(z_3,0,0)。证明 TT 有平方根,但 TT 不是正的。

3. 求 A=(106610)A=\begin{pmatrix}10&6\\6&10\end{pmatrix} 的正平方根。

4. 设 TT 正。证明 null⁡T=null⁡T\operatorname{null}\sqrt T=\operatorname{null}T。

5. 判断 B=(2332)B=\begin{pmatrix}2&3\\3&2\end{pmatrix} 是否正算子;若不是,给出一个具体的 vv 使 ⟨Bv,v⟩<0\langle Bv,v\rangle<0。

答案(点开)

1. 设 Tv=λvTv=\lambda v,v≠0v\ne0。则 0≤⟨Tv,v⟩=λ⟨v,v⟩=λ∥v∥20\le\langle Tv,v\rangle=\lambda\langle v,v\rangle=\lambda\|v\|^2,
故 λ≥0\lambda\ge0;又 TT 自伴故 λ∈R\lambda\in\mathbb{R}。
反例:(100−1)\begin{pmatrix}1&0\\0&-1\end{pmatrix},自伴、特征值 1,−11,-1,取 v=(0,1)v=(0,1) 得内积 −1-1。

2. 取 R(z1,z2,z3)=(z2,z3,0)R(z_1,z_2,z_3)=(z_2,z_3,0),则 R2(z1,z2,z3)=R(z2,z3,0)=(z3,0,0)=TR^2(z_1,z_2,z_3)=R(z_2,z_3,0)=(z_3,0,0)=T,故 TT 有平方根。
但 TT 不自伴:关于标准基,TT 的矩阵是 (001000000)\begin{pmatrix}0&0&1\\0&0&0\\0&0&0\end{pmatrix},
其伴随(转置)是 (000000100)≠T\begin{pmatrix}0&0&0\\0&0&0\\1&0&0\end{pmatrix}\ne T。
所以 TT 不正。(不自伴已经足够否定,不必再去算内积。)
这正对应 3.1 中「(d) 要求平方根为正」——光有平方根不够。

3. AA 对称,tr⁡=20\operatorname{tr}=20,det⁡=100−36=64\det=100-36=64,特征值 λ2−20λ+64=0\lambda^2-20\lambda+64=0 得 λ=16,4\lambda=16,4,均非负。
设 A=(abbd)\sqrt A=\begin{pmatrix}a&b\\b&d\end{pmatrix},则 a2+b2=10a^2+b^2=10、b(a+d)=6b(a+d)=6、b2+d2=10b^2+d^2=10。
由首末两式 a=±da=\pm d;a=−da=-d 与第二式矛盾,故 a=da=d,于是 ab=3ab=3、a2+b2=10a^2+b^2=10,
得 (a+b)2=16(a+b)^2=16,(a−b)2=4(a-b)^2=4,取非负支:a=3,b=1a=3,b=1。

A=(3113),(3113)2=(106610) ✓\sqrt A=\begin{pmatrix}3&1\\1&3\end{pmatrix},\qquad \begin{pmatrix}3&1\\1&3\end{pmatrix}^2=\begin{pmatrix}10&6\\6&10\end{pmatrix}\ ✓

4. 「⊆\subseteq」:若 Tv=0\sqrt Tv=0,则 Tv=T(Tv)=0Tv=\sqrt T(\sqrt Tv)=0。
「⊇\supseteq」:若 Tv=0Tv=0,则 ⟨Tv,v⟩=0\langle Tv,v\rangle=0,由 3.5 的证明 ∥Tv∥2=⟨Tv,v⟩=0\|\sqrt Tv\|^2=\langle Tv,v\rangle=0,故 Tv=0\sqrt Tv=0。
(注意这里是「T\sqrt T 自伴 + 3.5」的直接应用,不必再做谱分解。)

5. 不正。tr⁡=4\operatorname{tr}=4,det⁡=4−9=−5\det=4-9=-5,特征值满足 λ2−4λ−5=0\lambda^2-4\lambda-5=0,即 55 与 −1-1。
取 v=(1,−1)v=(1,-1):Bv=(2−3, 3−2)=(−1,1)Bv=(2-3,\ 3-2)=(-1,1),⟨Bv,v⟩=(−1)(1)+(1)(−1)=−2<0\langle Bv,v\rangle=(-1)(1)+(1)(-1)=-2<0。
注意它对称且元素全为正,这正是误区 1 的陷阱。


八、复习卡

项 内容
定义 TT 正   ⟺  \iff T=T∗T=T^* 且 ∀v: ⟨Tv,v⟩≥0\forall v:\ \langle Tv,v\rangle\ge0
复空间 「自伴」可省(⟨Tv,v⟩\langle Tv,v\rangle 恒实时自动自伴);实空间不可省
术语 又叫「非负算子」「正半定算子」;「正定」是更强的 ⟨Tv,v⟩>0 (v≠0)\langle Tv,v\rangle>0\ (v\ne0)
最小例子 00;II;cI (c≥0)cI\ (c\ge0);正交投影 PUP_U
(b) 特征值 正   ⟺  \iff 自伴且特征值全 ≥0\ge0
(c) 对角化 正   ⟺  \iff 某规范正交基下矩阵为对角且对角元非负
(d) 正平方根 正   ⟺  \iff 存在正算子 RR:R2=TR^2=T
(e) 自伴平方根 正   ⟺  \iff 存在自伴 RR:T=R2T=R^2
(f) 乘积形式 正   ⟺  \iff 存在 RR:T=R∗RT=R^*R
唯一性 正平方根唯一;T\sqrt T 专指它
唯一性证明骨架 在 TT 的特征向量上 RR 被逼成 λ v\sqrt\lambda\,v;特征基 ⇒\Rightarrow 处处确定
T∗TT^*T 对任意 T:V→WT:V\to W,T∗TT^*T 恒正;∥Tv∥2=⟨T∗Tv,v⟩\|Tv\|^2=\langle T^*Tv,v\rangle
零空间 null⁡(T∗T)=null⁡T\operatorname{null}(T^*T)=\operatorname{null}T;null⁡T=null⁡T\operatorname{null}\sqrt T=\operatorname{null}T
关键小引理 正算子 + ⟨Tv,v⟩=0⇒Tv=0\langle Tv,v\rangle=0\Rightarrow Tv=0(用 T\sqrt T 拆)
可视化 T\sqrt T:各正交方向的伸缩倍数逐项开方,方向不动
判据 特征值全 ≥0\ge0(对称阵);或配方成平方和;或找 R∗RR^*R
常见陷阱 元素全正 ≠\ne 正;有平方根 ≠\ne 正;平方根不唯一

九、出处与延伸(可选,不影响阅读)

  • 主线:这一讲的框架来自《线性代数应该这样学》(第 4 版)第 7 章 C 节「正算子」,对应 PDF 页码 223–227、印刷页码 209–213。该节的编排是:定义 → 三个例子 → 平方根定义与例子 → 六条等价刻画 → 正平方根唯一 → 记号 → 例 → 「⟨Tv,v⟩=0⇒Tv=0\langle Tv,v\rangle=0\Rightarrow Tv=0」→ 习题。
  • 定义与编号对照:正算子定义 7.34;例子 7.35;平方根定义 7.36、例 7.37;六条等价刻画 7.38;唯一性 7.39;记号 7.40;例 7.41;引理 7.43。本讲把这些合并重排成「定义 → 例子 → 主定理 → 唯一性 → 应用」,并把证明拆成更细的步。
  • 术语别名:正算子 = positive operator;有些文献写「正半定算子 / positive semidefinite operator」,含义相同。「正定算子」通常指更强的情形(⟨Tv,v⟩>0\langle Tv,v\rangle>0 对一切 v≠0v\ne0),与之对应的是「半正定」。
  • 为什么叫「正」而不是「非负」:这是一个历史遗留的用词。正算子对应的是非负数而非正数,严格说术语名不副实;算子理论界沿用了「正」这个叫法,本课也随大流。读到「正的」时,脑子里请自动替换为「非负的」。
  • 想深挖时可以查:
    • 算子单调性与 Löwner 序(  \sqrt{\;} 在正算子序下单调);
    • 矩阵函数与谱映射(T=p(T)\sqrt T=p(T)、eTe^{T}、log⁡T\log T);
    • Cholesky 分解与半正定判定(数值线代教材,如 Golub & Van Loan 第 4 章);
    • 量子力学中的可观测量:自伴算子;密度矩阵:迹为 11 的正算子——这是「正算子」在物理里最重要的一次出场。
  • 与后续讲次的接口:28 讲用 T∗T\sqrt{T^*T} 做极分解;29 讲用它定义奇异值;35 讲把 ⟨Tv,v⟩\langle Tv,v\rangle 直接看成二次型,与本讲的符号判据完全重合。