Assuming 101 is the row number in the mpc.bus array for bus 101:

define_constants;
mpopt = mpoption('out.all', 1, 'verbose', 2);
mpopt = mpoption(mpopt, 'cpf.stop_at', 'FULL', 'cpf.step', 
0.2,’cpf.plot.bus’,101);
mpopt = mpoption(mpopt, 'cpf.plot.level', 2);
mpcb = loadcase('Case_04'); % load base case
mpct = mpcb; % set up target case with
mpct.bus(101, [PD QD]) = mpcb.bus(101, [PD QD]) * 2.5 ; % and increased load
results = runcpf(mpcb, mpct, mpopt);
Shri


From: "[email protected]<mailto:[email protected]>" 
<[email protected]<mailto:[email protected]>>
Reply-To: MATPOWER discussion forum 
<[email protected]<mailto:[email protected]>>
Date: Saturday, May 13, 2017 at 10:31 PM
To: "[email protected]<mailto:[email protected]>" 
<[email protected]<mailto:[email protected]>>
Subject: Load increase in one busbar using runcpf

Dear all

Can anyone tell me what command I can use to increase a load in a specific 
busbar (let’s say bus 101) and get the P-V curve on that busbar please ?

I have these commands to increase overall load 2.5 times increased:
define_constants;
mpopt = mpoption('out.all', 1, 'verbose', 2);
mpopt = mpoption(mpopt, 'cpf.stop_at', 'FULL', 'cpf.step', 0.2);
mpopt = mpoption(mpopt, 'cpf.plot.level', 2);
mpcb = loadcase('Case_04'); % load base case
mpct = mpcb; % set up target case with
mpct.bus(:, [PD QD]) = mpcb.bus(:, [PD QD]) * 2.5 ; % and increased load
results = runcpf(mpcb, mpct, mpopt);

Best Regards
Kawsar

Reply via email to