﻿ Fast Fourier Transform Based Computation of American Options under Economic Recession Induced Volatility Uncertainty

Journal of Mathematical Finance
Vol.09 No.03(2019), Article ID:94538,28 pages
10.4236/jmf.2019.93026

Fast Fourier Transform Based Computation of American Options under Economic Recession Induced Volatility Uncertainty

Philip Ajibola Bankole, Olabisi O. Ugbebor

Received: May 16, 2019; Accepted: August 19, 2019; Published: August 22, 2019

ABSTRACT

The menace of Economic recession to uncertainty in the payoff of investments and standard of living cannot be over emphasized. This paper presents fast Fourier transform method for the valuation of American style options under the exposure of Economic recession. A multi-factor affine Exponential jump model with Recession induced Stochastic volatility and Intensity, which is a partial Integro-Differential Equation (PIDE) is presented. We show how to determine the characteristic function of the model via generating function. A close form characteristic formula for a financial claim satisfying the PIDE in pricing both European style and American style options in Fourier based transform was done. A numerical based Fourier transform algorithm FFT for European call option valuation was extended to the model under study. The algorithm was further extended to American call options valuation by adding premium price to the European call options price. Numerical result was presented to reflect the effect of economic recession induced volatility on options prices and that of the usual volatility. The result shows some significant vicissitudes in the options values in the two states of the Economy. The result output indicated that the model is effective and reliable compared to other existing models. The fast Fourier transform (FFT) approach gave better option value and compared to both Black-Scholes Merton (BSM) and American Option solver as shown in the table under numerical result section. We used Nigerian Flourmill Stock (NFS) prices for data calibration and reported the stock performance during the first Nigerian recession and recovery year in the Appendix section.

Keywords:

Economic Recession, Volatility Change, Characteristic Function in Affine Form, Partial Integro-Differential Equation PIDE, American Dividend Paying Options, Fourier Transform, Fast Fourier Transform (FFT)

1. Introduction

Financial mathematics encompasses many relevant areas in which option valuation is one of them and connects many different fields of study in mathematics. It relies on the application of various mathematical concepts and tools for survival. Such concepts are Probability theory, optimization theory, numerical analysis, partial differential equations, Ordinary Differential Equations, integral representations, transformation and many more. Among various transformation techniques obtainable in Mathematics such as Hankel transform, Mellin transform, Hilbert transform, Laplace transform, z-transform, and Fourier transform, this paper applies fast Fourier transform based computation of Options focusing on uncertainty in the options value induced by Economic Recession.

Instability of Economy has created a wide disparity between Economic activities and the values of Assets found in the economy. It is very obvious that Assets’ values are been affected by the state of the Economy. The effect of Economic recession to the falling of the payoff of investments and standard of living is enormous and should not be taken for granted. This is evident as the data from Nigeria Stock Exchange (NSE) indices prior to recession outbreak, during recession and after recession revealed huge lower performance of Nigeria stocks during the Nigeria economy recession of 2016. The National Bureau of Economic Research (NBER) has been recognised in the United States as an official body saddled with the responsibility of providing an accurate information on economic recession dates and business cycle publishing since 1929. The definition of recession according to NBER cited in [1] is a significant decline in economic activity spread across the economy, lasting more than a few months, normally visible in production, employment, real income, and other indicators. Some other definitions of economic recession obtainable in financial press emphasised that recession begins with two consecutive quarters of decline in Gross Domestic Product (GDP). NBER ponders on GDP as the single best measure of total economic activity and reflects on the GDP definition to be too narrow in measuring economic activity and to reliably date economic recessions. There exists scenario whereby recession may not include two consecutive quarters of negative growth such as United States recession of 2001 but according to Thomas Hsu [1] , declines in GDP are closely correlated to recession periods. Nigerian economic recession outbreak in year 2016 was based on decline in GDP and rise in some other macroeconomic indicators such as high inflation rate and unemployment rate after two consecutive quarters. Some measures were invented by axiomatic methods: for example, as cited in [2] . Probability Measure was invented by A.N. Kolmogoroff in 1933, Possibility Measure by L.A. Zadeh in 1978 and Uncertain Measure (B. Liu) in 2007 (see [2] [3] [4] [5] [6] ). Probability theory has been seen as the vehicle for dealing with uncertainty in finance and insurance risks. Probability theory as a mathematical theory is useful in describing and analysing situations where randomness or uncertainty are present [7] . Definition of uncertainty had been given by different authors in various scenarios. According to [8] (Dungan et al., 2002), Uncertainty is a multi-faceted characterization about data or predictions made from data that may include several concepts including error, accuracy, validity, quality, noise and confidence and reliability. Uncertainty connotes difficulty in predicting the outcome of an event due to inexact knowledge of information concerning the outcome of the event in a particular time. Uncertainty appears to endogenously increase during recessions, as lower economic growth induces greater micro and macro uncertainty [9] . It was stressed in [10] that investment growth experiences downward trend during the period of recession.

The following authors [11] - [21] have considered uncertainty with respect to models. The aspect of uncertainty in terms of economic recession was not discussed in their formulations. It was recognized that prices of insurance and financial risks are not solely determined by mathematics or economic theories but are also impacted by activities taking place in financial markets. During the period of economic crisis, prices of assets tend to change stochastically in the financial market (especially while dealing with stocks) since there is much tendency for fluctuations of asset prices due to fluctuation in the underlying term structures typically, the volatility variation. The study of volatility of stock prices is very essential due to its effect on the prices of stocks in the financial market. Pfante, O. and Bertschinger, N. (2018), [21] consider the uncertainty of volatility estimates from Heston Greeks. It was stressed in [22] that volatility has been the cause for several Statistical properties of observed stock prices processes. It was further emphasized in [22] that volatility clustering is commonly accompanied by other large fluctuations and similar for small changes. As volatility is not observed it has to be estimated from market prices, i.e., as the implied volatility from option prices. There is strong tendency that predictions of the volatility of asset prices induced by economic recession factor among other market risk factors may not really be accurate. This becomes a challenge for investors that do not like taking risk, that is, the risk averters. There are some existing methods of calculating the volatility of asset prices in financial market especially while dealing with instantaneous prices (i.e. the real time prices) of assets in the financial market. Among the methods known to us are 1) Standard deviation approach, 2) using Historical value of the volatility and 3) implied volatility. The Standard deviation approach tells us how tightly the stock price is grouped around the mean or moving average. When the prices are spread apart, the implication is that one has a relatively large standard deviation but if the prices are relatively closed to each other or bunched together, this connotes that the standard deviation is negligible.

The causal effect of volatility on stock prices is discussed briefly as follows. The stock market prices rise when volatility decreases and increase in volatility causes fall in stock market value (prices). Increase in volatility leads to increase in market risk but decrease in returns of the market. In the case of options on stocks market, the causal effect of volatility changes on options depends on the type. For example, increase in volatility leads to increase in call options value but for put options, increase in volatility leads to decline in the payoff of the put options vice versa. An intuition behind the introduction of recession induced volatility uncertainty is revealed by huge volatility fluctuations during the period of Economic recession compared to the period of normalcy (recession-free). This in turn affects investors’ prediction in the market. Since Economic recession induces a high level of uncertainty on investors activities to include decision making and the payoffs of stocks in general, we then proposes giving a close attention to volatility changes in relation to economic recession and financial models for options valuation. Nigeria economic recession outbreak in 2016 and its effects on the payoffs uncertainty of Nigeria Stocks Exchange (NSE) among other investments is among the motivating factors for proposing economic recession induced volatility formulation in Options pricing. A good knowledge of the behavior or the level of uncertainty in economic recession induced volatility by investor’s will help in decision making during recession period.

The rest part of the paper is organized as follows: Preliminaries on financial modeling is discussed in Section 2 in addition to some other subsections to include Fourier transform, uncertainty and uncertain measure. Section 3 deals with accounting for jumps in the asset price linked with recession, Section 4 shows the model formulation. Numerical Fourier based transform of options is presented in Section 5 while Section 6 is conclusion.

2. Preliminaries on Financial Modeling

2.1. Some Literature Review on Uncertainty

The theory of uncertainty in financial market could be traced back to the research work of [2] [3] [4] . Without uncertainty, the probabilities of risky events are known and frictionless markets can precisely price contracts contingent on risky events broadly. The volatility of the stock market or GDP is often used as a measure of uncertainty because when a data series becomes more volatile it is harder to forecast. Other common measures of uncertainty include forecaster disagreement, mentions of “uncertainty” in news, and the dispersion of productivity shocks to firms [9] . It was further stressed by Nicholas Bloom [9] that the volatility of stock markets, bond markets, exchange rates, and GDP growth all rise suddenly during economy recessions. In Philip et al. [10] , it was stressed that investment under the exposure of economic recession tends to have negative growth. This in turn poses a lot of challenges on investors in a financial or money market as investors’ decision in the market also depends on the daily information on the state of the economy.

In the history of options pricing in financial market, the famous Black - Scholes model [23] with the assumption of constant rate of return and volatility have been criticized by many researchers due to the fact that it does not reflect the stochastic nature of financial markets wholly. As a result of this deficiency in the B-S model, some other realistic models have been formulated and come to stay in the sense that the models shows better random movement of financial market. Example of such models are: Double exponential jump model, Regime switching model, Stochastic volatility models of Heston [24] , etc. In [25] , it was asserted that regime-switching behaviour captures changing preferences and beliefs of investors concerning asset prices as the state of financial market changes. Regime-Switching formulation in valuation of financial market was known to be introduced by J. Hamilton [26] . Since then, there have been some development on considering his suggestion. A closed form formula for the valuation of European call option in a two-state economy under the assumption of no-arbitrage was presented in [27] by Guo, X. European and American-style derivatives pricing was presented in [28] by Bollen using lattice-based approach.

Fourier transform approach of pricing options especially European-style options have been looked into by various authors. The one with great popularity known to us among others was that of Carr and Madan (1999) [29] in pricing European call option. Some other authors have also shown other mathematical formulation of Fourier transform of European options such as in [30] . As far as we know, the issue of economy recession factor in option pricing has not been incorporated into fast Fourier Transform of both European and American-style options. Very few authors have worked on Fourier Transform of American-style options compare to European style options. The notable one known to us as at the time of this research was the work of Oleksandr, Z. in (2010) [31] . Geske-Johnson scheme with Richardson Extrapolation was adopted to numerically extend fast Fourier Transform algorithm to American-style security, which features a continum of potential exercise times up to expiration.

Our contribution to the existing literature on options valuation is highlighted in sequel:

1) We incorporate economic recession induces volatility uncertainty into exponential jump model with stochastic volatility and intensity. The economic recession induce volatility is assumed to be an uncertain variable.

2) Derivation of Characteristic function of the affine model with recession uncertainty effect is presented.

3) The Fourier transform of the affine model is performed.

4) The Characteristic function of the affine model is extended to fast Fourier transform algorithm of Carr & Madam [29] in pricing European options. It was extended to pricing American style options by adding time premium such that the limit value tends to zero at the expiry time.

5) Volatility Surface of the affine model based on the recession induced uncertainty is presented.

2.2. Uncertainty

Definition 1: [5] [6] . An uncertain variable $\xi$ is a measurable function from an uncertainty space to the set of real numbers.

Definition 2: O.O. Ugbebor [6] . Let $\Gamma$ be a non-empty set. A collection $\mathcal{L}$ consisting of subsets of $\Gamma$ is called an algebra over $\Gamma$ if the following three conditions hold:

