|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
function [x,XA]=GaussXQByOrder(A,b) q' X: S% e2 f# |: G& ]! Z1 t# Z* E
N = size(A);
2 I. y7 N' C6 L% }" Gn = N(1);
+ a+ D# ]* G* ]3 A/ |
$ L3 Y" P4 w. y6 u( M8 j" kfor i=1:(n-1)
4 }: ]0 B7 P* H- D7 l for j=(i+1):n
' y1 q: s, H/ B Y. k if(A(i,i)==0). X8 J3 i) x/ U3 v
disp('对角元素为0!'); %防止对角元素为0
& a( B8 t Z1 H' S7 y return;
3 M+ T7 g3 J+ a! f. T end, ?* }% S) N5 F$ _
l = A(j,i);
: y6 F+ v! D/ Z- C+ G: _ m = A(i,i);0 K- U5 l3 p/ J' d! Q& Y
A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程, f- p' U; F# K- a( S* O
b(j)=b(j)-l*b(i)/m;6 W, F* C f4 u1 s6 n) f, n
end9 d! M; U# c# J/ m' K1 M
end
' O7 d: n# P" w; ?
2 U x" t$ F9 r4 ix=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数
# G2 ?1 K) k6 [, p8 l# uXA = A; %消元后的系数矩阵
' R$ x7 _+ a% _6 Y R$ I5 ^, K9 h8 U5 Q! r1 ]. ^+ \
K5 i$ H( b) X" K/ V) l在命令窗口进行运行已知出现# j4 ^7 E1 }; h) Y) e) z; e
未定义函数或变量 'SolveUpTriangle'。
6 c. j c0 Q6 t* f/ U( `
) o$ o7 f& E" O出错 GaussXQByOrder (line 18)
6 E) ]' V/ O3 P+ q: B0 gx=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数: B7 c) @9 t4 Z# f" `" k
8 R9 ?& o, d4 \3 `# P, |, |5 p
|
|