arxiv: v1 [math.pr] 4 Feb 2016

Transcription

1 Mittag-Leffler Lévy Processes Arun Kumar an N. S. Upahye *Inian Statistical Institute, Chennai Center, Taramani, Chennai-636, Inia an **Department of Mathematics, Inian Institute of Technology Maras, Chennai 636, Inia arxiv: v1 [math.pr] 4 Feb 216 Abstract In this article, we introuce Mittag-Leffler Lévy process an provie two alternative representations of this process. First, in terms of Laplace transform of the marginal ensities an next as a suborinate stochastic process. Both these representations are useful in analyzing the properties of the process. Since integer orer moments for this process are not finite, we obtain fractional orer moments. Its ensity function an corresponing Lévy measure ensity is also obtaine. Further, heavy taile behavior of ensities an stochastic self-similarity of the process is emonstrate. Our results generalize an complement the results available on Mittag-Leffler istribution in several irections. Key wors: Mittag-Leffler istribution, suborinate stochastic processes, Lévy ensities. 1 Introuction In recent years there is an increase attention on Mittag-Leffler (ML) function as well as on ML probability istribution. ML istribution is a natural generalization of exponential istribution. ML waiting times are use in efining a fractional version of stanar Poisson process that is also calle fractional Poisson process (see e.g. Meerschaert et al. 21; Repin an Saichev, 2; Laskin, 23; Beghin an Orsingher, 29). Let X be a ML istribute ranom variable (rv) with parameters < α 1 an >. Then the Laplace transform (LT) of X is given by (see e.g. Cahoy et al. 21; Pillai, 199) Ee ux + u α. (1.1) Pillai (199) establishe the infinite ivisibility an geometric infinite ivisibility for X, along with other interesting properties. The ML ensity function can be written in terms of ML function as f X (x;, α) 1 E α, ( x α ), x >, where E α, (z) 1 t α t, (1.2) 2πi t α z C 1

2 with contour C starts an ens at an circles aroun the singularities an branch points of the integran (see e.g. Gorenflo et al. 214). The infinite series representation of E α, (z) is given by E α, (z) k z k Γ(αk + ). (1.3) Note that ML istributions are heavy taile an in the omain of attraction of stable laws (see e.g. Feller, 1971), which can be seen as follows. Let X 1, X 2,..., X n be ii ML istribute ranom variables with LT as in (1.1), then LT of the rescale rv T n 1/α n 1/α n i1 X i is given by which is the LT of a stable istribution. E(e utn ) (1 + u α /n) n e uα, as n (1.4) 2 Mittag-Leffler Lévy Process (MLLP) As mentione earlier, ML istribution is infinitely ivisible an hence we can efine a Lévy process corresponing to this istribution. We characterize the MLLP by efining the LT of its marginal ensity as follows. Definition 2.1 (MLLP). A stochastic process M α, (t) is a MLLP if it is a Lévy process an its marginal ensity has the LT given by Ee um α,(t) ( ) t. (2.1) + u α Recently, suborinate stochastic processes have attracte much attention ue to their applications in fractional partial ifferential equations, stochastic volatility moeling in finance, an for interesting probabilistic properties. Next, we erive suborinate stochastic representation of MLLP. Let G,β (t) be a gamma Lévy process such that G,β (t) gamma(, βt) with ensity f G,β (t)(y) Further, the LT of G,β (t) is given by βt Γ(βt) yβt 1 e y, y >, t >. Ee ug,β(t) ( ) βt. (2.2) + u Let S α (t) be a stable Lévy process with stability parameter α (, 1). Note that stable process has strictly increasing sample paths an its tail probability ecays polynomially such 2

