threephasese-onlyenquationconstraint-gn/func_deltF.m

function deltF=func_deltF(wVolt,wPD,wQD,PD0,PD,QD,QD0,Volt,mVolt,Busnum,Loadi)
% t3=2*wPD.*(PD-PD0);
% t4=2*wQD.*(QD-QD0);
% t5=2*wVolt.*(Volt-mVolt);
% deltF=[sparse(t3);
%     sparse(t4);
%     sparse(t5);
%     sparse(1*Busnum*3,1);
% ];
totalM=length(PD)+length(QD)+length(Volt);
deltF=[
    eye(length(PD)),sparse(length(PD),totalM-length(PD));
    sparse(length(QD),length(PD)),eye(length(PD)),sparse(length(QD),totalM-length(PD)-length(QD));
    sparse(length(Volt),length(PD)+length(QD)),eye(length(Volt));
    sparse(length(Volt),totalM);
]';
deltF=sparse(deltF);
end