설명

몇 가지 유용한 기능을 갖춘 매우 간단한 MPEC 모델:
 a) 독특한 해결책이 있다
 b) 일치하는 모든 변수는 상한입니다.

기고자: Steve Dirkse, 2013년 9월

소형 모델 유형 :MPEC

카테고리 : 슬롯 게임 테스트 라이브러리

메인 파일 : mpec06.gms

$title 단순 MPEC 고유 솔루션 상한 일치(MPEC06,SEQ=617)

$onText
  몇 가지 유용한 기능을 갖춘 매우 간단한 MPEC 모델:
   a) 독특한 해결책이 있다
   b) 일치하는 모든 변수는 상한입니다.

기고자: Steve Dirkse, 2013년 9월
$offText

$설정되지 않은 경우 TESTTOL $set TESTTOL 1e-4

스칼라
  톨 / %TESTTOL% /
  c1 / 5 /
  c2 / 8 /
  ;
자유 변수 x1, x2, y;
x1.up = -1;
x2.up = -2;
y.up = 0;
자유 변수 z;
방정식 f1, f2, o;

f1 .. 2*x1 + x2 + y + c1 =L= 0;
f2 .. x1 + 2*x2 + y + c2 =N= 0;
o .. sqr(x1+5) + sqr(x2-2) + sqr(y) =E= z;

모델 m / f1.x1, f2.x2, o /;

mpec min z를 사용하여 m을 해결합니다.
if (m.solvestat = %solveStat.capabilityProblems%),
  abort$[m.modelstat <> %modelStat.noSolutionReturned%] '잘못된 상태 코드',
    m.solvestat, m.modelstat;
  abort.noerror '해결 방법이 없습니다. 더 이상 확인할 필요가 없습니다.';
그렇지 않으면
  abort$[not((m.solvestat = %solveStat.normalCompletion%) 및 (m.modelstat = %modelStat.optimal% 또는
              m.modelstat = %modelStat.locallyOptimal% 또는 m.modelstat = %modelStat.feasibleSolution%))] '잘못된 상태 코드',
    m.solvestat, m.modelstat;
;

스칼라 r1, r2;
r1 = 2*x1.l + x2.l + y.l + c1;
r2 = x1.l + 2*x2.l + y.l + c2;

파일 로그 / '' /;
putclose 로그
 ' ' /
 'F1 :=', r1:6:2, ' perp x1 ', x1.lo:8:2, ' <= ', x1.L:6:2 ' <= ', x1.up:6:2 /
 'F2 :=', r2:6:2, ' perp x2 ', x2.lo:8:2, ' <= ', x2.L:6:2 ' <= ', x2.up:6:2 /
 ' y = ', y.L:6:2 /
 ' ' / ;

abort$(abs(x1.l+1.0) > tol) 'var x1.l은 -1.0이어야 합니다.', x1.l;
abort$(abs(x2.l+2.4) > tol) 'var x2.l은 -2.4여야 합니다.', x2.l;
abort$(abs( y.l+2.2) > tol) 'var y.l은 -2.2여야 합니다.', y.l;
abort$(abs( z.l-40.2)> tol) 'var z.l은 40.2여야 합니다.', z.l;

abort$(abs(f1.l+6.6) > tol) 'equ f1.l은 -6.6이어야 합니다.', f1.l;
abort$(abs(f2.l+0.0) > tol) 'equ f2.l은 0.0이어야 합니다.', f2.l;
abort$(abs( o.l-0.0) > tol) 'equ o.l은 0.0이어야 합니다', o.l;