Re: Question regarding scaling the loads and generation in MATPower
I am extremely sorry about this. Thank you!
On Mon, Apr 13, 2015 at 4:48 PM, Ray Zimmerman wrote:
> I almost started to become superstitious …
>
> … then I noticed a typo in the last line before the call to runpf().
> Left-hand side has QD, right-hand side has PD.
>
>Ray
>
>
> On Apr 13, 2015, at 3:58 PM, Shruti Rao wrote:
>
> Hi all,
>
> I was trying to scale the loads and generation in MATPower using the
> following lines:
>
> clc;clear all;
> define_constants;
> mpcb = loadcase('case118'); % load base case
> mpct = mpcb; % set up target case with
> mpct.bus(:, VM) = 1;
> mpct.bus(:, VA) = 0;
> mpct.gen(:, VG) = 1;
> mpopt = mpoption('pf.tol', 1e-3,'pf.nr.max_it',
> 100,'pf.enforce_q_lims',1,'out.all', 1, 'verbose', 3);
> mpct.gen(:,PG)=mpct.gen(:,PG)*1.82736;
> mpct.gen(:,QG)=mpct.gen(:,QG)*1.82736;
> mpct.bus(:, PD) = mpcb.bus(:, PD)*1.82736; % and increased load
> mpct.bus(:, QD) = mpcb.bus(:, PD)*1.82736; % and increased load
> results = runpf(mpct, mpopt);
>
> The above code said that the power flow did not converge at that scaled
> loading level. However, on changing the order of the scaling (with the rest
> of the code remaining same) as follows, I got a solution:
>
> mpct.bus(:,PD)=mpcb.bus(:,PD)*1.82736;
> mpct.gen(:,PG)=mpcb.gen(:,PG)*1.82736;
> mpct.bus(:,QD)=mpcb.bus(:,QD)*1.82736;
> mpct.gen(:,QG)=mpcb.gen(:,QG)*1.82736;
>
> I have attached my published MATLAB code with the results.
>
> Why should the two codes give different results? Is there a reason behind
> having to follow a particular order? Which order should we follow?
>
> Thank you!
> --
> Best Regards,
> Shruti Dwarkanath Rao
>
> Graduate Research Assistant
> School of Electrical, Computer and Energy Engineering
> Arizona State University
> Tempe, AZ, 85281
> 650 996 0116
>
>
>
>
--
Best Regards,
Shruti Dwarkanath Rao
Graduate Research Assistant
School of Electrical, Computer and Energy Engineering
Arizona State University
Tempe, AZ, 85281
650 996 0116
Re: Question regarding scaling the loads and generation in MATPower
I almost started to become superstitious …
… then I noticed a typo in the last line before the call to runpf(). Left-hand
side has QD, right-hand side has PD.
Ray
> On Apr 13, 2015, at 3:58 PM, Shruti Rao wrote:
>
> Hi all,
>
> I was trying to scale the loads and generation in MATPower using the
> following lines:
>
> clc;clear all;
> define_constants;
> mpcb = loadcase('case118'); % load base case
> mpct = mpcb; % set up target case with
> mpct.bus(:, VM) = 1;
> mpct.bus(:, VA) = 0;
> mpct.gen(:, VG) = 1;
> mpopt = mpoption('pf.tol', 1e-3,'pf.nr.max_it',
> 100,'pf.enforce_q_lims',1,'out.all', 1, 'verbose', 3);
> mpct.gen(:,PG)=mpct.gen(:,PG)*1.82736;
> mpct.gen(:,QG)=mpct.gen(:,QG)*1.82736;
> mpct.bus(:, PD) = mpcb.bus(:, PD)*1.82736; % and increased load
> mpct.bus(:, QD) = mpcb.bus(:, PD)*1.82736; % and increased load
> results = runpf(mpct, mpopt);
>
> The above code said that the power flow did not converge at that scaled
> loading level. However, on changing the order of the scaling (with the rest
> of the code remaining same) as follows, I got a solution:
>
> mpct.bus(:,PD)=mpcb.bus(:,PD)*1.82736;
> mpct.gen(:,PG)=mpcb.gen(:,PG)*1.82736;
> mpct.bus(:,QD)=mpcb.bus(:,QD)*1.82736;
> mpct.gen(:,QG)=mpcb.gen(:,QG)*1.82736;
>
> I have attached my published MATLAB code with the results.
>
> Why should the two codes give different results? Is there a reason behind
> having to follow a particular order? Which order should we follow?
>
> Thank you!
> --
> Best Regards,
> Shruti Dwarkanath Rao
>
> Graduate Research Assistant
> School of Electrical, Computer and Energy Engineering
> Arizona State University
> Tempe, AZ, 85281
> 650 996 0116
>
Question regarding scaling the loads and generation in MATPower
Hi all,
I was trying to scale the loads and generation in MATPower using the
following lines:
clc;clear all;
define_constants;
mpcb = loadcase('case118'); % load base case
mpct = mpcb; % set up target case with
mpct.bus(:, VM) = 1;
mpct.bus(:, VA) = 0;
mpct.gen(:, VG) = 1;
mpopt = mpoption('pf.tol', 1e-3,'pf.nr.max_it',
100,'pf.enforce_q_lims',1,'out.all', 1, 'verbose', 3);
mpct.gen(:,PG)=mpct.gen(:,PG)*1.82736;
mpct.gen(:,QG)=mpct.gen(:,QG)*1.82736;
mpct.bus(:, PD) = mpcb.bus(:, PD)*1.82736; % and increased load
mpct.bus(:, QD) = mpcb.bus(:, PD)*1.82736; % and increased load
results = runpf(mpct, mpopt);
The above code said that the power flow did not converge at that scaled
loading level. However, on changing the order of the scaling (with the rest
of the code remaining same) as follows, I got a solution:
mpct.bus(:,PD)=mpcb.bus(:,PD)*1.82736;
mpct.gen(:,PG)=mpcb.gen(:,PG)*1.82736;
mpct.bus(:,QD)=mpcb.bus(:,QD)*1.82736;
mpct.gen(:,QG)=mpcb.gen(:,QG)*1.82736;
I have attached my published MATLAB code with the results.
Why should the two codes give different results? Is there a reason behind
having to follow a particular order? Which order should we follow?
Thank you!
--
Best Regards,
Shruti Dwarkanath Rao
Graduate Research Assistant
School of Electrical, Computer and Energy Engineering
Arizona State University
Tempe, AZ, 85281
650 996 0116
MATPower_Conflicting_solutions.pdf
Description: Adobe PDF document
