The DC OPF in MATPOWER neglects losses. If you want to include losses, you will
have to use the AC OPF. Or you could estimate the losses from the DC OPF
solution by computing the average current in each branch and multiplying by the
corresponding square of resistance in the line.
Ray
> On Apr 23, 2018, at 11:47 AM, Mohammed Alhajri <[email protected]> wrote:
>
> Hello all
>
> how i can include the losses in my code?
>
> my code is
>
> ++++++++++++++++++++++++++++++++++++++++++++++
> clear all;
> clc;
> define_constants;
> factor = [0.827315297246503; 0.870524914020136; 0.89225313150181;
> 0.874238893371013; 0.841447907563717; 0.84108393381471;
> 0.79110106954204; 0.797588753768103; 0.801421029392904;
> 0.822018440318329; 0.864051234184459; 0.876248439711015;
> 0.888800014120353; 0.922571691519356; 0.978924730305068; 1;
> 0.975195764982014; 0.903702013699042; 0.8424313185939;
> 0.786194728460529; 0.846555864529582; 0.851944597303249;
> 0.834213349185084; 0.84694392020701;];
> mpc0 = loadcase('OETC_OPF');
> nt = length(factor);
> diary R;
> for t = 1:nt
> mpc(t) = mpc0;
> [mpc(t).bus] = ...
> scale_load(factor(t), mpc(t).bus);
> results(t) = rundcopf(mpc(t));
> end
>
> ++++++++++++++++++++++++++++++++++++++++++++++
>
> but it does not calculate the loses!
>
> On 23 April 2018 at 18:09, Mohammed Alhajri <[email protected]
> <mailto:[email protected]>> wrote:
> thank you doctor
>
> بتاريخ ٢٠١٨/٠٤/٢٣ ٦:٠٧ م، كتب "Ray Zimmerman" <[email protected]
> <mailto:[email protected]>>:
> The function most_summary() is only for summarizing the output of a call to
> most(). It does not apply to runopf().
>
> Ray
>
>> On Apr 21, 2018, at 6:27 AM, Mohammed Alhajri <[email protected]
>> <mailto:[email protected]>> wrote:
>>
>> thank you, i did it and it works
>>
>> ++++++++++++++++++++++++++++++++++++++++++++
>>
>> clear all;
>> clc;
>> define_constants;
>> factor = [0.827315297246503; 0.870524914020136; 0.89225313150181;
>> 0.874238893371013; 0.841447907563717; 0.84108393381471;
>> 0.79110106954204; 0.797588753768103; 0.801421029392904;
>> 0.822018440318329; 0.864051234184459; 0.876248439711015;
>> 0.888800014120353; 0.922571691519356; 0.978924730305068; 1;
>> 0.975195764982014; 0.903702013699042; 0.8424313185939;
>> 0.786194728460529; 0.846555864529582; 0.851944597303249;
>> 0.834213349185084; 0.84694392020701;];
>> mpc0 = loadcase('OETC_OPF');
>> nt = length(factor);
>> for t = 1:nt
>> mpc(t) = mpc0;
>> [mpc(t).bus] = ...
>> scale_load(factor(t), mpc(t).bus);
>> results(t) = runopf(mpc(t));
>> end
>>
>> ++++++++++++++++++++++++++++++++++++++++++++
>>
>> put how i can do the same job but using the function
>>
>> most_summary()
>>
>> can you write the code, because i tried a lot but it did not work using the
>> function most_summary()
>>
>> On 19 April 2018 at 23:56, Ray Zimmerman <[email protected]
>> <mailto:[email protected]>> wrote:
>> Use the scale_load()
>> <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/scale_load.html>
>> function in a loop, something like shown here:
>> https://www.mail-archive.com/[email protected]/msg06477.html
>> <https://www.mail-archive.com/[email protected]/msg06477.html>
>>
>> Ray
>>
>>
>>> On Apr 17, 2018, at 12:27 PM, Mohammed Alhajri <[email protected]
>>> <mailto:[email protected]>> wrote:
>>>
>>> Hello all,
>>>
>>> i want to run dc opf for 24 hr
>>>
>>> i have the load scale for 24hr, and i want to multiply each load bus with
>>> this load scale to get the load for 24hr in MW
>>>
>>> how can i do that?
>>>
>>> i have attached:
>>>
>>> - case file
>>> -load scale
>>> -dcopf code
>>>
>>>
>>>
>>> <load scale.xlsx><OETC_DCOPF_2.m><OETC_OPF.m>
>>
>>
>
>
