?? sor.m
字號:
function [k,w,x]=ssor(a,b,w,x0,e,N)
%用逐次超松弛法解線性方程組Ax=b;
%格式[k,w,x]=ssor(a,b,w,x0,e,N),a為系數矩陣,b為右端向量,x返回解向量,k迭代次數,w松弛因子,x0初始向量,e容取誤差,N迭代次數上限,;
n=length(b);
x=x0;x0=x+2*e;
k=0;
while max(abs(x0-x))>e&k<N
k=k+1;x0=x;
for i=1:n
s=0;m=0;
for j=1:i-1
s=s+a(i,j)*x(j);
end
for j=i+1:n
m=m+a(i,j)*x0(j);
end
x(i)=(1-w)*x0(i)+w/a(i,i)*(b(i)-s-m);
end
end
if k==N error('已超過迭代次數');end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -