1. Look at function makeJac() 2. That looks like the question of power system sensitivity analysis. You may have to program it yourself.
regards, Gorazd ________________________________ Od: [email protected] [[email protected]] v imenu Ameeratul Balqis [[email protected]] Poslano: 11. julij 2017 4:39 Za: MATPOWER discussion forum Zadeva: Re: To make changes in runpf.m 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 + ,jxij= 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]<mailto:[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]<mailto:[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?
