Hi Paulo,
Please see the code below:
mpc = case33bw;
c = mpc.branch(:,11); % 0,1 branch status
mpc2 = mpc;
mpc2.branch = mpc.branch(c~=0,:);
mpopt = mpoption('pf.alg', 'PQSUM');
result = radial_pf(mpc2,mpopt);
I remember that in order to use radial_pf, the network should be radial,
which means that even if there are open branches which create a loop, they
should be removed.
Hope this helps.
Kind regards,
Ilias
*Dr Ilias Sarantakos*
Research Associate
School of Engineering
Newcastle University
Στις Παρ, 23 Οκτ 2020 στις 6:48 π.μ., ο/η Paulo Moisés <
[email protected]> έγραψε:
> Hi,
>
> I have a difficulty in running a radial power flow by using the radia_pf
> function in the matpower.
> Can you help me in this regard?
>
> Thanks
> Paulo
>