설명
이 모델은 외부 리소스를 사용하여 가속화하는 방법을 탐구합니다.
수입을 늘리기 위해 추가 자원을 공급함으로써 개발
그리고 투자. 터미널 조건은 텍스트에서 재현할 수 없습니다.
모든 단위는 수십억 1965루피입니다.
대형 모델 유형 :LP
메인 파일 : pak.gms
$title 성장과 원조의 최적 패턴 (PAK,SEQ=34)

$onText
이 모델은 가속화를 위해 외부 리소스의 사용을 탐구합니다.
수입을 늘리기 위해 추가 자원을 공급함으로써 개발
그리고 투자. 터미널 조건은 텍스트에서 재현할 수 없습니다.
모든 단위는 수십억 1965루피입니다.

Chenery, H B 및 Macewan, A, 9장: 최적의 성장 패턴
그리고 원조. Chenery, H B, Ed, 구조 변화 및 개발
정책. 옥스포드 대학 출판부, 뉴욕 및 옥스퍼드, 1979.

키워드: 선형 계획법, 성장 패턴, 경제 발전, 자본
          투자, 경제 성장 모델, 자본 유입, 금융 투자,
          경제적 자원
$offText

$sTitle 기본 데이터
세트
   테 '계획 기간 연장' / 1962*1985 /
   t(te) '계획기간' / 1963*1985 /
   j '섹터' / 비거래, 거래 /;

스칼라
   fbb '외국 원조 1962' / 1.183 /
   sb '구해내는 1962' / 3.381 /
   tib '총 투자액 1962' / 4.564 /
   mb '수입 1962' / 3.743 /
   eb '수출 1962' / 2.559 /
   gnpb 'gnp 1962' / 37.380 /
   cb '소비 1962' / 33.999 /
   rho '할인율' / .08 /
   r '포스트 플랜 할인' / .10 /
   g '계획 이후 성장' / .073 /
   가마 '외국자본비용' / 2.0 /
   d '계획 후 가중치' / 1.0 /
   alpha '한계 저축률' / .24 /
   mgnp 'gnp의 한계수입율' / .10 /
   mi '투자 한계수입율' / .35 /
   p '인구 증가' / .025 /
   베타 '투자의 최대 성장' / .13 /
   ee '수출 증가' / .049 /
   q '구호 배급' / .5 /
   num '지원 없이 보낸 기간' / 4 /;

매개변수
   k(j) '자본생산비율' / 비거래 3.0, 거래 4.5 /
   delt(t) '할인 요소'
   dis '포스트 호라이즌 단점에 대한 할인'
   vb(j) '기준 연도 출력'
   e(t) '수출';

e(t) = eb*(1 + ee)**ord(t);
디스플레이 e;

delt(t) = (1 + rho)**(-ord(t));
디스플레이 델타;

dis = (1 + r)**(-카드(t))*(1 - 알파)*(1 + g)/(r - g);
디스플레이 디스플레이;

vb("거래되지 않음") = gnpb;

$sTitle 모델 정의
변수
   gnp(t) '국민 총생산'
   v(t,j) '순 출력'
   ti(te) '총 투자'
   i(te,j) '투자'
   ks(te,j) '자본금'
   s(t) '총 절감액'
   f(t) '순자본유입'
   FB '총 할인 지원'
   m(t) '전통적인 수입품'
   c(te) '소비'
   w '복지';

양수 변수 v, i, s;

방정식
   gnpd(t) 'gnp 정의'
   invd(t) '투자 정의'
   invt(te) '투자 총액'
   tgap(t) '무역 격차'
   incd(t) '국민소득 정의'
   capb(t,j) '용량 균형'
   kbal(te,j) '자본잔고'
   savl(t) '최대 절감액'
   impl(t) '최소 수입'
   invu(te) '투자 상한'
   invl(te) '투자 하한값'
   conl(te) '소비 하한'
   fup (t) 'f(t)의 상한'
   taid '총 원조 정의'
   wdef '복지 정의';

gnpd(t).. gnp(t) =e= sum(j, v(t,j));

invd(t)..ti(t) =e= s(t) + f(t);

invt(te).. ti(te) =e= sum(j, i(te,j));

tgap(t).. f(t) =e= m(t) - e(t) - v(t,"거래됨");

incd(t).. gnp(t) =e= c(t) + ti(t) - f(t);

capb(t,j).. v(t,j) =l= vb(j) + 1/k(j)*ks(t,j);

kbal(te+1,j)..ks(te+1,j) =e= ks(te,j) + i(te,j);

savl(t).. s(t) =l= sb + 알파*(gnp(t) - gnpb);

impl(t).. m(t) =g= mb + mgnp*(gnp(t) - gnpb) + mi*(ti(t) - tib);

invu(te+1)..ti(te+1) =l= (1 + 베타)*ti(te);

invl(te+1)..ti(te+1) =g= ti(te);

conl(te+1).. c(te+1) =g= (1 + p)*c(te);

fup(t).. f(t) =l= q*gnp(t);

taid.. fb =e= sum(t, delt(t)*f(t));

wdef.. w =e= sum(t, delt(t)*c(t)) - gama*fb + d*dis*gnp("1985");

ks.fx("1962",j) = 0;
i.fx("1962","비거래됨") = tib;
i.fx("1962","거래됨") = 0;
c.fx("1962") = cb;
f.up(t) = inf$(card(t) - ord(t) >= num);

모델 pak1 / 모두 /;

lp를 사용하여 w를 최대화하는 pak1을 해결합니다.

매개변수 담당자 '요약 보고서(십억 루피)';
담당자(t,"f") = f.l(t);
담당자(t,"gnp") = gnp.l(t);
담당자(t,"ti") = ti.l(t);
담당자(t,"s") = s.l(t);
담당자(t,"c") = c.l(t);

디스플레이 담당자;