3 that P(S α (t) > x) x α, as x for some constant >. Stable processes are self-similar such that S α (ct) c 1/α S α (t), c >, an the LT of S α (t) is given by Ee usα(t) e tuα. (2.3) Definition 2.2 (Alternative representation of MLLP). Let S α (t) be a stable Lévy process an G,β (t) be gamma Lévy process inepenent of S α (t). A three parameter MLLP can be efine as a suborinate stochastic process by using gamma Lévy process as a suborinator such that Note that M α,,β (t) : S α (G,β (t)). (2.4) Ee umα,,β(t) Ee usα(g,β(t)) ( ) E Ee usα(g,β(t)) G,β (t) E ( e ) uα G,β (t) (using (2.3)) ( ) βt (using (2.2)), + u α which is same as in (2.1) for β 1. For simplicity, throughout the article, we assume β 1 an enote M α, (t) : M α,,1 (t). We have, M α, (t) : S α (G,1 (t)) G,1 (t) 1/α S α (1), (2.5) using self-similar property of stable processes. Hence ML istribution is a mixture of gamma an stable istributions. It is well known that given a ranom variable U, uniformly istribute on [, 1], the transforme ranom variable ln U / is exponentially istribute with mean. Thus, (2.5) can be use in generating ML ranom numbers. Let U 1, U 2 an U 3 be inepenent an uniformly istribute over [,1]. Note that for t 1, G,1 (1) exp(). Thus by using Kanter (1975) result on stable ranom number generation an (2.5), we have M α, (1) sin(απu 1)[sin((1 π)u 2 )] 1/α 1 [sin(πu 1 )] 1/α ln U 2 1/α 1 ln U 3 1/α Similar to (1.4), we have convergence in istribution for MLLP process. 1/α. (2.6) Proposition /α M α, (t) lim S t t 1/α α (1). 3

4 Proof. We have 1/α M α, (t) t 1/α 1/α ( ) 1/α G,1 (t) S α (1) (using (2.5)) t 1/α (EG,1 (1)) 1/α S α (1) S α (1) as t, where the last convergence follows from the fact that for a suborinator T t, lim t T t /t a.s. ET 1 (see e.g. Bertion, 1996, p. 92). Further, EG,1 (1) 1/ an the result that if X n X an Y n c for some constant c, then X n Y n cx as n. 3 MLLP ensity We next obtain the infinite series representation of ensity function of MLLP. The main iea is to invert the LT given in (2.1). Proposition 3.1. MLLP has marginal ensities in following infinite series form f Mα, (t)(x) ( 1) k t+k Γ(t + k) x α(t+k) 1, x >. (3.1) Γ(t)Γ(k + 1) Γ(α(t + k)) k Proof. We have ( ) t f Mα, (t)(u) Ee umα,(t) + u α ( t u αt 1 + ) t u α k k k ( 1) k t+k Γ(t + k) Γ(t)Γ(k + 1) u α(t+k). ( ) t t+k u α(t+k) k Now, by inverting the LT in both sies an using L 1 (1/u a ) x a 1 /Γ(a), a >, where L 1 enotes the inverse LT, we have the esire result. Using the marginal ensity function of MLLP, we obtain the Lévy ensity corresponing to Lévy-Khintchine representation of this process. Proposition 3.2. The Lévy ensity in Lévy-Khintchine representation of the process is given by ν(x) α x E α,1( x α ), x >. (3.2) 4

