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
Visio-Drawing1.pdf
Description: Adobe PDF document
