电力系统算法

zibuyu_2017

function [X]=Gauss(A,B)
; N4 @1 U0 h* J$ G- y' I( @/ |[m,n]=size(A);
1 k' v& u! v" b7 o3 Dl=length(B);
) n3 b) n0 R8 I- o' Mif n~=m
9 z) b% X$ k% J* J1 ~$ H. u8 K    fprintf('wrong');
    return;end
- l  ?0 ~# K/ Q+ Wif m~=l
7 k1 a8 C$ E6 K& k    fprintf('wrong');
$ ^9 n5 Y1 b' w- l+ ]" a& I    return;end
( K* l/ W& T' |; B5 n* X; ufprintf('go to calculate');
for k=1:n-1
* S5 r0 `* S! U! b7 l% D    max=0;
( o* u6 g+ S) N) Q    for j=k:n
        if abs(A(j,k))>max
            max=abs(A(j,k));mark=j;
. a) \# [9 p$ f& `/ p        end
4 [% {0 |9 D- n1 r: ]# U: p    end
    if abs(max)<1e-5
7 ?/ e; o4 B9 x4 q% T$ d4 `: d  ?       fprintf('wrong');
       return;end
# w) X  f' k, m2 Q- t1 f       for i=k:n
          a=A(k,i);A(k,i)=A(mark,i);A(mark,i)=a;
0 ^9 l, v: o( c' ~       end
! X* M$ L0 V6 J       b=B(k);B(k)=B(mark);B(mark)=b;
    for j=k+1:n
& Z* d$ H! K3 g4 p       m=A(j,k)/A(k,k);
2 g4 M; V. C! Z. B       for i=k:n
" a$ [1 h, Z; V1 S           A(j,i)=A(j,i)-m*A(k,i);
- @. V" `1 C% J8 M       end
. [* m/ f0 l8 N' E* ^) O  Q       B(j)=B(j)-m*B(k);
7 K; a, h5 ^$ \' k- P   end
. u# V* u: F2 S3 [End
& y. V( t5 m' ~0 M. _' s  `for j=n:-1:1
0 P: G+ X6 M' k9 O+ P* x1 L    for i=j+1:n
    B(j)=B(j)-A(j,i)*X(i);
    end
2 Q' H# I2 y6 A# G: }3 Z7 ?    X(j)=B(j)/A(j,j);
end
, _8 @: L6 R5 ]+ Q$ S/ J+ O/ g
6 P8 j. u$ o! T3 i

zibuyu_2017

这个帖子不错，大家快来顶起来！