power_flow.m

蟻拓演算法應用在電力系統的無效功率潮流分析

%   Power flow solution by Newton-Raphson method
%clear;
%data30;
%Y_bus;
basemva=100;   accuracy=0.001;   maxiter=20;   flag=0;
nbs=length(busdata(:,1));  ngr=length(find(busdata(:,2) == 2));  nld=nbs-ngr-1;
Vm=busdata(:,3);   Theta=busdata(:, 4)*pi/180;
P0=(busdata(:,7)-busdata(:,5))/basemva;
Q0=(busdata(:,8)-busdata(:,6)+busdata(:,11))/basemva;
Gy=real(Ybus);     By=imag(Ybus);
XX=[Theta(2:nbs);   Vm(ngr+2:nbs)];     % (10.31)  size: (2*nbs-ngr-2) x 1
nl =(linedata(:,1));     nr = linedata(:,2);        R = linedata(:,3);
X = linedata(:,4);     Bc = j*(linedata(:,5))/2;  a = linedata(:,6);
%Index = find( a<= 0);   a(Index) = 1;
for  iter=1:maxiter
     THIJ = Theta*ones(1,nbs)-ones(nbs,1)*Theta.';   % 璸衡à