Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Ramsey Growth in Discrete and Continuous Time

Authors
Affiliations
Johns Hopkins University
Econ-ARK
Johns Hopkins University
Econ-ARK

This section solves a continuous-time version of the Ramsey/Cass-Koopmans (RCK) model using the Hamiltonian method, and shows the relationship between that method and the discrete-time approach.

The problem is to choose a path of consumption per capita c\cons from the present moment (arbitrarily called time 0) into the infinite future, {c}0∞\{c\}_{0}^{\infty}, that solves the problem

max⁡{c}0∞∫0∞u(c)e−ϑt\max_{\{c\}_{0}^{\infty}} \int_{0}^{\infty} \uFunc(\cons) e^{-\timeRate t}

subject to

k˙=f(k)−c−(ξ+δ)kk>0 ∀ t\begin{gathered}\begin{aligned} \dot{\kap} & = \fFunc(\kap) - \cons - (\popGro+\depr)\kap \\ \kap & > 0 ~\forall~t \end{aligned}\end{gathered}

where ϑ\timeRate is the time preference rate, ξ\popGro is the population growth rate, and δ\depr is the depreciation rate. (In continuous time, we think of all variables as implicitly being a function of time, but it is cumbersome to write, e.g., c(t)\cons(t) everywhere, so the time argument is omitted; we are also thinking of the initial value of capital at date 0 as being a ‘given’ in the problem, so that k(0)=∙\kap(0)=\bullet for some specific value of ∙\bullet).

To emphasize the similarity between the continuous-time and the discrete-time solutions where we have typically used the roman VV to denote value, for the continuous-time problem we define ‘curly’ value as a function of the initial level of capital as V(k)\mathcal{V}(\kap).

The current-value (discounted) Hamiltonian is

H(k,c,λ)=u(c)+(f(k)−c−(ξ+δ)k)λ\mathcal{H}(\kap,\cons,\lambda) = \uFunc(\cons) + (\fFunc(\kap)-c-(\popGro+\depr)\kap)\lambda

where k\kap is the state variable, c\cons is the control variable, and λ\lambda is the costate variable.

λ\lambda is the continuous-time equivalent of a Lagrange multiplier, so its value should be equivalent to the value of relaxing the corresponding constraint by an infinitesimal amount. But the constraint in question is the capital-accumulation constraint. Thus λ\lambda should be equal to the value of having a tiny bit more capital, V(k+Δk)−V(k)Δk\frac{\mathcal{V}(\kap+\Delta k)-\mathcal{V}(\kap)}{\Delta k}. In other words, you can think of λ=V′(k)\lambda=\mathcal{V}'(\kap).

The first necessary Hamiltonian condition for optimality is

∂H∂c=0u′(c)=λu′(c)=V′(k).\begin{gathered}\begin{aligned} \frac{\partial \mathcal{H}}{\partial c} & = 0 \\ \uFunc^{\prime}(\cons) & = \lambda \\ \uFunc^{\prime}(\cons) & = \mathcal{V}'(\kap) . \end{aligned}\end{gathered}

Note the similarity between (4) and the result we usually obtain by using the Envelope theorem in the discrete-time problem,

u′(ct)=V′(kt).\begin{gathered}\begin{aligned} \uFunc^{\prime}(\cons_{t}) & = \VFunc^{\prime}(\kap_{t}) . \end{aligned}\end{gathered}

Thus, you can use the intuition you (should have) developed by now about why the marginal utility of consumption should be equal to the marginal value of extra resources to understand this Hamiltonian optimality condition.

The second necessary condition is

λ˙=ϑλ−λ(f′(k)−(ξ+δ))(λ˙λ)=ϑ−(f′(k)−(ξ+δ))\begin{gathered}\begin{aligned} \dot{\lambda} & = \timeRate\lambda - \lambda(\fFunc^{\prime}(\kap)-(\popGro+\depr)) \\ \left(\frac{\dot{\lambda}}{\lambda}\right) & = \timeRate-(\fFunc^{\prime}(\kap)-(\popGro+\depr)) \end{aligned}\end{gathered}

which expresses the growth rate of λ\lambda at an annual rate (because the interest rate r\rfree and time preference rate ϑ\timeRate are measured at an annual rate).

To interpret this in terms of our discrete-time model, begin with the condition

V′(kt)=RβV′(kt+1).\begin{gathered}\begin{aligned} \VFunc^{\prime}(\kap_{t}) & = \Rfree\Discount \VFunc^{\prime}(\kap_{t+1}). \end{aligned}\end{gathered}

The final necessary condition is just that the accumulation equation for capital is satisfied,

k˙=f(k)−c−(ξ+δ)k.\dot{\kap} = \fFunc(\kap)-c-(\popGro+\depr)\kap.

This is the continuous-time equivalent of what we have previously called the Dynamic Budget Constraint.

