If you look at equation (1) it states that the cost function (not the marginal 
cost) of the generators is given by:

C(P) = a + b*P + c*P^2

Thus the marginal cost function (which enters the objective function as 
parameters in the f and H matrix) is:

Mc = b + 2*c*P

Therefore I assume that you have to double all the values in column 5 of the 
gencost matrix to get correct results:

gencost = [
    2   0   0   3       0.008   10  0;
    2   0   0   3       0.012   15  0;
    2   0   0   3       0.016   50  0;
    2   0   0   3       0.01   12  0;
    2   0   0   3       0.012   15.5    0;
    2   0   0   3       0.014   15.5    0;
    2   0   0   3       0.016   21.5    0;
    2   0   0   3       0.012   16  0;
    2   0   0   3       0.01   14  0;
    2   0   0   3       0.008   13  0;
    2   0   0   3       0.012   16  0;
    2   0   0   3       0.018   31  0;
];

I haven't tried it, but I will assume it works. Look at the lambdas after 
simulating. It should also be safe to set all loads in the Pd column of the bus 
matrix.

From: [email protected] 
[mailto:[email protected]] On Behalf Of xiwen wang
Sent: 3. mars 2010 01:05
To: MATPOWER discussion forum; [email protected]
Subject: simple AC OPF problem

Dear Ray,
Thanks very much for the wonderful tool.
I am repeating the 11-bus case as in the attached paper, but I can not get the 
results for the base case (table 2 and 3).
Also attached the mat case file. Please let me know what I have done wrong.

If anyone happened done this before, please share the case file with me.

I really appreciate your help.

Cheers,
X


Reply via email to