存放近似解.由式可以看出,每迭代一步只需计算一次矩阵与向量的乘法.
二.具体的算例和操作步骤
2.1. Gauss-Seidel迭代法的伪代码 1.输入问题的参数A,b 2.分解A为D,L,U. 3.计算迭代方程G,f.
4.开始迭代,随机设定一个初值. 5.以迭代方程更新x的值.
6.如果到达迭代次数,则进入步骤7;否则,回到步骤5. 7.输出x,结束.
2.2.具体的算例验证算法的有效性 求解如下的线性方程组
8x1-3x2+2x3=20
4x1+11x2-x3=33 6x+3x+12x=36
23 1
这个方程的真实解为(3,2,1). 程序运行结果: 情况1:输入GS(A,b) GS(A,b) xhis =
0 0 0 2.5000 2.0909 1.2273 2.9773 2.0289 1.0041 3.0098 1.9968 0.9959 2.9998 1.9997 1.0002 2.9998 2.0001 1.0001 3.0000 2.0000 1.0000