문헌 : 이론적 부분 은 Chang (1988) 및 Achdou et al. (2015) 는 각각 숫자 부분입니다.
모델
1 인당 표기법에서 다음과 같은 확률 적 최적 성장 문제를 고려하십시오.
모든 것이 표준 인 dz 를 제외하고 표준 입니다. 표준 위너 과정의 증가, 즉, Z (t) \ SIM \ mathcal {N} (0, t) . 모집단 성장률은 평균 n 및 분산 \ sigma ^ 2 입니다.
성 최대씨∫∞0이자형− ρ tu ( c ) d티디k = [ f( k ) − ( n − σ2) k - c ] dt - σk d지c ∈ [ 0 , f( k ) ]k ( 0 ) = k0
디지지( t ) ~ N( 0 , t )엔σ2
분석 솔루션
Cobb-Douglas 기술은
에프( k ) = kα,α ∈ ( 0 , 1 )
및 CRRA 유틸리티
u ( c ) = c1 - γ1 - γ,γ> 1.
Hamilton-Jacobi 설정 -벨만 방정식 (HJB-e)
ρ v ( k ) = 최대씨{ c1 - γ1 - γ+ v'( k ) ( kα− ( n − σ2) k − c ) + v′ ′( k ) k2σ22}
첫 번째 주문 조건 (FOC)은
c = v'( k )− 1γ= : π( k )
여기서
π( ⋅ ) 는 정책 기능을 나타냅니다.
FOC를 HJB-e로 교체하십시오.
ρ v ( k ) = v'( k )γ− 1γ1 - γ+ v'( k ) kα− v'( k ) ( n − σ2) k - v'( k )γ− 1γ+ v′ ′( k ) k2σ22.
우리의 기능적 형태 추측 가진 ( POSCH (2009 당량. 41) )
v ( k ) = Ψv ( k )
v ( k ) = Ψ k1 - α γ1 - α γ
여기서 는 상수입니다. 의 제 1 및 제 2 차 미분 값 주어진다
vΨV
V'( k )V′ ′( k )= Ψ k− α γ= - α γΨ k− 1 − α γ.
그러면 HJB-e는 를 읽습니다.
⟺ρ Ψ k1 - α γ1 - α γ= Ψγ− 1γ케이α(1−γ)1 - γ+ Ψkα (1−γ)− ( n −σ2) Ψ k1 - α γ− Ψγ− 1γ케이α ( 1 −γ)− αγΨ k1 - α γσ22케이1 - α γ( ρ1 - α γ+ n − σ2( 1 - α γ2) ) = kα ( 1 −γ)[ 1 + Ψ− 1γγ1 - γ]
다음 조건이 인 경우 최대화 된 HJB-e는 true입니다.
ρ = ( − n + σ2( 1 - α γ2))(1−αγ)∧Ψ=(γ−1γ)−γ
를 로 대체 하여 진정한 가치 함수
v v ( k ) =Ψv
V(k)=(γ−1γ)−γk1 −αγ1−αγ.
따라서 결정론과 확률 론적 가치 함수는 동일해야합니다. 그런 다음 정책 기능을 쉽게 제공 할 수 있습니다 (FOC 및 파생 기능 사용)
π( k ) = ( 1 − 1γ) kα.
이 함수는 에도 의존하지 않습니다 .σ
수치 근사
나는 upwind scheme으로 HJB-e를 풀었다. 오차 허용 오차 . 아래 그림에서 다양한 대한 정책 기능을 플로팅합니다 . 를 들어 나는 진정한 솔루션 (보라색)에 도착한다. 그러나 경우 대략적인 정책 기능은 실제 정책 기능과 다릅니다. 가 의존하지 않기 때문에 어느 것이 맞지 않아야 합니까? ϵ = 1 e - 10σσ→ 0σ> 0π( k )σ
- 근사화 된 정책 기능은 어떤에 대해 동일해야 누구 확인할 수 있습니다 진정한 하나의 독립적이기 때문에, ?σσ