facat
/
distributionnetwork-power2current
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

修复了目标函数二次导数中一个系数的bug。不过这个系数的取值还需要考虑一下。 

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk 2014-12-04 16:33:34 +08:00
parent 1c2b53d2b5
commit 0b8015e9d0
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
func_deltdeltF.m
View File

@ -3,7 +3,7 @@ function deltdeltF=func_deltdeltF(busNum,fsY11,Loadi,wV1r,wV1i,wI1r,wI1i,V1measu
deltdeltF=[ deltdeltF=[
%zeros(busNum*2,busNum*2+length(Loadi)*2); %zeros(busNum*2,busNum*2+length(Loadi)*2);
sparse(Loadi,Loadi,-4*( -2*V1r(Loadi).*V1r(Loadi)+ abs(V1measurement).^2-(V1r(Loadi).^2+ V1i(Loadi).^2) )./[wV1r.*wV1r],busNum,busNum),sparse(Loadi,Loadi+busNum,-4*V1r(Loadi).*( -2*V1i(Loadi))./[wV1r.*wV1r],busNum,busNum+length(Loadi)*2); sparse(Loadi,Loadi,-4*( -2*V1r(Loadi).*V1r(Loadi)+ abs(V1measurement).^2-(V1r(Loadi).^2+ V1i(Loadi).^2) )./[wV1r.*wV1r],busNum,busNum),sparse(Loadi,Loadi+busNum,-4*V1r(Loadi).*( -2*V1i(Loadi))./[wV1r.*wV1r],busNum,busNum+length(Loadi)*2);
sparse(Loadi,Loadi,-4*V1i(Loadi).*( -2*V1r(Loadi))./[wV1i.*wV1i],busNum,busNum),sparse(Loadi,Loadi+busNum,-4*( -2*V1i(Loadi).*V1i(Loadi)+ abs(V1measurement).^2-(V1r(Loadi).^2+ V1i(Loadi).^2) )./[wV1r.*wV1r],busNum,busNum+length(Loadi)*2); sparse(Loadi,Loadi,-4*V1i(Loadi).*( -2*V1r(Loadi))./[wV1i.*wV1i],busNum,busNum),sparse(Loadi,Loadi+busNum,-4*( -2*V1i(Loadi).*V1i(Loadi)+ abs(V1measurement).^2-(V1r(Loadi).^2+ V1i(Loadi).^2) )./[wV1i.*wV1i],busNum,busNum+length(Loadi)*2);
zeros(length(Loadi)*2,busNum*2),2*eye(length(Loadi)*2)*diag(1./[wI1r.^2;wI1i.^2]); zeros(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); %sparse(Loadi,busNum*2+Loadi,2,busNum,busNum*2+length(Loadi)*2);
]; ];