Try using a smaller continuation step size. mpopt = mpoption(mpopt, 'cpf.stop_at', 'full', 'cpf.step', 0.05);
Shri > On Feb 12, 2015, at 4:03 AM, "Natakorn Thasnas" <[email protected]> > wrote: > > Hi All > > I am trying to run a continuation power flow with the following code > > define_constants; > mpopt = mpoption('out.all', 0, 'verbose', 2); > mpopt = mpoption(mpopt, 'cpf.stop_at', 'full', 'cpf.step', 0.2); > mpopt = mpoption(mpopt, 'cpf.plot.level', 2 ,'cpf.plot.bus', 14); > mpcb = loadcase('case14result'); % load base case > mpct = mpcb; % set up target case with > mpct.gen(:, [PG QG]) = mpcb.gen(:, [PG QG])*2.5; > mpct.bus(:,[PD QG]) = mpcb.bus(:,[PD QG])*2.5; > results = runcpf(mpcb, mpct, mpopt); > > After running it, I got the warning message "corrector did not converge in > 1000 iterations" and the PV curve is not full nose curve as showing in > attached file (c14.jpg) > > <c14.jpg> > > Then I increased "max_it" ​from 1000 to 10000. I got the warning messages as > following > > Warning: Matrix is close to singular or badly scaled. > Results may be inaccurate. RCOND = 1.944067e-018. > > In cpf_corrector at 143 > In runcpf at 268 > In cpfrun at 55 > step 48 : lambda = -7.770, corrector did not converge in 10000 iterations > > > Could you please suggest me how to fix it to get "the full nose curve"? > > Thank you in advance. > > Best Regards, > > Natakorn Thasnas > Ph.D Student > Khonkaen University, Thailand