5 Proof. For positive Lévy process (i.e. suborinator) with probability ensity f(x, t), Lévy measure ensity ν(x) is given by (see e.g. Barnorff-Nielsen, 21) Thus, using (3.1), we obtain ν(x) lim f(x, t). t t 1 ν(x) ( 1) k k xαk 1 Γ(αk) lim Γ(t + k) t tγ(t)γ(k + 1) k ( 1) k k xαk 1 kγ(αk) k α ( 1) k (x α ) k x Γ(αk + 1) α x E α,1( x α ) (using (1.3)). k Remark 3.1. For α 1, (3.2) reuces to ν(x) x 1 e x an which is the Lévy ensity corresponing to the gamma Lévy procees G,1 (t) (see e.g. Applebaum, 29, p. 55). Remark 3.2. Since, E α,1 ( x α ) 1 as x, we have ν(x) 1/x as x, where f(x) g(x) as x x means that lim x x f(x)/g(x) 1. Hence ν(x)x. Thus sample paths of MLLP are strictly increasing by an application of Theorem 21.3 of Sato (1999). Similar to gamma Lévy process the Lévy measure for MLLP is concentrate at origin an hence this process has an infinite arrival rate of jumps, most of which are small. For asymptotic behavior of ensity function, we use Tauberian theorem. First, we L(ct) recall that a function L(t) is slowly varying at some t, if for all fixe c >, lim 1. t t L(t) For the sake of convenience, we state Tauberian theorem here (see e.g. Bertoin, 1996, p. 1). Theorem 3.1 (Tauberian Theorem). Let L : (, ) (, ) be a slowly varying function at (respectively ) an let q. Then, for a function U : (, ) (, ), with corresponing LT Ũ the following are equivalent: (i) U(t) t q L(t)/Γ(1 + q), t (respectively t ). (ii) Ũ(u) u q 1 L(1/q), q (respectively q ). Proposition 3.3. The ensity function of MLLP has following asymptotic behaviors t f Mα, (t)(x) Γ(αt) xαt 1, as x, 1 ( + Γ( αt) xα ) t x αt 1, as x. 5

