Fast Fourier Transform of Multi-Assets Options under Economic Recession Induced Uncertainties

A Fast Fourier transform approach has been presented by Carr & Madan (2009) on a single underlying asset. In this current research paper, we present fast Fourier transform algorithm for the valuation of Multi-asset Options under Economic Recession Induced Uncertainties. The issue of multi-dimension in both finite and infinite case of Options is part of the focus of this research. The notion of economic recession was incorporated. An intuition behind the introduction of recession induced volatility uncertainty is revealed by huge volatility variation during the period of economic recession compared to the period of recession-free. Nigeria economic recession outbreak in 2016 and its effects on the uncertainty of the payoffs of Nigeria Stocks Exchange (NSE) among other investments was among the motivating factors for proposing economic recession induced volatility in options pricing. The application of the proposed Fast Fourier Transform algorithm in handling multi-assets options was shown. A new result on options pricing was achieved and capable of yielding efficient option prices during and out of recession. Numerical results were presented on assets in 3-dimensions as an illustration taking Black Scholes prices as a bench mark for method effectiveness comparison. The key findings of this research paper among other crucial contributions could be seen in computational procedure of options valuation in multi-dimensions and uncertainties in options payoffs under the exposure of economic recession.


Introduction
In the history of options valuation, the work of Black-Scholes [1]  be forgotten no matter the shortcomings we may find in their work such as the issue of constant values assumed by the volatility and interest rates. Divers financial assets are traded in financial market. Some are risky assets while some are riskfree. Multi-asset options generally are a set of options whose pay-offs are governed by more than one underlying assets. A typical example of multi-asset option is a basket option. Basket option is based on two or more underlying assets or stocks. There is European type basket call or put option whereby the overall value of a basket of assets plays the role of the price of the underlying asset in single-asset European. Likewise, we have American style basket options as well. Other types of multi-assets options include: spread option, exchange options, Quanto options, Rainbow options, etc. Various types of multi-asset options could be described based on their payoff function. Beom, J. K. et al. (2013) [2] presents a simple numerical approach for American put options valuation.
The focus of most researchers is limited to two dimensions in the valuation of multi-asset options. Some authors such as Yuwei Chen [3] gave their formulation in terms of Black Scholes Partial Differential Equation (PDE) and considered numerical method for obtaining the prices. The authors, Carol Alexander et al., (2009) [4] did not relent in their effort as they presented an analytic approximation for Multi-asset Options Pricing but their approach was different from the one we gave here. In recent time, an existence result of solutions for a multi-dimensional stochastic control problem with singular state constraint was presented by Ulrich Horst and Xiaonyu Xia, (2018) [5]. The problem at hand is not centered on the specialty of the type of multi-asset but how to handle the nature of multi-assets amidst other conditional forces of uncertainty and dimensional challenges in the valuation process. According to Yuwei Chen [3], most multi-asset European options and all American options have no analytical formula giving their prices, therefore, numerical methods are employed to approximate such prices.
We consider challenges of uncertainty posed by instability of economic factors on the payoff of assets. The keen economic factor to us is recession. We consider volatility variation imposed on options prices during economic recession which we termed economic recession induced volatility in addition to other volatility uncertainties of options price from other sources. valuating an option on a single underlying asset to valuation of multi-assets options. Our major contribution is seen in valuation of multi-assets option in both finite and infinite dimensions and incorporation of economic recession factor on the volatility of assets during the state of economic recession.
The rest part of this paper is organized as follows. In Section 2, we gave preliminaries on Fast Fourier Transform. Section 3 is the Main result, which consists of subsections Fast Fourier Transform of multi-assets in finite dimension and Fast Fourier Transform of multi-assets in infinite dimension. Section 4 is the application on three assets in 3-dimensional case while Section 5 shows Numerical results on multi-assets in 3-dimensions. The last section is the conclusion part.

Some Definitions on Fourier Transform
Definition 1: Let : be a real-value function for x R ∈ , then the Fourier transform of f is defined as which belongs to L 2 -spaces.
Some certain axioms on Fourier transform exist in literature. Among such axioms are: convolution axiom, Translation axiom, Linearity axiom, Differentiation axiom, Symmetric axiom etc.

Some Axioms of Fourier Transform
• Convolution axiom: This property is considered very crucial as convolution of functions is mapped into multiplication in Fourier space. Let Then their convolution is defined by The Fourier transform of the convolution of the functions f and g with convolution operator  is given by The proof of Equation (4) is given in what follows. Proof: To this end, one will have . Similarly, ; The notion of Fourier transform is related to characteristic function in the following sense.
Then a characteristic function ( ), IR φ ω ω ∈ is defined as the Fourier transform of ( ) One can recover the probability distribution function ( ) P x by taking the inverse Fourier transform of the characteristics function just as the definition of inverse Fourier transform presented in Equation (2). This is given by The reader can see the reference ( [10]) and other standard texts on Fourier transform for more information.
There are available models for predicting economic recession probabilities especially on US recession (See [11]). The application of probability theory is found useful in such estimation [12]. Philip et al. [13] gave an intuition in predicting price dynamics in the presence of economic recession via recurrence relation technique and further presented a simple model showing the future of an economy under the exposure of recession if strong strategies and policies were not put in place to revive the economy. Fabio, Fornari and Antonio, Mele in (2009) [14] showed some connections between Financial Volatility and Economic Activity.
In sequel, we examine uncertainty on options value pose by recession on the volatility parameter of an option.
Assumption: Supposed that there are two sources of volatility variation on the assets. These are classified into economic recession induced volatility, rec v , and the usual volatility from other sources, * v , all assumed to be stochastic in nature. The volatility of the options is then given by where rec v is economy recession induced volatility in the stock market which arises due to information flow from the recessed Economy while * v is the volatility in the market from other sources. The study of parameters uncertainties in some financial models has received attention in the past few years. Two of those authors are [15] and [16]. Unfortunantely, many authors are yet to consider uncertaianties parameters that focus on economic recession in most Financial models use especially in derivative pricing. This present research is believed to stand out in that regards as it is hoped that more attention of authors will be drawn to economic recession and its effects in options pricing.

