2014年10月13日
11:25
求解例1
Clear
clc
i=1;
a=[8 -3 2;4 11 -1;6 3 12];%输入系数矩阵
d=diag(diag(a));%求对角矩阵D
l=d-tril(a);%求下三角矩阵L
u=d-triu(a);%求上三角矩阵U
b=[20;33;36];%输入b
x0=zeros(3,1);%输入迭代初始值
G=inv(d-l)*u;%迭代矩阵
f=inv(d-l)*b;%求f
x=G*x0+f;%迭代方程
while norm(x-x0,inf)>=1e-6%设置精度（迭代停止条件）
x0=x;
x=G*x0+f;
i=i+1;
end
x%输出x
i%输出迭代次数
11:25
求解例1
Clear
clc
i=1;
a=[8 -3 2;4 11 -1;6 3 12];%输入系数矩阵
d=diag(diag(a));%求对角矩阵D
l=d-tril(a);%求下三角矩阵L
u=d-triu(a);%求上三角矩阵U
b=[20;33;36];%输入b
x0=zeros(3,1);%输入迭代初始值
G=inv(d-l)*u;%迭代矩阵
f=inv(d-l)*b;%求f
x=G*x0+f;%迭代方程
while norm(x-x0,inf)>=1e-6%设置精度（迭代停止条件）
x0=x;
x=G*x0+f;
i=i+1;
end
x%输出x
i%输出迭代次数
