|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
function [x,XA]=GaussXQByOrder(A,b)7 n, R* v& \8 }" l
N = size(A);) V+ B8 T) i1 ^2 S
n = N(1);# A, ~, Q9 e. h& ~( F- M E
* W0 |- D4 ?7 i# ?: ~for i=1:(n-1)
. \& D3 P' X1 q; x1 r1 {' z for j=(i+1):n
6 a7 O4 _* Q" B/ R; }; }3 g g( `0 o7 Y2 G2 I if(A(i,i)==0)# P) j) \( H. l5 T5 j$ p
disp('对角元素为0!'); %防止对角元素为05 I3 X7 a# ^' G: h
return;8 C l* u9 y; p+ n( V9 O7 ~
end
: T9 T! j: B$ e8 u; t* x l = A(j,i);
# j6 B/ n. g4 G3 |6 L z2 g& Q m = A(i,i);" N" M7 F* M& a9 M5 I% o
A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程' l6 \" ]' g9 ^$ W
b(j)=b(j)-l*b(i)/m;1 V2 K( \; p9 D$ i2 y
end8 V0 \2 k, ?( _$ ]( u1 `! c
end+ m5 W6 ?" J: C; T* j3 h. H( m5 q
9 L* t' W- W! e `" R6 D0 ]1 Px=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数) C6 |4 b0 f- O7 C) O( i, J9 X: r
XA = A; %消元后的系数矩阵
% h5 A& M$ ?4 M0 [2 z7 p
' l1 e7 \. S! u& u% T2 g9 R
9 T# x# G+ X) g$ T在命令窗口进行运行已知出现 O5 ?; O9 A& W
未定义函数或变量 'SolveUpTriangle'。
4 m' y; x2 n, p( l! o/ |
& ?3 I( m* }. {! q$ t出错 GaussXQByOrder (line 18)4 b1 z1 n1 E0 j' C+ d) O6 @
x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数( [0 C& c Q2 g" ?* h; t
4 g+ m& T1 m( W8 R2 E% ] |
|