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-11-10 22:33:39 +08:00
parent d411972751
commit 3db4877bb7
4 changed files with 14 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
IPMLoop.m
View File

@ -39,10 +39,10 @@ while(abs(Gap)>0.00001)
deltG=func_deltG(busNum,Loadi,I1r,I1i); deltG=func_deltG(busNum,Loadi,I1r,I1i);
end end
%% %%
L_1Z=diag(Init_Z./Init_L); % L_1Z=diag(Init_Z./Init_L);
U_1W=diag(Init_W./Init_U); % U_1W=diag(Init_W./Init_U);
%% 形成海森阵 %% 形成海森阵
deltdeltF=func_deltdeltF(busNum,fsY1,Loadi,wV1r,wV1i,wI1r,wI1i); deltdeltF=func_deltdeltF(busNum,fsY1,Loadi,wV1r,wV1i,wI1r,wI1i,V1measurement,V1r,V1i);
% deltdeltF=0; % deltdeltF=0;
%% 形成ddHy %% 形成ddHy
% ddh=func_ddh(busNum,Loadi,Init_Z,Init_W); % ddh=func_ddh(busNum,Loadi,Init_Z,Init_W);

4
func_deltF.m
View File

@ -9,8 +9,8 @@ function deltF=func_deltF(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,busNum
%% %%
deltF=[ deltF=[
%zeros(busNum*2,1); %zeros(busNum*2,1);
sparse(Loadi,1,-2*(real(V1measurement)-V1r(Loadi))./wV1r./wV1r,busNum,1); sparse(Loadi,1,2*( abs(V1measurement).^2- (V1r(Loadi).^2+ V1i(Loadi).^2) ).*(-2*V1r(Loadi))./wV1r./wV1r,busNum,1);
sparse(Loadi,1,-2*(imag(V1measurement)-V1i(Loadi))./wV1i./wV1i,busNum,1); sparse(Loadi,1,2*( abs(V1measurement).^2- (V1r(Loadi).^2+ V1i(Loadi).^2) ).*(-2*V1i(Loadi))./wV1i./wV1i,busNum,1);
%-2*( [real(Imeasurement);imag(Imeasurement)]-[I1r;I1i]); %-2*( [real(Imeasurement);imag(Imeasurement)]-[I1r;I1i]);
sparse(1:length(Loadi),1,-2*(real(I1measurement)-I1r(Loadi))./wI1r./wI1r,length(Loadi),1); sparse(1:length(Loadi),1,-2*(real(I1measurement)-I1r(Loadi))./wI1r./wI1r,length(Loadi),1);
sparse(1:length(Loadi),1,-2*(imag(I1measurement)-I1i(Loadi))./wI1i./wI1i,length(Loadi),1); sparse(1:length(Loadi),1,-2*(imag(I1measurement)-I1i(Loadi))./wI1i./wI1i,length(Loadi),1);

6
func_deltdeltF.m
View File

@ -1,9 +1,9 @@
function deltdeltF=func_deltdeltF(busNum,fsY11,Loadi,wV1r,wV1i,wI1r,wI1i) function deltdeltF=func_deltdeltF(busNum,fsY11,Loadi,wV1r,wV1i,wI1r,wI1i,V1measurement,V1r,V1i)
deltdeltF=[ deltdeltF=[
%zeros(busNum*2,busNum*2+length(Loadi)*2); %zeros(busNum*2,busNum*2+length(Loadi)*2);
sparse(Loadi,Loadi,2./[wV1r.*wV1r],busNum,busNum*2+length(Loadi)*2); sparse(Loadi,Loadi,-4*( -2*wV1r.*wV1r+ abs(V1measurement).^2-(V1r(Loadi).^2+ V1i(Loadi).^2) )./[wV1r.*wV1r],busNum,busNum),sparse(Loadi,Loadi+busNum,-4*wV1r.*( -2*wV1i)./[wV1r.*wV1r],busNum,busNum+length(Loadi)*2);
sparse(Loadi,busNum+Loadi,2./[wV1i.*wV1i],busNum,busNum*2+length(Loadi)*2); sparse(Loadi,Loadi,-4*wV1i.*( -2*wV1r)./[wV1i.*wV1i],busNum,busNum),sparse(Loadi,Loadi+busNum,-4*( -2*wV1i.*wV1i+ abs(V1measurement).^2-(V1r(Loadi).^2+ V1i(Loadi).^2) )./[wV1r.*wV1r],busNum,busNum+length(Loadi)*2);
zeros(length(Loadi)*2,busNum*2),2*eye(length(Loadi)*2)*diag(1./[wI1r.^2;wI1i.^2]); zeros(length(Loadi)*2,busNum*2),2*eye(length(Loadi)*2)*diag(1./[wI1r.^2;wI1i.^2]);
%sparse(Loadi,busNum*2+Loadi,2,busNum,busNum*2+length(Loadi)*2); %sparse(Loadi,busNum*2+Loadi,2,busNum,busNum*2+length(Loadi)*2);
]; ];

12
run.m
View File

@ -146,13 +146,13 @@ mphaseASpotLoadQ=phaseASpotLoadQ.*(1+normrnd(0,sigma,length(phaseASpotLoadQ),1))
mphaseBSpotLoadQ=phaseBSpotLoadQ.*(1+normrnd(0,sigma,length(phaseBSpotLoadQ),1)); mphaseBSpotLoadQ=phaseBSpotLoadQ.*(1+normrnd(0,sigma,length(phaseBSpotLoadQ),1));
mphaseCSpotLoadQ=phaseCSpotLoadQ.*(1+normrnd(0,sigma,length(phaseCSpotLoadQ),1)); mphaseCSpotLoadQ=phaseCSpotLoadQ.*(1+normrnd(0,sigma,length(phaseCSpotLoadQ),1));
%µçѹÁ¿²âÁ¿ %µçѹÁ¿²âÁ¿
% mVoltpA=VoltpA.*(1+normrnd(0,sigma,length(VoltpA),1)); mVoltpA=VoltpA.*(1+normrnd(0,sigma,length(VoltpA),1));
% mVoltpB=VoltpB.*(1+normrnd(0,sigma,length(VoltpB),1)); mVoltpB=VoltpB.*(1+normrnd(0,sigma,length(VoltpB),1));
% mVoltpC=VoltpC.*(1+normrnd(0,sigma,length(VoltpC),1)); mVoltpC=VoltpC.*(1+normrnd(0,sigma,length(VoltpC),1));
% %
mVoltpA=abs(VoltpA).*(1+normrnd(0,sigma,length(VoltpA),1)); % mVoltpA=abs(VoltpA).*(1+normrnd(0,sigma,length(VoltpA),1));
mVoltpB=abs(VoltpB).*(1+normrnd(0,sigma,length(VoltpB),1)).*exp(1j*-120/180*pi); % mVoltpB=abs(VoltpB).*(1+normrnd(0,sigma,length(VoltpB),1)).*exp(1j*-120/180*pi);
mVoltpC=abs(VoltpC).*(1+normrnd(0,sigma,length(VoltpC),1)).*exp(1j*+120/180*pi); % mVoltpC=abs(VoltpC).*(1+normrnd(0,sigma,length(VoltpC),1)).*exp(1j*+120/180*pi);
% %
% mVoltpA=sparse(ones(busNum,1)); % mVoltpA=sparse(ones(busNum,1));
% mVoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi); % mVoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi);