Stochastisches Wachstum in kontinuierlicher Zeit

Literatur: Siehe Chang (1988) für den theoretischen Teil und Achdou et al. (2015) jeweils für den numerischen Teil.

Modell

Betrachten Sie das folgende stochastische Problem des optimalen Wachstums in der Pro-Kopf-Notation.

\begin{aligned} max_{c} \int_{0}^{\infty} e^{- ρ t} u (c) d t \\ s.t. & d k = [f (k) - (n - σ^{2}) k - c] d t - σ k d z \\ c \in [0, f (k)] \\ k (0) = k_{0} \end{aligned}

$\begin{align} &\max_{c}\int^\infty_0 e^{-\rho t}u(c)dt\\ \text{s.t.}~~~& dk = [f(k) - (n-\sigma^2) k - c]dt - \sigma kdz\\ &c\in[0,f(k)]\\ &k(0) = k_0 \end{align}$

d z

$dz$

z (t) \sim N (0, t)

$z(t)\sim\mathcal{N}(0,t)$

n

$n$

σ^{2}

$\sigma^2$

Analytische Lösung

Wir gehen davon aus, dass die Cobb-Douglas-Technologie

\begin{aligned} f (k) = k^{α}, α \in (0, 1) \end{aligned}

$\begin{align} f(k) = k^\alpha,\quad \alpha\in(0,1) \end{align}$

und CRRA-Dienstprogramm Richten Sie den Hamilton-Jacobi ein -Bellman-Gleichung (HJB-e)

\begin{aligned} u (c) = \frac{c^{1 - γ}}{1 - γ}, γ > 1. \end{aligned}

$\begin{align} u(c) = \frac{c^{1-\gamma}}{1-\gamma},\quad \gamma > 1. \end{align}$

