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

乱搞了一通，发现没有错。 

Signed-off-by: facat <facat@ipso.laptop>
This commit is contained in:
facat 2014-04-21 10:33:35 +08:00
parent b906ac299c
commit 0561f915c9
1 changed files with 22 additions and 21 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

43
run.m
View File

@ -2,7 +2,7 @@ clc
clear clear
lineZ=readLineZ('.\feeder13\lineParameter.txt'); lineZ=readLineZ('.\feeder13\lineParameter.txt');
[ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ... [ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ...
phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY]=dataRead(lineZ,'.\feeder13\data1.txt'); phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY]=dataRead(lineZ,'.\feeder13\data2.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;
@ -12,16 +12,19 @@ Tf2p=inv(Tp2f);
fsY1amp=abs(fsY1); fsY1amp=abs(fsY1);
[r,c,fsY1ang]=find(fsY1); [r,c,fsY1ang]=find(fsY1);
fsY1ang=angle(fsY1ang); fsY1ang=angle(fsY1ang);
%%参考电压
refV0=1;
refV1=0;
refV2=0;
%refV1=1;
Pabc=phaseASpotLoadP+phaseBSpotLoadP+phaseCSpotLoadP; Pabc=phaseASpotLoadP+phaseBSpotLoadP+phaseCSpotLoadP;
Qabc=phaseASpotLoadQ+phaseBSpotLoadQ+phaseCSpotLoadQ; Qabc=phaseASpotLoadQ+phaseBSpotLoadQ+phaseCSpotLoadQ;
busNum=length(phaseASpotLoadP); busNum=length(phaseASpotLoadP);
%给序电压赋初值 %给序电压赋初值
Vmf1=sparse(ones(busNum,1)); Vmf1=sparse(ones(busNum,1));
Vmf2=sparse(ones(busNum,1)); Vaf1=sparse(zeros(busNum,1));
Vmf0=sparse(busNum,1); Vf1=Vmf1.*exp(1j*Vaf1);
Vaf1=sparse(busNum,1); %Vaf1=sparse(zeros(busNum,1));
Vaf2=sparse(busNum,1);
Vaf0=sparse(busNum,1);
%先求解正序的 %先求解正序的
PQi=nodeNum; PQi=nodeNum;
PG=sparse(busNum,1); PG=sparse(busNum,1);
@ -38,6 +41,7 @@ fsY00=fsY0;
fsY22=fsY2; fsY22=fsY2;
Vf2=sparse(busNum,1); Vf2=sparse(busNum,1);
If2=sparse(busNum,1); If2=sparse(busNum,1);
% Vf0=sparse(busNum,1);
Vf0=sparse(busNum,1); Vf0=sparse(busNum,1);
If0=sparse(busNum,1); If0=sparse(busNum,1);
%准备序矩阵 %准备序矩阵
@ -49,14 +53,6 @@ fsY0(:,Balance)=0;
fsY0=fsY0+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum); fsY0=fsY0+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum);
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(Balance,busNum, ...
PQi,PG,QG,QGi,PD,QD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量
maxD=max(abs([dP;dQ;]));
jaco=Jacobi(Balance,busNum,QGi,Vmf1,YdotSinVolt,YdotCosVolt,diag_Volt_YdotSin,diag_Volt_YdotCos);%雅克比矩阵
[dV, dVangle]=Solv(busNum,jaco,dP,dQ);%解出修正量
[Vmf1, Vaf1]=Modify(Vmf1,Vaf1,dV,dVangle,1);
fprintf('第 %d 次迭代， 最大不平衡量为 %f\n\n',k,full(maxD));
%%
%转换为三相电压 %转换为三相电压
VoltpABC=Tf2p*conj([ Vf0'; (Vmf1.*exp(1j*Vaf1))'; Vf2']); VoltpABC=Tf2p*conj([ Vf0'; (Vmf1.*exp(1j*Vaf1))'; Vf2']);
VoltpA=conj(VoltpABC(1,:)'); VoltpA=conj(VoltpABC(1,:)');
@ -73,16 +69,21 @@ while(k<=kmax && maxD> EPS)
If2(Balance)=0; If2(Balance)=0;
Vf0=fsY0\If0; Vf0=fsY0\If0;
Vf2=fsY2\If2; Vf2=fsY2\If2;
% Vf2=Vf2+refV2;
% Vf0=Vf0+refV0;
[dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ...
PQi,PG,QG,QGi,PD,QD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量
maxD=max(abs([dP;dQ;]));
jaco=Jacobi(Balance,busNum,QGi,Vmf1,YdotSinVolt,YdotCosVolt,diag_Volt_YdotSin,diag_Volt_YdotCos);%雅克比矩阵
[dV, dVangle]=Solv(busNum,jaco,dP,dQ);%解出修正量
[Vmf1, Vaf1]=Modify(Vmf1,Vaf1,dV,dVangle,1);
fprintf('第 %d 次迭代， 最大不平衡量为 %f\n\n',k,full(maxD));
end end
Vf1=Vmf1.*exp(1j*Vaf1); Vf1=Vmf1.*exp(1j*Vaf1);
%% %%
% V=[1.01*exp(1j*0);1.01*exp(-1j*pi/360*120);1.01*exp(-1j*pi/360*240)]; V=[1.0*exp(1j*0);1.0*exp(-1j*pi/180*120);1.0*exp(1j*pi/180*120)];
% abs(Tp2f*V) abs(Tp2f*V);
% Tf2p*Tp2f*V Tf2p*Tp2f*V;
% Vf0=Vf0+refV0;
% Vf2=Vf2+refV2;
% Vf0=Vf0*exp(-1j*(-1.0471975511966));
% Vf1=Vf1*exp(-1j*(1.0471975511966));
%% %%
(Vf0.*conj(fsY00*Vf0)+Vf1.*conj(fsY11*Vf1)+Vf2.*conj(fsY22*Vf2))*3; (Vf0.*conj(fsY00*Vf0)+Vf1.*conj(fsY11*Vf1)+Vf2.*conj(fsY22*Vf2))*3;
conj(Tf2p*[If0(2);If1(2);If2(2)]).*(Tf2p*[Vf0(2);Vf1(2);Vf2(2)]); conj(Tf2p*[If0(2);If1(2);If2(2)]).*(Tf2p*[Vf0(2);Vf1(2);Vf2(2)]);