八、(本題10分) 設 $A,B,C$ 均為 $n$ 階半正定實對稱陣, 使得 $ABC$ 是對稱陣, 即滿足 $ABC=CBA$. 證明: $ABC$ 也是半正定陣.
證明 我們先引用如下引理, 它是前 4 種證法的基礎, 這也是本題中兩個半正定陣的特殊情形.
引理 1 (高代白皮書例 9.55) 設 $A,B$ 都是 $n$ 階半正定實對稱陣, 則 $AB$ 的所有特征值都是非負實數. 特別地, 若 $AB$ 是對稱陣, 即 $AB=BA$, 則 $AB$ 也是半正定陣.
證法 1 (攝動法+數學歸納法) Case 1: 若 $C$ 為正定陣, 則 $ABC$ 合同於 $C^{-\frac{1}{2}}(ABC)C^{-\frac{1}{2}}=(C^{-\frac{1}{2}}AC^{-\frac{1}{2}})(C^{\frac{1}{2}}BC^{\frac{1}{2}})$. 注意到 $C^{-\frac{1}{2}}AC^{-\frac{1}{2}}$ 和 $C^{\frac{1}{2}}BC^{\frac{1}{2}}$ 都是半正定實對稱陣並且它們的乘積還是對稱陣, 故由引理 1 可知它們的乘積是半正定陣, 從而 $ABC$ 也是半正定陣.
Case 2: 若 $Ax=Cx=0$ 沒有公共非零解, 則對任意的正實數 $t$, $C+tA$ 都是正定陣. (事實上, $C+tA$ 是半正定陣. 若 $r(C+tA)<n$, 則存在非零實列向量 $\alpha$, 使得 $(C+tA)\alpha=0$, 從而 $0=\alpha'(C+tA)\alpha=\alpha'C\alpha+t\alpha'A\alpha$. 由 $A,C$ 的半正定性即得 $\alpha' A\alpha=\alpha' C\alpha=0$, 再由白皮書例 8.44 可知 $A\alpha=C\alpha=0$, 這與假設矛盾. 因此 $r(C+tA)=n$, 從而 $C+tA$ 為正定陣. 我們也可用白皮書例 9.73 來證明這一結論.) 注意到 $AB(C+tA)$ 仍為對稱陣, 故由 Case 1 可得 $AB(C+tA)$ 是半正定陣, 令 $t\to 0+$, 則有 $ABC$ 也是半正定陣.
Case 3: 若 $Ax=Cx=0$ 有公共非零解, 不妨設為單位特征向量 $\alpha_1$, 擴張為 $\mathbb{R}^n$ (取標准內積) 的一組標准正交基 $\{\alpha_1,\alpha_2,\cdots,\alpha_n\}$. 令 $P=(\alpha_1,\alpha_2,\cdots,\alpha_n)$, 則 $P$ 為正交陣且有 $$P'AP=\begin{pmatrix} 0 & 0 \\ 0 & A_1 \\ \end{pmatrix},\,\,\,\,P'BP=\begin{pmatrix} * & * \\ * & B_1 \\ \end{pmatrix},\,\,\,\,P'CP=\begin{pmatrix} 0 & 0 \\ 0 & C_1 \\ \end{pmatrix},$$ 於是 $P'ABCP=\begin{pmatrix} 0 & 0 \\ 0 & A_1B_1C_1 \\ \end{pmatrix}$. 由 $ABC=CBA$ 可知 $A_1B_1C_1=C_1B_1A_1$, 又 $A_1,B_1,C_1$ 為半正定陣, 故若證明開始處對階數 $n$ 進行歸納, 則由歸納假設可知 $A_1B_1C_1$ 為半正定陣, 從而 $ABC$ 也是半正定陣.
引理 2 (高代白皮書例 9.72) 設 $A,C$ 都是 $n$ 階半正定實對稱陣, 則存在非異實方陣 $D$, 使得 $$D'AD=\mathrm{diag}\{I_r,0\},\,\,\,\,D'CD=\mathrm{diag}\{\lambda_1,\lambda_2,\cdots,\lambda_n\}.$$
證法 2 (同時合同對角化) 由引理 2 可知兩個半正定陣可同時合同對角化. 由於題目的條件和結論在合同變換 $A\to D'AD$, $B\to D^{-1}B(D^{-1})'$, $C\to D'CD$ 下不改變, 故不妨從一開始就假設 $A,C$ 就是引理 2 中的兩個對角陣. 將 $B,C$ 按照 $A$ 的分塊方式進行對應的分塊: $B=\begin{pmatrix} B_{11} & B_{12} \\ B_{12}' & B_{22} \\ \end{pmatrix}$, $C=\begin{pmatrix} \Lambda_1 & 0 \\ 0 & \Lambda_2 \\ \end{pmatrix}$, 其中 $\Lambda_1,\Lambda_2$ 是半正定對角陣. 於是 $ABC=\begin{pmatrix} B_{11}\Lambda_1 & B_{12}\Lambda_2 \\ 0 & 0 \\ \end{pmatrix}$, 由 $ABC$ 的對稱性可得 $B_{12}\Lambda_2=0$, 並且 $B_{11}\Lambda_1$ 是對稱陣. 注意到 $B_{11},\Lambda_1$ 都是半正定陣, 故由引理 1 知 $B_{11}\Lambda_1$ 是半正定陣, 從而 $ABC=\mathrm{diag}\{B_{11}\Lambda_1,0\}$ 也是半正定陣.
證法 3 (半正定陣的性質--主對角塊占優) 由於題目的條件和結論在合同變換 $A\to D'AD$, $B\to D^{-1}B(D^{-1})'$, $C\to D'CD$ 下不改變, 故不妨從一開始就假設 $B=\begin{pmatrix} I_r & 0 \\ 0 & 0 \\ \end{pmatrix}$ 為合同標准型. 將 $A,C$ 進行對應的分塊: $A=\begin{pmatrix} A_{11} & A_{12} \\ A_{12}' & A_{22} \\ \end{pmatrix}$, $C=\begin{pmatrix} C_{11} & C_{12} \\ C_{12}' & C_{22} \\ \end{pmatrix}$, 則 $ABC=\begin{pmatrix} A_{11}C_{11} & A_{11}C_{12} \\ A_{12}'C_{11} & A_{12}'C_{12} \\ \end{pmatrix}$. 注意到 $C$ 是半正定陣, 故由白皮書例 8.45 可知 $r(C_{11},C_{12})=r(C_{11})$, 再由線性方程組的求解理論可知, 存在實矩陣 $M$, 使得 $C_{11}M=C_{12}$. 注意到 $ABC$ 是對稱陣, 對其實施分塊對稱初等變換可化為: $$\begin{pmatrix} I_r & 0 \\ -M' & I_{n-r} \\ \end{pmatrix}\begin{pmatrix} A_{11}C_{11} & A_{11}C_{12} \\ A_{12}'C_{11} & A_{12}'C_{12} \\ \end{pmatrix}\begin{pmatrix} I_r & -M \\ 0 & I_{n-r} \\ \end{pmatrix}=\begin{pmatrix} A_{11}C_{11} & 0 \\ 0 & 0 \\ \end{pmatrix}.$$ 於是 $A_{11}C_{11}$ 也是對稱陣, 又 $A_{11},C_{11}$ 都是半正定陣, 故由引理 1 知 $A_{11}C_{11}$ 是半正定陣, 從而 $ABC$ 也是半正定陣.
證法 4 由於題目的條件和結論在合同變換 $A\to D'AD$, $B\to D^{-1}B(D^{-1})'$, $C\to D'CD$ 下不改變, 故不妨從一開始就假設 $A=\begin{pmatrix} I_r & 0 \\ 0 & 0 \\ \end{pmatrix}$ 為合同標准型. 設 $BC=\begin{pmatrix} M_{11} & M_{12} \\ M_{21} & M_{22} \\ \end{pmatrix}$, $CB=\begin{pmatrix} N_{11} & N_{12} \\ N_{21} & N_{22} \\ \end{pmatrix}$ 為對應的分塊, 則由 $A(BC)=(CB)A$ 可得 $$\begin{pmatrix} M_{11} & M_{12} \\ 0 & 0 \\ \end{pmatrix}=\begin{pmatrix} N_{11} & 0 \\ N_{21} & 0 \\ \end{pmatrix},$$ 於是 $M_{11}=N_{11}$, $M_{12}=0$, $N_{21}=0$. 特別地, $ABC=\mathrm{diag}\{M_{11},0\}$ 為對稱陣, 並且 $BC=\begin{pmatrix} M_{11} & 0 \\ M_{21} & M_{22} \\ \end{pmatrix}$ 為分塊下三角陣. 由引理 1 可知, $BC$ 的特征值全部大於等於零, 於是實對稱陣 $M_{11}$ 的特征值也全部大於等於零, 故 $M_{11}$ 是半正定陣, 從而 $ABC$ 也是半正定陣.
引理 3 (高代白皮書例 9.107) 設 $A,C$ 是兩個實對稱陣且滿足 $AC=CA$, 則存在正交陣 $P$, 使得 $P'AP$ 和 $P'CP$ 都是對角陣.
證法 5 ($B$ 正定的特殊情形, 同時正交對角化) 由 $ABC=CBA$ 可得 $(AB^{\frac{1}{2}})(B^{\frac{1}{2}}C)=(CB^{\frac{1}{2}})(B^{\frac{1}{2}}A)$, 進一步可得 $(B^{\frac{1}{2}}AB^{\frac{1}{2}}) (B^{\frac{1}{2}}CB^{\frac{1}{2}})=(B^{\frac{1}{2}}CB^{\frac{1}{2}})(B^{\frac{1}{2}}AB^{\frac{1}{2}})$. 由引理 3 可知, 存在正交陣 $P$, 使得 $$P'(B^{\frac{1}{2}}AB^{\frac{1}{2}})P=\mathrm{diag}\{\lambda_1,\lambda_2,\cdots,\lambda_n\},\,\,\,\,P'(B^{\frac{1}{2}}CB^{\frac{1}{2}})P=\mathrm{diag}\{\mu_1,\mu_2,\cdots,\mu_n\},$$ 由 $A,C$ 的半正定性可知 $\lambda_i\geq 0$, $\mu_i\geq 0\,(1\leq i\leq n)$. 於是 $$P'B^{\frac{1}{2}}(ABC)B^{\frac{1}{2}}P=\mathrm{diag}\{\lambda_1\mu_1,\lambda_2\mu_2,\cdots,\lambda_n\mu_n\}$$ 為半正定陣, 從而 $B^{\frac{1}{2}}(ABC)B^{\frac{1}{2}}$ 也是半正定陣. 當 $B$ 是正定陣時, 可推出 $ABC$ 也是半正定陣. $\Box$
注 1 先給兩個反例. 引理 1 告訴我們: 兩個半正定陣乘積 (未必對稱陣) 的所有特征值都是非負實數, 但 $m\,(m\geq 3)$ 個半正定陣的乘積就不滿足這個性質. 例如: $A_1=\begin{pmatrix} 2 & -1 \\ -1 & 1 \\ \end{pmatrix}$, $A_2=\begin{pmatrix} 1 & 1 \\ 1 & 1 \\ \end{pmatrix}$, $A_3=\begin{pmatrix} 1 & -2 \\ -2 & 5 \\ \end{pmatrix}$ 都是半正定陣, 但 $A_1A_2A_3=\begin{pmatrix} -1 & 3 \\ 0 & 0 \\ \end{pmatrix}$ 的特征值中有負數. 引理 1 和本題結論告訴我們: 兩個或三個半正定陣的乘積, 如果是對稱陣, 那么一定也是半正定陣, 但 $m\,(m\geq 4)$ 個半正定陣的乘積, 即使它是對稱陣, 也不一定是半正定陣. 例如: $A_1,A_2,A_3$ 如上, $A_4=\begin{pmatrix} 10 & 3 \\ 3 & 1 \\ \end{pmatrix}$ 半正定, $A_1A_2A_3A_4=\begin{pmatrix} -1 & 0 \\ 0 & 0 \\ \end{pmatrix}$ 是對稱陣, 但它不是半正定陣.
注 2 上述反例告訴我們: 引理 1 的證明方法 (利用相似變換或降階公式把 $AB$ 的特征值轉化為 $A^{\frac{1}{2}}BA^{\frac{1}{2}}$ 的特征值, 再利用 $B$ 的半正定性即得) 很難推廣到三個半正定陣的情形, 主要原因是合同變換只保持對稱陣的半正定性, 但並不保持非對稱陣的特征值不變 (注意我們只有 $ABC$ 的對稱性, 並沒有 $AB$ 或 $BC$ 的對稱性). 上述反例也告訴我們: 本題的五種證明方法, 即使它們用遍了半正定陣幾乎所有重要的性質和結論, 但也很難推廣到四個或四個以上半正定陣的情形.
注 3 本題的正定版本是: 設 $A,B,C$ 均為 $n$ 階正定實對稱陣, 使得 $ABC$ 是對稱陣, 即滿足 $ABC=CBA$, 則 $ABC$ 也是正定陣. 正定版本可由半正定版本以及矩陣 $ABC$ 的非異性得到, 也可由本題五種證法的正定版本直接得到, 讀者可以做一嘗試.
注 4 本題完全做出 (得分在8分以上) 的同學為: 丁思成 (證法 2), 周爍星 (證法 2), 顧天翊 (證法 2), 王捷翔 (17級, 證法 3). 另外, 封清給出了證法 3 的前半部分, 林潔 (17級) 給出了證法 4 的前半部分, 陳在遠 (17級) 和陳欽品 (17級) 在證法 1 (攝動法) 方面做了嘗試, 馬珩元 (17級) 和王祝斌 (17級) 給出了證法 5 ($B$ 正定的特殊情形).