Affine processes on positive semidefinite matrices

Size: px
Start display at page:

Download "Affine processes on positive semidefinite matrices"

Transcription

1 Affine processes on positive semiefinite matrices Josef Teichmann (joint work with C. Cuchiero, D. Filipović an E. Mayerhofer) ETH Zürich Zürich, June 2010 C. Cuchiero, D. Filipović, E. Mayerhofer, J. Teichmann Affine processes on positive semiefinite matrices

2 1 Introuction Overview Motivation: Multivariate stochastic covariance moels Literature 2 Definition of affine processes on S + Affine processes on S + are Feller an regular 3 Main theorem Amissible parameters Ieas of the proof 4 1 / 69

3 Overview Motivation: Multivariate stochastic covariance moels Literature Introuction Affine processes on S +, the cone of positive semiefinite matrices, are stochastically continuous time-homogeneous Markov processes with state space S +, whose Laplace transform have exponential-affine epenence on the initial state, E x [e Tr(uXt)] = e φ(t,u) Tr(ψ(t,u)x). 2 / 69

4 Overview Motivation: Multivariate stochastic covariance moels Literature Introuction Affine processes on S +, the cone of positive semiefinite matrices, are stochastically continuous time-homogeneous Markov processes with state space S +, whose Laplace transform have exponential-affine epenence on the initial state, E x [e Tr(uXt)] = e φ(t,u) Tr(ψ(t,u)x). Aim: Full characterization of this class of processes 3 / 69

5 Overview Motivation: Multivariate stochastic covariance moels Literature Introuction Affine processes on S +, the cone of positive semiefinite matrices, are stochastically continuous time-homogeneous Markov processes with state space S +, whose Laplace transform have exponential-affine epenence on the initial state, E x [e Tr(uXt)] = e φ(t,u) Tr(ψ(t,u)x). Aim: Full characterization of this class of processes Necessary conitions on the parameters of the infinitesimal generator implie by the efinition of an affine process. 4 / 69

6 Overview Motivation: Multivariate stochastic covariance moels Literature Introuction Affine processes on S +, the cone of positive semiefinite matrices, are stochastically continuous time-homogeneous Markov processes with state space S +, whose Laplace transform have exponential-affine epenence on the initial state, E x [e Tr(uXt)] = e φ(t,u) Tr(ψ(t,u)x). Aim: Full characterization of this class of processes Necessary conitions on the parameters of the infinitesimal generator implie by the efinition of an affine process. Sufficient conitions for the existence of affine processes on S +. 5 / 69

7 Overview Motivation: Multivariate stochastic covariance moels Literature Introuction Affine processes on S +, the cone of positive semiefinite matrices, are stochastically continuous time-homogeneous Markov processes with state space S +, whose Laplace transform have exponential-affine epenence on the initial state, E x [e Tr(uXt)] = e φ(t,u) Tr(ψ(t,u)x). Aim: Full characterization of this class of processes Necessary conitions on the parameters of the infinitesimal generator implie by the efinition of an affine process. Sufficient conitions for the existence of affine processes on S +. Probabilistic properties like infinite ivisibility, excursion theory, etc. 6 / 69

8 Overview Motivation: Multivariate stochastic covariance moels Literature Motivation: Affine stochastic covariance moels One-imensional affine stochastic (co)variance moels (Heston [15], Barnorff-Nielsen Shepar moel [1], etc.). 7 / 69

