function deltdeltF=func_deltdeltF(wVolt,wPD,wQD,ContrlCount,Loadi) %ContrlCount=size(PVi,1)+size(PGi,1)+size(wD,1)+Busnum*2; %P,Q,Volt theta这些控制变量数 C=[wPD' wQD']; CVolt=wVolt; sizeC=length(C); diagC=sparse(1:sizeC,1:sizeC,C,sizeC,sizeC); deltdeltF=[ diagC*2,sparse(sizeC,ContrlCount-sizeC); %sparse(ContrlCount-sizeC,ContrlCount); diag(CVolt)*2,sparse(length(CVolt),ContrlCount-length(CVolt)); sparse(length(CVolt),ContrlCount); ]; end