function deltdeltF=func_deltdeltF(Busnum,GenC,PVi,PGi) ContrlCount=size(PVi,1)+size(PGi,1)+Busnum*2; %P,Q,Volt theta这些控制变量数 deltdeltF=[diag(GenC(:,2))*2,zeros(size(GenC,1),ContrlCount-size(GenC,1)); zeros(ContrlCount-size(GenC,1),ContrlCount); ]; %#ok sizeGenC=size(GenC(:,2),1); diagC=sparse(1:sizeGenC,1:sizeGenC,GenC(:,2),sizeGenC,sizeGenC); deltdeltF=[ diagC*2,sparse(sizeGenC,ContrlCount-sizeGenC); sparse(ContrlCount-sizeGenC,ContrlCount); ]; end