|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
function [x,XA]=GaussXQByOrder(A,b)
8 | Z Y0 V( HN = size(A);2 e9 m& p. ]. ]" S! `
n = N(1);
; k' e; l" {' G" x0 {7 Z/ K$ g9 r
: V) f) y# m+ @. D# m7 Y0 E! _for i=1:(n-1) f L( j* N: y! [2 y9 s+ w9 }
for j=(i+1):n
5 I0 A: C9 M% b% q' n& _8 O if(A(i,i)==0)1 C2 \2 H) \. b6 n6 j
disp('对角元素为0!'); %防止对角元素为01 r3 [& Y. a0 B$ E- P r
return;4 i& h6 x& z* v6 j& c7 R( L
end2 q" y& O3 Z9 u1 }4 n
l = A(j,i);$ v0 q; ^2 N/ F+ J. o
m = A(i,i);) t2 D' M. @! v: b1 R6 ?
A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程& i' d, Y7 d0 _- X% O) R0 s
b(j)=b(j)-l*b(i)/m;/ s6 ?4 Y3 }* z' l6 E4 J
end, w6 x' o7 ]7 o- P) N* o
end
" {5 P: Q- \( `! w7 h' w/ R! Y7 p! r3 n- n
x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数
' P6 c4 g" K: C8 b) |* B) M( ^XA = A; %消元后的系数矩阵2 ]3 U7 s( J, ?1 W2 z9 E
6 u3 ]2 t9 D8 Z7 n. F/ ]6 p. Z+ S. c
/ d9 {7 ^% B1 ^9 o G: |; m: V
在命令窗口进行运行已知出现
9 w/ d- o( {) ?, ?, g. n) ^# r未定义函数或变量 'SolveUpTriangle'。) U1 t0 B1 ? C/ U( R- Z2 t3 t
( |: u, d! {( x, t# t7 c出错 GaussXQByOrder (line 18), |8 q% X: A. U6 J: Q$ t
x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数
: ?2 f) J2 y3 |7 [ l4 y/ [5 l: |0 a0 i0 c; x
|
|