2012-05-25 15:38:31 +08:00
|
|
|
|
function [Init_Z,Init_L,Init_W,Init_U,Init_Y,PG,QG,Volt,UAngel,PD]=Modification(Init_Z,Init_L,Init_W,Init_U,Init_Y,deltZ,deltL,deltW,deltU,deltX,deltY,PG,QG,Volt,UAngel,PVi,ContrlCount,Balance,Busnum,PGi,PD,Loadi)
|
2012-05-22 11:33:21 +08:00
|
|
|
|
AlphaP=FormAlphaP(Init_L,deltL,Init_U,deltU);
|
2012-07-25 15:08:54 +08:00
|
|
|
|
%fprintf('AlphaP %f\n',full(AlphaP));
|
2012-05-22 11:33:21 +08:00
|
|
|
|
AlphaD=FormAlphaD(Init_Z,deltZ,Init_W,deltW);
|
2012-07-25 15:08:54 +08:00
|
|
|
|
%fprintf('AlphaD %f\n',full(AlphaD));
|
2012-05-22 11:33:21 +08:00
|
|
|
|
|
|
|
|
|
|
Init_Z=Init_Z+AlphaD*deltZ';
|
|
|
|
|
|
Init_L=Init_L+AlphaP*deltL';
|
|
|
|
|
|
Init_W=Init_W+AlphaD*deltW';
|
|
|
|
|
|
Init_U=Init_U+AlphaP*deltU';
|
|
|
|
|
|
Init_Y=Init_Y+AlphaD*deltY';
|
|
|
|
|
|
%PG(PGi)=PG(PGi)+deltX(1:size(PGi,1));
|
|
|
|
|
|
PG(PGi)=PG(PGi)+AlphaP*deltX(1:size(PGi,1));
|
|
|
|
|
|
%QG(PVi)=QG(PVi)+deltX(size(PGi,1)+1:size(PVi,1)+size(PGi,1) );
|
|
|
|
|
|
QG(PVi)=QG(PVi)+AlphaP*deltX(size(PGi,1)+1:size(PVi,1)+size(PGi,1) );
|
2012-05-25 15:38:31 +08:00
|
|
|
|
t=deltX(size(PVi,1)+size(PGi,1)+1:size(PVi,1)+size(PGi,1)+size(Loadi,1));
|
|
|
|
|
|
PD(Loadi)=PD(Loadi)+AlphaP*t;
|
|
|
|
|
|
t=deltX(size(PVi,1)+size(PGi,1)+size(Loadi,1)+1:ContrlCount)';
|
2012-05-22 11:33:21 +08:00
|
|
|
|
t(Busnum+Balance)=0;
|
|
|
|
|
|
%Volt=Volt+AlphaP*t(2:2:2*Busnum);<EFBFBD><EFBFBD>ʱ<EFBFBD><EFBFBD>һ<EFBFBD><EFBFBD>20111227
|
|
|
|
|
|
%UAngel=UAngel+AlphaP*t(1:2:2*Busnum);<EFBFBD><EFBFBD>ʱ<EFBFBD><EFBFBD>һ<EFBFBD><EFBFBD>20111227
|
2012-05-22 11:41:52 +08:00
|
|
|
|
balVolt=Volt(Balance);
|
2012-05-22 11:33:21 +08:00
|
|
|
|
Volt=Volt+AlphaP*t(1:Busnum);
|
2012-05-22 11:41:52 +08:00
|
|
|
|
Volt(Balance)=balVolt;
|
2012-05-22 11:33:21 +08:00
|
|
|
|
UAngel=UAngel+AlphaP*t(Busnum+1:2*Busnum);
|
|
|
|
|
|
end
|
