数值分析——方程组求根
1.问题
AX=B
A=[a11,a12,a13;a21,a22,a23;a31,a32,a33]
U=[0,a12,a13;0,0,a23;0,0,0]
D=[a11,0,0;0,a22,0;0,0,a33]
L=[0,0,0;a21,0,0;a31,a32,0]
A=U+D+L
b=[b1;b2;b3]
2.Jocobi方式迭代
2.1方程组表示
-
分离变量
a11x1=-a12x2+a13x3+b1
a22x2=-a21x1-a23x3+b2
a33x3=-a31x1-a32x2+b3 -
系数化为1
x1=(-a12x2+a13x3)/a11+b1/a11
x2=(-a21x1-a23x3)/a22+b2/a22
x3=(-a31x1-a32x2)/a33+b3/a33x=Bx+f(x加粗表示向量)
x=(x1;x2;x3)
B=(0,-a12/a11,-a13/a11;-a21/a22,0,-a23/a22;-a31/a33,-a32/33,0);
f=(b1/a11;b2/a22;b3/a33) -
构造迭代格式
x(k+1)=B x(k)+f
2.2 矩阵表示
- 分离
D