|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
function [x,XA]=GaussXQByOrder(A,b)1 K7 Z& |# Q4 C* g
N = size(A);
. R% X% |- Y' G7 g- S5 ?1 u; Ln = N(1);
% `# G f5 V2 M: V9 }% [* G
% h# p( d- i1 A: A- N. G Ifor i=1:(n-1)
1 }0 _8 B9 [" E1 a$ `5 d" w for j=(i+1):n
2 i& H* ?% e. E if(A(i,i)==0)
6 S$ A2 j9 V! A$ I disp('对角元素为0!'); %防止对角元素为09 L0 B v: Z4 P5 ?4 `2 `
return;
' X3 I! k: |6 {* g end
6 D3 e8 @' ]/ `3 f% Q7 j- P i l = A(j,i);2 f& q1 A" o1 X V. ?0 ]$ j1 l
m = A(i,i);& K1 a1 x( H! ^: [/ q
A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程
3 q# X& J, u- |$ }# g& ~/ g b(j)=b(j)-l*b(i)/m;
% A+ H- s. R6 \, q# x0 `: X6 u6 l end6 `7 r, F9 X8 M' h
end; L5 d0 m" R _
5 i9 l# B e; T' A% F* B) s
x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数2 N/ A# j+ c$ m
XA = A; %消元后的系数矩阵
, s7 E h# }) C6 f6 ^6 E/ G: q8 S# O0 w4 d& x( ^
% W7 Z o. j0 E3 {* m* F# c9 l
在命令窗口进行运行已知出现
/ E, u" d% i5 i' W, L# L( t未定义函数或变量 'SolveUpTriangle'。' K& m* ~! D4 ?# [+ N: X7 ~. {
1 R+ P- b& T4 e& B* B% r: r
出错 GaussXQByOrder (line 18)
, I) v6 }( m0 i8 }" _x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数
o! I6 [$ j. c, B, Q; ?$ S
/ B& s8 i i Q/ @$ a* _& m |
|