|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
function [x,XA]=GaussXQByOrder(A,b)6 {# r [7 {* r
N = size(A);
4 ^) y8 F) ~: Jn = N(1); d( i* n2 \9 f/ A9 E, q
+ H3 z' O7 L% ?' y, \ U' cfor i=1:(n-1)
- ^) w7 P @ b) G$ z' R for j=(i+1):n; D8 h/ O1 H. W1 G7 |- J
if(A(i,i)==0)
, S( e! L. y2 F+ j# r0 S/ B* {& X disp('对角元素为0!'); %防止对角元素为0( G- I' @. q5 Y6 Z- q
return;
8 l) E/ r2 `2 q* B end% i4 c3 |0 l; x+ f% d
l = A(j,i);
5 g0 m% x: Z5 q/ \. V# V: f m = A(i,i);' A& n1 X" Q5 M3 |! I5 [
A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程
. z3 O z0 }7 k) |) M5 y' q- ^' L b(j)=b(j)-l*b(i)/m;4 p- m) r7 O+ F: y& t. t, I
end
6 p5 s; h+ n1 U m- S$ rend
/ b% h3 u: g- s1 d# ]' Y
5 ?. N: X2 U+ M$ V0 W/ \2 {x=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数$ l d3 Q1 p1 S) |+ b( n
XA = A; %消元后的系数矩阵
( V5 n& A, \ J, |$ h' n8 U4 V' M' M
7 T; ^9 G, O/ d- M
在命令窗口进行运行已知出现
, H6 m9 }! X% a+ u未定义函数或变量 'SolveUpTriangle'。 D7 ?/ \* R0 K6 s/ X/ S9 q
; M) `- d3 C$ d; k' t: B出错 GaussXQByOrder (line 18)
4 k; L3 J) P! v" N1 Rx=SolveUpTriangle(A,b); %通用的求上三角系数矩阵线性方程组的函数' _% r) h" F+ F; J8 n
& J3 S4 G6 o, W! d2 w6 E% y v |
|