From f4d463540d0a15d844f41151feca090bb8cc6126 Mon Sep 17 00:00:00 2001 From: "dugg@lab-desk" Date: Sat, 24 May 2014 14:05:28 +0800 Subject: [PATCH] =?UTF-8?q?14=E6=AC=A1=E6=94=B6=E6=95=9B=20Mat=5FG=3D[=20?= =?UTF-8?q?=20=20=20=20sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;=20=20?= =?UTF-8?q?=20=20=20sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;=20=20=20?= =?UTF-8?q?=20=20sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;=20=20=20=20?= =?UTF-8?q?=20sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;=20=20=20=20=20Vo?= =?UTF-8?q?lt'-mVolt'-bigM*Vbi-0.03;=20=20=20=20=20Volt'-mVolt'+bigM*Vbi+0?= =?UTF-8?q?.03;=20=20=20=20=20Vbi;=20=20=20=20=20PDbi;=20=20=20=20=20QDbi;?= =?UTF-8?q?=20=20=20=20=20(1-Vbi).*Vbi;=20=20=20=20=20(1-PDbi).*PDbi;=20%?= =?UTF-8?q?=20=20=20=20=20(1-QDbi).*QDbi;=20=20=20=20=20];?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: dugg@lab-desk --- FormG.m | 2 +- FormLw.m | 2 +- FormLz.m | 2 +- OPF_Init.m | 2 +- func_ddg.m | 4 ++-- func_deltG.m | 12 ++++++------ 6 files changed, 12 insertions(+), 12 deletions(-) diff --git a/FormG.m b/FormG.m index 1070e2b..78ab555 100644 --- a/FormG.m +++ b/FormG.m @@ -11,7 +11,7 @@ Mat_G=[ PDbi; QDbi; (1-Vbi).*Vbi; -% (1-PDbi).*PDbi; + (1-PDbi).*PDbi; % (1-QDbi).*QDbi; ]; end \ No newline at end of file diff --git a/FormLw.m b/FormLw.m index 49b60b7..f4970c1 100644 --- a/FormLw.m +++ b/FormLw.m @@ -20,7 +20,7 @@ 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)); % sparse(ones(length(Loadi),1)); ]; t2=Mat_G+Init_U'-t1; diff --git a/FormLz.m b/FormLz.m index 8401046..753dd74 100644 --- a/FormLz.m +++ b/FormLz.m @@ -24,7 +24,7 @@ QDL(indQD(9:12:end))=0.95*realQD(indQD(9:12:end)); 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)); % sparse(0*ones(length(Loadi),1)); ]; t2=Mat_G-Init_L'-t1; diff --git a/OPF_Init.m b/OPF_Init.m index 7476b68..a471100 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)*1; 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 b4ddd9c..96d1fc2 100644 --- a/func_ddg.m +++ b/func_ddg.m @@ -1,6 +1,6 @@ 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)); +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); @@ -11,7 +11,7 @@ CCVbiZW=Init_Z(length(Loadi)*4+Busnum*2+Busnum+length(Loadi)*2+1:length(Loadi)*4 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(-0*eye(length(Loadi))),sparse(length(Loadi),length(Loadi)); +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(-0*eye(length(Loadi))); ]; end \ No newline at end of file diff --git a/func_deltG.m b/func_deltG.m index 0278149..e8c9af3 100644 --- a/func_deltG.m +++ b/func_deltG.m @@ -97,11 +97,11 @@ dgCCQDbi_dQDbi=sparse(diag(1-2*QDbi)); % 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,dgCCVbi_dPD; - dg3_dQD,dg32_dQD,dg4_dQD,dg42_dQD,dg5_dQD,dg6_dQD,dg7_dQD,dgPDbi_dQD,dgQDbi_dQD,dgCCVbi_dQD; - dg3_dx,dg32_dx,dg4_dx,dg42_dx,dg5_dx,dg6_dx,dg7_dx,dgPDbi_dx,dgQDbi_dx,dgCCVbi_dx; - dg3_dvbi,dg32_dvbi,dg4_dvbi,dg42_dvbi,dg5_dvbi,dg6_dvbi,dg7_dvbi,dgPDbi_dvbi,dgQDbi_dvbi,dgCCVbi_dvbi; - dg3_dPDbi,dg32_dPDbi,dg4_dPDbi,dg42_dPDbi,dg5_dPDbi,dg6_dPDbi,dg7_dPDbi,dgPDbi_dPDbi,dgQDbi_dPDbi,dgCCVbi_dPDbi; - dg3_dQDbi,dg32_dQDbi,dg4_dQDbi,dg42_dQDbi,dg5_dQDbi,dg6_dQDbi,dg7_dQDbi,dgPDbi_dQDbi,dgQDbi_dQDbi,dgCCVbi_dQDbi; +deltG=[dg3_dPD,dg32_dPD,dg4_dPD,dg42_dPD,dg5_dPD,dg6_dPD,dg7_dPD,dgPDbi_dPD,dgQDbi_dPD,dgCCVbi_dPD,dgCCPDbi_dPD; + dg3_dQD,dg32_dQD,dg4_dQD,dg42_dQD,dg5_dQD,dg6_dQD,dg7_dQD,dgPDbi_dQD,dgQDbi_dQD,dgCCVbi_dQD,dgCCPDbi_dQD; + dg3_dx,dg32_dx,dg4_dx,dg42_dx,dg5_dx,dg6_dx,dg7_dx,dgPDbi_dx,dgQDbi_dx,dgCCVbi_dx,dgCCPDbi_dx; + dg3_dvbi,dg32_dvbi,dg4_dvbi,dg42_dvbi,dg5_dvbi,dg6_dvbi,dg7_dvbi,dgPDbi_dvbi,dgQDbi_dvbi,dgCCVbi_dvbi,dgCCPDbi_dvbi; + dg3_dPDbi,dg32_dPDbi,dg4_dPDbi,dg42_dPDbi,dg5_dPDbi,dg6_dPDbi,dg7_dPDbi,dgPDbi_dPDbi,dgQDbi_dPDbi,dgCCVbi_dPDbi,dgCCPDbi_dPDbi; + dg3_dQDbi,dg32_dQDbi,dg4_dQDbi,dg42_dQDbi,dg5_dQDbi,dg6_dQDbi,dg7_dQDbi,dgPDbi_dQDbi,dgQDbi_dQDbi,dgCCVbi_dQDbi,dgCCPDbi_dQDbi; ]; end \ No newline at end of file