Up to now in this course we haven’t thought very much about what the time period is. Generally, we have expressed things in terms of yearly rates, so that for example we might choose R=1.04\Rfree=1.04 and β=1/(1+ϑ)=1/(1.04)\Discount=1/(1+\timeRate)=1/(1.04) to represent an interest rate of 4 percent and a discount rate of 4 percent.

One of the attractive features of the time-consistent model we have been using is that it generates self-similar behavior as the time interval is changed. Thus if we wanted to solve a quarterly version of the model we would choose R=1.01\Rfree=1.01 and β=1/1.01\Discount=1/1.01 and it would imply consumption of almost exactly 1/4 of the amount implied by the annual model, so that four quarters of such behavior would aggregate to the prediction of the annual model.

To put this in the most general form, suppose R\Rfree and β\Discount correspond to ‘annual rate’ values and we want to divide the year into mm periods. Then the appropriate interest rate and discount factor on a per-period basis would be R1/m\Rfree^{1/m} and β1/m\Discount^{1/m}. Thus the discrete-time equation could be rewritten

V′(kt)=R1/mβ1/mV′(kt+1)\begin{gathered}\begin{aligned} \VFunc^{\prime}(\kap_{t}) & = \Rfree^{1/m}\Discount^{1/m} \VFunc^{\prime}(\kap_{t+1}) \end{aligned}\end{gathered}

where the time interval is now 1/m1/mth of a year (e.g. if mm=52, we’re talking weekly, so that period t+1t+1 is one week after period tt). Now we can use our old friend, the fact that ez≈1+ze^{z} \approx 1+z, to note that this is approximately

V′(kt)≈[er]1/m[e−ϑ]1/mV′(kt+1)=e(1/m)re−(1/m)ϑV′(kt+1)=e(1/m)(r−ϑ)V′(kt+1)=e(1/m)(r−ϑ)(V′(kt)+ΔV′(kt+1))1=e(1/m)(r−ϑ)(V′(kt)+ΔV′(kt+1)V′(kt))e(1/m)(ϑ−r)=(V′(kt)+ΔV′(kt+1)V′(kt))V′(kt)(e(1/m)(ϑ−r)−1⏟≈1+(1/m)(ϑ−r)−1)=ΔV′(kt+1)ΔV′(kt+1)V′(kt)≈(1/m)(ϑ−r)mΔV′(kt+1)V′(kt)≈(ϑ−r).\begin{gathered}\begin{aligned} \VFunc^{\prime}(\kap_{t}) & \approx [e^{\rfree}]^{1/m}[e^{-\timeRate}]^{1/m} \VFunc^{\prime}(\kap_{t+1}) \\ & = e^{(1/m)\rfree}e^{-(1/m)\timeRate} \VFunc^{\prime}(\kap_{t+1}) \\ & = e^{(1/m)(\rfree-\timeRate)} \VFunc^{\prime}(\kap_{t+1}) \\ & = e^{(1/m)(\rfree-\timeRate)} (\VFunc^{\prime}(\kap_{t})+\Delta \VFunc^{\prime}(\kap_{t+1})) \\ 1 & = e^{(1/m)(\rfree-\timeRate)} \left(\frac{\VFunc^{\prime}(\kap_{t})+\Delta \VFunc^{\prime}(\kap_{t+1})}{\VFunc^{\prime}(\kap_{t})}\right) \\ e^{(1/m)(\timeRate-\rfree)} & = \left(\frac{\VFunc^{\prime}(\kap_{t})+\Delta \VFunc^{\prime}(\kap_{t+1})}{\VFunc^{\prime}(\kap_{t})}\right) \\ \VFunc^{\prime}(\kap_{t})(\underbrace{e^{(1/m)(\timeRate - \rfree)}-1}_{\approx 1+(1/m)(\timeRate-\rfree)-1}) & = \Delta \VFunc^{\prime}(\kap_{t+1}) \\ \frac{\Delta \VFunc^{\prime}(\kap_{t+1})}{\VFunc^{\prime}(\kap_{t})} & \approx (1/m)(\timeRate - \rfree) \\ \frac{m\Delta \VFunc^{\prime}(\kap_{t+1})}{\VFunc^{\prime}(\kap_{t})} & \approx (\timeRate - \rfree) . \end{aligned}\end{gathered}

We defined the interest rate and time preference rate on an annual basis, but the time interval between tt and t+1t+1 is only (1/m)(1/m)th of a year. Thus mΔV′(kt+1)m \Delta \VFunc^{\prime}(\kap_{t+1}) expresses the speed of change in V′(kt)\VFunc^{\prime}(\kap_{t}) at an annual rate.

