|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
function [x,XA]=GaussXQByOrder(A,b)
$ E% V+ y- p& o$ ON = size(A);+ b8 g; [6 D2 h7 U: k8 y8 _
n = N(1);
; {. r; [$ k; j* ~8 X8 t" [) o
7 ]1 B6 g: H( }0 P+ f* o+ Efor i=1:(n-1)$ Y0 K+ U7 f0 A+ ]
for j=(i+1):n; N. e& e X" r6 R; Q
if(A(i,i)==0)( h3 ?, f! h7 Y. ?- {
disp('对角元素为0!'); %防止对角元素为02 Q @% l+ e' {$ F7 m9 k
return;& L4 C1 W. q: l4 w2 g. j
end
& c5 x8 h! c6 @) p l = A(j,i);
" V1 c4 l9 |: t+ ] m = A(i,i);3 M$ g$ Y/ k* a
A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程5 @" m$ ?" o$ ?% N/ W
b(j)=b(j)-l*b(i)/m;, F; J6 u+ S @ r+ S( t' T- C9 C
end- R4 b. H4 q7 x1 B0 S9 h% _
end
. e- y( o- b/ R4 W' O0 I) Y) G; f# X
9 ^- A7 I6 {/ s; q9 @/ mx=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数8 Y( f! t$ q2 H+ t
XA = A; %消元后的系数矩阵
$ ]/ o( J6 Y0 z9 A
+ O3 \4 k$ [! Y* X3 c$ @) \; A; Y, g; E0 Q& D1 G! P
在命令窗口进行运行已知出现* |6 O# n& V b1 c
未定义函数或变量 'SolveUpTriangle'。
$ \; j1 F, c8 C- j: a, U$ B/ y+ h! O, ^5 J2 l5 s/ B; M" I7 b) k
出错 GaussXQByOrder (line 18)
! _. l9 K9 t* d* F$ Q; c- J9 g* Px=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数8 B0 j" N7 x0 @0 x2 f3 H! n5 \
1 ?8 b M/ c, ~! ^1 W& P* w
|
|