The question about extend OPF and the optimization tool

[neuqxhlinwei]

Hi, 

    I see. Can we add the power loss into objective function in runopf() ?

 

---------------------------------------------------------

Best regards,

Wei Lin 林巍

Zhejiang University, 

Zheda Road 38#, Hangzhou, P.R.China, 310027

+86-13675865210

 

发件人: Ray Zimmerman [mailto:r...@cornell.edu] 
发送时间: 2018年1月11日 23:56
收件人: neuqxhlin...@163.com <mailto:neuqxhlin...@163.com> 
主题: Re: The question about extend OPF and the optimization tool

 

(FYI, in the future, please address such questions to the  
<http://www.pserc.cornell.edu/matpower/mailinglists.html#discusslist> 
MATPOWER-L mailing list, not to me directly)

 

Yes, it is possible to add user-defined costs to the OPF. Please see Sections 
6.3 and 6.4 and Chapter 7 in the MATPOWER User 
<http://www.pserc.cornell.edu/matpower/docs/MATPOWER-manual-6.0.pdf> ’s Manual. 
If you need general non-linear costs, you will need to use the latest 
development version of MATPOWER on GitHub 
<https://github.com/MATPOWER/matpower/>   You can find the corresponding 
version of the manual here 
<https://www.dropbox.com/s/3qkji1cuyi075uj/MATPOWER-manual.pdf?dl=0>   
describing how to use the non-linear user-defined costs and constraints.

 

Regarding YALMIP, MATPOWER does not use it except when using SDP_PF, so you 
cannot use it to modify the problem normally solved by runopf().

 

   Ray

 

 

On Jan 10, 2018, at 9:54 PM, <neuqxhlin...@163.com 
<mailto:neuqxhlin...@163.com> > <neuqxhlin...@163.com 
<mailto:neuqxhlin...@163.com> > wrote:

 

Dear matpower:

Thanks for matpower solving lots of problem in my power simulation.

I have a question about the objective function in OPF. Can the system add some 
user-defined function expecting the generators cost, Like power loss , voltage 
deviation ?

And the matpower can be with some optimization tools like yalmip ? When I use 
yalmip , I define some variable but it can not be solve in matpower?

Code like:

      clc;clear;

x = sdpvar(14,1);

assign(x,ones(14,1));

mpc = loadcase('case14');

mpc.bus(:,4) = mpc.bus(:,4)  + value(x);

mpopt = mpoption('verbose',0,'out.all',0);

result = runpf(mpc,mpopt);

f = sum(abs(get_losses(result)));

F= [sum(x)<=50];

sdp_result = solvesdp(F,f); 

       and the result show the yalmip is successfully sovled , but the result 
is not right.

Hope you can help me.

 

---------------------------------------------------------

Best regards,

Wei Lin 林巍

Zhejiang University, 

Zheda Road 38#, Hangzhou, P.R.China, 310027

+86-13675865210