电力系统算法

zibuyu_2017

function [X]=Gauss(A,B)
  `6 ]4 {3 ]9 L7 h1 X  N[m,n]=size(A);
l=length(B);
if n~=m
( L4 ]7 ?$ d$ w% f2 R1 ~    fprintf('wrong');
    return;end
. l- h+ [6 T2 }5 L' C7 V- D9 yif m~=l
% S+ i) ^8 e& {5 Z# C$ M1 ?    fprintf('wrong');
. F! B* e1 _; \) A$ ]" c    return;end
fprintf('go to calculate');
) V* }! H9 s2 ]# j# a0 z: w* D" bfor k=1:n-1
  X/ `$ W  ?* l1 V" s2 g    max=0;
    for j=k:n
: l# O7 e+ _: L( z9 d+ ^! ^0 z        if abs(A(j,k))>max
4 b9 S' ^' `& x8 V' G% o! f1 u            max=abs(A(j,k));mark=j;
; k/ L- ]- U- p' {- r3 H        end
. |2 F% Z8 E  l/ {. N- [2 F0 K5 B4 c    end
0 Y2 C0 E8 p2 v7 n    if abs(max)<1e-5
. n! B* F' t5 O& U+ P  e3 d+ j       fprintf('wrong');
/ T7 \5 p6 v2 T  l0 P% T       return;end
- f' b5 f* j4 L5 q  K- T       for i=k:n
9 X, ?" V/ T; B% u          a=A(k,i);A(k,i)=A(mark,i);A(mark,i)=a;
4 x" B+ c/ `& u3 g" A       end
       b=B(k);B(k)=B(mark);B(mark)=b;
    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);
8 ^6 l+ t4 l- y9 k1 t; {       end
       B(j)=B(j)-m*B(k);
" C& j( w6 u' K$ A( e   end
5 W2 C  k' S  l' }7 NEnd
7 v/ X5 G7 S- q7 L5 Gfor j=n:-1:1
    for i=j+1:n
: m2 n% U* z. Z( a/ [    B(j)=B(j)-A(j,i)*X(i);
    end
! t2 P/ v1 F7 k1 s( I! E  H    X(j)=B(j)/A(j,j);
end
; z1 b) M5 Z8 P1 r& F" i, e

zibuyu_2017

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