facat
/
ipmthreephasese
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

做了点小测试 

Signed-off-by: facat@lab.com <facat@lab.com>
This commit is contained in:
facat@lab.com 2014-11-04 22:13:43 +08:00
parent acab2b5d17
commit f255d7285d
1 changed files with 4 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
OPF.m
View File

@ -221,12 +221,9 @@ mQD3P=rPD3P.*QD3PSigma;
wVolt=1./(abs(mVoltABCV*sigma).^2); wVolt=1./(abs(mVoltABCV*sigma).^2);
wPD=1./(abs(mPD3P*sigma).^2); wPD=1./(abs(mPD3P*sigma).^2);
wQD=1./(abs(mQD3P*sigma).^2); wQD=1./(abs(mQD3P*sigma).^2);
wVolt(1:3:end)=0; wVolt(5)=0;
wPD(1:3:end)=0; wPD(5)=0;
wQD(1:3:end)=0; wQD(5)=0;
% wVolt=1./abs(VoltSigma).^2;
% wPD=1./abs(PD3PSigma).^2;
% wQD=1./abs(QD3PSigma).^2;
%% %%
RestraintCount=size(Loadi,1)*2; %约束条件数,放开所有QD RestraintCount=size(Loadi,1)*2; %约束条件数,放开所有QD
Init_Z=sparse(ones(1,RestraintCount)); Init_Z=sparse(ones(1,RestraintCount));
@ -237,7 +234,7 @@ Init_Y=sparse(1,2*Busnum*3);%
Gap=(Init_L*Init_Z'-Init_U*Init_W'); Gap=(Init_L*Init_Z'-Init_U*Init_W');
KK=0; KK=0;
ContrlCount=size(Loadi,1)*2+Busnum*6; ContrlCount=size(Loadi,1)*2+Busnum*6;
kmax=300; kmax=900;
%% %%
Precision=1e-5; Precision=1e-5;
CenterA=0.1; CenterA=0.1;