6 Proof. We have f Mα, (t)(u) Ee um α,(t) ( ) t + u α t u αt (1 + u α ) t t u αt as u. Now, by the application of Tauberian theorem, we have f Mα, (t)(x) Further, we have ) t ) t f Mα, (t)(u) (1 + uα (1 + [(1 u α + u α ) t ] t (1 + [1 u α 1 uα u α ) t [ ( (1 + u α t u αt ) ] t u α + u 2α ) t (1 + u α t u αt as u L(1/u)u ( αt 1) 1 t Γ(αt) xαt 1, as x. ) )] t (1 + uα Note that the function L(u) ( + u α ) t is slowly varying at. Hence again by the application of Tauberian Theorem, we have f Mα, (t)(x) 1 Γ( αt) ( + xα ) t x αt 1, as x. 4 MLLP fractional orer moments Using (2.4), we have E(M α, (t)) E (E(S α (G,β (t)) G,β (t))), since mean is infinite for α-stable istribution. Thus, we nee to fin out the fractional moments of the MLLP. Here we use LT efinition of MLLP to obtain the fractional orer moments. For a positive ranom variable X, the fractional orer moments from LT φ(u) can be obtaine as follows. u [φ(u)]up 1 u u [Ee ux ]u p 1 s [ ] E u [e ux ]u p 1 u (By Fubini Theorem) [ ] ( 1)E X e ux u p 1 u ( 1)Γ(p)EX 1 p. Thus (1 p)-th orer moment, where p is a real number such that p (, 1) is given by E(X 1 p ) ( 1) Γ(p) 6 u [φ(u)]up 1 u. (4.1)

7 Corollary 4.1. The fractional orer moments of MLLP are given by E(M α, (t) q t ) (1 q Γ(1 q) B q α, t + q ), < q < α. (4.2) α Proof. By taking p 1 q in (4.1) an using (2.1), we have E(M α, (t) q ) 1 [ ] t u q u Γ(1 q) u + u α αtt u α q 1 u Γ(1 q) ( + u α ) t+1 q α t Γ(1 q) t y q α (1 + y) t+1 y (by substituting uα y) y (1 q α ) 1 q α Γ(1 q) (1 + y) y (1 q α )+(t+ q α ) t (1 q α Γ(1 q) B q α, t + q ). α Remark 4.1. Alternatively representation (2.5) of MLLP can also be use in fining the fractional orer moments. Note that E (M α, (t)) q E[G,1 (t)] q/α E[S α (1)] q Γ( q + t) α q α Γ(t) Γ(1 q α ) Γ(1 q) t q α Γ(1 q) B (1 q α, t + q α Remark 4.2. for α 1 an t 1 in (4.2), we have EM 1, (1) q 1 ) q Γ(1 q) B(1 q, 1 + q) Γ(1+q) q, which is the q-th orer moment of exponential istribution with parameter. 5 Stochastic self-similarity We use the efinition of stochastic self-similarity introuce by Kozubowski et al. (26). For completeness purpose we recall their efinition. Definition 5.1 (Stochastic self-similarity). Let X(t) be a stochastic process. Let T c (t), c (1, ) be a family of stochastic processes inepenent of X(t) such that T c () a.s. with non-ecreasing sample paths an ET c (t) ct. The process X(t) is stochastically self-similar with inex η > with respect to the family T c (t) if {X(T c (t)), t } {c η X(t), t }, c (1, ), (5.1) where enotes the equivalence of finite imensional istributions. 7

8 Let N p (t) be a negative binomial Lévy process with rift efine as N p (t) : t + NB p (t), where NB p (t) is a Lévy process such that ( ) t + j 1 P(NB p (t) j) p t (1 p) j, j, 1, t >. (5.2) j We have the following result for the stochastic self-similarity of the MLLP. Proposition 5.1. MLLP is stochastic self-similar with inex 1/α with respect to the process N 1/c (t). That is {M α, (N 1/c (t)), t } {c 1/α M α, (t), t }. (5.3) Proof. Since M α, (N 1/c (t)) an c 1/α M α, (t) are both Lévy processes to prove (5.3), it is sufficient to show the equivalence of marginal istributions at t 1. Thus, we have E ( e um α,(n 1/c (1)) ) E ( e usα(g,1(n 1/c (1))) ) E ( E ( e usα(g,1(n 1/c (1))) G,1 (N 1/c (1)) )) E ( e uα G,1 (N 1/c (1)) ) E ( E ( e uα G,1 (N 1/c (1)) N 1/c (1) )) ( ) N1/c (1) ( ) 1+NB1/c (1) E E + u α + u α ( ) 1+j ( ) j (using (5.2)) + u α c c j + cu α, which is equal to E(e c1/α um α, (1) ) an hence the result. 6 Simulation The sample paths of MLLP can be simulate by suborinating a iscretize gamma Lévy process with stable process. We first simulate a iscretize gamma process on equally space intervals. Then we simulate the increments of MLLP process Y k M α, (t k ) M α, (t k 1 ) which conitionally on the values of G,1 (t k ) an G,1 (t k 1 ) is a stationary stable process. Note that Y k M α, (t k ) M α, (t k 1 ) S α (G,1 (t k )) S α (G,1 (t k 1 )) S α (G,1 (t k ) G,1 (t k 1 )) S α (G,1 (t k t k 1 )) (G,1 (t k t k 1 )) 1/α S α (1). 8

9 Using the above ieas, we have following algorithm to simulate the sample paths of MLLP. We use Python 3.5.1, matplotlib an numpy for sample paths simulation. (i) Choose a finite interval [, t]. Suppose our aim is to have n-simulate values at t 1 t/n,, t n 1 (n 1)t/n, an t n t. (ii) Generate a vector G of size n of gamma variates such that G (G 1, G 2,, G n ), with G i gamma(, t i t i 1 ). (iii) Now generate an inepenent vector of size n of stanar α-stable ranom numbers S (S 1, S 2,, S n ). (iv) Compute Y k G 1/α k S k, k 1 an M i i k1 Y k. Then M 1, M 2,, M n enote n- simulate values from MLLP. (a) α.9 an 1. (b) α.4 an 1. Figure 1: Sample Paths of MLLP 7 Generalization As the integer orer moments of MLLP are not finite, we can generalize MLLP to a tempere MLLP by taking suborination with respect to a tempere stable suborinator (hereafter 9

10 TSS). TSS are infinitely ivisible, have exponentially ecaying tail probabilities an have all moments finite. These properties are obtaine on the cost of self-similarity. Let S α,µ (t) be the TSS with inex α (, 1) an tempering parameter µ. TSS are obtaine by exponential tempering in the istribution of stable processes (see e.g. Rosinski, 27). Further, the LT of ensity of TSS (see Meerschaert et al. 213) is We efine, tempere MLLP as follows. E(e usα,µ(t) ) e t ( (u+µ) α µ α ). (7.1) Definition 7.1 (Tempere MLLP). Tempere MLLP is efine by suborinating TSS with inepenent gamma process such that M α,,µ(t) : S α,µ (G,1 (t)), α (, 1), >, µ >, t. (7.2) By using the stanar conitioning argument with (7.1), we obtain ( ) t E(e um α,,µ (t) ). (7.3) µ α + (µ + u) α Without proof, we state the following result for marginal ensity an Lévy ensity for tempere MLLP. The proof is analogous to Propositions 3.1 an 3.2. Proposition 7.1. The marginal ensity function an Lévy ensity for tempere MLLP are given by f M α,,µ (t)(x) t e µx Further, the Lévy ensity is k (µ α ) k Γ(t + k) x α(t+k) 1 Γ(k + 1)Γ(t) Γ(α(t + k)). ν (x) αe µx x E α,1((µ α )x α ), x >. We have EG,1 (t) t/ an var(g,1 (t)) t/ 2. Further, E(S α,µ (t)) αt an µ 1 α var(s α,µ (t)) α(1 α)t. Using the information above an with the help of conitioning argument, we have E(Mα,,µ (t)) αt an var(m α(1 α)t µ 2 α µ 1 α α,,µ (t)) + α2 t. Moreover, for µ 2 α 2 µ 2 α tempere MLLP all orer moments are finite an tail probability ecays exponentially. References Applebaum, D. 29. Lévy Processes an Stochastic Calculus. 2n e., Cambrige University Press, Cambrige, U.K. 1

11 Barnorff-Nielsen, O. E. 2. Probability ensities an Lévy ensities. Research 18, Aarhus Univ., Centre for Mathematical Physics an Stochastics (MaPhySto). Beghin, L. an Orsingher, E. 29. Fractional Poisson processes an relate ranom motions. Electron. J. Probab., Bertoin, J Lévy Processes. Cambrige University Press, Cambrige. Cahoy, D. O., Uchaikin, V. V. an Woyczynski, W. A. 21. Parameter estimation for fractional Poisson processes. J. Statist. Plann. Inf Feller, W Introuction to Probability Theory an its Applications. Vol. II. John Wiley, New York. Gorenflo, R., Kilbas, A. A., Mainari, F. an Rogosin, S. V Mittag-Leffler Functions, Relate Topics an Applications. Springer, New York. Kanter, M Stable ensities uner change of scale an total variation inequlities. Annals of Probability 3, Kozubowski, T. J., Meerschaert, M. M. an Pogorski, K. 26. motion. Av. Appl. Prob. 38, Fractional Laplace Laskin, N. 23. Fractional Poisson process. Commun. Nonlinear Sci. Numer. Simul. 8, Meerschaert, M. M., Nane, E., Vellaisamy, P Transient anamolous subiffusions on boune omains. Proc. Amer. Math. Soc., 141, Meerschaert, M. M., Nane, E., Vellaisamy, P The fractional Poisson process an the inverse stable suborinator. Electron. J. Probab., 16, Pillai, R. N On Mittag-Leffler functions an relate istributions. Ann. Inst. Statist. Math. 42, Repin, O.N. an Saichev, A.I. 2 Fractional Poisson law. Raiophys. an Quantum Electronics. 43, Rosiński, J., 27. Tempering stable processes. Stochastic Process Appl. 117, Sato, K.-I Lévy Processes an Infinitely Divisible Distributions. Cambrige University Press. 11