Использование лямбда для функции ограничения

Question

import numpy 
from numpy import asarray 

Initial = numpy.asarray [2.0, 4.0, 5.0, 3.0, 5.0, 6.0]  # Initial values to start with 


bounds = [(1, 5000), (1, 6000), (2, 100000), (1, 50000), (1.0, 5000), (2, 1000000)] 

# actual passed bounds 

b1 = lambda x: numpy.asarray([1.4*x[0] - x[0]]) 
b2 = lambda x: numpy.asarray([1.4*x[1] - x[1]]) 
b3 = lambda x: numpy.asarray([x[2] - x[3]])  
constraints = numpy.asarray([b1, b2, b3]) 

opt= optimize.fmin_slsqp(func,Initial,ieqcons=constraints,bounds=bounds, full_output=True,iter=200,iprint=2, acc=0.01) 

Проблема: Я хочу передать в неравенство ограничений. Считают, что у меня есть 6 параметров

[ a, b, c, d, e, f] 

в Initial ценности, и мои ограничения:

a<=e<=1.4*a ('e' varies from a to 1.4*a) 
b<=f<=1.4*b ('f' varies from b to 1.4*b) 
c>d   ('c' must always be greater than d) 

Но это не работает должным образом. Я не знаю, что это за ошибка. Есть ли лучший способ передать мои ограничения как функцию? Пожалуйста, помогите мне.

Answer 1

Основываясь на комментарии Роберта Керна, я удалил свой предыдущий ответ. Вот ограничения как непрерывные функции:

b1 = lambda x: x[4]-x[0] if x[4]<1.2*x[0] else 1.4*x[0]-x[4] 
b2 = lambda x: x[5]-x[1] if x[5]<1.2*x[1] else 1.4*x[1]-x[5] 
b3 = lambda x: x[2]-x[3] 

Примечание: для этого синтаксиса требуется Python 2.5 или более.

Чтобы получить ограничение a<=e<=1.4*a, обратите внимание, что 1.2*a это точка на полпути между a и 1.4*a.

Под этим пунктом, то есть все e<1.2*a, мы используем непрерывную функцию e-a. Таким образом, общая функция ограничения отрицательна, когда e<a, обрабатывая условие нижнего запредельного состояния, ноль на нижней границе e==a, а затем положительно для e>a до половины.

Выше половины, то есть всего e>1.2*a, мы используем вместо этого непрерывную функцию 1.4*a-e. Это означает, что общая функция ограничения отрицательна, когда e>1.4*a, обрабатывая условие верхнего вне поля, ноль на верхней границе e==1.4*a, а затем положительно, когда e<1.4*a, вплоть до половины точки.

В середине пути, где e==1.2*a, обе функции имеют одинаковое значение. Это означает, что общая функция непрерывна.

Код: documentation for ieqcons.

^{- Здесь предварительно Python 2.5 Синтаксис: b1 = lambda x: (1.4*x[0]-x[4], x[4]-x[0])[x[4]<1.2*x[0]]}