참조
카테고리 : GAMS NOA 라이브러리
메인파일 : fiat.gms
$onText
스파크 점화 엔진의 안정성 마진 분석
피아트 데드라.
참고자료:
B.R. Barmish, 선형 시스템의 견고성을 위한 새로운 도구.
McMillan 출판사, 뉴욕, 1994.
M. Abate, B. Barmish, C. Murillo-Sanchez, R. Tempo, 응용 프로그램
스파크 점화 엔진의 강력한 안정성 분석을 위한 몇 가지 새로운 도구:
사례 연구. IEEE 트랜스. 대조 시스템 tech., vol.2, 1994, pp. 22.
Neculai Andrei, "모델, 테스트 문제 및 응용 프로그램
수학적 프로그래밍". 기술 언론, 부쿠레슈티, 2003.
신청서 A41, 407페이지.
Floudas, C.A., Pardalos, P.M., et al. "테스트 문제 핸드북
로컬 및 글로벌 최적화". Kluwer 학술 출판사, Dordrecht,
1999.
문제 7.3.6. 103페이지, 문제 16번을 테스트하세요.
$offText
변수
q1, q2, q3, q4, q5, q6, q7
w 주파수
k 안정성 마진
a0, a1, a2, a3, a4, a5, a6, a7
objval 목적 함수 변수;
무료 변수 objval;
방정식
f 목적 함수
g1
g2
b1l, b1u
비투엘, 비투유
b3l, b3u
b4l, b4u
b5l, b5u
b6l, b6u
b7l, b7u
ga0, ga1, ga2, ga3, ga4, ga5, ga6, ga7;
f .. objval =e=k;
g1 .. -a6*POWER(w,6) + a4*POWER(w,4) - a2*POWER(w,2) + a0 =e= 0;
g2 .. a7*POWER(w,6) - a5*POWER(w,4) + a3*POWER(w,2) - a1 =e= 0;
b1l .. 3.4329-1.02721*k =l= q1;
b1u .. q1 =l= 3.4320+1.02721*k;
b2l .. 0.1627-0.06*k =l= q2;
b2u .. q2 =l= 0.1627+0.06*k;
b3l .. 0.1139-0.0782*k =l= q3;
b3u .. q3 =l= 0.1139+0.0782*k;
b4l .. 1.2539-0.3068*k =l= q4;
b4u .. q4 =l= 1.2539+0.3068*k;
b5l .. 0.0208-0.0108*k =l= q5;
b5u .. q5 =l= 0.0208+0.08*k;
b6l .. 5.0247-2.4715*k =l= q6;
b6u .. q6 =l= 5.0247+2.4715*k;
b7l .. 1.0-2*k =l= q7;
b7u .. q7 =l= 1.0+2*k;
ga0 .. a0 =e= 6.82079e-05슬롯 무료체험1슬롯 무료체험3*POWER(q4,2) + 6.82079e-05슬롯 무료체험1슬롯 무료체험2슬롯 무료체험4슬롯 무료체험5;
ga1 .. a1 =e= 0.00076176*전력(q2,2)*전력(q5,2) +
0.00076176*파워(q3,2)*파워(q4,2) +
0.000402141슬롯 무료체험1슬롯 무료체험2*POWER(q5,2) +
0.00337606슬롯 무료체험1슬롯 무료체험3*POWER(q4,2) +
6.82079e-05슬롯 무료체험1슬롯 무료체험4슬롯 무료체험5 + 0.00051612*POWER(q2,2)슬롯 무료체험5슬롯 무료체험6 +
0.00337606슬롯 무료체험1슬롯 무료체험2슬롯 무료체험4슬롯 무료체험5 + 6.82079e-05슬롯 무료체험1슬롯 무료체험2슬롯 무료체험4슬롯 무료체험7 +
6.28987e-05슬롯 무료체험1슬롯 무료체험2슬롯 무료체험5슬롯 무료체험6 + 0.000402141슬롯 무료체험1슬롯 무료체험3슬롯 무료체험4슬롯 무료체험5 +
6.28987e-05슬롯 무료체험1슬롯 무료체험3슬롯 무료체험4슬롯 무료체험6 + 0.00152352슬롯 무료체험2슬롯 무료체험3슬롯 무료체험4슬롯 무료체험5 +
0.00051612슬롯 무료체험2슬롯 무료체험3슬롯 무료체험4슬롯 무료체험6;
ga2 .. a2 =e= 0.000402141슬롯 무료체험1*POWER(q5,2) + 0.00152352슬롯 무료체험2*POWER(q5,2) +
0.0552*전력(q2,2)*전력(q5,2) +
0.0552*전력(q3,2)*전력(q4,2) +
0.0189477슬롯 무료체험1슬롯 무료체험2*POWER(q5,2) + 0.034862슬롯 무료체험1슬롯 무료체험3*POWER(q4,2) +
0.00336706슬롯 무료체험1슬롯 무료체험4슬롯 무료체험5 + 6.82079e-05슬롯 무료체험1슬롯 무료체험4슬롯 무료체험7 +
6.28987e-05슬롯 무료체험1슬롯 무료체험5슬롯 무료체험6 + 0.00152352슬롯 무료체험3슬롯 무료체험4슬롯 무료체험5 +
0.00051612슬롯 무료체험3슬롯 무료체험4슬롯 무료체험6 - 0.00234048*POWER(q3,2)슬롯 무료체험4슬롯 무료체험6 +
0.034862슬롯 무료체험1슬롯 무료체험2슬롯 무료체험4슬롯 무료체험5 + 0.0237398*전력(q2,2)슬롯 무료체험5슬롯 무료체험6 +
0.00152352*파워(q2,2)슬롯 무료체험5슬롯 무료체험7 +
0.00051612*파워(q2,2)슬롯 무료체험6슬롯 무료체험7 +
0.00336706슬롯 무료체험1슬롯 무료체험2슬롯 무료체험4슬롯 무료체험7 + 0.00287416슬롯 무료체험1슬롯 무료체험2슬롯 무료체험5슬롯 무료체험6 +
0.000804282슬롯 무료체험1슬롯 무료체험2슬롯 무료체험5슬롯 무료체험7 + 6.28987e-05슬롯 무료체험1슬롯 무료체험2슬롯 무료체험6슬롯 무료체험7 +
0.0189477슬롯 무료체험1슬롯 무료체험3슬롯 무료체험4슬롯 무료체험5 + 0.00287416슬롯 무료체험1슬롯 무료체험3슬롯 무료체험4슬롯 무료체험6 +
0.000402141슬롯 무료체험1슬롯 무료체험3슬롯 무료체험4슬롯 무료체험7 + 0.1104슬롯 무료체험2슬롯 무료체험3슬롯 무료체험4슬롯 무료체험5 +
0.0237398슬롯 무료체험2슬롯 무료체험3슬롯 무료체험4슬롯 무료체험6 + 0.00152352슬롯 무료체험2슬롯 무료체험3슬롯 무료체험4슬롯 무료체험7 -
0.00234048슬롯 무료체험2슬롯 무료체험3슬롯 무료체험5슬롯 무료체험6 + 0.00103224슬롯 무료체험2슬롯 무료체험5슬롯 무료체험6;
ga3 .. a3 =e= 0.189477슬롯 무료체험1*POWER(q5,2) + 0.1104슬롯 무료체험2*POWER(q5,2) +
0.00051612슬롯 무료체험5슬롯 무료체험6 + 파워(q2,2)*파워(q5,2) +
0.00076176*파워(q2,2)*파워(q7,2) +
전력(q3,2)*전력(q4,2) +
0.1586슬롯 무료체험1슬롯 무료체험2*POWER(q5,2) + 0.000402141슬롯 무료체험1슬롯 무료체험2*POWER(q7,2) +
0.0872슬롯 무료체험1슬롯 무료체험3*POWER(q4,2) + 0.034862슬롯 무료체험1슬롯 무료체험4슬롯 무료체험5 +
0.00336706슬롯 무료체험1슬롯 무료체험4슬롯 무료체험7 + 0.00287416슬롯 무료체험1슬롯 무료체험5슬롯 무료체험6 +
6.28987e-05슬롯 무료체험1슬롯 무료체험6슬롯 무료체험7 + 0.00103224슬롯 무료체험2슬롯 무료체험6슬롯 무료체험7 +
0.1104슬롯 무료체험3슬롯 무료체험4슬롯 무료체험5 +
0.0237398슬롯 무료체험3슬롯 무료체험4슬롯 무료체험6 + 0.00152352슬롯 무료체험3슬롯 무료체험4슬롯 무료체험7 -
0.00234048슬롯 무료체험3슬롯 무료체험5슬롯 무료체험6 +
0.1826*전력(q2,2)슬롯 무료체험5슬롯 무료체험6 + 0.1104*전력(q2,2)슬롯 무료체험5슬롯 무료체험7 +
0.0237398*파워(q2,2)슬롯 무료체험6슬롯 무료체험7 - 0.0848*파워(q3,2)슬롯 무료체험4슬롯 무료체험6 +
0.0872슬롯 무료체험1슬롯 무료체험2슬롯 무료체험4슬롯 무료체험5 + 0.034862슬롯 무료체험1슬롯 무료체험2슬롯 무료체험4슬롯 무료체험7 +
0.0215658슬롯 무료체험1슬롯 무료체험2슬롯 무료체험5슬롯 무료체험6 + 0.0378954슬롯 무료체험1슬롯 무료체험2슬롯 무료체험5슬롯 무료체험7 +
0.00287416슬롯 무료체험1슬롯 무료체험2슬롯 무료체험6슬롯 무료체험7 + 0.1586슬롯 무료체험1슬롯 무료체험3슬롯 무료체험4슬롯 무료체험5 +
0.0215658슬롯 무료체험1슬롯 무료체험3슬롯 무료체험4슬롯 무료체험6 + 0.0189477슬롯 무료체험1슬롯 무료체험3슬롯 무료체험4슬롯 무료체험7 +
2슬롯 무료체험2슬롯 무료체험3슬롯 무료체험4슬롯 무료체험5 + 0.1826슬롯 무료체험2슬롯 무료체험3슬롯 무료체험4슬롯 무료체험6 + 0.1104슬롯 무료체험2슬롯 무료체험3슬롯 무료체험4슬롯 무료체험7 -
0.0848슬롯 무료체험2슬롯 무료체험3슬롯 무료체험5슬롯 무료체험6 - 0.00234048슬롯 무료체험2슬롯 무료체험3슬롯 무료체험6슬롯 무료체험7 +
0.00076176*파워(q5,2) + 0.0474795슬롯 무료체험2슬롯 무료체험5슬롯 무료체험6 +
0.000804282슬롯 무료체험1슬롯 무료체험5슬롯 무료체험7 + 0.00304704슬롯 무료체험2슬롯 무료체험5슬롯 무료체험7;
ga4 .. a4 =e= 0.1586슬롯 무료체험1*POWER(q5,2) + 0.000402141슬롯 무료체험1*POWER(q7,2) +
2슬롯 무료체험2*전력(q5,2) + 0.00152352슬롯 무료체험2*전력(q7,2) +
0.0237398슬롯 무료체험5슬롯 무료체험6 +
0.00152352슬롯 무료체험5슬롯 무료체험7 + 0.00051612슬롯 무료체험6슬롯 무료체험7 +
0.0552*전력(q2,2)*전력(q7,2) +
0.0189477슬롯 무료체험1슬롯 무료체험2*POWER(q7,2) +
0.0872슬롯 무료체험1슬롯 무료체험4슬롯 무료체험5 + 0.034862슬롯 무료체험1슬롯 무료체험4슬롯 무료체험7 + 0.0215658슬롯 무료체험1슬롯 무료체험5슬롯 무료체험6 +
0.00287416슬롯 무료체험1슬롯 무료체험6슬롯 무료체험7 + 0.0474795슬롯 무료체험2슬롯 무료체험6슬롯 무료체험7 + 2슬롯 무료체험3슬롯 무료체험4슬롯 무료체험5 +
0.1826슬롯 무료체험3슬롯 무료체험4슬롯 무료체험6 + 0.1104슬롯 무료체험3슬롯 무료체험4슬롯 무료체험7 - 0.0848슬롯 무료체험3슬롯 무료체험5슬롯 무료체험6 -
0.00234048슬롯 무료체험3슬롯 무료체험6슬롯 무료체험7 + 2*전력(q2,2)슬롯 무료체험5슬롯 무료체험7 +
0.1826*전력(q2,2)슬롯 무료체험6슬롯 무료체험7 + 0.0872슬롯 무료체험1슬롯 무료체험2슬롯 무료체험4슬롯 무료체험7 +
0.3172슬롯 무료체험1슬롯 무료체험2슬롯 무료체험5슬롯 무료체험7 + 0.0215658슬롯 무료체험1슬롯 무료체험2슬롯 무료체험6슬롯 무료체험7 +
0.1586슬롯 무료체험1슬롯 무료체험3슬롯 무료체험4슬롯 무료체험7 + 2슬롯 무료체험2슬롯 무료체험3슬롯 무료체험4슬롯 무료체험7 - 0.0848슬롯 무료체험2슬롯 무료체험3슬롯 무료체험6슬롯 무료체험7 +
0.0552*POWER(q5,2) + 0.3652슬롯 무료체험2슬롯 무료체험5슬롯 무료체험6 + 0.0378954슬롯 무료체험1슬롯 무료체험5슬롯 무료체험7 +
0.2208슬롯 무료체험2슬롯 무료체험5슬롯 무료체험7;
ga5 .. a5 =e= 0.0189477슬롯 무료체험1*POWER(q7,2) + 0.1104슬롯 무료체험2*POWER(q7,2) +
0.1826슬롯 무료체험5슬롯 무료체험6 + 0.1104슬롯 무료체험5슬롯 무료체험7 + 0.0237398슬롯 무료체험6슬롯 무료체험7 +
파워(q2,2)*파워(q7,2) + 0.1586슬롯 무료체험1슬롯 무료체험2*파워(q7,2) +
0.0872슬롯 무료체험1슬롯 무료체험4슬롯 무료체험7 + 0.0215658슬롯 무료체험1슬롯 무료체험6슬롯 무료체험7 + 0.3652슬롯 무료체험2슬롯 무료체험6슬롯 무료체험7 +
2슬롯 무료체험3슬롯 무료체험4슬롯 무료체험7 - 0.0848슬롯 무료체험3슬롯 무료체험6슬롯 무료체험7 + POWER(q5,2) +
0.00076176*POWER(q7,2) + 0.3172슬롯 무료체험1슬롯 무료체험5슬롯 무료체험7 + 4슬롯 무료체험2슬롯 무료체험5슬롯 무료체험7;
ga6 .. a6 =e= 0.1586슬롯 무료체험1*POWER(q7,2) + 2슬롯 무료체험2*POWER(q7,2) + 2슬롯 무료체험5슬롯 무료체험7 +
0.1826슬롯 무료체험6슬롯 무료체험7 + 0.0552*POWER(q7,2);
ga7 .. a7 =e= POWER(q7,2);
* 경계
슬롯 무료체험1.UP = 3.4329;
슬롯 무료체험2.UP = 0.1627;
슬롯 무료체험3.UP = 0.1139;
슬롯 무료체험4.LO = 0.2539;
슬롯 무료체험5.UP = 0.0208;
슬롯 무료체험6.LO = 2.0247;
슬롯 무료체험7.LO = 1;
w.LO = 0;
w.UP = 10;
k.LO = 0;
k.UP = 10;
* 초기점
q1.1=0.2;
q2.1=0.02;
q3.1=0.1;
q4.1=0.3;
q5.1=0;
q6.1=2;
q7.1=4.5;
w.l=0;
k.l=2;
모델 법정화폐 /ALL/;
NLP를 사용하여 법정화폐를 해결하세요. objval을 최소화하세요.
* 법정화폐 종료