|
回复 1# abababa
这个是不是错的啊,应该是$P^TAP=E$,其中$E$是单位矩阵。我一直没证明出$A=P^TP$。
下面是我的证明:
因为$A$正定实对称,所以能被正交对角化,即存在正交矩阵$C$使得$A=C^T{\rm{diag}(c_1^2,\cdots,c_n^2)}C$,因为正交矩阵转置也正交,所以这里不妨变换$C^T \to C$,即$A=C{\rm{diag}(c_1^2,\cdots,c_n^2)}C^T$。然后$A^{-1}=C{\rm{diag}}(1/c_1^2,\cdots,1/c_n^2)C^T$,记$\sqrt{A^{-1}}=C{\rm{diag}}(1/c_1,\cdots,1/c_n)C^T$.
因为$\sqrt{A^{-1}}^TB\sqrt{A^{-1}}$也是实对称矩阵,所以也能正交对角化为$M^T\sqrt{A^{-1}}^TB\sqrt{A^{-1}}M={\rm{diag}}(\lambda_1,\cdots,\lambda_n)$,其中$M$是正交矩阵。然后令$P=\sqrt{A^{-1}}M$,就有$P^TBP=\rm{diag}(\lambda_1,\cdots,\lambda_n)$。
这时
\begin{align*}
P^TAP &= \left(M^T\sqrt{A^{-1}}^T\right)\left(C{\rm{diag}}(c_1^2,\cdots,c_n^2)C^T\right)\left(\sqrt{A^{-1}}M\right)\\
&= M^T\left(C{\rm{diag}}(\frac{1}{c_1},\cdots,\frac{1}{c_n})C^T\right)\left(C{\rm{diag}}(c_1^2,\cdots,c_n^2)C^T\right)\left(C{\rm{diag}}(\frac{1}{c_1},\cdots,\frac{1}{c_n})C^T\right)M\\
&= M^T\left(C{\rm{diag}}(1,\cdots,1)C^T\right)M = M^TEM = E
\end{align*} |
|