And how do you change the status of a branch.Assuming you want to take a branch out of service(tripping it)
Sent from my iPhone > On 31 Mar 2018, at 19:41, Carlos E Murillo-Sanchez > <ce.murillosanc...@gmail.com> wrote: > > Assuming that the system data is in structure 'mpc', and that the buses are > numbered consecutively starting from 1, this should give you the injection > currents: > > % Run power flow > results = runpf(mpc); > % Get column pointers > define_constants; > % Compute vector of complex bus voltages > V = results.bus(:, VM) .* exp(sqrt(-1)*results.bus(:, VA)*pi/180); > % Compute nodal admittance matrix > Ybus = makeYbus(results); > % Compute vector of complex bus injection currents > Ibus = Ybus * V; > > Carlos. > > Mostafa Mohammadpourfard wrote: >> Thank you, Carlos. May I ask you to write the code here? I am not a power >> system student. Thanks >> >> >>> On Sat, Mar 31, 2018 at 9:41 AM, Carlos E Murillo-Sanchez >>> <ce.murillosanc...@gmail.com> wrote: >>> You have the bus voltages; now you just need to multiply the system's nodal >>> admittance matrix by the complex voltage vector to get the complex currents >>> vector. See the makeYbus function in MATTPOWER. >>> >>> Carlos. >>> >>> Arkan Arkan wrote: >>>> Thank you, Ilias. But, unfortunately, there is not any direct information >>>> about current in the Manpower like you mentioned and we should calculate >>>> based on the information obtained from runpf. >>>> >>>> On Fri, Mar 30, 2018 at 9:06 AM, Ilias Sarantakos <il.saranta...@gmail.com >>>> <mailto:il.saranta...@gmail.com>> wrote: >>>> >>>> Hi Arkan, >>>> >>>> If you type ''result = runpf(mpc)'', then you will get the power >>>> flow results of your network, stored in the ''result'' variable. >>>> Now that you have the voltages (magnitudes and angles) at all >>>> buses, you can calculate the current at each branch. I think it >>>> would be helpful to check Appendix B in MATPOWER manual, which can >>>> guide you how to access the results of your network. For example, >>>> ''result.bus(1,8)'' and ''result.bus(1,9)'', contain the magnitude >>>> and angle of voltage at bus 1. >>>> >>>> Hope this helps. >>>> >>>> Kind regards, >>>> >>>> Ilias Sarantakos >>>> >>>> 2018-03-29 2:25 GMT+01:00 Arkan Arkan <arkan.m2...@gmail.com >>>> <mailto:arkan.m2...@gmail.com>>: >>>> >>>> Hi everyone, >>>> >>>> I am wondering how can I get the magnitude and the phase angle >>>> of the current in the branches of any case file in the Matpower. >>>> >>>> Thanks for your help. >>>> >>>> >>>> >>> >>> >> >