|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
function [x,XA]=GaussXQByOrder(A,b)8 {' ]" D# I2 b) q9 R
N = size(A);) Y" F5 T; p) `: _0 k3 H
n = N(1);& o; n8 B) A$ |& ~
% J) R- A- n i" |- x) l6 X
for i=1:(n-1) a1 d1 ?. I i# m2 w1 \0 s$ _( \2 e
for j=(i+1):n N1 i/ [' r$ ~- _: B
if(A(i,i)==0)
, r5 U: [7 b1 }2 l' O. l disp('对角元素为0!'); %防止对角元素为04 X/ t7 Z+ e5 e2 v" {7 ]# C
return;
5 V9 O* T- m* {, }% c( q4 U& } end: u! C! L3 D: _+ B& d' x
l = A(j,i);
/ I: ^1 v, A* E1 Z m = A(i,i);
% q+ f9 U4 A3 r A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程0 d. b& P( k: Y
b(j)=b(j)-l*b(i)/m;
c* f8 w* g4 A end
0 T' i0 W: T0 E( U! iend* K. _2 d& d+ j: [# p$ Z
2 j& B7 y" e7 h7 }) h: wx=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数( {! i0 l! x7 x- g, n8 _' U" D7 X
XA = A; %消元后的系数矩阵
9 N& K1 j- W) Y) `! c J; t& k2 e- B& A: ~6 X7 y
) c+ u$ N, S' w( C在命令窗口进行运行已知出现
0 @9 k% P& f7 e5 a- a: ~* y# q未定义函数或变量 'SolveUpTriangle'。/ Y! Y" e2 O! J+ K
" |% }: Q8 o6 X9 Z5 x7 y/ h
出错 GaussXQByOrder (line 18)
$ M3 @, d3 I `! \" c& L* ~x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数& G, c9 M# [* F( f, P
/ E0 m' v# `% W7 m" l! e& I |
|