5

Param我正在尝试将三次样条拟合到下面的数据点，当我使用类似示例m.x = m.Param(value=np.linspace(-1, 6))或使用常量时我有点困惑Const。

import numpy as np
from gekko import GEKKO
xm = np.array([0, 1, 2, 3, 4, 5])
ym = np.array([0.1, 0.2, 0.3, 0.5, 1.0, 0.9])

m = GEKKO()
m.x = m.Param(value=np.linspace(-1, 6))
m.y = m.Var()
m.options.IMODE = 2
m.cspline(m.x, m.y, xm, ym)
m.solve(disp=False)

p = GEKKO()
p.x = p.Var(value=1, lb=0, ub=5)
p.y = p.Var()
p.cspline(p.x, p.y, xm, ym)
p.Obj(-p.y)
p.solve(disp=False)
4

1 回答 1

1

常数Const是预期不会改变的单个标量值。参数Param具有初始值，但用户可以通过数据更改它。固定值FV或操纵变量MV是两种特殊类型的参数，它们具有成为求解器决策变量的额外选项。FVan和之间的区别在于MVanFV具有一个值，而 an可以在数据 ( ) 或时间 ( ) 维度MV上具有不同的值。IMODE=2IMODE=4-9

您有一个很好的例子，将三次样条拟合到数据，然后求解范围内的最大值0 < x < 5。

三次样条最大化

import numpy as np
from gekko import GEKKO
xm = np.array([0, 1, 2, 3, 4, 5])
ym = np.array([0.1, 0.2, 0.3, 0.5, 1.0, 0.9])

m = GEKKO()
m.x = m.Param(value=np.linspace(-1, 6))
m.y = m.Var()
m.options.IMODE = 2
m.cspline(m.x, m.y, xm, ym)
m.solve(disp=False)

p = GEKKO()
p.x = p.Var(value=1, lb=0, ub=5)
p.y = p.Var()
p.cspline(p.x, p.y, xm, ym)
p.Maximize(p.y)
p.solve(disp=False)

import matplotlib.pyplot as plt
plt.plot(xm,ym,'rs',label='Data')
plt.plot(m.x,m.y,'r.-',label='Cubic Spline')
plt.plot(p.x,p.y,'bo',label='Maximize')
plt.xlabel('x'), plt.ylabel('y')
plt.legend()
plt.show()

在您的情况下， a Param（或MVwith STATUS=0）是适当的 gekko 对象。AConst给出一个错误：ValueError: Constant value must be scalar.

于 2019-11-22T11:55:04.770 回答