Thanks Dr. Zimmerman, I was able to get it working on my system Thai
On Tue, Jun 30, 2015 at 5:11 AM, Ray Zimmerman <[email protected]> wrote: > It sounds like you want set your generator costs to zero and put a > user-defined quadratic cost on the deviations of the voltage magnitudes > from 1. Something like the following should do it for case30 … > > define_constants; > mpopt = mpoption('verbose', 2, 'opf.ac.solver', 'MIPS'); > mpc = loadcase('case30'); > nb = size(mpc.bus, 1); > ng = size(mpc.gen, 1); > mpc.N = sparse(1:nb, nb+(1:nb), 1, nb, 2*nb+2*ng); %% select Vm > rhat = 1; %% offset 1 p.u. > k = 0; > d = 2 > m = 1; > mpc.fparm = ones(nb, 1) * [d rhat k m]; > mpc.Cw = 10000 * ones(nb, 1); > mpc.gencost(:, COST:end) = 0; > r = runopf(mpc, mpopt); > > Hope this helps, > > Ray > > > > > On Jun 29, 2015, at 2:22 PM, Thai Tran <[email protected]> wrote: > > Hi All, > > I want to run an opf where the objective is to minimize voltage > variability (i.e keep the voltage across the feeder as close to 1 p.u as > possible) how should I structure the cost matrix? > > Thanks > Thai > > >
