Thank you very much for the explanation. Firstly I am sorry, actually in my previous email, I was trying to make changes on the Jacobian Matrix, not the admittance matrices.
Therefore, 1. How to find the Jacobian Matrix in order to modify the set of linear equation so that it is relevant to my work? 2. And, from below equations PLosses = sum[αij (PiPj + QiQj) + βij(PiPj - QiQj)] (1) where,αij = rij/ ( ViVj) cos (ẟi - ẟj) and βij = rij/ ( ViVj) sin (ẟi - ẟj ) rij + ,*j*xij= Zij are the ijth element of [Zbus] matrix with [Zbus]=[Ybus] -1 The sensitivity factor of real power loss with respect to real power injection from DG is given by αi= (partial derivative of PLosses/ partial derivative of Pi) = 2 (sum) (α ij Pj - βijQj) where can i find the equation in the runpf.m file? I'm still with less experience using the MATPOWER code. Please help me. Thank you. On Thu, Jul 6, 2017 at 8:52 PM, Ray Zimmerman <[email protected]> wrote: > It is not clear what you are trying to accomplish, but MATPOWER code is > just like any other Matlab code. You can duplicate and rename runpf.m and > modify it to your heart’s content. However, multiplying an entire bus or > branch matrix by 2 is certainly not a valid thing to do. For one thing the > bus types and the branch from and to bus numbers will no longer be valid. > > Ray > > > On Jul 5, 2017, at 11:45 PM, Ameeratul Balqis <[email protected]> > wrote: > > > > Hello all, > > I would like to make a changes on the admittance matrices formation > whenever I run the code of runpf(edited).m. Thus, I would like to ask where > is the right place to edit on: > > > > Because when I've edit on > > TCL MERGE ERROR ( 07/06/2017 08:52:37 ): "invalid command name > "Ybus,"" > OutmailID: 121641754, List: 'matpower-l', MemberID: 77921808 > SCRIPT: "build admittance matrices > [Ybus, Yf, Yt] = > makeYbus(baseMVA, bus, branch); > > repeat = 1; > > as (example) > > >" > build admittance matrices > > [Ybus, Yf, Yt] = makeYbus(2*baseMVA, 2*bus, 2*branch); > > > > repeat = 1; > > > > it will shows that; > > > > Error in ==> runpf at 193 > > [Ybus, Yf, Yt] = makeYbus(2*baseMVA,2* bus, 2*branch); > > > > Is there any editing on the runpf.m can be made in order for me to apply > a changed runpf.m in my term so that I can applied it in the case study? > > > >
