交流社区 › 电力系统与电网技术 › 电力系统算法

zibuyu_2017 发表于 2017-8-7 08:40:42

电力系统算法

function =Gauss(A,B)
=size(A);
l=length(B);
if n~=m
    fprintf('wrong');
    return;end
if m~=l
    fprintf('wrong');
    return;end
fprintf('go to calculate');
for k=1:n-1
    max=0;
    for j=k:n
      if abs(A(j,k))>max
            max=abs(A(j,k));mark=j;
      end
    end
    if abs(max)<1e-5
       fprintf('wrong');
       return;end
       for i=k:n
          a=A(k,i);A(k,i)=A(mark,i);A(mark,i)=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);
       end
       B(j)=B(j)-m*B(k);
   end
End
for j=n:-1:1
    for i=j+1:n
    B(j)=B(j)-A(j,i)*X(i);
    end
    X(j)=B(j)/A(j,j);
end

zibuyu_2017 发表于 2017-8-7 08:42:58

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

查看完整版本: 电力系统算法