The Main Result
1) The Fast Fourier Transform (FFT) algorithm goes thus.
2) First get the characteristic function of the distribution. 2. Derive an analytic expression for the Fourier Transform of the modified options price wrt the characteristic function.
3) Find the Fourier inversion of the pricing function. 4) Perform Discrete Fourier transform in both step 2 and 3 above and apply Trapezoidal rule. 5) Choose uniform grid size for the computation of the FFT in step 4.
6) Make choices on the values assumed by other parameters such as N, the decay rate ( α ), etc. for optimality purpose.

Fast Fourier Transform of Multi-Assets in Finite Dimension
Consider an n-factor underlying asset , 1, 2, , , H S K of the options at maturity time on the given underlying multi-assets is defined as where the j S and j K are respectively asset and strike prices.
The parameter c is a control switch between call or put options. If c as defined in the Equations (19) and (20) where  is taken to be risk neutral measure while T q is the joint density function of each  . The integral representation of Equation (21) We gave the characteristic function corresponding to the joint density as In order to ensure the square integrability of the right hand side of the Equation (22), we take its product and a decaying term α over j k .
This we expressed as For the call option, the Fourier transform of the modified option prices is given by This is further simplified as follows For a known characteristic function or a calculated affine characteristic function T ϕ of a given process, it is possible to derive an analytic close form formula for ( ) Rewriting the integrals in series form using trapezoidal rule, equation (30)    , such that in an effort to capture the imaginary part of the complex form of j y , which is inclusive in Fourier transform, we have The options prices T V over log strikes from Equation (31)

Remarks:
1) The Equation (35) above is fast Fourier transform formula which can be adopted to calculate the value of European call option prices on finite -dimensional underlying multi-assets.
2) The corresponding European put option value can be obtained if we set the control parameter 1 c = − in the Equations (19)-(20) such that the sign attached to the log strikes and log assets changes in Equation (21) follow the same procedure through to arrive at (35). Alternatively, the put values can be obtained via put-call parity.

Fast Fourier Transform of Multi-Assets in an Infinite Dimension
Following similar algorithm presented in finite dimensional case, the above formulation can be extended to an infinite dimension as follows.

S K c S K c S K c S K c S K
where the j S and j K are respectively asset and strike prices. The parameter c is a control in switching between call or put options as discussed earlier.
where  still remain as risk neutral measure and T q is the joint density function for each We acknowledge Carr & Madan (1999) [6] for providing a fast Fourier Transform framework on single underlying asset which we extended to multi-assets in both finite dimension and infinite dimension in this paper. where the j S and j K are respectively asset and strike prices.

Application on Three Assets in 3-Dimensional Case
The option values computation in 3-dimension is given by where  is a risk neutral measure and T q is the joint density function of each The characteristic function corresponding to the joint density is given by , , e , 1, 2, 3.
Multiplying the right side of Equation (47) and decaying term say α over 1 2 3 , , k k k so as to have square integrable of the right side of (47).
We expressed the result as ( ) ( ) 1 1  2 2  3 3   1  2  3  1  2  3  1  2  3 , , e For the call option, the Fourier transform of the modified option prices is given by , , e , , d d d Whence, 3-dimensional version of , 1, 2, 3 j y j = , in (57) so as to capture the imaginary part of the complex form of j y , we have

Parameters Specification
Consider call option in European sense on an underlying vector of assets and Strikes prices given as For illustration purpose, we assumed the volatilities are stochastic and randn() function was used to generate 3 array of numbers between 1 and 3 for volatilities from the two sources explained by (18) in each case. The following results in the Tables were obtained.

Tables of Results
Using the same data for the parameters in Table 1 except that dividend parameter is set to zero. We obtain Table 2.

Conclusion
Fast Fourier Transform algorithm for valuation of multi-assets in both finite and infinite dimension was presented in this research work. The notion of economic recession induced volatility uncertainty was introduced in our formulation. An intuition behind the introduction of recession induced volatility uncertainty is revealed by huge volatility variation during the period of Economic recession compared to the period of normalcy (recession-free). Nigeria economic recession outbreak in 2016 and its effects on the payoffs uncertainty of Nigeria Stocks Exchange (NSE) among other invetments were among the motivating factors for proposing economic recession induced volatility formulation in Options pricing. This work generalizes the Fast Fourier transform algorithm of Carr & Madan (2009) [6] in the valuation of European call options on a single underlying asset. The algorithm is extended further to include some other steps, then applied to value multi-assets European options in finite and infinite dimensions. Numerical experiment was presented taking Black Scholes model as a bench mark. It is believed that this research work will enhance better understanding of Fast Fourier Transform valuation of multi-assets options in multi-dimensions. We gave the Options output metrics between the two methods. The result obtained is still comparable with little errors hinged on parameters uncertainties.