Hi Arun, Since MATPOWER’s continuation power flow transfer is set up in a very general fashion as the difference between a base and target case, which may involve changes at multiple buses, there isn’t a single “power” number associated with each step. However, you can recover any power of interest from the results.cpf.lam_p and/or results.cpf.lam_c values and transfer (target - base) vector from equation (5.3) in the manual. The lambda values are used to interpolate between the base injections (lambda = 0) and the target injections (lambda = 1).
-- Ray Zimmerman Senior Research Associate B30 Warren Hall, Cornell University, Ithaca, NY 14853 USA phone: (607) 255-9645 > On Feb 3, 2015, at 12:47 PM, arun s nair <[email protected]> wrote: > > Hi All > > I am trying to use matpower to evaluate IEEE 39 bus system. I am able to > plot the PV curve of load buses and was also able to get the critical power > and voltage values. > > I would like to know if its possible to get the power values at each > iterations along with the voltage values ? > > > I am using runcpf, and used the below code. > > mpopt = mpoption('out.all',0,'verbose',2,'out.bus',1); > mpopt = mpoption(mpopt,'cpf.stop_at','nose','cpf.step',0.2); > mpopt = mpoption(mpopt,'cpf.plot.bus',3,'cpf.plot.level',2); > > mpcb = loadcase('case39'); % load base case > mpct = mpcb; % set up target case with > %mpct.gen(:,[PG QG]) = mpcb.gen(:,[PG QG])*2.5; > mpct.bus(3,PD) = mpcb.bus(3,PD)*2.5; > results = runcpf(mpcb, mpct, mpopt); > > CriticalPower = results.bus(3,3) > CriticalVoltage = results.bus(3,8) > > > Thanking in advance > > With Regards. > > Arun Nair
