|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
function [x,XA]=GaussXQByOrder(A,b)
! |0 c" R0 B3 f5 Z! vN = size(A);
: a+ Y/ ?0 ]9 Mn = N(1);
: g: o+ g0 R/ i1 h- w+ Q5 h
, n- w; g! b; ~7 v; gfor i=1:(n-1)4 w" \8 \2 T+ H: v& j
for j=(i+1):n5 \/ e# i9 e0 k
if(A(i,i)==0)+ D R' ~9 \, G2 L# x1 U. |/ c
disp('对角元素为0!'); %防止对角元素为0) Q+ c$ ]$ J8 J' v- t- A
return;
6 q+ {+ G8 p+ E( p/ X7 ~' t+ } end: _% }$ x5 u% H2 M6 f* s
l = A(j,i);; r* C; H4 D7 S. R' v! o7 V
m = A(i,i);; L+ @5 k) I4 B7 c+ ~
A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程3 g" N1 D* P! H6 Y, _
b(j)=b(j)-l*b(i)/m;* v, L1 |* p5 e+ A& O5 w& Y- R
end0 k3 b" r4 S9 L4 E6 u- Y$ V9 E* w
end
8 q$ T' Y5 W7 g2 R% S
/ K' F1 L1 g! @x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数
+ `! ^7 t2 e5 d# q( I+ EXA = A; %消元后的系数矩阵+ k1 J8 [! H# _4 A# M( }
4 ` b2 V& p; d" M Z" V9 U, ? b1 R* y! f9 |& b5 C
在命令窗口进行运行已知出现* x# Q) @3 P3 X I' a/ k
未定义函数或变量 'SolveUpTriangle'。
) e" ~; j4 ?' P& Z6 [% \! t, Q- e, k
出错 GaussXQByOrder (line 18)2 ^' _; _3 M8 k2 @
x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数
4 H% G! k" A7 X4 Y( P
1 P" z: w+ G# t+ c5 D5 h |
|