Dear Ray, Thank you for your response. For calculating DC state estimation I use AC formula but put all voltage to 1.
V = ones(size(mpc.bus_name)) [dSbus_dVm, dSbus_dVa] = dSbus_dV(Ybus, V) j = full(real(dSbus_dVa)) The problem is that my Jacobian size becomes to 14*14 but my measurement is 19 and state is 14 for IEEE 14 bus system. it means that I need 19*14 Jacobin. how can I change the measurement? Kind Regards, Naime On Sat, Apr 21, 2018 at 12:46 AM, Ray Zimmerman <r...@cornell.edu> wrote: > 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 <naimeahm...@gmail.com> 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 > > >