Hi Fred, I greatly appreciate! I understand what the problem was. Would you please help with a similar system as follows. This is the actual system which I am dealing with and the one I've already given you was a reduced version of it. It's a 5 bus system with 4 branches. I have attached a figure of the system.
mpc.baseMVA = 1;
mpc.bus = [
1 3 0 0 0 0 1 1.06 0 0 1 1 1;
2 2 0 0 0 0 1 1.045 -4.98 0 1 1 1;
3 1 0 0 0 0 1 1.01 -12.72 0 1 1.06
0;
4 1 0 0 0 0 1 1.019 -10.33 0 1 1.06
0;
5 1 0.5 0 0 0 1 1.06 0 0 1 1.06 0;
];
mpc.gen = [
1 2 0 100 -100 1.02 100 1 318 0 0 0
0 0 0 0 0 0 0 0 0;
2 2 0 100 -100 1.02 100 1 318 0 0 0
0 0 0 0 0 0 0 0 0;
];
mpc.branch = [
1 3 .2 0 0 9900 0 0 0 0 1 -360 360;
2 3 .2 0 0 9900 0 0 0 0 1 -360 360;
3 4 .2 0 0 9900 0 0 0 0 1 -360 360;
4 5 .2 0 0 9900 0 0 0 0 1 -360 360;
];
mpc.branch(:,6)=inf;%no limits on line's power
mpc.gencost = [
2 0 0 3 0.0430293 20 0;
2 0 0 3 0.0430293 20 0;
];
And, this is the manual solution:
x0=[
0
0
0
0
0
1
1
0.9
0.7
0.5
0.5
0.5
0
0
];
Best Regards,
Arash
On Mon, Dec 16, 2013 at 1:14 AM, Fred Kanjelesa <[email protected]>wrote:
> Hi Arash,
> It appears to me you may have omiited the generator data in your 4-bus
> case.
>
> i tried your case with some simple definition of gen. data and it run for
> both power flow and optimal power flow. generator data is as follows below
>
>
> %% generator data
> % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max
> Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf
> mpc.gen = [
> 1 2 0 100 -100 1.02 100 1 318 0 0 0 0 0 0 0 0 0 0 0 0;
> ];
>
> Please find attached the case modeling and the results,
>
> best regards
>
>
> On 16 December 2013 10:48, Arash Alimardani <[email protected]>wrote:
>
>> Hello all,
>>
>> I have the simplest OPF problem, but I can't get matpower to do it for
>> me. I would appreciate if someone could help me.
>>
>> I have attached a figure of the system. Here is the description of the
>> system:
>>
>> 4 buses including 1 generator and 3 load buses:
>>
>> mpc.bus = [
>> 1 3 0 0 0 0 1 1.06 0 0 1 1.06
>> 0.94;
>> 2 1 0 0 0 0 1 1.045 -4.98 0 1
>> 1.06 0.94;
>> 3 1 0 0 0 0 1 1.01 -12.72 0 1
>> 1.06 0.94;
>> 4 1 .25 0 0 0 1 1.019 -10.33 0 1
>> 1.06 0.94;
>> ];
>>
>> as you can see the load of bus 4 is P=0.25, Q=0. All other loads are zero.
>>
>> As for branches:
>>
>> mpc.branch = [
>> 1 2 1/3 0 0 9900 0 0 0 0 1 -360
>> 360;
>> 2 3 1/3 0 0 9900 0 0 0 0 1 -360
>> 360;
>> 3 4 1/3 0 0 9900 0 0 0 0 1 -360
>> 360;
>> ];
>>
>> As you can see, only a resistance with the value of 0.333 is the
>> impedance of each line.
>>
>> The data is in per unit and:
>> mpc.baseMVA = 1;
>>
>> And the cost of production:
>> mpc.gencost = [
>> 2 0 0 3 0.0430293 20 0;
>> ];
>>
>>
>> The solution which I have by solving it manually is:
>> x0=[0;0;0;0;1;.833333;.666666;.5;.5;0];
>>
>> However, matpower doesn't converge to it even when I give it as the
>> initial value (x0).
>>
>> Please, help me figure out the problem with my work.
>>
>> Regards,
>> Arash Alimardani
>>
>
>
>
> --
>
> *Fred Kanjelesa.*NCC-Trading/ZESCO
> Zambia/Lusaka.
>
> Mobile: +260967/955-880912
> Extension: 5124
>
> "I can do all things through Christ who strengthens me." *Philippians
> 4:13*
>
Visio-Drawing1.pdf
Description: Adobe PDF document
