I not sure I follow completely, but if you are talking about line flows, rather
than bus injections then you want to use Bf rather than Bbus from the return
values of makeBdc().
Ray
> On Apr 19, 2018, at 10:16 PM, naime ahmadi <[email protected]> wrote:
>
> Dear Ray,
> Thank you. I want to calculate this Jacobin.
> P(i,j)=B(i,j)*(theta(i)-theta(j))
> PP=[P(1,2);P(1,4);P(1,5);P(2,3);P(2,4);P(4,5)] (measurement vector)
> J=double(jacobian(PP,theta)) getting derivatives of power flow between line
> i and j in respect to all other theta
> J = J(:,2:end); omit first column as a slack bus.
>
> Kind Regards,
> Naime