\begin{aligned} ρ v (k) = max_{c} {\frac{c^{1 - γ}}{1 - γ} + v^{'} (k) (k^{α} - (n - σ^{2}) k - c) + v^{″} (k) \frac{k^{2} σ^{2}}{2}} \end{aligned}

$\begin{align} \rho v(k) = \max_c\left\{\frac{c^{1-\gamma}}{1-\gamma} + v'(k)(k^\alpha - (n - \sigma^2)k - c) + v''(k)\frac{k^2\sigma^2}{2}\right\} \end{align}$

Die Bedingung erster Ordnung (FOC) lautet

\begin{aligned} c = v^{'} (k)^{- \frac{1}{γ}} =: π (k) \end{aligned}

$\begin{align} c = v'(k)^{-\frac{1}{\gamma}}=:\pi(k) \end{align}$ wobei

π (\cdot)

$\pi(\cdot)$ bezeichnet die Richtlinienfunktion.

Resubstituiere FOC in HJB-e

\begin{aligned} ρ v (k) = \frac{v^{'} (k)^{\frac{γ - 1}{γ}}}{1 - γ} + v^{'} (k) k^{α} - v^{'} (k) (n - σ^{2}) k - v^{'} (k)^{\frac{γ - 1}{γ}} + v^{″} (k) \frac{k^{2} σ^{2}}{2} . \end{aligned}

$\begin{align} \rho v(k) = \frac{v'(k)^{\frac{\gamma-1}{\gamma}}}{1-\gamma} + v'(k)k^\alpha - v'(k)(n - \sigma^2)k - v'(k)^{\frac{\gamma-1}{\gamma}} + v''(k)\frac{k^2\sigma^2}{2}. \end{align}$

Wir schätzen eine funktionale Form von mit ( Posch (2009, Gl. 41) ) $v(k)$

\begin{aligned} v (k) = Ψ \frac{k^{1 - α γ}}{1 - α γ} \end{aligned}

$\begin{align} v(k) = \Psi \frac{k^{1-\alpha\gamma}}{1-\alpha\gamma} \end{align}$

Dabei ist eine Konstante. Die Ableitung erster und zweiter Ordnung von ist gegeben durch $\Psi$ $v$

\begin{aligned} v^{'} (k) & = Ψ k^{- α γ} \\ v^{″} (k) & = - α γ Ψ k^{- 1 - α γ} . \end{aligned}

$\begin{align} v'(k) &= \Psi k^{-\alpha\gamma}\\ v''(k) &= -\alpha\gamma\Psi k^{-1-\alpha\gamma}. \end{align}$

Das HJB-e liest dann

\begin{aligned} ρ Ψ \frac{k^{1 - α γ}}{1 - α γ} = \frac{Ψ^{\frac{γ - 1}{γ}} k^{α (1 - γ)}}{1 - γ} + Ψ k^{α (1 - γ)} - (n - σ^{2}) Ψ k^{1 - α γ} - Ψ^{\frac{γ - 1}{γ}} k^{α (1 - γ)} - α γ Ψ k^{1 - α γ} \frac{σ^{2}}{2} \\ ⟺ & k^{1 - α γ} (\frac{ρ}{1 - α γ} + n - σ^{2} (1 - \frac{α γ}{2})) = k^{α (1 - γ)} [1 + Ψ^{- \frac{1}{γ}} \frac{γ}{1 - γ}] \end{aligned}

$\begin{align} &\rho \Psi\frac{k^{1-\alpha\gamma}}{1-\alpha\gamma} = \frac{\Psi^{\frac{\gamma-1}{\gamma}}k^{\alpha(1-\gamma)}}{1-\gamma} + \Psi k^{\alpha(1-\gamma)} - (n-\sigma^2) \Psi k^{1-\alpha\gamma} - \Psi^{\frac{\gamma-1}{\gamma}} k^{\alpha(1-\gamma)} - \alpha\gamma\Psi k^{1-\alpha\gamma}\frac{\sigma^2}{2}\\[2mm] \Longleftrightarrow \quad & k^{1-\alpha\gamma}\left(\frac{\rho}{1-\alpha\gamma} + n - \sigma^2\left(1 - \frac{\alpha\gamma}{2}\right)\right) = k^{\alpha(1-\gamma)}\left[1+\Psi^{-\frac{1}{\gamma}}\frac{\gamma}{1-\gamma}\right] \end{align}$

Das maximierte HJB-e ist wahr, wenn die folgenden Bedingungen

\begin{aligned} ρ = (- n + σ^{2} (1 - \frac{α γ}{2})) (1 - α γ) \land Ψ = {(\frac{γ - 1}{γ})}^{- γ} \end{aligned}

$\begin{align} \rho = \left(-n + \sigma^2\left(1 - \frac{\alpha\gamma}{2}\right)\right)(1-\alpha\gamma)\quad \wedge \quad \Psi = \left(\frac{\gamma-1}{\gamma}\right)^{-\gamma} \end{align}$

Resubstituiere in was schließlich die wahre Wertefunktion ergibt. $\Psi$ $v$

\begin{aligned} v (k) = {(\frac{γ - 1}{γ})}^{- γ} \frac{k^{1 - α γ}}{1 - α γ} . \end{aligned}

$\begin{align} v(k) = \left(\frac{\gamma-1}{\gamma}\right)^{-\gamma} \frac{k^{1-\alpha\gamma}}{1-\alpha\gamma}. \end{align}$

Wie kommt es, dass nicht von abhängt ? $v$ $\sigma$

Die deterministische und die stochastische Wertfunktion müssen also gleich sein. Die Richtlinienfunktion ist dann leicht gegeben durch (benutze FOC und Ableitung der Wertfunktion)

\begin{aligned} π (k) = (1 - \frac{1}{γ}) k^{α} . \end{aligned}

$\begin{align} \pi(k) = \left(1-\frac{1}{\gamma}\right)k^\alpha. \end{align}$

Beachten Sie, dass diese Funktion auch nicht von abhängt . $\sigma$

Numerische Approximation

Ich habe den HJB-e mit einem Gegenwind gelöst. Fehlertoleranz . In der folgenden Abbildung zeichne ich die Richtlinienfunktion für das Variieren von $\epsilon=1e-10$ $\sigma$ . Für ich die wahre Lösung (lila). Aber für weicht die angenäherte Richtlinienfunktion von der wahren ab. Was sollte nicht der Fall sein, da nicht von abhängt , oder? $\sigma\to 0$ $\sigma>0$ $\pi(k)$ $\sigma$

Kann jemand bestätigen, dass die ungefähren Richtlinienfunktionen für alle dieselben sein sollten? $\sigma$ , da das wahre von unabhängig ist ? $\sigma$

academic-graduate continuous-time optimal-control numerical-methods ahnungslos
quelle

Was mich hier stört, ist die erste "iff" -Bedingung nach dem Schreiben "das maximierte HJB-e ist wahr, wenn die folgenden Bedingungen gelten": Dies ist eine sehr spezifische Gleichheitsrelation , die zwischen allen Parametern der Modell- Präferenz-Parameter gelten muss. Bevölkerungswachstum, Kapitalproduktivität und Volatilität. Ich frage mich: Können wir wirklich mit erratenen Funktionen arbeiten, deren Gültigkeit von einer so engen Bedingung abhängt, die von den Parametern abhängt?

Alecos Papadopoulos

Nun, hier lege ich tatsächlich als Funktion der vier verbleibenden Parameter fest. Die Gleichung ist also immer wahr, wenn zusätzlich gilt. Ich frage mich: Gibt es eine Regel, nach der das Erraten einer Funktion nicht zulässig ist? Ich meine, wir sind daran interessiert, die wahre Lösung zu finden, und unter bestimmten Bedingungen erhalten wir die wahre Lösung. Ich bin nicht sicher, was Sie hier aus theoretischer Sicht stört? Sicher, es kann die empirische Arbeit einschränken, aber darum geht es hier nicht. Wir sind eher daran interessiert, die HJBe zu lösen, und das kann getan werden. Wenn ein Empiriker (1/2)

ρ = ρ (α, γ, n, σ)

$\rho = \rho(\alpha, \gamma, n, \sigma)$

ρ > 0

$\rho > 0$

ahnungslos

schätzt und wir stellen fest, dass die Bedingung verletzt ist, dann können wir das Modell ablehnen. Die Lösung bleibt jedoch grundsätzlich wahr. (2/2)

{α, γ, n, ρ, σ}

$\{\alpha, \gamma,n,\rho,\sigma\}$

ρ = . . . .

$\rho = ....$

ahnungslos

Mein Anliegen ist nicht die empirische Validität. Ich frage mich, inwieweit die genaue Einschätzung der Funktionsform der Wertfunktion von dieser Beziehung zwischen den Parametern abhängt. Wenn wir ohne Bezugnahme auf empirische Daten annehmen, dass die Beziehung nicht zutrifft, was dann? Sollten wir eine Wertefunktion erraten, die in nicht einmal exponentiell ist , oder würde es ausreichen, die exponentielle Struktur beizubehalten, aber verschiedene Methoden zu versuchen, um die Parameter darin einzuschließen? (Ich beschäftige mich übrigens auch mit Ihrer Hauptfrage, da diese Diskussion wahrscheinlich eher am Rande stattfindet.)

k

$k$

Alecos Papadopoulos,

Sind Sie sicher, dass das Optimierungsproblem korrekt angegeben wurde? Es gibt zum Beispiel keine Erwartung, die an ? Wie es jetzt angegeben ist, nehmen und daher wahrscheinlich irgendeinen Wert an, wenn der Wiener-Prozess .

f (k)

$f(k)$

k

$k$

f (k)

$f(k)$

z

$z$

Hans

Stochastisches Wachstum in kontinuierlicher Zeit

Modell

Analytische Lösung

Numerische Approximation

Antworten: