School of Computer Science and Engineering, University of Kurdistan Hewler, 40m St., Erbil, Kurdistan, Iraq
5.
Department of Electronics and Communication Engineering, Vemu Institute of Technology, Chithoor, India
6.
Department of Electronics and Communications Engineering and University Centre for Research & Development, Chandigarh University, Mohali-140413, Punjab
7.
Health Technology Research Institute, Amirkabir University of Technology (Tehran Polytechnic), Iran
8.
School of Automation and Electronic Information, Xiangtan University, Xiangtan 411105, China
Academic Editor: Yang Kuang
Received:
29 July 2023
Revised:
04 September 2023
Accepted:
11 September 2023
Published:
18 September 2023
The significance of discrete neural models lies in their mathematical simplicity and computational ease. This research focuses on enhancing a neural map model by incorporating a hyperbolic tangent-based memristor. The study extensively explores the impact of magnetic induction strength on the model's dynamics, analyzing bifurcation diagrams and the presence of multistability. Moreover, the investigation extends to the collective behavior of coupled memristive neural maps with electrical, chemical, and magnetic connections. The synchronization of these coupled memristive maps is examined, revealing that chemical coupling exhibits a broader synchronization area. Additionally, diverse chimera states and cluster synchronized states are identified and discussed.
Citation: Prasina Alexander, Fatemeh Parastesh, Ibrahim Ismael Hamarash, Anitha Karthikeyan, Sajad Jafari, Shaobo He. Effect of the electromagnetic induction on a modified memristive neural map model[J]. Mathematical Biosciences and Engineering, 2023, 20(10): 17849-17865. doi: 10.3934/mbe.2023793
Related Papers:
[1]
YongKyung Oh, JiIn Kwak, Sungil Kim .
Time delay estimation of traffic congestion propagation due to accidents based on statistical causality. Electronic Research Archive, 2023, 31(2): 691-707.
doi: 10.3934/era.2023034
[2]
Ming Wei, Shaopeng Zhang, Bo Sun .
Airport passenger flow, urban development and nearby airport capacity dynamic correlation: 2006-2019 time-series data analysis for Tianjin city, China. Electronic Research Archive, 2022, 30(12): 4447-4468.
doi: 10.3934/era.2022226
[3]
Yuxia Liu, Qi Zhang, Wei Xiao, Tianguang Chu .
Characteristic period analysis of the Chinese stock market using successive one-sided HP filter. Electronic Research Archive, 2023, 31(10): 6120-6133.
doi: 10.3934/era.2023311
[4]
Jiping Xing, Xiaohong Jiang, Yu Yuan, Wei Liu .
Incorporating mobile phone data-based travel mobility analysis of metro ridership in aboveground and underground layers. Electronic Research Archive, 2024, 32(7): 4472-4494.
doi: 10.3934/era.2024202
[5]
Moutian Liu, Lixia Duan .
In-phase and anti-phase spikes synchronization within mixed Bursters of the pre-Bözinger complex. Electronic Research Archive, 2022, 30(3): 961-977.
doi: 10.3934/era.2022050
[6]
Zhihui Wang, Yanying Yang, Lixia Duan .
Dynamic mechanism of epileptic seizures induced by excitatory pyramidal neuronal population. Electronic Research Archive, 2023, 31(8): 4427-4442.
doi: 10.3934/era.2023226
[7]
Li Li, Zhiguo Zhao .
Inhibitory autapse with time delay induces mixed-mode oscillations related to unstable dynamical behaviors near subcritical Hopf bifurcation. Electronic Research Archive, 2022, 30(5): 1898-1917.
doi: 10.3934/era.2022096
[8]
Jianmin Hou, Quansheng Liu, Hongwei Yang, Lixin Wang, Yuanhong Bi .
Stability and bifurcation analyses of p53 gene regulatory network with time delay. Electronic Research Archive, 2022, 30(3): 850-873.
doi: 10.3934/era.2022045
[9]
Jieqiong Yang, Panzhu Luo .
Study on the spatial correlation network structure of agricultural carbon emission efficiency in China. Electronic Research Archive, 2023, 31(12): 7256-7283.
doi: 10.3934/era.2023368
[10]
Ariel Leslie, Jianzhong Su .
Modeling and simulation of a network of neurons regarding Glucose Transporter Deficiency induced epileptic seizures. Electronic Research Archive, 2022, 30(5): 1813-1835.
doi: 10.3934/era.2022092
Abstract
The significance of discrete neural models lies in their mathematical simplicity and computational ease. This research focuses on enhancing a neural map model by incorporating a hyperbolic tangent-based memristor. The study extensively explores the impact of magnetic induction strength on the model's dynamics, analyzing bifurcation diagrams and the presence of multistability. Moreover, the investigation extends to the collective behavior of coupled memristive neural maps with electrical, chemical, and magnetic connections. The synchronization of these coupled memristive maps is examined, revealing that chemical coupling exhibits a broader synchronization area. Additionally, diverse chimera states and cluster synchronized states are identified and discussed.
1.
Introduction
In quantitative finance, volatility is a crucial factor in pricing derivatives, dynamic hedging, and portfolio management within the financial market. For FX (foreign exchange) options, prices are quoted in volatility. Given its significance in evaluating financial derivatives, volatility has been a focal point in both academia and practice.
For many years, academics have been studying how to model volatility. One of the most famous volatility models is the Black-Scholes model [1]. They described the dynamic of the asset price in terms of geometric Brownian model (GBM) in which volatility of the underlying asset is assumed to be a constant and the stock price has a log-normal distribution. Because of this, the Black-Scholes model has the advantage of being tractable, but has the disadvantage of not reflecting the time-varying volatility in the real market. In fact, in the real market, there is high probability that excess skewness and leptokurtosis occur as properties of the risk-neutral probability distribution unlike the log-normal distribution. Since then, to overcome these shortcomings of the Black-Scholes model, two major types of volatility models have been developed: Stochastic volatility (SV) models and local volatility models.
In SV models, taking account of many extraordinary volatility behaviors in the financial market after 1987 Financial Crash, the existence of a nonflat implied volatility has become remarkable. Therefore, participants interested in financial transactions have taken notice of the models that can predict the movement of financial assets, noting that the volatility of an underlying asset follows a stochastic process, and then the (pure) SV model has been proposed for describing and reflecting real situations in financial markets. In fact, the Heston model (cf. Heston [2]) and the fast mean-reverting SV model given by Fouque et al. [3] have become representative SV models, which are designed to capture the phenomenon of the mean-reversion of volatility in the real market. In addition, the Hull and White model [4] adopted the instantaneous variance process as a geometric Brownian motion. The mean-reverting Ornstein-Uhlenbeck process was also used in several models as a stochastic process (see Scott [5]; Chesney and Scott [6]; Schöbel and Zhu [7]). In addition, Heston model [2] assumed that the volatility is a Cox-Ingersoll-Ross (CIR) process and the Heston model is the one of most popular stochastic models because of the tractability.
The local volatility models were developed by Dupire [8] and Derman and Kani [9] for the continuous and discrete cases, respectively, which are also called non-parametric local volatility models. In this model, the volatility depends on the both asset price and time, addressing the importance of the correlation between the change in the underlying asset price and the randomness of volatility to price options. In addition, Cox and Ross [10] proposed the constant elasticity of variance (CEV) model as the parametric local volatility models. Especially, as seen in Tian et al. [11], the CEV model has been known to generate the U-shaped implied volatility contrary to the Black-Scholes model. However, in the CEV model mentioned by Cox and Ross [10], the movement in volatility and underlying asset price has perfect correlation, either positively or negatively, relying on the elasticity parameter. However, in the empirical study, there is definite correlation at all times, showing volatility is time-varying as shown in Ghysels et al. [12].
Recently, considering that the local volatility models or the SV models fail to capture the empirical evidence demonstrating that the implied volatility of equity options exhibits smile and skew curves simultaneously. A model that combines these two volatility models, called a hybrid model, has been proposed by many researchers, concentrating that the mixed models were designed to have the advantages of these two type of volatility models. First, one hybrid model was developed based on the non-parametric local volatility model and the stochastic model (see Van der Stoep et al. [13]; Tian et al. [11]; Cui et al. [14]), namely stochastic local volatility models. Second, another hybrid model that combines the CEV model and the SV model has also been developed (e.g., Andersen and Piterbarg [15]; Lord et al. [16]; Choi et al. [17]; Cui et al. [18]).
As volatility has a direct impact on pricing and hedging performance, several volatility derivatives have been introduced. In particular, these products enable investors to manage the volatility risk, that is, they are used to hedge volatility risk. Among the volatility derivatives, variance swap is the most famous [19,20]. The payoff of variance swap is calculated from the difference between the realized variance and a strike price of the contract. Starting with the studies of Carr and Madan [19] and Demeterfi et al. [21], who evaluated variance swap by replicating the swap payoff with a portfolio of European options, various studies have been introduced. Zhu and Lian [22] priced the variance swap under the Heston model. Zheng and Kwok [23] developed a probabilistic approach to evaluate the variance swap under an SV model with jumps in the underlying asset price. A closed-form pricing formula for the discrete-time variance swap under general affine GARCH (generalized autoregressive conditional heteroskedasticity) type models was introduced by Badescu et al. [20]. Issaka [24] evaluate the variance swap by assuming that the underlying assets follow the multiple SV models. Xi and Wong [25] priced the variance swap in a rough Heston model [26] and investigated the sensitivity to roughness on the variance swap.
In addition, as another volatility derivative, there is a timer option. The payoff of the timer option is similar to that of the vanilla option, but it has a random maturity contrary to the vanilla option, whose payoff is determined by only the fixed maturity embedded in the option. The timer option expires when the accumulated realized variance of the underlying asset reaches a predetermined level. In particular, timer options were first studied by Neuberger [27] as the "mileage option", another name for timer options, in academia and launched by the Société Générale Corporate and Investment Banking (SG CIB) in 2007. Since April 2007, timer options have been traded on the market (see [28]). After the trade began, several studies have been conducted on the evaluation of time options. Bick [29] considered the continuous version of timer options, and Li [30] implemented the Monte Carlo simulation to price timer options. Saunders [31] investigated the pricing of timer options under the fast mean-reverting SV and provided the closed-form approximation formulas. Liang et al. [32] obtained closed-form pricing formulas for timer options under the 3/2 and Heston models by using the path-integral techniques. Bernard and Cui [33] developed a one-dimensional problem for pricing of a timer call option under a general SV model. In addition, they showed the empirical results in the Hull-White and Heston models. Li [34] studied perpetual timer options under general SV models. Thus, they obtained pricing formulas for the Heston model and the 3/2 model. Ma et al. [35] derived a fast approximate analytic method to evaluate timer options under the Vasicek interest rate model. Li [36] applied the Heston model to price the timer options. In addition, as mentioned in Zhang et al. [37], the perpetual timer option was examined based on the Hull-White SV model. [38] developed concise pricing formulas for timer options using a probabilistic method to explicitly derive probability densities in stochastic volatility models, including Heston and 3/2. Recently, Wang et al. [39] employed multiscale volatility models to price timer options, stressing that the factor of the volatility is assumed to follow both fast mean-reverting and slowly varying factors. In particular, [40] introduced the return timer option, a financial derivative with a random expiry triggered by the first occurrence of a return exceeding upper or lower barriers.
In this study, we investigate the timer option prices under the CEV model. In previous studies, the pricing of the time options has been only evaluated using various stochastic volatility models. To the best of our knowledge, there has been no research about the timer option pricing based upon CEV model in which the volatility process is considered to be the SV given by Heston [2] or Fouque et al. [3]. In fact, since the model dynamics of the underlying asset for the classical timer option is a form of the SV model, including the fast-mean reverting process given by Fouque et al. [3] or CIR process described by Heston [2], the pricing procedure of the timer options based upon CEV model is very similar to that of European options under the stochastic and local volatility (SVCEV) model introduced by Choi et al. [17], whose model is driven by the fast-mean reverting volatility process and constant elasticity of variance process. According to Choi et al. [17], the SVCEV model shows some improvements over the traditional CEV model. The SVCEV model has been widely used to evaluate various contingent claims. For example, Bock et al. [41] developed the pricing formula of a European barrier option under the SVCEV model. In addition, Kim et al. [42] applied the SVCEV model to evaluate the real options. Furthermore, several extended models based on SVCEV have been developed. For instance, Choi et al. [43] introduced a multiscale hybrid model consisting of fast and slow factors and evaluated an equity-linked annuity under the multiscale hybrid model. The previous SVCEV model [17] derived an approximation solution using up to the first order correction. Recently, an extended model has been developed using up to the second order correction term [44].
In this article, by using asymptotic analysis on the partial differential equations (PDEs) for timer options, we derive the approximation formulas of timer options with the CEV model. Later, we call the timer options under the CEV model TO-CEVs. To demonstrate that our approximated solution for the timer options has been obtained accurately and effectively, we compare the option value with a Monte-Carlo price. Moreover, in the numerical implications, we analyze the price sensitivity of the SV for the given elasticity value on the timer options. Above all else, the most remarkable result of our paper is that the standard European option under CEV may or may not be overvalued compared with timer option with CEV, depending on the elasticity value, which is contrary to the results of Sawyer [28] emphasizing that the price of the European vanilla call options can be quite overvalued compared to the price of timer call options in the region of in-the-money (ITM).
The remainder of this paper is organized as follows: In Section 2, we model stochastic differential equations (SDEs) to price a timer option under CEV models. Next, in Section 3, we derive the approximation solution of the option value taking advantage of asymptotic analysis. In Section 4, we examine the pricing accuracy of the approximated formula for the timer options through a Monte-Carlo simulation and analyze the parameter sensitivity on the first approximated option price. Finally, Section 5 provides concluding remarks.
2.
Model formulation
In this section, we provide the stochastic dynamics for the price of a timer call option based on hybrid stochastic and local volatility. According to Choi et al. [17], we assume that the underlying asset price model, Xt, and its instantaneous variance process, Yt, satisfy the following stochastic differential equation (SDE) model under a real probability measure P:
dXt=μXtdt+f(Yt)Xθ/2tdWXt,
(2.1)
dYt=α(m−Yt)dt+βdWYt,
(2.2)
where μ is the expected return rate of the underlying asset process, Xt, f(⋅) is a smooth function satisfying 0<c1≤f≤c2<∞ for some constants c1 and c2, θ is an elasticity parameter, α and β are positive constants, and m is the long-run mean level of Y. In addition, WXt and WYt are standard Brownian motions such that d⟨WX,WY⟩t=ρdt, where |ρ|≤1.
Herein, based on the fact that the volatility of the underlying asset has been demonstrated to be mean-reverting through data-driven study [45], we construct the model of the volatility as a fast mean-reverting Ornstein-Uhlenbeck process, as depicted in SDE (2.2). Since the ergodic process Yt is explicitly presented by Yt=m+(Y0−m)e−αt+β∫t0e−α(t−s)dWYt, the process Yt follows the normal distribution. i.e., Y∼N(m+(Y0−m)e−αt,β22α(1−e−2αt)). Interestingly, as t goes to ∞, its invariant distribution is given by N(m,β22α), which is independent of the initial volatility level, Y0. Nevertheless, empirical studies grounded in financial market data (refer to Choi et al. [17] and Choi et al. [44]) indicate that a (pure) stochastic volatility process (in the case of θ=2) may not accurately capture the geometric properties of the implied volatility, specifically, the implied-volatility-smile and implied-volatility-skew phenomena. Therefore, to overcome the aforementioned limitations observed in the SV model, we consider the timer options based upon the CEV diffusion, which is very similar to that of European options under stochastic and local volatility (in brief, SVCEV), proposed by Choi et al. [17].
The most distinctive feature of the SVCEV model proposed by Choi et al. [17] lies in the significant impact of the elasticity parameter θ. This economic parameter, when set to θ=2, represents the classical Black-Scholes market, whereas for θ>2, it depicts the phenomenon noticed in commodity markets (especially, coal, copper, or gold [46]) where an increase in commodity prices leads to a concurrent rise in volatility, known as the inverse leverage effect. On the other hand, when θ<2, it captures the leverage effect commonly observed in stock markets, where a fall in stock prices results in an increase in the underlying asset's volatility. Therefore, the SVCEV model of volatility for the underlying asset, as adopted from Choi et al. [17], can be considered as a suitable and realistic volatility model, efficiently capturing the random structures across various financial markets.
Then, by the Girsanov theorem in Øksendal [47], under the risk-neutral measure Q, the stochastic dynamics in (2.1) and (2.2) are transformed into
dXt=rXtdt+f(Yt)Xθ/2td¯WXt,
(2.3)
dYt={1ϵ(m−Yt)−ν√2√ϵΛ(Yt)}dt+ν√2√ϵd¯WYt,
(2.4)
where r is risk-free interest rate, ν is the standard deviation of the invariant distribution of Y, and Λ(Yt) is the market price of volatility risk. In addition, ¯WX and ¯WY are transformed Brownian motions such that ⟨¯WX,¯WY⟩t=ρdt.
We now demonstrate the timer call option and derive its PDE. Let us define the accumulated variance process at time t as
It=∫t0f2(Yu)Xθ−2udu.
We note that if θ=2, then it reduces to the form of the accumulated variance process found in standard timer options. Next, if we denote a pre-determined variance budget as B, depending on the option buyer's choice, then the stopping timer or random maturity, τB, is given by
τB=inf{t>0,It=B}.
(2.5)
In other words, the expiration of the timer option coincides with the initial moment at which the accumulated variance budget, It, reaches the predetermined variance budget level, B.
Then, the payoff of a timer call option is max(XτB−K,0), where K is the strike price. Therefore, we obtain the price of the timer call option at time t as follows:
where EQ is the expectation with respect to risk-neutral measure Q.
Now, by using the well-known Feynman-Kac formula and employing the well-known fact, as described in Li [36], that the price of timer options does not depend on the time variable t, we can transform this conditional expectation (2.6) into the following partial differential equation:
3.
Asymptotic analysis and first-order approximated solution
In this section, we aim to build an approximate solution utilizing the asymptotic analysis proposed by Fouque et al. [45]. Referring to Fouque et al. [45], for small parameter ϵ, we expand the original solution P(x,y,I) with respect to √ϵ in the following form:
P(x,y,I)=P0(x,y,I)+√ϵP1(x,y,I)+ϵP2(x,y,I)+⋯,
(3.1)
where P0,P1,⋯ are functions of (x,y,I) satisfying P0(x,B)=(x−K)+ and Pn(x,B)=0 for all n≥1. Now, if we substitute the expansion in (3.1) into the singular perturbed PDE presented in (2.7), we obtain
Now, we present two useful and fundamental results, known as the growth condition and the centering condition, for asymptotic analysis.
Lemma 3.1.(Growth condition) As y→∞, if the leading-order price P0(x,y,I) and the correction term P1(x,y,I) do not grow as much as ∂P0∂y∼ey2/2 and ∂P1∂y∼ey2/2, respectively, then P0 and P1 are independent of the unobservable variable y.
Proof. The proof is similar to that in Theorem 3.1 of Ha et al. [48] or Lemma 1 and Theorem 4.1 in the work of Choi et al. [17]. From (3.2), we have the following PDEs:
L0P0=0,L0P1+L1P0=0.
(3.3)
Next, the ODE L0P0=0 in (3.3) yields
P0(x,I)=c1(x,I)∫y0e(m−u)22ν2du+c2(x,I)
for some functions c1(x,I) and c2(x,I). Then, by using the assumption, as y→∞, if the leading-order price P0 does not grow as much as ∂P0∂y∼ey2/2, then c1(x,I)=0. Therefore, P0 is independent of the unobservable y. Similarly, if we apply this assumption to L0P1+L1P0=0, the correction term P1 is independent of y. □
Lemma 3.2.(Centering condition) If the ODE L0G(y)+G=0, a Poisson equation for the function G(y), has a unique solution, then the following centering condition must be satisfied:
⟨G⟩≡∫RG(y)Φ(y)dy=0,
where Φ(y)=1√2πνe−(y−m)22ν2 is the invariant density distribution of the ergodic process Y.
Proof. Refer to Section 5.2 in Fouque et al. [45]. □
Proposition 3.1.Under the growth condition presented in Lemma 3.1, the leading-order price P0(x,I) satisfies the following PDE:
Proof. According to Lemma 3.2, since L0P2+L2P0=0 is a Poisson equation of P2 with respect to y and P0 is independent of the variable y, the following centering condition must be hold:
⟨L2P0⟩=⟨L2⟩P0=0,
where
⟨L2⟩=ˆσ2xθ−2∂∂I+ˆσ22xθ∂2∂x2+r(x∂∂x−⋅)≡LTO−CEV(ˆσ)
and the effective volatility ˆσ is given by
ˆσ=√⟨f2⟩.
Therefore, we can obtain the PDE (3.4) with (3.5). □
Now, we provide the analytic form solution to the PDE presented in Proposition 3.1.
Theorem 3.1.The leading-order price, P0(x,I), is given by
where the modified Bessel function of the first kind of order q, denoted as Bq(x), is given by
Bq(x)=∞∑n=0(x2)2n+qr!Γ(n+1+q).
Proof. If we define the new state variable ξ as ξ=Iˆσ2xθ−2, then ∂P0∂I=1ˆσ2xθ−2∂P0∂ξ holds. Therefore, the PDE problem in (3.4) and (3.5) can be transformed into
∂P0∂ξ+ˆσ22xθ∂2P0∂x2+r(x∂P0∂x−P0)=0,
(3.7)
P0(x,Bˆσ2xθ−2)=max(x−K,0).
(3.8)
Surprisingly, this PDE is related to the pricing formula of a European call option under the CEV diffusion, except that the current time t and expiry date T are now replaced by the initial accumulated variance budget ξ and variance budget B, respectively. Therefore, referring to Lipton [49], the solution of PDE (3.7) with boundary condition (3.8) is given by
Finally, if we replace ξ with I and use the change of variables ξ=Iˆσ2xθ−2, then we can obtain the leading-order price P0 as in (3.6). □
Next, we continue our asymptotic approaches to obtain an explicit expression for correction term P1, which plays a key role in this research and will reflect the impact of stochastic volatility on timer prices.
Proposition 3.2.Under the growth condition given in Lemma 3.1, the first-order correction term P1(x,I) satisfies the following non-homogeneous PDE:
LTO−CEV(ˆσ)P1(x,I)=H(x,I;θ),
(3.9)
P1(x,B)=0,
(3.10)
where the non-homogeneous term H(x,I;θ) is given by
Proof. As shown in Theorem 3.1, if we define the new state variable ξ as ξ=Iˆσ2xθ−2, then ∂P1∂I=1ˆσ2xθ−2∂P1∂ξ holds. Therefore, the PDE problem in (3.9) and (3.10) can be transformed into
∂P1∂ξ+ˆσ22xθ∂2P1∂x2+r(x∂P1∂x−P1)=0,
(3.13)
P1(x,Bˆσ2xθ−2)=0.
(3.14)
Now, using the Feynman-Kac formula and referring to Choi et al. [17], the solution to the PDE (3.13) with the source term, H(x,I;θ), as defined in Proposition 3.2, and with a zero boundary condition (3.14), is given by
P1(x,ξ)=EQ[−∫Bˆσ2xθ−2ξe−r(s−ξ)H(x,s)ds|Xt=x].
Finally, if we replace ξ with I, then we can obtain the analytic-form solution P1(x,I) as described in (3.12). □
4.
Accuracy and numerical experiments
This section investigates the error of the approximation of the timer call option price P(x,y,I) given by (2.6) under the CEV model and examines the price sensitivities on the TO-CEV with regard to the model parameters.
4.1. Accuracy of the approximation
If we combine the results of leading-order term P0(x,I) presented in Theorem 3.1 and the correction term P1(x,I) described in Theorem 3.2, then we can obtain a first-order approximation to the TO-CEV price P(x,y,I), which is denoted by ˜Pϵ, given by
P(x,y,I)≈P0(x,I)+√ϵP1(x,I)≡˜Pϵ(x,I).
(4.1)
Then, one may want to examine the accuracy of the above approximation ˜Pϵ(x,I). In this regard, by referring to Fouque and Lorig [50], if we assume that the payoff function and its derivatives are continuously differentiable and bounded, utilizing the asymptotic expansions methods, then the accuracy of ˜Pϵ(x,I) is given by
|P(x,y,I)−˜Pϵ(x,I)|≤O(ϵ).
(4.2)
However, for a non-smooth payoff, that is, the payoff function is not continuously differentiable, it is required to introduce a regularization argument because of the singularity at x=K. In this research, we numerically investigate the accuracy of the analytic form solution to the PDE (2.7) via Monte-Carlo simulations (refer to Table 1 in Section 4.2.1), instead of providing the mathematical proof for the accuracy.
Table 1.
For given θ=1.8,1.9,2.1 and θ=2.2, solutions ˜Pϵ and errors are given against ϵ, where "R.Error" is the relative percentage error. Note that the baseline parameters are given by X0=100,T=1.0,r=0.03,K=100,B=0.3,v=0.05,⟨Λ(y)ψ′(y)⟩=0.0141,⟨f(y)ψ′(y)⟩=0.0023,ρ=−0.1,m=0.5,u=1.0,f(y)=ey and ⟨f(y)2⟩=0.0272.
In this section, we provide numerical implications by using the asymptotic approximations (4.1) of the TO-CEV prices. In addition, we compare the option price calculated by the closed-form with that obtained from the Monte-Carlo method. The parameters we have chosen in the numerical experiments are described by X0=100,T=1.0,r=0.03,K=100,B=0.3,v=0.05,⟨Λ(y)ψ′(y)⟩=0.0141,⟨f(y)ψ′(y)⟩=0.0023,ρ=−0.1,m=0.5,u=1.0,f(y)=ey and ⟨f(y)2⟩=0.0272, referring to Ha et al. [48].
4.2.1. Accuracy and Monte-Carlo simulation
We show the accuracy of the model formula by comparing the price of the timer option, presented by the model formula and the option value through Monte Carlo (MC) simulation. Table 1 demonstrates a comparison of the results. In the table, PMC represents the Monte Carlo price and ˜Pϵ is the approximate option price given by (4.2). In addition, we calculate the relative percentage error given by 100×|PMC−˜Pϵ(x,I)PMC| for the numerical implications. According to Table 1, by performing simulations 10000, 50000, or 100000 times, we obtain stable MC simulation results. Furthermore, we confirm that the value of relative error approaches zero as ϵ goes to zero. The average CPU times (in seconds) required for the execution of PMC are 27.5863, 4894.1, and 34745 for simulations conducted 10000, 50000, and 100000 times, respectively. In contrast, the CPU time for ˜Pϵ is 0.6485. Consequently, our options' formula verifies the excellence of the model not only with regards to the computational accuracy but also in terms of the computational efficiency.
4.2.2. Sensitivity to the parameters
Figure 1 displays the changes in the correction term price P1 with respect to the accumulated variance under the variance budget for a given θ value. As seen in both figures, as the value of the accumulated variance is less than the variance budget (B), the impact of the correction price is significant, but when the v is greater than the B, the values of the correction term are zero and have no influence on the option price. In fact, in the case of timer options, as the accumulated variance reaches the variance budget, the options get exercised. Therefore, when the aggregated variance is larger than the budget level, the effect of the stochastic volatility (SV) must be zero, and then the correction term price remains zero. Moreover, it can be seen that the greater the variance budget level, the lower the price impact of the correction term is for the accumulated variance. It implies that as the time-to-maturity of the timer options is higher (as the budget levels get bigger), the price sensitivities become smaller, which ultimately decreases the influence of the stochastic volatility (SV) on the timer option. What's more remarkable is that, as shown in Figure 1, the more the elasticity value increases, the more the price of the correction term is sensitive to the variance budget. It means that the effect of the SV on the timer option tends to rise as the elasticity parameter gets larger, and in particular, it becomes more prominent as the expiration is longer for the larger elasticity value.
Figure 1.
Value of the correction term Pϵ1 in terms of the accumulated variance v for the given θ. Note: X0=100,T=1.0,r=0.03,K=100,⟨Λ(y)ψ′(y)⟩=0.0141, ⟨f(y)ψ′(y)⟩=0.0023, ρ=−0.1,ϵ=0.001,m=log(0.1),u=1.0,f(y)=exp(y), ⟨f(y)2⟩=0.0272.
Figure 2 shows the changes in the correction term P1 of the timer option with constant elasticity of variance (CEV) with regard to the underlying asset for each strike price and for a given elasticity parameter. As seen in Figure 2(a)–(d), the correction prices have a tendency to decrease sharply as the underlying asset is near the strike price (K) regardless of the value of the elasticity θ. It implies that the SV has a more significant impact on the time option if the value of the underlying asset is adjacent to the strike value (K) for the elasticity parameter θ, especially, exhibiting a hump phenomenon for the graph of the correction term for θ>2.
Figure 2.
Value of the correction term Pϵ1 with respect to the underlying asset x for the given θ. Note: B=0.03,v=0.01. The values of the remaining parameters are the same as in Figure 1.
Figure 3 presents the comparison of the price of the timer option (TO) and the price of the timer CEV options (TO-CEVs). Taking account of the fact that the classical timer option has the elasticity θ=2, note that the price of TO-CEVs tends to decrease as the elasticity value gets bigger, and the sensitivities of the option value get larger when the elasticity decreases.
Figure 3.
Comparison of the values of timer CEV option (TO-CEV) with respect to the underlying asset x for a given θ. Note: B=0.3,v=0.05. The values of the remaining parameters are the same as in Figure 1.
Figure 4 exhibits the changes of the TO-CEV prices against the underlying asset for a given elasticity value, comparing them with those of the European CEV options (EO-CEVs). As shown in Figure 4(a)–(f), one can observe that the price gap between the TO-CEVs and EO-CEVs depends on the value of the elasticity. We can see that, for θ>2 or θ=2, definitely, the prices of EO-CEVs tend to be bigger than that of the TO-CEVs, and as the elasticity parameter increases, the phenomenon is more pronounced, demonstrating that our findings that TO-CEVs are noticeably underestimated compared with EO-CEVs are consistent with those of Sawyer [28]. However, in the case of θ<2, the price difference between TO-CEVs and EO-CEVs becomes smaller as the elasticity value grows down-after all, the TO-CEVs tend to be overpriced compared to the EO-CEVs in terms of underlying asset price, contrary to the results of Sawyer [28]. In particular, as seen in Figure 4(c), the EO-CEVs and the TO-CEVs have almost the same values in the region of deep ITM. In fact, referring to the empirical evidence verifying that the elasticity value calibrated from the historical data analysis of the volatility of the S&P 500 index is between 1.9 and 2, given by Choi et al. [17] and Figure 4(c), (d), it can be observed that the price of the EO-CEVs tend to get larger than that of the TO-CEVs in the real financial situation, especially, in the area of ITM. It implies that, even in the case of θ<2, if we take account of the elasticities obtained in the real situation, the value of TO-CEVs becomes undervalued compared to that of EO-CEVs in the region of ITM, similar to the findings described by Sawyer [28].
Figure 4.
Value of Vanilla and timer option prices by CEV models with respect to the underlying asset x for the given θ. Note: B=0.022,v=0.02. The values of the remaining parameters are the same as in Figure 1.
In this paper, we investigated the valuing problem of timer options under the constant elasticity of variance (CEV). In the real financial market, timer options are one of the financial derivatives that take account of the level of volatility so that the options are expired as the realized variance arrives at the variance budget. The purpose of this research is to analyze the price sensitivities between the timer options with the CEV (TO-CEVs) and the standard European CEV options (EO-CEVs).
First, under the partial differential equations (PDEs) obtained from the underlying asset model and Feynman-Kac formula, we derive the approximated analytic solutions for the options by making use of the method of asymptotic analysis given by Fouque et al. [3]. Second, we implemented the numerical experiments to verify the pricing accuracy of our analytic-form formulas for the timer CEV options, compared to Monte Carlo prices. Third, we observed the pricing sensitivities of the options in terms of some model parameters, emphasizing the impact of the SV on the timer options in terms of the elasticity parameter. Fourth, we found that the timer CEV options tend to be underpriced compared with the standard European CEV options for θ>2 or the elasticity values that reflect the real market situation, in line with the results of [28]. Finally, our research can be more widely applied to several challengeable volatility managing problems or risk management strategies, including the problem for the timer options under diverse underlying asset models.
Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
Acknowledgments
The work of S. Y. Choi was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2021R1F1A1046138) and the Gachon University research fund of 2023 (2023039990001). The research by J. H. Yoon was supported by the the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2022R1A5A1033624 and No. 2023R1A2C1006600).
Conflict of interest
The authors have no conflicts of interest to declare.
References
[1]
M. Courbage, V. I. Nekorkin, Map based models in neurodynamics, Int. J. Bifurcation Chaos, 20 (2010), 1631–1651. https://doi.org/10.1142/S0218127410026733 doi: 10.1142/S0218127410026733
[2]
B. Ibarz, J. M. Casado, M. A. Sanjuán, Map-based models in neuronal dynamics, Phys. Rep., 501 (2011), 1–74. https://doi.org/10.1016/j.physrep.2010.12.003 doi: 10.1016/j.physrep.2010.12.003
[3]
A. Holko, M. Mȩdrek, Z. Pastuszak, K. Phusavat, Epidemiological modeling with a population density map-based cellular automata simulation system, Expert Syst. Appl., 48 (2016), 1–8. https://doi.org/10.1016/j.eswa.2015.08.018 doi: 10.1016/j.eswa.2015.08.018
[4]
V. V. Tarasova, V. E. Tarasov, Logistic map with memory from economic model, Chaos Solitons Fractals, 95 (2017), 84–91. https://doi.org/10.1016/j.chaos.2016.12.012 doi: 10.1016/j.chaos.2016.12.012
[5]
K. Tanaka, T. Kinkyo, S. Hamori, Financial hazard map: Financial vulnerability predicted by a random forests classification model, Sustainability, 10 (2018), 1530. https://doi.org/10.3390/su10051530 doi: 10.3390/su10051530
[6]
Z. Hua, Y. Zhou, Exponential chaotic model for generating robust chaos, IEEE Trans. Syst. Man Cybern., 51 (2019), 3713–3724. https://doi.org/10.1109/TSMC.2019.2932616 doi: 10.1109/TSMC.2019.2932616
[7]
M. Ausloos, The Logistic Map and the Route to Chaos: From the Beginnings to Modern Applications, Springer Science & Business Media, 2006.
[8]
B. Cessac, B. Doyon, M. Quoy, M. Samuelides, Mean-field equations, bifurcation map and route to chaos in discrete time neural networks, Physica D, 74 (1994), 24–44. https://doi.org/10.1016/0167-2789(94)90024-8 doi: 10.1016/0167-2789(94)90024-8
[9]
Q. Xu, T. Liu, S. Ding, H. Bao, Z. Li, B. Chen, Extreme multistability and phase synchronization in a heterogeneous bi-neuron Rulkov network with memristive electromagnetic induction, Cognit. Neurodyn., 17 (2023), 755–766. https://doi.org/10.1007/s11571-022-09866-3 doi: 10.1007/s11571-022-09866-3
[10]
I. Bashkirtseva, L. Ryashko, J. M. Seoane, M. A. Sanjuán, Noise-induced complex dynamics and synchronization in the map-based Chialvo neuron model, Commun. Nonlinear Sci. Numer. Simul., 116 (2023), 106867. https://doi.org/10.1016/j.cnsns.2022.106867 doi: 10.1016/j.cnsns.2022.106867
[11]
J. Sausedo-Solorio, A. Pisarchik, Synchronization of map-based neurons with memory and synaptic delay, Phys. Lett. A, 378 (2014), 2108–2112. https://doi.org/10.1016/j.physleta.2014.05.003 doi: 10.1016/j.physleta.2014.05.003
[12]
I. Franović, V. Miljković, Power law behavior related to mutual synchronization of chemically coupled map neurons, Eur. Phys. J. B, 76 (2010), 613–624. https://doi.org/10.1140/epjb/e2010-00205-4 doi: 10.1140/epjb/e2010-00205-4
[13]
Q. Xu, T. Liu, C. T. Feng, H. Bao, H. G. Wu, B. C. Bao, Continuous non-autonomous memristive Rulkov model with extreme multistability, Chin. Phys. B, 30 (2021), 128702. https://doi.org/10.1088/1674-1056/ac2f30 doi: 10.1088/1674-1056/ac2f30
[14]
S. He, K. Rajagopal, A. Karthikeyan, A. Srinivasan, A discrete Huber-Braun neuron model: From nodal properties to network performance, Cognit. Neurodyn., 17 (2023), 301–310. https://doi.org/10.1007/s11571-022-09806-1 doi: 10.1007/s11571-022-09806-1
[15]
F. Wang, H. Cao, Mode locking and quasiperiodicity in a discrete-time Chialvo neuron model, Commun. Nonlinear Sci. Numer. Simul., 56 (2018), 481–489. https://doi.org/10.1016/j.cnsns.2017.08.027 doi: 10.1016/j.cnsns.2017.08.027
X. Zhang, C. Li, E. Dong, Y. Zhao, Z. Liu, A conservative memristive system with amplitude control and offset boosting, Int. J. Bifurcation Chaos, 32 (2022), 2250057. https://doi.org/10.1142/S0218127422500572 doi: 10.1142/S0218127422500572
[18]
Y. Jiang, C. Li, C. Zhang, Y. Zhao, H. Zang, A double-memristor hyperchaotic oscillator with complete amplitude control, IEEE Trans. Circuits Syst. I, 68 (2021), 4935–4944. https://doi.org/10.1109/TCSI.2021.3121499 doi: 10.1109/TCSI.2021.3121499
[19]
H. Jeong, L. Shi, Memristor devices for neural networks, J. Phys. D, 52 (2018), 023003. https://doi.org/10.1088/1361-6463/aae223 doi: 10.1088/1361-6463/aae223
[20]
D. B. Strukov, G. S. Snider, D. R. Stewart, R. S. Williams, The missing memristor found, Nature, 453 (2008), 80–83. https://doi.org/10.1038/nature06932 doi: 10.1038/nature06932
[21]
H. Wang, C. L. Li, S. Banerjee, S. B. He, Novel memristor and memristor-based applications, Eur. Phys. J. Spec. Top., 231 (2022), 2973–2977. https://doi.org/10.1140/epjs/s11734-022-00697-1 doi: 10.1140/epjs/s11734-022-00697-1
[22]
M. Ge, Y. Jia, Y. Xu, L. Yang, Mode transition in electrical activities of neuron driven by high and low frequency stimulus in the presence of electromagnetic induction and radiation, Nonlinear Dyn., 91 (2018), 515–523. https://doi.org/10.1007/s11071-017-3886-2 doi: 10.1007/s11071-017-3886-2
[23]
S. P. Adhikari, C. Yang, H. Kim, L. O. Chua, Memristor bridge synapse-based neural network and its learning, IEEE Trans. Neural Networks Learn. Syst., 23 (2012), 1426–1435. https://doi.org/10.1109/TNNLS.2012.2204770 doi: 10.1109/TNNLS.2012.2204770
[24]
J. C. Magee, C. Grienberger, Synaptic plasticity forms and functions, Ann. Rev. Neurosci., 43 (2020), 95–117. https://doi.org/10.1146/annurev-neuro-090919-022842 doi: 10.1146/annurev-neuro-090919-022842
[25]
C. Chen, F. Min, Y. Zhang, B. Bao, Memristive electromagnetic induction effects on Hopfield neural network, Nonlinear Dyn., 106 (2021), 2559–2576. https://doi.org/10.1007/s11071-021-06910-5 doi: 10.1007/s11071-021-06910-5
[26]
S. Hu, Y. Liu, Z. Liu, T. Chen, J. Wang, Q. Yu, et al., Associative memory realized by a reconfigurable memristive Hopfield neural network, Nat. Commun., 6 (2015), 7522. https://doi.org/10.1038/ncomms8522 doi: 10.1038/ncomms8522
[27]
K. Rajagopal, A. Karthikeyan, S. Jafari, F. Parastesh, C. Volos, I. Hussain, Wave propagation and spiral wave formation in a Hindmarsh-Rose neuron model with fractional-order threshold memristor synapse, Int. J. Mod. Phys. B, 34 (2020), 2050157. https://doi.org/10.1142/S021797922050157X doi: 10.1142/S021797922050157X
[28]
Q. Xu, Z. Ju, S. Ding, C. Feng, M. Chen, B. Bao, Electromagnetic induction effects on electrical activity within a memristive Wilson neuron model, Cognit. Neurodyn., 16 (2022), 1221–1231. https://doi.org/10.1007/s11571-021-09764-0 doi: 10.1007/s11571-021-09764-0
[29]
S. Qiao, C. Gao, Complex dynamics of a non-smooth temperature-sensitive memristive Wilson neuron model, Commun. Nonlinear Sci. Numer. Simul., 125 (2023), 107410. https://doi.org/10.1016/j.cnsns.2023.107410 doi: 10.1016/j.cnsns.2023.107410
[30]
C. Li, Y. Yang, X. Yang, Y. Lu, Application of discrete memristors in logistic map and Hindmarsh-Rose neuron, Eur. Phys. J. Spec. Top., 231 (2022), 3209–3224. https://doi.org/10.1140/epjs/s11734-022-00645-z doi: 10.1140/epjs/s11734-022-00645-z
[31]
Q. Lai, C. Lai, H. Zhang, C. Li, Hidden coexisting hyperchaos of new memristive neuron model and its application in image encryption, Chaos Solitons Fractals, 158 (2022), 112017. https://doi.org/10.1016/j.chaos.2022.112017 doi: 10.1016/j.chaos.2022.112017
[32]
C. Gao, S. Qiao, X. An, Global multistability and mechanisms of a memristive autapse-based Filippov Hindmash-Rose neuron model, Chaos Solitons Fractals, 160 (2022), 112281. https://doi.org/10.1016/j.chaos.2022.112281 doi: 10.1016/j.chaos.2022.112281
[33]
S. Qiao, C. Gao, X. An, Hidden dynamics and control of a Filippov memristive hybrid neuron model, Nonlinear Dyn., 111 (2023), 10529–10557. https://doi.org/10.1016/j.chaos.2022.112017 doi: 10.1016/j.chaos.2022.112017
[34]
A. Arenas, A. Díaz-Guilera, J. Kurths, Y. Moreno, C. Zhou, Synchronization in complex networks, Phys. Rep., 469 (2008), 93–153. https://doi.org/10.1016/j.physrep.2008.09.002 doi: 10.1016/j.physrep.2008.09.002
[35]
Z. Wang, R. Ramamoorthy, X. Xi, H. Namazi, Synchronization of the neurons coupled with sequential developing electrical and chemical synapses, Math. Biosci. Eng., 19 (2022), 1877–1890. https://doi.org/10.3934/mbe.2022088 doi: 10.3934/mbe.2022088
[36]
S. N. Chowdhury, S. Majhi, M. Ozer, D. Ghosh, M. Perc, Synchronization to extreme events in moving agents, New J. Phys., 21 (2019), 073048. https://doi.org/10.1088/1367-2630/ab2a1f doi: 10.1088/1367-2630/ab2a1f
[37]
G. Vivekanandhan, H. Natiq, Y. Merrikhi, K. Rajagopal, S. Jafari, Dynamical analysis and synchronization of a new memristive Chialvo neuron model, Electronics, 12 (2023), 545. https://doi.org/10.3390/electronics12030545 doi: 10.3390/electronics12030545
[38]
S. He, Complexity and chimera states in a ring-coupled fractional-order memristor neural network, Front. Appl. Math. Stat., 6 (2020), 24. https://doi.org/10.3389/fams.2020.00024 doi: 10.3389/fams.2020.00024
[39]
Z. Wang, H. Tian, O. Krejcar, H. Namazi, Synchronization in a network of map-based neurons with memristive synapse, Eur. Phys. J. Spec. Top., 231 (2022), 4057–4064. https://doi.org/10.1140/epjs/s11734-022-00691-7 doi: 10.1140/epjs/s11734-022-00691-7
[40]
F. Parastesh, S. Jafari, H. Azarnoush, Z. Shahriari, Z. Wang, S. Boccaletti, et al., Chimeras, Phys. Rep., 898 (2021), 1–114. https://doi.org/10.1016/j.physrep.2020.10.003 doi: 10.1016/j.physrep.2020.10.003
[41]
S. Majhi, B. K. Bera, D. Ghosh, M. Perc, Chimera states in neuronal networks: A review, Phys. Life Rev., 28 (2019), 100–121. https://doi.org/10.1016/j.plrev.2018.09.003 doi: 10.1016/j.plrev.2018.09.003
[42]
I. Franović, K. Todorović, N. Vasović, N. Burić, Cluster synchronization of spiking induced by noise and interaction delays in homogenous neuronal ensembles, Chaos, 22 (2012). https://doi.org/10.1063/1.4753919
[43]
I. Franović, K. Todorović, N. Vasović, N. Burić, Spontaneous formation of synchronization clusters in homogenous neuronal ensembles induced by noise and interaction delays, Phys. Rev. Lett., 108 (2012), 094101. https://doi.org/10.1103/PhysRevLett.108.094101 doi: 10.1103/PhysRevLett.108.094101
[44]
M. Mehrabbeik, F. Parastesh, J. Ramadoss, K. Rajagopal, H. Namazi, S. Jafari, Synchronization and chimera states in the network of electrochemically coupled memristive Rulkov neuron maps, Math. Biosci. Eng., 18 (2021), 9394–9409. https://doi.org/10.3934/mbe.2021462 doi: 10.3934/mbe.2021462
[45]
K. Li, B. Bao, J. Ma, M. Chen, H. Bao, Synchronization transitions in a discrete memristor-coupled bi-neuron model, Chaos Solitons Fractals, 165 (2022), 112861. https://doi.org/10.1016/j.chaos.2022.112861 doi: 10.1016/j.chaos.2022.112861
[46]
B. Ramakrishnan, M. Mehrabbeik, F. Parastesh, K. Rajagopal, S. Jafari, A new memristive neuron map model and its network's dynamics under electrochemical coupling, Electronics, 11 (2022), 153. https://doi.org/10.3390/electronics11010153
[47]
S. Wang, Z. Wei, Synchronization of coupled memristive Hindmarsh–Rose maps under different coupling conditions, AEU-Int. J. Electron. Commun., 161 (2023), 154561. https://doi.org/10.1016/j.aeue.2023.154561 doi: 10.1016/j.aeue.2023.154561
[48]
G. Baghdadi, S. Jafari, J. C. Sprott, F. Towhidkhah, M. H. Golpayegani, A chaotic model of sustaining attention problem in attention deficit disorder, Commun. Nonlinear Sci. Numer. Simul., 20 (2015), 174–185. https://doi.org/10.1016/j.cnsns.2014.05.015 doi: 10.1016/j.cnsns.2014.05.015
[49]
B. Bao, H. Qian, Q. Xu, M. Chen, J. Wang, Y. Yu, Coexisting behaviors of asymmetric attractors in hyperbolic-type memristor based Hopfield neural network, Front. Comput. Neurosci., 11 (2017), 81. https://doi.org/10.3389/fncom.2017.00081 doi: 10.3389/fncom.2017.00081
[50]
D. Ding, J. Luo, X. Shan, Y. Hu, Z. Yang, L. Ding, Coexisting behaviors of a fraction-order novel hyperbolic-type memristor Hopfield neuron network based on three neurons, Int. J. Mod. Phys. B, 34 (2020), 2050302. https://doi.org/10.1142/S0217979220503026 doi: 10.1142/S0217979220503026
[51]
R. Li, Z. Wang, E. Dong, A new locally active memristive synapse-coupled neuron model, Nonlinear Dyn., 104 (2021), 4459–4475. https://doi.org/10.1007/s11071-021-06574-1 doi: 10.1007/s11071-021-06574-1
[52]
K. Li, H. Bao, H. Li, J. Ma, Z. Hua, B. Bao, Memristive Rulkov neuron model with magnetic induction effects, IEEE Trans. Ind. Inf., 18 (2021), 1726–1736. https://doi.org/10.1109/TⅡ.2021.3086819 doi: 10.1109/TⅡ.2021.3086819
[53]
S. Garai, S. Karmakar, S. Jafari, N. Pal, Coexistence of triple, quadruple attractors and Wada basin boundaries in a predator-prey model with additional food for predators, Commun. Nonlinear Sci. Numer. Simul., 121 (2023), 107208. https://doi.org/10.1016/j.cnsns.2023.107208 doi: 10.1016/j.cnsns.2023.107208
[54]
Y. Zhou, J. Gao, K. D. White, I. Merk, K. Yao, Perceptual dominance time distributions in multistable visual perception, Biol. Cybern., 90 (2004), 256–263. https://doi.org/10.1007/s00422-004-0472-8 doi: 10.1007/s00422-004-0472-8
[55]
D. Durstewitz, G. Deco, Computational significance of transient dynamics in cortical networks, Eur. J. Neurosci., 27 (2008), 217–227. https://doi.org/10.1111/j.1460-9568.2007.05976.x doi: 10.1111/j.1460-9568.2007.05976.x
[56]
T. Ionescu, Exploring the nature of cognitive flexibility, New Ideas Psychol., 30 (2012), 190–200. https://doi.org/10.1016/j.newideapsych.2011.11.001 doi: 10.1016/j.newideapsych.2011.11.001
[57]
Y. Xu, Y. Jia, J. Ma, A. Alsaedi, B. Ahmad, Synchronization between neurons coupled by memristor, Chaos Solitons Fractals, 104 (2017), 435–442. https://doi.org/10.1016/j.chaos.2017.09.002 doi: 10.1016/j.chaos.2017.09.002
[58]
P. Zhou, Y. Xu, J. Ma, Dynamical and coherence resonance in a photoelectric neuron under autaptic regulation, Physica A, 620 (2023), 128746. https://doi.org/10.1016/j.physa.2023.128746 doi: 10.1016/j.physa.2023.128746
[59]
A. N. Pisarchik, A. E. Hramov, Coherence resonance in neural networks: Theory and experiments, Phys. Rep., 1000 (2023), 1–57. https://doi.org/10.1016/j.physrep.2022.11.004 doi: 10.1016/j.physrep.2022.11.004
Table 1.
For given θ=1.8,1.9,2.1 and θ=2.2, solutions ˜Pϵ and errors are given against ϵ, where "R.Error" is the relative percentage error. Note that the baseline parameters are given by X0=100,T=1.0,r=0.03,K=100,B=0.3,v=0.05,⟨Λ(y)ψ′(y)⟩=0.0141,⟨f(y)ψ′(y)⟩=0.0023,ρ=−0.1,m=0.5,u=1.0,f(y)=ey and ⟨f(y)2⟩=0.0272.
Figure 1. Bifurcation diagrams of the original neural map model (a), and its memristive version (b–d) for different magnetic induction effects as a function of A. (a) μ=0, (b) μ=0.1,(c)μ=0.2,(d)μ=0.3. The bifurcations are obtained using attractor-following method (forward: red color and backward: blue color) using the initial condition [1,0]
Figure 2. Bifurcation diagram of the memristive neural map model for A=8 as a function of μ. (a) The blue, red and cyan colors correspond to the initial conditions [−1,0] and [1,0], and [−5,4],respectively. (b) The forward continuation is used for the initial condition [1,0] (red) and [-1, 0] (blue)
Figure 3. Dynamics of the neural map model for (a) A=8 and (b) A=9 by adding the memristor with different strengths. (a1, b1) μ=0, (a2, b2) μ=0.02, (a3, b3) μ=0.05, (a4, b4) μ=0.1. The initial condition is [−1,0]
Figure 4. Basin of attraction of coexisting periodic (blue) and chaotic (red) dynamics in memristive neural map model for A=8. (a) μ=0.1, (b) μ=0.25
Figure 5. Coexisting dynamics of the memristive neural map model for A = 8. (a) μ=0.1, (b) μ=0.2, (c) μ=0.51. The left and right panels show the time series corresponding to the initial conditions [1,0] and[−1,0], respectively
Figure 6. (a–c) The synchronization error of coupled memristive neural maps as a function of μ and coupling strength for (a) electrical coupling (gm=gc=0), (b) magnetic coupling (ge=gc=0), (c) chemical coupling (ge=gm=0). The first and second rows correspond to two different network structures where in the first row P=10 and in the second row P=25. (d–f) The synchronization error (E) and the synchronization factor (R) as a function of coupling strength for μ=0.5. (d) electrical coupling, (e) magnetic coupling, (f) chemical coupling
Figure 7. Dynamical behavior of electrically coupled neural map models (a) spatiotemporal pattern and the time snapshot for μ=0.4 and ge=0.007, (b)the corresponding time series of part (a). (c) spatiotemporal pattern and the time snapshot for μ=0.4 and ge=0.009, (d)the corresponding time series of part (c)
Figure 8. Dynamical behavior of magnetically coupled neural map models (a) spatiotemporal pattern and the time snapshot for μ=0.1 and gm=0.03, (b)the corresponding time series of part (a). (c) spatiotemporal pattern and the time snapshot for μ=0.1 and gm=0.08, (d)the corresponding time series of part (c)
Figure 9. Dynamical behavior of magnetically coupled neural map models (a) spatiotemporal pattern and the time snapshot for μ=0.5 and gm=0.04, (b)the corresponding time series of part (a). (c) spatiotemporal pattern and the time snapshot for μ=0.5 and gm=0.07, (d)the corresponding time series of part (c)
Figure 10. Dynamical behavior of chemically coupled neural map models (a) spatiotemporal pattern and the time snapshot for μ=0.1 and gc=0.002, (b)the corresponding time series of part (a). (c) spatiotemporal pattern and the time snapshot for μ=0.1 and gc=0.005, (d)the corresponding time series of part (c)
Figure 11. Dynamical behavior of chemically coupled neural map models (a) spatiotemporal pattern and the time snapshot for μ=0.1 and gc=0.04, (b)the corresponding time series of part (a). (c) spatiotemporal pattern and the time snapshot for μ=0.4 and gc=0.03, (d)the corresponding time series of part (c)