本文通过一个具体的数学例子，展示了如何使用赛德尔迭代法求解线性方程组。其中包括了迭代过程的具体实现代码，以及如何设定迭代精度来获得满意的解。
%计算方法书上p62 例3.5
%赛德尔迭代
%程序：张建波
%2005-11-18 PM23:00
a=[10 -2 -1;-2 10 -1;-1 -2 5;-1 -2 5];
b=[3 15 10];
x=[0;0;0];
N=3;%N 是解得个数，X1,X2,X3
es=0.0001; %精度控制
while(0<1)
es0=0;
%开始迭代
for i=1:N
s=x(i);
x(i)=0;
x(i)=(b(i)-a(i,:)*x)/a(i,i);
es0=max(abs(s-x(i)),es0);
end
x
if es0<es break;
end
end
%赛德尔迭代
%程序：张建波
%2005-11-18 PM23:00
a=[10 -2 -1;-2 10 -1;-1 -2 5;-1 -2 5];
b=[3 15 10];
x=[0;0;0];
N=3;%N 是解得个数，X1,X2,X3
es=0.0001; %精度控制
while(0<1)
es0=0;
%开始迭代
for i=1:N
s=x(i);
x(i)=0;
x(i)=(b(i)-a(i,:)*x)/a(i,i);
es0=max(abs(s-x(i)),es0);
end
x
if es0<es break;
end
end
扫一扫
1233
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
