我正在尝试解决非线性优化问题。我通过创建下面的代码重复了我的问题。 Python 返回 TypeError: object of type 'int' has no len()
。如何在约束函数中包含 IF 语句?
控制台打印以下内容:
File "<ipython-input-196-8d29d410dcea>", line 1, in <module>
runfile('C:/Users/***/Documents/***/Project/untitled.py', wdir='C:/Users/***/Documents/***/***/Project')
File "C:\Users\***\Anaconda3\lib\site-packages\spyder_kernels\customize\spydercustomize.py", line 704, in runfile
execfile(filename, namespace)
File "C:\Users\***\Anaconda3\lib\site-packages\spyder_kernels\customize\spydercustomize.py", line 108, in execfile
exec(compile(f.read(), filename, 'exec'), namespace)
File "C:/Users/***/Documents/***/***/Project/untitled.py", line 27, in <module>
m.Equation(Cx(x1,x2,x3,x4) < 0)
File "C:/Users/***/Documents/***/***/Project/untitled.py", line 17, in Cx
if K > 15:
File "C:\Users\***\Anaconda3\lib\site-packages\gekko\gk_operators.py", line 25, in __len__
return len(self.value)
File "C:\Users\***\Anaconda3\lib\site-packages\gekko\gk_operators.py", line 134, in __len__
return len(self.value)
TypeError: object of type 'int' has no len()
-
from gekko import GEKKO
m = GEKKO()
def Cr(x1,x2,x3,x4):
return (x1*x4*(x1+x2+x3)+x3**2)
def Cw(x1,x2,x3,x4):
return x1*x2*x3*x4
def Ck(x1,x2,x3,x4):
return x1*x2*x3*x4+1
def Cx(x1,x2,x3,x4):
K = Ck(x1,x2,x3,x4)
if K > 15: #Issue here
K = 15
return x1**2+x2**2+x3**2+x4**2 - K
x1 = m.Var(value=1,lb=-5000,ub=5000)
x2 = m.Var(value=1,lb=-5000,ub=5000)
x3 = m.Var(value=1,lb=-5000,ub=5000)
x4 = m.Var(value=1,lb=-5000,ub=5000)
m.Equation(Cw(x1,x2,x3,x4) >= 14)
m.Equation(Cx(x1,x2,x3,x4) < 0)
m.Obj(Cr(x1,x2,x3,x4))
m.solve(disp=False)
print(x1.value)
print(x2.value)
print(x3.value)
print(x4.value)
-
我希望让 GEKKO 在约束中使用 IF 语句运行,我不关心代码中的优化问题是否有解决方案。提前谢谢你。
最佳答案
(免责声明:我不知道这个库或它能为你做什么)
if 语句使这个问题不可微,这使 NLP 求解器(如 Ipopt)的假设无效。
就 MINLP 求解器(Bonmin、Couenne)而言,这可以通过重新表述来实现(并且当所需的辅助二进制变量已放宽时,所产生的问题是可微分的)。期待 lib 为您做这件事是毫无疑问的。
因此,您似乎需要遵守某些 MINLP 模型的规则,例如 Bonmin here 所描述的.没有“基于 if 的分支”的概念。
要么像 MIP 世界中常见的那样引入指示变量,请参阅 here .忽略开销这个想法是这样的:
K_ = Ck(x1,x2,x3,x4)
I = K_ > 15 (binary variable; see link for formulation idea)
return x1**2+x2**2+x3**2+x4**2 - I*15 - (1-I) * K_
那么这是一个MINLP。
在解释方程式时,您可能无需使用额外的二进制变量(并接触 MINLP)就可以逃脱,例如:
return x1**2+x2**2+x3**2+x4**2 - min(Ck(x1,x2,x3,x4), 15)
这也是不可微分的,但可以很容易地重新表述(有一个怪癖),例如:
return x1**2+x2**2+x3**2+x4**2 - A
# extra constraints
A <= Ck(x1,x2,x3,x4)
A <= 15
如果我们能够强制朝着尽可能大的 A 迈进。这意味着,它必须是目标的一部分:
m.Obj(Cr(x1,x2,x3,x4) + c * A) (if it's a maximization problem)
这将是一个 NLP,但是 c 的值需要注意(它必须足够大)。
关于python - Gekko 非线性优化,约束函数评估 if 语句时出现对象类型错误,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/55436296/