thank you so much llias On Mon, Apr 2, 2018 at 11:06 AM, Ilias Sarantakos <[email protected]> wrote:
> Hi Winnie, > > For example, ''mpc.branch(1,11) = 0'', would take branch 1 out of service. > For more information, you can check Appendix B in MATPOWER manual (page > 124, Branch Data). > > Hope this helps. > > Kind regards, > > Ilias > > 2018-03-31 21:16 GMT+01:00 Winnie apiyo <[email protected]>: > >> 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 < >> [email protected]> 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 < >> <[email protected]>[email protected]> 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 < >>>> <[email protected]>[email protected] <mailto: >>>> [email protected]>> 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 < <[email protected]> >>>> [email protected] >>>> <mailto:[email protected]>>: >>>> >>>> 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. >>>> >>>> >>>> >>>> >>> >>> >> >> >
