6.1 用高斯消去法求解线性方程组 2X1- X2+3X3 = 1 4X1+2X2+5X3 = 4 -3X1+4X2-3X3 = 17
解:x1=-10.4117;x2=2.5294;x3=8.1176。
6.2 用列主元素消元法求解线性方程组 2X1- X2+3X3 = 1 4X1+2X2+5X3 = 4 -3X1+4X2-3X3 = 17
解: x1=-10.4117;x2=2.5294;x3=8.1176。
6.3 用全主元素消元法求解线性方程组 X1+2X2+3X3 = 5 2X1– X2+9X3 = 3 -3X1+ X2–3X3 = 17
解:x1=-6.96;x2=2.71;x3=2.18。
6.4 比较1,2,、3题的消元过程,三种消元的不同点和优点是什么?能消元的条件是什么?
解:高斯消去法:只要对角线上的元素非零就可以进行消元,否则需对换行。
列主元素消元法 :对角线上的元素非零且其绝对值在本列对角线以下元素中最大时就可以进行消元,否则需对换行。
全主元素消元法:对角线上的元素非零且其绝对值在本列(行)对角线以下(后)的所有元素中最大时就可以进行消元,否则需对换行(列)。
用计算机实现解方程组时,列主元素消元法和全主元素消元法可减少计算过程中的舍入误差。
6.5 用矩阵的三角分解算法求解线性方程组 X1+2X2+3X3 = 0 2X1+2X2+8X3 = -4 -3X1-10X2-2X3 = -11
解:分解矩阵ALR有:
231001231210022
28 23102321003解方程组Lyb:
100y100 210y24得y4 321y3113解方程组Rxy:
123x10 022x24得方程组的解x003x33
11 110021210,求A,A。 6.6 对矩阵A101210012解:A为实对称矩阵,行范数A4,列范数A14。
6.7 利用矩阵范数的性质证明:I1,A11(I为单位矩阵) A证明:由于||A||||AI||||A||||I||,等式两边同除以||A||可得I1;又因为
1||I||||A1A||||A1||||A||,不等式1|||A1||||A||两边同除以||A||可得A11。 A
6.8 对于线性方程组 KX1+X2 = 1
X1+KX2+X3 = 2
X2+2X3= 3 K≠0
写出Jacobi迭代公式, K取何值时Jacobi迭代法收敛?说明理由。选择一个合适的参数K,
(0)T
选择初始向量X =(0,0,0),迭代三步。
解: Jacobi迭代公式xBxg为:
0x11 x2kx30由于迭代矩阵B的特征方程为: (21k12010x1k12x2 kkx3302112)0 2kk所以使迭代收敛的k应满足条件:
1121 2kk(0)
因此,当k=2时,迭代收敛。选择初始向量X
=(0,0,0),三步迭代的结果如下:
T
x(1)112021(2)(3)1 x0 x
3213226.9 对于线性方程组
KX1+ X2 = 1
X1+ KX2+ X3 = 2 K≠0 X2+ KX3 = 3
(1) 写出Jacobi迭代法和Seidel迭代法矩阵形式的迭代公式。
(0)(1)(2)
(2)选择一个合适的K的值,使Seidel迭代法收敛,选择初始向量X ,计算出XX(3) X。
解: Jacobi迭代公式xBxg为:
0x11 x2kx30Seidel迭代公式为:
1k1k010x1k12x2 kkx330k00010x11/k001x2/k1/k2 0221/k000x33/k2/kx11/k2 x21/kx30 (1/k1/k2由于Seidel迭代矩阵B的特征方程为:
12)0 2k所以使迭代收敛的k应满足条件: |k|1
因此,当k=2时,迭代收敛。选择初始向量X如下:
(0)
=(0,0,0),三步Seidel迭代的结果
T
x(1)11178322173(2)(3) x x 4165214161
6.10 已知方程组Axb,其中A121,b 0.312(1) 试讨论用雅Jacobi迭代法和Seidel迭代法求解此方程组的收敛性。
(2) 若有迭代公式x(k1)x(k)(Ax(k)b),试确定使该迭代公式收敛的的取值范围。 解:Jacobi迭代式为
02(k)1x(k1)x 00.32IBJ其(BJ)20.320.6
0.61,故Jacobi迭代收敛。
Seidel迭代式为
0x(k1)02(k)1x 0.61.7IBG2(0.6) 00.6其(BG)0.61,故Seidel迭代收敛。 对于以下迭代式
x(k1)x(k)(Ax(k)b)(IA)x(k)b
IA2(10.6)(10.6)
0.311故IA的特征值为1(10.6),1(10.6)。 当5(10.6)0时,有(IA)1,从而迭代收敛。
因篇幅问题不能全部显示,请点此查看更多更全内容
Copyright © 2019- haog.cn 版权所有 赣ICP备2024042798号-2
违法及侵权请联系:TEL:199 1889 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务