I tried to look for 'makeJac' function as told and search in the runpf.m file, however, the only 'make' that I found were 'makeBdc' , 'makeSbus', 'makeYbus', and 'makeB'.
So, is the problem only in my MATPOWER file (I need to download new matpower files again) since I can't found the 'makeJac()' function anywhere in my runpf.m file? On Tue, Jul 11, 2017 at 3:29 PM, Bone, Gorazd <[email protected]> wrote: > 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 + ,*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? >> >> >> >> >
