|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
function [x,XA]=GaussXQByOrder(A,b)4 K0 {5 G; [- H! b
N = size(A);& v0 E, ]" F: }! F
n = N(1);
( M/ \5 |8 n/ R, s \$ Z
3 [1 z4 h1 `! X7 c1 o$ U2 s8 ffor i=1:(n-1)4 d3 x" C8 y& j3 ]2 j
for j=(i+1):n ?, f/ e2 L" W# p& h
if(A(i,i)==0)
3 ~1 t& n8 [+ a* p disp('对角元素为0!'); %防止对角元素为0, N" I+ {8 {+ y( b9 s: {7 N
return;
8 A4 T. q& `# y; C end
; g# G$ e" J* W. k: h l = A(j,i);
, E+ X9 ?5 t6 g. p7 j! l! g m = A(i,i);
- A) [9 _+ f& V l( J7 U6 ~2 J* O A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程
* j9 a, Y$ Y6 v b(j)=b(j)-l*b(i)/m;
8 C9 a# M9 p- N. o end6 M2 I+ O$ n% n6 ~5 b
end4 Q; r6 t! P5 X& t
, f9 r% w2 o' B7 j! [5 Z) A3 `
x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数
- L$ x) O/ Z" M( S5 r* c9 VXA = A; %消元后的系数矩阵 \- n o. {, ~: h& y, ~8 [+ ~# A
8 D4 O4 T) Y* q3 K, B% [
; O! p/ C; F0 \! u& b在命令窗口进行运行已知出现
. ~) @1 X' U/ p$ [未定义函数或变量 'SolveUpTriangle'。
! M |+ C, Q# P2 L* i' F4 \
0 x$ T$ W, n; Z3 T出错 GaussXQByOrder (line 18)
. D$ G% `' Y9 \ ]' V) Ax=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数; K) w+ V9 B" j% }2 o
* ] I9 Q* j1 B: t z0 y |
|