From 431b0c8f780c3d11ab1d98f5d67d0347fe8e592a Mon Sep 17 00:00:00 2001 From: "dugg@lab-desk" Date: Sat, 24 May 2014 12:04:03 +0800 Subject: [PATCH] =?UTF-8?q?1.=E5=88=A0=E6=8E=89=E4=BA=92=E8=A1=A5=E7=BA=A6?= =?UTF-8?q?=E6=9D=9F=EF=BC=8C=E7=9C=8B=E7=9C=8B=E6=98=AF=E4=BB=80=E4=B9=88?= =?UTF-8?q?=E5=BD=B1=E5=93=8D=E4=BA=86=E6=94=B6=E6=95=9B=E6=80=A7=E3=80=82?= =?UTF-8?q?=E5=BD=93=E5=89=8D=E6=94=B6=E6=95=9B=E6=AC=A1=E6=95=B0=E4=B8=BA?= =?UTF-8?q?13=20Mat=5FG=3D[=20=20=20=20=20sparse(PD(Loadi))-mPD(Loadi)-big?= =?UTF-8?q?M*PDbi-0.1;=20=20=20=20=20sparse(PD(Loadi))-mPD(Loadi)+bigM*PDb?= =?UTF-8?q?i+0.1;=20=20=20=20=20sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1?= =?UTF-8?q?;=20=20=20=20=20sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;=20?= =?UTF-8?q?=20=20=20=20Volt'-mVolt'-bigM*Vbi-0.03;=20=20=20=20=20Volt'-mVo?= =?UTF-8?q?lt'+bigM*Vbi+0.03;=20=20=20=20=20Vbi;=20=20=20=20=20PDbi;=20=20?= =?UTF-8?q?=20=20=20QDbi;=20%=20=20=20=20=20(1-Vbi).*Vbi;=20%=20=20=20=20?= =?UTF-8?q?=20(1-PDbi).*PDbi;=20%=20=20=20=20=20(1-QDbi).*QDbi;=20=20=20?= =?UTF-8?q?=20=20];=202.=E6=9B=B4=E6=96=B0=E4=BA=86.gitignore?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: dugg@lab-desk --- .gitignore | 3 ++- FormG.m | 6 +++--- FormLw.m | 10 +++++----- FormLz.m | 10 +++++----- OPF_Init.m | 2 +- func_ddg.m | 19 ++++++++++--------- func_deltG.m | 8 +++++++- 7 files changed, 33 insertions(+), 25 deletions(-) diff --git a/.gitignore b/.gitignore index 21229f7..43eb260 100644 --- a/.gitignore +++ b/.gitignore @@ -1 +1,2 @@ -*.asv \ No newline at end of file +*.asv +*.7z \ No newline at end of file diff --git a/FormG.m b/FormG.m index 9add28b..c10690a 100644 --- a/FormG.m +++ b/FormG.m @@ -10,8 +10,8 @@ Mat_G=[ Vbi; PDbi; QDbi; - (1-Vbi).*Vbi; - (1-PDbi).*PDbi; - (1-QDbi).*QDbi; +% (1-Vbi).*Vbi; +% (1-PDbi).*PDbi; +% (1-QDbi).*QDbi; ]; end \ No newline at end of file diff --git a/FormLw.m b/FormLw.m index 4484aa1..9bbbce2 100644 --- a/FormLw.m +++ b/FormLw.m @@ -18,11 +18,11 @@ indQD=find(realQD>0); QDU(indQD(3:12:end))=1.55*realQD(indQD(3:12:end)); QDU(indQD(9:12:end))=1.05*realQD(indQD(9:12:end)); t1=([0*PDU',100*ones(1,length(PDU)),0*QDU',100*ones(1,length(QDU)),0*VoltU,100*ones(1,length(VoltU)),1*ones(1,Busnum+length(Loadi)*2)])'; -t1=[t1; - sparse(ones(Busnum,1)); - sparse(ones(length(Loadi),1)); - sparse(ones(length(Loadi),1)); - ]; +% t1=[t1; +% sparse(ones(Busnum,1)); +% sparse(ones(length(Loadi),1)); +% sparse(ones(length(Loadi),1)); +% ]; t2=Mat_G+Init_U'-t1; Lw=t2; diff --git a/FormLz.m b/FormLz.m index 4df63fb..e091443 100644 --- a/FormLz.m +++ b/FormLz.m @@ -22,11 +22,11 @@ QDL(indQD(3:12:end))=0.95*realQD(indQD(3:12:end)); QDL(indQD(9:12:end))=0.95*realQD(indQD(9:12:end)); % QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF; t1=([-100*ones(1,length(PDL)),0*PDL',-100*ones(1,length(QDL)),0*QDL',-100*ones(1,length(VoltL)),0*VoltL,0*ones(1,Busnum+length(Loadi)*2)])'; -t1=[t1; - sparse(0*ones(Busnum,1)); - sparse(0*ones(length(Loadi),1)); - sparse(0*ones(length(Loadi),1)); -]; +% t1=[t1; +% sparse(0*ones(Busnum,1)); +% sparse(0*ones(length(Loadi),1)); +% sparse(0*ones(length(Loadi),1)); +% ]; t2=Mat_G-Init_L'-t1; Lz=t2; diff --git a/OPF_Init.m b/OPF_Init.m index 7f9ddba..da090d6 100644 --- a/OPF_Init.m +++ b/OPF_Init.m @@ -3,7 +3,7 @@ Loadi=find(QD~=0 | PD~=0); %Loadi=[1:Busnum]'; RestraintCount=size(Loadi,1)*4+Busnum*2+Busnum+length(Loadi)*2; %约束条件数,放开所有QD %互补约束数量 -RestraintCount=RestraintCount+Busnum+length(Loadi)*2; +RestraintCount=RestraintCount;%+Busnum+length(Loadi)*2; t_Bal_volt=Volt(Balance); Volt=sparse(1*ones(1,Busnum)); Volt(Balance)=t_Bal_volt; diff --git a/func_ddg.m b/func_ddg.m index a75d304..5142d31 100644 --- a/func_ddg.m +++ b/func_ddg.m @@ -1,11 +1,12 @@ function ddgzw=func_ddg(Busnum,ContrlCount,Loadi,Init_Z,Init_W) -CCVbiZW=Init_Z(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+1:length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum)+Init_W(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+1:length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum); -CCPDbiZW=Init_Z(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+1:length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+length(Loadi))+Init_W(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+1:length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+length(Loadi)); -CCQDbiZW=Init_Z(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+length(Loadi)+1:end)+Init_W(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+length(Loadi)+1:end); -ddgzw=[ -sparse(length(Loadi)*2+Busnum*2,ContrlCount); -sparse(Busnum,length(Loadi)*2+Busnum*2),sparse(-2*eye(Busnum))*diag(CCVbiZW),sparse(Busnum,length(Loadi)*2); -sparse(length(Loadi),length(Loadi)*2+Busnum*2+Busnum),sparse(-2*eye(length(Loadi)))*diag(CCPDbiZW),sparse(length(Loadi),length(Loadi)); -sparse(length(Loadi),length(Loadi)*2+Busnum*2+Busnum+length(Loadi)),sparse(-2*eye(length(Loadi)))*diag(CCQDbiZW); -]; +% CCVbiZW=Init_Z(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+1:length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum)+Init_W(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+1:length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum); +% CCPDbiZW=Init_Z(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+1:length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+length(Loadi))+Init_W(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+1:length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+length(Loadi)); +% CCQDbiZW=Init_Z(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+length(Loadi)+1:end)+Init_W(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+Busnum+length(Loadi)+1:end); +% ddgzw=[ +% sparse(length(Loadi)*2+Busnum*2,ContrlCount); +% sparse(Busnum,length(Loadi)*2+Busnum*2),sparse(-2*eye(Busnum))*diag(CCVbiZW),sparse(Busnum,length(Loadi)*2); +% sparse(length(Loadi),length(Loadi)*2+Busnum*2+Busnum),sparse(-2*eye(length(Loadi)))*diag(CCPDbiZW),sparse(length(Loadi),length(Loadi)); +% sparse(length(Loadi),length(Loadi)*2+Busnum*2+Busnum+length(Loadi)),sparse(-2*eye(length(Loadi)))*diag(CCQDbiZW); +% ]; +ddgzw=0; end \ No newline at end of file diff --git a/func_deltG.m b/func_deltG.m index d1dcb3e..0ca1725 100644 --- a/func_deltG.m +++ b/func_deltG.m @@ -95,5 +95,11 @@ deltG=[dg3_dPD,dg32_dPD,dg4_dPD,dg42_dPD,dg5_dPD,dg6_dPD,dg7_dPD,dgPDbi_dPD,dgQD dg3_dPDbi,dg32_dPDbi,dg4_dPDbi,dg42_dPDbi,dg5_dPDbi,dg6_dPDbi,dg7_dPDbi,dgPDbi_dPDbi,dgQDbi_dPDbi,dgCCVbi_dPDbi,dgCCPDbi_dPDbi,dgCCQDbi_dPDbi; dg3_dQDbi,dg32_dQDbi,dg4_dQDbi,dg42_dQDbi,dg5_dQDbi,dg6_dQDbi,dg7_dQDbi,dgPDbi_dQDbi,dgQDbi_dQDbi,dgCCVbi_dQDbi,dgCCPDbi_dQDbi,dgCCQDbi_dQDbi; ]; - +deltG=[dg3_dPD,dg32_dPD,dg4_dPD,dg42_dPD,dg5_dPD,dg6_dPD,dg7_dPD,dgPDbi_dPD,dgQDbi_dPD; + dg3_dQD,dg32_dQD,dg4_dQD,dg42_dQD,dg5_dQD,dg6_dQD,dg7_dQD,dgPDbi_dQD,dgQDbi_dQD; + dg3_dx,dg32_dx,dg4_dx,dg42_dx,dg5_dx,dg6_dx,dg7_dx,dgPDbi_dx,dgQDbi_dx; + dg3_dvbi,dg32_dvbi,dg4_dvbi,dg42_dvbi,dg5_dvbi,dg6_dvbi,dg7_dvbi,dgPDbi_dvbi,dgQDbi_dvbi; + dg3_dPDbi,dg32_dPDbi,dg4_dPDbi,dg42_dPDbi,dg5_dPDbi,dg6_dPDbi,dg7_dPDbi,dgPDbi_dPDbi,dgQDbi_dPDbi; + dg3_dQDbi,dg32_dQDbi,dg4_dQDbi,dg42_dQDbi,dg5_dQDbi,dg6_dQDbi,dg7_dQDbi,dgPDbi_dQDbi,dgQDbi_dQDbi; +]; end \ No newline at end of file