矢量编程潮流及最优潮流算法.zip_optimal power flow_power system_最优潮流_潮流计算

clc;clear all [FileName,PathName]=uigetfile('*.dat'); tic; path=strcat(PathName,FileName); data=textread(path); [nodeNum,branchNum,balanceNode,base,Kmax,Precision,LinePara,GroundPara,TransPara,PowerPara,PvPara,P,Q,PVV,PVI]=ReadData(data); [Y]=FormY(nodeNum,LinePara,GroundPara,TransPara); V=ones(nodeNum,1); %[I,J]=find(V); %V=sparse(I,J,V(I),nodeNum,1); format long V=V+PVV; [I,J]=find(V); V=sparse(I,J,V(I),nodeNum,1); d=zeros(nodeNum,1); Y1=abs(Y); disp(' 迭代次数 最大不平衡量 ') for k=1:Kmax di=repmat(d,1,nodeNum); dj=repmat(d',nodeNum,1); aij=angle(Y); dij=di-dj-aij; [DP,DQ]=culation(V,Y1,dij,P,Q,PVI,balanceNode); DPQ=[DP;DQ]; s=abs(DPQ); if max(s)<Precision break; end %keyboard JJ=jacobian(V,Y1,dij,nodeNum,PVI,balanceNode); DdDV=-JJ\DPQ; d=d+DdDV(1:nodeNum); DV=(DdDV(nodeNum+1:2*nodeNum)); V=(V)+inv(diag(V))*(DV); disp([k-1,max(s)]) k=k+1; end toc k; V;