电力系统算法

zibuyu_2017

function [X]=Gauss(A,B)
[m,n]=size(A);
2 y0 K* }3 Q. u# T" H% _* [l=length(B);
if n~=m
: i! q. C' B2 i: B( a  Y    fprintf('wrong');
    return;end
/ e- O2 N# \9 k' [# ~if m~=l
    fprintf('wrong');
    return;end
fprintf('go to calculate');
6 S+ A' s  e. Y& D5 n: j+ @+ Gfor k=1:n-1
    max=0;
) s9 {& Z0 B/ o& q/ |    for j=k:n
2 K% O( b% m8 U7 G! _- m0 B8 q( H        if abs(A(j,k))>max
5 J$ Z/ B) e3 A3 k- o3 j- y1 g            max=abs(A(j,k));mark=j;
) ]# r8 }5 T8 R# s- e6 O( r2 a! Y6 ]        end
    end
* t/ X4 F; I' _9 D+ [# K    if abs(max)<1e-5
% s0 l3 V. p. I- S3 d2 J8 ]6 a! h       fprintf('wrong');
       return;end
       for i=k:n
          a=A(k,i);A(k,i)=A(mark,i);A(mark,i)=a;
! f% Y/ R7 g' Q       end
       b=B(k);B(k)=B(mark);B(mark)=b;
7 X; P% Q! @# Q: }9 v. f* H& P    for j=k+1:n
       m=A(j,k)/A(k,k);
       for i=k:n
           A(j,i)=A(j,i)-m*A(k,i);
       end
       B(j)=B(j)-m*B(k);
   end
/ {3 }* [1 p- [+ jEnd
for j=n:-1:1
    for i=j+1:n
    B(j)=B(j)-A(j,i)*X(i);
    end
5 V) @8 c3 ]4 f0 m5 C    X(j)=B(j)/A(j,j);
end
6 }+ T; J/ \7 V; D

zibuyu_2017

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