|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
function [x,XA]=GaussXQByOrder(A,b)
( ~; O6 x" j6 h, r; Z! nN = size(A);
' d& K( r) n6 A( gn = N(1);
) Z7 e! d3 d* z! n( H/ j, Y9 A$ e0 m+ y* D5 h6 B
for i=1:(n-1)
: J% V, E. ?! D6 |1 p6 l ~9 O for j=(i+1):n
/ }& }8 W. u( t( O' L3 A* y if(A(i,i)==0)
" \4 [: e; [3 e7 k& c; p9 B$ b0 i disp('对角元素为0!'); %防止对角元素为01 n! _7 X4 j! w
return;3 ]+ O! ?7 O( O3 [, I M
end" l! M9 y2 _& _1 u3 G
l = A(j,i);! A$ B( E7 h ~& j; w K
m = A(i,i);
# ~3 a- K) t3 G$ S6 b$ Y A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程
$ s: V2 \* [" J8 O, w! @- r b(j)=b(j)-l*b(i)/m;
" l$ E' O: A; }# k, _0 N end
v; _ D& }/ ~2 @% nend
3 K# i1 ^3 U( c5 k* d4 q
/ y$ a1 y, l v5 U1 ?* I# kx=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数8 t$ q4 W' g' G' i
XA = A; %消元后的系数矩阵
. H; v4 E8 ~) ?/ c( T; y
# ` ^; f5 }: e, y# y$ Y$ G
0 b- v: N# u9 {; l9 ~" O在命令窗口进行运行已知出现3 {% z8 Q! z x1 H9 h4 _) `7 p
未定义函数或变量 'SolveUpTriangle'。
" H2 d4 v0 S/ b# g M: |3 z5 m _9 V8 a4 q, ~% V# A* D8 }4 i
出错 GaussXQByOrder (line 18)
m' d! U2 s5 ax=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数- V# W. k& K9 W! Q
( c2 `5 \3 E3 N x: D |
|