Now, note that since the effective interest rate in this model is f′(k)−(ξ+δ)\fFunc^{\prime}(\kap)-(\popGro+\depr), equation (10) is basically the same as (6) since λ=V′(k)\lambda = \mathcal{V}'(\kap) and mΔV′(kt+1)=(d/dt)V′(k)=λ˙m \Delta \VFunc^{\prime}(\kap_{t+1}) = (d/dt)\mathcal{V}'(\kap) = \dot{\lambda}. Hence, the second optimality condition in the Hamiltonian optimization method is basically equivalent to the condition V′(kt)=RβV′(kt+1)\VFunc^{\prime}(\kap_{t})=\Rfree\Discount \VFunc^{\prime}(\kap_{t+1}) from the discrete-time optimization method!

The final required condition (the transversality constraint) is

lim⁡t→∞λke−ϑt=0\begin{gathered}\begin{aligned} \lim_{t \rightarrow \infty} \lambda \kap e^{-\timeRate t}& = 0 \end{aligned}\end{gathered}

The translation of this into the discrete-time model is

lim⁡t→∞βtu′(ct)kt=0.\lim_{t \rightarrow \infty} \Discount^{t} \uFunc^{\prime}(\cons_{t}) \kap_{t} = 0.

Consider the simple model with a constant gross interest rate R\Rfree and CRRA utility. In that model, recall that ct+1=(Rβ)1/ρct\cons_{t+1}= (\Rfree\Discount)^{1/\CRRA} \cons_{t}. Thus considered from time zero (12) becomes

lim⁡t→∞βt(c0((Rβ)1/ρ)t)−ρkt=0=c0−ρβt[(Rβ)t/ρ]−ρkt=c0−ρβtβ−tR−tkt=c0−ρR−tkt→lim⁡t→∞R−tkt=0.\begin{gathered}\begin{aligned} \lim_{t \rightarrow \infty} \Discount^{t} (\cons_{0} ((\Rfree\Discount)^{1/\CRRA})^{t})^{-\CRRA} \kap_{t} & = 0 \\ & = \cons_{0}^{-\CRRA} \Discount^{t}[(\Rfree\Discount)^{t/\CRRA}]^{-\CRRA} \kap_{t} \\ & = \cons_{0}^{-\CRRA} \Discount^{t} \Discount^{-t} \Rfree^{-t} \kap_{t} \\ & = \cons_{0}^{-\CRRA} \Rfree^{-t} \kap_{t} \\ \rightarrow \lim_{t \rightarrow \infty} \Rfree^{-t} \kap_{t} & = 0 . \end{aligned}\end{gathered}

What this says is that you cannot behave in such a way that you expect kt\kap_{t} to grow faster than the interest rate forever.[1]

This is the infinite-horizon version of the intertemporal budget constraint.

Among the infinite number of time paths of c\cons and k\kap that will satisfy the first order conditions above, only one will also satisfy this transversality constraint, because all the others imply a violation of the intertemporal budget constraint.

Now differentiate (4) with respect to time

c˙u′′(c)=λ˙\begin{gathered}\begin{aligned} \dot{\cons} \uFunc^{\prime\prime}(\cons) & = \dot{\lambda} \end{aligned}\end{gathered}

and substitute this into equation (6) to get

c˙u′′(c)u′(c)=(ϑ−(f′(k)−(ξ+δ)))c˙=−u′(c)u′′(c)(f′(k)−(ξ+δ)−ϑ)\begin{gathered}\begin{aligned} \frac{\dot{\cons} \uFunc^{\prime\prime}(\cons)}{\uFunc^{\prime}(\cons)} & = (\timeRate-(\fFunc^{\prime}(\kap)-(\popGro+\depr))) \\ \dot{\cons} & = -\frac{\uFunc^{\prime}(\cons)}{\uFunc^{\prime\prime}(\cons)} (\fFunc^{\prime}(\kap)-(\popGro+\depr)-\timeRate) \end{aligned}\end{gathered}

using the fact derived earlier that for a CRRA utility function u(c)=c1−ρ/(1−ρ),−u′′(c)c/u′(c)=ρ\uFunc(\cons)=c^{1-\CRRA}/(1-\CRRA), -\uFunc^{\prime\prime}(\cons) c/\uFunc^{\prime}(\cons) = \CRRA,, this becomes

=(c/ρ)(f′(k)−(ξ+δ)−ϑ)c˙/c=ρ−1(f′(k)−(ξ+δ)−ϑ)\begin{gathered}\begin{aligned} & = (\cons/\CRRA) (\fFunc^{\prime}(\kap)-(\popGro+\depr)-\timeRate) \\ \dot{\cons}/\cons & = \CRRA^{-1}(\fFunc^{\prime}(\kap)-(\popGro+\depr)-\timeRate) \end{aligned}\end{gathered}
Footnotes
  1. Note that this also rules out negative kt\kap_{t} values that grow faster than the interest rate.