facat
/
stateestimateionyalmip-lu9-forcomparision
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

加上噪音 

Signed-off-by: facat <dugg@21cn.com>
This commit is contained in:
facat 2013-04-06 21:53:24 +08:00
parent 8990b8a278
commit 4db769846c
2 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
OPF_Init.m
View File

@ -33,8 +33,8 @@ wQG=1*ones(size(PVi,1),1);
%randInt=randperm(size(Loadi,1)); %randInt=randperm(size(Loadi,1));
%randPDind=randInt(1:10); %randPDind=randInt(1:10);
randPDind=0; randPDind=0;
wPD=1*ones(sum(PD~=0),1); wPD=1/0.05^2*ones(sum(PD~=0),1);
wQD=1.5*ones(sum(QD~=0),1); wQD=1/0.05^2*ones(sum(QD~=0),1);
%wD(randPDind)=0;%一些负荷不约束 %wD(randPDind)=0;%一些负荷不约束
%wD(7)=0; %wD(7)=0;
% wD(11)=0; % wD(11)=0;

3
Run_YALMIP.m
View File

@ -41,6 +41,9 @@ UAngel=sdpvar(Busnum,1);
PD=sdpvar(Busnum,1); PD=sdpvar(Busnum,1);
QD=sdpvar(Busnum,1); QD=sdpvar(Busnum,1);
AngleIJ=sdpvar(Busnum,Busnum,'full'); AngleIJ=sdpvar(Busnum,Busnum,'full');
%% ¼ÓÔëÒô
PD0=(1+normrnd(0,0.05,length(PD0),1)).*PD0;
QD0=(1+normrnd(0,0.05,length(QD0),1)).*QD0;
%% 目标函数 %% 目标函数
Objective=ObjectiveFun(PG,PG0,PGi,QG,QG0,PVi,PD,PD0,QD,QD0,wPG,wQG,wPD,wQD,Loadi); Objective=ObjectiveFun(PG,PG0,PGi,QG,QG0,PVi,PD,PD0,QD,QD0,wPG,wQG,wPD,wQD,Loadi);
%AngleIJ=sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum); %AngleIJ=sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum);