runcpf returns a MATPOWER case struct that has the loading/genration and the
voltages at its last continuation step. By using the option 'cpf_stop_at_nose',
the continuation power flow will stop when it is near the nose point. You can
then extract the voltages and power from the MATPOWER case struct.
For example,
define_constants;
mpopt = mpoption('out.all',0,'verbose',2);
mpopt = mpoption('cpf.stop_at','nose','cpf.step',0.2);
mpopt = mpoption(mpopt,'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(:,[PD QG]) = mpcb.bus(:,[PD QG])*2.5
results = runcpf(mpcb, mpct, mpopt);
Vm = results.bus(:,8); % Voltage magnitude at the nose point for all the buses.
Other values can be extracted from the case struct by accessing the appropriate
fields and their corresponding columns. The data format for the case struct is
defined in caseformat.m.
Shri
From: arun s nair <arunsnair...@gmail.com<mailto:arunsnair...@gmail.com>>
Reply-To: MATPOWER discussion forum
<matpowe...@list.cornell.edu<mailto:matpowe...@list.cornell.edu>>
Date: Fri, 16 Jan 2015 12:06:50 -0700
To: MATPOWER discussion forum
<matpowe...@list.cornell.edu<mailto:matpowe...@list.cornell.edu>>
Subject: Re: Matpower- Continuation Power Flow
Hi
Thank you Abhyanker for your input, the example in the manual worked.
I am really new to Matpower and I am learning it from scratch. My intention is
to get the critical value of power and voltage values(Nose point of the PV
curve) of a particular bus, Can anyone please tell me is it possible to get
the values using Matpower using continuation power flow ?
I am trying to modify the example code from the manual (I am working on IEEE 39
bus system)
define_constants;
mpopt = mpoption('out.all',0,'verbose',2);
mpopt = mpoption('cpf.stop_at','nose','cpf.step',0.2);
mpopt = mpoption(mpopt,'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(:,[PD QG]) = mpcb.bus(:,[PD QG])*2.5
results = runcpf(mpcb, mpct, mpopt);
results.cpf
Thanking in advance
Regards
Nair
On Thu, Jan 15, 2015 at 3:25 PM, Abhyankar, Shrirang G.
<abhy...@mcs.anl.gov<mailto:abhy...@mcs.anl.gov>> wrote:
Arun,
You have the same generation/load for both the base (mpcb) and target (mpct)
cases. What are you seeing is that runcpf is solving the base case over and
over again since there is zero power transfer. Use different generation/load
for mpcb and mpct and run the continuation power flow again. See the MATPOWER
User's Manual<http://www.pserc.cornell.edu//matpower/manual.pdf> for a CPF
example on a test 9-bus case.
Shri
From: arun s nair <arunsnair...@gmail.com<mailto:arunsnair...@gmail.com>>
Reply-To: MATPOWER discussion forum
<matpowe...@list.cornell.edu<mailto:matpowe...@list.cornell.edu>>
Date: Thu, 15 Jan 2015 13:02:57 -0700
To: <MATPOWER-L@cornell.edu<mailto:MATPOWER-L@cornell.edu>>
Subject: Matpower- Continuation Power Flow
Hi All
I am trying to use matpower to evaluate IEEE 39 bus system and I tried to run
an example code for runcpf (Continuation Power Flow)
define_constants;
mpopt = mpoption('out.all',0,'verbose',2);
mpopt = mpoption('cpf.stop_at','nose','cpf.step',0.2);
mpopt = mpoption(mpopt,'cpf.plot.level',2);
mpcb = loadcase(t_case9_pfv2); % load base case
mpct = mpcb; % set up target case with
mpct.gen(:, [PG QG]) = mpcb.gen(:, [PG QG])
mpct.bus(2,[PD]) = mpcb.bus(2, [PD])
results = runcpf(mpcb, mpct, mpopt);
results.cpf
The program is keep on running and is not converging. Can anybody please help
on this.
If anyone can send me a running example for runcpf of matpower , it will be of
great help and if anybody has worked on IEEE 39 bus using runcpf, please let me
know
Thanking you all in advance.
Regards
Arun Nair
