참조
카테고리 : 슬롯 커뮤니티 NOA 라이브러리
메인파일 : surface.gms
$onText
최소한의 표면 문제.
일부 조건에 따라 그래프의 배열을 최소화하는 함수 f를 찾습니다.
f 도메인의 경계에 대한 제약 조건.
Boyd, S., Vandenberghe, L., 볼록 최적화, Cambridge University Press,
케임브리지, 2004.
$offText
세트 X /I1*I21/;
설정 Y /J1*J21/;
SET 내부(X,Y);
* 내부에서 i1, i21 제외
inside(X,Y)$(not((ord(X)=1) and (ord(X)=card(X)))) = yes;
내부에 표시;
스칼라 K /10/;
변수 obj, f(x,y);
양수 변수 f(X,Y);
* 변수의 한계, 초기 조건, 고정 조건:
f.up(x,y)=1;
f.l(x,y) =1.0;
f.fx(X,Y)$((ord(X)=1) 또는 (ord(X)=card(X))) = 1;
방정식 objfun;
objfun..obj =E= (1/sqr(K)) *
합계((X,Y) $(내부(X,Y)),
sqrt( sqr((F(X+1,Y)-F(X,Y))/K) +
sqr((F(X,Y+1)-F(X,Y))/K) + 1) ) ;
모델 표면 /all/;
obj를 최소화하는 nlp를 사용하여 표면을 해결합니다.
$ifThenI x%mode%==xbook
파일 res1 /surf1.dat/
res1을 넣어;
"배열 표면 ="obj.l; 놓다 /;
loop(Y, put Y.tl:6; loop(X, put F.l(X,Y):6:2 ); put /;) put /;
$endIf
* 끝면