threephasese-onlyenquationconstraint-gn/jacobian_M.m

function [Jacob]=jacobian_M3(Busnum,Volt,Y,Angle,UAngel,r,c)   
%**************************************************************************
%  <EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD> : <EFBFBD>Ӻ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>γ<EFBFBD><EFBFBD>ſɱȾ<EFBFBD><EFBFBD><EFBFBD>Jacobian
%  <EFBFBD><EFBFBD>    <EFBFBD>ߣ<EFBFBD> 
%  <EFBFBD><EFBFBD><EFBFBD><EFBFBD>ʱ<EFBFBD>䣺2010.12
%**************************************************************************
%%<EFBFBD><EFBFBD><EFBFBD><EFBFBD>ͼ<EFBFBD><EFBFBD><EFBFBD><EFBFBD>6¥<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>д
%% <EFBFBD>ֱ<EFBFBD><EFBFBD><EFBFBD><EFBFBD>ſ˱Ⱦ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>H,L,N,J<EFBFBD><EFBFBD><EFBFBD>й<EFBFBD><EFBFBD>޹<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>P,Q
AngleIJ=UAngel(r)-UAngel(c)-Angle;
mat_AngleIJ=sparse(r,c,AngleIJ,Busnum*3,Busnum*3);
H=diag(Volt)*Y.*sin(mat_AngleIJ)*diag(Volt)-diag(Y.*sin(mat_AngleIJ)*Volt)*diag(Volt);
N=-diag(Volt)*Y.*cos(mat_AngleIJ)*diag(Volt)+diag(Y.*cos(mat_AngleIJ)*Volt)*diag(Volt);
J=diag(Y.*cos(mat_AngleIJ)*Volt)+Y.*cos(mat_AngleIJ)*diag(Volt);
L=diag(Y.*sin(mat_AngleIJ)*Volt)+Y.*sin(mat_AngleIJ)*diag(Volt);
t1=[J,H;
    L,N;
];
Jacob=t1;
end