18 lines
662 B
Mathematica
18 lines
662 B
Mathematica
|
function [Init_Y,PG,QG,Volt,UAngel,PD,QD]=Modification(Init_Y,deltX,deltY,PG,QG,Volt,UAngel,ContrlCount,Balance,Busnum,PD,QD,Loadi)
|
|||
|
|
Init_Y=Init_Y+deltY;
|
|||
|
|
t=deltX(1:size(Loadi,1)*2);
|
|||
|
|
PD=PD+t(1:length(Loadi));
|
|||
|
|
QD=QD+t(length(Loadi)+1:length(Loadi)*2);
|
|||
|
|
t=deltX(1:Busnum*3*2)';
|
|||
|
|
%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
|
|||
|
|
% balVolt=Volt( (Balance-1)*3+1 );
|
|||
|
|
Volt=Volt+t(1:Busnum*3)';
|
|||
|
|
Volt( (Balance-1)*3+1 )=1;
|
|||
|
|
Volt( (Balance-1)*3+2 )=1;
|
|||
|
|
Volt( (Balance-1)*3+3 )=1;
|
|||
|
|
t(Busnum*3+(Balance-1)*3+1)=0;
|
|||
|
|
t(Busnum*3+(Balance-1)*3+2)=0;
|
|||
|
|
t(Busnum*3+(Balance-1)*3+3)=0;
|
|||
|
|
UAngel=UAngel+t(Busnum*3+1:2*Busnum*3)';
|
|||
|
|
end
|