Мне нужно сгенерировать все точки парето, но я получил эту ошибку.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)
Что происходит, если заменить все экземпляры 'f1' с' x'? Моя кишка говорит мне, что вам запрещены вызовы произвольных функций в '@ NLobjective' –
Да, но я бы использовал произвольную функцию для повторного использования этой функции для создания точек парето. – sparkle
Для этого варианта использования, поскольку ваши функции являются просто выражениями замкнутой формы, вы должны просто использовать '@ NLexpression' вместо пользовательских функций. Это, как правило, будет работать лучше, потому что JuMP в настоящее время не вычисляет производные второго порядка, если есть определенные пользователем функции. – mlubin