Tempered Fractional Calculus
|
|
- Amos Underwood
- 5 years ago
- Views:
Transcription
1 Tempered Fractional Calculus International Symposium on Fractional PDEs: Theory, Numerics and Applications 3 5 June 2013 Mark M. Meerschaert Department of Statistics and Probability Michigan State University mcubed@stt.msu.edu Partially supported by NSF grant DMS and NIH grant R01-EB
2 Abstract Fractional Calculus has a close relation with Probability. Random walks with heavy tails converge to stable stochastic processes, whose probability densities solve space-fractional diffusion equations. Continuous time random walks, with heavy tailed waiting times between particle jumps, converge to non-markovian stochastic limits, whose probability densities solve time-fractional diffusion equations. Time-fractional derivatives and integrals of Brownian motion produce fractional Brownian motion, a useful model in many applications. Fractional derivatives and integrals are convolutions with a power law. Including an exponential term leads to tempered fractional derivatives and integrals. Tempered stable processes are the limits of random walk models where the power law probability of long jumps is tempered by an exponential factor. These random walks converge to tempered stable stochastic process limits, whose probability densities solve tempered fractional diffusion equations. Tempered power law waiting times lead to tempered fractional time derivatives, which have proven useful in geophysics. Applying this idea to Brownian motion leads to tempered fractional Brownian motion, a new stochastic process that can exhibit semi-long range dependence. The increments of this process, called tempered fractional Gaussian noise, provide a useful new stochastic model for wind speed data.
3 Acknowledgments Farzad Sabzikar, Statistics and Probability, Michigan State Boris Baeumer, Maths & Stats, University of Otago, New Zealand Anna K. Panorska, Mathematics and Statistics, U of Nevada Parthanil Roy, Indian Statistical Institute, Kolkata Hans-Peter Scheffler, Math, Universität Siegen, Germany Qin Shao, Mathematics, University of Toledo, Ohio Alla Sikorskii, Statistics and Probability, Michigan State Yong Zhang, Desert Research Institute, Las Vegas Nevada
4 New Books Stochastic Models for Fractional Calculus Mark M. Meerschaert and Alla Sikorskii De Gruyter Studies in Mathematics 43, 2012 Advanced graduate textbook ISBN Mathematical Modeling, 4th Edition Mark M. Meerschaert Academic Press, Elsevier, 2013 Advanced undergraduate / beginning graduate textbook (new sections on particle tracking and anomalous diffusion) ISBN
5 Continuous time random walks R d X 2 X 3 0 X 1 W 1 W 2 W 3 W 4 W 5 X 4 X 5 t A random particle arrives at location S(n) = X X n at time T n = W 1 + +W n. After N t = max{n 0 : T n t} jumps, particle location is S(N t ).
6 CTRW limit theory If P(X n > x) x α and E(X n ) = 0 then n 1/α S(nu) A(u). The α-stable limit A(u) has pdf p(x,u) with Fourier transform ˆp(k,u) = e uψ A(k) and Fourier symbol ψ A (k) = (ik) α for 1 < α 2. If P(W n > t) t β then n 1/β T nu D u. The β-stable limit has pdf g(t,u) with Laplace transform g(s,u) = e uψ D(s) and Laplace symbol ψ D (s) = s β for 0 < β < 1. Inverse process: n β N nt E t = inf{u > 0 : D u > t}, and then n β/α S(N nt ) = n β/α S(n β n β N nt ) (n β ) 1/α S(n β E t ) A(E t ). Since E t has a pdf h(u,t) with LT h(u,s) = s 1 ψ D (s)e uψd(s), the CTRW limit A(E t ) has pdf q(x,t) = 0 p(x,u)h(u,t)du u P(A(u) = x E t = u)p(e t = u).
7 Probability and fractional calculus CTRW limit pdf has Fourier-Laplace transform q(k,s) = = 0 ˆp(k,u) h(u,s)du 0 eu(ik)α s β 1 e usβ du = sβ 1 s β (ik) α Then s β q(k,s) s β 1 = (ik) α q(k,s) which inverts to β t q(x,t) = Dα x q(x,t), using Caputo on the LHS and Riemann-Liouville on the RHS. Note β t codes power law waiting times, D α x power law jumps. Reduces to traditional diffusion for α = 2 and β = 1.
8 Anomalous diffusion processes Here α = 2.0,1.7 (left/right) and β = 1.0,0.9 (top/bottom). At At t t AEt AEt t t
9 Tempered fractional calculus Tempered stable D t has pdf with LT g(t,u) = e uψ D(s) where the Laplace symbol ψ D (s) = (λ+s) β λ β for some λ > 0 (small). Tempered stable A(u) has pdf with FT ˆp(k,u) = e uψ A(k) where the Fourier symbol ψ A (k) = (λ+ik) α λ α ikαλ α 1 with λ > 0. Again E t has pdf with LT h(u,s) = s 1 ψ D (s)e uψ D(s), and again A(E t ) has pdf q(x,t) = 0 p(x,u)h(u,t)du with FLT q(k,s) = s 1 ψ D (s) ψ D (s) ψ A (k) Then ψ D (s) q(k,s) s 1 ψ D (s) = ψ A (k) q(k,s). Invert to get β,λ t q(x,t) = D α,λ x q(x,t). Tempered jumps P(X n > x) x α e λx and waiting times.
10 Space fractional tempered stable Tempered stable Lévy motion with α = λ = λ = λ = λ =
11 Tempered power laws in finance AMZN stock price changes fit a tempered power law model P(X > x) x 0.6 e 0.3x for x large ln(p(x> x)) oo o o o ooooo ooooooooooo oooo oooooo oooo ooooo oooo o o o ooo o o o o o o o o o o ln(x)
12 Tempered power laws in hydrology Tempered power law model P(X > x) x 0.6 e 5.2x for increments in hydraulic conductivity at the MADE site. ln(p(x>x)) ln(x)
13 Tempered power laws in atmospheric science Tempered power law model P(X > x) x 0.2 e 0.01x for daily precipitation data at Tombstone AZ. ln(p(x>x)) ln(x)
14 Tempered stable pdf in macroeconomics Density data One-step BARMA forecast errors for annual inflation rates, and a tempered stable with α = 1.1 and λ = 12 (rough fit).
15 Tempered time-fractional diffusion model Fitted concentration data from a 3-D supercomputer simulation. ADE fit uses α = 2,β = 1. Without cutoff uses λ = 0.
16 Numerical codes for fractional diffusion For α > 0 we have D α x f(x) = lim h 0h α α f(x) where ( ) ( ) α α f(x) = ( 1) j α f(x jh), = Γ(α+1) j j j!γ(α j +1) j=0 Use to construct explicit/implicit finite difference codes. These codes are mass-preserving since j=0 ( ) α ( 1) j = 0. j Change jh to (j α )h for unconditional stability (implicit Euler). Use operator splitting for 2-d, 3-d, or reaction term. Finite element and finite volume methods also available.
17 Numerical example Exact solution p(x,t) = e t x 3 and numerical approximation to p(x, t) = d(x) 1.8 p(x,t) t x 1.8 +r(x,t) on 0 < x < 1 with p(x,0) = x 3, p(0,t) = 0, p(1,t) = e t, r(x,t) = (1+x)e t x 3, d(x) = Γ(2.2)x 2.8 / Exact CN
18 Tempered fractional diffusion codes For 0 < α < 1, lim h 0 h α α,λ f(x) = D α,λ x f(x) where ( ) α,λ α f(x) = ( 1) j e λjh f(x jh) (1 e λh ) α f(x). j j=0 Codes are mass-preserving since (by the Binomial formula) ( ) α ( 1) j e λjh = (1 e λh ) α j j=0 Stable, consistent implicit Euler codes: shift jh (j α )h. Crank-Nicolson codes are O( t 2 + x). Apply Richardson extrapolation to get O( t 2 + x 2 ). For 1 < α < 2, h α α,λ f(x) D α,λ x f(x)+αλ α 1 f (x).
19 Numerical example Exact solution u(x,t) = x β e λx t /Γ(1+β) and Crank-Nicolson solution (extrapolated) to t p(x,t) = c(x)d α,λ x p(x,t)+r(x,t) on 0 < x < 1 with α = 1,6, λ = 2, β = 2.8, p(x,0) = x β e λx /Γ(β +1), p(0,t) = 0, p(1,t) = e λ t /Γ(β +1), c(x) = x α Γ(1+β α)/γ(β +1), and r(x,t) = r 1 (x,t)e λx t Γ(1+β α)/γ(β +1) where r 1 (x,t) = (1 α)λα x α+β Γ(β +1) + αβλα 1 x α+β 1 Γ(β) 2x β Γ(1+β α). t x CN Max Error Rate CNX Max Error Rate 1/10 1/ /20 1/ /40 1/ /80 1/
20 Fractional derivatives: Integral forms In the simplest case 0 < α < 1 the generator form is D α x f(x) = α ( ) f(x) f(x y) y α 1 dy. Γ(1 α) 0 Check: Since f(x y) has FT e ikyˆf(k), then the RHS has FT α Γ(1 α) 0 ( 1 e iky ) ˆf(k)y α 1 dy. For λ > 0 it is not hard to compute [Prop in MS 2012] α Γ(1 α) 0 ( 1 e (λ+ik)y ) y α 1 dy = (λ+ik) α and then (let λ 0) D α x f(x) has FT (ik)αˆf(k). For 1 < α < 2 D α α(α 1) xf(x) = Γ(2 α) and again the RHS has FT (ik) αˆf(k). 0 ( f(x y) f(x)+yf (x) ) y α 1 dy
21 Tempered fractional derivatives In the simplest case 0 < α < 1 the generator form is D α,λ x f(x) = α ( ) f(x) f(x y) e λy y α 1 dy. Γ(1 α) 0 Check: The integral on the RHS has Fourier transform α ( 1 e iky ) ˆf(k)e λy y α 1 dy Γ(1 α) 0 α ( = (1 e (λ+ik)y ) (1 e λy ) ) ˆf(k)y α 1 dy Γ(1 α) 0 which reduces to [(λ+ik) α λ α ]ˆf(k). For 1 < α < 2 we have D α,λ α(α 1) x f(x) = Γ(2 α) 0 ( f(x y) f(x)+yf (x) ) e λy y α 1 dy and the RHS has FT [(λ+ik) α λ α ikαλ α 1 ]ˆf(k).
22 Laplace transform approach Since e λt f(t) has Laplace transform f(s λ) and since D α t LT s α f(s), we see that [ e λt f(t) ] dt = s α f(s λ). 0 e st D α t Using the shift property one more time, we see that 0 e st e λt D α [ t e λt f(t) ] dt = (s+λ) α f(s). Then for 0 < α 1 we have f(t) has D α,λ t f(t) = e λt D α t and for 1 < α 2 we have [ e λt f(t) ] λ α f(t) D α,λ x f(x) = e λx D α x[ e λx f(x) ] λ α f(x) αλ α 1 f (x).
23 Fractional derivatives (other integral forms) Recall that in the simplest case 0 < α < 1 the generator form D α x f(x) = 1 Γ(1 α) 0 ( f(x) f(x y) ) αy α 1 dy. Integrate by parts udv = uv vdu with u = f(x) f(x y) and dv = αy α 1 dy to get the Caputo form 1 Γ(1 α) 0 f (x y)y α dy = 1 Γ(1 α) f (u)(x u) α + du where (x u) + = (x u) for x > u and (x u) + = 0 otherwise. Move the derivative outside to get the Riemann-Liouville form 1 d Γ(1 α) dx Caputo is I 1 α x 0 f(x y)y α dy = 1 d Γ(1 α) dx f(u)(x u) α + dy. D 1 x f(x) and Riemann-Liouville is D1 x I1 α f(x) where I α x f(x) = 1 Γ(α) f(u)(x u)α 1 + du. x
24 Fractional Brownian motion Given Z n iid with mean zero and finite variance, the time series ( ) X t = α α Z n = ( 1) j Z j t j j=0 is long range dependent: E(X t X t+j ) j 2H 2 with H = 1/2 α. Then n H (X 1 + +X [nt] ) B H (t) fractional Brownian motion. Here B H (t) = α t B(t) α t B(0) where B(t) is a Brownian motion: B H (t) = 1 Γ(1 α) + [ ] (t u) α + (0 u) α + B(du) with B(du) = B (u)du in the distributional sense (Caputo). Hence FBM is the fractional integral of the white noise B(du).
25 Tempered fractional Brownian motion Tempered fractional Brownian motion with 1/2 < α < 1/2: B α,λ (t) = + [ e λ(t u) + (t u) α + e λ(0 u) +(0 u) α ] + B(du). TFBM is the tempered fractional integral of white noise B(du): I α,λ t f(t) = 1 Γ(α) f(u)e λ(t u) +(t u) α 1 + du. Tempered fractional Gaussian noise: X t = B α,λ (t) B α,λ (t 1). Semi-long range dependence: E(X t X t+j ) j 2H 2 for j small to moderate, then E(X t X t+j ) j 2 as j, where H = 1/2 α. For any Z n iid with mean zero and finite variance, the time series X t = α,λ Z n also exhibits semi-long range dependence.
26 FBM and tempered FBM FBM (thin line) and TFBM (thick line) with same B(du) for H = 0.3,λ = 0.03 (left) and H = 0.7,λ = 0.01 (right). X(t) X(t) t t Sample paths are Hölder continuous of order H. Scaling: B H (ct) d = c H B H (t) and B α,λ (ct) d = c H B α,cλ (t).
27 Semi-long range dependence The covariance functions E(X t X t+j ) for FGN and TFGN with H = 0.7 (and λ = 0.001) are quite similar until j gets large lag
28 Spectral density FGN spectral density blows up at low frequencies for H > 1/2: f H (ω) = 1 2π j= e iωj E(X t X t+j ) ω 1 2H as ω 0. For TFGN with λ 0 we get a similar result f α,λ (ω) ω 2 [λ 2 +ω 2 ] H+1/2 as ω 0. Hence f α,λ (ω) ω 1 2H for moderate frequencies, but remains bounded for very low frequencies (Davenport spectrum).
29 Spectral density comparison Spectral density for FGN and TFGN with H = 0.7,λ = frequency
30 Davenport spectrum Kolmogorov invented FBM to model turbulence in the inertial subrange. The Davenport spectrum f(ω) ω 2 /[1+ω 2 ] H+1/2 extends the model to include production and dissipation. Since TFBM has the same spectral density, it provides a comprehensive stochastic model for turbulence. Figure reproduced from Beaupuits et al. (2004).
31 Davenport spectrum for wind gusts Spectral density of wind gustiness from Davenport (1961). TFBM provides a stochastic model for the Davenport spectrum.
32 Wind power study Spectral density of wind speed at the Chajnantor radio telescope site in Chile.
33 Summary Tempered power laws Tempered fractional derivatives Numerical methods Tempered fractional Brownian motion Davenport model for wind speed
34 References 1. I.B. Aban, M.M. Meerschaert, and A.K. Panorska (2006) Parameter Estimation for the Truncated Pareto Distribution. Journal of the American Statistical Association: Theory and Methods. 101(473), P. Abry, P. Goncalves and P. Flandrin (1995) Wavelets, spectrum analysis and 1/f processes. In Wavelets and statistics, Springer New York. 3. B. Baeumer and M.M. Meerschaert (2010) Tempered stable Levy motion and transient super-diffusion. Journal of Computational and Applied Mathematics 233, J. P. Pérez Beaupuits, A. Otárola, F. T. Rantakyrö, R. C. Rivera, S. J. E. Radford, and L.-Å. Nyman (2004) Analysis of wind data gathered at Chajnantor. ALMA Memo Á. Cartea and D. del Castillo-Negrete (2007) Fluid limit of the continuous-time random walk with general Lévy jump distribution functions. Phys. Rev. E 76, A. Chakrabarty and M. M. Meerschaert (2011) Tempered stable laws as random walk limits. Statistics and Probability Letters 81(8), A. V. Chechkin, V. Yu. Gonchar, J. Klafter and R. Metzler (2005) Natural cutoff in Lévy flights caused by dissipative nonlinearity. Phys. Rev. E 72, S. Cohen and J. Rosiński (2007) Gaussian approximation of multivariate Lévy processes with applications to simulation of tempered stable processes, Bernoulli 13, A. G. Davenport (1961) The spectrum of horizontal gustiness near the ground in high winds. Quarterly Journal of the Royal Meteorological Society 87, P. Flandrin (1989) On the spectrum of fractional Brownian motions. IEEE Trans. on Info. Theory IT-35,
35 11. R. N. Mantegna and H. E. Stanley (1994) Stochastic process with ultraslow convergence to a Gaussian: The truncated Lévy flight. Phys. Rev. Lett. 73, M.M. Meerschaert and H.P. Scheffler (2001) Limit Theorems for Sums of Independent Random Vectors: Heavy Tails in Theory and Practice. Wiley Interscience, New York. 13. M.M. Meerschaert and H.P. Scheffler (2004) Limit theorems for continuous-time random walks with infinite mean waiting times. J. Appl. Probab. 41(3), Meerschaert, M.M., Scheffler, H.P. (2008) Triangular array limits for continuous time random walks. Stoch. Proc. Appl. 118(9), Meerschaert, M.M., Y. Zhang and B. Baeumer (2008) Tempered anomalous diffusion in heterogeneous systems. Geophys. Res. Lett. 35, L M.M. Meerschaert and A. Sikorskii (2012) Stochastic Models For Fractional Calculus. De Gruyter, Berlin/Boston. 17. Meerschaert, M.M., P. Roy and Q. Shao (2012) Parameter estimation for tempered power law distributions. Communications in Statistics Theory and Methods 41(10), M.M. Meerschaert (2013) Mathematical Modeling. 4th Edition, Academic Press, Boston. 19. M.M. Meerschaert and F. Sabzikar (2013) Tempered fractional Brownian motion. Preprint at Rosiński, J. (2007), Tempering stable processes. Stoch. Proc. Appl. 117, C. Tadjeran, M.M. Meerschaert, H.P. Scheffler (2006) A second order accurate numerical approximation for the fractional diffusion equation. Journal of Computational Physics 213(1),
36 Simulating tempered stable laws Simulation codes for stable random variates are widely available. If X > 0 has stable density density f(x), TS density is e λx f(x) f λ (x) = 0 e λu f(u)du Take Y exp(λ) independent of X, (X i,y i ) IID with (X,Y). Let N = min{n : X n Y n }. Then X N f λ (x). Proof: Compute P(X N x) = P(X x X Y) by conditioning, then take d/dx to verify.
37 Triangular array scheme (SPL 2011) Take P(X > x) Cx α with 1 < α < 2. Triangular array limit [nt] k=1 n 1/α X k b (n) t is stable. Define tempering variables: P(Z > u) = u α A(t) u r α 1 e λr dr Replace n 1/α X k by Z k if n 1/α X k > Z k. Triangular array limit is tempered stable. Exponential tempering: sum of α and α 1 tempered stables.
38 Tail estimation (CIS 2012) Hill-type estimator: Assume P(X > x) Cx α e λx for x large, use order statistics X (1) X (2) X (n). Conditional MLE given X (n k+1) > L X (n k) : T 1 : = T 2 : = 1 = k i=1 k i=1 k i=1 (logx (n i+1) logl) (X (n i+1) L) x (n i+1) kx (n i+1) + ˆα(T 2 T 1 x (n i+1) ) ˆλ = (k ˆαT 1 )/T 2 Ĉ = k n Lˆα eˆλl R code available at
39 Testing for pure power law tail (JASA 06) Null hypthesis H 0 : P(X > x) = Cx α Pareto for x > L. Test based on extreme value theory rejects H 0 if X (1) < ( nc lnq ) 1/α where α, C can be estimated using Hill s estimator ˆα H = k 1 k {lnx (n i+1) lnx (n k) } i=1 Ĉ H = (k/n)(x (n k) )ˆα H 1 Simple p-value formula p = exp{ ncx α (n) }.
Covariance structure of continuous time random walk limit processes
Covariance structure of continuous time random walk limit processes Alla Sikorskii Department of Statistics and Probability Michigan State University Fractional Calculus, Probability and Non-local Operators:
More informationContinuous Time Random Walk Limits in Bounded Domains
Continuous Time Random Walk Limits in Bounded Domains Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY May 16, 2012 Acknowledgement Boris Baeumer, University of Otago, New Zealand.
More informationCauchy Problems in Bounded Domains and Iterated Processes
Cauchy Problems in Bounded Domains and Iterated Processes Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 18, 21 Continuous time random walks R d X 2 X 3 X 1 J 1 J 2 J 3 J
More informationCTRW Limits: Governing Equations and Fractal Dimensions
CTRW Limits: Governing Equations and Fractal Dimensions Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 19-23, 2013 Joint work with Z-Q. Chen, M. D Ovidio, M.M. Meerschaert,
More informationA tempered multiscaling stable model to simulate transport in regional scale fractured media
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl043609, 2010 A tempered multiscaling stable model to simulate transport in regional scale fractured media Yong Zhang,
More informationAnomalous diffusion and turbulence
Anomalous diffusion and turbulence MURI Project Review Aberdeen Proving Ground 4 December 2018 Mark M. Meerschaert University Distinguished Professsor Department of Statistics and Probability Michigan
More informationCauchy Problems Solved by Running Subordinate Processes
Cauchy Problems Solved by Running Subordinate Processes Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 18, 21 Continuous time random walks R d X 2 X 3 X 1 J 1 J 2 J 3 J 4
More informationReflected spectrally negative stable processes and fractional Cauchy problems
Reflected spectrally negative stable processes and fractional Cauchy problems 7th International Conference on Lévy Processes: Theory and Applications 15 19 July 2013 Mark M. Meerschaert Department of Statistics
More informationSemi-Markov approach to continuous time random walk limit processes
Semi-Markov approach to continuous time random walk limit processes Joint Mathematics Meetings San Antonio, TX January 13, 2015 Mark M. Meerschaert Michigan State University East Lansing, Michigan, USA
More informationJournal of Computational and Applied Mathematics. Tempered stable Lévy motion and transient super-diffusion
Journal of Computational and Applied Mathematics 33 1 438 448 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam Tempered
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationCorrelation Structure of Fractional Pearson Diffusions
Correlation Structure of Fractional Pearson Diffusions Alla Sikorskii Department of Statistics and Probability Michigan State University International Symposium on Fractional PDEs: Theory, Numerics and
More informationBoundary Conditions for Tempered Fractional Diffusion
Boundary Conditions for Tempered Fractional Diffusion Anna Lischke a,, James F. Kelly b, Mark M. Meerschaert b a Division of Applied Mathematics, Brown University, Providence, RI 02912 b Department of
More informationStatistics and Probability Letters
Statistics and Probability Letters 83 (23) 83 93 Contents lists available at SciVerse ScienceDirect Statistics and Probability Letters journal homepage: www.elsevier.com/locate/stapro Fractal dimension
More informationFractional Cauchy Problems for time-changed Processes
Fractional Cauchy Problems for time-changed Processes Erkan Nane Department of Mathematics and Statistics Auburn University February 17, 211 Erkan Nane (Auburn University) Fractional Cauchy Problems February
More informationUCSD ECE250 Handout #27 Prof. Young-Han Kim Friday, June 8, Practice Final Examination (Winter 2017)
UCSD ECE250 Handout #27 Prof. Young-Han Kim Friday, June 8, 208 Practice Final Examination (Winter 207) There are 6 problems, each problem with multiple parts. Your answer should be as clear and readable
More informationTheory of fractional Lévy diffusion of cold atoms in optical lattices
Theory of fractional Lévy diffusion of cold atoms in optical lattices, Erez Aghion, David Kessler Bar-Ilan Univ. PRL, 108 230602 (2012) PRX, 4 011022 (2014) Fractional Calculus, Leibniz (1695) L Hospital:
More informationSpace-Time Duality and Medical Ultrasound
Space-Time Duality and Medical Ultrasound James F. Kelly and Mark M. Meerschaert A Workshop on Future Directions in Fractional Calculus Research and Applications Michigan State University October 21, 2016
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationCDA5530: Performance Models of Computers and Networks. Chapter 3: Review of Practical
CDA5530: Performance Models of Computers and Networks Chapter 3: Review of Practical Stochastic Processes Definition Stochastic ti process X = {X(t), t T} is a collection of random variables (rvs); one
More informationPhysica A. Extremal behavior of a coupled continuous time random walk
Physica A 39 (211) 55 511 Contents lists available at ScienceDirect Physica A journal homepage: www.elsevier.com/locate/physa Extremal behavior of a coupled continuous time random walk Rina Schumer a,,
More informationStationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos
N Cufaro Petroni: CYCLOTRONS 2007 Giardini Naxos, 1 5 October, 2007 1 Stationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos CYCLOTRONS 2007 Giardini Naxos, 1 5 October Nicola
More informationThe Laplace driven moving average a non-gaussian stationary process
The Laplace driven moving average a non-gaussian stationary process 1, Krzysztof Podgórski 2, Igor Rychlik 1 1 Mathematical Sciences, Mathematical Statistics, Chalmers 2 Centre for Mathematical Sciences,
More informationLecture 25: Large Steps and Long Waiting Times
Lecture 25: Large Steps and Long Waiting Times Scribe: Geraint Jones (and Martin Z. Bazant) Department of Economics, MIT Proofreader: Sahand Jamal Rahi Department of Physics, MIT scribed: May 10, 2005,
More informationON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP
Dedicated to Professor Gheorghe Bucur on the occasion of his 7th birthday ON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP EMIL POPESCU Starting from the usual Cauchy problem, we give
More informationLocal vs. Nonlocal Diffusions A Tale of Two Laplacians
Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2
More informationReview 1: STAT Mark Carpenter, Ph.D. Professor of Statistics Department of Mathematics and Statistics. August 25, 2015
Review : STAT 36 Mark Carpenter, Ph.D. Professor of Statistics Department of Mathematics and Statistics August 25, 25 Support of a Random Variable The support of a random variable, which is usually denoted
More informationSUBDYNAMICS OF FINANCIAL DATA FROM FRACTIONAL FOKKER PLANCK EQUATION
Vol. 4 (29) ACTA PHYSICA POLONICA B No 5 SUBDYNAMICS OF FINANCIAL DATA FROM FRACTIONAL FOKKER PLANCK EQUATION Joanna Janczura, Agnieszka Wyłomańska Institute of Mathematics and Computer Science Hugo Steinhaus
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationTopics in fractional Brownian motion
Topics in fractional Brownian motion Esko Valkeila Spring School, Jena 25.3. 2011 We plan to discuss the following items during these lectures: Fractional Brownian motion and its properties. Topics in
More informationCORRELATED CONTINUOUS TIME RANDOM WALKS
CORRELATED CONTINUOUS TIME RANDOM WALKS MARK M. MEERSCHAERT, ERKAN NANE, AND YIMIN XIAO Abstract. Continuous time random walks impose a random waiting time before each particle jump. Scaling limits of
More informationProduct measure and Fubini s theorem
Chapter 7 Product measure and Fubini s theorem This is based on [Billingsley, Section 18]. 1. Product spaces Suppose (Ω 1, F 1 ) and (Ω 2, F 2 ) are two probability spaces. In a product space Ω = Ω 1 Ω
More informationDefinition: Lévy Process. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 2: Lévy Processes. Theorem
Definition: Lévy Process Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 2: Lévy Processes David Applebaum Probability and Statistics Department, University of Sheffield, UK July
More informationA Class of Fractional Stochastic Differential Equations
Vietnam Journal of Mathematics 36:38) 71 79 Vietnam Journal of MATHEMATICS VAST 8 A Class of Fractional Stochastic Differential Equations Nguyen Tien Dung Department of Mathematics, Vietnam National University,
More informationFractional Quantum Mechanics and Lévy Path Integrals
arxiv:hep-ph/9910419v2 22 Oct 1999 Fractional Quantum Mechanics and Lévy Path Integrals Nikolai Laskin Isotrace Laboratory, University of Toronto 60 St. George Street, Toronto, ON M5S 1A7 Canada Abstract
More informationFigure 10.1: Recording when the event E occurs
10 Poisson Processes Let T R be an interval. A family of random variables {X(t) ; t T} is called a continuous time stochastic process. We often consider T = [0, 1] and T = [0, ). As X(t) is a random variable
More informationInvariance Properties of the Negative Binomial Lévy Process and Stochastic Self-similarity
International Mathematical Forum, 2, 2007, no. 30, 1457-1468 Invariance Properties of the Negative Binomial Lévy Process and Stochastic Self-similarity Tomasz J. Kozubowski Department of Mathematics &
More informationAnomalous diffusion equations
Anomalous diffusion equations Ercília Sousa Universidade de Coimbra. p.1/34 Outline 1. What is diffusion? 2. Standard diffusion model 3. Finite difference schemes 4. Convergence of the finite difference
More informationEstimation of the long Memory parameter using an Infinite Source Poisson model applied to transmission rate measurements
of the long Memory parameter using an Infinite Source Poisson model applied to transmission rate measurements François Roueff Ecole Nat. Sup. des Télécommunications 46 rue Barrault, 75634 Paris cedex 13,
More informationON THE NUMERICAL SOLUTION FOR THE FRACTIONAL WAVE EQUATION USING LEGENDRE PSEUDOSPECTRAL METHOD
International Journal of Pure and Applied Mathematics Volume 84 No. 4 2013, 307-319 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v84i4.1
More informationA New Fractional Process: A Fractional Non-homogeneous Poisson Process Enrico Scalas
A New Fractional Process: A Fractional Non-homogeneous Poisson Process Enrico Scalas University of Sussex Joint work with Nikolai Leonenko (Cardiff University) and Mailan Trinh Fractional PDEs: Theory,
More informationAnomalous Transport and Fluctuation Relations: From Theory to Biology
Anomalous Transport and Fluctuation Relations: From Theory to Biology Aleksei V. Chechkin 1, Peter Dieterich 2, Rainer Klages 3 1 Institute for Theoretical Physics, Kharkov, Ukraine 2 Institute for Physiology,
More informationQualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama
Qualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama Instructions This exam has 7 pages in total, numbered 1 to 7. Make sure your exam has all the pages. This exam will be 2 hours
More informationTime fractional Schrödinger equation
Time fractional Schrödinger equation Mark Naber a) Department of Mathematics Monroe County Community College Monroe, Michigan, 48161-9746 The Schrödinger equation is considered with the first order time
More informationChapter 5. Chapter 5 sections
1 / 43 sections Discrete univariate distributions: 5.2 Bernoulli and Binomial distributions Just skim 5.3 Hypergeometric distributions 5.4 Poisson distributions Just skim 5.5 Negative Binomial distributions
More informationµ X (A) = P ( X 1 (A) )
1 STOCHASTIC PROCESSES This appendix provides a very basic introduction to the language of probability theory and stochastic processes. We assume the reader is familiar with the general measure and integration
More informationShape of the return probability density function and extreme value statistics
Shape of the return probability density function and extreme value statistics 13/09/03 Int. Workshop on Risk and Regulation, Budapest Overview I aim to elucidate a relation between one field of research
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationCoupled continuous time random walks in finance
Physica A 37 (26) 114 118 www.elsevier.com/locate/physa Coupled continuous time random walks in finance Mark M. Meerschaert a,,1, Enrico Scalas b,2 a Department of Mathematics and Statistics, University
More informationVolatility and Returns in Korean Futures Exchange Markets
Volatility and Returns in Korean Futures Exchange Markets Kyungsik Kim*,, Seong-Min Yoon and Jum Soo Choi Department of Physics, Pukyong National University, Pusan 608-737, Korea Division of Economics,
More informationIsogeometric Analysis for the Fractional Laplacian
Isogeometric Analysis for the Fractional Laplacian Kailai Xu & Eric Darve CME500 May 14, 2018 1 / 26 Outline Fractional Calculus Isogeometric Analysis Numerical Experiment Regularity for the Fractional
More informationA Note on Certain Stability and Limiting Properties of ν-infinitely divisible distributions
Int. J. Contemp. Math. Sci., Vol. 1, 2006, no. 4, 155-161 A Note on Certain Stability and Limiting Properties of ν-infinitely divisible distributions Tomasz J. Kozubowski 1 Department of Mathematics &
More informationwhere r n = dn+1 x(t)
Random Variables Overview Probability Random variables Transforms of pdfs Moments and cumulants Useful distributions Random vectors Linear transformations of random vectors The multivariate normal distribution
More informationBMIR Lecture Series on Probability and Statistics Fall, 2015 Uniform Distribution
Lecture #5 BMIR Lecture Series on Probability and Statistics Fall, 2015 Department of Biomedical Engineering and Environmental Sciences National Tsing Hua University s 5.1 Definition ( ) A continuous random
More informationHW7 Solutions. f(x) = 0 otherwise. 0 otherwise. The density function looks like this: = 20 if x [10, 90) if x [90, 100]
HW7 Solutions. 5 pts.) James Bond James Bond, my favorite hero, has again jumped off a plane. The plane is traveling from from base A to base B, distance km apart. Now suppose the plane takes off from
More informationGeneralised Fractional-Black-Scholes Equation: pricing and hedging
Generalised Fractional-Black-Scholes Equation: pricing and hedging Álvaro Cartea Birkbeck College, University of London April 2004 Outline Lévy processes Fractional calculus Fractional-Black-Scholes 1
More informationFrom discrete time random walks to numerical methods for fractional order differential equations
From discrete time random walks to numerical methods for fractional order differential equations Click to edit Present s Name Dr. Christopher Angstmann SANUM 2016 Never Stand Still Faculty of Science School
More informationAn Efficient Numerical Method for Solving. the Fractional Diffusion Equation
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 1-12 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 An Efficient Numerical Method for Solving the Fractional
More informationSPACE-TIME FRACTIONAL DERIVATIVE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 133, Number 8, Pages 2273 2282 S 0002-9939(05)07949-9 Article electronically published on March 14, 2005 SPACE-TIME FRACTIONAL DERIVATIVE OPERATORS
More informationReview For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: = 0 : homogeneous equation.
Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: y y x y2 = 0 : homogeneous equation. x2 v = y dy, y = vx, and x v + x dv dx = v + v2. dx =
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationStochastic Stokes drift of a flexible dumbbell
Stochastic Stokes drift of a flexible dumbbell Kalvis M. Jansons arxiv:math/0607797v4 [math.pr] 22 Mar 2007 Department of Mathematics, University College London, Gower Street, London WC1E 6BT, UK January
More informationCDA6530: Performance Models of Computers and Networks. Chapter 3: Review of Practical Stochastic Processes
CDA6530: Performance Models of Computers and Networks Chapter 3: Review of Practical Stochastic Processes Definition Stochastic process X = {X(t), t2 T} is a collection of random variables (rvs); one rv
More informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
More informationPoisson Processes. Particles arriving over time at a particle detector. Several ways to describe most common model.
Poisson Processes Particles arriving over time at a particle detector. Several ways to describe most common model. Approach 1: a) numbers of particles arriving in an interval has Poisson distribution,
More informationSystems Driven by Alpha-Stable Noises
Engineering Mechanics:A Force for the 21 st Century Proceedings of the 12 th Engineering Mechanics Conference La Jolla, California, May 17-20, 1998 H. Murakami and J. E. Luco (Editors) @ASCE, Reston, VA,
More informationA Short Introduction to Diffusion Processes and Ito Calculus
A Short Introduction to Diffusion Processes and Ito Calculus Cédric Archambeau University College, London Center for Computational Statistics and Machine Learning c.archambeau@cs.ucl.ac.uk January 24,
More informationLower Tail Probabilities and Related Problems
Lower Tail Probabilities and Related Problems Qi-Man Shao National University of Singapore and University of Oregon qmshao@darkwing.uoregon.edu . Lower Tail Probabilities Let {X t, t T } be a real valued
More informationRandom variables and transform methods
Chapter Random variables and transform methods. Discrete random variables Suppose X is a random variable whose range is {,,..., } and set p k = P (X = k) for k =,,..., so that its mean and variance are
More informationVARIATIONAL METHODS FOR KIRCHHOFF TYPE PROBLEMS WITH TEMPERED FRACTIONAL DERIVATIVE
Electronic Journal of Differential Equations, Vol. 2018 (2018, No. 34, pp. 1 13. ISSN: 1072-6691. UL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu VAIATIONAL METHODS FO KICHHOFF TYPE POBLEMS
More informationMultivariate Generalized Ornstein-Uhlenbeck Processes
Multivariate Generalized Ornstein-Uhlenbeck Processes Anita Behme TU München Alexander Lindner TU Braunschweig 7th International Conference on Lévy Processes: Theory and Applications Wroclaw, July 15 19,
More information2 Functions of random variables
2 Functions of random variables A basic statistical model for sample data is a collection of random variables X 1,..., X n. The data are summarised in terms of certain sample statistics, calculated as
More informationMODELS FOR COMPUTER NETWORK TRAFFIC
MODELS FOR COMPUTER NETWORK TRAFFIC Murad S. Taqqu Boston University Joint work with Walter Willinger, Joshua Levy and Vladas Pipiras,... Web Site http://math.bu.edu/people/murad OUTLINE Background: 1)
More informationWeak Ergodicity Breaking
Weak Ergodicity Breaking Eli Barkai Bar-Ilan University 215 Ergodicity Ergodicity: time averages = ensemble averages. x = lim t t x(τ)dτ t. x = xp eq (x)dx. Ergodicity out of equilibrium δ 2 (, t) = t
More informationChing-Han Hsu, BMES, National Tsing Hua University c 2015 by Ching-Han Hsu, Ph.D., BMIR Lab. = a + b 2. b a. x a b a = 12
Lecture 5 Continuous Random Variables BMIR Lecture Series in Probability and Statistics Ching-Han Hsu, BMES, National Tsing Hua University c 215 by Ching-Han Hsu, Ph.D., BMIR Lab 5.1 1 Uniform Distribution
More informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More informationFractional Pearson diffusions
Fractional Pearson diffusions N.N. Leonenko, School of Mathematics, Cardiff University, UK joint work with Alla Sikorskii (Michigan State University, US), Ivan Papić and Nenad Šuvak (J.J. Strossmayer University
More informationNumerical Solution of Space-Time Fractional Convection-Diffusion Equations with Variable Coefficients Using Haar Wavelets
Copyright 22 Tech Science Press CMES, vol.89, no.6, pp.48-495, 22 Numerical Solution of Space-Time Fractional Convection-Diffusion Equations with Variable Coefficients Using Haar Wavelets Jinxia Wei, Yiming
More informationFractional Diffusion Theory and Applications Part II
Fractional Diffusion Theory and Applications Part II p. 1/2 Fractional Diffusion Theory and Applications Part II 22nd Canberra International Physics Summer School 28 Bruce Henry (Trevor Langlands, Peter
More informationMaximum Likelihood Drift Estimation for Gaussian Process with Stationary Increments
Austrian Journal of Statistics April 27, Volume 46, 67 78. AJS http://www.ajs.or.at/ doi:.773/ajs.v46i3-4.672 Maximum Likelihood Drift Estimation for Gaussian Process with Stationary Increments Yuliya
More informationDetermining and Forecasting High-Frequency Value-at-Risk by Using Lévy Processes
Determining and Forecasting High-Frequency Value-at-Risk by Using Lévy Processes W ei Sun 1, Svetlozar Rachev 1,2, F rank J. F abozzi 3 1 Institute of Statistics and Mathematical Economics, University
More informationTHREE-DIMENSIONAL HAUSDORFF DERIVATIVE DIFFUSION MODEL FOR ISOTROPIC/ANISOTROPIC FRACTAL POROUS MEDIA
Cai, W., et al.: Three-Dimensional Hausdorff Derivative Diffusion Model... THERMAL SCIENCE: Year 08, Vol., Suppl., pp. S-S6 S THREE-DIMENSIONAL HAUSDORFF DERIVATIVE DIFFUSION MODEL FOR ISOTROPIC/ANISOTROPIC
More informationSpace-Time Duality and High-Order Fractional Diffusion. Abstract
Space-Time Duality and High-Order Fractional Diffusion James F. Kelly and Mark M. Meerschaert Department of Statistics and Probability, arxiv:188.161v2 [cond-mat.stat-mech] 26 Jan 219 Michigan State University,
More informationStat410 Probability and Statistics II (F16)
Stat4 Probability and Statistics II (F6 Exponential, Poisson and Gamma Suppose on average every /λ hours, a Stochastic train arrives at the Random station. Further we assume the waiting time between two
More informationContinuous Random Variables and Continuous Distributions
Continuous Random Variables and Continuous Distributions Continuous Random Variables and Continuous Distributions Expectation & Variance of Continuous Random Variables ( 5.2) The Uniform Random Variable
More informationSimulation of stationary fractional Ornstein-Uhlenbeck process and application to turbulence
Simulation of stationary fractional Ornstein-Uhlenbeck process and application to turbulence Georg Schoechtel TU Darmstadt, IRTG 1529 Japanese-German International Workshop on Mathematical Fluid Dynamics,
More informationMOMENTS AND CUMULANTS OF THE SUBORDINATED LEVY PROCESSES
MOMENTS AND CUMULANTS OF THE SUBORDINATED LEVY PROCESSES PENKA MAYSTER ISET de Rades UNIVERSITY OF TUNIS 1 The formula of Faa di Bruno represents the n-th derivative of the composition of two functions
More informationMULTIDIMENSIONAL WICK-ITÔ FORMULA FOR GAUSSIAN PROCESSES
MULTIDIMENSIONAL WICK-ITÔ FORMULA FOR GAUSSIAN PROCESSES D. NUALART Department of Mathematics, University of Kansas Lawrence, KS 6645, USA E-mail: nualart@math.ku.edu S. ORTIZ-LATORRE Departament de Probabilitat,
More informationMATH Solutions to Probability Exercises
MATH 5 9 MATH 5 9 Problem. Suppose we flip a fair coin once and observe either T for tails or H for heads. Let X denote the random variable that equals when we observe tails and equals when we observe
More informationarxiv: v1 [math.ap] 26 Mar 2013
Analytic solutions of fractional differential equations by operational methods arxiv:134.156v1 [math.ap] 26 Mar 213 Roberto Garra 1 & Federico Polito 2 (1) Dipartimento di Scienze di Base e Applicate per
More informationStochastic solutions for fractional wave equations
Nonlinear Dyn DOI.7/s7-4-99-z ORIGINAL PAPER Stochastic solutions for fractional wave equations Mark M. Meerschaert René L. Schilling Alla Sikorskii Received: 6 November 3 / Accepted: 5 February 4 Springer
More informationName of the Student: Problems on Discrete & Continuous R.Vs
Engineering Mathematics 05 SUBJECT NAME : Probability & Random Process SUBJECT CODE : MA6 MATERIAL NAME : University Questions MATERIAL CODE : JM08AM004 REGULATION : R008 UPDATED ON : Nov-Dec 04 (Scan
More informationThe exponential distribution and the Poisson process
The exponential distribution and the Poisson process 1-1 Exponential Distribution: Basic Facts PDF f(t) = { λe λt, t 0 0, t < 0 CDF Pr{T t) = 0 t λe λu du = 1 e λt (t 0) Mean E[T] = 1 λ Variance Var[T]
More informationMATH 56A SPRING 2008 STOCHASTIC PROCESSES 197
MATH 56A SPRING 8 STOCHASTIC PROCESSES 197 9.3. Itô s formula. First I stated the theorem. Then I did a simple example to make sure we understand what it says. Then I proved it. The key point is Lévy s
More informationMean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations
Mean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations Ram Sharan Adhikari Assistant Professor Of Mathematics Rogers State University Mathematical
More informationFinite element approximation of the stochastic heat equation with additive noise
p. 1/32 Finite element approximation of the stochastic heat equation with additive noise Stig Larsson p. 2/32 Outline Stochastic heat equation with additive noise du u dt = dw, x D, t > u =, x D, t > u()
More informationSOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD
SOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD R. C. Mittal 1 and Ruchi Nigam 2 1 Department of Mathematics, I.I.T. Roorkee, Roorkee, India-247667. Email: rcmmmfma@iitr.ernet.in
More informationDepartment of Mathematics
Department of Mathematics Ma 3/13 KC Border Introduction to Probability and Statistics Winter 217 Lecture 13: The Poisson Process Relevant textbook passages: Sections 2.4,3.8, 4.2 Larsen Marx [8]: Sections
More informationExponential Distribution and Poisson Process
Exponential Distribution and Poisson Process Stochastic Processes - Lecture Notes Fatih Cavdur to accompany Introduction to Probability Models by Sheldon M. Ross Fall 215 Outline Introduction Exponential
More informationLeast Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises
Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,
More information