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

先不考虑电容 

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk 2014-10-29 21:56:38 +08:00
parent 7409723aae
commit 4c2f16c8c7
1 changed files with 12 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
run.m
View File

@ -65,6 +65,10 @@ while(k<=kmax && maxD> EPS)
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));
%先不要电容
SA=0;
SB=0;
SC=0;
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);
@ -231,7 +235,6 @@ mIf2=-If2;
% fsY11(Balance,:)=0; % fsY11(Balance,:)=0;
% fsY11=fsY11+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum); % fsY11=fsY11+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum);
% mIf1(3)=1; % mIf1(3)=1;
%% 先算正序的 %% 先算正序的
%平衡节点电流 %平衡节点电流
fprintf('正序\n'); fprintf('正序\n');
@ -252,10 +255,8 @@ sigmaI1=normrnd(0,sigma,length(Loadi),1);
I1measurement=mIf1(Loadi).*(1+sigmaI1);%测量值是等效发电机电流 I1measurement=mIf1(Loadi).*(1+sigmaI1);%测量值是等效发电机电流
wI1r=abs( real(mIf1(Loadi).*sigmaI1) ); wI1r=abs( real(mIf1(Loadi).*sigmaI1) );
wI1i=abs( imag(mIf1(Loadi).*sigmaI1) ); wI1i=abs( imag(mIf1(Loadi).*sigmaI1) );
[ V1r,V1i,I1r,I1i ]=IPMLoop(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,BalI1r,BalI1i,busNum,Loadi,fsY11,Balance,1 ); [ V1r,V1i,I1r,I1i ]=IPMLoop(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,BalI1r,BalI1i,busNum,Loadi,fsY11,Balance,1 );
% f=sum(([real(Imeasurement);imag(Imeasurement)]-[-I1r;-I1i]).^2)+sum((real(Vmeasurement)-V1r(Loadi)).^2)+sum((imag(Vmeasurement)-V1i(Loadi)).^2); f=sum(([real(I1measurement);imag(I1measurement)]-[I1r;I1i]).^2)+sum((real(rV1)-V1r).^2)+sum((imag(rV1)-V1i).^2);
% f=sum(([real(-mIf1(Loadi));imag(-mIf1(Loadi))]-[-I1r;-I1i]).^2)+sum((real(rV1(Loadi))-V1r(Loadi)).^2)+sum((imag(rV1(Loadi))-V1i(Loadi)).^2);
% fprintf('目标值 %f\n',full(f)); % fprintf('目标值 %f\n',full(f));
%% 算负序的 %% 算负序的
fprintf('负序\n'); fprintf('负序\n');
@ -277,7 +278,7 @@ I2measurement=mIf2(Loadi).*(1+sigmaI2);%
wI2r=abs( real(mIf2(Loadi).*sigmaI2) ); wI2r=abs( real(mIf2(Loadi).*sigmaI2) );
wI2i=abs( imag(mIf2(Loadi).*sigmaI2) ); wI2i=abs( imag(mIf2(Loadi).*sigmaI2) );
[ V2r,V2i,I2r,I2i ]=IPMLoop(V2measurement,wV2r,wV2i,I2measurement,wI2r,wI2i,BalI2r,BalI2i,busNum,Loadi,fsY22,Balance,0 ); [ V2r,V2i,I2r,I2i ]=IPMLoop(V2measurement,wV2r,wV2i,I2measurement,wI2r,wI2i,BalI2r,BalI2i,busNum,Loadi,fsY22,Balance,0 );
% f=sum(([real(Imeasurement);imag(Imeasurement)]-[-I1r;-I1i]).^2)+sum((real(Vmeasurement)-V2r(Loadi)).^2)+sum((imag(Vmeasurement)-V2i(Loadi)).^2); f=sum(([real(I2measurement);imag(I2measurement)]-[I2r;I2i]).^2)+sum((real(rV2)-V2r).^2)+sum((imag(rV2)-V2i).^2);
% fprintf('目标值 %f\n',full(f)); % fprintf('目标值 %f\n',full(f));
%% 算零序 %% 算零序
fprintf('零序\n'); fprintf('零序\n');
@ -299,9 +300,9 @@ I0measurement=mIf0(Loadi).*(1+sigmaI0);%
wI0r=abs( real(mIf0(Loadi).*sigmaI0) ); wI0r=abs( real(mIf0(Loadi).*sigmaI0) );
wI0i=abs( imag(mIf0(Loadi).*sigmaI0) ); wI0i=abs( imag(mIf0(Loadi).*sigmaI0) );
[ V0r,V0i,I0r,I0i ]=IPMLoop(V0measurement,wV0r,wV0i,I0measurement,wI0r,wI0i,BalI0r,BalI0i,busNum,Loadi,fsY00,Balance,0 ); [ V0r,V0i,I0r,I0i ]=IPMLoop(V0measurement,wV0r,wV0i,I0measurement,wI0r,wI0i,BalI0r,BalI0i,busNum,Loadi,fsY00,Balance,0 );
%% f=sum(([real(I0measurement);imag(I0measurement)]-[I0r;I0i]).^2)+sum((real(rV0)-V0r).^2)+sum((imag(rV0)-V0i).^2);
% f=sum(([real(Imeasurement);imag(Imeasurement)]-[-I1r;-I1i]).^2)+sum((real(Vmeasurement)-V0r(Loadi)).^2)+sum((imag(Vmeasurement)-V0i(Loadi)).^2);
% fprintf('目标值 %f\n',full(f)); % fprintf('目标值 %f\n',full(f));
%状态量 %状态量
% SEVoltpA=sparse(ones(busNum,1)); % SEVoltpA=sparse(ones(busNum,1));
% SEVoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi); % SEVoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi);
@ -325,8 +326,11 @@ SEVf0=V0r+1j*V0i;
SEVf1=V1r+1j*V1i; SEVf1=V1r+1j*V1i;
SEVf2=V2r+1j*V2i; SEVf2=V2r+1j*V2i;
SEVoltpABC=Tf2p*conj([ SEVf0'; SEVf1'; SEVf2']); SEVoltpABC=Tf2p*conj([ SEVf0'; SEVf1'; SEVf2']);
% SEVoltpABC2=Tf2p*conj([ rV0'; rV1'; rV2']);
% 三序电流 % 三序电流
SEIf0=I0r+1j*I0r; SEIf0=I0r+1j*I0r;
SEIf1=I1r+1j*I1r; SEIf1=I1r+1j*I1r;
SEIf2=I2r+1j*I2r; SEIf2=I2r+1j*I2r;
SEIpABC=Tf2p*conj([SEIf0';SEIf1';SEIf2']); SEIpABC=Tf2p*conj([SEIf0';SEIf1';SEIf2']);
%看一下差多少
(full(abs(VoltpABC))-abs(SEVoltpABC))./abs(SEVoltpABC)*100