Next Article in Journal
Risk Mutualization in Central Clearing: An Answer to the Cross-Guarantee Phenomenon from the Financial Stability Viewpoint
Next Article in Special Issue
A Novel Implementation of Siamese Type Neural Networks in Predicting Rare Fluctuations in Financial Time Series
Previous Article in Journal
Earnings Management, Related Party Transactions and Corporate Performance: The Moderating Role of Internal Control
Previous Article in Special Issue
Merton Investment Problems in Finance and Insurance for the Hawkes-Based Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transformations of Telegraph Processes and Their Financial Applications

by
Anatoliy A. Pogorui
1,
Anatoliy Swishchuk
2 and
Ramón M. Rodríguez-Dagnino
3,*
1
Department of Mathematics, Zhytomyr Ivan Franko State University, Valyka Berdychivska St. 40, 10008 Zhytomyr, Ukraine
2
Department of Mathematics & Statistics, University of Calgary, Calgary, AB T2N 1N4, Canada
3
School of Engineering and Sciences, Tecnologico de Monterrey, Av. Eugenio Garza Sada 2501 Sur. C.P., 64849 Monterrey, Mexico
*
Author to whom correspondence should be addressed.
Submission received: 2 July 2021 / Revised: 2 August 2021 / Accepted: 5 August 2021 / Published: 17 August 2021
(This article belongs to the Special Issue Stochastic Modelling in Financial Mathematics)

Abstract

:
In this paper, we consider non-linear transformations of classical telegraph process. The main results consist of deriving a general partial differential Equation (PDE) for the probability density (pdf) of the transformed telegraph process, and then presenting the limiting PDE under Kac’s conditions, which may be interpreted as the equation for a diffusion process on a circle. This general case includes, for example, classical cases, such as limiting diffusion and geometric Brownian motion under some specifications of non-linear transformations (i.e., linear, exponential, etc.). We also give three applications of non-linear transformed telegraph process in finance: (1) application of classical telegraph process in the case of balance, (2) application of classical telegraph process in the case of dis-balance, and (3) application of asymmetric telegraph process. For these three cases, we present European call and put option prices. The novelty of the paper consists of new results for non-linear transformed classical telegraph process, new models for stock prices based on transformed telegraph process, and new applications of these models to option pricing.

1. Introduction

In 1951, Goldstein and Kac (see Goldstein (1951); Kac (1974) and also Kac (1950)) proposed an interesting random motion model for the movement of a particle on the line (or one dimension). The particle was moving at a constant velocity v in any of the two directions and traveling a random distance drawn from an exponential probability distribution with parameter λ . Therefore, this is a random motion driven by a Poisson process with intensity λ . After one movement, the particle changes its direction of motion in the opposite direction under the same stochastic conditions. This particle motion can be modeled as a random motion governed by a switching Poisson process with alternating directions and having exponentially distributed holding times.
In an independent manner, Goldstein and Kac solved this problem and they found that the solution satisfies the one-dimensional telegraph equation, which has a similar mathematical form as the Heaviside telegraph equation appearing in deterministic problems of wave propagation in electrical transmission lines, namely:
2 f ( t , x ) t 2 + 2 λ f ( t , x ) t = v 2 2 f ( t , x ) x 2 .
In the case of the random motion probabilistic model this equation is also called the Goldstein–-Kac telegraph equation or classical telegraph equation.
This seminal work have been extended in many publications worldwide by introducing variations of this basic idea. For instance, we should mention applications in financial market theory of the one-dimensional jump telegraph process, which is a generalization of the telegraph process, Ratanov (2007, 2010), López and Ratanov (2012, 2014), and Ratanov and Melnikov (2008). Explicit formula for the distribution of the integrated telegraph process (or Kac’s process) first appeared in Janssen and Siebert (1981). Its proof was presented in Steutel (1985) (see also Orgingher (1990) for more details). Some connections between telegraph equation and heat equation may be found in Janssen (1990) along with asymptotic properties of integrated telegraph process. Distributions of the integrated telegraph process were obtained in Di Masi et al. (1994) in both symmetric (intensities of transitions are same) and asymmetric cases (intensities of transitions are different). In the hydrodynamic limit, this process approximates the diffusion process on the line. Some probabilistic analysis of the telegrapher’s process with drift by means of relativistic transformations were considered in Beghin et al. (2001). Variety of transformations of telegraph process and its association with many areas were studied by Orsingher; see Orsingher (1985); Orsingher and Beghin (2009); Orsingher and De Gregorio (2007); Orsingher and Ratanov (2002); Orsingher and Somella (2004); Orgingher (1990). They include hyperbolic equations, fractional diffusion equations, random flights, planar and cyclic random motions, among others. Some recent works consider a telegraph equation with time-dependent coefficients Angelani and Garra (2019), Markov-modulated Lévy processes with two different regimes of restarting Ratanov (2020), some generalizations of the classical Black–Scholes models in finance Stoynov (2019), jump-telegraph process with exponentially distributed interarrival times Di Crescenzo and Meoli (2018), and a model to describe the vertical motions in the Campi Flegrei volcanic region consisting of a Brownian motion process driven by a generalized telegraph process Travaglino et al. (2018).
The application of the telegraph process for option prices was studied in Ratanov (2007); Ratanov and Melnikov (2008); Ratanov (2010); Kolesnik and Ratanov (2013). Some applications of classical telegraph process in finance were considered in Pogorui et al. (2021b).
We note that the classical telegraph process is the simplest case of one-dimensional random evolutions (REs). A good introduction to RE may be found in Pinsky (1991). Random evolutions driven by the hyper-parabolic operators were considered in Kolesnik and Pinsky (2011). Many ideas, methods, classifications, applications, and examples of REs are presented in the handbook Swishchuk (1997). Many applications of REs in finance and insurance are considered in Swishchuk (2000).
In this paper, we consider transformations of classical telegraph process. We also give three applications of transform telegraph process in finance: (1) application of classical telegraph process ( h ( x ) = x ) in the case of balance; (2) application of classical telegraph process ( h ( x ) = x ) in the case of dis-balance; and (3) application of asymmetric telegraph process ( h ( x ) has a special form presented below) in finance. The novelty of the paper consists of new results for transformed classical telegraph process, new models for stock prices, and new applications of these models to option pricing. Function h ( x ) is used to generalize classical telegraph equation to obtain, e.g., asymmetric telegraph process, diffusion process on a circle, etc.
The main idea of application of telegraph process in finance is the following one: Instead of the geometric Brownian motion (GBM) we propose the following model for the price S t of a stock at time t :
S t ϵ = S 0 exp x ( t ) ,
where x ( t ) = 0 t v ( s ) d s , v ( t ) is a continuous-time Markov chain with state space ( v 1 , v 2 ) , and with λ i being the rates of the exponential waiting times, i = 1 , 2 .
To satisfy Kac’s conditions, we consider the scaled model for stock price:
S t ϵ = S 0 exp ε 1 0 t v s ε 2 d s ,
and then taking a limit when ϵ 0 . In this manner, process S t ϵ converges weakly to GBM S ^ t = S 0 e μ t + σ W ( t ) , with specific constants μ and σ , and { W ( t ) } is the Wiener process. We use the latter model to calculate European call and put options pricing.
For example, the modeling of cash flow in high-frequency and algorithmic trading can be fitted by this model. If transactions happen in a short time (i.e., every millisecond), then because of the scale 1 / ϵ 2 , v ( t / ϵ 2 ) is switched quickly between two states { v 1 , v 2 } according to Markov chain ξ ( t / ϵ 2 ) . Hence, S t ϵ is changing quickly as well. If we use a longer time interval t / ϵ 2 , instead of time t in milliseconds, we can fit this model for different purposes such as market making, liquidation, acquisition, etc., purposes. Thus, we can apply our asymptotic results considered above. Our modeling approach is an alternative to the well-known Black–Scholes model based on geometric Brownian motion. It is well-known that Brownian motion (Wiener process) has some mathematical difficulties that make it difficult to fit to real data. For instance, it has trajectories continuous everywhere but differentiable nowhere, it is fractal with Hausdorff dimension equals 1.5, it has zero length of free path, and infinite velocity at any point of time.
The paper is organized as follows: Section 2 considers transformations of classical telegraph process. Section 3 presents financial applications of transformed telegraph process. Section 4 concludes the paper.

