Method Documentation
Expected path
12min
this method implements a mechanism for incorporating upper and lower bounds into the forecast of a saving plan, thus effectively capturing the inherent volatility of the respective plan the calculation is based on the solution bounds derived from stochastic differential equations throughout this documentation, we assume that the initial expected return and risk of the portfolio are denoted as ยต and ฯ, respectively, and are provided as inputs in this context, ฯ
(t) represents the future price projection of a portfolio, characterized by an expected return following a normal distribution with parameters (ยต, ฯ) the variable t denotes the time of projection, while ฯ
0 denotes the initial value of the portfolio at t = 0 here, ฯ
is a stochastic process which satisfies itรดโs lemma \[1] \begin{align } d\upsilon(t) & =\frac{\partial\upsilon}{\partial t}dt+\frac{\partial\upsilon}{\partial w}dw+\frac{1}{2}\frac{\partial^{2}\upsilon}{\partial w^{2}}(dw)^{2} \end{align } where w is a wiener process \[2], also known as a brownian motion, characterized by a normal distribution \mathcal{n}(0, \sqrt{dt}) in this case \left(dw\right)^2 = dt projection with cashflow for a portfolio with initial investment ฯ
0, the solution by \nu(t) = \nu {0} e^{\left(\mu \frac{1}{2}\sigma^{2}\right)t + \sigma w(t)} using itรดโs lemma follows geometric brownian motion d\nu = \mu \nu dt + \sigma \nu dw the differential form dฯ
describes the change of investment value for a given confidence interval z, using the solution, the bounds to ฯ
are given by \begin{align } \upsilon^{u}(t) & =\upsilon {0}e^{\left(\mu \frac{1}{2}\sigma^{2}\right)t+\sigma z\sqrt{t}}\\\\ \upsilon^{l}(t) & =\upsilon {0}e^{\left(\mu \frac{1}{2}\sigma^{2}\right)t \sigma z\sqrt{t}} \end{align } where \begin{align } \upsilon^{u}(t) & \text{the upper bound of the confidence interval for } \upsilon \text{ at time } t, \\\\ \upsilon^{l}(t) & \text{the lower bound of the confidence interval for } \upsilon \text{ at time } t \end{align } projection without cashflow in this section, we consider a portfolio with initial investment ฯ
0 and future cashflows given by \left\\{ (t {c,1}, s {1}), (t {c,2}, s {2}), \ldots, (t {c,n}, s {n}) \right\\} where s i >=0 is the amount of cashflow at each instant time t c,i>=0 to represent these discrete values in continuous time, we use the continuous function of cashflows given by c i(t) = s i \delta {t i}(t) where ฮด is the dirac delta distribution \[3] for this kind of portfolio with initial investment and cashflows, to calculate the evolution of portfolio value, we adapt the method in \[4] where the solution given by \upsilon(t) = \upsilon {0}e^{(\mu \frac{1}{2}\sigma^{2})t+\sigma w(t)} + \sum {i}\int {t {c,i}}^{t}c {i}(\tau)e^{(\mu \frac{1}{2}\sigma^{2})(t \tau)+\sigma(w(t) w(\tau))}d\tau follows d\upsilon = (\mu\upsilon + \sum {i}c {i})dt + \sigma\upsilon dw to calculate the lower and upper bound of the solution for a given confidence interval z, we use numerical approximation as follow the numerical time step is denoted by โt and the index of step is given by c {i}(t) = s {i}\delta {t {i}}(t) the dirac delta is approximated by \delta {t {c,i}}(k\delta t) \approx \begin{cases} \frac{1}{\delta t}, & \text{for } k {c,i} \leq k \leq k {c,i}+1 \\\\ 0, & \text{otherwise} \end{cases} where k {c,i} = \frac{t {c,i}}{\delta t} here, we can denote the set which contains the step index where cashflow occurs \omega {c} = \left\\{ \frac{t {c,1}}{\delta t}, \frac{t {c,2}}{\delta t}, \ldots, \frac{t {c,n}}{\delta t} \right\\} the approximation to the lower and upper bound of the solution can be described by \begin{align } \upsilon^{u}(k\delta t) \approx & \upsilon {0}e^{\left(\mu \frac{1}{2}\sigma^{2}\right)k\delta t+\sigma z\sqrt{k\delta t}} + \sum {k {c,i} \in \omega {k}} s {i}e^{\left(\mu \frac{1}{2}\sigma^{2}\right)(k k {c,i})\delta t+\sigma z\left(\sqrt{k\delta t} \sqrt{k {c,i}\delta t}\right)} \\\\ \upsilon^{l}(k\delta t) \approx & \upsilon {0}e^{\left(\mu \frac{1}{2}\sigma^{2}\right)k\delta t \sigma z\sqrt{k\delta t}} + \sum {k {c,i} \in \omega {k}} s {i}e^{\left(\mu \frac{1}{2}\sigma^{2}\right)(k k {c,i})\delta t \sigma z\left(\sqrt{k\delta t} \sqrt{k {c,i}\delta t}\right)} \end{align } where ฯ
l and ฯ
u are the lower and upper bound of future projection, respec tively, given by the model the set ฯk contains the step index where the effect of its corresponding cashflow is added to the solution, i e \omega {k} = \left\\{ k {c} \in \omega {c} \\,|\\, k {c} \leq k \right\\} note this section generalize the method in the previous section (without cashflow) the method in previous section is identical to the method in this section when all the cashflows are set to zero references \[1] bagchi, a (1993) optimal control of stochastic systems international series in systems and control engineering prentice hall international \[2] wiener, n (1981) collected works, edited by masani, p vol 3, mit press \[3] dirac, paul (1930) the principles of quantum mechanics (1st ed ), oxford university press \[4] shreve, s e and soner, h m (1994) optimal investment and consumption with transaction costs the annals of applied probability