83 lines
3.2 KiB
Matlab
83 lines
3.2 KiB
Matlab
function deltG=func_deltG(Busnum,PVi,PGi,Loadi,PD,QD)
|
|
sizeLoadi=size(Loadi,1);
|
|
%dg3 PD-M*b-t-PD0
|
|
dg3_dPD=sparse(1:size(Loadi,1),1:size(Loadi,1),ones(size(Loadi,1),1),size(Loadi,1),size(Loadi,1));
|
|
%dg32 PD+M*B+t+PD0
|
|
dg32_dPD=dg3_dPD;
|
|
%dg4 QD-M*B-t-QD0
|
|
dg4_dPD=sparse(size(Loadi,1),length(Loadi));
|
|
%dg42 QD+M*B+t+QD0
|
|
dg42_dPD=sparse(size(Loadi,1),length(Loadi));
|
|
dg5_dPD=sparse(size(Loadi,1),Busnum);
|
|
dg6_dPD=dg5_dPD;
|
|
%dg7 ÊÇ VbiµÄÔ¼Êø
|
|
dg7_dPD=sparse(sizeLoadi,Busnum);
|
|
dpdbi_dPD=sparse(sizeLoadi,sizeLoadi);
|
|
dqdbi_dPD=sparse(sizeLoadi,sizeLoadi);
|
|
%%
|
|
dg3_dQD=sparse(length(Loadi),length(Loadi));
|
|
dg32_dQD=sparse(sizeLoadi,sizeLoadi);
|
|
dg4_dQD=sparse(1:size(Loadi,1),1:size(Loadi,1),ones(size(Loadi,1),1),size(Loadi,1),size(Loadi,1));
|
|
dg42_dQD=sparse(1:size(Loadi,1),1:size(Loadi,1),ones(size(Loadi,1),1),size(Loadi,1),size(Loadi,1));
|
|
dg5_dQD=sparse(size(Loadi,1),Busnum);
|
|
dg6_dQD=dg5_dQD;
|
|
dg7_dQD=sparse(sizeLoadi,Busnum);
|
|
dpdbi_dQD=sparse(sizeLoadi,sizeLoadi);
|
|
dqdbi_dQD=sparse(sizeLoadi,sizeLoadi);
|
|
%%
|
|
dg3_dx=sparse(2*Busnum,sizeLoadi);
|
|
dg32_dx=sparse(2*Busnum,sizeLoadi);
|
|
dg4_dx=sparse(2*Busnum,length(Loadi));
|
|
dg42_dx=sparse(2*Busnum,sizeLoadi);
|
|
dg5_dx=[sparse(1:Busnum,1:Busnum,ones(Busnum,1),Busnum,Busnum);
|
|
sparse(Busnum,Busnum);
|
|
];
|
|
dg6_dx=dg5_dx;
|
|
dg7_dx=sparse(2*Busnum,Busnum);
|
|
dpdbi_dx=sparse(2*Busnum,sizeLoadi);
|
|
dqdbi_dx=sparse(2*Busnum,sizeLoadi);
|
|
%%
|
|
dg3_dvbi=sparse(Busnum,sizeLoadi);
|
|
dg32_dvbi=sparse(Busnum,sizeLoadi);
|
|
dg4_dvbi=sparse(Busnum,length(Loadi));
|
|
dg42_dvbi=sparse(Busnum,sizeLoadi);
|
|
dg5_dvbi=sparse(-eye(Busnum,Busnum));
|
|
dg6_dvbi=sparse(eye(Busnum,Busnum));
|
|
dg7_dvbi=sparse(eye(Busnum,Busnum));
|
|
dpdbi_dvbi=sparse(Busnum,sizeLoadi);
|
|
dqdbi_dvbi=sparse(Busnum,sizeLoadi);
|
|
%%
|
|
dg3_dpdbi=sparse(-eye(sizeLoadi,sizeLoadi));
|
|
dg32_dpdbi=sparse(eye(sizeLoadi,sizeLoadi));
|
|
dg4_dpdbi=sparse(sizeLoadi,sizeLoadi);
|
|
dg42_dpdbi=sparse(sizeLoadi,sizeLoadi);
|
|
dg5_dpdbi=sparse(sizeLoadi,Busnum);
|
|
dg6_dpdbi=sparse(sizeLoadi,Busnum);
|
|
dg7_dpdbi=sparse(sizeLoadi,Busnum);
|
|
dpdbi_dpdbi=sparse(eye(sizeLoadi,sizeLoadi));
|
|
dqdbi_dpdbi=sparse(sizeLoadi,sizeLoadi);
|
|
%%
|
|
dg3_dqdbi=sparse(sizeLoadi,sizeLoadi);
|
|
dg32_dqdbi=sparse(sizeLoadi,sizeLoadi);
|
|
dg4_dqdbi=sparse(-eye(sizeLoadi,sizeLoadi));
|
|
dg42_dqdbi=sparse(eye(sizeLoadi,sizeLoadi));
|
|
dg5_dqdbi=sparse(sizeLoadi,Busnum);
|
|
dg6_dqdbi=sparse(sizeLoadi,Busnum);
|
|
dg7_dqdbi=sparse(sizeLoadi,Busnum);
|
|
dpdbi_dqdbi=sparse(sizeLoadi,sizeLoadi);
|
|
dqdbi_dqdbi=sparse(eye(sizeLoadi,sizeLoadi));
|
|
%%
|
|
deltG=[dg3_dPD,dg32_dPD,dg4_dPD,dg42_dPD,dg5_dPD,dg6_dPD,dg7_dPD,dpdbi_dPD,dqdbi_dPD;
|
|
dg3_dQD,dg32_dQD,dg4_dQD,dg42_dQD,dg5_dQD,dg6_dQD,dg7_dQD,dpdbi_dQD,dqdbi_dQD;
|
|
dg3_dx,dg32_dx,dg4_dx,dg42_dx,dg5_dx,dg6_dx,dg7_dx,dpdbi_dx,dqdbi_dx;
|
|
dg3_dvbi,dg32_dvbi,dg4_dvbi,dg42_dvbi,dg5_dvbi,dg6_dvbi,dg7_dvbi,dpdbi_dvbi,dqdbi_dvbi;
|
|
dg3_dpdbi,dg32_dpdbi,dg4_dpdbi,dg42_dpdbi,dg5_dpdbi,dg6_dpdbi,dg7_dpdbi,dpdbi_dpdbi,dqdbi_dpdbi;
|
|
dg3_dqdbi,dg32_dqdbi,dg4_dqdbi,dg42_dqdbi,dg5_dqdbi,dg6_dqdbi,dg7_dqdbi,dpdbi_dqdbi,dqdbi_dqdbi;
|
|
];
|
|
% deltG=[dg3_dPD,dg32_dPD,dg4_dPD,dg5_dPD,dg6_dPD,dg7_dPD;
|
|
% dg3_dQD,dg32_dQD,dg4_dQD,dg5_dQD,dg6_dQD,dg7_dQD;
|
|
% dg3_dx,dg32_dx,dg4_dx,dg5_dx,dg6_dx,dg7_dx;
|
|
% dg3_dvbi,dg32_dvbi,dg4_dvbi,dg5_dvbi,dg6_dvbi,dg7_dvbi;
|
|
% dg3_dpdbi,dg32_dpdbi,dg4_dpdbi,dg5_dpdbi,dg6_dpdbi,dg7_dpdbi;
|
|
% ];
|
|
end |