这讲讲什么 :把「对角元之和」这个看起来最依赖坐标的数,证明成算子本身的不变量;再由这一条免费换来三个结论——迹等于特征值之和、tr ( S T ) = tr ( T S ) \operatorname{tr}(ST)=\operatorname{tr}(TS) tr ( S T ) = tr ( T S ) 、以及「S T − T S = I ST-TS=I S T − T S = I 不可能」。前置 :19 讲(上三角矩阵:对角元就是全部特征值)、17 讲(特征值与重数)、20 讲(可对角化判据,自测第 5 题用)、10 讲(换基公式 A = C − 1 B C A=C^{-1}BC A = C − 1 B C )、22–23 讲(内积与规范正交基)一句话 :迹是「既线性、又循环不变、又把恒等算子映成维数」的唯一 泛函,所以它必然等于特征值之和。
一、为什么需要它 #
1.1 一个便宜的数,能替你省下大量验算 #
求特征值有多贵,你已经有体会:n n n 阶矩阵要解一个 n n n 次方程,n ≥ 5 n\ge 5 n ≥ 5 时连根式解都没有。19 讲那条「上三角之后对角元即特征值」是好消息,但它只在你已经 拿到上三角形状时才免费。
现在换个角度问:有没有一个数,能让我在算完特征值之后一秒钟自查?
有。把对角线上的数加起来。因为第三节要证的定理会说:
tr T = λ 1 + λ 2 + ⋯ + λ n ( 各特征值按其重数出现 ) (34.1) \operatorname{tr}T=\lambda_1+\lambda_2+\cdots+\lambda_n\qquad(\text{各特征值按其重数出现}) \tag{34.1} tr T = λ 1 + λ 2 + ⋯ + λ n ( 各特征值按其重数出现 ) ( 34.1 )
你辛苦求出三个特征值,加起来对不上对角元之和,那就一定算错了——不是「可能有错」,是一定 错。这种一行字的校验,比再算一遍特征多项式便宜太多。
1.2 更深的理由:基无关的量太稀有 #
回想你手里已有的算子不变量。可逆性算一个,但它只给「是/否」。秩算一个,但它是整数。要一个数 ,而且要对所有算子都有定义、还随算子连续变化——这样的东西非常少。
原因不难理解:算子换一组基,矩阵变成 C − 1 B C C^{-1}BC C − 1 B C ,各个元素全被打乱。想让某个「元素算出来的数」在被打乱后保持不变,等于要求这个数对所有 可逆 C C C 都不动。这是一条极强的约束,绝大多数「元素的和」都过不了。
迹过了。行列式也过了。这两个恰好是特征多项式最高次的两个系数 (第 3.5 节会说清)。所以它们组成一对,不是巧合。
1.3 还有一层:迹是「被迫存在」的 #
设 τ \tau τ 是从 L ( V ) \mathcal{L}(V) L ( V ) 到 F \mathbb{F} F 的线性泛函,满足两条:
τ ( S T ) = τ ( T S ) 对一切 S , T ; τ ( I ) = dim V (34.2) \tau(ST)=\tau(TS)\quad\text{对一切 }S,T;\qquad \tau(I)=\dim V \tag{34.2} τ ( S T ) = τ ( T S ) 对一切 S , T ; τ ( I ) = dim V ( 34.2 )
结论是:这样的 τ \tau τ 只有一个 ,就是迹。这件事的意义不只是「唯一性漂亮」,而是:只要你想要一个「不管乘法顺序」的线性读数,你手里就自动只剩下迹,没有别的选择。
下面这条随之而来的结论,是本讲最值得记住的应用——它证明了一个「不存在」 :
在有限维空间上,找不到算子 S , T S,T S , T 使 S T − T S = I ST-TS=I S T − T S = I 。
证明只有一行:左边取迹是 0 0 0 ,右边取迹是 dim V ≠ 0 \dim V\ne 0 dim V = 0 。一个等式两边迹不同,等式自然不成立。这种「用不变量否证存在性」的手法,在数学里叫不变量论证 ,迹是它最干净的例子。
读完这一讲你应该能回答:
对角元之和明明随基而变,凭什么说它是算子的属性?
tr ( A B ) = tr ( B A ) \operatorname{tr}(AB)=\operatorname{tr}(BA) tr ( A B ) = tr ( B A ) 到底允许 A , B A,B A , B 不同形状吗?允许换几次顺序?
为什么「迹等于特征值之和」必须假设在复空间上(或特征多项式完全分解)?
迹是 L ( V ) \mathcal{L}(V) L ( V ) 上的线性泛函——这句话能推出什么具体的算法?
二、定义与最小例子 #
2.1 矩阵的迹 #
定义 :设 A A A 是元素属于 F \mathbb{F} F 的 n × n n\times n n × n 方阵。A A A 的迹 记作 tr A \operatorname{tr}A tr A ,定义为对角线上各元素之和:
tr A = ∑ j = 1 n A j , j = A 1 , 1 + A 2 , 2 + ⋯ + A n , n (34.3) \operatorname{tr}A=\sum_{j=1}^{n}A_{j,j}=A_{1,1}+A_{2,2}+\cdots+A_{n,n} \tag{34.3} tr A = j = 1 ∑ n A j , j = A 1 , 1 + A 2 , 2 + ⋯ + A n , n ( 34.3 )
字面意思就是「trace,痕迹」:只看从左上角走到右下角那一条线上留下了什么。
对角元要求行号等于列号,只有方阵才有完整的一套。2 × 3 2\times 3 2 × 3 矩阵谈不上「对角元之和」。 但要留意一个细节:第 3.1 节的 tr ( A B ) = tr ( B A ) \operatorname{tr}(AB)=\operatorname{tr}(BA) tr ( A B ) = tr ( B A ) 里,A A A 是 m × n m\times n m × n 、B B B 是 n × m n\times m n × m ,两个都不是方阵 (除非 m = n m=n m = n ),而乘积 A B AB A B 是 m × m m\times m m × m 、B A BA B A 是 n × n n\times n n × n ,两个都是方的。 「迹只对方阵有定义」和「迹公式里可以出现非方阵」并不冲突。
2.2 最小例子:一个 3 × 3 3\times 3 3 × 3 #
A = ( 3 − 1 − 2 5 2 − 3 1 6 0 ) A=\begin{pmatrix}3&-1&-2\\ 5&2&-3\\ 1&6&0\end{pmatrix} A = 3 5 1 − 1 2 6 − 2 − 3 0
对角线是 3 , 2 , 0 3,2,0 3 , 2 , 0 ,所以
tr A = 3 + 2 + 0 = 5 \operatorname{tr}A=3+2+0=5 tr A = 3 + 2 + 0 = 5
就这么简单。真正的内容全在后面:这个 5 5 5 不随基改变。
2.3 算子的迹:先记账,再还账 #
我们真正想要的是算子 的迹,而不是某个矩阵的迹。自然的想法是:取 V V V 的一组基,算 T T T 关于这组基的矩阵,再取迹。
tr T : = tr M ( T , ( v 1 , … , v n ) ) \operatorname{tr}T\;:=\;\operatorname{tr}\mathcal{M}\big(T,(v_1,\dots,v_n)\big) tr T := tr M ( T , ( v 1 , … , v n ) )
但这个写法有个漏洞:换一组基,右边会不会变? 不变,这个定义才立得住;会变,那它只是「某个矩阵的迹」,不配叫 T T T 的迹。
答案是不变(第 3.2 节证)。所以定义是合法的,而且在定义里可以写「v 1 , … , v n v_1,\dots,v_n v 1 , … , v n 是 V V V 的任意 一组基」——「任意」这两个字,是靠一条定理买来的,不是许愿。
2.4 一个 2 × 2 2\times 2 2 × 2 的全集:四种典型算子 #
把 R 2 \mathbb{R}^2 R 2 上最简单的几个算子排一排,先建立手感(都为标准基下的矩阵):
算子
矩阵
对角元
迹
特征值
恒等 I I I
( 1 0 0 1 ) \begin{pmatrix}1&0\\ 0&1\end{pmatrix} ( 1 0 0 1 )
1 , 1 1,1 1 , 1
2 2 2
1 , 1 1,1 1 , 1
旋转 R θ R_\theta R θ
( cos θ − sin θ sin θ cos θ ) \begin{pmatrix}\cos\theta&-\sin\theta\\ \sin\theta&\cos\theta\end{pmatrix} ( cos θ sin θ − sin θ cos θ )
cos θ , cos θ \cos\theta,\cos\theta cos θ , cos θ
2 cos θ 2\cos\theta 2 cos θ
e i θ , e − i θ \mathrm{e}^{\mathrm{i}\theta},\mathrm{e}^{-\mathrm{i}\theta} e i θ , e − i θ (和 = 2 cos θ =2\cos\theta = 2 cos θ )
投影到 x x x 轴
( 1 0 0 0 ) \begin{pmatrix}1&0\\ 0&0\end{pmatrix} ( 1 0 0 0 )
1 , 0 1,0 1 , 0
1 1 1
1 , 0 1,0 1 , 0
幂零 N N N
( 0 1 0 0 ) \begin{pmatrix}0&1\\ 0&0\end{pmatrix} ( 0 0 1 0 )
0 , 0 0,0 0 , 0
0 0 0
0 , 0 0,0 0 , 0
四行对照着看,能读出三件事:
旋转的迹是 2 cos θ 2\cos\theta 2 cos θ 。 复特征值成对出现 e ± i θ \mathrm{e}^{\pm\mathrm{i}\theta} e ± i θ ,和是实数 2 cos θ 2\cos\theta 2 cos θ ,正好等于对角元之和。这条关系在物理里到处冒头。
投影的迹是 1 1 1 ,等于它「留下的维数」(像是 1 1 1 维,零空间 1 1 1 维)。这个观察第七节自测第 5 题会升级成定理。
幂零矩阵的迹是 0 0 0 ,但它不是零算子。 所以「迹为零」远不等于「算子为零」。这一点第 3.10 节会给出一个例外情形(正算子),对照着记。
2.5 一个复空间的例子 #
定义 T ∈ L ( C 3 ) T\in\mathcal{L}(\mathbb{C}^3) T ∈ L ( C 3 ) 为
T ( z 1 , z 2 , z 3 ) = ( 3 z 1 − z 2 − 2 z 3 , 3 z 1 + 2 z 2 − 3 z 3 , z 1 + 2 z 2 ) T(z_1,z_2,z_3)=(3z_1-z_2-2z_3,\ 3z_1+2z_2-3z_3,\ z_1+2z_2) T ( z 1 , z 2 , z 3 ) = ( 3 z 1 − z 2 − 2 z 3 , 3 z 1 + 2 z 2 − 3 z 3 , z 1 + 2 z 2 )
关于标准基的矩阵是
( 3 − 1 − 2 3 2 − 3 1 2 0 ) \begin{pmatrix}3&-1&-2\\ 3&2&-3\\ 1&2&0\end{pmatrix} 3 3 1 − 1 2 2 − 2 − 3 0
对角元 3 , 2 , 0 3,2,0 3 , 2 , 0 ,故 tr T = 5 \operatorname{tr}T=5 tr T = 5 。而 T T T 的特征值是 1 , 2 + 3 i , 2 − 3 i 1,\ 2+3\mathrm{i},\ 2-3\mathrm{i} 1 , 2 + 3 i , 2 − 3 i (你可以自己验,重数都是 1 1 1 ),它们的和
1 + ( 2 + 3 i ) + ( 2 − 3 i ) = 5 1+(2+3\mathrm{i})+(2-3\mathrm{i})=5 1 + ( 2 + 3 i ) + ( 2 − 3 i ) = 5
与迹一致。这个例子值得亲手算一遍:你从两条完全不同的路得到同一个数 ,一条只读对角元,一条解三次方程。两条路的交点,就是这一讲的全部内容。
三、几个必须自己算一遍的结论 #
3.1 tr ( A B ) = tr ( B A ) \operatorname{tr}(AB)=\operatorname{tr}(BA) tr ( A B ) = tr ( B A ) #
这是全讲的发动机。注意它的适用范围比想象中宽:A A A 是 m × n m\times n m × n ,B B B 是 n × m n\times m n × m ,不要求两者是方阵,也不要求同形状 。
证明 :A B AB A B 是 m × m m\times m m × m 矩阵,它的第 j j j 个对角元是
( A B ) j , j = ∑ k = 1 n A j , k B k , j (AB)_{j,j}=\sum_{k=1}^{n}A_{j,k}B_{k,j} ( A B ) j , j = k = 1 ∑ n A j , k B k , j
于是
tr ( A B ) = ∑ j = 1 m ∑ k = 1 n A j , k B k , j (34.4) \operatorname{tr}(AB)=\sum_{j=1}^{m}\sum_{k=1}^{n}A_{j,k}B_{k,j} \tag{34.4} tr ( A B ) = j = 1 ∑ m k = 1 ∑ n A j , k B k , j ( 34.4 )
到这里停一下,看清楚要干什么:把两个求和的次序换过来。有限和换序永远合法(这与无穷级数不同,那里换序会出事)。换完:
tr ( A B ) = ∑ k = 1 n ∑ j = 1 m B k , j A j , k = ∑ k = 1 n ( B A ) k , k = tr ( B A ) (34.5) \operatorname{tr}(AB)=\sum_{k=1}^{n}\sum_{j=1}^{m}B_{k,j}A_{j,k}=\sum_{k=1}^{n}(BA)_{k,k}=\operatorname{tr}(BA) \tag{34.5} tr ( A B ) = k = 1 ∑ n j = 1 ∑ m B k , j A j , k = k = 1 ∑ n ( B A ) k , k = tr ( B A ) ( 34.5 )
最后一步:∑ j B k , j A j , k \sum_j B_{k,j}A_{j,k} ∑ j B k , j A j , k 正是 B A BA B A 的第 k k k 个对角元,而 B A BA B A 是 n × n n\times n n × n ,对 k k k 求和就得到它的迹。■ \blacksquare ■
这个证明里只有一件事在起作用: 双重和的指标可以互换。没有别的技巧。所以凡是「求和 + 双线性」的场合,这套论证都能原样搬过去——第 3.9 节的迹内积、第六节量子力学里的期望值,剥到最里层都是它。
① tr ( A B ) = tr ( B A ) \operatorname{tr}(AB)=\operatorname{tr}(BA) tr ( A B ) = tr ( B A ) 不 意味着 A B = B A AB=BA A B = B A 。取 A = ( 0 1 0 0 ) A=\begin{pmatrix}0&1\\ 0&0\end{pmatrix} A = ( 0 0 1 0 ) 、B = ( 0 0 1 0 ) B=\begin{pmatrix}0&0\\ 1&0\end{pmatrix} B = ( 0 1 0 0 ) :A B = ( 1 0 0 0 ) AB=\begin{pmatrix}1&0\\ 0&0\end{pmatrix} A B = ( 1 0 0 0 ) ,B A = ( 0 0 0 1 ) BA=\begin{pmatrix}0&0\\ 0&1\end{pmatrix} B A = ( 0 0 0 1 ) ,两个矩阵不同,但迹都是 1 1 1 。 ② 它不 意味着 tr ( A B ) = tr A ⋅ tr B \operatorname{tr}(AB)=\operatorname{tr}A\cdot\operatorname{tr}B tr ( A B ) = tr A ⋅ tr B (见误区 1)。 ③ 它允许循环移位 ,不允许任意重排。tr ( A B C ) = tr ( C A B ) = tr ( B C A ) \operatorname{tr}(ABC)=\operatorname{tr}(CAB)=\operatorname{tr}(BCA) tr ( A B C ) = tr ( C A B ) = tr ( B C A ) 都成立,但一般 tr ( A B C ) ≠ tr ( A C B ) \operatorname{tr}(ABC)\ne\operatorname{tr}(ACB) tr ( A B C ) = tr ( A C B ) ——把 B B B 和 C C C 直接对调不是循环移位(见第 4.3 节的「环上珠子」)。
3.2 算子的矩阵的迹不依赖于基的选取 #
结论 :设 T ∈ L ( V ) T\in\mathcal{L}(V) T ∈ L ( V ) ,u 1 , … , u n u_1,\dots,u_n u 1 , … , u n 与 v 1 , … , v n v_1,\dots,v_n v 1 , … , v n 是 V V V 的两组基,则
tr M ( T , ( u 1 , … , u n ) ) = tr M ( T , ( v 1 , … , v n ) ) (34.6) \operatorname{tr}\mathcal{M}\big(T,(u_1,\dots,u_n)\big)=\operatorname{tr}\mathcal{M}\big(T,(v_1,\dots,v_n)\big) \tag{34.6} tr M ( T , ( u 1 , … , u n ) ) = tr M ( T , ( v 1 , … , v n ) ) ( 34.6 )
证明 :记 A = M ( T , ( u ∙ ) ) A=\mathcal{M}(T,(u_\bullet)) A = M ( T , ( u ∙ )) 、B = M ( T , ( v ∙ ) ) B=\mathcal{M}(T,(v_\bullet)) B = M ( T , ( v ∙ )) 。10 讲的换基公式说,存在可逆矩阵 C C C 使
A = C − 1 B C A=C^{-1}BC A = C − 1 B C
(C C C 就是从新基到旧基的过渡矩阵。)于是
tr A = tr ( C − 1 ( B C ) ) = tr ( ( B C ) C − 1 ) = tr ( B ( C C − 1 ) ) = tr B (34.7) \begin{aligned}\operatorname{tr}A&=\operatorname{tr}\big(C^{-1}(BC)\big)\\ &=\operatorname{tr}\big((BC)\,C^{-1}\big)\\ &=\operatorname{tr}\big(B\,(CC^{-1})\big)\\ &=\operatorname{tr}B\end{aligned} \tag{34.7} tr A = tr ( C − 1 ( B C ) ) = tr ( ( B C ) C − 1 ) = tr ( B ( C C − 1 ) ) = tr B ( 34.7 )
第二步就是把 C − 1 C^{-1} C − 1 从最前面循环移到后面去,用的正是 3.1。第三步用结合律,第四步用 C C − 1 = I CC^{-1}=I C C − 1 = I 与 tr ( B I ) = tr B \operatorname{tr}(BI)=\operatorname{tr}B tr ( B I ) = tr B 。■ \blacksquare ■
同一个套路可以证明「tr ( C − 1 B C ) \operatorname{tr}(C^{-1}BC) tr ( C − 1 B C ) 与 C C C 无关」这一类命题:把 C − 1 C^{-1} C − 1 一路往右推,推到它和 C C C 碰头,消掉。 行列式那边是另一个套路(乘性:det ( C − 1 B C ) = det C − 1 det B det C \det(C^{-1}BC)=\det C^{-1}\det B\det C det ( C − 1 B C ) = det C − 1 det B det C ),两者不要混。
于是定义站得住:
tr T : = tr M ( T , ( v 1 , … , v n ) ) , v 1 , … , v n 为 V 的任意一组基 (34.8) \operatorname{tr}T:=\operatorname{tr}\mathcal{M}\big(T,(v_1,\dots,v_n)\big),\qquad v_1,\dots,v_n\ \text{为}\ V\ \text{的任意一组基} \tag{34.8} tr T := tr M ( T , ( v 1 , … , v n ) ) , v 1 , … , v n 为 V 的任意一组基 ( 34.8 )
3.3 手算一次:对角元会变,和不变 #
空口说「不变」没有说服力,动手换一组基看看。
取 B = ( 1 1 0 2 ) B=\begin{pmatrix}1&1\\ 0&2\end{pmatrix} B = ( 1 0 1 2 ) (对角元 1 , 2 1,2 1 , 2 ,tr B = 3 \operatorname{tr}B=3 tr B = 3 ),再取 C = ( 1 2 3 5 ) C=\begin{pmatrix}1&2\\ 3&5\end{pmatrix} C = ( 1 3 2 5 ) ,det C = − 1 \det C=-1 det C = − 1 ,C − 1 = ( − 5 2 3 − 1 ) C^{-1}=\begin{pmatrix}-5&2\\ 3&-1\end{pmatrix} C − 1 = ( − 5 3 2 − 1 ) 。算 C − 1 B C C^{-1}BC C − 1 B C :
B C = ( 1 1 0 2 ) ( 1 2 3 5 ) = ( 4 7 6 10 ) BC=\begin{pmatrix}1&1\\ 0&2\end{pmatrix}\begin{pmatrix}1&2\\ 3&5\end{pmatrix}=\begin{pmatrix}4&7\\ 6&10\end{pmatrix} B C = ( 1 0 1 2 ) ( 1 3 2 5 ) = ( 4 6 7 10 )
C − 1 ( B C ) = ( − 5 2 3 − 1 ) ( 4 7 6 10 ) = ( − 8 − 15 6 11 ) C^{-1}(BC)=\begin{pmatrix}-5&2\\ 3&-1\end{pmatrix}\begin{pmatrix}4&7\\ 6&10\end{pmatrix}=\begin{pmatrix}-8&-15\\ 6&11\end{pmatrix} C − 1 ( B C ) = ( − 5 3 2 − 1 ) ( 4 6 7 10 ) = ( − 8 6 − 15 11 )
新矩阵的对角元是 − 8 -8 − 8 和 11 11 11 ——离原来的 1 , 2 1,2 1 , 2 差得远,甚至符号都变了。但
− 8 + 11 = 3 = 1 + 2 -8+11=3=1+2 − 8 + 11 = 3 = 1 + 2
对角元是坐标的产物,和是算子的属性。 这句话可以当作整讲的注脚。顺便,− 8 -8 − 8 和 11 11 11 都不是这个算子的特征值(特征值是 1 , 2 1,2 1 , 2 )——所以「对角元是特征值」这条只对上三角矩阵成立,绝不能推广到一般矩阵。
3.4 迹等于特征值之和 #
结论 :设 F = C \mathbb{F}=\mathbb{C} F = C ,T ∈ L ( V ) T\in\mathcal{L}(V) T ∈ L ( V ) ,dim V = n \dim V=n dim V = n 。则
tr T = λ 1 + ⋯ + λ n \operatorname{tr}T=\lambda_1+\cdots+\lambda_n tr T = λ 1 + ⋯ + λ n
其中各特征值出现的次数等于它的(代数)重数。
证明 :按 19 讲,复空间上存在一组基使 T T T 关于它有上三角矩阵,且对角元恰是 T T T 的全部特征值(按重数计)。于是
tr T = tr M ( T , 该基 ) = 对角元之和 = 特征值之和(含重数) \operatorname{tr}T=\operatorname{tr}\mathcal{M}(T,\text{该基})=\text{对角元之和}=\text{特征值之和(含重数)} tr T = tr M ( T , 该基 ) = 对角元之和 = 特征值之和(含重数)
第一步用的是 3.2 给的自由:想看迹,就挑最方便的基 。挑上三角那组,对角元一目了然。■ \blacksquare ■
这个证明短得惊人,但它的关键不在计算,在 3.2:正因为迹与基无关,「挑一组好基来看」才是合法操作。如果没有 3.2,上三角那组基算出来的只是「那个矩阵的迹」,跟 T T T 没关系。
在 R \mathbb{R} R 上,旋转 R θ R_\theta R θ (θ \theta θ 不是 π \pi π 的倍数)一个实特征值都没有 ,谈不上「特征值之和」。 但它的迹 2 cos θ 2\cos\theta 2 cos θ 有定义。所以等式 (34.1) 必须在「特征值齐全」的场合才有意义—— 要么 F = C \mathbb{F}=\mathbb{C} F = C ,要么特征多项式在 F \mathbb{F} F 上完全分解。 这时也可以说:迹是特征值之和在实数域上的「残影」。
3.5 迹与特征多项式 #
把特征多项式按「特征值」写出来(复空间,λ 1 , … , λ n \lambda_1,\dots,\lambda_n λ 1 , … , λ n 按重数计):
p ( z ) = ( z − λ 1 ) ( z − λ 2 ) ⋯ ( z − λ n ) p(z)=(z-\lambda_1)(z-\lambda_2)\cdots(z-\lambda_n) p ( z ) = ( z − λ 1 ) ( z − λ 2 ) ⋯ ( z − λ n )
展开,看最高两次的系数:
p ( z ) = z n − ( λ 1 + ⋯ + λ n ) z n − 1 + ⋯ + ( − 1 ) n λ 1 λ 2 ⋯ λ n (34.9) p(z)=z^n-(\lambda_1+\cdots+\lambda_n)z^{n-1}+\cdots+(-1)^n\lambda_1\lambda_2\cdots\lambda_n \tag{34.9} p ( z ) = z n − ( λ 1 + ⋯ + λ n ) z n − 1 + ⋯ + ( − 1 ) n λ 1 λ 2 ⋯ λ n ( 34.9 )
于是直接读出两条:
tr T = − [ z n − 1 的系数 ] , det T = ( − 1 ) n [ 常数项 ] (34.10) \boxed{\ \text{tr}\,T=-[z^{n-1}\text{ 的系数}],\qquad \det T=(-1)^n[\text{常数项}]\ } \tag{34.10} tr T = − [ z n − 1 的系数 ] , det T = ( − 1 ) n [ 常数项 ] ( 34.10 )
这也是为什么第 3.3 节那个 B = ( 1 1 0 2 ) B=\begin{pmatrix}1&1\\ 0&2\end{pmatrix} B = ( 1 0 1 2 ) 的特征多项式 z 2 − 3 z + 2 z^2-3z+2 z 2 − 3 z + 2 里,中项系数 − 3 -3 − 3 正好是 − tr B -\operatorname{tr}B − tr B 。
一个记忆框架:迹和行列式是初等对称多项式的两个端点 n n n 个数的初等对称多项式有 n n n 个:e 1 = ∑ λ i e_1=\sum\lambda_i e 1 = ∑ λ i 、e 2 = ∑ i < j λ i λ j e_2=\sum_{i<j}\lambda_i\lambda_j e 2 = ∑ i < j λ i λ j 、…、e n = ∏ λ i e_n=\prod\lambda_i e n = ∏ λ i 。 特征多项式把它们全部 打包了:p ( z ) = z n − e 1 z n − 1 + e 2 z n − 2 − ⋯ + ( − 1 ) n e n p(z)=z^n-e_1z^{n-1}+e_2z^{n-2}-\cdots+(-1)^ne_n p ( z ) = z n − e 1 z n − 1 + e 2 z n − 2 − ⋯ + ( − 1 ) n e n 。 迹是 e 1 e_1 e 1 ,行列式是 e n e_n e n 。中间那些 e 2 , … , e n − 1 e_2,\dots,e_{n-1} e 2 , … , e n − 1 也各有矩阵表达式(所有 k k k 阶主子式之和),只是没有专门的记号与名字。所以「迹 + 行列式」是你能免费读到的两个端点,中间那些要另外花力气。
3.6 迹是线性泛函 #
结论 :tr : L ( V ) → F \operatorname{tr}:\mathcal{L}(V)\to\mathbb{F} tr : L ( V ) → F 是线性的,即
tr ( S + T ) = tr S + tr T , tr ( λ T ) = λ tr T (34.11) \operatorname{tr}(S+T)=\operatorname{tr}S+\operatorname{tr}T,\qquad \operatorname{tr}(\lambda T)=\lambda\operatorname{tr}T \tag{34.11} tr ( S + T ) = tr S + tr T , tr ( λ T ) = λ tr T ( 34.11 )
并且
tr ( S T ) = tr ( T S ) , tr I = dim V (34.12) \operatorname{tr}(ST)=\operatorname{tr}(TS),\qquad \operatorname{tr}I=\dim V \tag{34.12} tr ( S T ) = tr ( T S ) , tr I = dim V ( 34.12 )
证明 :固定一组基,把算子换成矩阵。矩阵的迹显然满足前两条(对角元逐个相加、逐个乘 λ \lambda λ ,而矩阵的加法与标量乘是逐元素的)。第三条就是 3.1 搬到算子层面(用 M ( S T ) = M ( S ) M ( T ) \mathcal{M}(ST)=\mathcal{M}(S)\mathcal{M}(T) M ( S T ) = M ( S ) M ( T ) )。第四条:I I I 的矩阵是单位矩阵,对角元 n n n 个 1 1 1 。■ \blacksquare ■
「tr \operatorname{tr} tr 是 L ( V ) \mathcal{L}(V) L ( V ) 上的线性泛函」这句话在 12 讲的对偶空间语言里是什么意思?L ( V ) \mathcal{L}(V) L ( V ) 是 n 2 n^2 n 2 维向量空间,线性泛函是它的对偶空间里的一个元素。迹是其中一个具体、可算、且与基无关 的元素。
任何线性泛函 τ : L ( V ) → F \tau:\mathcal{L}(V)\to\mathbb{F} τ : L ( V ) → F ,只要满足 τ ( S T ) = τ ( T S ) \tau(ST)=\tau(TS) τ ( S T ) = τ ( T S ) 与 τ ( I ) = dim V \tau(I)=\dim V τ ( I ) = dim V ,就必有 τ = tr \tau=\operatorname{tr} τ = tr 。 证明思路(不必背,但要能复述三步):取基 v 1 , … , v n v_1,\dots,v_n v 1 , … , v n ,令 P j , k P_{j,k} P j , k 是「把第 k k k 个坐标搬到第 j j j 位」的算子,即 P j , k ( a 1 v 1 + ⋯ + a n v n ) = a k v j P_{j,k}(a_1v_1+\cdots+a_nv_n)=a_kv_j P j , k ( a 1 v 1 + ⋯ + a n v n ) = a k v j 。则 ① P j , k P k , k = P j , k P_{j,k}P_{k,k}=P_{j,k} P j , k P k , k = P j , k ,故 τ ( P j , k ) = τ ( P k , k P j , k ) \tau(P_{j,k})=\tau(P_{k,k}P_{j,k}) τ ( P j , k ) = τ ( P k , k P j , k ) ;而当 j ≠ k j\ne k j = k 时 P k , k P j , k = 0 P_{k,k}P_{j,k}=0 P k , k P j , k = 0 ,于是 τ ( P j , k ) = 0 \tau(P_{j,k})=0 τ ( P j , k ) = 0 。 ② P j , k P k , j = P j , j P_{j,k}P_{k,j}=P_{j,j} P j , k P k , j = P j , j 与 P k , j P j , k = P k , k P_{k,j}P_{j,k}=P_{k,k} P k , j P j , k = P k , k 合起来(配合①)给出 τ ( P j , j ) = τ ( P k , k ) \tau(P_{j,j})=\tau(P_{k,k}) τ ( P j , j ) = τ ( P k , k ) :所有对角元的 τ \tau τ 值相同,记作 c c c 。 ③ τ ( I ) = ∑ j τ ( P j , j ) = n c \tau(I)=\sum_j\tau(P_{j,j})=n\,c τ ( I ) = ∑ j τ ( P j , j ) = n c ,而 τ ( I ) = dim V = n \tau(I)=\dim V=n τ ( I ) = dim V = n ,故 c = 1 c=1 c = 1 。 于是任一算子 T = ∑ j , k M ( T ) j , k P j , k T=\sum_{j,k}\mathcal{M}(T)_{j,k}P_{j,k} T = ∑ j , k M ( T ) j , k P j , k 满足 τ ( T ) = ∑ j M ( T ) j , j = tr T \tau(T)=\sum_j\mathcal{M}(T)_{j,j}=\operatorname{tr}T τ ( T ) = ∑ j M ( T ) j , j = tr T 。意义 :这解释了为什么 tr ( S T ) = tr ( T S ) \operatorname{tr}(ST)=\operatorname{tr}(TS) tr ( S T ) = tr ( T S ) 不是一条孤立的技术性引理,而是刻画 了迹——它是「乘法顺序无关的线性读数」这件事的唯一解。
3.7 用迹证「S T − T S = I ST-TS=I S T − T S = I 不存在」 #
结论 :设 dim V < ∞ \dim V<\infty dim V < ∞ ,则不存在 S , T ∈ L ( V ) S,T\in\mathcal{L}(V) S , T ∈ L ( V ) 使 S T − T S = I ST-TS=I S T − T S = I 。
证明 :
tr ( S T − T S ) = tr ( S T ) − tr ( T S ) = 0 \operatorname{tr}(ST-TS)=\operatorname{tr}(ST)-\operatorname{tr}(TS)=0 tr ( S T − T S ) = tr ( S T ) − tr ( T S ) = 0
而 tr I = dim V ≠ 0 \operatorname{tr}I=\dim V\ne 0 tr I = dim V = 0 。两个算子迹不同,故不相等。■ \blacksquare ■
它是「不变量论证 」的样板:你不需要研究 S T − T S ST-TS S T − T S 长什么样,只需要找一个「两边必然不同」的读数。 迹恰好就是这样一个读数,而且它便宜到只数对角元。 对比一下:你可以用行列式,但 det ( S T − T S ) \det(ST-TS) det ( S T − T S ) 与 det I = 1 \det I=1 det I = 1 之间推不出矛盾(det \det det 不可加,S T − T S ST-TS S T − T S 的长相也读不出来)。能承担这个论证的,必须是可加的不变量——这就把范围缩到了迹附近。
顺便说:这条在无限维不成立 。在全体多项式构成的空间上,取 T T T 为「乘 x x x 」、S S S 为求导,则 S T − T S = I ST-TS=I S T − T S = I 。验一下:作用在 x k x^k x k 上,S T ST S T 是先乘再导,得 ( k + 1 ) x k (k+1)x^k ( k + 1 ) x k ;T S TS T S 是先导再乘,得 k x k kx^k k x k ;两者相减恰是 x k x^k x k 。迹在那里根本没有定义——无限维的那个和不收敛。定理的边界正好落在「迹有定义」这条线上 ,这个观察比定理本身更有意思。
3.8 内积空间上的迹:把对角元看成投影 #
结论 :设 V V V 是内积空间,e 1 , … , e n e_1,\dots,e_n e 1 , … , e n 是一组规范正交基 。则
tr T = ⟨ T e 1 , e 1 ⟩ + ⟨ T e 2 , e 2 ⟩ + ⋯ + ⟨ T e n , e n ⟩ (34.13) \operatorname{tr}T=\langle Te_1,e_1\rangle+\langle Te_2,e_2\rangle+\cdots+\langle Te_n,e_n\rangle \tag{34.13} tr T = ⟨ T e 1 , e 1 ⟩ + ⟨ T e 2 , e 2 ⟩ + ⋯ + ⟨ T e n , e n ⟩ ( 34.13 )
证明 :矩阵的第 k k k 个对角元是 T T T 在 e k e_k e k 上的坐标的第 k k k 个分量。而 23 讲给过:在规范正交基下,向量的第 k k k 个坐标就是它和 e k e_k e k 的内积。取向量 v = T e k v=Te_k v = T e k ,得 M ( T ) k , k = ⟨ T e k , e k ⟩ \mathcal{M}(T)_{k,k}=\langle Te_k,e_k\rangle M ( T ) k , k = ⟨ T e k , e k ⟩ ,对 k k k 求和即得。■ \blacksquare ■
这个公式给了一个非常好的直觉,第 4.2 节展开。这里先记账。
一般基下 M ( T ) k , k ≠ ⟨ T e k , e k ⟩ \mathcal{M}(T)_{k,k}\ne\langle Te_k,e_k\rangle M ( T ) k , k = ⟨ T e k , e k ⟩ 。 例:R 2 \mathbb{R}^2 R 2 取 f 1 = ( 1 , 1 ) f_1=(1,1) f 1 = ( 1 , 1 ) (未单位化),T = I T=I T = I 。则 ⟨ T f 1 , f 1 ⟩ = 2 \langle Tf_1,f_1\rangle=2 ⟨ T f 1 , f 1 ⟩ = 2 ,但 I I I 关于 { f 1 , … } \{f_1,\dots\} { f 1 , … } 的矩阵对角元是 1 1 1 。 差的正是 ∥ f 1 ∥ 2 = 2 \|f_1\|^2=2 ∥ f 1 ∥ 2 = 2 这个因子。内积公式只对单位正交的基成立。
3.9 迹内积:tr ( T ∗ S ) \operatorname{tr}(T^*S) tr ( T ∗ S ) 是 L ( V , W ) \mathcal{L}(V,W) L ( V , W ) 上的内积 #
把 3.8 与 3.1 合起来,会冒出一个很好用的东西:对有限维内积空间 V , W V,W V , W ,定义
⟨ S , T ⟩ : = tr ( T ∗ S ) (34.14) \langle S,T\rangle:=\operatorname{tr}(T^*S) \tag{34.14} ⟨ S , T ⟩ := tr ( T ∗ S ) ( 34.14 )
则这确实 是 L ( V , W ) \mathcal{L}(V,W) L ( V , W ) 上的一个内积。要验三件事:
共轭对称 :tr ( T ∗ S ) ‾ = tr ( S ∗ T ) \overline{\operatorname{tr}(T^*S)}=\operatorname{tr}(S^*T) tr ( T ∗ S ) = tr ( S ∗ T ) ,这靠 tr ( M ∗ ) = tr M ‾ \operatorname{tr}(M^*)=\overline{\operatorname{tr}M} tr ( M ∗ ) = tr M (见下方)。
对第一个变量线性 :T ∗ T^* T ∗ 与 tr \operatorname{tr} tr 各自线性。
正定性 :tr ( T ∗ T ) ≥ 0 \operatorname{tr}(T^*T)\ge 0 tr ( T ∗ T ) ≥ 0 ,且等号只在 T = 0 T=0 T = 0 时成立。这条要算一下:取 V V V 的规范正交基 e 1 , … , e n e_1,\dots,e_n e 1 , … , e n ,由 3.8(对算子 T ∗ T T^*T T ∗ T 用它),
tr ( T ∗ T ) = ∑ k = 1 n ⟨ T ∗ T e k , e k ⟩ = ∑ k = 1 n ⟨ T e k , T e k ⟩ = ∑ k = 1 n ∥ T e k ∥ 2 (34.15) \operatorname{tr}(T^*T)=\sum_{k=1}^{n}\langle T^*Te_k,e_k\rangle=\sum_{k=1}^{n}\langle Te_k,Te_k\rangle=\sum_{k=1}^{n}\|Te_k\|^2 \tag{34.15} tr ( T ∗ T ) = k = 1 ∑ n ⟨ T ∗ T e k , e k ⟩ = k = 1 ∑ n ⟨ T e k , T e k ⟩ = k = 1 ∑ n ∥ T e k ∥ 2 ( 34.15 )
这是若干非负数之和,为零当且仅当每个 T e k = 0 Te_k=0 T e k = 0 ,即 T = 0 T=0 T = 0 。
它长什么样 :若再把每个 T e k Te_k T e k 在 f 1 , … , f m f_1,\dots,f_m f 1 , … , f m 下展开,则
tr ( T ∗ T ) = ∑ k = 1 n ∑ j = 1 m ∣ ⟨ T e k , f j ⟩ ∣ 2 (34.16) \operatorname{tr}(T^*T)=\sum_{k=1}^{n}\sum_{j=1}^{m}\big|\langle Te_k,f_j\rangle\big|^2 \tag{34.16} tr ( T ∗ T ) = k = 1 ∑ n j = 1 ∑ m ⟨ T e k , f j ⟩ 2 ( 34.16 )
而 ⟨ T e k , f j ⟩ \langle Te_k,f_j\rangle ⟨ T e k , f j ⟩ 正是 T T T 的矩阵的第 j j j 行 k k k 列元素。所以:
tr ( T ∗ T ) \operatorname{tr}(T^*T) tr ( T ∗ T ) = 矩阵全部元素绝对值的平方和。
这就是线性代数里那个最常见的范数 ∥ A ∥ F = ∑ j , k ∣ a j , k ∣ 2 \|A\|_F=\sqrt{\sum_{j,k}|a_{j,k}|^2} ∥ A ∥ F = ∑ j , k ∣ a j , k ∣ 2 的坐标无关解释。它有个名字:Frobenius 范数 (也叫 Hilbert–Schmidt 范数)。
同一个矩阵,T ∗ T T^*T T ∗ T 的迹给出「元素平方和」,而算子范数(22 讲)给出「最坏方向上的拉伸倍数」。两者不同,别混。 关系是 ∥ T ∥ ≤ ∥ T ∥ F ≤ n ∥ T ∥ \|T\|\le\|T\|_F\le\sqrt{n}\,\|T\| ∥ T ∥ ≤ ∥ T ∥ F ≤ n ∥ T ∥ (n = dim V n=\dim V n = dim V )。Frobenius 范数好算(数元素即可),算子范数难算(要优化)。 数值线性代数里用前者做估计、用后者做保证,就是这个原因。
同样地,tr ( T ∗ ) = tr T ‾ \operatorname{tr}(T^*)=\overline{\operatorname{tr}T} tr ( T ∗ ) = tr T :转置共轭把矩阵的元素换成共轭再转置,对角元变成 a j , j ‾ \overline{a_{j,j}} a j , j ,和不就是 tr A ‾ \overline{\operatorname{tr}A} tr A 。这个细节在做自伴算子(25 讲)的迹时会用到:T = T ∗ T=T^* T = T ∗ 时 tr T \operatorname{tr}T tr T 必为实数。
3.10 两个「迹为零却有意思」的对照 #
由 3.4,tr T = ∑ λ i \operatorname{tr}T=\sum\lambda_i tr T = ∑ λ i ,于是:
T T T 幂零(某次方为零)⇒ \Rightarrow ⇒ tr T = 0 \operatorname{tr}T=0 tr T = 0 。 因为幂零算子的特征值全是 0 0 0 (19 讲:上三角后对角元全零)。反过来不成立:( 0 1 0 0 ) \begin{pmatrix}0&1\\ 0&0\end{pmatrix} ( 0 0 1 0 ) 幂零,但 ( 1 0 0 − 1 ) \begin{pmatrix}1&0\\ 0&-1\end{pmatrix} ( 1 0 0 − 1 ) 迹为 0 0 0 而特征值 ± 1 \pm1 ± 1 ,不是幂零。
T T T 正算子(27 讲)且 tr T = 0 \operatorname{tr}T=0 tr T = 0 ⇒ \Rightarrow ⇒ T = 0 T=0 T = 0 。 因为正算子的特征值都 ≥ 0 \ge 0 ≥ 0 (谱定理),而迹是它们之和;若干个非负数加起来为 0 0 0 ,只能全为 0 0 0 ,于是 T T T 的所有特征值为零,正算子因而 T = 0 T=0 T = 0 。
这两条并排看特别有教益:同样是「迹为零」,一边只说得出「特征值正负抵消」,另一边能说「就是零」。 差别全在「特征值的符号是否受控」。这也是为什么量子力学里 tr ( ρ A ) \operatorname{tr}(\rho A) tr ( ρ A ) 能当成一个期望值:ρ \rho ρ 正、tr ρ = 1 \operatorname{tr}\rho=1 tr ρ = 1 ,这两个条件把迹钉成了一个平均值。
四、几何直观 / 记忆钩子 #
4.1 钩子一:迹是「平均缩放率」乘维数 #
tr T = ∑ λ i \operatorname{tr}T=\sum\lambda_i tr T = ∑ λ i ,所以
tr T n = 特征值的算术平均 \frac{\operatorname{tr}T}{n}=\text{特征值的算术平均} n tr T = 特征值的算术平均
对旋转 R θ R_\theta R θ :平均是 cos θ \cos\theta cos θ ,迹是 2 cos θ 2\cos\theta 2 cos θ 。θ = 0 \theta=0 θ = 0 时迹 = 2 =2 = 2 (不动),θ = π \theta=\pi θ = π 时迹 = − 2 =-2 = − 2 (翻面),θ = π / 2 \theta=\pi/2 θ = π /2 时迹 = 0 =0 = 0 (转 90 ∘ 90^\circ 9 0 ∘ ,把每个方向都转走)。迹的符号在告诉你「总体上是在放大、缩小,还是转得厉害」。
4.2 钩子二:迹是「各方向上自己投影的总和」 #
由 3.8,tr T = ∑ ⟨ T e k , e k ⟩ \operatorname{tr}T=\sum\langle Te_k,e_k\rangle tr T = ∑ ⟨ T e k , e k ⟩ 。每一项 ⟨ T e k , e k ⟩ \langle Te_k,e_k\rangle ⟨ T e k , e k ⟩ 的意思是:
把 e k e_k e k 用 T T T 推一下,再看它在原方向 上留了多少。
若 e k e_k e k 是特征向量(T e k = λ e k Te_k=\lambda e_k T e k = λ e k ),这一项就是 λ \lambda λ ;
若 T T T 把 e k e_k e k 转走了(垂直),这一项就是 0 0 0 。
所以迹只看得见「不把方向转跑的那部分」 。这也解释了为什么幂零矩阵迹为零:它把每个方向都往更深一层推,没有留下原方向的影子。
4.3 钩子三:只有循环,没有交换 #
tr ( A B C ) = tr ( C A B ) = tr ( B C A ) \operatorname{tr}(ABC)=\operatorname{tr}(CAB)=\operatorname{tr}(BCA) tr ( A B C ) = tr ( C A B ) = tr ( B C A )
记忆法:把 A , B , C A,B,C A , B , C 想象成串在一条环上的珠子,只能整体转动环,不能拔下来重穿。 一旦你把 B B B 和 C C C 对调(不经过旋转),等式一般就断。
反例 (自己验一遍,两行乘法):取
A = ( 1 0 0 0 ) , B = ( 0 1 0 0 ) , C = ( 0 0 1 0 ) A=\begin{pmatrix}1&0\\ 0&0\end{pmatrix},\quad B=\begin{pmatrix}0&1\\ 0&0\end{pmatrix},\quad C=\begin{pmatrix}0&0\\ 1&0\end{pmatrix} A = ( 1 0 0 0 ) , B = ( 0 0 1 0 ) , C = ( 0 1 0 0 )
则 B C = ( 1 0 0 0 ) BC=\begin{pmatrix}1&0\\ 0&0\end{pmatrix} B C = ( 1 0 0 0 ) ,A ( B C ) = ( 1 0 0 0 ) A(BC)=\begin{pmatrix}1&0\\ 0&0\end{pmatrix} A ( B C ) = ( 1 0 0 0 ) ,迹 = 1 =1 = 1 ; 而 A C = 0 AC=0 A C = 0 (A A A 只留第一行,C C C 只留第二列,第一行第二列是 0 0 0 ),故 A C B = 0 ACB=0 A C B = 0 ,迹 = 0 =0 = 0 。1 ≠ 0 1\ne0 1 = 0 ,循环性到此为止。
因为全部工具只有一条 :tr ( X Y ) = tr ( Y X ) \operatorname{tr}(XY)=\operatorname{tr}(YX) tr ( X Y ) = tr ( Y X ) ,它对任意的「切开点」都成立。tr ( A B C ) → tr ( B C A ) \operatorname{tr}(ABC)\to\operatorname{tr}(BCA) tr ( A B C ) → tr ( B C A ) 就是取 X = A X=A X = A 、Y = B C Y=BC Y = B C ;tr ( A B C ) → tr ( C A B ) \operatorname{tr}(ABC)\to\operatorname{tr}(CAB) tr ( A B C ) → tr ( C A B ) 就是取 X = A B X=AB X = A B 、Y = C Y=C Y = C 。 这两种切法,恰好给出三次循环。而「把 B , C B,C B , C 对调」这个排列没有对应的切法—— 它需要同时把 B B B 挪到 C C C 后面、把 C C C 挪到 B B B 前面,方向相反,一次 X Y → Y X XY\to YX X Y → Y X 办不到。结论 :能被单次或多次「切一刀再交换」造出来的重排,只有循环那一族。
4.4 一张对照表:迹 vs 行列式 #
项
迹 tr \operatorname{tr} tr
行列式 det \det det
特征值层面
和 ∑ λ i \sum\lambda_i ∑ λ i
积 ∏ λ i \prod\lambda_i ∏ λ i
特征多项式
z n − 1 z^{n-1} z n − 1 的系数取负
常数项乘 ( − 1 ) n (-1)^n ( − 1 ) n
基无关
是(靠循环性)
是(靠乘性 det ( C − 1 B C ) = det B \det(C^{-1}BC)=\det B det ( C − 1 B C ) = det B )
循环性
tr ( A B C ) = tr ( C A B ) \operatorname{tr}(ABC)=\operatorname{tr}(CAB) tr ( A B C ) = tr ( C A B )
det ( A B C ) = det ( C A B ) \det(ABC)=\det(CAB) det ( A B C ) = det ( C A B ) (因为都可交换地乘)
可逆判据
不 能判(tr = 0 \operatorname{tr}=0 tr = 0 也可逆)
能判(det = 0 ⟺ \det=0\iff det = 0 ⟺ 不可逆)
计算成本
数 n n n 个元素
需要消元,O ( n 3 ) O(n^3) O ( n 3 )
几何含义
各方向缩放率的和
体积缩放倍数
迹是加法性的,行列式是乘法性的。 所以「可逆」这种与乘有关的问题归行列式;「求和、线性、期望」这种与加有关的问题归迹。
4.5 钩子四:迹是「最便宜的基无关数」 #
需要校验时先算迹,需要判可逆时先行列式,需要精确的谱信息时才老实求特征多项式。这是「成本从低到高」的三级 ,考试和工程里都按这个顺序用。
五、常见误区 #
误区 1:以为 tr ( A B ) = tr A ⋅ tr B \operatorname{tr}(AB)=\operatorname{tr}A\cdot\operatorname{tr}B tr ( A B ) = tr A ⋅ tr B 反例:A = B = ( 1 0 0 1 ) A=B=\begin{pmatrix}1&0\\ 0&1\end{pmatrix} A = B = ( 1 0 0 1 ) ,左边 = 2 =2 = 2 ,右边 = 4 =4 = 4 。tr \operatorname{tr} tr 是 线性泛函(保加法、保标量乘),但不是 乘法同态。行列式才是乘法同态,别把两个的性质对调。 这条错误的根源是:把「线性」误读成「和乘法相容」。线性只关加法。
误区 2:以为 tr ( A 2 ) = ( tr A ) 2 \operatorname{tr}(A^2)=(\operatorname{tr}A)^2 tr ( A 2 ) = ( tr A ) 2 同上,取 A = ( 1 0 0 2 ) A=\begin{pmatrix}1&0\\ 0&2\end{pmatrix} A = ( 1 0 0 2 ) :A 2 = ( 1 0 0 4 ) A^2=\begin{pmatrix}1&0\\ 0&4\end{pmatrix} A 2 = ( 1 0 0 4 ) ,迹 = 5 =5 = 5 ;而 ( tr A ) 2 = 9 (\operatorname{tr}A)^2=9 ( tr A ) 2 = 9 。 正确的算法是 tr ( A 2 ) = λ 1 2 + λ 2 2 \operatorname{tr}(A^2)=\lambda_1^2+\lambda_2^2 tr ( A 2 ) = λ 1 2 + λ 2 2 (由 3.4,因为 A 2 A^2 A 2 的特征值是 λ i 2 \lambda_i^2 λ i 2 )——平方进到特征值里去,不是进到和里去 。
误区 3:以为「迹与基无关」可以推广到别的对角元函数 迹无关,det \det det 也无关,但「对角元的平方和」「对角线乘积」「第一大对角元」全都有赖于基。 更细的一层:∑ j , k ∣ a j , k ∣ 2 \sum_{j,k}|a_{j,k}|^2 ∑ j , k ∣ a j , k ∣ 2 (全部元素平方和)在规范正交 换基下不变(3.9),在一般换基下变 。 判据:问「这个量是不是某个循环不变的表达式的迹」,是则不变,否则别指望。
误区 4:以为只有同阶方阵才能用 tr ( A B ) = tr ( B A ) \operatorname{tr}(AB)=\operatorname{tr}(BA) tr ( A B ) = tr ( B A ) 恰恰相反,这条最有用的时候正是 A , B A,B A , B 不同形状 。取 A A A 为 2 × 3 2\times3 2 × 3 、B B B 为 3 × 2 3\times2 3 × 2 :
A = ( 1 2 0 0 1 1 ) , B = ( 1 0 0 1 2 3 ) A=\begin{pmatrix}1&2&0\\ 0&1&1\end{pmatrix},\qquad B=\begin{pmatrix}1&0\\ 0&1\\ 2&3\end{pmatrix} A = ( 1 0 2 1 0 1 ) , B = 1 0 2 0 1 3
A B AB A B 是 2 × 2 2\times2 2 × 2 :A B = ( 1 2 2 4 ) AB=\begin{pmatrix}1&2\\ 2&4\end{pmatrix} A B = ( 1 2 2 4 ) ,tr ( A B ) = 1 + 4 = 5 \operatorname{tr}(AB)=1+4=5 tr ( A B ) = 1 + 4 = 5 。B A BA B A 是 3 × 3 3\times3 3 × 3 :B A = ( 1 2 0 0 1 1 2 7 3 ) BA=\begin{pmatrix}1&2&0\\ 0&1&1\\ 2&7&3\end{pmatrix} B A = 1 0 2 2 1 7 0 1 3 ,tr ( B A ) = 1 + 1 + 3 = 5 \operatorname{tr}(BA)=1+1+3=5 tr ( B A ) = 1 + 1 + 3 = 5 。
两者阶数都不同 ,当然不是同一个矩阵;但迹一样。这正说明迹是一条比「矩阵相等」弱得多、因而适用面宽得多的信息。 顺便留意 A B AB A B 的对角元 1 , 4 1,4 1 , 4 与 B A BA B A 的 1 , 1 , 3 1,1,3 1 , 1 , 3 ——连对角元本身都对不上,只有和相等 。这一点别指望更多。
旋转 R 90 ∘ = ( 0 − 1 1 0 ) R_{90^\circ}=\begin{pmatrix}0&-1\\ 1&0\end{pmatrix} R 9 0 ∘ = ( 0 1 − 1 0 ) 在 R \mathbb{R} R 上没有实特征值,在 C \mathbb{C} C 上才有 ± i \pm\mathrm{i} ± i ,和为零,与迹 = 0 {}=0 = 0 一致。 但如果你在 R \mathbb{R} R 上硬说「实特征值之和」,那句话连主语都没有。 正确表述:F = C \mathbb{F}=\mathbb{C} F = C 时成立;一般情形下,需假设特征多项式在 F \mathbb{F} F 上完全分解。 (19 讲那条等价条件。)
幂零矩阵 ( 0 1 0 0 ) \begin{pmatrix}0&1\\ 0&0\end{pmatrix} ( 0 0 1 0 ) 迹为零但显然非零。 只有加上「正算子」这个前提(3.10),迹为零才能推出算子为零。永远先问:我这两个特征值有没有符号约束?
六、拓展(可选,可跳过) #
6.1 迹在量子力学里是「期望值」 #
量子态是一个正算子 ρ \rho ρ (密度矩阵),满足 tr ρ = 1 \operatorname{tr}\rho=1 tr ρ = 1 ——因为「概率之和为一」。可观测量 A A A 的期望值写成
⟨ A ⟩ = tr ( ρ A ) \langle A\rangle=\operatorname{tr}(\rho A) ⟨ A ⟩ = tr ( ρ A )
为什么必须是迹?因为 tr ( ρ A ) = tr ( A ρ ) \operatorname{tr}(\rho A)=\operatorname{tr}(A\rho) tr ( ρ A ) = tr ( A ρ ) ,所以「先算态再测」与「先测算子再配态」给出同一个数——这个次序无关性 正是期望值该有的性质。而 3.1 是它唯一的来源。
6.2 无限维的那条边界 #
3.7 在有限维成立,无限维失败(求导与乘 x x x 的例子)。这不是巧合:只要给 S , T S,T S , T 加上「迹类」之类的紧性限制,结论在无限维仍能救回来,而那套语言正是量子场论里用的。你的定理能走多远,取决于你定义「迹」时那个求和收敛不收敛。
6.3 迹范数(本讲之外的一句预告) #
如果你手上有个算子 T T T ,它的奇异值是 σ 1 ≥ ⋯ ≥ σ n ≥ 0 \sigma_1\ge\cdots\ge\sigma_n\ge 0 σ 1 ≥ ⋯ ≥ σ n ≥ 0 (29 讲),那么
∥ T ∥ tr : = tr T ∗ T = σ 1 + ⋯ + σ n \|T\|_{\text{tr}}:=\operatorname{tr}\sqrt{T^*T}=\sigma_1+\cdots+\sigma_n ∥ T ∥ tr := tr T ∗ T = σ 1 + ⋯ + σ n
也是一个「大小」,叫迹范数 (核范数)。它和 3.9 的 Frobenius 范数的区别是:F 范数是 σ 1 2 + ⋯ + σ n 2 \sqrt{\sigma_1^2+\cdots+\sigma_n^2} σ 1 2 + ⋯ + σ n 2 (平方和开根),迹范数是 σ i \sigma_i σ i 直接相加。前者惩罚大奇异值,后者一视同仁。压缩感知与矩阵补全的理论里,迹范数扮演「秩的凸替代」的角色。
上面的迹范数不在本讲主线的原文范围内,是从「T ∗ T \sqrt{T^*T} T ∗ T 的迹」这个自然延伸推出来的。 你只要记住 tr \operatorname{tr} tr 在这里的作用:把「对角元」换成「奇异值」再求和,就得到一种新的范数。
6.4 一句群表示论 #
若 G G G 是群,ρ : G → G L ( V ) \rho:G\to GL(V) ρ : G → G L ( V ) 是表示,则 χ ( g ) = tr ρ ( g ) \chi(g)=\operatorname{tr}\rho(g) χ ( g ) = tr ρ ( g ) 叫特征标 (character)。它满足 χ ( h g h − 1 ) = χ ( g ) \chi(hgh^{-1})=\chi(g) χ ( h g h − 1 ) = χ ( g ) ——又是因为循环性。整个表示论的第一层分类,就是靠这个「循环不变的迹」把共轭类上的值记下来。迹之所以到处都是,是因为「循环不变」正是「与基无关」的代数版。
七、自测(先做,再看答案) #
1. 设 A = ( 1 2 0 2 1 0 0 0 5 ) A=\begin{pmatrix}1&2&0\\ 2&1&0\\ 0&0&5\end{pmatrix} A = 1 2 0 2 1 0 0 0 5 。(a) 求 tr A \operatorname{tr}A tr A ;(b) 已知 A A A 有两个特征值是 3 3 3 与 − 1 -1 − 1 ,不用算特征多项式,求出第三个。
2. 设 A A A 为 2 × 3 2\times3 2 × 3 矩阵、B B B 为 3 × 2 3\times2 3 × 2 矩阵。证明 tr ( A B ) = tr ( B A ) \operatorname{tr}(AB)=\operatorname{tr}(BA) tr ( A B ) = tr ( B A ) ,并指出此时 A B AB A B 与 B A BA B A 的阶数。
3. 设 T ∈ L ( C 3 ) T\in\mathcal{L}(\mathbb{C}^3) T ∈ L ( C 3 ) 的矩阵是 ( 51 − 12 − 21 60 − 40 − 28 57 − 68 1 ) \begin{pmatrix}51&-12&-21\\ 60&-40&-28\\ 57&-68&1\end{pmatrix} 51 60 57 − 12 − 40 − 68 − 21 − 28 1 。有人告诉你 − 48 -48 − 48 与 24 24 24 是它的特征值。不许动笔算出第三个。
4. 证明:不存在 S , T ∈ L ( C 4 ) S,T\in\mathcal{L}(\mathbb{C}^4) S , T ∈ L ( C 4 ) 使 S T − T S = I ST-TS=I S T − T S = I 。
5. 设 P ∈ L ( V ) P\in\mathcal{L}(V) P ∈ L ( V ) 满足 P 2 = P P^2=P P 2 = P (幂等)。证明 tr P = dim range P \operatorname{tr}P=\dim\operatorname{range}P tr P = dim range P 。
答案(点开) 1. (a) 对角元是 1 , 1 , 5 1,1,5 1 , 1 , 5 ,故 tr A = 1 + 1 + 5 = 7 \operatorname{tr}A=1+1+5=7 tr A = 1 + 1 + 5 = 7 。 (b) 由 3.4,三个特征值之和等于迹:3 + ( − 1 ) + λ = 7 3+(-1)+\lambda=7 3 + ( − 1 ) + λ = 7 ,得 λ = 5 \lambda=5 λ = 5 。自己验一遍 :A A A 是分块上三角的 ( B 0 0 5 ) \begin{pmatrix}B&0\\ 0&5\end{pmatrix} ( B 0 0 5 ) (B = ( 1 2 2 1 ) B=\begin{pmatrix}1&2\\ 2&1\end{pmatrix} B = ( 1 2 2 1 ) ), 故特征多项式 = det ( B − z I ) ⋅ ( 5 − z ) =\det(B-zI)\cdot(5-z) = det ( B − z I ) ⋅ ( 5 − z ) 。而 det ( B − z I ) = ( 1 − z ) 2 − 4 = z 2 − 2 z − 3 = ( z − 3 ) ( z + 1 ) \det(B-zI)=(1-z)^2-4=z^2-2z-3=(z-3)(z+1) det ( B − z I ) = ( 1 − z ) 2 − 4 = z 2 − 2 z − 3 = ( z − 3 ) ( z + 1 ) , 全谱是 3 , − 1 , 5 3,-1,5 3 , − 1 , 5 ,和 3 − 1 + 5 = 7 3-1+5=7 3 − 1 + 5 = 7 ,与迹一致 ✓。这条技巧在考试里极常用 :给了 n − 1 n-1 n − 1 个特征值,第 n n n 个白送,不必再碰特征多项式。 顺带一个反例提醒:( 1 2 2 1 ) \begin{pmatrix}1&2\\ 2&1\end{pmatrix} ( 1 2 2 1 ) 的迹是 2 2 2 ,特征值是 3 , − 1 3,-1 3 , − 1 (和也是 2 2 2 ); 若你把它误算成 ± 1 \pm1 ± 1 (那是 ( 0 1 1 0 ) \begin{pmatrix}0&1\\ 1&0\end{pmatrix} ( 0 1 1 0 ) 的谱),和就变成 0 ≠ 2 0\ne2 0 = 2 ——迹一算不等,错在哪立刻暴露。这就是第 1.1 节说的「一行字校验」。
2. A B AB A B 是 2 × 2 2\times2 2 × 2 ,B A BA B A 是 3 × 3 3\times3 3 × 3 。设 A = ( a j , k ) A=(a_{j,k}) A = ( a j , k ) (1 ≤ j ≤ 2 1\le j\le2 1 ≤ j ≤ 2 ,1 ≤ k ≤ 3 1\le k\le3 1 ≤ k ≤ 3 )、B = ( b k , j ) B=(b_{k,j}) B = ( b k , j ) 。则tr ( A B ) = ∑ j = 1 2 ∑ k = 1 3 a j , k b k , j \operatorname{tr}(AB)=\sum_{j=1}^{2}\sum_{k=1}^{3}a_{j,k}b_{k,j} tr ( A B ) = ∑ j = 1 2 ∑ k = 1 3 a j , k b k , j ,而 tr ( B A ) = ∑ k = 1 3 ∑ j = 1 2 b k , j a j , k \operatorname{tr}(BA)=\sum_{k=1}^{3}\sum_{j=1}^{2}b_{k,j}a_{j,k} tr ( B A ) = ∑ k = 1 3 ∑ j = 1 2 b k , j a j , k 。 有限双重和可以换序,故两者相等。注意:两个矩阵阶数不同,不可能相等 ,但它们的迹可以相等——这两件事互不矛盾。
3. 迹是 51 + ( − 40 ) + 1 = 12 51+(-40)+1=12 51 + ( − 40 ) + 1 = 12 。三个特征值之和为 12 12 12 ,故第三个是 12 − ( − 48 ) − 24 = 36 12-(-48)-24=36 12 − ( − 48 ) − 24 = 36 。 (这条正是「迹是免费校验」的实战:特征值 − 48 , 24 , 36 -48,24,36 − 48 , 24 , 36 都不好猜,但它们的和一眼可读。)
4. 若存在,则取迹:tr ( S T − T S ) = tr ( S T ) − tr ( T S ) = 0 \operatorname{tr}(ST-TS)=\operatorname{tr}(ST)-\operatorname{tr}(TS)=0 tr ( S T − T S ) = tr ( S T ) − tr ( T S ) = 0 , 而 tr I = dim C 4 = 4 ≠ 0 \operatorname{tr}I=\dim\mathbb{C}^4=4\ne0 tr I = dim C 4 = 4 = 0 。一个算子迹为 0 0 0 、另一个迹为 4 4 4 ,不可能相等(同一个算子只有一个迹)。注意结论对维数没有要求 ,任何有限正维数都行;关键是 dim V ≠ 0 \dim V\ne 0 dim V = 0 。
5. 思路:先证 P P P 可对角化(于是对角元就是特征值),再数有几个 1 1 1 。 ① P 2 = P P^2=P P 2 = P 即 P ( P − I ) = 0 P(P-I)=0 P ( P − I ) = 0 ,说明极小多项式整除 z ( z − 1 ) z(z-1) z ( z − 1 ) 。它无重根,故由 20 讲的极小多项式判据,P P P 可对角化。 于是存在一组基使 P P P 的矩阵是对角阵,对角元 λ 1 , … , λ n \lambda_1,\dots,\lambda_n λ 1 , … , λ n 各自只能取 0 0 0 或 1 1 1 (因为它们是特征值,而特征值必是极小多项式的零点)。 ② 由 3.4(迹与基无关),tr P = λ 1 + ⋯ + λ n = # { j : λ j = 1 } \operatorname{tr}P=\lambda_1+\cdots+\lambda_n=\#\{j:\lambda_j=1\} tr P = λ 1 + ⋯ + λ n = # { j : λ j = 1 } 。 ③ 另一方面 P P P 在 range P \operatorname{range}P range P 上是恒等、在 null P \operatorname{null}P null P 上是零: 前者因为 range P \operatorname{range}P range P 中每个向量都形如 P v Pv P v ,而 P ( P v ) = P 2 v = P v P(Pv)=P^2v=Pv P ( P v ) = P 2 v = P v ; 后者按定义。再注意 V = null P ⊕ range P V=\operatorname{null}P\oplus\operatorname{range}P V = null P ⊕ range P (把 v v v 拆成 ( v − P v ) + P v (v-Pv)+Pv ( v − P v ) + P v , 前一项被 P P P 打成 P v − P 2 v = 0 Pv-P^2v=0 P v − P 2 v = 0 )总成立。所以 1 1 1 这个特征值的特征空间恰是 range P \operatorname{range}P range P , 其维数 dim range P \dim\operatorname{range}P dim range P 就等于「λ j = 1 \lambda_j=1 λ j = 1 的个数」。 合起来 tr P = dim range P \operatorname{tr}P=\dim\operatorname{range}P tr P = dim range P 。■ \blacksquare ■ (投影的迹 = 它留下的维数——正是 2.4 那张表里第三行的观察。)
八、复习卡 #
项
内容
矩阵的迹
方阵对角元之和 tr A = ∑ j A j , j \operatorname{tr}A=\sum_j A_{j,j} tr A = ∑ j A j , j
算子的迹
取任意一组基,算矩阵的迹;由 3.2 与基无关
循环性
tr ( A B ) = tr ( B A ) \operatorname{tr}(AB)=\operatorname{tr}(BA) tr ( A B ) = tr ( B A ) ,A A A 为 m × n m\times n m × n 、B B B 为 n × m n\times m n × m (不必方、不必同形)
允许的变形
循环移位:tr ( A B C ) = tr ( C A B ) = tr ( B C A ) \operatorname{tr}(ABC)=\operatorname{tr}(CAB)=\operatorname{tr}(BCA) tr ( A B C ) = tr ( C A B ) = tr ( B C A )
不允许的
任意对调(tr ( A B C ) ≠ tr ( A C B ) \operatorname{tr}(ABC)\ne\operatorname{tr}(ACB) tr ( A B C ) = tr ( A C B ) 一般);tr ( A B ) ≠ tr A ⋅ tr B \operatorname{tr}(AB)\ne\operatorname{tr}A\cdot\operatorname{tr}B tr ( A B ) = tr A ⋅ tr B
基无关的证明
A = C − 1 B C A=C^{-1}BC A = C − 1 B C ,把 C − 1 C^{-1} C − 1 一路右推到 C C C 旁边消掉
迹 = 特征值和
tr T = ∑ λ i \operatorname{tr}T=\sum\lambda_i tr T = ∑ λ i (含重数),需 F = C \mathbb{F}=\mathbb{C} F = C 或特征多项式完全分解
特征多项式
tr T = − [ z n − 1 的系数 ] \operatorname{tr}T=-[z^{n-1}\text{ 的系数}] tr T = − [ z n − 1 的系数 ] ;det T = ( − 1 ) n [ 常数项 ] \det T=(-1)^n[\text{常数项}] det T = ( − 1 ) n [ 常数项 ]
线性泛函
tr ( S + T ) = tr S + tr T \operatorname{tr}(S+T)=\operatorname{tr}S+\operatorname{tr}T tr ( S + T ) = tr S + tr T ,tr ( λ T ) = λ tr T \operatorname{tr}(\lambda T)=\lambda\operatorname{tr}T tr ( λ T ) = λ tr T
唯一刻画
循环不变 + tr I = dim V \operatorname{tr}I=\dim V tr I = dim V ⇒ \Rightarrow ⇒ 只能是迹
内积公式
规范正交基下 tr T = ∑ ⟨ T e k , e k ⟩ \operatorname{tr}T=\sum\langle Te_k,e_k\rangle tr T = ∑ ⟨ T e k , e k ⟩
不可能性
有限维无解 S T − T S = I ST-TS=I S T − T S = I (迹左 0 0 0 、右 dim V \dim V dim V )
迹内积
⟨ S , T ⟩ = tr ( T ∗ S ) \langle S,T\rangle=\operatorname{tr}(T^*S) ⟨ S , T ⟩ = tr ( T ∗ S ) 是 L ( V , W ) \mathcal{L}(V,W) L ( V , W ) 上的内积
Frobenius 范数
tr ( T ∗ T ) = ∑ j , k ∣ a j , k ∣ 2 = ∑ σ i 2 \operatorname{tr}(T^*T)=\sum_{j,k}|a_{j,k}|^2=\sum\sigma_i^2 tr ( T ∗ T ) = ∑ j , k ∣ a j , k ∣ 2 = ∑ σ i 2
转置共轭
tr ( T ∗ ) = tr T ‾ \operatorname{tr}(T^*)=\overline{\operatorname{tr}T} tr ( T ∗ ) = tr T ;故自伴算子的迹为实数
零迹的两种命运
幂零 ⇒ \Rightarrow ⇒ 迹零(反之不然);正算子且迹零 ⇒ \Rightarrow ⇒ 算子零
幂等的迹
P 2 = P P^2=P P 2 = P ⇒ \Rightarrow ⇒ tr P = dim range P \operatorname{tr}P=\dim\operatorname{range}P tr P = dim range P
与行列式的区别
迹可加不可乘;行列式可乘。可逆性归行列式,求和与期望归迹
九、出处与延伸(可选,不影响阅读) #
主线出处 :上述定义与结论系统对应《线性代数应该这样学》(第 4 版)第 8 章 8D 节「联系矩阵与算子的桥梁——迹」,PDF 页码 p286–290(印刷页码 p272–276)。原文编号:8.47(矩阵的迹)、8.48(3 × 3 3\times3 3 × 3 例)、8.49(tr ( A B ) = tr ( B A ) \operatorname{tr}(AB)=\operatorname{tr}(BA) tr ( A B ) = tr ( B A ) )、8.50(迹与基无关)、8.51(算子的迹)、8.52(复空间上迹等于特征值之和)、8.53(C 3 \mathbb{C}^3 C 3 例)、8.54(迹与特征多项式)、8.55(内积空间上的迹)、8.56(迹是线性的)、8.57(S T − T S ≠ I ST-TS\ne I S T − T S = I )。
本讲的顺序调整 :原文先讲矩阵、再证基无关、最后转向算子;本讲把「基无关」的动机(第 1.1、1.2 节)提前,并把「为什么必须先证基无关才能定义」放在 2.3 明说,因为这是初学者最容易糊过去的一步。逻辑关系未变。
习题来源 :自测第 3 题、第 5 题以及 3.9 的迹内积、3.10 的正算子结论,分别对应原文 8D 节习题 7、2、12、5。第 3 题的三个特征值是 − 48 , 24 , 36 -48,24,36 − 48 , 24 , 36 。
迹与行列式并列的展开 :第 9 章(p291 起)用交错多重线性型重讲行列式,并证明「行列式 = 特征值之积」。想对照「和 vs 积」这条对称性,那里是最完整的版本。
数值视角 :任何线性代数库都有 trace,成本 O ( n ) O(n) O ( n ) 。用它校验特征值分解(eig 的输出求和与 trace 比对)是工程上最常见的一行自查。
术语别名 :「迹」= trace = Spur(德)= след(俄);「线性泛函」= linear functional;「特征标」= character;「Frobenius 范数」= Hilbert–Schmidt 范数。
想深挖时可查 :Horn & Johnson《Matrix Analysis》里关于迹的不等式(如 ∣ tr ( A ∗ B ) ∣ 2 ≤ tr ( A ∗ A ) tr ( B ∗ B ) \big|\operatorname{tr}(A^*B)\big|^2\le\operatorname{tr}(A^*A)\operatorname{tr}(B^*B) tr ( A ∗ B ) 2 ≤ tr ( A ∗ A ) tr ( B ∗ B ) ,它就是 3.9 那个内积下的 Cauchy–Schwarz 不等式);量子力学教材里 tr ( ρ A ) \operatorname{tr}(\rho A) tr ( ρ A ) 的推导。
迹的全部魔力集中在一条等式上:tr ( A B ) = tr ( B A ) \operatorname{tr}(AB)=\operatorname{tr}(BA) tr ( A B ) = tr ( B A ) 。 后面你遇到的每个「迹的妙用」,剥到最后一层都是它。 若只记一件事,记这条,并且记住它的证明只有一步——双重和换序 。