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 


  

Reply via email to