求解线性方程组(2 & -1 & 1 1 & 1 & 1 1 & 1 & -2)的Jacobi方法迭代矩阵的谱半径等于() A (3)/(2) B (sqrt(3))/(2) C (5)/(2) D (sqrt(5))/(2)
求解线性方程组$\left(\begin{array}{cccccc}2 & -1 & 1 \\ 1 & 1 & 1 \\ 1 & 1 & -2\end{array}\right)\left(\begin{array}{cccccc}x_1 \\ x_2 \\ x_3\end{array}\right)=\left(\begin{array}{cccccc}1 \\ 4 \\ 5\end{array}\right)$的Jacobi方法迭代矩阵的谱半径等于()
A $\frac{3}{2}$
B $\frac{\sqrt{3}}{2}$
C $\frac{5}{2}$
D $\frac{\sqrt{5}}{2}$
题目解答
答案
解析
考查要点:本题主要考查Jacobi迭代法迭代矩阵的构造及其谱半径的计算。
解题思路:
- 分解系数矩阵:将系数矩阵分解为对角矩阵$D$、严格下三角矩阵$L$和严格上三角矩阵$U$。
- 构造迭代矩阵:根据公式$B = -D^{-1}(L + U)$,计算Jacobi迭代矩阵。
- 求特征值:通过求解特征方程$\det(B - \lambda I) = 0$,得到迭代矩阵的特征值。
- 计算谱半径:取特征值模的最大值作为谱半径。
关键点:
- 矩阵分解的正确性直接影响迭代矩阵的构造。
- 特征值的计算需注意行列式的展开和化简,避免符号错误。
- 谱半径的定义是特征值模的最大值,需特别注意虚数特征值的模。
1. 分解系数矩阵
系数矩阵$A$为:
$A = \begin{pmatrix}2 & -1 & 1 \\1 & 1 & 1 \\1 & 1 & -2\end{pmatrix}$
分解为:
- 对角矩阵$D$:
$D = \begin{pmatrix}2 & 0 & 0 \\0 & 1 & 0 \\0 & 0 & -2\end{pmatrix}$ - 严格下三角矩阵$L$:
$L = \begin{pmatrix}0 & 0 & 0 \\1 & 0 & 0 \\1 & 1 & 0\end{pmatrix}$ - 严格上三角矩阵$U$:
$U = \begin{pmatrix}0 & -1 & 1 \\0 & 0 & 1 \\0 & 0 & 0\end{pmatrix}$
2. 构造Jacobi迭代矩阵
计算$D^{-1}$:
$D^{-1} = \begin{pmatrix}\frac{1}{2} & 0 & 0 \\0 & 1 & 0 \\0 & 0 & -\frac{1}{2}\end{pmatrix}$
计算$L + U$:
$L + U = \begin{pmatrix}0 & -1 & 1 \\1 & 0 & 1 \\1 & 1 & 0\end{pmatrix}$
计算$D^{-1}(L + U)$:
$D^{-1}(L + U) = \begin{pmatrix}0 & -\frac{1}{2} & \frac{1}{2} \\1 & 0 & 1 \\-\frac{1}{2} & -\frac{1}{2} & 0\end{pmatrix}$
迭代矩阵$B$为:
$B = -D^{-1}(L + U) = \begin{pmatrix}0 & \frac{1}{2} & -\frac{1}{2} \\-1 & 0 & -1 \\\frac{1}{2} & \frac{1}{2} & 0\end{pmatrix}$
3. 求特征值
构造特征方程$\det(B - \lambda I) = 0$:
$\det\begin{pmatrix}-\lambda & \frac{1}{2} & -\frac{1}{2} \\-1 & -\lambda & -1 \\\frac{1}{2} & \frac{1}{2} & -\lambda\end{pmatrix} = 0$
展开行列式并化简得:
$-\lambda^3 - \frac{5}{4}\lambda = 0 \implies \lambda(\lambda^2 + \frac{5}{4}) = 0$
解得特征值:
$\lambda_1 = 0, \quad \lambda_{2,3} = \pm \frac{\sqrt{5}}{2}i$
4. 计算谱半径
特征值的模为:
$|\lambda_1| = 0, \quad |\lambda_{2,3}| = \frac{\sqrt{5}}{2}$
谱半径为$\frac{\sqrt{5}}{2}$,对应选项D。