2. Transformations of Classical Telegraph Process

Suppose that x t is the telegraph process, such that
d d t x t = 1 ξ t v ,
where v > 0 , ξ t is an alternating Markov process with phase space X = 0 , 1 and infinitesimal generator matrix (or intensity matrix)
Q = q 00 q 01 q 10 q 11 = λ λ λ λ .
Consider a differentiable function h : R R , such that there exists the inverse function h 1 . Let us introduce the following process,
z t = h x t .
Then,
d d t z t = d d t h x t = h x t 1 ξ t v = h h 1 z t 1 ξ t v .
The bivariate process z t , ξ t is a Markov process Korolyuk and Swishchuk (1995), Swishchuk (1997, 2000), and its infinitesimal generator L is given by
L φ x , n = h h 1 x 1 n v x φ x , n + Q φ x , n ,
where Q φ x , n = q n 0 φ x , 0 + q n 1 φ x , 1 , n = 1 , 2 .
Or in more details
L φ x , 0 = v h h 1 x x φ x , 0 λ φ x , 0 + λ φ x , 1
L φ x , 1 = v h h 1 x x φ x , 1 λ φ x , 1 + λ φ x , 0 .
Denote by f x , n , t , x R , t 0 the pdf of the process z t , ξ t in the case when ξ t = n , n 0 , 1 . Then,
t f x , n , t = L f x , n , t
or in more details
t f x , 0 , t = v h h 1 x x f x , 0 , t λ f x , 0 , t + λ f x , 1 , t
t f x , 1 , t = v h h 1 x x f x , 1 , t λ f x , 1 , t + λ f x , 0 , t .
Equivalently, these equations can be expressed in matrix form as follows:
t v h h 1 x x + λ λ λ t + v h h 1 x x + λ f x , 0 , t f x , 1 , t = 0 0 .
Let us define the following notation,
f x , t : = f x , 0 , t + f x , 1 , t .
It is easily seen that f x , t d x is the pdf of z t at x, i.e., f x , t d x = P z t x , x + d x .
Then, f x , t satisfies the following equation,
det t v h h 1 x x + λ λ λ t + v h h 1 x x + λ f x , t = 0 ,
or equivalently,
2 t 2 v 2 h h 1 x 2 2 x 2 + 2 λ t f x , t = 0 .
It is well-known under Kac’s conditions Kac (1950) that the telegraph process x ( t ) converges weakly to Wiener process W ( t ) and, hence, h ( x ( t ) ) converges weakly to h ( W ( t ) ) . Passing in the last equation to the Kac’s limit Kac (1950), i.e., when λ + and ν + in such a way that ν 2 / λ σ 2 > 0 , we obtain
t f h ( w ) ( x , t ) = σ 2 2 ( h ( h 1 ( x ) ) ) 2 2 x 2 f h ( w ) ( x , t ) ,
where f h ( w ) ( x , t ) is the pdf of the process h ( W ( t ) ) .
On the other hand, the pdf f x , t of the telegraph process with the initial density distribution g x and equally probable velocities v and v satisfies the telegraph equation
2 t 2 v 2 2 x 2 + 2 λ t f x , t = 0 , t 0 , x R ,
and initial conditions
f x , 0 = g x , t f x , t t = 0 = 0 .
It is well-known that the unique solution of Cauchy problem (5) is given by the following formula
f x , t = e λ t 2 g v t x + g v t + x + e λ t 2 v v t v t g y + x p y , t d y ,
where p y , t = λ I 0 λ v v t 2 y 2 + t I 0 λ v v t 2 y 2 .
In the particular case, when the telegraph process starts from x = 0 , with equally probable velocities v and v , its pdf f x , t satisfies the telegraph equation
2 t 2 v 2 2 x 2 + 2 λ t f x , t = 0
with initial conditions f x , 0 = δ x , t f x , t t = 0 = 0 , and is given by
f x , t = e λ t 2 δ v t x + δ v t + x
+ e λ t 2 v λ I 0 λ v v t 2 x 2 + t I 0 λ v v t 2 x 2 ,
where v t 2 x 2 .
It is also well-known that if a random variable X has the pdf f X x and Y = h X and there exists the inverse function h 1 , then the pdf f Y y of Y is as follows:
f Y y = f X h 1 y d h 1 d x y .
Therefore, the solution f h x , t of Equation () with initial conditions f h x , 0 = g h 1 x , t f h x , t t = 0 = 0 is given by the following formula:
f h x , t = e λ t 2 d h 1 d x x
× { g v t h 1 x + g v t + h 1 x
+ 1 v v t v t g y + h 1 x λ I 0 λ v v t 2 y 2 + t I 0 λ v v t 2 y 2 d y ,
where v t 2 h 1 x 2 .
In particular, the solution f h x , t of Equation () with initial conditions f h x , 0 = δ h 1 x , t f h x , t t = 0 = 0 is given by the following formula:
f h x , t = e λ t 2 d h 1 d x x
× δ v t h 1 x + δ v t + h 1 x + 1 v λ I 0 λ v v t 2 h 1 x 2
+ t I 0 λ v v t 2 h 1 x 2 ,
where v t 2 h 1 x 2 .
Now, let us consider a case where h is a differentiable mapping h : R R or ( C ) and the inverse function does not necessarily exist for it.
In particular, in the case where h x = e i x , we have
d d t z t = d d t e i x t = i e i x t d d t x t = i z t 1 ξ t v .
Hence,
d d t cos x ( t ) = v 1 ξ t + 1 sin x ( t ) d d t sin x ( t ) = v 1 ξ t cos x ( t ) .
By assuming α ( t ) = cos x ( t ) and β ( t ) = sin x ( t ) , we have
d d t α ( t ) = v 1 ξ ( t ) + 1 β ( t )
d d t β ( t ) = v 1 ξ ( t ) α ( t ) .
The process α ( t ) , β ( t ) , ξ ( t ) is Markov and its infinitesimal generator L is given by the formula Korolyuk and Swishchuk (1995)
L φ α , β , ξ t
= v 1 ξ t + 1 β α φ α , β , ξ t + v 1 ξ t α β φ u , v , ξ t + Q φ α , β , ξ t .
Denoting by f α , β , n , t , x R , t 0 the pdf of the process α t , β t , ξ t in the special case when ξ t = n , n 0 , 1 , we have
t f α , β , n , t = L f α , β , n , t
or in more details
t f α , β , 0 , t =
v β α f α , β , 0 , t + v α β f α , β , 0 , t λ f α , β , 0 , t + λ f α , β , 1 , t t f α , β , 1 , t
= v β α f α , β , 1 , t v α β f α , β , 2 , t λ f α , β , 1 , t + λ f α , β , 0 , t .
Passing to the polar coordinate system α = r cos φ , β = r sin φ , we have
α f = r f cos φ φ f sin φ r β f = r f sin φ φ f cos φ r .
Taking into account that f does not depend on r and substituting (9) into (8), we have
t f φ , 0 , t = v sin 2 φ φ f φ , 0 , t v cos 2 φ φ f φ , 0 , t + λ f φ , 1 , t
t f φ , 1 , t = v sin 2 φ φ f φ , 1 , t + v cos 2 φ φ f φ , 0 , t λ f φ , 1 , t + λ f φ , 0 , t .
Denote by
f φ , t : = f φ , 0 , t + f φ , 1 , t .
In much the same manner as we obtained (3), we have
det a 11 λ λ a 22 f φ , t = 0 ,
where
a 11 = t v sin 2 φ φ + v cos 2 φ φ + λ
a 22 = t + v sin 2 φ φ v cos 2 φ φ + λ .
It is easily verified that
det a 11 λ λ a 22 = 2 t 2 + 2 λ t cos 2 2 φ 2 φ 2 + sin 4 φ 2 φ .
Therefore, the pdf f φ , t satisfies the following equation:
2 t 2 + 2 λ t v 2 cos 2 2 φ 2 φ 2 + v 2 sin 4 φ 2 φ f φ , t = 0 ,
with initial conditions f φ , 0 = δ φ , φ f φ , t φ = 0 = 0 .
Passing to the Kac’s limit in (10) as λ and v , such that v 2 λ σ 2 > 0 we obtain
t f φ , t = σ 2 2 cos 2 2 φ 2 φ 2 f φ , t sin 4 φ 2 φ f φ , t .
By analogy with a diffusion process on a line, this equation can be interpreted as the equation for a diffusion process on a circle De Gregorio and Iafrate (2020).
Remark 1.
Let us consider the case where h ( x ) = e x , i.e., z ( t ) = e x ( t ) . Then we have h ( h 1 ( x ) ) = x . That is,
2 t 2 v 2 x 2 2 x 2 + 2 λ t f x , t = 0 , t 0 , x R .
The solution to the last equation with the initial conditions f h ( x , 0 ) = δ ( ln ( x ) ) = δ ( x 1 ) , and t f h ( x , t ) | t = 0 = 0 , is given by the following formula:
f x , t = e λ t 2 | x | δ v t ln ( x ) + δ v t + ln ( x )
+ 1 v λ I 0 λ v v t 2 ( ln x ) 2 + t I 0 λ v v t 2 ( ln x ) 2 ,
where v t 2 ( ln x ) 2 .
Remark 2.
We note that if λ and v , such that v 2 λ σ 2 > 0 , i.e., under the Kac’s condition Kac (1974), the process e x ( t ) weakly converges to the geometric Brownian motion e W ( t ) that is used for many authors for modeling of a stock price in the Black–Scholes formula.

3. Financial Applications of Transformed Telegraph Process

In this section, we give three applications of transformations of telegraph process in finance: (1) application of classical telegraph process ( h ( x ) = x ) in the case of balance; (2) application of classical telegraph process ( h ( x ) = x ) in the case of dis-balance; and (3) application of asymmetric telegraph process ( h ( x ) has a special form presented below) in finance.
We note that the asymmetric telegraph process is a telegraph process where the particle is allowed to move forward or backward with two different velocities, v 1 , v 2 , Beghin et al. (2001), De Gregorio and Iafrate (2020), and López and Ratanov (2014). Furthermore, two different velocity switching rates λ 1 and λ 2 , are allowed in this process. Thus, the underlying telegraph signal can be modeled as a continuous-time Markov chain { X ( t ) , t 0 } with state space ( v 1 , v 2 ) , where λ i is the rate of the exponential sojourn time when the telegraph signal is in state ( 1 ) i + 1 v i , i = 1 , 2 , and X ( 0 ) is uniformly distributed on ( v 1 , v 2 ) . Then, the asymmetric telegraph process is defined as x ( t ) = 0 t X ( s ) d s .
It is straightforward to recover the classical telegraph process as a particular case when v 1 = v 2 = v , and λ 1 = λ 2 = λ .
Thus, in our case, the balance condition for classical telegraph process means ( v 1 + v 2 ) / 2 = 0 , i.e., v 2 = v 1 , and this is the case because v 1 = v 2 = v ; thus, v + ( v ) = 0 . Dis-balance condition for classical telegraph process means that ( v 1 + v 2 ) / 2 0 , i.e., v 2 v 1 . Thus, we have different velocities v 1 , v 2 , that satisfy the dis-balance condition.
We note that in the case of asymmetric telegraph process the transformed telegraph process may be presented by the following function h ( x ) :
h ( x ( t ) ) = 0 t v 1 v 2 2 + v 1 + v 2 2 ( 1 X ( 0 ) = v 1 1 X ( 0 ) = v 2 ) ( 1 ) N ( s ) d s ,
where N ( t ) is a counting process which has intensities λ i and λ i + ( 1 ) i + 1 for switching times S 2 i 1 and S 2 i , respectively, i = 1 , 2 .

3.1. Application of Classical Telegraph Process in Finance: Balance Case

It is possible to derive an analog model of the Black–Scholes formula after an application of the asymptotic results presented in Pogorui and Rodríguez-Dagnino (2008, 2009); Pogorui et al. (2021a).
The classical symmetric telegraph process is defined as follows: a particle is allowed to move forward or backward with velocities v , v , in an alternate manner, and the process has a velocity switching rate λ . Hence, the underlying telegraph signal can be modeled as a continuous-time Markov chain { X ( t ) , t 0 } with state space ( v , v ) , where λ is the rate of the exponential sojourn time in the interval when the telegraph signal is in state ( 1 ) i + 1 v . Therefore, the classical symmetric telegraph process is defined as x ( t ) = 0 t X ( s ) d s .
The probability law of the asymmetric telegraph process has an absolutely continuous component f that satisfies the following hyperbolic Equation (see Beghin et al. (2001)):
2 f ( t , x ) t 2 = v 2 2 f ( t , x ) x 2 2 λ f ( t , x ) t .
By taking the limits for λ , v + , such that
v λ σ ,
one obtains the governing equation of a Brownian motion (Kac-type condition). it is also possible to show that the classical symmetric telegraph process converges in distribution to a Brownian motion, i.e.,
x ( t ) d σ W ( t ) ,
where W ( t ) is a Wiener process (standard Brownian motion) and
σ : = v λ .
Let us consider the following model for a stock price:
S t = S 0 e x ( t ) ,
where x ( t ) is a classical symmetric telegraph process. Under above-mentioned Kac’s conditions we can state that
S t = S 0 e x ( t ) d S ^ t : = S 0 e σ W ( t ) .
After applying Itô’s formula (Shreve (2004)), we found S ^ t satisfies the following stochastic differential Equation (SDE):
d S ^ t = b S ^ t d t + σ S ^ t d W ( t ) ,
where
b : = σ 2 2 .
Now, let us define the following process:
Z ( t ) : = exp b r σ W t 1 2 b r σ 2 t ,
where r > 0 is the interest rate, and W t = W ( t ) .
Then, it is not difficult to see that E P [ Z ( t ) ] = 1 , hence according to Novikov’s result Novikov (1980), this process is a positive martingale. Now, let us define the probability measure Q (recall E P [ Z ( t ) ] = 1 ) on a complete probability space ( Ω , F ) :
Q ( A ) : = E P [ Z ( T ) 1 A ] ,
where 1 A is the indicator operator of the set A F , and the process Z ( t ) is defined above.
In a similar manner, we can define the following process:
W ^ t : = W t b r σ t ,
where b is defined above.
We can find that the stochastic process W ^ t is a standard Wiener process after applying Girsanov’s theorem Shreve (2004), under the probability measure Q. After this fact, measure Q is called a risk-neutral or martingale measure. Then, our stock price S ^ t in (14), under the risk-neutral measure Q, satisfies the following SDE:
d S ^ t = r S ^ t d t + σ S ^ t d W ^ t , σ = v λ .
Therefore, we can write the equivalent Black-Scholes formula for European call option price C ( t ) for our model in (15):
C ( t ) = S 0 N ( d 1 ) K e r ( T t ) N ( d 2 ) ,
where
d 1 : = ln ( S 0 / K ) + ( r + v 2 / 2 λ ) ( T t ) ( v / λ ) T t ,
d 2 : = d 1 ( v / λ ) T t ,
and
N ( x ) = Φ ( x ) = 1 2 π x e u 2 2 d u
is the cdf (cumulative distribution function) of a standard normal random variable with zero mean and unit variance, K is a strike price, and T is the maturity.
Example 1
(European call option for limiting telegraph process in balance case). Let us suppose the following numerical values: S 0 = $ 50 , K = $ 50 , r = 0.01, v = 0.01, λ = 0.01, T = 1 . Then, applying Formulas (16) and (17), we can obtain the following European call option price at time t = 0 :
C ( 0 ) = 50 × 0.559617 50 × e 0.01× 1 × 0.519938= $ 2.2426
or C ( 0 ) = 224.26 cents.
Below we show the numerical values of time evolution of C ( t ) dependent on λ (upon fixing v), on v (upon fixing λ ), see Figure 1, and of C ( t ) as a function of v and λ after fixing t, see Figure 2.
Remark 3.
In Figure 3 we show the cost behavior according to Equation (16) as a function of σ. The Black–Scholes limit case occurs by letting σ = v / λ = 0.01/ 0.01 = 0.1 .
Now, we can say that on longer time interval the BS formula works better but on shorter time interval our formula produces a better performance. The same for volatility: If volatility is bigger than 0.1, then C ( t ) is bigger.

3.2. Application of Classical Telegraph Process in Finance: Dis-Balance Case

Now, we study the one-dimensional transport process in the case of dis-balance. We consider first the scaled telegraph process and its limiting case (Section 3.2.1). Then, we applied the limiting process to option pricing; the stock price in this case is modeled as a geometric limiting telegraph process (Section 3.2.2).

3.2.1. Asymptotic Results for Scaled Telegraph Process

Consider a Markov process ξ t , t 0 with two states 0 , 1 and the generator matrix
Q = λ 1 1 1 1 .
Let us introduce the following random evolution or transport process
x t = x + 0 t v ( s ) d s ,
where
v t = v 0 , when ξ t = 0 ; v 1 , when ξ t = 1 .
The generator A of the bi-variate process
ς t = x t , ξ t , t 0
is as follows
A φ x , 0 = v 0 x φ x , 0 + λ φ x , 1 λ φ x , 0
A φ x , 1 = v 1 x φ x , 1 + λ φ x , 0 λ φ x , 1 ,
where φ D ( A ) is the domain of A, and x R .
The generator A can be interpreted in the following equivalent manner: Denote by Z = R × 0 , 1 and
T t φ x , i = R φ ( z ) P ς t d z ς t d z ς 0 = ( x , i ) ς 0 = ( x , i ) , i 0 , 1 .
Then,
A φ x , i = lim t 0 + T t φ x , i φ x , i t .
Considering u i ( x , t ) = T t φ x , i , we have
t u i x , t = A T t φ x , i = A u i x , t .
Let us introduce the scaled telegraph process
x ε t = x 0 + 1 ε 0 t v s ε 2 d s ,
with velocities v i / ε , i 0 , 1 . It is easily seen that v t ε 2 is Markovian and its generator is of the form 1 ε 2 Q .
Thus, we have the system of Kolmogorov differential equations:
t u 0 ε ( x , t ) = v 0 ε x u 0 ε ( x , t ) + λ ε 2 u 1 ε ( x , t ) λ ε 2 u 0 ε ( x , t )
t u 1 ε x , t = v 1 ε x u 1 ε x , t + λ ε 2 u 0 ε x , t λ ε 2 u 1 ε x , t ,
or in matrix form
t u ε ( x , t ) = 1 ε V u ε ( x , t ) + 1 ε 2 Q u ε ( x , t )
where
u ε ( x , t ) = u 0 ε ( x , t ) u 1 ε ( x , t ) , V = v 0 0 0 v 1 ,
and
V = v 0 x 0 0 v 1 x , Q = λ 1 1 1 1 .
Now, consider the potential matrix of ξ t , Korolyuk and Swishchuk (1995), Swishchuk (1997, 2000):
R 0 = 0 Π e Q t d t ,
where e Q t = p i j t ; i , j 0 , 1 are transition probabilities, and
Π = 1 2 1 2 1 2 1 2
is the projector matrix on N Q = ker Q . It is easily verified that R 0 = Π Q + Π 1 .
The balance condition implies v 0 + v 1 2 = 0 , i.e., v 1 = v 0 . It is easily verified that the balance condition can be also written as Π V Π = 0 (see Korolyuk and Swishchuk (1995), Swishchuk (1997, 2000)).
We are interested in the following dis-balance case: v 0 = v + Δ 1 and v 1 = v Δ 2 , where Δ i = ε a i , i = 1 , 2 , ε > 0 . It is easy to see that the generator of ζ ε t = x ε ( t ) , ξ t ε 2 has the following form:
A ε = 1 ε V + A + 1 ε 2 Q ,
where
A = a 1 0 0 a 2 .
Denote as u i ε ( x , t ) = T t ε φ x , i = Z φ z P ζ ε t d z / ζ ε 0 = x , i .
Then, much in the same way we obtain the following matrix equation
t u ε ( x , t ) = 1 ε V u ε ( x , t ) + A u ε ( x , t ) + 1 ε 2 Q u ε ( x , t ) .
By applying asymptotic expansion, Korolyuk and Swishchuk (1995), we have:
u ε x , t = u ( 0 ) x , t + n = 1 ε n u ( n ) x , t + v ( n ) x , t / ε 2 ,
where u ( n ) x , t , n = 0 , 1 , 2 , are the regular terms of the expansion whereas v ( n ) x , t / ε 2 , n = 1 , 2 , are the singular ones.
Then, by substituting (23) into (22), we obtain
Q u ( 0 ) x , t = 0 , Q u ( 1 ) x , t + V u ( 0 ) ( x , t ) = 0 ,
Q u ( 2 ) x , t + V u ( 1 ) ( x , t ) + A u ( 0 ) ( x , t ) t u ( 0 ) ( x , t ) = 0 ,
Q u ( k + 2 ) x , t + V u ( k + 1 ) ( x , t ) + A u ( k ) ( x , t ) t u ( k ) ( x , t ) = 0 ,
for k 0 .
Thus, u ( 0 ) x , t N Q , i.e., Π u ( 0 ) x , t = u ( 0 ) x , t .
From (24) it follows that
u ( 1 ) x , t = R 0 V u ( 0 ) ( x , t ) + c 1 ( t ) ,
where c 1 ( t ) N Q .
Much in the same way, we have
Q u ( 2 ) x , t = t u ( 0 ) ( x , t ) V u ( 1 ) ( x , t )
= t u ( 0 ) ( x , t ) V R 0 V u ( 0 ) ( x , t ) V c 1 ( t ) A u ( 0 ) ( x , t ) .
From the properties of Π it follows that
Π Q u ( 2 ) x , t = 0
= t u ( 0 ) ( x , t ) Π V R 0 V Π u ( 0 ) ( x , t ) Π A Π u ( 0 ) ( x , t ) .
Hence, the first term u ( 0 ) x , t of the Expansion (23) satisfies the diffusion Equation (27).
The matrix Equation (22) can be written as follows:
t v ε x a 1 x + λ ε 2 λ ε 2 λ ε 2 t + v ε x + a 2 x + λ ε 2 u ( ε ) x , t
= Ψ u ( ε ) x , t = 0 .
It is easily seen that the function f ε ( x , t ) = u 0 ε ( x , t ) + u 1 ε ( x , t ) is the solution of the following equation:
det Ψ f ε ( x , t ) = 0 .
Let us write (28) in more detail as follows:
2 λ t v 2 2 x 2 + λ a 2 a 1 x
+ ε 2 2 t 2 + a 2 a 1 2 x t v a 2 + a 1 ε 2 x 2 a 1 a 2 2 x 2
× f ε ( x , t ) = 0 .
Let us define the notation u ( 0 ) ( x , t ) : = u 0 ( x , t ) , u 1 ( x , t ) and f 0 ( x , t ) : = u 0 ( x , t ) + u 1 ( x , t ) . Since u ε x , t > u ( 0 ) ( x , t ) as ε > 0 Korolyuk and Swishchuk (1995), we have lim ε 0 f ε ( x , t ) = f 0 ( x , t ) .
From (29) it follows that
t v 2 2 λ 2 x 2 + a 2 a 1 2 x f 0 ( x , t ) = 0 .
Hence, if a 2 a 1 , then f 0 ( x , t ) satisfies the diffusion Equation (30) with drift a 1 a 2 2 and diffusion v 2 λ .
Remark 4.
It follows from (30) that scaled telegraph process x ϵ ( t ) in (19) weakly converges to a diffusion process with drift coefficient a 1 a 2 2 and diffusion coefficient v 2 λ .

3.2.2. Application in Finance: Black-Scholes Formula for Geometric Limiting Telegraph Process

The well-known geometric Brownian motion (GBM) is used to model a price S t of a stock (Shreve (2004)) at time t, such that
S t = S 0 exp μ t + σ W t ,
where μ and σ are the drift and volatility of the stock, respectively, and { W t } is the Wiener process.
After substituting μ = a 1 a 2 2 and σ = v λ in (31) the diffusion equation for the process μ t + σ W t can be obtained. Therefore, the Black–Scholes formula is obtained by considering an exponential Brownian motion for the share price S t .
As a consequence, we propose the following formula for the price S t of a stock at time t (see (19)):
S t ϵ = S 0 exp ε 1 0 t v s ε 2 d s ,
where v ( s ) is defined in (19) above. This formula represents an alternative to the formulation based on GBM.
This new formulation can be used to model cash flow in high-frequency and algorithmic trading. In many financial applications transactions happen in short periods of time (every millisecond), thus the stochastic process v ( t / ϵ 2 ) is switched quickly between two states { v 0 , v 1 } according to an underlying Markov chain ξ ( t / ϵ 2 ) , because of the scale 1 / ϵ 2 , and it means that S t ϵ is changing quickly as well. In some cases, we can assume longer periods of time t (instead of time in milliseconds), thus t / ϵ 2 can be used and it might be needed in applications, such as liquidation, acquisition, market making, etc. Therefore, we can apply our asymptotic results considered above for a better model fitting in these cases. For instance, below we show how to obtain an option pricing formula that is analogue to the Black–Scholes formula, for our model of a stock price. Of course, our modeling approach and results may be applied to other problems in mathematical finance, such as portfolio optimization, optimal control, etc.
Thus, by applying the results from the previous subsection we can state the following weak convergence
S t ϵ ϵ 0 S ^ t ,
in Skorokhod topology, where
S ^ t = S 0 exp a 1 a 2 2 t + v λ W t .
Now, applying Itô’s formula (Shreve (2004)), then S ^ t satisfies the stochastic differential equation:
d S ^ t = μ ^ S ^ t d t + v λ S ^ t d W t ,
where
μ ^ : = a 1 a 2 2 + v 2 2 λ .
Let us define the following stochastic process:
Z ( t ) : = exp μ ^ r v / λ W t 1 2 μ ^ r v / λ 2 t ,
where r > 0 is the interest rate.
Hence, E P [ Z ( t ) ] = 1 , and using Novikov’s result (Novikov (1980)), we conclude that this process is a positive martingale.
Let us define the new probability measure Q (recall E P [ Z ( t ) ] = 1 ) on a complete probability space ( Ω , F ) :
Q ( A ) : = E P [ Z ( T ) 1 A ] ,
where 1 A is an indicator operator of the set A F , and Z ( t ) is the stochastic process defined above.
We also define the stochastic process:
W ^ t : = W t μ ^ r v / λ t ,
where μ ^ is defined above.
After applying Girsanov’s theorem (Shreve (2004)), under the probability measure Q , the stochastic process we conclude that W ^ t is a standard Wiener process. We call measure Q a risk-neutral or martingale measure. Then, under the risk-neutral measure Q our stock price S ^ t satisfies the following SDE:
d S ^ t = r S ^ t d t + v λ S ^ t d W ^ t .
Therefore, we can write the Black–Scholes formula for European put option price P ( t ) for our model:
P ( t ) = K e r ( T t ) N ( d 2 ) S 0 N ( d 1 )
where
d 1 : = ln ( S 0 / K ) + ( r + v 2 / 2 λ ) ( T t ) ( v / λ ) T t ,
d 2 : = d 1 ( v / 2 λ ) T t ,
and N ( x ) is the cdf of a standard normal random variable with zero mean and unit variance, K is a strike price, and T is the maturity.
We note that
d 1 : = ln ( K / S 0 ) ( r + v 2 / 2 λ ) ( T t ) ( v / λ ) T t ,
and
d 2 = ( v / 2 λ ) T t d 1 .
We also note that European call option price is:
C ( t ) = S 0 N ( d 1 ) K e r ( T t ) N ( d 2 ) ,
where d 1 , d 2 are defined in (33).
Example 2
(European Put Option for Limiting Telegraph Process). Suppose the numerical values S 0 = $ 50 , K = $ 50 , r = 0.01, v = 0.02, λ = 0.02, T = 1 . Therefore, applying Formulas (32) and (33), we obtain the following European put option price at time t = 0 :
P ( 0 ) = 50 × e 0.01 × 1 × 0.471814 50 × 0.44376854= $ 1.167542
or P ( 0 ) = 116.7542 cents.
Below we present some graphs of the time evolution of P ( t ) as a function of λ (upon fixing v), and as a function of v (upon fixing λ ); see Figure 4, and on v and λ after fixing t; see Figure 5.
Remark 5.
In Figure 6 we show the cost behavior according to Equation (32) as a function of σ. The Black–Scholes limit case occurs by letting σ = 0.14 .
Now, we can say that on longer time interval the BS formula works better but on shorter time interval our formula produces a better performance. The same for volatility: If volatility is bigger than 0.14, then P ( t ) is bigger.

3.3. Asymmetric Telegraph Process and Its Financial Application

The asymmetric telegraph process is a telegraph process where the particle is allowed to move forward or backward with two different velocities, v 1 , v 2 ; and it has been studied in Beghin et al. (2001) and López and Ratanov (2014). Furthermore, this stochastic process can have two different velocity switching rates λ 1 and λ 2 . Hence, the underlying telegraph stochastic signal is modeled as a continuous-time Markov process ( X ( t ) , t 0 ) with state space ( v 1 , v 2 ) , where λ i is the rate of the exponential sojourn time when the telegraph signal is in state ( 1 ) i + 1 v i , i = 1 , 2 , and X ( 0 ) is uniformly distributed on ( v 1 , v 2 ) . Hence, the asymmetric telegraph process is defined as x ( t ) = 0 t X ( s ) d s . It is clear that the classical telegraph process can be recovered as a particular case when v 1 = v 2 = c , λ 1 = λ 2 = λ . The probability law of the asymmetric telegraph process has an absolutely continuous component f and it satisfies the following hyperbolic Equation (see Beghin et al. (2001)):
2 f ( t , x ) t 2 = ( v 1 v 2 ) 2 f ( t , x ) x 2 + ( v 2 v 1 ) 2 f ( t , x ) t x ( λ 1 + λ 2 ) f ( t , x ) t + 1 2 [ ( v 2 v 1 ) ( λ 1 + λ 2 ) ( λ 2 λ 1 ) ( v 1 + v 2 ) ] f ( t , x ) t .
Under Kac-type conditions we can take the limit in the first equation, and we can obtain the governing equation of a Brownian motion with drift. Hence, after taking the limits for λ i , v i + in such a manner that
λ 1 λ 2 ν 2 > 0 , v i λ i σ i , i = 1 , 2 , λ 2 v 1 λ 1 v 2 λ 1 + λ 2 δ ,
it is not difficult to show that the marginal distributions of the asymmetric telegraph process converges to a drifted Brownian motion
x ( t ) d δ t + σ W ( t ) ,
where W ( t ) is a standard Brownian motion and
σ : = σ 1 σ 2 ( σ 1 2 + σ 2 2 ) / 2 .
Remark 6.
We note that for the symmetric case, v 1 = v 2 = c , λ 1 = λ 2 = λ , the symmetric telegraph process x ( t ) under Kac’s conditions converges to standard Wiener process with volatility σ = c / λ :
x ( t ) d ( c / λ ) W ( t ) .
Let us consider the following model for a stock price:
S t = S 0 e x ( t ) ,
where x ( t ) is a telegraph process. Under above-mentioned Kac’s conditions we can state that
S t = S 0 e x ( t ) d S ^ t : = S 0 e δ t + σ W ( t ) .
Applying Itô’s formula (Shreve (2004)), the stochastic process S ^ t satisfies the following SDE:
d S ^ t = a S ^ t d t + σ S ^ t d W t ,
where
a : = δ + σ 2 2 .
Let us define the following stochastic process:
Z ( t ) : = exp a r σ W t 1 2 a r σ 2 t ,
where r > 0 is the interest rate.
Hence, E P [ Z ( t ) ] = 1 , and applying Novikov’s result (Novikov (1980)), we conclude that this process is a positive martingale. Now, let us define the new probability measure Q (recall E P [ Z ( t ) ] = 1 ) on a complete probability space ( Ω , F ) :
Q ( A ) : = E P [ Z ( T ) 1 A ] ,
where 1 A is an indicator operator of the set A F , and Z ( t ) is the stochastic process defined above.
We can also define the following process:
W ^ t : = W t a r σ t ,
where a is defined above.
By applying Girsanov’s theorem (Shreve (2004)), under the probability measure Q , then we conclude that the stochastic process W ^ t is a standard Wiener process. We call measure Q a risk-neutral or martingale measure. Thus, under the risk-neutral measure Q our stock price S ^ t satisfies the following SDE:
d S ^ t = r S ^ t d t + σ S ^ t d W ^ t .
Therefore, we can write the alternative form of Black–Scholes formula for European call option price C ( t ) for our model:
C ( t ) = S 0 N ( d 1 ) K e r ( T t ) N ( d 2 ) ,
where
d 1 : = ln ( S 0 / K ) + ( r + σ 2 / 2 ) ( T t ) σ T t ,
d 2 : = d 1 σ T t ,
σ : = σ 1 σ 2 ( σ 1 2 + σ 2 2 ) / 2 , σ i : = v i λ i ,
and N ( x ) is the cdf of a zero mean normal random variable with unit variance, K is a strike price, and T is the maturity.
Example 3
(European Call Option for Asymmetric Limiting telegraph Process). Suppose the numerical values S 0 = $ 50 , K = $ 50 , r = 0.01, v 1 = 0.01, v 2 = 0.02, λ 1 = 0.01, λ 2 = 0.02, T = 1 . Therefore, after applying formulas (34)-(36), we have σ 1 = 0.1 , σ 2 = 0.141421356, thus σ = 0.1155, and that European call option price at time t = 0 is:
C ( 0 ) = 50 × 0.557383 50 × e 0.01× 1 × 0.51151487= $ 2.547891
or C ( 0 ) = 254.7891 cents.
Now, below we present some graphs of the time evolution of C ( t ) dependent on λ 1 , λ 2 (upon fixing v i ), on v 1 , v 2 (upon fixing λ i ), see Figure 7, and on v 1 and λ 1 , see Figure 8.
Remark 7.
In Figure 9 we show the cost behavior according to Equation (34) as a function of σ. The Black–Scholes limit case occurs by letting σ = 0.115079291 (according to (36)).
Now, we can say that on longer time interval the BS formula works better but on shorter time interval our formula produces a better performance. The same for volatility: If volatility is bigger than 0.115079291, then C ( t ) is bigger.

4. Conclusions and Future Work

In this paper, we considered transformations of classical telegraph process. We also presented three applications of transform telegraph process in finance: (1) application of classical telegraph process in the case of balance, (2) application of classical telegraph process in the case of dis-balance, and (3) application of asymmetric telegraph process in finance. For these three cases, we presented European call and put option prices. The novelty of the paper consists of new results for transformed classical telegraph process, new models for stock prices and new applications of these models to option pricing.
As for the future work we could consider other problems in mathematical finance, such as portfolio optimization, optimal control, etc. Furthermore, we will calibrate λ , v and σ according to high-frequency and algorithmic trading (HFT) real data to see a better fit of our model. We will also try to apply our models of a stock price to optimization problems in HFT, such as optimal liquidation, acquisition, and market making. We will perform a comparative analysis of different models, including ours, in finance based on real data, as well.

Author Contributions

Conceptualization, A.A.P., A.S. and R.M.R.-D.; Formal analysis, A.A.P., A.S. and R.M.R.-D.; Investigation, A.A.P., A.S. and R.M.R.-D.; Software, R.M.R.-D.; Writing—original draft, A.A.P.; Writing—review & editing, A.A.P., A.S. and R.M.R.-D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The second author thanks NSERC for continuing support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Angelani, Luisa, and Roberto Garra. 2019. Probability distributions for the run-and-tumble models with variable speed and tumbling rate. Modern Stochastics: Theory and Applications 6: 3–12. [Google Scholar] [CrossRef]
  2. Beghin, Luisa, Luciano Nieddu, and Enzo Orsingher. 2001. Probabilistic analysis of the telegrapher’s process with drift by means of relativistic transformations. Special issue: Advances in applied stochastics. Journal of Applied Mathematics and Stochastic Analysis 14: 11–25. [Google Scholar] [CrossRef] [Green Version]
  3. De Gregorio, Alessandro, and Francesco Iafrate. 2020. Telegraph random evolutions on a circle. arXiv arXiv:2011.03025v2. [Google Scholar]
  4. Di Crescenzo, Antonio, and Alessandro Meoli. 2018. On a jump-telegraph process driven by an alternating fractional Poisson process. Journal of Applied Probability 55: 94–111. [Google Scholar] [CrossRef]
  5. Di Masi, Giovanni B., Yuri Kabanov, and J. Runggaldier. 1994. Mean-variance hedging of options on stocks with Markov volatilities. Theory of Probability and Its Applications 39: 211–22. [Google Scholar] [CrossRef]
  6. Goldstein, Sydney. 1951. On diffusion by discontinuous movements and on the telegraph equation. The Quarterly Journal of Mechanics and Applied Mathematics 4: 129–56. [Google Scholar] [CrossRef]
  7. Janssen, Arnold. 1990. The distance between the Kac process and the Wiener process with application to generalized telegraph equations. Journal of Theoretical Probability 3: 349–60. [Google Scholar] [CrossRef]
  8. Janssen, Arnold, and Eberhord Siebert. 1981. Convolution semigroups and generalized telegraph equations. Mathematische Zeitschrift 177: 519–32. [Google Scholar] [CrossRef]
  9. Kac, Mark. 1950. On some connections between probability theory and differential and integral equations. Paper presented at the Fifth Berkeley Symposium on Mathematical Statistics and Probability, Berkeley, LA, USA, July 31–August 12; Volume 12, pp. 189–215. [Google Scholar]
  10. Kac, Mark. 1974. A stochastic model related to the telegrapher’s equation. Rocky Mountain Journal of Mathematics 4: 497–509. [Google Scholar] [CrossRef]
  11. Kolesnik, Alexander D., and Mark A. Pinsky. 2011. Random evolutions are driven by the hyper-parabolic operators. Journal of Statistical Physics 142: 828–46. [Google Scholar] [CrossRef]
  12. Kolesnik, Alexander D., and Nikita Ratanov. 2013. Telegraph Processes and Option Pricing. Heidelberg: Springer. [Google Scholar]
  13. Korolyuk, Volodymyr Semonovych, and Anatoliy V. Swishchuk. 1995. Evolution of Systems in Random Media. Boca Raton: CRC Press. [Google Scholar]
  14. López, Oscar, and Nikita Ratanov. 2012. Option pricing under jump-telegraph model with random jumps. Journal of Applied Probability 49: 838–49. [Google Scholar] [CrossRef] [Green Version]
  15. López, Oscar, and Nikita Ratanov. 2014. On the asymmetric telegraph processes. Journal of Applied Probability 51: 569–89. [Google Scholar] [CrossRef] [Green Version]
  16. Novikov, Aleksandr A. 1980. On conditions for uniform integrability for continuous exponential martingales. Stochastic Differential Systems 25: 304–10. [Google Scholar]
  17. Orsingher, Enzo. 1985. Hyperbolic equations arising in random models. Stochastic Processes and their Applications 21: 93–106. [Google Scholar] [CrossRef] [Green Version]
  18. Orsingher, Enzo, and Luisa Beghin. 2009. Fractional diffusion equations and processes with randomly-varying time. The Annals of Probability 37: 206–49. [Google Scholar] [CrossRef]
  19. Orsingher, Enzo, and Alessandro De Gregorio. 2007. Random flights in higher spaces. Journal of Theoretical Probability 20: 769–806. [Google Scholar] [CrossRef]
  20. Orsingher, Enzo, and Ratanov Ratanov. 2002. Planar random motions with drift. Journal of Applied Mathematics and Stochastic Analysis 15: 205–21. [Google Scholar] [CrossRef] [Green Version]
  21. Orsingher, Enzo, and Andrea M. Somella. 2004. A cyclic random motion in R3 with four directions and finite velocity. Stochastics and Stochastics Reports 76: 113–33. [Google Scholar] [CrossRef]
  22. Orgingher, Enzo. 1990. Probability law, flow functions, maximum distribution of wave-governed random motions and their connections with Kirchoff’s laws. Stochastic Processes and Their Applications 34: 49–66. [Google Scholar] [CrossRef] [Green Version]
  23. Pinsky, Mark. 1991. Lectures on Random Evolutions. Singapore: World Scientific. [Google Scholar]
  24. Pogorui, Anatoliy A, and Ramón M. Rodríguez-Dagnino. 2008. Evolution process as an alternative to diffusion process and Black-Scholes formula. Paper presented at the 5th Conference in Actuarial Science & Finance on Samos, Samos, Greece, September 4–7. [Google Scholar]
  25. Pogorui, Anatoliy A., and Ramón M. Rodríguez-Dagnino. 2009. Evolution process as an alternative to diffusion process and Black-Scholes formula. Random Operators and Stochastic Equations 17: 61–68. [Google Scholar] [CrossRef] [Green Version]
  26. Pogorui, Anatoliy A., Anatoliy V. Swishchuk, and Ramón M. Rodríguez-Dagnino. 2021a. Random Motion in Markov and Semi-Markov Random Environment 1: Homogeneous and Inhomogeneous Random Motions. Volume 1. New York: ISTE Ltd. & Wiley. [Google Scholar]
  27. Pogorui, Anatoliy A., Anatoliy V. Swishchuk, and Ramón M. Rodríguez-Dagnino. 2021b. Random Motion in Markov and Semi-Markov Random Environment 2: High-dimensional Random Motions and Financial Applications. Volume 2. New York: ISTE Ltd. & Wiley. [Google Scholar]
  28. Ratanov, Nikita. 2007. A jump telegraph model for option pricing. Quantitative Finance 7: 575–83. [Google Scholar] [CrossRef]
  29. Ratanov, Nikita. 2010. Option pricing model based on a Markov-modulated diffusion with jumps. The Brazilian Journal of Probability and Statistics 24: 413–31. [Google Scholar] [CrossRef]
  30. Ratanov, Nikita. 2020. Kac–Lévy Processes. Journal of Theoretical Probability 33: 239–67. [Google Scholar] [CrossRef]
  31. Ratanov, Nikita, and Alexander Melnikov. 2008. On financial markets based on telegraph processes. Stochastics 80: 247–68. [Google Scholar] [CrossRef]
  32. Shreve, Steven. 2004. Stochastic Calculus for Finance: Continuous-Time Models. New York: Springer. [Google Scholar]
  33. Steutel, Fred W. 1985. Poisson processes and a Bessel function integral. SIAM Review 27: 73–77. [Google Scholar] [CrossRef] [Green Version]
  34. Stoynov, Pavel. 2019. Financial models beyond the classical Black-Scholes. AIP Conference Proceedings 2352: 030035. [Google Scholar]
  35. Swishchuk, Anatoliy V. 1997. Random Evolutions and their Applications. Dordrecht: Kluwer Academic Publishers. [Google Scholar]
  36. Swishchuk, Anatoliy V. 2000. Random Evolutions and their Applications: New Trends. Dordrecht: Kluwer Academic Publishers. [Google Scholar]
  37. Travaglino, Fabio, Antonio Di Crescenzo, Barbara Martinucci, and Roberto Scarpa. 2018. A New Model of Campi Flegrei Inflation and Deflation Episodes Based on Brownian Motion Driven by the Telegraph Process. Mathematical Geosciences 50: 961–75. [Google Scholar] [CrossRef]
Figure 1. Time evolution of European call option price C ( t ) as a function of v (a) and λ (b).
Figure 1. Time evolution of European call option price C ( t ) as a function of v (a) and λ (b).
Risks 09 00147 g001
Figure 2. Time evolution of European call option price C ( t ) as a function of v and λ for t = 0 (a) and t = 0.6 (b).
Figure 2. Time evolution of European call option price C ( t ) as a function of v and λ for t = 0 (a) and t = 0.6 (b).
Risks 09 00147 g002
Figure 3. Time evolution of European call option price C ( t ) as a function of σ .
Figure 3. Time evolution of European call option price C ( t ) as a function of σ .
Risks 09 00147 g003
Figure 4. Time evolution of European put option price P ( t ) and dependent on v (a) and λ (b).
Figure 4. Time evolution of European put option price P ( t ) and dependent on v (a) and λ (b).
Risks 09 00147 g004
Figure 5. Time evolution of European put option price P ( t ) and dependent on v and λ for t = 0 (a) and t = 0.6 (b).
Figure 5. Time evolution of European put option price P ( t ) and dependent on v and λ for t = 0 (a) and t = 0.6 (b).
Risks 09 00147 g005
Figure 6. Time evolution of European put option price P ( t ) as a function of σ .
Figure 6. Time evolution of European put option price P ( t ) as a function of σ .
Risks 09 00147 g006
Figure 7. Time evolution of European call option price C ( t ) by fix λ , t = 0 (a) and t = 0.6 (b) in the asymmetric case.
Figure 7. Time evolution of European call option price C ( t ) by fix λ , t = 0 (a) and t = 0.6 (b) in the asymmetric case.
Risks 09 00147 g007
Figure 8. Time evolution of European call option price C ( t ) for fix v, t = 0 (a) and t = 0.6 (b) in the asymmetric case.
Figure 8. Time evolution of European call option price C ( t ) for fix v, t = 0 (a) and t = 0.6 (b) in the asymmetric case.
Risks 09 00147 g008
Figure 9. Time evolution of European call option price C ( t ) as a function of σ .
Figure 9. Time evolution of European call option price C ( t ) as a function of σ .
Risks 09 00147 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pogorui, A.A.; Swishchuk, A.; Rodríguez-Dagnino, R.M. Transformations of Telegraph Processes and Their Financial Applications. Risks 2021, 9, 147. https://0-doi-org.brum.beds.ac.uk/10.3390/risks9080147

AMA Style

Pogorui AA, Swishchuk A, Rodríguez-Dagnino RM. Transformations of Telegraph Processes and Their Financial Applications. Risks. 2021; 9(8):147. https://0-doi-org.brum.beds.ac.uk/10.3390/risks9080147

Chicago/Turabian Style

Pogorui, Anatoliy A., Anatoliy Swishchuk, and Ramón M. Rodríguez-Dagnino. 2021. "Transformations of Telegraph Processes and Their Financial Applications" Risks 9, no. 8: 147. https://0-doi-org.brum.beds.ac.uk/10.3390/risks9080147

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop