电力系统算法

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

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

您需要登录才可以下载或查看，没有账号？立即加入

×

function [X]=Gauss(A,B)
[m,n]=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

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

账号		自动登录	找回密码
密码			立即加入

电力系统算法

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

相关帖子