facat
/
distributionnetwork-power2current
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

1.重新排列了代码 

2.考虑了数据文件中没有电容器项的情况

Signed-off-by: facat <facat@ipso.laptop>
This commit is contained in:
facat 2014-04-23 11:33:34 +08:00
parent 356a3d0fca
commit 276db0e00e
2 changed files with 28 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
dataRead.m
View File

@ -45,6 +45,7 @@ phaseBSpotLoadQ( ismember(setIJ,spotloads(:,1) ) )=spotloads(:,5);
phaseCSpotLoadQ( ismember(setIJ,spotloads(:,1) ) )=spotloads(:,7); phaseCSpotLoadQ( ismember(setIJ,spotloads(:,1) ) )=spotloads(:,7);
%% 补偿电容 %% 补偿电容
cap=data(zeroEntries(2)+1:zeroEntries(3)-1,:); cap=data(zeroEntries(2)+1:zeroEntries(3)-1,:);
if ~isempty(cap)
capNode=nodeNum(cap(:,1)==setIJ); capNode=nodeNum(cap(:,1)==setIJ);
%数据文件中的补偿电容给的是补偿容量，要将其转换为导纳。 %数据文件中的补偿电容给的是补偿容量，要将其转换为导纳。
capB=cap(:,2:4)/1000; capB=cap(:,2:4)/1000;
@ -58,5 +59,10 @@ offSet=(capNode-1)*3+1;
phaseABCY=phaseABCY+sparse(offSet,offSet,1j*capB(:,2),busNum*3,busNum*3); phaseABCY=phaseABCY+sparse(offSet,offSet,1j*capB(:,2),busNum*3,busNum*3);
offSet=(capNode-1)*3+1; offSet=(capNode-1)*3+1;
phaseABCY=phaseABCY+sparse(offSet,offSet,1j*capB(:,3),busNum*3,busNum*3); phaseABCY=phaseABCY+sparse(offSet,offSet,1j*capB(:,3),busNum*3,busNum*3);
else
cap=struct();
cap.capNode=[];
cap.capB=sparse(1,3);
end
end end

42
run.m
View File

@ -1,9 +1,9 @@
clc clc
clear clear
lineZ=readLineZ('.\feeder13\lineParameter.txt'); lineZ=readLineZ('D:\share\feeder123\lineParameter.txt');
[ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ... [ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ...
phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY ... phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY ...
cap]=dataRead(lineZ,'.\feeder13\data1.txt'); cap]=dataRead(lineZ,'D:\share\feeder123\data.txt');
a=exp(1j*2*pi/3); a=exp(1j*2*pi/3);
Tp2f=1/3*[1 1 1; Tp2f=1/3*[1 1 1;
1 a a^2; 1 a a^2;
@ -55,18 +55,18 @@ VoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi);
VoltpC=sparse(ones(busNum,1)).*exp(1j*+120/180*pi); VoltpC=sparse(ones(busNum,1)).*exp(1j*+120/180*pi);
while(k<=kmax && maxD> EPS) while(k<=kmax && maxD> EPS)
k=k+1; k=k+1;
%把补偿电容看作负荷 %把补偿电容看作负荷
SA=VoltpA.*conj(VoltpA.*sparse(cap.capNode,1,1j*cap.capB(:,1),busNum,1)); SA=VoltpA.*conj(VoltpA.*sparse(cap.capNode,1,1j*cap.capB(:,1),busNum,1));
SB=VoltpB.*conj(VoltpA.*sparse(cap.capNode,1,1j*cap.capB(:,2),busNum,1)); SB=VoltpB.*conj(VoltpA.*sparse(cap.capNode,1,1j*cap.capB(:,2),busNum,1));
SC=VoltpC.*conj(VoltpA.*sparse(cap.capNode,1,1j*cap.capB(:,3),busNum,1)); SC=VoltpC.*conj(VoltpA.*sparse(cap.capNode,1,1j*cap.capB(:,3),busNum,1));
iterPD=PD+real(SA+SB+SC)/3; iterPD=PD+real(SA+SB+SC)/3;
iterQD=QD+imag(SA+SB+SC)/3; iterQD=QD+imag(SA+SB+SC)/3;
iterPhaseASpotLoadP=phaseASpotLoadP+real(SA); iterPhaseASpotLoadP=phaseASpotLoadP+real(SA);
iterPhaseBSpotLoadP=phaseBSpotLoadP+real(SB); iterPhaseBSpotLoadP=phaseBSpotLoadP+real(SB);
iterPhaseCSpotLoadP=phaseCSpotLoadP+real(SC); iterPhaseCSpotLoadP=phaseCSpotLoadP+real(SC);
iterPhaseASpotLoadQ=phaseASpotLoadQ+imag(SA); iterPhaseASpotLoadQ=phaseASpotLoadQ+imag(SA);
iterPhaseBSpotLoadQ=phaseBSpotLoadQ+imag(SB); iterPhaseBSpotLoadQ=phaseBSpotLoadQ+imag(SB);
iterPhaseCSpotLoadQ=phaseCSpotLoadQ+imag(SC); iterPhaseCSpotLoadQ=phaseCSpotLoadQ+imag(SC);
[dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ... [dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ...
PQi,PG,QG,QGi,iterPD,iterQD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量 PQi,PG,QG,QGi,iterPD,iterQD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量
maxD=max(abs([dP;dQ;])); maxD=max(abs([dP;dQ;]));
@ -154,14 +154,14 @@ maxD=10000;
tic tic
while(k<=kmax && maxD> EPS) while(k<=kmax && maxD> EPS)
k=k+1; k=k+1;
[dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance3P,busNum*3, ... [dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance3P,busNum*3, ...
PQi3P,PG,QG,QGi3P,PD3P,QD3P,Vp3m,Vp3a,Y,Yangle,r,c,0,0,0,0); PQi3P,PG,QG,QGi3P,PD3P,QD3P,Vp3m,Vp3a,Y,Yangle,r,c,0,0,0,0);
maxD=max(abs([dP;dQ;])); maxD=max(abs([dP;dQ;]));
jaco=Jacobi(Balance3P,busNum*3,QGi3P,Vp3m,YdotSinVolt,YdotCosVolt,diag_Volt_YdotSin,diag_Volt_YdotCos);%雅克比矩阵 jaco=Jacobi(Balance3P,busNum*3,QGi3P,Vp3m,YdotSinVolt,YdotCosVolt,diag_Volt_YdotSin,diag_Volt_YdotCos);%雅克比矩阵
[dV, dVangle]=Solv(busNum*3,jaco,dP,dQ);%解出修正量 [dV, dVangle]=Solv(busNum*3,jaco,dP,dQ);%解出修正量
[Vp3m, Vp3a]=Modify(Vp3m,Vp3a,dV,dVangle,1); [Vp3m, Vp3a]=Modify(Vp3m,Vp3a,dV,dVangle,1);
fprintf('第 %d 次迭代， 最大不平衡量为 %f\n',k,full(maxD)); fprintf('第 %d 次迭代， 最大不平衡量为 %f\n',k,full(maxD));
fprintf('迭代时间%f\n',toc); fprintf('迭代时间%f\n',toc);
end end
NewtonToc=toc; NewtonToc=toc;
fprintf('牛顿法计算时间 %f\n',NewtonToc); fprintf('牛顿法计算时间 %f\n',NewtonToc);