From 1c2b53d2b59558480c4a262914be93b65adcc030 Mon Sep 17 00:00:00 2001 From: "dugg@lab-desk" Date: Thu, 4 Dec 2014 11:19:08 +0800 Subject: [PATCH] =?UTF-8?q?=E4=BF=AE=E5=A4=8D=E4=BA=86=E7=9B=AE=E6=A0=87?= =?UTF-8?q?=E5=87=BD=E6=95=B0=E6=B1=82=E5=AF=BC=E7=9A=84=E9=94=99=E8=AF=AF?= =?UTF-8?q?=E3=80=82?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: dugg@lab-desk --- FormLwstate1.m | 8 ++++---- FormLzstate1.m | 8 ++++---- IPMLoop.m | 4 ++-- feeder13/data1.txt | 6 +++--- func_deltdeltF.m | 4 ++-- run.m | 6 +++--- 6 files changed, 18 insertions(+), 18 deletions(-) diff --git a/FormLwstate1.m b/FormLwstate1.m index d552ff4..3274aec 100644 --- a/FormLwstate1.m +++ b/FormLwstate1.m @@ -14,12 +14,12 @@ upper(pIi+length(Ir))=1.2*Ii(pIi); upper(nIi+length(Ir))=0.8*Ii(nIi); %太小的数都要放宽一些 -tooSmall=find(abs(Ir)<0.05); +tooSmall=find(abs(Ir)<0.0005); % upper(tooSmall)=0.99*abs(Ir(tooSmall)); -upper(tooSmall)=0.2; -tooSmall=find(abs(Ii)<0.05); +upper(tooSmall)=0.002; +tooSmall=find(abs(Ii)<0.0005); % upper(tooSmall+length(Ir))=0.99*abs(Ii(tooSmall)); -upper(tooSmall+length(Ir))=0.2; +upper(tooSmall+length(Ir))=0.002; % upper([4,5,6])=[1;1;1]; % upper([4,5,6])=0.2*Ii(nIi); diff --git a/FormLzstate1.m b/FormLzstate1.m index 2adaea9..6a2e11e 100644 --- a/FormLzstate1.m +++ b/FormLzstate1.m @@ -13,12 +13,12 @@ lower(pIi+length(Ir))=0.8*Ii(pIi); lower(nIi+length(Ir))=1.2*Ii(nIi); %太小的数都要放宽一些 -tooSmall=find(abs(Ir)<0.05); +tooSmall=find(abs(Ir)<0.0005); % lower(tooSmall)=-0.99*abs(Ir(tooSmall)); -lower(tooSmall)=-0.2; -tooSmall=find(abs(Ii)<0.05); +lower(tooSmall)=-0.002; +tooSmall=find(abs(Ii)<0.0005); % lower(tooSmall+length(Ir))=-.99*abs(Ii(tooSmall)); -lower(tooSmall+length(Ir))=-0.2; +lower(tooSmall+length(Ir))=-0.002; % lower=-ones(length(Ir)*2,1); % lower([4,5,6])=[-1;-1;-1]; diff --git a/IPMLoop.m b/IPMLoop.m index da6deee..364f926 100644 --- a/IPMLoop.m +++ b/IPMLoop.m @@ -43,7 +43,7 @@ while(abs(Gap)>0.00001) % U_1W=diag(Init_W./Init_U); %% 形成海森阵 deltdeltF=func_deltdeltF(busNum,fsY1,Loadi,wV1r,wV1i,wI1r,wI1i,V1measurement,V1r,V1i); - deltdeltF=0; +% deltdeltF=0; %% 形成ddHy % ddh=func_ddh(busNum,Loadi,Init_Z,Init_W); ddh=0; @@ -55,7 +55,7 @@ while(abs(Gap)>0.00001) end %% 开始构建deltF deltF=func_deltF(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,busNum,fsY1,Loadi,V1r,V1i,I1r,I1i); - deltF=0; +% deltF=0; %% Luu=Init_U.*Init_W+Init_u*ones(RestraintCount,1); Lul=Init_L.*Init_Z-Init_u*ones(RestraintCount,1); diff --git a/feeder13/data1.txt b/feeder13/data1.txt index 1c37398..843b67a 100644 --- a/feeder13/data1.txt +++ b/feeder13/data1.txt @@ -6,9 +6,9 @@ 0 632 0 0 0 0 -630 80 70 80 60 80 70 -631 40 10 40 10 42 10 -632 120 120 120 120 120 120 +630 90 70 80 40 30 10 +631 40 10 60 30 100 80 +632 12 12 120 120 11 11 0 diff --git a/func_deltdeltF.m b/func_deltdeltF.m index 1a9ed34..9aae298 100644 --- a/func_deltdeltF.m +++ b/func_deltdeltF.m @@ -2,8 +2,8 @@ function deltdeltF=func_deltdeltF(busNum,fsY11,Loadi,wV1r,wV1i,wI1r,wI1i,V1measu deltdeltF=[ %zeros(busNum*2,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); + sparse(Loadi,Loadi,-4*( -2*V1r(Loadi).*V1r(Loadi)+ abs(V1measurement).^2-(V1r(Loadi).^2+ V1i(Loadi).^2) )./[wV1r.*wV1r],busNum,busNum),sparse(Loadi,Loadi+busNum,-4*V1r(Loadi).*( -2*V1i(Loadi))./[wV1r.*wV1r],busNum,busNum+length(Loadi)*2); + sparse(Loadi,Loadi,-4*V1i(Loadi).*( -2*V1r(Loadi))./[wV1i.*wV1i],busNum,busNum),sparse(Loadi,Loadi+busNum,-4*( -2*V1i(Loadi).*V1i(Loadi)+ 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 05910a9..d69405e 100644 --- a/run.m +++ b/run.m @@ -150,9 +150,9 @@ 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);