pinene.gms : 알파-피넨 COPS 2.0의 등각화 #5

설명

열적 등거리화에서 반응 계수를 결정합니다.
알파-피넨.

이 모델은 COPS 벤치마킹 제품군에서 나온 것입니다.
참조http://www-unix.mcs.anl.gov/~more/cops/.이산화 지점 수는 다음 명령을 사용하여 지정할 수 있습니다
라인 매개변수 --nh. NH에 대한 COPS 성능 테스트가 보고되었습니다.
= 25, 50, 100, 200

대형 모델 유형 :NLP

카테고리 : 무료 슬롯 게임 모델 라이브러리

메인 파일 : pinene.gms   포함: copspart.inc

$title 알파-피넨 COPS 2.0 #5의 등거리화 (PINENE,SEQ=233)

$onText
열적 등거리화에서 반응 계수를 결정합니다.
알파-피넨.

이 모델은 COPS 벤치마킹 제품군에서 나온 것입니다.
http://www-unix.mcs.anl.gov/~more/cops/를 참조하세요.

이산화 지점 수는 다음 명령을 사용하여 지정할 수 있습니다.
라인 매개변수 --nh. NH에 대한 COPS 성능 테스트가 보고되었습니다.
= 25, 50, 100, 200

Dolan, E D 등, JJ, 벤치마킹 최적화
COPS가 포함된 소프트웨어. 기술. 대표, 수학과 컴퓨터
과학부, 2000.

Box, G E P, Hunter, W G, MacGregor, J F 및 Erjavec, J,
분석과 관련된 몇 가지 문제
다중응답 데이터. 기술 15(1973), 33-51.

Tjoa, I B 및 Biegler, L T, 동시 솔루션 및
매개변수 추정을 위한 최적화 전략
미분-대수 방정식 시스템. 인디애나
화학. 결의안. 30(1991), 376-385.

Averick, B M, Carter, R G, More, J J 및 Xue, G L,
MINPACK-2 테스트 문제 모음. 기술. 대표,
수학과 컴퓨터 과학부,
아르곤 국립 연구소, 1992.

Ascher, U M, Matteij, R M M 및 Russell, R D,
경계값 문제의 수치해
상미분 방정식의 경우. 1995년 시암.

키워드: 비선형 계획법, 화학 공학, 열 등거리 측정,
          화학 동역학, 반응 속도 계수
$offText

$설정되지 않은 경우 nh $set nh 50

세트
   무료 슬롯 게임 '미분방정식' / ne1*ne5 /
   np 'ODE 매개변수' / np1*np5 /
   nc '배치 포인트' / nc1*nc3 /
   nh '파티션 간격' / nh1*nh%nh% /
   nm '측정값' / 1*8 /;

매개변수
   bc(무료 슬롯 게임) 'ODE 초기 조건' / ne1 100, (ne2*ne5) 0 /
   tau(nm) '관찰이 이루어진 시간' / 1 1230.0, 2 3060.0
                                                 3 4920.0, 4 7800.0
                                                 5 10680.0, 6 15030.0
                                                 7 22620.0, 8 36420.0 /;

테이블 z(nm,무료 슬롯 게임) '관찰'
         ne1 ne2 ne3 ne4 ne5
   1 88.35 7.3 2.3 0.4 1.75
   2 76.4 15.6 4.5 0.7 2.8
   3 65.1 23.1 5.3 1.1 5.8
   4 50.4 32.9 6.0 1.5 9.3
   5 37.5 42.7 6.0 1.9 12.0
   6 25.9 49.1 5.9 2.2 17.0
   7 14.0 57.4 5.1 2.6 21.0
   8 4.5 63.1 3.8 2.9 25.7 ;

$batInclude copspart.inc nc3 8

양의 변수 theta(np) 'ODE 매개변수';

방정식
   collocation_eqn1(nh,nc)
   collocation_eqn2(nh,nc)
   collocation_eqn3(nh,nc)
   collocation_eqn4(nh,nc)
   collocation_eqn5(nh,nc);

collocation_eqn1(i,j).. Duc[i,j,'ne1'] =e= - (theta['np1'] + theta['np2'])*uc[i,j,'ne1'];

collocation_eqn2(i,j).. Duc[i,j,'ne2'] =e= theta['np1']*uc[i,j,'ne1'];

collocation_eqn3(i,j).. Duc[i,j,'ne3'] =e= theta['np2']*uc[i,j,'ne1']
                                        - (세타['np3'] + 세타['np4'])*uc[i,j,'ne3']
                                        + 세타['np5']*uc[i,j,'ne5'];

collocation_eqn4(i,j).. Duc[i,j,'ne4'] =e= theta['np3']*uc[i,j,'ne3'];

collocation_eqn5(i,j).. Duc[i,j,'ne5'] =e= theta['np4']*uc[i,j,'ne3']
                                        - 세타['np5']*uc[i,j,'ne5'];

v.fx['nh1',s] = bc(s);

모델 피넨/모두/;

$if set workSpace pinene.workSpace = %workSpace%

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