|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
function [x,XA]=GaussXQByOrder(A,b)
3 Z7 O2 D$ y' K! n7 n% r( XN = size(A);! ]6 l8 A& Q( Y3 t
n = N(1);
$ ?+ D8 I% X, L. q9 W2 b$ }5 F, a" Q4 L+ }. |
for i=1:(n-1)
\2 D, ^4 v3 i/ g3 j for j=(i+1):n
6 E0 q+ g& s) U( z3 ~# Y( o if(A(i,i)==0)! U0 @3 J. Z9 T- q0 }
disp('对角元素为0!'); %防止对角元素为0
3 [" ]' }, m f3 A& v return;8 G0 F; w0 G7 D9 `- c, I
end
$ x! r: x# R, | l = A(j,i);
1 h8 g+ k \% y' h7 ^7 D, N m = A(i,i);
" T& d% l6 J I v0 |" O A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程
/ Q& }! o7 S( e. I$ _; C b(j)=b(j)-l*b(i)/m;
* s3 H8 P. q0 E; E end9 ^% p$ ~. L& q O, [; z- `; u1 U
end
7 G. k) w8 u9 f- @4 o9 h4 `7 F/ l9 O3 h
x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数
; g( u% Z& K8 j# u% B, J: HXA = A; %消元后的系数矩阵: l8 @/ R. n5 i+ J7 h! `) c- |) V
( n( ~3 b4 u3 o8 n/ x C( [: f
/ s0 O* ~* a) N! J/ C' m在命令窗口进行运行已知出现
6 b1 W5 q4 g: @- m; d# o未定义函数或变量 'SolveUpTriangle'。5 z- A! u9 b' A+ a p
7 N" m5 X# ], T3 [1 _/ c* V6 [: D
出错 GaussXQByOrder (line 18)* x2 O$ n N9 x& Q4 c$ _
x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数# K7 F* G( o+ J; K% ^
* ?7 V5 a3 H( ~
|
|