9 Overview Motivation: Multivariate stochastic covariance moels Literature Motivation: Affine stochastic covariance moels One-imensional affine stochastic (co)variance moels (Heston [15], Barnorff-Nielsen Shepar moel [1], etc.). Risk neutral ynamics for the log-price process Y t an the R + -value variance process X t : X t = (b + βx t t + σ X t W t + J t, X 0 = x, ( Y t = r X ) t t + X t B t, Y 0 = y. 2 B, W : correlate Brownian motions, J: pure jump process, r: constant interest rate. 8 / 69

10 Overview Motivation: Multivariate stochastic covariance moels Literature Motivation: Affine stochastic covariance moels One-imensional affine stochastic (co)variance moels (Heston [15], Barnorff-Nielsen Shepar moel [1], etc.). Risk neutral ynamics for the log-price process Y t an the R + -value variance process X t : X t = (b + βx t t + σ X t W t + J t, X 0 = x, ( Y t = r X ) t t + X t B t, Y 0 = y. 2 B, W : correlate Brownian motions, J: pure jump process, r: constant interest rate. Efficient valuation of European options since the moment generating function is explicitly known (up to the solution of an ODE ) an is of the following form E x,y [ e zx t+vy t ] = e Φ(t,z,v)+Ψ(t,z,v)x+vy. 9 / 69

11 Overview Motivation: Multivariate stochastic covariance moels Literature Multivariate affine stochastic covariance moels Extension to multivariate stochastic covariance moels with the aim to / 69

12 Overview Motivation: Multivariate stochastic covariance moels Literature Multivariate affine stochastic covariance moels Extension to multivariate stochastic covariance moels with the aim to......capture the epenence structure between ifferent assets, 11 / 69

13 Overview Motivation: Multivariate stochastic covariance moels Literature Multivariate affine stochastic covariance moels Extension to multivariate stochastic covariance moels with the aim to......capture the epenence structure between ifferent assets,...obtain a consistent pricing framework for multi-asset options such as basket options, 12 / 69

14 Overview Motivation: Multivariate stochastic covariance moels Literature Multivariate affine stochastic covariance moels Extension to multivariate stochastic covariance moels with the aim to......capture the epenence structure between ifferent assets,...obtain a consistent pricing framework for multi-asset options such as basket options,...use them for moel-base ecision-making in the area of portfolio optimization an heging of correlation risk. 13 / 69

15 Overview Motivation: Multivariate stochastic covariance moels Literature Multivariate affine stochastic covariance moels Multivariate stochastic covariance moels consist of a -imensional logarithmic price process with risk-neutral ynamics ( Y t = r1 1 ) 2 X t iag t + X t B t, Y 0 = y, an stochastic covariation process. 0 X = Y, Y. B: -imensional Brownian motion. r: constant interest rate 1: the vector whose entries are all equal to one. X iag : the vector containing the iagonal entries of X. 14 / 69

16 Overview Motivation: Multivariate stochastic covariance moels Literature Multivariate affine stochastic covariance moels Multivariate stochastic covariance moels consist of a -imensional logarithmic price process with risk-neutral ynamics ( Y t = r1 1 ) 2 X t iag t + X t B t, Y 0 = y, an stochastic covariation process. 0 X = Y, Y. B: -imensional Brownian motion. r: constant interest rate 1: the vector whose entries are all equal to one. X iag : the vector containing the iagonal entries of X. In orer to qualify for a covariation process, X must be specifie as a process in S +. An affine ynamics for X is interesting out of several reasons. 15 / 69

17 Overview Motivation: Multivariate stochastic covariance moels Literature Multivariate affine stochastic covariance moels The following affine ynamics for X have been propose in the literature: X t = (b + HX t + X t H )t + X t W t Σ + Σ Wt Xt + J t, X 0 = x S +, b: suitably chosen matrix in S +, H, Σ: invertible matrices, W a stanar -matrix of Brownian motions possibly correlate with B, J a pure jump process whose compensator is an affine function of X. 16 / 69

18 Overview Motivation: Multivariate stochastic covariance moels Literature Multivariate affine stochastic covariance moels Analytic tractability: Uner some (mil) technical conitions, the following affine transform formula hols: [ ] E x,y e Tr(zXt)+v Y t = e Φ(t,z,v)+Tr(Ψ(t,z,v)x)+v y for appropriate arguments z S is an v C. The functions Φ an Ψ solve a system of generalize Riccati ODEs. Option pricing via Fourier methos 17 / 69

19 Overview Motivation: Multivariate stochastic covariance moels Literature Multivariate affine stochastic covariance moels Analytic tractability: Uner some (mil) technical conitions, the following affine transform formula hols: [ ] E x,y e Tr(zXt)+v Y t = e Φ(t,z,v)+Tr(Ψ(t,z,v)x)+v y for appropriate arguments z S is an v C. The functions Φ an Ψ solve a system of generalize Riccati ODEs. Option pricing via Fourier methos v y in the moment generating function correspons to a homogeneity assumption an implies that the covariation process X is a Markov process in its own filtration, which motivates the analysis of affine processes on S / 69

20 Overview Motivation: Multivariate stochastic covariance moels Literature Relate Literature Affine processes on R m + R n : D. Duffie, D. Filipović an W. Schachermayer [11]: Characterization of affine processes on R m + R n. D. Filipović an E. Mayerhofer [12]: Characterization of affine Diffusion processes. M. Keller-Ressel, W. Schachermayer an J. Teichmann [16]: Regularity of affine processes.... many other authors. 19 / 69

21 Overview Motivation: Multivariate stochastic covariance moels Literature Relate Literature Affine processes on R m + R n : D. Duffie, D. Filipović an W. Schachermayer [11]: Characterization of affine processes on R m + R n. D. Filipović an E. Mayerhofer [12]: Characterization of affine Diffusion processes. M. Keller-Ressel, W. Schachermayer an J. Teichmann [16]: Regularity of affine processes.... many other authors. Affine processes on S + - Theory of Wishart processes: M. Bru [3]: Existence an uniqueness of Wishart processes of type X t = (δi )t + X t W t + Wt Xt, X 0 S +, for δ 1. C. Donati-Martin, Y. Doumerc, H. Matsumoto an M. Yor [10]: Properties of Wishart processes. 20 / 69

22 Overview Motivation: Multivariate stochastic covariance moels Literature Relate Literature Affine processes on S + Applications from mathematical Finance: O. Barnorff-Nielsen an R. Stelzer [2]: Matrix-value Lévy riven Ornstein-Uhlenbeck processes of finite variation. B. Buraschi et al. [5, 4]: Correlation risk an portfolio optimization. J. Da Fonseca et al. [6, 7, 8, 9]: Multivariate stochastic covariance an option pricing. C. Gourieroux an R. Sufana [13, 14]: Wishart quaratic term structure moels. M. Leippol an F. Trojani [17]: S + -value affine jump iffusions an financial applications (multivariate option pricing, interest rate moels, etc.) 21 / 69

23 Definition of affine processes on S + Affine processes on S + are Feller an regular Setting an notation S : symmetric -matrices equippe with scalar prouct x, y = Tr(xy). 22 / 69

24 Definition of affine processes on S + Affine processes on S + are Feller an regular Setting an notation S : symmetric -matrices equippe with scalar prouct x, y = Tr(xy). S + : cone of symmetric -positive semiefinite matrices. 23 / 69

25 Definition of affine processes on S + Affine processes on S + are Feller an regular Setting an notation S : symmetric -matrices equippe with scalar prouct x, y = Tr(xy). S + : cone of symmetric -positive semiefinite matrices. S ++ : interior of S + in S, the cone of strictly positive efinite matrices. 24 / 69

26 Definition of affine processes on S + Affine processes on S + are Feller an regular Setting an notation S : symmetric -matrices equippe with scalar prouct x, y = Tr(xy). S + : cone of symmetric -positive semiefinite matrices. S ++ : interior of S + in S, the cone of strictly positive efinite matrices. Partial orer inuce by the cones: x y y x S + x y y x S ++. an 25 / 69

27 Definition of affine processes on S + Affine processes on S + are Feller an regular Setting an notation S : symmetric -matrices equippe with scalar prouct x, y = Tr(xy). S + : cone of symmetric -positive semiefinite matrices. S ++ : interior of S + in S, the cone of strictly positive efinite matrices. Partial orer inuce by the cones: x y y x S + an x y y x S ++. X time-homogeneous Markov process with state space S / 69

28 Definition of affine processes on S + Affine processes on S + are Feller an regular Setting an notation S : symmetric -matrices equippe with scalar prouct x, y = Tr(xy). S + : cone of symmetric -positive semiefinite matrices. S ++ : interior of S + in S, the cone of strictly positive efinite matrices. Partial orer inuce by the cones: x y y x S + an x y y x S ++. X time-homogeneous Markov process with state space S +. (P t ) t 0 : associate semigroup acting on boune measurable functions f : S + R, P t f (x) := E x [f (X t )] = f (ξ)p t (x, ξ), x S +. S + 27 / 69

29 Definition of affine processes on S + Affine processes on S + are Feller an regular Setting an notation S : symmetric -matrices equippe with scalar prouct x, y = Tr(xy). S + : cone of symmetric -positive semiefinite matrices. S ++ : interior of S + in S, the cone of strictly positive efinite matrices. Partial orer inuce by the cones: x y y x S + an x y y x S ++. X time-homogeneous Markov process with state space S +. (P t ) t 0 : associate semigroup acting on boune measurable functions f : S + R, P t f (x) := E x [f (X t )] = f (ξ)p t (x, ξ), x S +. S + X might be not conservative. Stanar extension of p t to S + with cemetery. { } 28 / 69

30 Definition of affine processes on S + Affine processes on S + are Feller an regular Definition Definition The Markov process X is calle affine if 1 it is stochastically continuous, that is, lim s t p s (x, ) = p t (x, ) weakly on S + x S +, an for every t an 2 its Laplace transform has exponential-affine epenence on the initial state: P t e u,x = e u,ξ p t (x, ξ) = e φ(t,u) ψ(t,u),x, S + for all t an u, x S +, for some functions φ : R + S + R + an ψ : R + S + S / 69

31 Definition of affine processes on S + Affine processes on S + are Feller an regular Regularity an Feller property Definition The affine process X is calle regular if the erivatives φ(t, u) ψ(t, u) F (u) =, R(u) = (1) t t t=0+ t=0+ exist an are continuous at u = / 69

32 Definition of affine processes on S + Affine processes on S + are Feller an regular Regularity an Feller property Definition The affine process X is calle regular if the erivatives φ(t, u) ψ(t, u) F (u) =, R(u) = (1) t t t=0+ t=0+ exist an are continuous at u = 0. Theorem Let X be an affine process with state space S +. Then, we have: 1 X is regular. 2 X is a Feller process. 31 / 69

33 Definition of affine processes on S + Affine processes on S + are Feller an regular Remarks on the proof The proof relies on several properties of the functions φ an ψ: Semi-flow property: For all t, s R + φ(t + s, u) = φ(t, u) + φ(s, ψ(t, u)), (2) ψ(t + s, u) = ψ(s, ψ(t, u)). (3) 32 / 69

34 Definition of affine processes on S + Affine processes on S + are Feller an regular Remarks on the proof The proof relies on several properties of the functions φ an ψ: Semi-flow property: For all t, s R + φ(t + s, u) = φ(t, u) + φ(s, ψ(t, u)), (2) ψ(t + s, u) = ψ(s, ψ(t, u)). (3) Orer relations: For all u, v S + with v u an for all t 0, φ(t, v) φ(t, u) an ψ(t, v) ψ(t, u). 33 / 69

35 Definition of affine processes on S + Affine processes on S + are Feller an regular Remarks on the proof The proof relies on several properties of the functions φ an ψ: Semi-flow property: For all t, s R + φ(t + s, u) = φ(t, u) + φ(s, ψ(t, u)), (2) ψ(t + s, u) = ψ(s, ψ(t, u)). (3) Orer relations: For all u, v S + with v u an for all t 0, φ(t, v) φ(t, u) an ψ(t, v) ψ(t, u). Analyticity: u φ(t, u) an u ψ(t, u) are analytic on S / 69

36 Definition of affine processes on S + Affine processes on S + are Feller an regular Remarks on the proof The proof relies on several properties of the functions φ an ψ: Semi-flow property: For all t, s R + φ(t + s, u) = φ(t, u) + φ(s, ψ(t, u)), (2) ψ(t + s, u) = ψ(s, ψ(t, u)). (3) Orer relations: For all u, v S + with v u an for all t 0, φ(t, v) φ(t, u) an ψ(t, v) ψ(t, u). Analyticity: u φ(t, u) an u ψ(t, u) are analytic on S ++. This implies Lemma The function ψ(t, u) S ++ for all u S ++ an t / 69

37 Definition of affine processes on S + Affine processes on S + are Feller an regular Remarks on the proof Remark on the proof of the above theorem: Regularity can be prove by methos on the regularity of (semi-)flows from Montgommery an Zippin. 36 / 69

38 Definition of affine processes on S + Affine processes on S + are Feller an regular Remarks on the proof Remark on the proof of the above theorem: Regularity can be prove by methos on the regularity of (semi-)flows from Montgommery an Zippin. The { Feller property follows from regularity an the fact that e u,x u S ++ } is ense in C0 (S + ) (Stone-Weierstrass). 37 / 69

39 Definition of affine processes on S + Affine processes on S + are Feller an regular Consequences of regularity: From regularity we can conclue local characteristics: Since the following limit exists, Ae u,x = e u,x 1 lim t 0 t E( e u,xt x 1 ), we know that the limit is of Lévy-Khintchine form on the cone S + Rx. Due to the affine property the previous equation is also an equation for the erivatives of φ an ψ. By the Markov property (semigroup property) these equations turn into ODEs for φ an ψ. 38 / 69

40 Definition of affine processes on S + Affine processes on S + are Feller an regular The function φ an ψ are solutions of ODEs: Hence we obtain φ(t, u) = F (ψ(t, u)), t φ(0, u) = 0, (4) ψ(t, u) = R(ψ(t, u)), t ψ(0, u) = u S +, (5) where F an R are efine in (1). Due to the particular form of F an R we shall call them generalize Riccati equations. 39 / 69

41 Definition of affine processes on S + Affine processes on S + are Feller an regular The functions F an R the question of existence of affine processes Theorem Moreover, φ(t, u) an ψ(t, u) solve the ifferential equations (4) an (5), where F an R have the following form F (u) = b, u + c (e u,ξ 1)m(ξ), S + \{0} R(u) = 2uαu + B (u) + γ ( ) µ(ξ) e u,ξ 1 + χ (ξ), u ξ 2 1. S + \{0} Conversely, let (α, b, β ij, c, γ, m, µ) be an amissible parameter set. Then there exists a unique affine process on S + with infinitesimal generator (6). 40 / 69

42 Main theorem Amissible parameters Ieas of the proof Infinitesimal generator Theorem If X is an affine process on S +, then its infinitesimal generator is of affine form: Af (x) = 1 A ijkl (x) 2 f (x) + b + B(x), f (x) (c + γ, x )f (x) 2 x ij x kl i,j,k,l + (f (x + ξ) f (x)) m(ξ) (6) S + \{0} + (f (x + ξ) f (x) χ(ξ), f (x) ) M(x, ξ), S + \{0} for some truncation function χ an amissible parameters ( α, b, B(x) = i,j β ij x ij, c, γ, m(ξ), M(x, ξ) = where A ijkl (x) = x ik α jl + x il α jk + x jk α il + x jl α ik. ) x, µ, ξ / 69

43 Main theorem Amissible parameters Ieas of the proof Relation to semimartingales Corollary Let X be a conservative affine process on S +. Then X is a semimartingale. Furthermore, there exists, possibly on an enlargement of the probability space, a -matrix of stanar Brownian motions W such that X amits the following representation ( t ) X t = x + b + χ(ξ)m(ξ) + B(X s) s, 0 S + \{0} t ( ) + XsWsΣ + Σ W s Xs 0 t ( ) + χ(ξ) µ X (t, ξ) (m(ξ) + M(X t, ξ)) t 0 S + \{0} t + (ξ χ(ξ))µ X (t, ξ), 0 S + \{0} where Σ is a matrix satisfying Σ Σ = α an µ X enotes the ranom measure associate with the jumps of X. 42 / 69

44 Main theorem Amissible parameters Ieas of the proof Amissible parameters linear iffusion coefficient: α S +, 43 / 69

45 Main theorem Amissible parameters Ieas of the proof Amissible parameters linear iffusion coefficient: α S +, linear jump coefficient: -matrix µ = (µ ij ) of finite signe measures on S + \ {0} with µ(e) S + for all E B(S + \ {0}) M(x, ξ) := x,µ(ξ) satisfies ξ 2 1 S + \{0} χ(ξ), u M(x, ξ) < for all x, u S + with x, u = 0, 44 / 69

46 Main theorem Amissible parameters Ieas of the proof Amissible parameters linear iffusion coefficient: α S +, linear jump coefficient: -matrix µ = (µ ij ) of finite signe measures on S + \ {0} with µ(e) S + for all E B(S + \ {0}) M(x, ξ) := x,µ(ξ) satisfies ξ 2 1 S + \{0} χ(ξ), u M(x, ξ) < for all x, u S + with x, u = 0, linear rift coefficient: the linear map B(x) = i,j βij x ij with β ij = β ji S satisfies B(x), u χ(ξ), u M(x, ξ) 0 S + \{0} for all x, u S + with x, u = 0, 45 / 69

47 Main theorem Amissible parameters Ieas of the proof Amissible parameters linear iffusion coefficient: α S +, linear jump coefficient: -matrix µ = (µ ij ) of finite signe measures on S + \ {0} with µ(e) S + for all E B(S + \ {0}) M(x, ξ) := x,µ(ξ) satisfies ξ 2 1 S + \{0} χ(ξ), u M(x, ξ) < for all x, u S + with x, u = 0, linear rift coefficient: the linear map B(x) = i,j βij x ij with β ij = β ji S satisfies B(x), u χ(ξ), u M(x, ξ) 0 S + \{0} for all x, u S + with x, u = 0, linear killing rate coefficient: γ S +, 46 / 69

48 Main theorem Amissible parameters Ieas of the proof Amissible parameters linear iffusion coefficient: α S +, linear jump coefficient: -matrix µ = (µ ij ) of finite signe measures on S + \ {0} with µ(e) S + for all E B(S + \ {0}) M(x, ξ) := x,µ(ξ) satisfies ξ 2 1 S + \{0} χ(ξ), u M(x, ξ) < for all x, u S + with x, u = 0, linear rift coefficient: the linear map B(x) = i,j βij x ij with β ij = β ji S satisfies B(x), u χ(ξ), u M(x, ξ) 0 S + \{0} for all x, u S + with x, u = 0, linear killing rate coefficient: γ S +, constant rift term: b ( 1)α, 47 / 69

49 Main theorem Amissible parameters Ieas of the proof Amissible parameters linear iffusion coefficient: α S +, linear jump coefficient: -matrix µ = (µ ij ) of finite signe measures on S + \ {0} with µ(e) S + for all E B(S + \ {0}) M(x, ξ) := x,µ(ξ) satisfies ξ 2 1 S + \{0} χ(ξ), u M(x, ξ) < for all x, u S + with x, u = 0, linear rift coefficient: the linear map B(x) = i,j βij x ij with β ij = β ji S satisfies B(x), u χ(ξ), u M(x, ξ) 0 S + \{0} for all x, u S + with x, u = 0, linear killing rate coefficient: γ S +, constant rift term: b ( 1)α, constant jump term: Borel measure m on S + ( ξ 1) m(ξ) <, S + \{0} \ {0} with 48 / 69

50 Main theorem Amissible parameters Ieas of the proof Amissible parameters linear iffusion coefficient: α S +, linear jump coefficient: -matrix µ = (µ ij ) of finite signe measures on S + \ {0} with µ(e) S + for all E B(S + \ {0}) M(x, ξ) := x,µ(ξ) satisfies ξ 2 1 S + \{0} χ(ξ), u M(x, ξ) < for all x, u S + with x, u = 0, linear rift coefficient: the linear map B(x) = i,j βij x ij with β ij = β ji S satisfies B(x), u χ(ξ), u M(x, ξ) 0 S + \{0} for all x, u S + with x, u = 0, linear killing rate coefficient: γ S +, constant rift term: b ( 1)α, constant jump term: Borel measure m on S + ( ξ 1) m(ξ) <, S + \{0} constant killing rate term: c R +. \ {0} with 49 / 69

51 Main theorem Amissible parameters Ieas of the proof Remark on the amissible parameters No constant iffusion part, linear part is of very specific form u, A(x)u = 4 x, uαu. 50 / 69

52 Main theorem Amissible parameters Ieas of the proof Remark on the amissible parameters No constant iffusion part, linear part is of very specific form u, A(x)u = 4 x, uαu. Very remarkable constant rift conition: b ( 1)α. For 2 the bounary of S + is curve an implies this relation between linear iffusion coefficient α an rift part b. 51 / 69

53 Main theorem Amissible parameters Ieas of the proof Remark on the amissible parameters No constant iffusion part, linear part is of very specific form u, A(x)u = 4 x, uαu. Very remarkable constant rift conition: b ( 1)α. For 2 the bounary of S + is curve an implies this relation between linear iffusion coefficient α an rift part b. Jumps escribe by m are of finite variation, for the linear jump part we have finite variation for the irections orthogonal to the bounary while parallel to the bounary general jump behavior is allowe. 52 / 69

54 Main theorem Amissible parameters Ieas of the proof Remark on the amissible parameters Parameters of prototype equation are amissible: X t = (b + HX t + X t H )t + X t W t Σ + Σ W t Xt + J t, 53 / 69

55 Main theorem Amissible parameters Ieas of the proof Remark on the amissible parameters Parameters of prototype equation are amissible: X t = (b + HX t + X t H )t + X t W t Σ + Σ W t Xt + J t, Σ R : α = Σ Σ S +, 54 / 69

56 Main theorem Amissible parameters Ieas of the proof Remark on the amissible parameters Parameters of prototype equation are amissible: X t = (b + HX t + X t H )t + X t W t Σ + Σ W t Xt + J t, Σ R : α = Σ Σ S +, b ( 1)Σ Σ, 55 / 69

57 Main theorem Amissible parameters Ieas of the proof Remark on the amissible parameters Parameters of prototype equation are amissible: X t = (b + HX t + X t H )t + X t W t Σ + Σ W t Xt + J t, Σ R : α = Σ Σ S +, b ( 1)Σ Σ, B(x) = Hx + H x: Hx + H x, u = 0 if u, x = 0, 56 / 69

58 Main theorem Amissible parameters Ieas of the proof Remark on the amissible parameters Parameters of prototype equation are amissible: X t = (b + HX t + X t H )t + X t W t Σ + Σ W t Xt + J t, Σ R : α = Σ Σ S +, b ( 1)Σ Σ, B(x) = Hx + H x: Hx + H x, u = 0 if u, x = 0, c = 0, γ = 0, 57 / 69

59 Main theorem Amissible parameters Ieas of the proof Remark on the amissible parameters Parameters of prototype equation are amissible: X t = (b + HX t + X t H )t + X t W t Σ + Σ W t Xt + J t, Σ R : α = Σ Σ S +, b ( 1)Σ Σ, B(x) = Hx + H x: Hx + H x, u = 0 if u, x = 0, c = 0, γ = 0, J a compoun Poisson process with intensity l + λ, x an jump istribution ν supporte on S + : m(ξ) = lν(ξ), M(x, ξ) = λ, x ν(ξ). 58 / 69

60 Main theorem Amissible parameters Ieas of the proof What are the new results Full characterization an exact assumptions uner which affine processes on S + actually exist. 59 / 69

61 Main theorem Amissible parameters Ieas of the proof What are the new results Full characterization an exact assumptions uner which affine processes on S + actually exist. Necessity an sufficiency of the remarkable rift conition b ( 1)α. 60 / 69

62 Main theorem Amissible parameters Ieas of the proof What are the new results Full characterization an exact assumptions uner which affine processes on S + actually exist. Necessity an sufficiency of the remarkable rift conition b ( 1)α. Extension of the moel class. 61 / 69

63 Main theorem Amissible parameters Ieas of the proof What are the new results Full characterization an exact assumptions uner which affine processes on S + actually exist. Necessity an sufficiency of the remarkable rift conition b ( 1)α. Extension of the moel class. General linear rift part B(x) = ij βij x ij. This allows epenency of the volatility of one asset on the other ones which is not possible for B(x) = Hx + xh. Example: = 2 an ( ) x22 x B(x) = 12. x 12 x / 69

64 Main theorem Amissible parameters Ieas of the proof What are the new results Full characterization an exact assumptions uner which affine processes on S + actually exist. Necessity an sufficiency of the remarkable rift conition b ( 1)α. Extension of the moel class. General linear rift part B(x) = ij βij x ij. This allows epenency of the volatility of one asset on the other ones which is not possible for B(x) = Hx + xh. Example: = 2 an ( ) x22 x B(x) = 12. x 12 x 11 Full generality of jumps (quaratic variation jumps parallel to the bounary). 63 / 69

65 Main theorem Amissible parameters Ieas of the proof Ieas an methos applie in the proof Necessary conitions: Since the limit exists Ae u,x = e u,x 1 lim t 0 t E( e u,xt x 1 ), we know that the limit is of Lévy-Khintchine form on the cone S + Rx. The structure of infinitely ivisible ranom variables on cones implies the form of F an R an it etermines the form of the infinitesimal generator since Ae u,x = ( F (u) R(u), x )e u,x. In orer to erive the conition on b we consier et(x t ). Invariance conitions for R + -value process imply then b ( 1)α. 64 / 69

66 Main theorem Amissible parameters Ieas of the proof Ieas an methos applie in the proof Sufficient conitions Problems: Unboune coefficients, x is not Lipschitz continuous at the bounary, infinitely active jumps. Regularization of the coefficients (replace x by x + ɛi ɛi ). Through stochastic invariance, we establish existence of an S + -value solution of the martingale problem for the generator of the regularize process Existence of an S + -value solution of the martingale problem for A. Uniqueness an existence of global solutions φ an ψ of the generalize Riccati equations (4) an (5) yiel uniqueness of the solution of the martingale problem Existence of an affine (Markov) process on S / 69

67 On S +, the laws of X t = δi t + X t W t + Wt Xt are those of non central Wishart istributions an are not infinitely ivisible. As a consequence of the conition on the constant rift, we obtain Corollary X is affine an its one-imensional marginal istributions are infinitely ivisible if an only if α = 0 or = / 69

68 Conclusion an Outlook Conclusion Characterization of S + -value affine processes. Exhaustive moel specification. Characterization of infinitely ivisible affine processes. Outlook Calibration an parameter estimation of multivariate stochastic covariance moels. Affine processes on symmetric cones. Calculation of stochastic quantities beyon marginals in the affine setting. Infinite imensional affine processes. 67 / 69

69 [1] Ole E. Barnorff-Nielsen an Neil Shephar. Moeling by Lévy processes for financial econometrics. In Lévy processes, pages Birkhäuser Boston, Boston, MA, [2] Ole Eiler Barnorff-Nielsen an Robert Stelzer. Positive-efinite matrix processes of finite variation. Probab. Math. Statist., 27(1):3 43, [3] Marie-France Bru. Wishart processes. Journal of Theoretical Probability, 4(4): , [4] B. Buraschi, A. Cieslak, an F. Trojani. Correlation risk an the term structure of interest rates. Working paper, University St.Gallen, [5] B. Buraschi, P. Porchia, an F. Trojani. Correlation risk an optimal portfolio choice. Working paper, University St.Gallen, [6] Jose Da Fonseca, Martino Grasseli, an Florian Ielpo. Heging (co)variance risk with variance swaps. Working paper, University St.Gallen, [7] Jose Da Fonseca, Martino Grasseli, an Florian Ielpo. Estimating the Wishart affine stochastic correlation moel using the empirical characteristic function. Working paper, [8] Jose Da Fonseca, Martino Grasseli, an Clauio Tebali. Option pricing when correlations are stochastic: an analytical framework. Review of Derivatives Research, 10(2): , [9] Jose Da Fonseca, Martino Grasseli, an Clauio Tebali. A multifactor volatility Heston moel. 68 / 69

70 J. Quant. Finance, 8(6): , [10] Catherine Donati-Martin, Yan Doumerc, Hiroyuki Matsumoto, an Marc Yor. Some properties of the Wishart processes an a matrix extension of the Hartman-Watson laws. Publ. Res. Inst. Math. Sci., 40(4): , [11] Darrel Duffie, Damir Filipović, an Walter Schachermayer. Affine processes an applications in finance. Ann. Appl. Prob., 13: , [12] D. Filipović an E. Mayerhofer. Affine iffusion processes: Theory an applications. In Avance Financial Moelling, volume 8 of Raon Ser. Comput. Appl. Math. Walter e Gruyter, Berlin, [13] Christian Gourieroux an Razvan Sufana. Wishart quaratic term structure moels. Working paper, CREST, CEPREMAP an University of Toronto, [14] Christian Gourieroux an Razvan Sufana. Derivative pricing with Wishart multivariate stochastic volatility: application to creit risk. Working paper, CREST, CEPREMAP an University of Toronto, [15] S. Heston. A close-form solution for options with stochastic volatility with appliactions to bon an currency options. Rev. of Financial Stuies. [16] Martin Keller-Ressel, Walter Schachermayer, an Josef Teichmann. Affine processes are regular. Preprint, [17] Markus Leippol an Fabio Trojani. Asset pricing with matrix affine jump iffusions. Working paper, / 69

Affine Processes are regular

Affine Processes are regular Affine Processes are regular Martin Keller-Ressel kemartin@math.ethz.ch ETH Zürich Based on joint work with Walter Schachermayer and Josef Teichmann (arxiv:0906.3392) Forschungsseminar Stochastische Analysis

More information

REGULARITY OF AFFINE PROCESSES ON GENERAL STATE SPACES

REGULARITY OF AFFINE PROCESSES ON GENERAL STATE SPACES REGULARITY OF AFFINE PROCESSES ON GENERAL STATE SPACES MARTIN KELLER-RESSEL, WALTER SCHACHERMAYER, AN JOSEF TEICHMANN Abstract. We consider a stochastically continuous, affine Markov process in the sense

More information

1 Math 285 Homework Problem List for S2016

1 Math 285 Homework Problem List for S2016 1 Math 85 Homework Problem List for S016 Note: solutions to Lawler Problems will appear after all of the Lecture Note Solutions. 1.1 Homework 1. Due Friay, April 8, 016 Look at from lecture note exercises:

More information

arxiv: v2 [math.pr] 23 Jun 2014

arxiv: v2 [math.pr] 23 Jun 2014 COMPUTATION OF COPULAS BY FOURIER METHODS ANTONIS PAPAPANTOLEON arxiv:08.26v2 [math.pr] 23 Jun 204 Abstract. We provide an integral representation for the (implied) copulas of dependent random variables

More information

Polynomial Preserving Jump-Diffusions on the Unit Interval

Polynomial Preserving Jump-Diffusions on the Unit Interval Polynomial Preserving Jump-Diffusions on the Unit Interval Martin Larsson Department of Mathematics, ETH Zürich joint work with Christa Cuchiero and Sara Svaluto-Ferro 7th General AMaMeF and Swissquote

More information

arxiv: v1 [math.pr] 4 Feb 2016

arxiv: v1 [math.pr] 4 Feb 2016 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

More information

Convergence of Random Walks

Convergence of Random Walks Chapter 16 Convergence of Ranom Walks This lecture examines the convergence of ranom walks to the Wiener process. This is very important both physically an statistically, an illustrates the utility of

More information

Some Examples. Uniform motion. Poisson processes on the real line

Some Examples. Uniform motion. Poisson processes on the real line Some Examples Our immeiate goal is to see some examples of Lévy processes, an/or infinitely-ivisible laws on. Uniform motion Choose an fix a nonranom an efine X := for all (1) Then, {X } is a [nonranom]

More information

AFFINE PROCESSES ARE REGULAR

AFFINE PROCESSES ARE REGULAR AFFINE PROCESSES ARE REGULAR MARTIN KELLER-RESSEL, WALTER SCHACHERMAYER, AND JOSEF TEICHMANN Abstract. We show that stochastically continuous, time-homogeneous affine processes on the canonical state space

More information

Optimal Control of Spatially Distributed Systems

Optimal Control of Spatially Distributed Systems Optimal Control of Spatially Distribute Systems Naer Motee an Ali Jababaie Abstract In this paper, we stuy the structural properties of optimal control of spatially istribute systems. Such systems consist

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Supermodular ordering of Poisson arrays

Supermodular ordering of Poisson arrays Supermodular ordering of Poisson arrays Bünyamin Kızıldemir Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological University 637371 Singapore

More information

Correction to: Yield curve shapes and the asymptotic short rate distribution in affine one-factor models

Correction to: Yield curve shapes and the asymptotic short rate distribution in affine one-factor models Finance Stoch (218) 22:53 51 https://doi.org/1.17/s78-18-359-5 CORRECTION Correction to: Yield curve shapes and the asymptotic short rate distribution in affine one-factor models Martin Keller-Ressel 1

More information

Levy Process and Infinitely Divisible Law

Levy Process and Infinitely Divisible Law Stat205B: Probability Theory (Spring 2003) Lecture: 26 Levy Process an Infinitely Divisible Law Lecturer: James W. Pitman Scribe: Bo Li boli@stat.berkeley.eu Levy Processes an Infinitely Divisible Law

More information

SOME LYAPUNOV TYPE POSITIVE OPERATORS ON ORDERED BANACH SPACES

SOME LYAPUNOV TYPE POSITIVE OPERATORS ON ORDERED BANACH SPACES Ann. Aca. Rom. Sci. Ser. Math. Appl. ISSN 2066-6594 Vol. 5, No. 1-2 / 2013 SOME LYAPUNOV TYPE POSITIVE OPERATORS ON ORDERED BANACH SPACES Vasile Dragan Toaer Morozan Viorica Ungureanu Abstract In this

More information

Introduction to Markov Processes

Introduction to Markov Processes Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav

More information

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine

More information

1. Aufgabenblatt zur Vorlesung Probability Theory

1. Aufgabenblatt zur Vorlesung Probability Theory 24.10.17 1. Aufgabenblatt zur Vorlesung By (Ω, A, P ) we always enote the unerlying probability space, unless state otherwise. 1. Let r > 0, an efine f(x) = 1 [0, [ (x) exp( r x), x R. a) Show that p f

More information

Energy behaviour of the Boris method for charged-particle dynamics

Energy behaviour of the Boris method for charged-particle dynamics Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion

More information

Convergence rates of moment-sum-of-squares hierarchies for optimal control problems

Convergence rates of moment-sum-of-squares hierarchies for optimal control problems Convergence rates of moment-sum-of-squares hierarchies for optimal control problems Milan Kora 1, Diier Henrion 2,3,4, Colin N. Jones 1 Draft of September 8, 2016 Abstract We stuy the convergence rate

More information

On a limit theorem for non-stationary branching processes.

On a limit theorem for non-stationary branching processes. On a limit theorem for non-stationary branching processes. TETSUYA HATTORI an HIROSHI WATANABE 0. Introuction. The purpose of this paper is to give a limit theorem for a certain class of iscrete-time multi-type

More information

CLARK-OCONE FORMULA BY THE S-TRANSFORM ON THE POISSON WHITE NOISE SPACE

CLARK-OCONE FORMULA BY THE S-TRANSFORM ON THE POISSON WHITE NOISE SPACE CLARK-OCONE FORMULA BY THE S-TRANSFORM ON THE POISSON WHITE NOISE SPACE YUH-JIA LEE*, NICOLAS PRIVAULT, AND HSIN-HUNG SHIH* Abstract. Given ϕ a square-integrable Poisson white noise functionals we show

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz A note on asymptotic formulae for one-imensional network flow problems Carlos F. Daganzo an Karen R. Smilowitz (to appear in Annals of Operations Research) Abstract This note evelops asymptotic formulae

More information

Discrete Operators in Canonical Domains

Discrete Operators in Canonical Domains Discrete Operators in Canonical Domains VLADIMIR VASILYEV Belgoro National Research University Chair of Differential Equations Stuencheskaya 14/1, 308007 Belgoro RUSSIA vlaimir.b.vasilyev@gmail.com Abstract:

More information

Jointly continuous distributions and the multivariate Normal

Jointly continuous distributions and the multivariate Normal Jointly continuous istributions an the multivariate Normal Márton alázs an álint Tóth October 3, 04 This little write-up is part of important founations of probability that were left out of the unit Probability

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Modelling and simulation of dependence structures in nonlife insurance with Bernstein copulas

Modelling and simulation of dependence structures in nonlife insurance with Bernstein copulas Moelling an simulation of epenence structures in nonlife insurance with Bernstein copulas Prof. Dr. Dietmar Pfeifer Dept. of Mathematics, University of Olenburg an AON Benfiel, Hamburg Dr. Doreen Straßburger

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

Optimal Control of Spatially Distributed Systems

Optimal Control of Spatially Distributed Systems Optimal Control of Spatially Distribute Systems Naer Motee an Ali Jababaie Abstract In this paper, we stuy the structural properties of optimal control of spatially istribute systems. Such systems consist

More information

Monotonicity for excited random walk in high dimensions

Monotonicity for excited random walk in high dimensions Monotonicity for excite ranom walk in high imensions Remco van er Hofsta Mark Holmes March, 2009 Abstract We prove that the rift θ, β) for excite ranom walk in imension is monotone in the excitement parameter

More information

A simple tranformation of copulas

A simple tranformation of copulas A simple tranformation of copulas V. Durrleman, A. Nikeghbali & T. Roncalli Groupe e Recherche Opérationnelle Créit Lyonnais France July 31, 2000 Abstract We stuy how copulas properties are moifie after

More information

arxiv: v1 [math-ph] 5 May 2014

arxiv: v1 [math-ph] 5 May 2014 DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic

More information

Markovian N-Server Queues (Birth & Death Models)

Markovian N-Server Queues (Birth & Death Models) Markovian -Server Queues (Birth & Death Moels) - Busy Perio Arrivals Poisson (λ) ; Services exp(µ) (E(S) = /µ) Servers statistically ientical, serving FCFS. Offere loa R = λ E(S) = λ/µ Erlangs Q(t) = number

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

NONLINEAR MARKOV SEMIGROUPS AND INTERACTING LÉVY TYPE PROCESSES

NONLINEAR MARKOV SEMIGROUPS AND INTERACTING LÉVY TYPE PROCESSES NONLINEAR MARKOV SEMIGROUPS AND INTERACTING LÉVY TYPE PROCESSES Vassili N. Kolokoltsov November 21, 2006 To appear in Journal of Statistical Physics. Abstract Semigroups of positivity preserving linear

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

Weighted Cramér-von Mises-type Statistics

Weighted Cramér-von Mises-type Statistics Weighte Cramér-von Mises-type Statistics LSTA, Université Paris VI 175 rue u Chevaleret F 7513 Paris, France (e-mail: p@ccr.jussieu.fr) Paul Deheuvels Abstract. We consier quaratic functionals of the multivariate

More information

ON THE RIEMANN EXTENSION OF THE SCHWARZSCHILD METRICS

ON THE RIEMANN EXTENSION OF THE SCHWARZSCHILD METRICS ON THE RIEANN EXTENSION OF THE SCHWARZSCHILD ETRICS Valerii Dryuma arxiv:gr-qc/040415v1 30 Apr 004 Institute of athematics an Informatics, AS R, 5 Acaemiei Street, 08 Chisinau, olova, e-mail: valery@ryuma.com;

More information

Modeling of Dependence Structures in Risk Management and Solvency

Modeling of Dependence Structures in Risk Management and Solvency Moeling of Depenence Structures in Risk Management an Solvency University of California, Santa Barbara 0. August 007 Doreen Straßburger Structure. Risk Measurement uner Solvency II. Copulas 3. Depenent

More information

AFFINE PROCESSES AND APPLICATIONS IN FINANCE

AFFINE PROCESSES AND APPLICATIONS IN FINANCE AFFINE PROCESSES AND APPLICATIONS IN FINANCE D. DUFFIE, D. FILIPOVIĆ, W. SCHACHERMAYER Abstract. We provide the definition and a complete characterization of regular affine processes. This type of process

More information

Lower bounds on Locality Sensitive Hashing

Lower bounds on Locality Sensitive Hashing Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,

More information

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics

More information

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

More information

REVERSIBILITY FOR DIFFUSIONS VIA QUASI-INVARIANCE. 1. Introduction We look at the problem of reversibility for operators of the form

REVERSIBILITY FOR DIFFUSIONS VIA QUASI-INVARIANCE. 1. Introduction We look at the problem of reversibility for operators of the form REVERSIBILITY FOR DIFFUSIONS VIA QUASI-INVARIANCE OMAR RIVASPLATA, JAN RYCHTÁŘ, AND BYRON SCHMULAND Abstract. Why is the rift coefficient b associate with a reversible iffusion on R given by a graient?

More information

Logarithmic spurious regressions

Logarithmic spurious regressions Logarithmic spurious regressions Robert M. e Jong Michigan State University February 5, 22 Abstract Spurious regressions, i.e. regressions in which an integrate process is regresse on another integrate

More information

PROBABILITY THEORY EXPONENTIAL FUNCTIONAL OF LÉVY PROCESSES: GENERALIZED WEIERSTRASS PRODUCTS AND WIENER-HOPF FACTORIZATION P. PATIE AND M.

PROBABILITY THEORY EXPONENTIAL FUNCTIONAL OF LÉVY PROCESSES: GENERALIZED WEIERSTRASS PRODUCTS AND WIENER-HOPF FACTORIZATION P. PATIE AND M. PROBABILITY THEORY EXPONENTIAL FUNCTIONAL OF LÉVY PROCESSES: GENERALIZED WEIERSTRASS PRODUCTS AND WIENER-HOPF FACTORIZATION P. PATIE AND M. SAVOV Abstract. In this note, we state a representation of the

More information

Permanent vs. Determinant

Permanent vs. Determinant Permanent vs. Determinant Frank Ban Introuction A major problem in theoretical computer science is the Permanent vs. Determinant problem. It asks: given an n by n matrix of ineterminates A = (a i,j ) an

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

7.1 Support Vector Machine

7.1 Support Vector Machine 67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to

More information

Introduction. A Dirichlet Form approach to MCMC Optimal Scaling. MCMC idea

Introduction. A Dirichlet Form approach to MCMC Optimal Scaling. MCMC idea Introuction A Dirichlet Form approach to MCMC Optimal Scaling Markov chain Monte Carlo (MCMC quotes: Metropolis et al. (1953, running coe on the Los Alamos MANIAC: a feasible approach to statistical mechanics

More information

Exact and high order discretization schemes. Wishart processes and their affine extensions

Exact and high order discretization schemes. Wishart processes and their affine extensions for Wishart processes and their affine extensions CERMICS, Ecole des Ponts Paris Tech - PRES Université Paris Est Modeling and Managing Financial Risks -2011 Plan 1 Motivation and notations 2 Splitting

More information

SEPARABLE TERM STRUCTURES AND THE MAXIMAL DEGREE PROBLEM. 1. Introduction This paper discusses arbitrage-free separable term structure (STS) models

SEPARABLE TERM STRUCTURES AND THE MAXIMAL DEGREE PROBLEM. 1. Introduction This paper discusses arbitrage-free separable term structure (STS) models SEPARABLE TERM STRUCTURES AND THE MAXIMAL DEGREE PROBLEM DAMIR FILIPOVIĆ Abstract. This paper discusses separable term structure diffusion models in an arbitrage-free environment. Using general consistency

More information

Multivariate Generalized Ornstein-Uhlenbeck Processes

Multivariate 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 information

Parameter estimation: A new approach to weighting a priori information

Parameter estimation: A new approach to weighting a priori information Parameter estimation: A new approach to weighting a priori information J.L. Mea Department of Mathematics, Boise State University, Boise, ID 83725-555 E-mail: jmea@boisestate.eu Abstract. We propose a

More information

Multi-View Clustering via Canonical Correlation Analysis

Multi-View Clustering via Canonical Correlation Analysis Technical Report TTI-TR-2008-5 Multi-View Clustering via Canonical Correlation Analysis Kamalika Chauhuri UC San Diego Sham M. Kakae Toyota Technological Institute at Chicago ABSTRACT Clustering ata in

More information

Markov Chains in Continuous Time

Markov Chains in Continuous Time Chapter 23 Markov Chains in Continuous Time Previously we looke at Markov chains, where the transitions betweenstatesoccurreatspecifietime- steps. That it, we mae time (a continuous variable) avance in

More information

A. Exclusive KL View of the MLE

A. Exclusive KL View of the MLE A. Exclusive KL View of the MLE Lets assume a change-of-variable moel p Z z on the ranom variable Z R m, such as the one use in Dinh et al. 2017: z 0 p 0 z 0 an z = ψz 0, where ψ is an invertible function

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

arxiv: v1 [math.dg] 3 Jun 2016

arxiv: v1 [math.dg] 3 Jun 2016 FOUR-DIMENSIONAL EINSTEIN MANIFOLDS WITH SECTIONAL CURVATURE BOUNDED FROM ABOVE arxiv:606.057v [math.dg] Jun 206 ZHUHONG ZHANG Abstract. Given an Einstein structure with positive scalar curvature on a

More information

On the number of isolated eigenvalues of a pair of particles in a quantum wire

On the number of isolated eigenvalues of a pair of particles in a quantum wire On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in

More information

Linear Algebra- Review And Beyond. Lecture 3

Linear Algebra- Review And Beyond. Lecture 3 Linear Algebra- Review An Beyon Lecture 3 This lecture gives a wie range of materials relate to matrix. Matrix is the core of linear algebra, an it s useful in many other fiels. 1 Matrix Matrix is the

More information

STABILITY ESTIMATES FOR SOLUTIONS OF A LINEAR NEUTRAL STOCHASTIC EQUATION

STABILITY ESTIMATES FOR SOLUTIONS OF A LINEAR NEUTRAL STOCHASTIC EQUATION TWMS J. Pure Appl. Math., V.4, N.1, 2013, pp.61-68 STABILITY ESTIMATES FOR SOLUTIONS OF A LINEAR NEUTRAL STOCHASTIC EQUATION IRADA A. DZHALLADOVA 1 Abstract. A linear stochastic functional ifferential

More information

On Decentralized Optimal Control and Information Structures

On Decentralized Optimal Control and Information Structures 2008 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June 11-13, 2008 FrC053 On Decentralize Optimal Control an Information Structures Naer Motee 1, Ali Jababaie 1 an Bassam

More information

Linear ODEs. Types of systems. Linear ODEs. Definition (Linear ODE) Linear ODEs. Existence of solutions to linear IVPs.

Linear ODEs. Types of systems. Linear ODEs. Definition (Linear ODE) Linear ODEs. Existence of solutions to linear IVPs. Linear ODEs Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems p. 1 Linear ODEs Types of systems Definition (Linear ODE) A linear ODE is a ifferential equation

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

Stable and compact finite difference schemes

Stable and compact finite difference schemes Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long

More information

arxiv: v2 [math.st] 15 Apr 2016

arxiv: v2 [math.st] 15 Apr 2016 Maximum Likelihoo Estimation for Wishart processes Aurélien Alfonsi, Ahme Kebaier, an Clément Rey CERMICS, École es Ponts, UPE, Projet MathRisk ENPC-INRIA-UMLV, Champs-sur-Marne, France Université Paris

More information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

More information

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device Close an Open Loop Optimal Control of Buffer an Energy of a Wireless Device V. S. Borkar School of Technology an Computer Science TIFR, umbai, Inia. borkar@tifr.res.in A. A. Kherani B. J. Prabhu INRIA

More information

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13) Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to

More information

Mathematical Review Problems

Mathematical Review Problems Fall 6 Louis Scuiero Mathematical Review Problems I. Polynomial Equations an Graphs (Barrante--Chap. ). First egree equation an graph y f() x mx b where m is the slope of the line an b is the line's intercept

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Quantile function expansion using regularly varying functions

Quantile function expansion using regularly varying functions Quantile function expansion using regularly varying functions arxiv:705.09494v [math.st] 9 Aug 07 Thomas Fung a, an Eugene Seneta b a Department of Statistics, Macquarie University, NSW 09, Australia b

More information

A PAC-Bayesian Approach to Spectrally-Normalized Margin Bounds for Neural Networks

A PAC-Bayesian Approach to Spectrally-Normalized Margin Bounds for Neural Networks A PAC-Bayesian Approach to Spectrally-Normalize Margin Bouns for Neural Networks Behnam Neyshabur, Srinah Bhojanapalli, Davi McAllester, Nathan Srebro Toyota Technological Institute at Chicago {bneyshabur,

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS

SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS GEORGE A HAGEDORN AND CAROLINE LASSER Abstract We investigate the iterate Kronecker prouct of a square matrix with itself an prove an invariance

More information

BOUNDEDNESS IN A THREE-DIMENSIONAL ATTRACTION-REPULSION CHEMOTAXIS SYSTEM WITH NONLINEAR DIFFUSION AND LOGISTIC SOURCE

BOUNDEDNESS IN A THREE-DIMENSIONAL ATTRACTION-REPULSION CHEMOTAXIS SYSTEM WITH NONLINEAR DIFFUSION AND LOGISTIC SOURCE Electronic Journal of Differential Equations, Vol. 016 (016, No. 176, pp. 1 1. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu BOUNDEDNESS IN A THREE-DIMENSIONAL ATTRACTION-REPULSION

More information

Analysis of tandem fluid queues

Analysis of tandem fluid queues Analysis of tanem flui queues Małgorzata M. O Reilly School of Mathematics University of Tasmania Tas 7001, Australia malgorzata.oreilly@utas.eu.au Werner Scheinhart Department of Applie Mathematics University

More information

EXISTENCE OF INVARIANT MANIFOLDS FOR STOCHASTIC EQUATIONS IN INFINITE DIMENSION

EXISTENCE OF INVARIANT MANIFOLDS FOR STOCHASTIC EQUATIONS IN INFINITE DIMENSION EXISTENCE OF INVARIANT MANIFOLDS FOR STOCHASTIC EQUATIONS IN INFINITE DIMENSION DAMIR FILIPOVIĆ AND JOSEF TEICHMANN Abstract. We provie a Frobenius type existence result for finite-imensional invariant

More information

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

Self-normalized Martingale Tail Inequality

Self-normalized Martingale Tail Inequality Online-to-Confience-Set Conversions an Application to Sparse Stochastic Banits A Self-normalize Martingale Tail Inequality The self-normalize martingale tail inequality that we present here is the scalar-value

More information

Numerical Methods with Lévy Processes

Numerical Methods with Lévy Processes Numerical Methods with Lévy Processes 1 Objective: i) Find models of asset returns, etc ii) Get numbers out of them. Why? VaR and risk management Valuing and hedging derivatives Why not? Usual assumption:

More information

Spurious Significance of Treatment Effects in Overfitted Fixed Effect Models Albrecht Ritschl 1 LSE and CEPR. March 2009

Spurious Significance of Treatment Effects in Overfitted Fixed Effect Models Albrecht Ritschl 1 LSE and CEPR. March 2009 Spurious Significance of reatment Effects in Overfitte Fixe Effect Moels Albrecht Ritschl LSE an CEPR March 2009 Introuction Evaluating subsample means across groups an time perios is common in panel stuies

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix

More information

arxiv: v4 [math.pr] 27 Jul 2016

arxiv: v4 [math.pr] 27 Jul 2016 The Asymptotic Distribution of the Determinant of a Ranom Correlation Matrix arxiv:309768v4 mathpr] 7 Jul 06 AM Hanea a, & GF Nane b a Centre of xcellence for Biosecurity Risk Analysis, University of Melbourne,

More information

Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs

Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs Preprints of the 8th IFAC Worl Congress Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs Guoong Shi ACCESS Linnaeus Centre, School of Electrical Engineering, Royal

More information

arxiv: v4 [q-fin.pr] 24 Jul 2012

arxiv: v4 [q-fin.pr] 24 Jul 2012 The explicit Laplace transform for the Wishart process Alessandro Gnoatto artino Grasselli January 31, 13 Abstract arxiv:117.748v4 [q-fin.pr 4 Jul 1 We derive the explicit formula for the joint Laplace

More information

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements Aaptive Gain-Scheule H Control of Linear Parameter-Varying Systems with ime-delaye Elements Yoshihiko Miyasato he Institute of Statistical Mathematics 4-6-7 Minami-Azabu, Minato-ku, okyo 6-8569, Japan

More information

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. Uniqueness for solutions of ifferential equations. We consier the system of ifferential equations given by x = v( x), () t with a given initial conition

More information

Tractability results for weighted Banach spaces of smooth functions

Tractability results for weighted Banach spaces of smooth functions Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March

More information