时间戳

首页 > 解决方案 > 在 Julia-JuMP 中实现非线性约束的问题

julia - 在 Julia-JuMP 中实现非线性约束的问题

问题描述

我正在尝试解决一个优化问题，该问题涉及在其约束之一中矩阵的逆。简化版如下：

a = [1 1e-3;1e-3 1]
n = 2
f(Q...) = inv(reshape(collect(Q), (n, n)))
g(x...) = reshape(collect(x), (n,n))
NSD = Model(with_optimizer(Ipopt.Optimizer))
@variable(NSD, Q[1:n, 1:n])
@variable(NSD, aux[1:n, 1:n])
register(NSD, :f, n*n, f, autodiff=true)
register(NSD, :g, n*n, g, autodiff=true)
@constraint(NSD, aux .== a)
@NLconstraint(NSD, invconst[i=1:n, j=1:n], f(Q...)[i,j] >= g(aux...)[i,j])
@objective(NSD, Max, 0)
optimize!(NSD)

但是，此代码不断向我显示以下我无法解决的错误消息： LoadError: Cannot divide a quadratic expression by an aff. expression

有人可以帮我解决这里的问题吗？

标签： juliajulia-jump

解决方案

用户定义的函数必须返回单个标量。他们不能返回一个矩阵。

您n的应用程序中有什么？如果n > 2，这几乎肯定不是您想要做的。

如果n = 2，只需明确地写出逆。

我正在添加一个对 JuMP 有帮助的示例：https ://github.com/jump-dev/JuMP.jl/pull/2379

推荐阅读