实对称矩阵的对角化

实对称矩阵的对角化,我们一般是找一个正交矩阵 \(Q\),使得 \(Q^{-1}AQ=D\),所以比普通的矩阵对角化多一个步骤:施密特正交化步骤。

我们先叙述实对称矩阵对角化的结论。

1,定理: \(A\) 是实对称矩阵,则一定存在一个正交矩阵 \(Q\),使得 \(Q^{-1}AQ=D\) 或者 \(Q^TAQ=D\)。

这个定理我们这里不证明了。

2,实对称矩阵对角化的步骤:

  • 求特征值;
  • 求特征向量;
  • 将同一个特征值所对应的不同特征向量正交化(施密特正交化方法);
  • 将所有正交特征向量规范化
  • 得到 \(Q\) 与 \(D\)。

我们来看例题。

例1:设 \(A=\begin{pmatrix}4&2&2\\2&4&2\\2&2&4\end{pmatrix}\),求正交矩阵 \(Q\),使得 \(Q^{-1}AQ=D\)。

解:(1)先求特征值。

\[\begin{align*}|A-\lambda I|&=\begin{vmatrix}4-\lambda&2&2\\2&4-\lambda&2\\2&2&4-\lambda\end{vmatrix}=\begin{vmatrix}8-\lambda&8-\lambda&8-\lambda\\2&4-\lambda&2\\2&2&4-\lambda\end{vmatrix}\\ &=(8-\lambda)\begin{vmatrix}1&1&1\\2&4-\lambda&2\\2&2&4-\lambda\end{vmatrix}=(8-\lambda)\begin{vmatrix}1&1&1\\0&2-\lambda&2\\0&0&2-\lambda\end{vmatrix}\\ &=(8-\lambda)(2-\lambda)^2=0\end{align*}\] 我们得到特征值 \( \lambda_1=8,\lambda_{2,3}=2\)。

(2)再来求特征向量。

当 \(\lambda=8\) 时,

\[\begin{align*}A-\lambda I&=\begin{pmatrix}-4&2&2\\2&-4&2\\2&2&-4\end{pmatrix}\sim\begin{pmatrix}-4&2&2\\2&-4&2\\0&0&0\end{pmatrix}\\ &\sim \begin{pmatrix}2&-4&2\\-4&2&2\\0&0&0\end{pmatrix}\sim \begin{pmatrix}2&-4&2\\0&-6&6\\0&0&0\end{pmatrix}\\ &\sim \begin{pmatrix}2&-4&2\\0&1&-1\\0&0&0\end{pmatrix}\sim \begin{pmatrix}1&0&-1\\0&1&-1\\0&0&0\end{pmatrix}\end{align*}\]

特征向量为 \(\vec{x}_1=\begin{pmatrix}1\\1\\1\end{pmatrix}\).

当 \(\lambda=2\) 时,

\[\begin{align*}A-\lambda I&=\begin{pmatrix}2&2&2\\2&2&2\\2&2&2\end{pmatrix}\sim \begin{pmatrix}1&1&1\\0&0&0\\0&0&0\end{pmatrix}\end{align*}\]

特征向量为 \(\vec{x}_2=\begin{pmatrix}-1\\1\\0\end{pmatrix}, \vec{x}_3=\begin{pmatrix}-1\\0\\1\end{pmatrix}\)。

(3)将 \(\vec{x}_1,\vec{x}_2, \vec{x}_3\) 规范正交化。 令\(\vec{b}_1=\vec{x}_1, \|\vec{b}_1\|=\sqrt3\),因为 \(\vec{x}_1\) 与 \(\vec{x}_2,\vec{x}_3\) 正交,我们只需要对它规范化就可以了。 \[\vec{e}_1=\frac{\vec{b}_1}{\|\vec{b}_1\|}=\begin{pmatrix}\frac{1}{\sqrt3}\\ \frac{1}{\sqrt3}\\ \frac{1}{\sqrt3}\end{pmatrix}\]

我们对 \(\vec{x}_2, \vec{x}_3\) 规范正交化。 令 \(\vec{b_2}=\vec{x}_2, \|\vec{b}_2\|=\sqrt2\),单位化可得

\[\vec{e}_2=\frac{\vec{b}_2}{\|\vec{b}_2\|}=\begin{pmatrix}-\frac{1}{\sqrt2}\\ \frac{1}{\sqrt2}\\ 0\end{pmatrix}\]

\(\vec{x}_3\cdot\vec{b}_2=1, \|\vec{b}_2\|^2=2\),所以

\[\vec{b_3}=\vec{x}_3-\frac{\vec{x}_3\cdot\vec{b}_2}{\|\vec{b}_2\|^2}\vec{b}_2=\begin{pmatrix}-1\\0\\1\end{pmatrix}-\frac{1}{2}\begin{pmatrix}-1\\1\\0\end{pmatrix}=\begin{pmatrix}-\frac{1}{2}\\-\frac{1}{2}\\1\end{pmatrix}\]

\(\|\vec{b}_3\|=\sqrt{\frac{3}{2}}\),将 \(\vec{b}_3\) 单位化。

\[\vec{e}^3=\frac{\vec{b}_3}{\|\vec{b}_3\|}=\begin{pmatrix}\frac{1}{\sqrt6}\\ -\frac{2}{\sqrt6}\\\frac{1}{\sqrt6}\end{pmatrix}\]

(4)对角化:我们令

\[Q=(\vec{e}_1,\vec{e}_2,\vec{e}_3)=\begin{pmatrix}-\frac{1}{\sqrt2}&\frac{1}{\sqrt3}&\frac{1}{\sqrt2}\\ \frac{1}{\sqrt2}&\frac{1}{\sqrt3}&-\frac{2}{\sqrt6}\\ 0&\frac{1}{\sqrt3}&\frac{1}{\sqrt6}\end{pmatrix},\quad D=\begin{pmatrix}8&&\\ &2&\\ &&2\end{pmatrix}\]

则 \(Q^{-1}AQ=D\)。