Dr. Zimmerman

I applied the change you have mentioned and the result is the following:


MATPOWER Version 6.0b1, 01-Jun-2016 -- AC Optimal Power Flow
MATLAB Interior Point Solver -- MIPS, Version 1.2.1, 01-Jun-2016
 (using built-in linear solver)
Numerically Failed

Did not converge in 9 iterations.

>>>>>  Did NOT converge (0.11 seconds)  <<<<<




2016-11-23 19:54 GMT+02:00 Ray Zimmerman <[email protected]>:

> I think you want to change the first column of your gencost entries to 2
> (polynomial) instead of 1 (piecewise linear). Hopefully, that fixes your
> problem and it will solve.
>
>    Ray
>
> > On Nov 23, 2016, at 10:30 AM, Alexandros Gesios <
> [email protected]> wrote:
> >
> > Hello to everyone
> >
> > I am new to MATPOWER and I am trying to operate an OPF analysis in a MV
> grid constisted of two feeders.
> > The first feeder has 10 inductive loads and the second feeder has 10
> photovoltaic units and their are connected to the grid by the slack bus.
> > What I want to achieve is loss minimization by allocating the reactive
> power among the photovoltaic units in a optimal way. The nominal active
> power of the PV units is 10.6 MW and the load is 6 MW so the extra active
> power is being absorbed by the HV/MV Grid.
> > Also I want the slack bus to operate at constant 1.05 pu voltage so in
> the mpc.bus slack bus has voltage limits Vmax=Vmin=1.05.
> > Also I want the PV units to operate in nominal power so in mpc.gen their
> active power limits are Pmax=Pmin=Pnom.
> >
> > Finally the mpc.gencost matrix constists of 11 rows and every row has
> these elements     1    0       0       2       2   4; because I read in
> the archive that in order to achieve loss minimization you need the same
> linear costs in the gencost matrix.
> >
> > When I run runopf(mycase) I take the following errors
> >
> > Index exceeds matrix dimensions.
> >
> > Error in opf_setup (line 84)
> >   x1 = mpc.gencost(pwl1, COST+2);
> >
> > Error in opf (line 198)
> > om = opf_setup(mpc, mpopt);
> >
> > Error in runopf (line 75)
> > [r, success] = opf(casedata, mpopt);
> >
> > Also, I would like to ask whether my problem can be formulated in
> MATPOWER.
> >
> >
> > Thank you in advance
> >
> >
> > Alexandros Gesios
>
>
>
>

Reply via email to