설명

이것은 해결책이 존재한다는 것을 보여주는 간단한 평형 문제입니다.
개별 최적화 문제에 대한 존재를 보장하지는 않습니다.
평형 솔루션.

개별 상담원 문제는 해결 가능한 LP입니다.
전반적인 문제에는 해결책이 없습니다.

변수 x2를 양수로 강제하면 평형 솔루션이 가능해집니다.

기고자: Michael Ferris, 2010년 3월

소형 모델 유형 :EQUIL

카테고리 : 크레이지 슬롯 EMP 라이브러리

메인 파일 : simpequil.gms

$title 단순 평형 (SIMPEQUIL,SEQ=49)

$onText
이는 해의 존재를 보여주는 간단한 평형 문제입니다.
개별 최적화 문제에 대한 존재를 보장하지는 않습니다.
평형 솔루션.

개별 상담원 문제는 해결 가능한 LP입니다.
전반적인 문제에는 해결책이 없습니다.

변수 x2를 양수로 강제하면 평형 솔루션이 가능해집니다.

기고자: Michael Ferris, 2010년 3월

$offText

옵션 limrow=0, limcol = 0;

변수 z1, z2;
양의 변수 x1, x2;
방정식 obj1, obj2, cons1, cons2;

obj1.. z1 =e= x1 + x2;
단점1.. x1*x2 =g= 1;

obj2.. z2 =e= x2;
단점2.. x1 + x2 =g= 1;

모델 구성품 /
  obj1, obj2,
  단점1, 단점2
      /;

x1.l = 1; x2.l = 1;

파일 정보 / '%emp.info%' /;
정보 입력 / '평형';
넣다 / '최소 z1 x1 obj1 cons1';
putclose / '최소 z2 x2 obj2 cons2';

emp를 사용하여 comp를 해결합니다.
* 실제로 실행 불가능한지 확인
중단$(comp.solvestat <> %solveStat.normalCompletion%
       또는 ((comp.modelstat <> %modelStat.infeasible%)
           and (comp.modelstat <> %modelStat.locallyInfeasible%)) ) '첫 번째 모델은 실현 불가능하지 않음';

* x2에 양의 하한을 추가하면 문제가 해결됩니다.
x2.lo = 0.5;

emp를 사용하여 comp를 해결합니다.
중단$(comp.solvestat <> %solveStat.normalCompletion%
       또는 comp.modelstat <> %modelStat.locallyOptimal%) '두 번째 모델이 국소적으로 최적이 아님';