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

1.上下界改为1.2到0.8 

2.用量测值确定上下界

Signed-off-by: dmy@lab <dmy@lab.lab>
This commit is contained in:
dmy@lab 2015-04-03 10:06:03 +08:00
parent 7652f913ea
commit cc7312d070
3 changed files with 9 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
FormLw.m
View File

@ -1,4 +1,4 @@
function Lw=FormLw(Mat_G,Init_U,Loadi,rPD3P,rQD3P) function Lw=FormLw(Mat_G,Init_U,Loadi,mPD3P,mQD3P)
% KK=999; % KK=999;
% %PU=GenU(:,2);%发电机有功上界 % %PU=GenU(:,2);%发电机有功上界
% %
@ -26,8 +26,8 @@ function Lw=FormLw(Mat_G,Init_U,Loadi,rPD3P,rQD3P)
% QDU=1.0*PD(Loadi).*sqrt(1 -PF.^2)./PF; % QDU=1.0*PD(Loadi).*sqrt(1 -PF.^2)./PF;
% PDU=1.8*ones(length(Loadi),1); % PDU=1.8*ones(length(Loadi),1);
% QDU=1.8*ones(length(Loadi),1); % QDU=1.8*ones(length(Loadi),1);
PDU=1.1*rPD3P; PDU=1.2*mPD3P;
QDU=1.1*rQD3P; QDU=1.2*mQD3P;
t1=([PDU',QDU'])'; t1=([PDU',QDU'])';
t2=Mat_G+Init_U'-t1; t2=Mat_G+Init_U'-t1;
Lw=t2; Lw=t2;

6
FormLz.m
View File

@ -1,4 +1,4 @@
function Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P) function Lz=FormLz(Mat_G,Init_L,Loadi,mPD3P,mQD3P)
% KK=999; % KK=999;
% %
% VoltL=(0.9)*ones(1,Busnum); % VoltL=(0.9)*ones(1,Busnum);
@ -23,8 +23,8 @@ function Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P)
% QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF; % QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF;
% PDL=0*ones(length(Loadi),1); % PDL=0*ones(length(Loadi),1);
% QDL=0*ones(length(Loadi),1); % QDL=0*ones(length(Loadi),1);
PDL=0.9*rPD3P; PDL=0.8*mPD3P;
QDL=0.9*rQD3P; QDL=0.8*mQD3P;
t1=([PDL',QDL'])'; t1=([PDL',QDL'])';
t2=Mat_G-Init_L'-t1; t2=Mat_G-Init_L'-t1;
Lz=t2; Lz=t2;

6
OPF.m
View File

@ -210,7 +210,7 @@ rQD3P(3:3:end)=phaseCSpotLoadQ;
rPD3P=rPD3P(Loadi); rPD3P=rPD3P(Loadi);
rQD3P=rQD3P(Loadi); rQD3P=rQD3P(Loadi);
%量测量 %量测量
sigma=0.03; sigma=0.01;
VoltSigma=(1+normrnd(0,sigma,length(rVoltABCV),1)); VoltSigma=(1+normrnd(0,sigma,length(rVoltABCV),1));
mVoltABCV=rVoltABCV.*VoltSigma; mVoltABCV=rVoltABCV.*VoltSigma;
PD3PSigma=(1+normrnd(0,sigma,length(rPD3P),1)); PD3PSigma=(1+normrnd(0,sigma,length(rPD3P),1));
@ -272,8 +272,8 @@ while(abs(Gap)>Precision)
Mat_G=FormG(PD3P,QD3P,Loadi); Mat_G=FormG(PD3P,QD3P,Loadi);
Mat_H=FormH(Busnum,Volt,PG3P,PD3P,QG3P,QD3P,Y,UAngel,r,c,Yangle,Loadi); Mat_H=FormH(Busnum,Volt,PG3P,PD3P,QG3P,QD3P,Y,UAngel,r,c,Yangle,Loadi);
Ly=Mat_H; Ly=Mat_H;
Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P); Lz=FormLz(Mat_G,Init_L,Loadi,mPD3P,mQD3P);
Lw=FormLw(Mat_G,Init_U,Loadi,rPD3P,rQD3P); Lw=FormLw(Mat_G,Init_U,Loadi,mPD3P,mQD3P);
Lx=FormLx(deltF,deltH,Init_Y,deltG,Init_Z,Init_W); Lx=FormLx(deltF,deltH,Init_Y,deltG,Init_Z,Init_W);
% YY=FormYY(Lul,Lz,Ly,Luu,Lw,Lx); % YY=FormYY(Lul,Lz,Ly,Luu,Lw,Lx);
%% 开始解方程 %% 开始解方程