참조

Neculai Andrei, 슬롯 사이트 기술을 사용한 비선형 최적화 애플리케이션, 스프링거 최적화 및 그 애플리케이션, 모델컨트롤3(12.25) 장최적 제어, 2013

메인파일 : control3.gms

$onText
비선형 동적 구속조건과 경계를 이용한 최적 제어 문제
일반 비선형 계획법 문제로 해결된 조건.

Divya Garg 외, 직접 궤적 최적화 및 비용 추정
유한수평과 무한수평 최적제어 문제
Radau 의사 스펙트럼 방법. 전산 최적화 및 응용,
vol.49, nr. 2011년 6월 2일, pp. 335-358.
$offText

n 상태 설정 / state1 /
     k 설정 /t1*t100/
     ku(k) 제어 범위
     ki(k) 초기 기간
     kt(k) 단말기 기간 ;

ku(k) = yes$(ord(k) lt 카드(k));
ki(k) = yes$(ord(k) eq 1);
kt(k) = ku(k)가 아님;

k, ki, kt, ku를 표시합니다.

매개변수 rk 페널티 제어 / 0.01 /
          xinit(n) 초기값 / state1 2 / ;

변수 x(n,k) 상태변수
          u(k) 제어변수
          j 기준

방정식 비용 기준 정의
          stateq(n,k) 상태 방정식 ;

비용..
j =e= .5*합((k,n), (x(n,k)) +
      .5*sum((ku), (u(ku))*rk*(u(ku))) );

상태q(n,k+1)..

x(n,k+1) =e= 2*x(n,k) + 2*u(k)*sqrt(x(n,k)) ;

모델 제어3 /all/;

$ifThenI x%mode%==xbook
x.l(n,k) = xinit(n);
x.fx(n,ki) = xinit(n);
x.fx(n,kt) = 2;
$endIf

nlp를 사용하여 j를 최소화하는 control3을 해결합니다.

x.l, ul을 표시합니다.

$ifThenI x%mode%==xbook
파일 res1 /control3.dat/;
res1을 넣어
루프(k, put x.l('state1',k):10:5,',', put/)

루프(k, u.l(k):10:5,',', put/)
$endIf

* 컨트롤 종료3