$\Gamma \in \mathcal{L}$ ;

• if $\wedge \in \mathcal{L}$, then ${\wedge }^{c}\in \mathcal{L}$ ; and

• if ${\wedge }_{1},{\wedge }_{2},\cdots ,{\wedge }_{n}\in \mathcal{L}$, then

$\underset{i=1}{\overset{n}{\cup }}{\wedge }_{i}\in \mathcal{L}$

Under the notion of countable union, if ${\wedge }_{1},{\wedge }_{2},\cdots ,{\wedge }_{n}\in \mathcal{L}$ such that

$\underset{i=1}{\overset{n}{\cup }}{\wedge }_{i}\in \mathcal{L}$

the collection $\mathcal{L}$ becomes a σ- algebra over $\Gamma$.

Example 1: Suppose $\wedge$ is a proper subset of $\Gamma$. Then the set $\left\{\varnothing ,\wedge ,{\wedge }^{c},\Gamma \right\}$ is a σ-algebra over $\Gamma$.

Remark 1: The collection $\Gamma$ is closed under countable union, countable intersection, difference and limit.

Definition 3: [5] [6] Let $\Gamma \ne \varnothing$ be a non-empty set and suppose $\mathcal{L}$ is a σ-algebra over $\Gamma$. Then $\left(\text{Γ},\mathcal{L}\right)$ is called a measurable space and any member in $\mathcal{L}$ is a measurable set.

Definition 4: [6] A function $\xi :=\left(\text{Γ},\mathcal{L}\right)\to R$ is said to be measurable if

${\xi }^{-1}\left(\mathcal{B}\right)=\left\{\gamma \in \Gamma |\xi \left(\gamma \right)\in \mathcal{B}\right\}\in \mathcal{L}$ (1)

where $\left(\text{Γ},\mathcal{L}\right)$ is a measurable space and $\mathcal{B}\left(R\right)$ is a Borel set of real numbers.

2.3. Uncertain Measure

B. Liu, [4] [5] highlighted the following three axioms for an uncertain process:

Axiom 1: (Normality Axiom); $\mathcal{M}\left\{\text{Γ}\right\}=1$ for the universal set $\Gamma$.

Axiom 2: (Duality Axiom): $\mathcal{M}\left(\wedge \right)+\mathcal{M}\left({\wedge }^{c}\right)=1$ for any even $\wedge$.

Axiom 3: (Subadditivity Axiom); $\mathcal{M}\left\{{\cup }_{i=1}^{\infty }{\wedge }_{i}\right\}\le {\sum }_{i=1}^{\infty }\mathcal{M}\left\{{\wedge }_{i}\right\}$ for every countable sequence of events ${\wedge }_{1},{\wedge }_{2},\cdots$.

Remark 2:

For some application of B. Liu (2007)’s uncertainty theory, the reader should see the research work of the following authors in [12] [13] [14] as cited earlier.

2.4. Fourier Transform

Fourier transform in mathematical finance was firstly used to determine the distribution of an underlying asset price under the stochastic volatility model by the inversion method, Stein et al. (1991) [18] . The following definitions and conditions are necessary for the understanding of Fourier transform technique.

Definition 5: Let $f\left(x\right)$ be a Lebesgue-measurable function of $x\in R$, then the ${L}^{2}$ -norm of f is defined by

$‖f‖={\left({\int }_{-\infty }^{\infty }{|f\left(x\right)|}^{2}\text{d}x\right)}^{\frac{1}{2}}$ (2)

where

${L}^{2}=\left\{f:‖f‖<\infty \right\}.$ (3)

Equivalently, $f\left(x\right)$ is a piecewise integrable real function over the entire real line satisfying the condition

${\int }_{-\infty }^{\infty }|f\left(x\right)|\text{d}x<\infty .$ (4)

For the Fourier transform and inverse Fourier transform of a function to exist, it is necessary that the function is absolutely integrable over R with finite value as defined in (3) above. An absolutely integrable functions on a given interval $\left(a,b\right)$ are said to be on the space of ${L}^{1}\left(a,b\right)$. The notion of Fourier transform can be extended to the square integrable functions.

Definition 6: A square integrable functions is defined by

${\int }_{-\infty }^{\infty }{|f\left(x\right)|}^{2}\text{d}x<\infty .$ (5)

The space of square integrable functions is represented by ${L}^{2}\left(R\right)$ over the real line or simply an interval $\left(a,b\right)$.

Definition 7: Let $f:R\to R$ be a real-value function for $x\in R$, then the Fourier transform of f is defined as

$\mathcal{F}\left(f\left(x\right);w\right)=\stackrel{^}{f}\left(w\right)={\int }_{-\infty }^{\infty }f\left(x\right){\text{e}}^{-i\omega x}\text{d}x,\text{\hspace{0.17em}}\omega \in R,$ (6)

where $i=\sqrt{-1}$ and $\omega$ is a parameter.

We can recover $f\left(x\right)$ from $f\left(\omega \right)$ by the inverse Fourier transform.

Definition 8: The inverse Fourier transform of the real valued function $f\left(x\right)$ is defined by

${\mathcal{F}}^{-1}\left(\stackrel{^}{f}\left(\omega \right);x\right)=f\left(x\right)=\frac{1}{2\text{π}}{\int }_{-\infty }^{\infty }\stackrel{^}{f}\left(\omega \right){\text{e}}^{i\omega x}\text{d}\omega ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\in R$ (7)

which belongs to either ${L}^{1}$ or ${L}^{2}$ - spaces.

3. Accounting for Jumps in the Asset Prices Linked with Recession

The major market parameters which most financial modeler will wish to consider while formulating a financial model relies on price, interest rate, dividend rate, volatility and time. In an economy threatened by recession, asset price tends to experience lots of fluctuation. Asset price especially stocks can never be stable. It is very necessary to consider interest rate and volatility to be stochastic in nature. Even if one considers interest rate to be constant, in reality volatility cannot be constant. The Stock volatility is seen as a measure of the uncertainty on the payoff or returns of the stock which investors look up to for decision making and taking. During economic recession, the rate of volatility variation is considered higher than the state of normalcy of an economy. As a result an economy recession factor is invoked in this paper as we shall see later.

Jumps in the Underlying Stock Price

Consider the dynamics of a stock price $S\left(t\right)$ given by

$\frac{\text{d}S\left(t\right)}{S\left(t\right)}=\left(r-q-\lambda \left(t\right)m\right)\text{d}t+\sqrt{\sigma \left(t\right)}\text{ }\text{d}{W}_{s}\left(t\right)+\left({\text{e}}^{\nu }-1\right)\text{d}N\left(t\right)$ (8)

where r is the interest rate, q is the dividend rate, $N\left(t\right)$ is a Poisson process with stochastic intensity $\lambda \left(t\right)$, m is an average jump amplitude given as $m={E}^{\stackrel{^}{P}}\left({\text{e}}^{\nu }-1\right)$ where jump size $\nu$ is a random variable. The equation consists of diffusion process and jump process. There are two sources of fluctuations in the stock price which we classified as “changes in economic state” and “changes due to supply and demand factors”. The inclusion of jump in the model is to cater for the arrival of useful information into the market that will have an abnormal consequence on the stock price. Among the factors that are responsible for jumps in the stock price highlighted by Matthew, S.M in [32] , we added economic recession factor. The factors are grouped into:

• Firm specific jumps: Caused by news inflow to the market on individual firm’s profit/loss report.

• Industry/sector specific jumps: Caused by news that can affect specific company or industry such as news on sudden declaration of Holiday.

• Market specific jumps: This jumps are caused by news such as oil prices, interest rates, credit spreads etc. that affect the market.

• Economic recession jumps: This is caused by the news inflow to the market on the state of the economy to include deep economy crisis with recession probability indicators.

In our model formulation of stock price to capture economy recession induced volatility uncertainty on stock in the next section, we further grouped the above highlighted factors “(a)-(c)” as other sources with volatility ${\sigma }_{t}^{*}$ while the factor (d) Economic Recession induced volatility source as ${\sigma }_{t}^{rec}$.

Poisson process $N\left(t\right)$ has probability density function given as

$P\left({N}_{t}=\nu \right)=\frac{\left(\lambda {t}^{\nu }\right)}{\nu !}{\text{e}}^{-\lambda t},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\nu =0,1,2,\cdots$ (9)

where $\lambda$ is the intensity. In line with Matthew [24] , if information on economic recession or other source of panic enters the financial market (stock) causing an instantaneous jump in the stock price such that the price jump change from $S\left(t\right)\to {\nu }_{t}S\left(t\right)$ where ${\nu }_{t}$ is taking as the absolute magnitude jump. Then the relative change in price is given by

$\frac{\text{d}S\left(t\right)}{S\left(t\right)}=\frac{{\nu }_{t}S\left(t\right)-S\left(t\right)}{S\left(t\right)}=\frac{S\left(t\right)\left({\nu }_{t}-1\right)}{S\left(t\right)}={\nu }_{t}-1.$ (10)

4. The Model Formulation

Let $\xi \left(\gamma \right)$ be an economy recession induced parameter variation define on an uncertain space $\left(\text{Γ},\mathcal{L},\mathcal{M}\right)$, with $\gamma =\left\{{\gamma }_{1},{\gamma }_{2}\right\}$ such that

$\mathcal{M}\left\{{\gamma }_{1}\right\}=\alpha ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathcal{M}\left\{{\gamma }_{2}\right\}=\beta .$

Then

