2012数值分析第二次作业及答案.docVIP

数值分析第二次作业及答案 用矩阵的直接三角分解法（LU分解）解方程组。 矩阵第一行乘以一数，成为,证明当时，有最小值。 3. 设有方程组，其中已知它有解。如果右端有小扰动，试估计由此引起的解的相对误差。 4. 设，已知方程组的精确解为。 计算条件数； 若近似解，计算剩余； 进行误差分析，此结果说明了什么？ 解：（1）于是 （2） （3）依误差公式，右端为 而左端为 这表明当为病态矩阵时，尽管剩余很小，误差估计仍然较大，因此，当为病态时，用大小作为检验解的准确度是不可靠的。 5. 对线性代数方程组 设法导出使雅可比（Jacobi）迭代法和高斯－赛德尔（G-S）迭代法均收敛的迭代格式，要求分别写出迭代格式，并说明收敛的理由。 解： 因其变换后为等价方程组，且严格对角占优，故雅可比和高斯－赛德尔迭代法均收敛。 雅可比迭代格式为： 高斯－赛德尔代格式为： 设方程组 试考察解此方程组的雅可比迭代法及高斯－赛德尔迭代法的收敛性。 用Jacobi、Gause-Seidel迭代法求解下列方程组 问是否收敛？为什么？若将原方程组变为 再用上述两种迭代法求解是否收敛？为什么？ 解：对第一个方程Jacobi迭代法的迭代阵是： 由故Jacobi迭代过程发散。 Gause-Seidel迭代法的迭代阵是： 由故Gause-Seidel迭代过程发散。 对第二个方程，因其系数矩阵是强对角占优的，故再用上述两种迭代均收敛。 8、设有. 证明解此方程组的Jacobi迭代法收敛的充要条件是 的根模 证明解此方程组的Gause-Seidel迭代法收敛的充要条件是 的根模 证：（1）由于Jacobi迭代法的迭代矩阵 于是 而由题设知 ，故Jacobi迭代收敛的充要条件：的根模，等价于的根模。 （2）Gauss-Seidel迭代法的迭代矩阵 于是 由题设知，故Gauss-Seidel迭代收敛的充要条件：的根模，等价于的根模。 9、 10．已知构造一个Householder变换矩阵H，使得。 解： 取， 而 其2－范数 所以 ， Householder变换矩阵H为 。 11．设，（1）把A相似化为拟上三角阵；（2）用QR算法求的（也是A的）全部特征值。（只需一次迭代） 解：（1）取 则有为拟上三角阵。 （2））取 ， 则有 。 再取 ， 则有 。 而 故有 7