2016-07-10 14 views
1

Мне нужно сгенерировать все точки парето, но я получил эту ошибку.JuMP: LoadError: непризнанная функция: f1 используется в нелинейном выражении

using JuMP 
using Gurobi 
using Gadfly 
using Ipopt 

m = Model(solver=IpoptSolver(print_level=0)) 
@variable(m, 0.1 <= x <= 1.0) 
@variable(m, 0.0 <= y <= 1.0) 

pareto_x = Float16[] 
pareto_y = Float16[] 

for i in 0.0:0.1:1.0 
    for j in 0.0:0.1:1.0 

     f1(x,y) = x 
     f2(x,y) = (2.0-exp(-((y-0.2)/0.004)^2)-0.8*exp(-((y-0.6)/0.4)^2))/x 

     @NLobjective(m, Min, i*f1(x,y) + j*f2(x,y)) ## <<-- ERROR HERE 

     status = solve(m) 

     println("Objective value: ", getobjectivevalue(m)) 
     x_opt = getvalue(x) 
     y_opt = getvalue(y) 
     println("x = ", x_opt) 
     println("y = ", y_opt) 

     push!(pareto_x,f1(x_opt,y_opt)) 
     push!(pareto_y,f2(x_opt,y_opt)) 
    end 
end 
plot(x=pareto_x, y=pareto_y) 
+0

Что происходит, если заменить все экземпляры 'f1' с' x'? Моя кишка говорит мне, что вам запрещены вызовы произвольных функций в '@ NLobjective' –

+0

Да, но я бы использовал произвольную функцию для повторного использования этой функции для создания точек парето. – sparkle

+0

Для этого варианта использования, поскольку ваши функции являются просто выражениями замкнутой формы, вы должны просто использовать '@ NLexpression' вместо пользовательских функций. Это, как правило, будет работать лучше, потому что JuMP в настоящее время не вычисляет производные второго порядка, если есть определенные пользователем функции. – mlubin

ответ

1

http://jump.readthedocs.io/en/latest/nlp.html#user-defined-functions

mysquare(x) = x^2 
myf(x,y) = (x-1)^2+(y-2)^2 

JuMP.register(:myf, 2, myf, autodiff=true) 
JuMP.register(:mysquare, 1, mysquare, autodiff=true) 

m = Model() 

@variable(m, x[1:2] >= 0.5) 
@NLobjective(m, Min, myf(x[1],mysquare(x[2]))) 

решаемые с этим

f1(x,y) = x 
f2(x,y) = (2.0-exp(-((y-0.2)/0.004)^2)-0.8*exp(-((y-0.6)/0.4)^2))/x 

isdefined(:f1) || JuMP.register(:f1, 2, f1, autodiff=true) 
isdefined(:f2) || JuMP.register(:f2, 2, f2, autodiff=true) 

@objective(m, Min, f1(x,y)) 
+0

Обратите внимание, что 'JuMP.register' определяет только функции для использования в нелинейных выражениях, например,' @ NLobjective' или '@ NLconstraint'. Вышеприведенный пример работает случайно, поскольку 'f1' является линейной функцией, используемой в' @ object'. – mlubin