|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
function [x,XA]=GaussXQByOrder(A,b)
4 b+ j$ W% e$ @4 }2 x3 Y$ FN = size(A);
' o/ a! L$ `- E3 Qn = N(1);
) C# C6 |! I# W2 b3 o" I; `( i8 {+ i8 {2 Y
for i=1:(n-1)
]4 ^' W5 {( a% b3 S' _& j. X for j=(i+1):n' L" d- e0 S0 A
if(A(i,i)==0)0 g# d1 H m& P) r S# N$ Z
disp('对角元素为0！'); %防止对角元素为0
}0 P0 [0 ]/ X/ X' ?' c return;. K/ d3 _; J) |" a) U4 O
end: s' b& h' ]1 z+ p2 G& ^8 U4 m1 O
l = A(j,i);1 s0 ]+ p) Q" k3 }9 Q ^ m- ^6 K
m = A(i,i);
) ?8 p* O3 g% l$ O. M# C$ G A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程
# ?$ M& S; S& k# V b(j)=b(j)-l*b(i)/m;- J" N m0 U2 _6 P; i: b3 L: _
end! _7 B6 Z/ A0 q
end
$ v+ L. u- s$ a# K6 z: O( V; Y D8 X: {. i+ A
x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数+ [ R+ H/ g: Y
XA = A; %消元后的系数矩阵& x8 l! a' S5 S
, w! y, S; Q( l# [* o) L" M
5 D7 v. \ H" R, i% `" ^2 [在命令窗口进行运行已知出现
3 }0 k m% w, j: |* T未定义函数或变量 'SolveUpTriangle'。# g8 N6 V2 w: x5 u
9 J7 U& p& R+ g/ {. G2 r
出错 GaussXQByOrder (line 18); |# `( s8 L0 |% d% \6 {
x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数
# L6 Q$ G7 z9 h6 c3 P2 N2 w1 _
) G3 K& w" m* L3 r; M4 @- Y) K& O6 c |
|
/1 
发表于 2020-2-27 15:15
收藏
淘帖
支持！
反对！