facat
/
pes2014-debug2-laplace
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

做电压上下界错误的测试，这个参数是最好的。 

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk 2015-02-07 14:47:20 +08:00
parent 1fd274a00d
commit a964aa453d
3 changed files with 17 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
OPF.m
View File

@ -17,7 +17,7 @@ for badDataNode=1:1
close close
[kmax,Precision,UAngel,Volt,Busnum,PVi,PVu,Balance,Y,Angle,P0,Q0,r,c,GB,Linei,Linej,Transfori,Transforj,GenU,GenL,GenC,PG,QG,PD,QD,CenterA,PGi,PVQU,PVQL,Liner,Linex,Lineb,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0]= ... [kmax,Precision,UAngel,Volt,Busnum,PVi,PVu,Balance,Y,Angle,P0,Q0,r,c,GB,Linei,Linej,Transfori,Transforj,GenU,GenL,GenC,PG,QG,PD,QD,CenterA,PGi,PVQU,PVQL,Liner,Linex,Lineb,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0]= ...
pf('E:\ËãÀý\feeder33\feeder33.txt'); pf('E:\ËãÀý\feeder33\feeder33.txt');
sigma=0.001; sigma=0.01;
RealPD=PD; RealPD=PD;
RealQD=QD; RealQD=QD;
rVolt=Volt; rVolt=Volt;
@ -39,7 +39,7 @@ for badDataNode=1:1
% mVolt(badDataNode)=rVolt(badDataNode)*(1-sigma*6); % mVolt(badDataNode)=rVolt(badDataNode)*(1-sigma*6);
%% »­Case AµÄͼ %% »­Case AµÄͼ
% figure('Color',[1 1 1]); % figure('Color',[1 1 1]);
[Busnum,Loadi,Volt,PD,QD,rVolt,UAngel,RealPD,RealQD,rUAngel,Vbi,PDbi,QDbi,plotGapA,isConverge]=subOPF([1:33],PD0,QD0,mVolt,sigma);%全部有 [Busnum,Loadi,Volt,PD,QD,rVolt,UAngel,RealPD,RealQD,rUAngel,Vbi,PDbi,QDbi,plotGapA,isConverge]=subOPF([22],PD0,QD0,mVolt,sigma);%全部有
% badDataResult(I,badDataNode)=sum(Vbi); % badDataResult(I,badDataNode)=sum(Vbi);
% badDataLocation(1:33,I)=Vbi; % badDataLocation(1:33,I)=Vbi;
% badDataLocation(34,I)=sum(abs((rVolt-Volt)./rVolt./length(rVolt)))+sum(abs( (UAngel(2:33)-rUAngel(2:33))./rUAngel(2:33)./length(rUAngel(2:33)))); % badDataLocation(34,I)=sum(abs((rVolt-Volt)./rVolt./length(rVolt)))+sum(abs( (UAngel(2:33)-rUAngel(2:33))./rUAngel(2:33)./length(rUAngel(2:33))));

6
openfile2.m
View File

@ -56,6 +56,12 @@ QD=QD/Base;
PD=sparse(PD); PD=sparse(PD);
QD=sparse(QD); QD=sparse(QD);
%QD=PD*sqrt(1-.85^2)/.85; %QD=PD*sqrt(1-.85^2)/.85;
%%
%增加大负荷
PD(22)=PD(22)*65;
%%
PG=sparse(PG); PG=sparse(PG);
QG=sparse(QG); QG=sparse(QG);
%% pv节点功率参数矩阵 %% pv节点功率参数矩阵

16
subOPF.m
View File

@ -46,7 +46,7 @@ Gap=(Init_L*Init_Z'-Init_U*Init_W');
KK=0; KK=0;
plotGap=zeros(1,60); plotGap=zeros(1,60);
ContrlCount=size(Loadi,1)*2+Busnum*2+Busnum+length(Loadi)*2; ContrlCount=size(Loadi,1)*2+Busnum*2+Busnum+length(Loadi)*2;
kmax=150; kmax=350;
Precision=Precision/1; Precision=Precision/1;
%% 加误差 %% 加误差
%找DG %找DG
@ -88,6 +88,8 @@ uQD(DGi)=abs(mQD(DGi)*3*sigma*0.3);
lVolt(DGi)=abs(mVolt(DGi)'*3*sigma*0.3); lVolt(DGi)=abs(mVolt(DGi)'*3*sigma*0.3);
uVolt(DGi)=abs(mVolt(DGi)'*3*sigma*0.3); uVolt(DGi)=abs(mVolt(DGi)'*3*sigma*0.3);
%% 不测量的数据 %% 不测量的数据
% mPD(22)=mPD(22)/65;
lPD(noMeasurei)=0.15*mPD(noMeasurei);%15% lPD(noMeasurei)=0.15*mPD(noMeasurei);%15%
uPD(noMeasurei)=0.15*mPD(noMeasurei); uPD(noMeasurei)=0.15*mPD(noMeasurei);
lQD(noMeasurei)=0.15*mQD(noMeasurei); lQD(noMeasurei)=0.15*mQD(noMeasurei);
@ -95,10 +97,10 @@ uQD(noMeasurei)=0.15*mQD(noMeasurei);
%电压伪量测错误 %电压伪量测错误
mVolt(noMeasurei)=0.95; mVolt(noMeasurei)=.95;
lVolt(noMeasurei)=0.7*mVolt(noMeasurei);%电压在0.93~1.07 lVolt(noMeasurei)=0.07*mVolt(noMeasurei);%电压在0.93~1.07
uVolt(noMeasurei)=0.7*mVolt(noMeasurei); uVolt(noMeasurei)=0.07*mVolt(noMeasurei);
@ -138,7 +140,7 @@ while(abs(Gap)>Precision*10)
%% 形成方程矩阵 %% 形成方程矩阵
Luu=Init_U'.*Init_W'+Init_u*ones(RestraintCount,1); Luu=Init_U'.*Init_W'+Init_u*ones(RestraintCount,1);
Lul=Init_L'.*Init_Z'-Init_u*ones(RestraintCount,1); Lul=Init_L'.*Init_Z'-Init_u*ones(RestraintCount,1);
bigM=1; bigM=0.4;
Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD,mQD,uPD,lPD,uQD,lQD,uVolt,lVolt); Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD,mQD,uPD,lPD,uQD,lQD,uVolt,lVolt);
Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi); Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi);
Ly=Mat_H; Ly=Mat_H;
@ -191,8 +193,8 @@ while(abs(Gap)>Precision*10)
%% 形成方程矩阵 %% 形成方程矩阵
Luu=Init_U'.*Init_W'+Init_u*ones(RestraintCount,1); Luu=Init_U'.*Init_W'+Init_u*ones(RestraintCount,1);
Lul=Init_L'.*Init_Z'-Init_u*ones(RestraintCount,1); Lul=Init_L'.*Init_Z'-Init_u*ones(RestraintCount,1);
bigM=1; bigM=0.4;
eps=Gap*0.1; eps=Gap*0.001;
Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD,mQD,uPD,lPD,uQD,lQD,uVolt,lVolt); Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD,mQD,uPD,lPD,uQD,lQD,uVolt,lVolt);
Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi); Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi);
Ly=Mat_H; Ly=Mat_H;