?? sor.m
字號:
%Successive Over Relaxation迭代程序
function [x,sp]=sor(a,b,n,x1,w)
%誤差
e=ones(n,1);
%迭代的解向量
x2=zeros(n,1);
%迭代的次數
k=0;
%當誤差沒有滿足要求時繼續迭代
while norm(e,2)>1e-6
%每隔5步顯示迭代結果
if (rem(k,5)==0)
str=sprintf('X%d= %g %g %g %g\n',k,x1);
disp(str)
end
k=k+1;
%計算x2(i)
for i=1:n
s1=0;
s2=0;
if i-1>0
for j=1:i-1
s1=s1+a(i,j)*x2(j);
end
end
for j=i:n
s2=s2+a(i,j)*x1(j);
end
x2(i)=x1(i)+w*(b(i)-s1-s2)/a(i,i);
end
%計算誤差向量
e=x2-x1;
%把x2作為新的迭代初值
x1=x2;
end
%返回迭代步數
sp=k;
%返回迭代值
x=x2;
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -