From 3db4877bb72a0da7a69d18918ac653768c171a52 Mon Sep 17 00:00:00 2001 From: "dugg@lab-desk" Date: Mon, 10 Nov 2014 22:33:39 +0800 Subject: [PATCH] =?UTF-8?q?=E6=94=B9=E4=BA=86=E4=B8=80=E4=B8=8B=E7=9B=AE?= =?UTF-8?q?=E6=A0=87=E5=87=BD=E6=95=B0=EF=BC=8C=E5=8F=88=E4=B8=8D=E6=94=B6?= =?UTF-8?q?=E6=95=9B=E4=BA=86=E3=80=82?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: dugg@lab-desk --- IPMLoop.m | 6 +++--- func_deltF.m | 4 ++-- func_deltdeltF.m | 6 +++--- run.m | 12 ++++++------ 4 files changed, 14 insertions(+), 14 deletions(-) diff --git a/IPMLoop.m b/IPMLoop.m index 939ca8c..364f926 100644 --- a/IPMLoop.m +++ b/IPMLoop.m @@ -39,10 +39,10 @@ while(abs(Gap)>0.00001) deltG=func_deltG(busNum,Loadi,I1r,I1i); end %% - L_1Z=diag(Init_Z./Init_L); - U_1W=diag(Init_W./Init_U); +% L_1Z=diag(Init_Z./Init_L); +% 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; %% 形成ddHy % ddh=func_ddh(busNum,Loadi,Init_Z,Init_W); diff --git a/func_deltF.m b/func_deltF.m index 8f3c15b..bb58902 100644 --- a/func_deltF.m +++ b/func_deltF.m @@ -9,8 +9,8 @@ function deltF=func_deltF(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,busNum %% deltF=[ %zeros(busNum*2,1); - sparse(Loadi,1,-2*(real(V1measurement)-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*V1r(Loadi))./wV1r./wV1r,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]); 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); diff --git a/func_deltdeltF.m b/func_deltdeltF.m index b034cf9..1a9ed34 100644 --- a/func_deltdeltF.m +++ b/func_deltdeltF.m @@ -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=[ %zeros(busNum*2,busNum*2+length(Loadi)*2); - sparse(Loadi,Loadi,2./[wV1r.*wV1r],busNum,busNum*2+length(Loadi)*2); - sparse(Loadi,busNum+Loadi,2./[wV1i.*wV1i],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,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]); %sparse(Loadi,busNum*2+Loadi,2,busNum,busNum*2+length(Loadi)*2); ]; diff --git a/run.m b/run.m index df9c356..22023b8 100644 --- a/run.m +++ b/run.m @@ -146,13 +146,13 @@ mphaseASpotLoadQ=phaseASpotLoadQ.*(1+normrnd(0,sigma,length(phaseASpotLoadQ),1)) mphaseBSpotLoadQ=phaseBSpotLoadQ.*(1+normrnd(0,sigma,length(phaseBSpotLoadQ),1)); mphaseCSpotLoadQ=phaseCSpotLoadQ.*(1+normrnd(0,sigma,length(phaseCSpotLoadQ),1)); %电压量测量 -% mVoltpA=VoltpA.*(1+normrnd(0,sigma,length(VoltpA),1)); -% mVoltpB=VoltpB.*(1+normrnd(0,sigma,length(VoltpB),1)); -% mVoltpC=VoltpC.*(1+normrnd(0,sigma,length(VoltpC),1)); +mVoltpA=VoltpA.*(1+normrnd(0,sigma,length(VoltpA),1)); +mVoltpB=VoltpB.*(1+normrnd(0,sigma,length(VoltpB),1)); +mVoltpC=VoltpC.*(1+normrnd(0,sigma,length(VoltpC),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); -mVoltpC=abs(VoltpC).*(1+normrnd(0,sigma,length(VoltpC),1)).*exp(1j*+120/180*pi); +% 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); +% mVoltpC=abs(VoltpC).*(1+normrnd(0,sigma,length(VoltpC),1)).*exp(1j*+120/180*pi); % % mVoltpA=sparse(ones(busNum,1)); % mVoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi);