FMIN_COBYLA에 대한 제약 함수 (scipy optimize)

Question

저는 fmin_cobyla 함수를 사용하여 Scipy Optimize를 사용 중입니다.

1. 모든 항목을 확인 모든 항목> = 0

사람이 있습니까 통찰력을 1

로 요약 확인 : 나는 제약 함수를 작성하기 위해 노력했습니다? 예를 들어 fmin_slsqp의 경우 함수가 있습니다.

def w_constraint(w, v, x0, x1): 
    return np.sum(w) - 1 

함수의 경계는 물론입니다.

하지만 이들은 구획에서 작동하지 않습니다.

Answer 1

불평등 쌍으로 평등을 공식화하십시오. 의미 :

• x >= y
• x <= y = -x >= -y

물론이 수치 문제를 (해결 자 많은 종류의 상처) 소개 할 수 있습니다,하지만 난 보았다

• x == y하면된다 이것은 COBYLA 인터페이스의 여러 구현에서 (equality-constraint support implem 그런 유혹).

  NLopt의 문서에서 주석 :

  (기본이 COBYLA 코드는 불평등 제약을 지원하는 평등 제약은 자동으로이 알고리즘의 경우 원인이없는 것 같다 불평등 제약의 쌍으로 변환됩니다. 문제가됩니다.)