$\xi \left(\gamma \right)=\left\{\begin{array}{ll}0,\hfill & \text{if}\text{\hspace{0.17em}}\gamma ={\gamma }_{1}\hfill \\ 1,\hfill & \text{if}\text{\hspace{0.17em}}\gamma ={\gamma }_{2}\hfill \end{array}$ (11)

is an uncertain variable. We have that

$\mathcal{M}\left\{\xi =0\right\}=\mathcal{M}\left\{\gamma |\xi \left(\gamma \right)=0\right\}=\mathcal{M}\left\{{\gamma }_{1}\right\}=\alpha$

and

$\mathcal{M}\left\{\xi =1\right\}=\mathcal{M}\left\{\gamma |\xi \left(\gamma \right)=1\right\}=\mathcal{M}\left\{{\gamma }_{2}\right\}=\beta$

such that $\mathcal{M}\left\{{\gamma }_{1}\right\}+\mathcal{M}\left\{{\gamma }_{2}\right\}\equiv \alpha +\beta =1$.

Example 2: Let the above economic recession induce parameter variation $\xi \left(\gamma \right)$ defined on an uncertain space $\left(\text{Γ},\mathcal{L},\mathcal{M}\right)$, with $\gamma =\left\{{\gamma }_{1},{\gamma }_{2}\right\}$ be defined such that

$\mathcal{M}\left\{{\gamma }_{1}\right\}=0.32,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathcal{M}\left\{{\gamma }_{2}\right\}=0.68$

satisfying Equation (11), then

$\mathcal{M}\left\{\xi =0\right\}=\mathcal{M}\left\{\gamma |\xi \left(\gamma \right)=0\right\}=\mathcal{M}\left\{{\gamma }_{1}\right\}=0.32$

and

$\mathcal{M}\left\{\xi =1\right\}=\mathcal{M}\left\{\gamma |\xi \left(\gamma \right)=1\right\}=\mathcal{M}\left\{{\gamma }_{2}\right\}=0.68$

Let $\left(\text{Ω},\mathcal{F},\left\{{\mathcal{F}}_{t}^{mket}\right\}\subseteq \left\{{\mathcal{F}}_{t}^{eco}\right\},\stackrel{^}{\mathcal{P}}\right)$ be a filtered probability space where ${\mathcal{F}}_{t}^{mket}$ and ${\mathcal{F}}_{t}^{eco}$ describes the filtration in the market and the economy where market activities takes place, respectively. The ${\mathcal{F}}_{t}^{mket}$ and ${\mathcal{F}}_{t}^{eco}$ are the market information and economy information available up to a time $t\in \left(0,T\right]$ respectively.

Definition 9: [33] An uncertain random variable is a measurable function $\xi \in {R}^{p}$ (resp. ${R}^{p×m}$ ) from an uncertainty probability space $\left(\Gamma ×\Omega ,\mathcal{L}\otimes \mathcal{F},\mathcal{M}×P\right)$ to the set in ${R}^{p}$ (resp. ${R}^{p×m}$ ), for any Borel set $A\in {R}^{p}$ (resp, ${R}^{p×m}$ ) the given set $\left\{\xi \in A\right\}=\left\{\left(\gamma ,\omega \right)\in \Gamma ×\Omega :\xi \in \left(\gamma ,\omega \right)\in A\right\}\in \mathcal{L}\otimes \mathcal{F}$.

The expected value of the uncertain random variable $\xi$ is defined as

$\begin{array}{c}E\left[\xi \right]={E}_{p}\left[{E}_{\mathcal{M}}\left[\xi \right]\right]\\ ={\int }_{\Omega }\left[{\int }_{0}^{\infty }\mathcal{M}\left\{\xi \ge r\right\}\text{d}r\right]P\left(\text{d}\omega \right)-{\int }_{\Omega }\left[{\int }_{0}^{\infty }\mathcal{M}\left\{\xi \le r\right\}\text{d}r\right]P\left(\text{d}\omega \right)\end{array}$ (12)

where ${E}_{p}$ and ${E}_{\mathcal{M}}$ represent the expected values under the uncertainty space and the probability space, respectively.

The above definition connects both the notion of probability and uncertainty such that the random variable is defined from a probability space to uncertainty space. It makes sense to introduce the notion of uncertainty to financial models especially while pricing during economy recession or strong financial crisis in the market. The works of the authors cited in this paper to have contributed to parameter uncertainty in financial models did not consider uncertainty with respect to economic recession. We extended the notion of uncertainty to the term structure of stochastic volatility to include economy recession.

Let an asset X be described as a two-state regime switching process which is free to jump between the states defined as

${X}_{t}=\left\{\begin{array}{ll}1,\hfill & \text{if the economy is in a state of Expansion};\hfill \\ 2,\hfill & \text{if the economy is in a state of Recession}\text{.}\hfill \end{array}$ (13)

Suppose further that there is room for transition between the two states which evolves as a Poisson process

$Prob\left({t}_{lm}^{\text{*}}>t\right)=\mathrm{exp}\left(-{\lambda }_{lm}t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}l,m=1,2$ (14)

where ${\lambda }_{lm}$ is the rate of transition from state l to state m and the times spent in state l before transiting to state m is ${t}_{lm}^{\text{*}}$.

Assumption 1: Assume there are two sources of volatility on the asset, ${X}_{t}$, i.e volatility ${\sigma }_{{X}_{t}}^{rec}$ and ${\sigma }_{{X}_{s}}^{othersources}$ representing market volatility triggered by economy recession and other sources (usual) respectively such that

${\sigma }_{{X}_{s}}^{market}={\sigma }_{t}^{rec}+{\sigma }_{s}^{usual}$, but ${\sigma }_{t}^{rec}\ne {\sigma }_{s}^{usual}$ (15)

Assumption 2: It is assumed that an economy cannot be in the two states at the same time. Therefore, when pricing in the financial (stock) market, the choice of the volatility parameter depends on the state of the economy. For example, if the stock is under the exposure of economy recession, then market volatility of the stock ${\sigma }^{market}={\sigma }_{t}^{rec}+{\sigma }_{s}^{othersources}$ since recession will have an influence on the market but if there is no recession, then market volatility ${\sigma }_{t}^{market}={\sigma }_{s}^{*}$ usual volatility, here $t=s$, since ${\sigma }_{s}^{rec}=0$.

Suppose further that the asset X price is defined on a filtered probability space $\left(\Omega ,{\mathcal{F}}_{t}^{mket},{\mathcal{F}}_{t}^{rec},\stackrel{^}{\mathcal{P}}\right)$ such that the market filtration is generated by the combination of Wiener process and jump process at a given time, $t\in \left[0,T\right]$ and $\stackrel{^}{P}$ is taken as a risk-neutral probability measure. The dynamics of the underlying stock price $S\left(t\right)$ (not necessarily a stock) is given by

$\left\{\begin{array}{l}\text{d}S\left(t\right)=\left(r-q-\lambda \left(t\right)m\right)S\left(t\right)\text{d}t+\sqrt{\sigma \left(t\right)}\text{ }\text{d}{W}_{s}\left(t\right)+\left({\text{e}}^{\nu }-1\right)S\left(t\right)\text{d}N\left(t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}S\left(0\right)={S}_{0}>0\\ \text{d}\sigma \left(t\right)={\kappa }_{\sigma }\left({\beta }^{*}+{\beta }^{rec}-\sigma \left(t\right)\right)\text{d}t+{\xi }_{\sigma }\sqrt{\sigma \left(t\right)}\text{ }\text{d}{W}_{\sigma }\left(t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sigma \left(0\right)={\sigma }_{0}>0\\ \text{d}\lambda \left(t\right)={\kappa }_{\lambda }\left(\theta -\lambda \left(t\right)\right)\text{d}t+{\xi }_{\lambda }\sqrt{\lambda \left(t\right)}\text{ }\text{d}{W}_{\lambda }\left(t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\lambda \left(0\right)={\lambda }_{0}>0\end{array}$ (16)

where r is the interest rate, q is the dividend rate, $N\left(t\right)$ is a Poisson process with stochastic intensity $\lambda \left(t\right)$, m is an average jump amplitude given as $m={E}^{\stackrel{^}{P}}\left({\text{e}}^{\nu }-1\right)$ where jump size $\nu$ is a random variable. The parameters ${\xi }_{\sigma }>0$ and ${\xi }_{\lambda }>0$ are constants with mean - reverting rates ${\kappa }_{\sigma }$ and ${\kappa }_{\lambda }$ which are also taken to be positive constants. The stock market stochastic volatility

$\sigma \left(t\right)=\left\{\begin{array}{ll}{\sigma }^{rec}+{\sigma }^{*},\hfill & \text{if the economy is in recession};\hfill \\ {\sigma }^{*},\hfill & \text{if economy is recession free}\text{.}\hfill \end{array}$ (17)

where ${\sigma }^{rec}$ is economy recession induced volatility in the stock market due to the information flow from the Recessed Economy into the stock market while ${\sigma }^{*}$ is the volatility in the market from other sources. Furthermore, the constants ${\beta }^{\text{*}}$ and ${\beta }^{rec}$ are respectively regarded as the usual long-term volatility, long-term volatility influenced by recession during the period of economy recession and $\theta$ is the intensity constant. The ${\beta }^{rec}$ as specified above is to account for the possibility of economy recession source of long-term volatility in addition to the usual long-term volatility that may arise during the life-span of the quoted option on the stock. Also, ${W}_{s}\left(t\right)$ and ${W}_{\sigma }\left(t\right)$ are Wiener processes but ${W}_{\lambda }\left(t\right)$ evolve independently of the correlated Wiener process ${W}_{s}\left(t\right)$ and ${W}_{\sigma }\left(t\right)$ in the above model. The correlation of ${W}_{s}\left(t\right)$ and ${W}_{\sigma }\left(t\right)$ is such that $〈{W}_{s}\left(t\right),{W}_{\sigma }\left(t\right)〉=\rho \text{d}t$ with $|\rho |$ is less than 1. As we take m to be an average jump amplitude since the stock price is expected to jump either upward or downward. Then setting

$m=\frac{p}{1-{\xi }_{u}}+\frac{q}{1+{\xi }_{d}}-1\equiv {\int }_{-\infty }^{\infty }{\text{e}}^{\nu }f\left(\nu \right)\text{d}\nu -1$ (18)

where p is the probability of upward move jump and q, the downward move jump such that $p+q=1$. In [16] , the mean positive jump and negative jumps

were given as $\frac{1}{{\xi }_{u}}$ and $\frac{1}{{\xi }_{d}}$, respectively. An analysis on double exponential

jump model with stochastic volatility and stochastic intensity was done by Jiexiang Huang et al in [34] and solved. The major distinction between our formulation here and that of [34] is the incorporation of economic recession induced uncertainty into the term structure of the volatility of the model. Consequently, we propose the model in Equation (16) to be referred to as an “Uncertain Affine Exponential jump model with Recession induced Stochastic volatility and stochastic Intensity”.

Remarks 3:

The stochastic volatility $\sigma {\left(t\right)}_{{x}_{t}}$ of the stock has been formulated to capture recession induced volatility.

4.1. Determination of Characteristics Function

The notion of characteristics function is indispensable in the study of random variables in the context of jump diffusion processes. It is a viable tool in Fourier transformation of options pricing. It is worth noting that the distribution function of jump diffusion processes in closed form may not be readily available or known but the characteristics function is explicitly known. However, the characteristics function of some certain stochastic processes or model may not be readily available in closed form but it can be determined.

Definition 10: The characteristics function of a random variable X is a function

${\phi }_{X}\left(\omega \right)=E\left({\text{e}}^{i\omega X}\right)={\int }_{-\infty }^{\infty }{f}_{X}\left(x\right){\text{e}}^{i\omega x}\text{d}x,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\omega \in R,$ (19)

where ${f}_{X}\left(x\right)$ is the probability density function of the random variable X.

This definition almost coincide with the definition of inverse Fourier Transform in (6) with the exception of $\frac{1}{\text{2π}}$. For a Poisson process say ${N}_{t}$ with

parameter $\lambda$, the characteristics function is given by

${\phi }_{N}\left(\omega \right)\equiv E\left({\text{e}}^{i\omega {N}_{t}}\right)=\mathrm{exp}\left\{\lambda t\left({\text{e}}^{i\omega }-1\right)\right\}$ (20)

There is relationship between Moment-generating function and characteristics function. Moment generating function say, M of a real valued random variable is very useful in the study of probability theory as an alternative specification of the probability distribution of the random variable $X\left(t\right)$.

Definition 11: Moment-generating function of a random variable, X, is defined as

${M}_{x}\left(t\right)=E\left({\text{e}}^{tX}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}t\in R,$ (21)

if the expectation exists.

Example 3: For a Poisson $Pois\left(\lambda \right)$ distribution, the moment-generating function ${M}_{x}\left(t\right)={\text{e}}^{\lambda \left({\text{e}}^{t}-1\right)}$ and its characteristic function $\phi \left(t\right)={\text{e}}^{\lambda \left({\text{e}}^{it}-1\right)}$.

For a continuous probability density function, the moment-generating function is

${M}_{x}\left(t\right)={\int }_{-\infty }^{\infty }{\text{e}}^{tx}f\left(x\right)\text{d}x$ (22)

where $f\left(x\right)$ is the probability distribution function.

To this end, let $x\left(\tau \right):=\mathrm{ln}S\left(\tau \right)$ with $\tau =T-t$ or $\tau ={t}^{\text{*}}-t$ as the case may be where $x\left(\tau \right)$ is the log stock price in the model (16) and $\tau$ is the time difference between the maturity time and the present time or the time difference between the optimal time and the present time for exercising the option as the case may be. The latter is for Early exercise of American option peradventure one is able to determine the optimal exercise time for the option. In order to derive the characteristics function for the proposed model in 16, we need to consider Feynman-Kac Formula in what follows.

4.2. The Feynman-Kac Formula Revisited

The formula states that a probabilistic expectation wrt some Ito-diffusion processes can be obtained as a solution of a related partial differential equation. For an example, for a 1-dimensional stochastic process $X\left(t\right)=\left({{X}^{\prime }}_{t},\cdots ,{X}_{t}^{d}\right)$ which is a solution of stochastic differential equation

$\text{d}{X}_{t}^{i}={\mu }_{i}\left(t,{X}_{t}\right)\text{d}t+{\sigma }_{i}\left(t,{X}_{t}\right)\text{d}{W}_{t}^{i}$

where ${W}_{t}^{i},i=1,\cdots ,d$ are Wieners process with correlation $〈\text{d}{W}_{t}^{i},\text{d}{W}_{t}^{j}〉={\rho }_{ij}\text{d}t$. Let $\mathcal{Z}\left(x\right)=\mathcal{Z}\left({x}_{1},\cdots ,{x}_{d}\right)$ be some payoff.

Then the function

$g\left(t,x\right):=E\left[\mathcal{Z}\left({X}_{T}\right)|{X}_{t}=x\right]$ (23)

is a solution of the PDE

$\frac{\partial g}{\partial t}+{\sum }_{i=1}^{d}\mu \left(t,x\right)\frac{\partial g}{\partial {x}_{i}}+\frac{1}{2}{\sum }_{i,j=1}^{d}{\rho }_{ij}{\sigma }_{i}\left(t,x\right){\sigma }_{j}\left(t,x\right)\frac{{\partial }^{2}g}{\partial {x}_{i}\partial {x}_{j}}=0$ (24)

with final condition $g\left(t,x\right)=\mathcal{Z}\left(x\right)$. For 1-dimensional stochastic process Xt which is a solution of the SDE

$\text{d}{X}_{t}=\mu \left(t,{X}_{t}\right)\text{d}{t}_{\sigma }\left(t,{X}_{t}\right)\text{d}{W}_{t}.$ (25)

By Feynman-Kac formula, for any function say

$g\left(t,{X}_{t}\right)=E\left[g\left(t,{X}_{t}\right)|X\left(t\right)\right]$ (26)

is the solution to the PDE

$\frac{\partial g}{\partial t}+\mu \left(t,{X}_{t}\right)\frac{\partial g}{\partial x}+\frac{1}{2}{\sigma }^{2}\left(t,{X}_{t}\right)\frac{{\partial }^{2}g}{\partial {x}^{2}}=0.$ (27)

The Equation (26) is useful as it can be extended to characteristics function as

$g\left(\tau ,{X}_{\tau }\right)=\phi \left(\omega ,{X}_{\tau },\tau \right)={\text{e}}^{i\omega {X}_{\tau }}.$ (28)

Similarly,

$g\left(\tau ,{X}_{t}\right)=\phi \left(\omega ,{X}_{t},t\right)={\text{e}}^{i\omega {X}_{t}}=E\left[{\text{e}}^{i\omega {X}_{T}}|{X}_{t}\right].$ (29)

Remarks 4:

1) $\phi \left(\omega ,{X}_{\tau },\tau \right)$ is the characteristics function of ${X}_{\tau }$.

2) The moment-generating function ${M}_{x}\left(\omega \right)={\phi }_{x}\left(-i\omega \right)$.

4.3. Solution to the Proposed Uncertain Affine Exponential Jump Model with Recession Induced Stochastic Volatility and Intensity

Consider the Equation (16) rewritten in the form

$\left\{\begin{array}{l}\text{d}\mathrm{ln}S\left(t\right)=\left(r-q-\lambda \left(t\right)m\right)\text{d}t+\sqrt{\sigma \left(t\right)}\text{ }\text{d}{W}_{s}\left(t\right)+\left({\text{e}}^{\nu }-1\right)\text{d}N\left(t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}S\left(0\right)={S}_{0}>0\\ \text{d}\sigma \left(t\right)={\kappa }_{\sigma }\left({\beta }^{*}+{\beta }^{rec}-\sigma \left(t\right)\right)\text{d}t+{\xi }_{\sigma }\sqrt{\sigma \left(t\right)}\text{ }\text{d}{W}_{\sigma }\left(t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sigma \left(0\right)={\sigma }_{0}>0\\ \text{d}\lambda \left(t\right)={\kappa }_{\lambda }\left(\theta -\lambda \left(t\right)\right)\text{d}t+{\xi }_{\lambda }\sqrt{\lambda \left(t\right)}\text{ }\text{d}{W}_{\lambda }\left(t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\lambda \left(0\right)={\lambda }_{0}>0\end{array}$ (30)

where $〈\text{d}{W}_{s}\left(t\right),\text{d}{W}_{\sigma }〉=\rho \text{d}t$ and $〈\text{d}{W}_{s}\left(t\right),\text{d}{W}_{\lambda }\left(t\right)〉=0=〈\text{d}{W}_{\sigma }\left(t\right),\text{d}{W}_{\lambda }\left(t\right)〉$ and the volatility is taken as $\sigma \left(t\right)={\sigma }^{*}+{\sigma }^{rec}$ as defined in Equation (17). The PIDE for the moment-generating function $M\left(\omega \right)$ of ${x}_{\tau }:=\text{ln}{S}_{\tau }$ for the model above is given by

$\begin{array}{l}{M}_{\tau }+\left(r-q-\frac{1}{2}\sigma \left(\tau \right)-\lambda \left(\tau \right)m\right){M}_{x}+\frac{1}{2}\sigma {M}_{xx}\\ \text{ }+{\kappa }_{\sigma }\left({\beta }^{\text{*}}+{\beta }^{rec}-\sigma \left(\tau \right)\right){M}_{\sigma }+\frac{1}{2}{\xi }_{\sigma }^{2}\sigma \left(\tau \right){M}_{\sigma \sigma }\\ \text{ }+\rho \xi \sigma {M}_{x\sigma }+{\kappa }_{\lambda }\left(\theta -\lambda \left(t\right)\right){M}_{\lambda }+\frac{1}{2}{\xi }_{\lambda }^{2}\lambda {M}_{\lambda \lambda }\\ \text{ }+\lambda {\int }_{-\infty }^{\infty }\left(M\left(x+\nu \right)-M\left(\nu \right)\right)f\left(\nu \right)\text{d}\nu =0\end{array}$ (31)

where ${M}_{\tau }={\text{e}}^{{x}_{\tau }\omega }$ is the final condition and $m={\int }_{-\infty }^{\infty }{\text{e}}^{\nu }f\left(\nu \right)\text{d}\nu -1$.

The solution to the PIDE (30) is speculated to be of the form

$M\left(\omega ,{x}_{\tau },{\sigma }_{\tau },{\lambda }_{\tau },\tau \right)={\text{e}}^{\left(r-q\right)\tau \omega +\omega {x}_{\tau }+A\left(\omega ,\tau \right)+B\left(\omega ,\tau \right){\sigma }_{\tau }+C\left(\omega ,\tau \right){\lambda }_{\tau }}.$ (32)

Substituting Equation (32) into (31) gives

$\begin{array}{l}-\left(\frac{\partial A}{\partial \tau }\left(\omega ,\tau \right)+\frac{\partial B}{\partial \tau }\left(\omega ,\tau \right){\sigma }_{\tau }+\frac{\partial C}{\partial \tau }\left(\omega ,\tau \right){\lambda }_{\tau }\right)\\ =\frac{1}{2}{\sigma }_{\tau }\left({\omega }^{2}-\omega \right)-{\lambda }_{\tau }\left(m\omega -{\int }_{-\infty }^{\infty }{\text{e}}^{\nu \omega }f\left(\nu \right)\text{d}\nu -1\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{ }+{\kappa }_{\sigma }\left({\beta }^{\text{*}}+{\beta }^{rec}-{\sigma }_{\tau }\right)B\left(\omega ,\tau \right)+\frac{1}{2}{\xi }_{\sigma }^{2}{\sigma }_{\tau }{B}^{2}\left(\omega ,\tau \right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{ }+\rho {\xi }_{\sigma }{\sigma }_{\tau }\omega B\left(\omega ,\tau \right)+{\kappa }_{\lambda }\left(\theta -{\lambda }_{\tau }\right)C\left(\omega ,\tau \right)+\frac{1}{2}{\xi }_{\lambda }^{2}{\lambda }_{\tau }{C}^{2}\left(\omega ,\tau \right)\end{array}$ (33)

Rearranging the equation in terms of the state variables volatility and intensity results to

$\begin{array}{l}\left[\frac{\partial B}{\partial \tau }+\frac{1}{2}{\xi }_{\sigma }^{2}{B}^{2}\left(\omega ,\tau \right)+\rho {\xi }_{\sigma }\omega B\left(\omega ,\tau \right)-{\kappa }_{\sigma }B\left(\omega ,\tau \right)-\frac{1}{2}\omega +\frac{1}{2}{\omega }^{2}\right]{\sigma }_{\tau }\\ \text{ }+\left[\frac{\partial C\left(\omega ,\tau \right)}{\partial \tau }+\frac{1}{2}{\xi }_{\lambda }^{2}{C}^{2}\left(\omega ,\tau \right)-{\kappa }_{\lambda }C\left(\omega ,\tau \right)+\left(m\omega -{\int }_{-\infty }^{\infty }{\text{e}}^{\nu \omega }f\left(\nu \right)\text{d}\nu -1\right)\right]{\lambda }_{\tau }\\ \text{ }+\frac{\partial A\left(\omega ,\tau \right)}{\partial \tau }+{\kappa }_{\sigma }\left({\beta }^{\text{*}}+{\beta }^{rec}\right)B\left(\omega ,\tau \right)+{\kappa }_{\lambda }\theta C\left(\omega ,\tau \right)=0\end{array}$ (34)

where

${\int }_{-\infty }^{\infty }\left({\text{e}}^{\nu \omega }-1\right)f\left(\nu \right)\text{d}\nu ={\phi }_{\text{jump}}\left(\omega \right)$ (35)

is the moment-generating function of the jump size distribution and m is an average jump amplitude given earlier by

$m={\int }_{-\infty }^{\infty }{\text{e}}^{\nu }f\left(\nu \right)\text{d}\nu -1.$ (36)

Equating the coefficients of the term structures (the stochastic volatility and the intensity) in Equation (35) to zero, the following ordinary differential equations were obtained:

$\frac{\partial A\left(\omega ,\tau \right)}{\partial \tau }=-\left({\kappa }_{\sigma }\left({\beta }^{*}+{\beta }^{rec}\right)B\left(\omega ,\tau \right)+{\kappa }_{\lambda }\theta C\left(\omega ,\tau \right)\right)$ (37)

$\frac{\partial B\left(\omega ,\tau \right)}{\partial \tau }=-\frac{1}{2}{\xi }_{\sigma }^{2}{B}^{2}\left(\omega ,\tau \right)+\left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)B\left(\omega ,\tau \right)+\frac{1}{2}\left(\omega -{\omega }^{2}\right)$ (38)

$\frac{\partial C\left(\omega ,\tau \right)}{\partial \tau }=-\frac{1}{2}{\xi }_{\lambda }^{2}{C}^{2}\left(\omega ,\tau \right)+{\kappa }_{\lambda }C\left(\omega ,\tau \right)-\left(m\omega -{\int }_{-\infty }^{\infty }{\text{e}}^{\nu \omega }f\left(\nu \right)\text{d}\nu -1\right)$ (39)

The above systems of solutions were solved as follows.

Starting with the Equation (38) first, we can see that Equation (38) is a Ricatti differential equation. Setting

$B\left(\omega ,\tau \right)=2\frac{{D}^{\prime }\left(\tau \right)}{{\xi }_{\sigma }^{2}D\left(\tau \right)}$ (40)

in Equation (38) and simplify further, one will have a second order differential equation given by

$\frac{{\partial }^{2}D}{\partial {\tau }^{2}}+\left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)\frac{\partial D}{\partial \tau }+\frac{1}{4}{\xi }_{\sigma }^{2}\left(\omega -{\omega }^{2}\right)D\left(\tau \right)=0.$ (41)

The general solution of the Equation (41) is given by

$D\left(\tau \right)={k}_{1}\text{exp}\left(-\frac{1}{2}{\psi }_{-\tau }\right)+{k}_{2}\text{exp}\left(\frac{1}{2}{\psi }_{+\tau }\right)$ (42)

where

${\psi }_{±}=\left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)+\sqrt{{\left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)}^{2}-{\xi }_{\sigma }\left(\omega -{\omega }^{2}\right)}.$ (43)

At $\tau =0$, the following boundary conditions hold:

$\begin{array}{l}D\left(0\right)={k}_{1}+{k}_{2}\\ {D}^{\prime }\left(0\right)=\frac{1}{2}{\psi }_{+}{k}_{2}-\frac{1}{2}{\psi }_{-}{k}_{1}=0.\end{array}$ (44)

The sum

${\psi }_{+}+{\psi }_{-}=2\sqrt{{\left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)}^{2}-{\xi }_{\sigma }\left(\omega -{\omega }^{2}\right)}\equiv 2\varsigma$

with

$\varsigma =\sqrt{{\left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)}^{2}-{\xi }_{\sigma }\left(\omega -{\omega }^{2}\right)}$

and the product

${\psi }_{-}\cdot {\psi }_{+}=-{\xi }_{\sigma }\left(\omega -{\omega }^{2}\right).$

The constants ${k}_{1}$ and ${k}_{2}$ has values $\frac{{\psi }_{+}D\left(0\right)}{2\varsigma }$ and $\frac{{\psi }_{-}D\left(0\right)}{2\varsigma }$ respectively using the boundary conditions in Equation (47) and initial conditions

$A\left(\omega ,0\right)=B\left(\omega ,0\right)=C\left(\omega ,0\right)=0$

Substituting for ${D}^{\prime }\left(\tau \right)$ and $D\left(\tau \right)$ in Equation (40)

$B\left(\omega ,\tau \right)=2\frac{{D}^{\prime }\left(\tau \right)}{{\xi }_{\sigma }^{2}D\left(\tau \right)}$

and after simplification gives

$B\left(\omega ,\tau \right)=-\left(\omega -{\omega }^{2}\right)\frac{{\text{e}}^{-\frac{1}{2}{\psi }_{+\tau }}-{\text{e}}^{-\frac{1}{2}{\psi }_{-\tau }}}{{\psi }_{-}{\text{e}}^{-\frac{1}{2}{\psi }_{+\tau }}+{\psi }_{+}{\text{e}}^{-\frac{1}{2}{\psi }_{-\tau }}}$ (45)

$B\left(\omega ,\tau \right)=-\left(\omega -{\omega }^{2}\right)\frac{1-{\text{e}}^{-{\varsigma }_{\tau }}}{{\psi }_{-}+{\psi }_{+}{\text{e}}^{-{\varsigma }_{\tau }}}.$ (46)

In line with Artur [35] . Letting $A\left(\omega ,\tau \right)\equiv E\left(\omega ,\tau \right)+F\left(\omega ,\tau \right)$.

We rewrite Equation (37) as

$\frac{\partial A}{\partial \tau }=\frac{\partial E}{\partial \tau }\left(\omega ,\tau \right)+\frac{\partial F}{\partial \tau }\left(\omega ,\tau \right)$ (47)

such that

$\frac{\partial E}{\partial \tau }\left(\omega ,\tau \right)={\kappa }_{\sigma }{\beta }^{\text{*}}B\left(\omega ,\tau \right)+{\kappa }_{\sigma }{\beta }^{rec}B\left(\omega ,\tau \right)$ (48)

$\frac{\partial F}{\partial \tau }\left(\omega ,\tau \right)={\kappa }_{\lambda }\theta C\left(\omega ,\tau \right).$ (49)

Integrating the Equation (48) firstly, the following equations emerged.

$E\left(\omega ,\tau \right)={\int }_{0}^{\tau }{\kappa }_{\sigma }{\beta }^{\text{*}}B\left(\omega ,s\right)\text{d}s+{\int }_{0}^{\tau }{\kappa }_{\sigma }{\beta }^{rec}B\left(\omega ,s\right)\text{d}s$ (50)

$E\left(\omega ,\tau \right)={\kappa }_{\sigma }{\beta }^{\text{*}}{\int }_{0}^{\tau }B\left(\omega ,s\right)\text{d}s+{\kappa }_{\sigma }{\beta }^{rec}{\int }_{0}^{\tau }B\left(\omega ,s\right)\text{d}s$ (51)

$E\left(\omega ,\tau \right)=\frac{-2\kappa {\beta }^{\text{*}}}{{\xi }_{\sigma }^{2}}{\int }_{0}^{\tau }\frac{{D}^{\prime }\left(s\right)}{D\left(s\right)}\text{d}s+\frac{-2\kappa {\beta }^{rec}}{{\xi }_{\sigma }^{2}}{\int }_{0}^{\tau }\frac{{D}^{\prime }\left(s\right)}{D\left(s\right)}\text{d}s$ (52)

$E\left(\omega ,\tau \right)=\frac{-2\kappa \left({\beta }^{*}+{\beta }^{rec}\right)}{{\xi }_{\sigma }^{2}}{\mathrm{ln}D\left(s\right)|}_{s=0}^{s=\tau }$ (53)

$E\left(\omega ,\tau \right)=\frac{-2\kappa \left({\beta }^{*}+{\beta }^{rec}\right)}{{\xi }_{\sigma }^{2}}\mathrm{ln}\left[\frac{D\left(\tau \right)}{D\left(0\right)}\right].$ (54)

$E\left(\omega ,\tau \right)=\frac{-2\kappa \left({\beta }^{*}+{\beta }^{rec}\right)}{{\xi }_{\sigma }^{2}}\mathrm{ln}\left[\frac{{\psi }_{-}{\text{e}}^{-\frac{1}{2}{\psi }_{+\tau }}+{\psi }_{+}{\text{e}}^{-\frac{1}{2}{\psi }_{-\tau }}}{2\varsigma }\right].$ (55)

The result is finally written as

$E\left(\omega ,\tau \right)=\frac{-2{\kappa }_{\sigma }\left({\beta }^{*}+{\beta }^{rec}\right)}{{\xi }_{\sigma }^{2}}\left[{\psi }_{+\tau }+2\mathrm{ln}\left(\frac{{\psi }_{-}+{\psi }_{+}{\text{e}}^{-{\varsigma }_{\tau }}}{2\varsigma }\right)\right].$ (56)

Similarly, integrating the ODE in the Equation (49), one will have

$F\left(\omega ,\tau \right)=\frac{-{\kappa }_{\lambda }\theta }{{\xi }_{\lambda }^{2}}\left[{\psi }_{+\tau }+2\mathrm{ln}\left(\frac{{\chi }_{-}+{\chi }_{+}{\text{e}}^{-{ϵ}_{\tau }}}{2ϵ}\right)\right].$ (57)

Writing an explicit solution for $A\left(\omega ,\tau \right)$ as the sum of the solutions $E\left(\omega ,\tau \right)$ and $F\left(\omega ,\tau \right)$ in (56) and (57) respectively yield

$\begin{array}{c}A\left(\omega ,\tau \right)=-\left\{\frac{{\kappa }_{\sigma }\left({\beta }^{\text{*}}+{\beta }^{rec}\right)}{{\xi }_{\sigma }^{2}}\left[{\psi }_{+\tau }+2\mathrm{ln}\left(\frac{{\psi }_{-}+{\psi }_{+}{\text{e}}^{-{\varsigma }_{\tau }}}{2\varsigma }\right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{{\kappa }_{\lambda }\theta }{{\xi }_{\lambda }^{2}}\left[{\psi }_{+\tau }+2\mathrm{ln}\left(\frac{{\chi }_{-}+{\chi }_{+}{\text{e}}^{-{ϵ}_{\tau }}}{2ϵ}\right)\right]\right\}.\end{array}$ (58)

Analogously,

$C\left(\omega ,\tau \right)=2\text{Λ}\left(\frac{1-{\text{e}}^{-{ϵ}_{\tau }}}{{\psi }_{-}+{\psi }_{+}+{\text{e}}^{-{\varsigma }_{\tau }}}\right)$ (59)

where

$\Lambda \left(\omega \right)={\int }_{-\infty }^{\infty }{\text{e}}^{\nu \omega }f\left(\nu \right)\text{d}\nu -1-\omega \left({\int }_{-\infty }^{\infty }{\text{e}}^{\nu }f\left(\nu \right)\text{d}\nu -1\right)$ (60)

$ϵ=\sqrt{{\kappa }_{\lambda }^{2}-2{\xi }_{\lambda }^{2}\text{Λ}\left(\omega \right)}$ (61)

${\chi }_{±}=ϵ\mp {\kappa }_{\lambda }$ (62)

${\psi }_{±}=\mp \left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)+\sqrt{{\left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)}^{2}-{\xi }_{\sigma }\left(\omega -{\omega }^{2}\right)}$ (63)

but $\varsigma =\sqrt{{\left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)}^{2}-{\xi }_{\sigma }\left(\omega -{\omega }^{2}\right)}.$

The summary of the solution to the model (16) in which we have incorporated economic recession induced volatility parameter on the stock market price is given by the characteristic function approach (the core for Fourier transform) in sequel as

$M\left(\omega ,{x}_{\tau },{\sigma }_{\tau },{\lambda }_{\tau },\tau \right)={\text{e}}^{\left(r-q\right)\tau \omega +\omega {x}_{\tau }+A\left(\omega ,\tau \right)+B\left(\omega ,\tau \right){\sigma }_{\tau }+C\left(\omega ,\tau \right){\lambda }_{\tau }}.$ (64)

where

$\begin{array}{c}A\left(\omega ,\tau \right)=-\left\{\frac{{\kappa }_{\sigma }\left({\beta }^{\text{*}}+{\beta }^{rec}\right)}{{\xi }_{\sigma }^{2}}\left[{\psi }_{+\tau }+2\mathrm{ln}\left(\frac{{\psi }_{-}+{\psi }_{+}{\text{e}}^{-{\varsigma }_{\tau }}}{2\varsigma }\right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{{\kappa }_{\lambda }\theta }{{\xi }_{\lambda }^{2}}\left[{\psi }_{+\tau }+2\mathrm{ln}\left(\frac{{\chi }_{-}+{\chi }_{+}{\text{e}}^{-{ϵ}_{\tau }}}{2ϵ}\right)\right]\right\}\end{array}$

$B\left(\omega ,\tau \right)=-\left(\omega -{\omega }^{2}\right)\frac{1-{\text{e}}^{-{\varsigma }_{\tau }}}{{\psi }_{-}+{\psi }_{+}{\text{e}}^{-{\varsigma }_{\tau }}}$

$C\left(\omega ,\tau \right)=2\text{Λ}\left(\frac{1-{\text{e}}^{-{ϵ}_{\tau }}}{{\psi }_{-}+{\psi }_{+}+{\text{e}}^{-{\varsigma }_{\tau }}}\right)$

$\text{Λ}\left(\omega \right)={\int }_{-\infty }^{\infty }{\text{e}}^{\nu \omega }f\left(\nu \right)\text{d}\nu -1-\omega \left({\int }_{-\infty }^{\infty }{\text{e}}^{\nu }f\left(\nu \right)\text{d}\nu -1\right)$

$ϵ=\sqrt{{\kappa }_{\lambda }^{2}-2{\xi }_{\lambda }^{2}\text{Λ}\left(\omega \right)}$

${\chi }_{±}=ϵ\mp {\kappa }_{\lambda }$

${\psi }_{±}=\mp \left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)+\sqrt{{\left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)}^{2}-{\xi }_{\sigma }\left(\omega -{\omega }^{2}\right)}$

$\varsigma =\sqrt{{\left({\kappa }_{\sigma }-\rho {\xi }_{\sigma }\omega \right)}^{2}-{\xi }_{\sigma }\left(\omega -{\omega }^{2}\right)}.$

For a financial claim $g\left(x,\sigma ,\lambda ,\tau \right)$ satisfying the Partial Integro-Differential Equation (PIDE) given in Equation (30) such that the payoff function say $f\left({\text{e}}^{x},K\right)$ is satisfied by the claim, i.e. $g\left(x,\sigma ,\lambda ,0\right)=f\left({\text{e}}^{x},K\right)$. The Fourier transform of $g\left(x\right)$ simply put defined on the PIDE in the light of the definition defined in Equation (5) is given as

$\stackrel{^}{g}\left(w\right)=\mathcal{F}\left(g\left(x\right);w\right)={\int }_{-\infty }^{+\infty }g\left(x\right){\text{e}}^{-iwx}\text{d}x$ (65)

and the corresponding inverse Fourier transform using definition in (6) is given by

$g\left(x\right)={\mathcal{F}}^{-1}\left(\stackrel{^}{g}\left(w\right);x\right)=\frac{1}{\text{2π}}{\int }_{-\infty }^{+\infty }\stackrel{^}{g}\left(w\right){\text{e}}^{iwx}\text{d}w$ (66)

In light of Theorem 3.1 given by Artur Sepp [35] , and Theorem 3.2 given by Lewis (2001) [36] on page 11, we gave another modified version as follow.

Theorem 5: (Characteristic formula)

Let an asset price ${x}_{\tau }=\mathrm{ln}{S}_{\tau }$ possesses an affine analytic characteristic function ${\phi }_{\tau }\left(w\right)$ at time $\tau \le T$. Define a regularity strip

$Reg\text{ }S:=\left\{w:\alpha <\Im \left(w\right)<\beta \right\}$

where $\Im \left(w\right)$ is the imaginary part of w lying between $\alpha$ and $\beta$ in the regularity strip. Suppose further that ${\text{e}}^{-\Im \left(w\right)x}g\left(x\right)$ is defined on space of ${L}^{1}\left(R\right)$ such that $\stackrel{^}{g}\left(w\right),\Im \left(w\right)\in {S}_{g}$ where ${S}_{g}$ is the payoff strip and $\stackrel{^}{g}\left(w\right)$ obeys the Fourier transform given in Equation (65). Then the option value is given by

$g\left(x\left(t\right)\right)=\frac{{\text{e}}^{-r\left(T-t\right)}}{2\text{π}}{\int }_{i\Im \left(w\right)-\infty }^{i\Im \left(w\right)+\infty }{\phi }_{T}\left(-w\right)\stackrel{^}{g}\left(w\right)\text{d}w$ (67)

where $w\in Reg\text{ }S\cap {S}_{g}$.

Proof:

In a risk-neutral world Q, we expressed

$g\left(x\left(t\right)\right)={E}_{Q}\left[{\text{e}}^{-r\left(T-t\right)}g\left(x\left(t\right)|{\mathcal{F}}_{T}\right)\right],\text{\hspace{0.17em}}\text{\hspace{0.17em}}t (68)

$g\left(x\left(t\right)\right)={\text{e}}^{-r\left(T-t\right)}I{E}_{Q}\left[g\left(x\left(T\right)\right)\right]$ (69)

(70)

$g\left(x\left(t\right)\right)=\frac{{\text{e}}^{-r\left(T-t\right)}}{2\text{π}}{\int }_{i\Im \left(w\right)-\infty }^{i\Im \left(w\right)+\infty }{\text{e}}^{-iwx\left(T\right)}{E}_{Q}\left({\text{e}}^{-iwx\left(T\right)}\right)\stackrel{^}{g}\left(w\right)$ (71)

$g\left(x\left(t\right)\right)=\frac{{\text{e}}^{-r\left(T-t\right)}}{2\text{π}}{\int }_{i\Im \left(w\right)-\infty }^{i\Im \left(w\right)+\infty }{\text{e}}^{-iwx\left(T\right)}{\phi }_{T}\left(-w\right)\stackrel{^}{g}\left(w\right)\text{d}w$ (72)

Corollary 6 : (Characteristic formula for early exercise Option)

For an American option, consider a range of time $t\le \tau \le T$, where t is the initial (starting) time, $\tau$ is an early exercise time and T is the expiry time. Suppose one is able to determine an optimal payoff time ${\tau }^{\text{*}} in a stopping region, then the early exercise payoff of the claim $g\left(x\left(\tau \right)\right)$ is given as

$g\left(x\left(\tau \right)\right)=\frac{{\text{e}}^{-r\tau }}{2\text{π}}{\int }_{i\Im \left(w\right)-\infty }^{i\Im \left(w\right)+\infty }{\text{e}}^{-iwx\left(\tau \right)}{\phi }_{\tau }\left(-w\right)\stackrel{^}{g}\left(w\right)\text{d}w$ (73)

whenever $\tau ={\tau }^{\text{*}}$.

The proof follows immediately from the above theorem in addition to optimality condition attached to early exercising of American options whenever the early exercise time $\tau ={\tau }^{\text{*}}$.

The justification of the corollary rely on the fact that an American option can be exercised at any time up to the expiry time. By symmetric property, we assume that there exists ${\phi }_{\tau }\left(-w\right)$ since $w\in Reg\text{ }S\cap {S}_{g}$, and by extension the entire integrand exists whenever $w\in Reg\text{ }S\cap {S}_{g}$.

5. Numerical Fourier Based Transform of

A good number of references such as Carr & Madan (1999) [29] among others has agreed that if the corresponding characteristic function of a risk-neutral density is given or derived, the Fourier transform representation can be written analytically.

Let the characteristics function ${\phi }_{x}\left(iw,{X}_{\tau },\tau \right)$ of the moment generating function ${M}_{x}\left(w,{X}_{\tau },\tau \right)$ derived above holds for the affine model. In general, we rewrite

${\phi }_{\tau }\left(w\right)=E\left(\mathrm{exp}\left(iw{S}_{\tau }\right)\right)={\int }_{-\infty }^{\infty }\mathrm{exp}\left(iw{S}_{\tau }\right){q}_{\tau }\left({S}_{\tau }\right)\text{d}{S}_{\tau }$ (74)

where $\tau$ := Optimal time in the time horizon $\left\{0={t}_{0}<{t}_{1}<{t}_{2}<\cdots <{t}_{n}=T\right\}$.

In the case of European call ${E}_{T}^{call}\left(k\right)$ maturing at time $t=T$ with exercising price $K={\text{e}}^{k}$ defined on an underlying stock S. Suppose further that ${q}_{T}\left(s\right)$ represent the probability distribution function pdf of ${s}_{T}=\mathrm{ln}{S}_{T}$. Then the fair price ${E}_{T}^{call}\left(k\right)$ is expressed as the present value of the expected payoff defined by

${E}_{T}^{call}\left(k\right)={\text{e}}^{-rT}{\int }_{k}^{\infty }\left({\text{e}}^{s}-{\text{e}}^{k}\right){q}_{T}\left(s\right)\text{d}s.$ (75)

With reference to Carr & Madan (1999) [29] , we rewrite modified call to accommodate damping factor say $\alpha \in {R}^{+}$ so that ${c}_{T}\left(k\right)={\text{e}}^{\alpha k}×{E}_{T}^{call}\left(k\right)$.

Define ${\psi }_{\tau }\left(u\right)$ as the Fourier transformation representation of ${c}_{T}\left(k\right)$ such that

${\psi }_{T}\left(k\right)={\int }_{-\infty }^{\infty }{\text{e}}^{iuk}{c}_{T}\left(k\right)\text{d}k$. (76)

Then the call option price function is written as

${E}_{T}^{call}\left(k\right)=\frac{\mathrm{exp}\left(-\alpha k\right)}{2\text{π}}{\int }_{-\infty }^{\infty }{\text{e}}^{-iuk}{\psi }_{T}\left(u\right)\text{d}u$ (77)

${E}_{T}^{call}\left(k\right)=\frac{\mathrm{exp}\left(-\alpha k\right)}{\text{π}}{\int }_{0}^{\infty }{\text{e}}^{-iuk}{\psi }_{T}\left(u\right)\text{d}u$ (78)

${\psi }_{\tau }\left(u\right)$ in (78) is the Fourier transform of the call price given by

${\psi }_{\tau }\left(u\right)={\int }_{-\infty }^{\infty }\mathrm{exp}\left(iuk\right){c}_{\tau }\left(k\right)\text{d}k$ (79)

${\psi }_{\tau }\left(u\right)={\int }_{-\infty }^{\infty }\mathrm{exp}\left(iuk\right)×\mathrm{exp}\left(-rT\right){\int }_{k}^{\infty }\mathrm{exp}\left(\alpha k\right)\left({\text{e}}^{s}-{\text{e}}^{k}\right){q}_{T}\left(s\right)\text{d}s\text{d}k$ (80)

${\psi }_{\tau }\left(u\right)={\int }_{-\infty }^{\infty }\mathrm{exp}\left(-rT\right){q}_{T}\left(s\right){\int }_{-\infty }^{s}\left({\text{e}}^{s+\alpha k+iuk}-{\text{e}}^{\left(\alpha +1+iu\right)k}\right)\text{d}k\text{d}s$ (81)

${\psi }_{\tau }\left(u\right)={\int }_{-\infty }^{\infty }\mathrm{exp}\left(-rT\right){q}_{T}\left(s\right){\int }_{-\infty }^{s}\left({\text{e}}^{s+\left(\alpha +iu\right)k}-{\text{e}}^{\left(1+\alpha +iu\right)k}\right)\text{d}k\text{d}s$ (82)

Integrating the term in the second integral wrt k gives

${\psi }_{\tau }\left(u\right)={\int }_{-\infty }^{\infty }\mathrm{exp}\left(-rT\right){q}_{T}\left(s\right)\left[\frac{{\text{e}}^{\left(\alpha +1+iu\right)s}}{\alpha +iu}-\frac{{\text{e}}^{\left(\alpha +1+iu\right)s}}{\alpha +1+iu}\right]\text{d}s$ (83)

Further simplification yields

${\psi }_{\tau }\left(u\right)=\mathrm{exp}\left(-rT\right)\left[\frac{{\phi }_{\tau }\left(u-\left(\alpha +1\right)i\right)}{{\alpha }^{2}+\alpha -{u}^{2}+i\left(2\alpha +1\right)u}\right]$ (84)

By rationalising the base yields

${\psi }_{\tau }\left(u\right)=\mathrm{exp}\left(-rT\right)\left[\frac{{\phi }_{\tau }\left(u-\left(\alpha +1\right)i\right)×\left(\left({\alpha }^{2}+\alpha -{u}^{2}\right)-i\left(2\alpha +1\right)u\right)}{\left(\underset{a}{\underset{︸}{{\alpha }^{2}+\alpha -{u}^{2}}}+\underset{b}{\underset{︸}{i\left(2\alpha +1\right)u}}\right)\left(\underset{a}{\underset{︸}{{\alpha }^{2}+\alpha -{u}^{2}}}-\underset{b}{\underset{︸}{i\left(2\alpha +1\right)u}}\right)}\right]$ (85)

${\psi }_{\tau }\left(u\right)=\mathrm{exp}\left(-rT\right)\left[\frac{{\phi }_{\tau }\left(u-\left(\alpha +1\right)i\right)×\left(\left({\alpha }^{2}+\alpha -{u}^{2}\right)-i\left(2\alpha +1\right)u\right)}{\left(\underset{\left(a+b\right)}{\underset{︸}{{\alpha }^{2}+\alpha -{u}^{2}+i\left(2\alpha +1\right)u}}\right)\left(\underset{\left(a-b\right)}{\underset{︸}{{\alpha }^{2}+\alpha -{u}^{2}-i\left(2\alpha +1\right)u}}\right)}\right]$ (86)

From difference of two’s squares, the denominator reduces to

${\psi }_{\tau }\left(u\right)=\mathrm{exp}\left(-rT\right)×{\phi }_{\tau }\left(u-\left(\alpha +1\right)i\right)\left[\frac{\left({\alpha }^{2}+\alpha -{u}^{2}\right)-i\left(2\alpha +1\right)u}{\underset{{a}^{2}}{\underset{︸}{{\left({\alpha }^{2}+\alpha -{u}^{2}\right)}^{2}}}-\underset{{b}^{2}}{\underset{︸}{\left({i}^{2}{\left(2\alpha +1\right)}^{2}{u}^{2}\right)}}}\right]$ (87)

Since ${i}^{2}=-1$, then negative sign in the denominator changed to positive,

Hence we have

${\psi }_{\tau }\left(u\right)=\mathrm{exp}\left(-rT\right){\phi }_{\tau }\left(u-\left(\alpha +1\right)i\right)×\left[\frac{\left({\alpha }^{2}+\alpha -{u}^{2}\right)-i\left(2\alpha +1\right)u}{\underset{{a}^{2}+{b}^{2}}{\underbrace{{\left({\alpha }^{2}+\alpha -{u}^{2}\right)}^{2}+\left({\left(2\alpha +1\right)}^{2}{u}^{2}\right)}}}\right]$ (88)

a and b as used above are for ease simplification purpose and henceforth suppressed.

$\begin{array}{c}{\psi }_{\tau }\left(u\right)={\text{e}}^{-rT}{\phi }_{\tau }\left(u-\left(\alpha +1\right)i\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}×\left[\frac{\left({\alpha }^{2}+\alpha -{u}^{2}\right)-i\left(2\alpha +1\right)u}{\left({\alpha }^{4}+2{\alpha }^{3}+{\alpha }^{2}-2\left({\alpha }^{2}+\alpha \right){u}^{2}+{u}^{4}\right)+\left(4{\alpha }^{2}+4\alpha +1\right){u}^{2}}\right]\end{array}$ (89)

Finally,

${\psi }_{\tau }\left(u\right)={\text{e}}^{-rT}{\phi }_{\tau }\left(u-\left(\alpha +1\right)i\right)×\left(\frac{{\alpha }^{2}+\alpha -{u}^{2}-i\left(2\alpha +1\right)u}{{\alpha }^{4}+2{\alpha }^{3}+{\alpha }^{2}+\left(2\left({\alpha }^{2}+\alpha \right)+1\right){u}^{2}+{u}^{4}}\right)$ (90)

Substituting Equation (90) into (79), we now have a complete analytical formula for the option price as

$\begin{array}{c}{E}_{T}^{call}\left(k\right)=\frac{\mathrm{exp}\left(-\alpha k\right)}{\text{π}}{\int }_{0}^{\infty }{\text{e}}^{-iuk}{\text{e}}^{-rT}{\phi }_{\tau }\left(u-\left(\alpha +1\right)i\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}×\left(\frac{{\alpha }^{2}+\alpha -{u}^{2}-i\left(2\alpha +1\right)u}{{\alpha }^{4}+2{\alpha }^{3}+{\alpha }^{2}+\left(2\left({\alpha }^{2}+\alpha \right)+1\right){u}^{2}+{u}^{4}}\right)\text{d}u.\end{array}$ (91)

Equivalently stated as

$\begin{array}{c}{E}_{T}^{call}\left(k\right)=\frac{\mathrm{exp}\left(-\alpha k\right)}{\text{π}}{\int }_{0}^{\infty }{\text{e}}^{-\left(rT+iuk\right)}{\phi }_{\tau }\left(u-\left(\alpha +1\right)i\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}×\left(\frac{{\alpha }^{2}+\alpha -{u}^{2}-i\left(2\alpha +1\right)u}{{\alpha }^{4}+2{\alpha }^{3}+{\alpha }^{2}+\left(2\left({\alpha }^{2}+\alpha \right)+1\right){u}^{2}+{u}^{4}}\right)\text{d}u.\end{array}$ (92)

Extending this formulation to American option, we express the value of American call price as the sum of European call Fourier prices and Early exercise premium price as

${A}_{\left(t\le T\right)}={E}_{T}^{call}\left(k\right)+{P}_{t},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{where}\text{\hspace{0.17em}}k={\mathrm{log}}_{e}K.$ (93)

${P}_{t}$ represent the early premium price of the American call option such that ${\text{lim}}_{t\to T}{P}_{t}=0$, since at the maturity time, continuation region varnishes and there is no waiting time any longer. In other words, once the option reaches the maturity date, if it has not been exercised at an early time $t, then the option must be exercised at $t=T$. The implication is that the American option prices coincides with that of European option prices at the maturity date, T.

Hence,

$\begin{array}{c}{A}_{t}\left(k\right)=\frac{\mathrm{exp}\left(-\alpha k\right)}{\text{π}}{\int }_{0}^{\infty }{\text{e}}^{-\left(rT+iuk\right)}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}×\left(\frac{{\phi }_{\tau }\left(u-\left(\alpha +1\right)i\right)\left({\alpha }^{2}+\alpha -{u}^{2}-i\left(2\alpha +1\right)u\right)}{{\alpha }^{4}+2{\alpha }^{3}+{\alpha }^{2}+\left(2\left({\alpha }^{2}+\alpha \right)+1\right){u}^{2}+{u}^{4}}\right)\text{d}u+{P}_{t}.\end{array}$ (94)

To extend FFT algorithm which is a numerical technique of Discrete Fourier Transform (DFT) to American options valuation for the models under study with economic recession induced uncertainty, we define Fast Fourier Transform (FFT) equation as

$H\left(k\right)=\underset{j=1}{\overset{N-1}{\sum }}{\text{e}}^{-i\frac{2\text{π}}{N}\left(j-1\right)\left(k-1\right)}{y}_{j},\text{\hspace{0.17em}}\text{\hspace{0.17em}}1\le k\le N.$ (95)

Applying FFT algorithm to (94), we have

${A}_{\tau }\left({k}_{u}\right)\approx \frac{\mathrm{exp}\left(-\alpha k\right)}{\text{π}}{\sum }_{j=1}^{N}{\text{e}}^{-i{u}_{j}\left(-\varpi +\zeta \left(u-1\right)\right)}{\psi }_{T}\left({u}_{j}\right)\eta +{P}_{t},\text{\hspace{0.17em}}\text{\hspace{0.17em}}1\le u\le N.$ (96)

where $1\le u\le N$,${k}_{u}=-\varpi +\zeta \left(u-1\right)$ and $2\varpi =\zeta N$.

$\zeta$ is the size of consistent spacing existing between N values of log strikes K. Setting ${u}_{j}=\left(j-1\right)\eta$, and substituting into Equation (96) yields

${A}_{\tau }\left({k}_{u}\right)\approx \frac{\mathrm{exp}\left(-\alpha k\right)}{\text{π}}{\sum }_{j=1}^{N}{\text{e}}^{-i{u}_{j}\zeta \eta \left(j-1\right)\left(u-1\right)}{\text{e}}^{i\varpi {u}_{j}}{\psi }_{T}\left({u}_{j}\right)\eta +{P}_{t},\text{\hspace{0.17em}}\text{\hspace{0.17em}}1\le u\le N.$ (97)

The size of the consistent spacing $\zeta$ between the N values of k has the following relation

$\zeta =\frac{2\text{π}}{\eta N}\stackrel{\text{iff}}{⇔}\zeta \eta =\frac{2\text{π}}{N}.$ (98)

The relation is observable if we make comparison between Equations (95) and (97).

The smaller the values assumed by $\eta$, the better the fineness of the integration grid and vice versa.

Numerical Experiment

Consider an American Stock with initial price ${S}_{0}=100$, Strike Price $K=80$, Risk free interest rate $r=0.04$, Dividend rate $q=0.002$, Time to Maturity $T=1$ year, Volatilities ${\sigma }^{*}=0.1,{\sigma }^{rec}=\left\{0.00,0.025\right\}$, Integrability Parameter $\alpha =0.25$, Fineness of integration grid point $N={2}^{12}$. The call option value on the underlying American stock is reported below.

Table 1 above shows the payoff of a dividend paying American style call options during economic recession period and recession free period in Nigeria which the proposed volatility change formulation presented in this paper stressed. The option prices obtained via Fast Fourier Transform (FFT) approach yields better results compared to the values obtained by BSM and American call options solver prices. This confirmed numerical performance of Fast Fourier Transform (FFT) technique among other methods of pricing options such as Analytical payoff via BSM and American call option solver.

Figure 1 shows the true representation of the call options value obtained during recession free period in Table 1 above while Figure 2 shows the payoff of the

Figure 1. FFT call options payoffs comparison among other methods without recession induced volatility.

Table 1. Numerical value for a dividend paying American call options.

Figure 2. FFT call options payoffs comparison among other methods with recession induced volatility.

American style call options during recession period. The inference we can draw from both the Table 1 and the two Figure 1 and Figure 2 is that during a very short maturity period of a month in the options life span, the options value in both economic states (recession and recession-free) coincides while as the maturity period increases beyond a month, the upward volatility change induced by recession causes slight increase in the call options prices. If further experiment is carried out on put options, reverse will be the case in terms of volatility change effect on the put options payoff.

6. Conclusions

The notion of economic recession and its effect on volatility uncertainty on the payoff of European and American Options based on some certain assumptions was presented in this paper. Economic recession is becoming a global issue and not too far to be recognized as a re-occuring incident. In the history of economic recession, US is at the forefront as the continent has experienced several economic recessions. An intuition behind the introduction of economic recession induced volatility uncertainty in this research could be traced to Nigeria economic recession outbreak in 2016. This has an enormous effect on investors, financial institutions and every other economic activities in the country and by extension affected some aspect of foreign transactions (exchange), import and export etc. by the individual within and outside the country. In relation to the options (or stock) market, volatility has been a reference point to the level of fluctuation in the market price of the underlying asset. It is seen as a metric for the speed and amount of changes stochastically in the underlying asset’s prices. An investor’s acquaintance with volatility gives a better comprehension of why option prices behave in certain ways and this will guide them in decision making.

The model we investigated is proposed to be referred to as Uncertain Affine Exponential jump model with Recession induced Stochastic volatility and Intensity. Fourier transform of the model was implemented having derived an affine close form of the characteristic function. A numerical based Fourier transform algorithm called Fast Fourier transform (FFT) which was a variant of Carr & Madan [29] FFT algorithm for European call option valuation was adopted to determine the European call Fourier prices. The algorithm was further extended to American call options valuation by adding premium price to the European call options prices.

We also reported Nigerian Flourmill stock performance, during recession and recovery year in this study. The Stock volatility is seen as a measure of the uncertainty on the payoff or returns of the stock which may require an estimation. The two common methods in practice for an estimation of volatility of stocks or assets generally are known to be an estimation of historical volatility or implied volatility on that stock. Estimation of Recession induced volatility may be difficult to be determined accurately. Nevertheless, the fact we try to establish here is that the stocks price experiences a high level of uncertainty during the period of recession. The volatility of stocks tends to increase during recession compared to the period of normalcy. One major point is that volatility is never constant in an ideal real life situation. The Flourmill stock prices data used for calibration purpose revealed that during Nigerian recession, the stock prices became more volatile compared to other periods. According to the Assumption 1 above in this paper, we suggest the use of historical data prior to recession period to determine the level of uncertainty posed by Economic recession which we refer to as Economic recession induced volatility. We hope the figures we generated using MATLAB was presented in the Appendix section. Figures A1-A5 showed more details in terms of volatility change effect on stock prices induced by economic recession.

Acknowledgements

The authors would like to thank the editors and the anonymous reviewers for their useful comments and suggestions.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Bankole, P.A. and Ugbebor, O.O. (2019) Fast Fourier Transform Based Computation of American Options under Economic Recession Induced Volatility Uncertainty. Journal of Mathematical Finance, 9, 494-521. https://doi.org/10.4236/jmf.2019.93026

References

1. 1. Hsu, T. (2016) U.S. Recession Forecasting Using Probit Models with Asset Index Predictor Variables. 1-41.

2. 2. Liu, B. (2007) Uncertainty Theory. 2nd Edition, Springer-Verlag, Berlin.

3. 3. Liu, B. (2008) Theory and Practice of Uncertain Programming (Studies in Fuzziness and Soft Computing). 2nd Edition, Volume 239, Springer, Berlin, 1-177.

4. 4. Liu, B. (2010) Uncertainty Theory. 3rd Edition, Springer, Berlin, 1-148. https://doi.org/10.1007/978-3-642-13959-8http://orsc.edu.cn/liu/ut.pdf

5. 5. Liu, B. (2016) Uncertainty Theory. 5th Edition, Springer, Berlin, 1-148. http://orsc.edu.cn/liu/ut.pdf

6. 6. Ugbebor, O.O. (2011) Unpublished M.Sc. Lecture Notes on Probability Theory, Stochastic Processes and Applications in Mathematical Finance. University of Ibadan, Ibadan.

7. 7. Peter, O. and Andersson, M. (2011) Probability, Statistics, and Stochastic Processes. John Wiley & Sons, Inc., Hoboken. https://doi.org/10.1002/9781118231296

8. 8. Dungan, J.L., Gao, D. and Pang, A.T. (2001) Definitions of Uncertainty. 1-3.

9. 9. Bloom, N. (2014) Fluctuations in Uncertainty. Journal of Economic Perspectives, 28, 153-176. https://doi.org/10.1257/jep.28.2.153

10. 10. Bankole, P.A., Ojo, E.K. and Odumosu, M.O. (2017) On Recurrence Relations and Application in Predicting Price Dynamics in the Presence of Economic Recession. International Journal of Discrete Mathematics, 2, 125-131.

11. 11. Patrick, S. (2011) Market Uncertainty and Market Instability. Proceedings of the IFC Conference on Initiatives to Address Data Gaps Revealed by the Financial Crisis, Vol. 34, Basel, 25-26 August 2010, 430-435.

12. 12. Chen, L., Peng, J., Zhang, B. and Rosyida, I. (2017) Diversified Models for Portfolio Selection Based on Uncertain Semivariance. International Journal of Systems Science, 48, 637-648. https://doi.org/10.1080/00207721.2016.1206985

13. 13. Chen, L., Peng, J., Liu, Z. and Zhao, R. (2017) Pricing and Effort Decisions for a Supply Chain with Uncertain Information. International Journal of Production Research, 55, 264-284. https://www.tandfonline.com/toc/tprs20/current https://doi.org/10.1080/00207543.2016.1204475

14. 14. Chen, L., Peng, J. and Zhang, B. (2017) Uncertain Goal Programming Models for Bicriteria Solid Transportation Problem. Applied Soft Computing, 51, 49-59. https://www.sciencedirect.com/journal/applied-soft-computing/vol/51/suppl/C https://doi.org/10.1016/j.asoc.2016.11.027

15. 15. Lindström, E. (2010) Implications of Parameter Uncertainty on Option Prices. Advances in Decision Sciences, 2010, Article ID: 598103. https://doi.org/10.1155/2010/598103

16. 16. Lindstrvm, E. and Strålfors, J. (2012) Model Uncertainty, Model Selection and Option Valuation. 34th Symposium I Anvendt Statistik, Köpenhamn, 23-25 January 2012, 229-238.

17. 17. Rigotti, L. and Shannon, C. (2004) Uncertainty and Risk in Financial Markets.

18. 18. Stein, E. and Stein, J. (1991) Stock Price Distribution with Stochastic Volatiltiy: An Analytic Approach. Review of Financial Studies, 4, 727-752. https://doi.org/10.1093/rfs/4.4.727

19. 19. Alexander, D. and Pietro, V. (1999) Option Pricing with Uncertain Fundamentals: Theory and Evidence on the Dynamics of Implied Volatilities. 1-66. https://doi.org/10.2139/ssrn.199332

20. 20. Rama, C. (2006) Model Uncertainty and Its Impact on the Pricing of Derivative Instruments. Mathematical Finance, 16, 519-547. https://doi.org/10.1111/j.1467-9965.2006.00281.x

21. 21. Pfante, O. and Bertschinger, N. (2018) Uncertainty of Volatility Estimates from Heston Greeks. Frontier Applied Mathematics and Statistics, 3, 27. https://doi.org/10.3389/fams.2017.00027

22. 22. Bouchaud, J.P. and Potters, M. (2003) Theory of Financial Risk and Derivative Pricing. Cambridge University Press, New York, 379. https://doi.org/10.1017/CBO9780511753893

23. 23. Black, F. and Scholes, M. (1973) Valuation of Options and Corporate Liabilities. Journal of Political Economy, 81, 637-654. https://doi.org/10.1086/260062

24. 24. Heston, S.L. (1993.) A Closed-Form Solution for Options with Stochastic Volatility with Applications to Bond and Currency Options. The Review of Financial Studies, 6, 327-343. https://doi.org/10.1093/rfs/6.2.327

25. 25. Zhu, S., Badran, A. and Lu, X. (2012) A New Exact Solution for Pricing European Options in a Two States Regime Switching Economy. Computers and Mathematics with Applications, 64, 2744-2755. https://doi.org/10.1016/j.camwa.2012.08.005

26. 26. Hamilton, J. (1989) A New Approach to the Economic Analysis of Nonstationary Time Series and the Business Cycle. Econometrica, 57, 357-384. https://doi.org/10.2307/1912559

27. 27. Guo, X. (2001) Information and Option Pricing. Journal of Quantitative Finance, 1, 38-44. https://doi.org/10.1088/1469-7688/1/1/302

28. 28. Bollen, N. (1998) Valuing Options in Regime-Switching Models. Journal of Derivatives, 6, 38-49. https://doi.org/10.3905/jod.1998.408011

29. 29. Carr, P. and Madan, D.B. (1999) Option Valuation Using Fast Fourier Transform. Journal of Computational Finance, 2, 61-73. https://doi.org/10.21314/JCF.1999.043

30. 30. Sajib, B., Ruppa, K.T. and Parimala, T. (2004) Fast Fourier Transform for Option Pricing: Improved Mathematical Modeling and Design of an Efficient Parallel Algorithm. International Conference on Computational Science and Its Applications, Assisi, 14-17 May 2004, 686-695. https://doi.org/10.1007/978-3-540-24767-8_72

31. 31. Zhylyevskyy, O. (2010) A Fast Fourier Transform Technique for Pricing American Options under Stochastic Volatility. Economics Publications, Iowa State University Digital Repository, 1-24. https://doi.org/10.1007/s11147-009-9041-6

32. 32. Matthew, S.M. (2007) A Two-Asset Jump Diffusion Model with Correlation. Mathematical Modeling and Scientific Computing, Michaelmas, 1-57.

33. 33. Fei, W.-Y. (2014) On Existence and Uniqueness of Solutions to Uncertain Backward Stochastic Differential Equations. Applied Mathematics: A Journal of Chinese Universities, 29, 53-66. https://doi.org/10.1007/s11766-014-3048-y

34. 34. Huang, J., et al. (2014) Option Pricing Using the Fast Fourier Transform under the Double Exponential Jump Model with Stochastic Volatility and Stochastic Intensity. Journal of Computational and Applied Mathematics, 263, 152-159. https://doi.org/10.1016/j.cam.2013.12.009

35. 35. Artur, S. (2003) Fourier Transform for Option Pricing under Affine Jump-Diffusions: An Overview. http://www.hot.ee/seppar

36. 36. Lewis, A.L. (2001) A Simple Option Formula for General Jump-Diffusion and Other Exponential Levy Processes. Envision Financial Systems and OptionCity.net, 1-24.

Appendix

The following Figures A1-A5 were generated using Nigerian Flour mill stock prices of year 2016-2017 through MATLAB.

Figure A1. Nigerian flourmill stock prices movement during recession 2016 and recovery year 2017.

Figure A2. Comparison of stock prices variance during Nigeria economic recession.

Figure A3. Bar charts of flourmill stock prices movement during Nigeria economic recession and recovery year.

Figure A4. Flourmill stock price metrics during economic recession period.

Figure A5. Rotated side view of flourmill stock prices metrics during Nigerian economic recession 2016.