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

1.修复了这个地方的bug 

lPD=abs(RealPD*3*sigma);
uPD=abs(RealPD*3*sigma);
lQD=abs(RealQD*3*sigma);
uQD=abs(RealQD*3*sigma);
lVolt=abs(rVolt'*3*sigma);
uVolt=abs(rVolt'*3*sigma);
2.加入了考虑发电机的情况，其实也就是加入负数负荷。

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk 2014-07-20 21:08:16 +08:00
parent a0a0d45e7d
commit 0e24eed4bf
1 changed files with 11 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
OPF.m
View File

@ -6,7 +6,7 @@ clear
%% %%
thesis=ForThesis(1,62); thesis=ForThesis(1,62);
[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:\ËãÀý\«¼Ê911_2751267_2012-09-05\newFIle9-2.txt'); pf('E:\ËãÀý\17\17.csv');
% pf('E:\算例\柳金Ⅰ926_21671693_2012-09-06\newFIle16.txt'); % pf('E:\算例\柳金Ⅰ926_21671693_2012-09-06\newFIle16.txt');
%pf('D:\Project\青秀降损项目\最小化潮流\最小潮流算例\原始\津头站津视922（3-1）_0.5_120%.txt'); %pf('D:\Project\青秀降损项目\最小化潮流\最小潮流算例\原始\津头站津视922（3-1）_0.5_120%.txt');
%pf('D:\Project\最小化潮流\最小潮流算例\仙海919.txt'); %pf('D:\Project\最小化潮流\最小潮流算例\仙海919.txt');
@ -56,12 +56,12 @@ mVolt=rVolt.*(1+normrnd(0,sigma,length(rVolt),1))';
mPD=PD0; mPD=PD0;
mQD=QD0; mQD=QD0;
%加估计上下界 %加估计上下界
lPD=RealPD*(1-3*sigma); lPD=abs(RealPD*3*sigma);
uPD=RealPD*(1+3*sigma); uPD=abs(RealPD*3*sigma);
lQD=RealQD*(1-3*sigma); lQD=abs(RealQD*3*sigma);
uQD=RealQD*(1+3*sigma); uQD=abs(RealQD*3*sigma);
lVolt=rVolt'*(1-3*sigma); lVolt=abs(rVolt'*3*sigma);
uVolt=rVolt'*(1+3*sigma); uVolt=abs(rVolt'*3*sigma);
%错误数据 %错误数据
%mVolt(2)=5; %mVolt(2)=5;
bigM=10; bigM=10;
@ -148,10 +148,11 @@ while(abs(Gap)>Precision*1)
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;
if KK>8 if KK>8
eps=eps*0.2; eps=eps*0.3;
if abs(eps)<1e-6 if abs(eps)<1e-6
eps=1e-6; eps=1e-5;
end end
eps;
end end
Lz=FormLz(Mat_G,Init_L,GenL,Busnum,PVQL,PD,RealPD,RealQD,Loadi,KK,PF,eps); Lz=FormLz(Mat_G,Init_L,GenL,Busnum,PVQL,PD,RealPD,RealQD,Loadi,KK,PF,eps);
Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD,PD0,QD0,Loadi,KK,PF,eps); Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD,PD0,QD0,Loadi,KK,PF,eps);
@ -167,6 +168,7 @@ while(abs(Gap)>Precision*1)
KK=KK+1; KK=KK+1;
end end
toc toc
eps
obj=sum(Vbi)+sum(PDbi)+sum(QDbi); obj=sum(Vbi)+sum(PDbi)+sum(QDbi);
fprintf('目标函数值 %.2f\n',full(obj)); fprintf('目标函数值 %.2f\n',full(obj));