29 lines
1.1 KiB
Matlab
29 lines
1.1 KiB
Matlab
function deltdeltF=func_deltdeltF(busNum,fsY11,Loadi,wV1r,wV1i,wI1r,wI1i,V1measurement,V1r,V1i)
|
|
% wV1r_=(1)^2*ones(length(wV1r),1);
|
|
% wV1i_=(1)^2*ones(length(wV1i),1);
|
|
% wI1r_=(1)^2*ones(length(wI1r),1);
|
|
% wI1i_=(1)^2*ones(length(wI1i),1);
|
|
% infWV=find(wV1r==Inf);
|
|
% infWI=find(wI1r==Inf);
|
|
% wV1r=wV1r_;
|
|
% wV1i=wV1i_;
|
|
% wI1r=wI1r_;
|
|
% wI1i=wI1i_;
|
|
% wV1r(infWV)=Inf;
|
|
% wV1i(infWV)=Inf;
|
|
% wI1r(infWI)=Inf;
|
|
% wI1i(infWI)=Inf;
|
|
|
|
W2=wV1r.*wV1r;
|
|
t1=-2*V1r(Loadi).*V1r(Loadi)+ abs(V1measurement).^2;
|
|
t2=-(V1r(Loadi).^2+ V1i(Loadi).^2);
|
|
ttt=-4*(t1+t2 )./W2;
|
|
t3=-4*V1r(Loadi).*( -2*V1i(Loadi))./W2;
|
|
deltdeltF=[
|
|
%zeros(busNum*2,busNum*2+length(Loadi)*2);
|
|
sparse(Loadi,Loadi,ttt,busNum,busNum),sparse(Loadi,Loadi+busNum,t3,busNum,busNum+length(Loadi)*2);
|
|
sparse(Loadi,Loadi,-4*V1i(Loadi).*( -2*V1r(Loadi))./W2,busNum,busNum),sparse(Loadi,Loadi+busNum,-4*( -2*V1i(Loadi).*V1i(Loadi)+ abs(V1measurement).^2-(V1r(Loadi).^2+ V1i(Loadi).^2) )./W2,busNum,busNum+length(Loadi)*2);
|
|
sparse(length(Loadi)*2,busNum*2),2*eye(length(Loadi)*2)*diag(1./[wI1r.^2;wI1i.^2]);
|
|
%sparse(Loadi,busNum*2+Loadi,2,busNum,busNum*2+length(Loadi)*2);
|
|
];
|
|
end |