3

私は PyOmo を使用して、最終的に Ipopt で解決される非線形モデルを生成しています。モデルは次のとおりです。

from pyomo.environ import *
from pyomo.dae import *

m = ConcreteModel()

m.t = ContinuousSet(bounds=(0,100)) 

m.T  = Param(default=100,mutable=True)
m.a  = Param(default=0.1)
m.kP = Param(default=20)
m.P  = Var(m.t, bounds=(0,None))
m.S  = Var(m.t, bounds=(0,None))
m.u  = Var(m.t, bounds=(0,1), initialize=0.5)

m.Pdot = DerivativeVar(m.P)
m.Sdot = DerivativeVar(m.S)

m.obj = Objective(expr=m.S[100],sense=maximize)

def _Pdot(M,i):
  if i == 0:
    return Constraint.Skip
  return M.Pdot[i] == (1-M.u[i])*(M.P[i]**0.75)

def _Sdot(M,i):
  if i == 0:
    return Constraint.Skip
  return M.Sdot[i] == M.u[i]*0.2*(M.P[i]**0.75)

def _init(M):
  yield M.P[0] == 2
  yield M.S[0] == 0
  yield ConstraintList.End

m.Pdotcon         = Constraint(m.t, rule=_Pdot)
m.Sdotcon         = Constraint(m.t, rule=_Sdot)
m.init_conditions = ConstraintList(rule=_init)

discretizer = TransformationFactory('dae.collocation')
discretizer.apply_to(m,nfe=100,ncp=3,scheme='LAGRANGE-RADAU')
discretizer.reduce_collocation_points(m,var=m.u,ncp=1,contset=m.t)

solver  = SolverFactory('ipopt')
results = solver.solve(m,tee=False)

モデルを実行すると、次のエラーが発生します。

Error evaluating constraint 1: can't evaluate pow'(0,0.75).

Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "/usr/local/lib/python3.5/dist-packages/pyomo/opt/base/solvers.py", line 577, in solve
    "Solver (%s) did not exit normally" % self.name)
pyutilib.common._exceptions.ApplicationError: Solver (asl) did not exit normally

エラーの最初の部分は Ipopt によるもので、2 番目の部分は PyOmo によるものです。M.P[i]**0.75明らかに、この問題は制約内の用語と関係がありませんが、電力を変更しても問題は解決しません (ただし、2.0うまくいきました)。

どうすればこれを解決できますか?

4

2 に答える 2

3

エラー メッセージには、pow'(0,0.75)評価できないことが示されています。この関数の'文字は、1 次導関数を示します (''は 2 次導関数を示します)。メッセージは事実上、一次導関数が存在しないか、ゼロで無限大になることを示しています。

この問題は簡単に解決できます。次のように、変数をゼロ以外の値にバインドします。

m.P  = Var(m.t, bounds=(1e-20,None))
m.S  = Var(m.t, bounds=(1e-20,None))
于 2016-06-17T00:38:13.990 に答える
2

リチャードの答えに追加します：

ipopt は、指定されていない場合は 0 を想定するため、変数の初期値を更新する必要がある場合もあります。そのため、最初の反復では変数が 0 で評価されます。

したがって：

m.P  = Var(m.t, bounds=(1e-20,None), initialize=1e-20)
m.S  = Var(m.t, bounds=(1e-20,None), initialize=1e-20)

初期化として 1e-20 の代わりに、問題により関連する値を使用できます

于 2020-06-18T11:41:36.310 に答える