scipy.optimize.shgo ValueError: Значение истинности массива с более чем одним элементом является неоднозначным. Используйте a.any () или a.all () - PullRequest
Новая
       7

scipy.optimize.shgo ValueError: Значение истинности массива с более чем одним элементом является неоднозначным. Используйте a.any () или a.all ()

0 голосов
/ 08 марта 2019

Я пытаюсь подогнать функцию y(x,T,p), чтобы получить коэффициенты a, b, c, d, e, f. Данные для y, x, T, p известны. С глобальным оптимизатором я хочу найти хорошую отправную точку. shgo кажется единственным, кто принимает constraints. 

import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import shgo

# test data
x = np.array([0.1,0.2,0.3,1])
T = np.array([300,300,300,300])
p = np.array([67.2,67.2,67.2,67.2])
y = np.array([30,50,55,67.2])

# function
def func(pars,x,T,p):
    a,b,c,d,e,f = pars
    return x*p+x*(1-x)*(a+b*T+c*T**2+d*x+e*x*T+f*x*T**2)*p

# residual
def resid(pars):
    return ((func(pars,x,T,p) - y) ** 2).sum()

# constraint: derivation is positive in every data point
def der(pars):
    a,b,c,d,e,f = pars
    return -p*((3*f*T**2+3*e*T+3*d)*x**2+((2*c-2*f)*T**2+(2*b-2*e)*T-2*d+2*a)*x-c*T**2-b*T-a-1)

con1 = ({'type':'ineq', 'fun':der})

# minimizer shgo
bounds = [(-1,1),(-1,1),(-1,1),(-1,1),(-1,1),(-1,1)]
res = shgo(resid, bounds, constraints=con1)
print("a = %f , b = %f, c = %f, d = %f, e = %f, f = %f" % (res[0], res[1], res[2], res[3], res[4], res[5]))

# plotting
x0 = np.linspace(0, 1, 100)
fig, ax = plt.subplots()
fig.dpi = 80
ax.plot(x,y,'ro',label='data')
for i,txt in enumerate(T):
    ax.annotate(txt,(x[i],y[i]))
ax.plot(x0, func(res.x, x0, 300,67.2), '-', label='fit1')
plt.xlabel('x')
plt.ylabel('y')
plt.legend()
plt.show()

С этим я получаю ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all() Я не знаю, что означает эта ошибка, и другие темы с той же ошибкой не очень помогают мне понять. Когда я использую локальный минимизатор (scipy.optimize.minimize с методом cobyla), ошибка не появляется.

Может ли кто-нибудь помочь мне понять мою проблему или даже помочь решить ее? Спасибо

EDIT: 

Traceback (most recent call last):
  File "C:\Users\...\Python\Python36\site-packages\scipy\optimize\_shgo_lib\triangulation.py", line 759, in __getitem__
    return self.cache[x]
KeyError: (0, 0, 0, 0, 0, 0)

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "C:/Users/.../test.py", line 70, in <module>
    res = shgo(resid, bounds, constraints=con1)
  File "C:\Users\...\Python\Python36\site-packages\scipy\optimize\_shgo.py", line 423, in shgo
    shc.construct_complex()
  File "C:\Users\...\Python\Python36\site-packages\scipy\optimize\_shgo.py", line 726, in construct_complex
    self.iterate()
  File "C:\Users\...\Python\Python36\site-packages\scipy\optimize\_shgo.py", line 869, in iterate
    self.iterate_complex()
  File "C:\Users\...\Python\Python36\site-packages\scipy\optimize\_shgo.py", line 890, in iterate_hypercube
    self.g_args)
  File "C:\Users\...\Python\Python36\site-packages\scipy\optimize\_shgo_lib\triangulation.py", line 121, in __init__
    self.n_cube(dim, symmetry=symmetry)
  File "C:\Users\...\Python\Python36\site-packages\scipy\optimize\_shgo_lib\triangulation.py", line 172, in n_cube
    self.C0.add_vertex(self.V[origintuple])
  File "C:\Users\...\Python\Python36\site-packages\scipy\optimize\_shgo_lib\triangulation.py", line 767, in __getitem__
    index=self.index)
  File "C:\Users\...\Python\Python36\site-packages\scipy\optimize\_shgo_lib\triangulation.py", line 681, in __init__
    if g(self.x_a, *args) < 0.0:
ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()
...