?? gs.m
字號:
function GS=GS(A,b,x0,eps)
%該函數文件是通過Gauss-Seidel選代格式求解線性方程組的近似解
%其中,A為方程組的系數矩陣,b為方程組的右端向量,x0為選代初值
if nargin==3
eps=1.0e-6;
elseif nargin<3
error
return
end
D=diag(diag(A)); %求解對角陣
L=-tril(A,-1); %求解出嚴格下對角陣
U=-triu(A,1); %求解出嚴格上對角陣
G=(D-L)\U; %求解出迭代矩陣
f=(D-L)\b;
GS=G*x0+f;
while norm(GS-x0)>=eps
x0=GS;
GS=G*x0+f;
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -