2014-05-14 10:00:41 +08:00
|
|
|
function Lz=FormLz(Mat_G,Init_L,Busnum,PD0,QD0,Loadi)
|
2012-11-16 22:10:18 +08:00
|
|
|
VoltL=(0.9)*ones(1,Busnum);
|
2012-05-25 15:38:31 +08:00
|
|
|
PDL=PD0(Loadi);
|
2013-03-05 10:43:18 +08:00
|
|
|
PDL(PDL>0)=0.800*PDL(PDL>0);
|
|
|
|
|
PDL(PDL<0)=1.200*PDL(PDL<0);
|
2012-06-06 10:19:08 +08:00
|
|
|
PDL(PDL==0)=-0.400;
|
2014-05-13 16:48:46 +08:00
|
|
|
realPD=PD0(Loadi);
|
|
|
|
|
indPD=find(realPD>0);
|
|
|
|
|
PDL(indPD(3:12:end))=0.95*realPD(indPD(3:12:end));
|
|
|
|
|
PDL(indPD(9:12:end))=0.45*realPD(indPD(9:12:end));
|
2013-01-21 18:28:16 +08:00
|
|
|
QDL=QD0(Loadi);
|
2013-03-05 10:43:18 +08:00
|
|
|
QDL(QDL>0)=0.800*QDL(QDL>0);
|
|
|
|
|
QDL(QDL<0)=1.200*QDL(QDL<0);
|
2013-01-21 21:47:54 +08:00
|
|
|
QDL(QDL==0)=-0.400;
|
2014-05-13 16:48:46 +08:00
|
|
|
realQD=QD0(Loadi);
|
|
|
|
|
indQD=find(realQD>0);
|
|
|
|
|
QDL(indQD(3:12:end))=0.95*realQD(indQD(3:12:end));
|
|
|
|
|
QDL(indQD(9:12:end))=0.95*realQD(indQD(9:12:end));
|
2014-05-14 10:12:34 +08:00
|
|
|
t1=([PDL',PDL',QDL',-5000*ones(1,length(VoltL)),0*VoltL,0*ones(1,Busnum)])';
|
2012-05-22 11:33:21 +08:00
|
|
|
t2=Mat_G-Init_L'-t1;
|
|
|
|
|
Lz=t2;
|
|
|
|
|
end
|
