Thank toy for your response. Here is the code. define_constants;mpcb = loadcase('case118');mpct = mpcb;mpct.bus(:, [PD QD]) = mpct.bus(:, [PD QD])*1.0;mpci = ext2int(mpct);Jacc = makeJac(mpci);Jac=full(Jacc); Jac_inv= inv(Jac);% inverse of Jsen_vpp=Jac_inv(1:117,118:181); % sensitivity of voltages to real power injection On Wednesday, September 6, 2023 at 10:46:30 AM EDT, Ray Daniel Zimmerman <r...@cornell.edu> wrote: Can you describe exactly what you are doing in a bit more detail (or maybe provide a few lines of code to illustrate)? I’m not sure where you are getting a matrix that is 117x64. Ray
On Sep 3, 2023, at 7:10 AM, Kaled Zareer <kaled...@yahoo.com> wrote: Dear MATPOWER team, I used inv(Jac) to obtain the voltage sensitivities of 118-bus. The size of "dv/dp" matrix is 117x64. I extracted the sensitivities related to load buses (i.e. 64x64 matrix). As it is known, the diagonal elements must be greater than the off-diagonal elements. However, the results do not coincide with this fact.Your comments are highly appreciated. Khaled