Multivariate Generalized Ornstein-Uhlenbeck Processes

Size: px
Start display at page:

Download "Multivariate Generalized Ornstein-Uhlenbeck Processes"

Transcription

1 Mulivariae Generalized Ornsein-Uhlenbeck Processes Ania Behme and Alexander Lindner Absrac De Haan and Karandikar [12] inroduced generalized Ornsein Uhlenbeck processes as one-dimensional processes (V 0 which are basically characerized by he fac ha for each h > 0 he equidisanly sampled process (V nh n N0 saisfies he random recurrence equaion V nh = A (nh,nh V (nh + B (nh,nh, n N, where (A (nh,nh, B (nh,nh n N is an i.i.d. sequence wih posiive A 0,h for each h > 0. We generalize his concep o a mulivariae seing and use i o define mulivariae generalized Ornsein Uhlenbeck (MGOU processes which occur o be characerized by a saring random variable and some Lévy process (X, Y in R m m R m. The sochasic differenial equaion an MGOU process saisfies is also derived. We furher sudy invarian subspaces and irreducibiliy of he models generaed by MGOU processes and use his o give necessary and sufficien condiions for he exisence of sricly saionary soluions of MGOU processes under some exra condiions Mahemaics subjec classificaion. 60G10, 60G51. Key words and phrases. generalized Ornsein-Uhlenbeck process, invarian subspace, irreducible model, Lévy process, muliplicaive Lévy process, sochasic exponenial. 1 Inroducion Le (ξ, η = (ξ, η 0 be a bivariae Lévy process and V 0 a random variable, independen of (ξ, η. Then, following De Haan and Karandikar [12] and Carmona e al. [6], he onedimensional process (V 0, given by V = e (V ξ 0 + e ξ s dη s, 0, (1.1 Technische Universiä Braunschweig, Insiu für Mahemaische Sochasik, Pockelssraße 14, D Braunschweig, Germany, a.behme@u-bs.de or a.lindner@u-bs.de, el.:+49/531/ or +49/531/ , fax:+49/531/

2 is called a generalized Ornsein Uhlenbeck (GOU process. We refer o Maller e al. [19] for furher informaion and references regarding GOU processes. A key feaure of hese processes is ha for any h > 0, he random sequence (V nh n N0 saisfies he random recurrence equaion V nh = A (nh,nh V (nh + B (nh,nh, n N, where (A (nh,nh, B (nh,nh n N is an i.i.d. (independen and idenically disribued sequence wih A 0,h > 0 almos surely. Wihou assuming independence of V 0 and (ξ, η, processes of he form (1.1 are he only processes having his propery for any h > 0 and which saisfy some naural exra condiions, as shown by De Haan and Karandikar [12]. In he presen paper we exend he seing of De Haan and Karandikar [12] o random marices wih real valued enries, i.e. we aim o consruc a process (V 0, wih V = (V (i,j 1 i m R m l 1 j l in coninuous ime which fulfills he random recurrence equaion V = A s, V s + B s, a.s., 0 s, (1.2 for random funcionals (A s, 0 s, (B s, 0 s such ha A s, R m m and B s, R m l, he A s, are supposed o be non-singular and (A (nh,nh, B (nh,nh, n N, are i.i.d. for all h > 0. We also aim o characerize all processes in coninuous ime which have his propery and saisfy some naural exra condiions. The obained soluions will be called mulivariae generalized Ornsein Uhlenbeck (MGOU processes since hey exend he key feaure of one-dimensional generalized Ornsein Uhlenbeck processes canonically. Observe ha he quesion of when a soluion of (1.2 exiss can be reaed separaely for each column of (V 0. Thus, if no saed oherwise, for simpliciy we se l = 1 hroughou his paper, hence V and B s, are elemens in R m. To moivae he menioned exra condiions, following he lines of De Haan and Karandikar [12] observe ha he condiion of (1.2 o hold for all 0 s yields A u, V u + B u, = V = A s, V s + B s, = A s, A u,s V u + A s, B u,s + B s,, 0 u s. Assuming ha (A s,, B s, 0 s is unique now leads o Assumpion 1(a given below while exending he i.i.d. propery of (A (nh,nh, B (nh,nh, n N, for all h > 0 ino he coninuous ime seing yields he requiremens 1(b and (c. Finally, i is naural o impose ha (A 0, 0 and (B 0, 0 are coninuous in probabiliy a 0 since his, ogeher wih 1(a,(b and (c, implies he exisence of càdlàg modificaions of he processes (A 0 := (A 0, 0 and (B 0 := (B 0, 0 as will be shown in Lemma 2.1 below. This moivaes Assumpion 1(d below. We denoe he se of all inverible real m m-marices by GL(R, m, he ideniy marix by I and by 0 he vecor (or marix having only zero enries. We wrie d = for equaliy in disribuion and P- lim for limis in probabiliy. Assumpion 1. Suppose he GL(R, m R m -valued random funcional (A s,, B s, 0 s wih A, = I and B, = 0 a.s. for all 0 saisfies he following four condiions. 2

3 (a For all 0 u s almos surely A u, = A s, A u,s and B u, = A s, B u,s + B s,. (1.3 (b For all 0 a b c d he families of random marices {(A s,, B s,, a s b} and {(A s,, B s,, c s d} are independen. (c For all 0 s i holds (A s,, B s, d = (A 0, s, B 0, s. (1.4 (d I holds P- lim 0 A 0, = I and P- lim 0 B 0, = 0. (1.5 The firs main resul of he paper will be a characerizaion of all random funcionals (A s,, B s, 0 s which saisfy Assumpion 1, in erms of appropriae driving Lévy processes. This will be achieved in Theorem 3.1 and hen be used o define MGOU processes as processes which saisfy (1.2 wih (A s,, B s, 0 s subjec o Assumpion 1. I will be also shown in Secion 3 ha MGOU processes saisfy he sochasic differenial equaion (SDE dv = du V + dl for appropriae Lévy processes U and L if he saring random variable V 0 is independen of (A 0,, B 0, 0, exending a corresponding one-dimensional resul of De Haan and Karandikar [12]. A new aspec compared o he one-dimensional GOU process is he possibiliy of he exisence of affine subspaces H of R m which are invarian under he model (1.2 in he sense ha A s, H + B s, H holds for all 0 s. In Secion 4 we give necessary and sufficien condiions for he exisence of an invarian affine subspace of he model (1.2 and show ha given he exisence of a d-dimensional invarian affine subspace H, afer an appropriae orhogonal ransformaion of he underlying space, he MGOU process wih V 0 H consiss of an (m d-dimensional consan process and an R d -valued MGOU process. Subsequenly in Secion 5 sricly saionary soluions of MGOU processes are reaed. Under some exra condiions we give necessary and sufficien condiions for heir exisence and deermine heir form, exending corresponding one-dimensional resuls of Behme e al. [3] and Lindner and Maller [18]. The proofs for he resuls of Secions 3 5 are given in Secions 6 8. A crucial ingredien for he derivaion of he necessary and sufficien condiions for saionariy are he resuls on saionary soluions of random recurrence equaions by Bougerol and Picard [4]. Secion 8 also conains several auxiliary resuls abou mulivariae sochasic exponenials. Some preliminary resuls are colleced in Secion 2, where we also se furher noaion used hroughou he paper. Random recurrence equaions have many applicaions in finance, biology or fracal images, o name jus a few, see e.g. Wong and Li [27], Tong [26], or Diaconis and Freedman [7]. Hence mulivariae generalized Ornsein Uhlenbeck processes as heir coninuous ime counerpars have considerable poenial for applicaions. In one dimension, various applicaions of he GOU process are known. For example, he volailiy of he COGARCH(1,1 process of Klüppelberg e al. [16] or he risk process of Paulsen [20] are one-dimensional GOU processes. As an example of an applicaion of he MGOU process 3

4 o finance, we presen in Example 3.6 he sae vecor process of he volailiy process of he COGARCH(q, p model of Brockwell e al. [5] as a special case of an MGOU process. Furher applicaions of MGOU processes as mulivariae volailiy models seem possible, bu we shall no pursue his opic furher in his paper bu leave i o fuure research. Finally, we menion ha major pars of he resuls of his paper have been obained in he firs named auhor s docoral hesis [2, Chaper 5]. 2 Preliminaries Throughou his paper for any marix M R m n we wrie M for is ranspose and le M (i,j denoe he componen in he ih row and jh column of M. Limis in disribuion will be denoed by d- lim or d, limis in probabiliy by P- lim or P, and almos surely will be abbreviaed by a.s.. The law of a random marix Y will be denoed by L(Y. We wrie N = {1, 2,...}, N 0 = {0, 1, 2,...} and log + (x := log max{x, 1} for x R. Jumps of a marix valued càdlàg process X = (X 0 will be denoed by X := X X wih X := lim s X s for > 0 and he convenion X 0 := 0. Muliplicaive Lévy processes Recall ha an (addiive Lévy process X = (X 0 wih values in R m l is a process wih saionary and independen (addiive incremens which has almos surely càdlàg pahs and sars a 0. Here, an incremen of X is given by X X s for s. We refer o Applebaum [1] or Sao [23] for furher informaion regarding Lévy processes. In he following i will be also necessary o consider muliplicaive Lévy processes wih values in he general linear group GL(R, m of order m, where he group operaion is marix muliplicaion. For ha, remark ha he group srucure allows us o define (muliplicaive lef incremens X Xs and (muliplicaive righ incremens Xs X for 0 s < of a GL(R, m-valued process. We say ha he process (X 0 in GL(R, m has independen lef incremens if for any n N, 0 < 1 <... < n, he random variables X 0, X 1 X0,..., X n X n are independen. The process has saionary lef incremens if X Xs d = X s X0 holds for all s <. Saionariy and independence of righ incremens is undersood analogously. Now following he noaions in he book of Liao [17] a càdlàg process (X 0 in GL(R, m, m 1, wih X 0 = I a.s. is called a (muliplicaive lef Lévy process, if i has independen and saionary righ incremens. Similarly, a càdlàg process (X 0 in GL(R, m, m 1, wih X 0 = I a.s. is called a (muliplicaive righ Lévy process, if i has independen and saionary lef incremens. Given a filraion F = (F 0, a lef Lévy process (X 0 in GL(R, m is called a lef F-Lévy process, if i is adaped o F and for any s < he righ incremen Xs X is independen of F s. Righ F-Lévy processes and (addiive F-Lévy processes are defined similarly. The following lemma gives he connecion beween he random funcionals A s, saisfying Assumpion 1 and muliplicaive Lévy processes. 4

5 Lemma 2.1. (a For any (A s, 0 s fulfilling Assumpion 1 he process (A 0 = (A 0, 0 has a càdlàg modificaion which is a righ Lévy process in GL(R, m. Conversely, if (A 0 is a righ Lévy process in GL(R, m, hen (A s, 0 s defined by A s, = A A s fulfills Assumpion 1. (b For any (A s,, B s, 0 s fulfilling Assumpion 1 he process (A, B 0 = (A 0,, B 0, 0 has a càdlàg modificaion. Proof. (a Since by Assumpion 1(a we have A A s = A s, i follows direcly from Assumpion 1(b and (c, ha (A 0 is a sochasic process in GL(R, m wih saionary and independen lef incremens. I is everywhere coninuous in probabiliy from he righ since by 1(a, (c and (d Similarly due o P- lim h 0 A +h = P- lim h 0 A,+h A = A, 0. P- lim h 0 A h = P- lim h 0 A A h, = A P- lim h 0 A h = A, 0, i is also coninuous in probabiliy from he lef such ha by [25, Theorem V.3] a càdlàg modificaion exiss which is a righ Lévy process in GL(R, m as specified above. The converse is rue by he definiion of righ Lévy processes. (b Since B +h = A +h A B + B,+h he process (B 0 is by Assumpion 1(c and (d everywhere coninuous in probabiliy from he righ and similarly from he lef. Hence i admis a càdlàg modificaion which can be shown by a simple exension of he proof in he one-dimensional case given in [12, Lemma 2.1]. Since every se (A s,, B s, 0 s of random funcionals saisfying Assumpion 1 admis a càdlàg modificaion (A, B 0 by he preceding lemma, we may and do resric aenion o such funcionals wih càdlàg pahs. Marix valued sochasic inegrals Given a filraion F = (F 0 saisfying he usual hypoheses (cf. [21, p. 3], a marixvalued sochasic process M = (M 0 is called an F-semimaringale or simply a semimaringale if every componen (M (i,j 0 is a semimaringale wih respec o he filraion F. For a semimaringale M in R m n and a locally bounded predicable process H in R l m he R l n -valued (lef sochasic inegral I = HdM is given by I (i,j = m k=1 H (i,k dm (k,j and in he same way for M R l m, H R m n, he R l n -valued sochasic (righ inegral J = dmh is given by J (i,j = m k=1 H (k,j dm (i,k. Sochasic inegrals of he form HdM H for locally bounded predicable processes H and H are defined similarly in he obvious way. Given wo semimaringales M and N in R l m and R m n he quadraic variaion [M, N] in R l n is defined by is componens via [M, N] (i,j = m k=1 [M (i,k, N (k,j ]. Similarly is 5

6 coninuous par [M, N] c is given by ([M, N] c (i,j = m k=1 [M (i,k, N (k,j ] c. Wih hese noaions, for wo semimaringales M and N in R m m and wo locally bounded predicable processes G and H in R m m we have he following a.s. equaliies as saed e.g. in Karandikar [15] [ ] G s dm s, dn s H s = G s d[m, N] s H s, 0, (2.1 (0, ] (0, ] [ M, (0, ] G s dn s ] = [ dm s G s, N (0, ] ], 0, (2.2 and he inegraion by pars formula akes he form (MN = M s dn s + dm s N s + [M, N], 0. (2.3 The mulivariae sochasic exponenial Sochasic exponenials of R m m -valued Lévy processes will play a crucial rule in our consideraions. We firs recall he definiion of lef and righ sochasic exponenials from [21, p ]. Definiion 2.2. Le (X 0 be a semimaringale in R m m. Then is lef sochasic exponenial E (X is defined as he unique R m m -valued, adaped, càdlàg soluion of he inegral equaion Z = I + Z s dx s, 0, (2.4 while he unique adaped, càdlàg soluion of Z = I + dx s Z s, 0, (2.5 will be called righ sochasic exponenial and denoed by E (X. Boh E (X and E (X are semimaringales. Unforunaely, unlike for one-dimensional sochasic exponenials as e.g. in [21, Theorem II.37], no closed form expression is available for general mulivariae sochasic exponenials, which makes heir reamen more difficul. The SDE of he sochasic exponenial for processes wih values in arbirary Lie groups has been sudied by Esrade [10]. Remark ha replacing Z and X by heir ransposes in (2.4 leads o he SDE (2.5 and vice versa. Hence we have E (X = E (X. (2.6 As has been observed by Karandikar [15] a necessary and sufficien condiion for nonsingulariy of he lef sochasic exponenial of an R m m -valued process X a ime, is o 6

7 claim ha (I + X s is inverible for all 0 < s. Due o he above saed relaionship beween lef and righ exponenial his resul holds rue also for righ exponenials and hence any sochasic exponenial is inverible for all 0 if and only if de(i + X 0 for all 0. (2.7 For GL(R, m-valued semimaringales, he sochasic logarihm is defined as follows. Definiion 2.3. Le (Z 0 be a GL(R, m-valued semimaringale wih Z 0 = I. Then he lef sochasic logarihm Log Z and righ sochasic logarihm Log Z of Z are defined by Log (Z = Zs dz s, and Log (Z = dz s Zs, 0, (2.8 respecively. I is clear from he defining SDE dz = Z dx for lef sochasic exponenials ha if X is a semimaringale saisfying (2.7 wih X 0 = 0, hen Log E (X = X and X is he unique semimaringale Y saisfying Y 0 = 0 and E (Y = E (X. The same is rue for righ sochasic exponenials and righ sochasic logarihms. The following one-o-one relaion beween muliplicaive Lévy processes and sochasic exponenials of addiive Lévy processes is a key observaion for he invesigaions in his paper. Proposiion 2.4. Le F = (F 0 be a filraion saisfying he usual hypoheses. Then for every F-Lévy process (X 0 in R m m fulfilling (2.7, he sochasic exponenial Z = E (X (resp. Z = E (X is a lef (resp. righ F-Lévy process in GL(R, m. Conversely, if Z = (Z 0 is a lef (resp. righ F-Lévy process in GL(R, m, hen Z is an F-semimaringale and Log Z (resp. Log Z is an addiive Lévy process in R m m saisfying (2.7. Skech of Proof. The firs par follows by simple calculaions using he Markov propery of X, and we refer o [2, Prop. 5.5] for a complee proof. The converse has been observed by Holevo [13] as a conclusion of resuls by Skorokhod [24]. Acually, here i is only observed ha Z is a semimaringale wih respec o is augmened naural filraion, H say, and ha Log Z and Log Z, resp., are H-Lévy processes, bu i is easy o see ha hen Log Z and Log Z are even F-Lévy processes, and since E ( Log Z = Z and E ( Log Z = Z, resp., i follows ha Z is an F-semimaringale. Again we refer o [2, Prop. 5.5] for deailed calculaions. Since he inverse and he ranspose of a lef Lévy process in GL(R, m are righ Lévy processes and vice versa, for any addiive Lévy process (X 0 fulfilling (2.7 he process ( E (X 0 is a righ Lévy process and hence by he above proposiion i is he righ 7

8 sochasic exponenial of anoher Lévy process (U 0. In fac (see [15, Theorem 1] if (X 0 is a semimaringale such ha (2.7 is fulfilled, hen i holds wih E (X = [ E (U ] = E (U, 0 U := X + [X, X] c + 0<s ( (I + Xs I + X s, 0. (2.9 Remark ha i follows from (2.9 by sandard calculaions ha he processes U and X fulfill he relaion U = X [X, U], 0, (2.10 and ha if X is a Lévy process, hen so is U and vice versa. 3 Mulivariae Generalized Ornsein-Uhlenbeck Processes In his secion we will characerize all families of random funcionals (A s,, B s, 0 s saisfying Assumpion 1 and hen will use his o define mulivariae generalized Ornsein- Uhlenbeck processes. Furher, we show ha every mulivariae generalized Ornsein- Uhlenbeck process (V 0 is a soluion of he SDE dv = du V + dl for a suiable R m m R m -valued Lévy process (U, L. Conversely, provided ha V 0 is F 0 -measurable for some filraion F = (F 0 saisfying he usual hypoheses such ha he Lévy process (U, L is a semimaringale wih respec o F, he soluion o his SDE is a mulivariae generalized Ornsein Uhlenbeck process. The proofs for he resuls of his secion are given in Secion 6. The following heorem characerizes all choices of random funcionals (A s,, B s, 0 s fulfilling Assumpion 1. Recall ha A = A 0,, B = B 0, and ha by Lemma 2.1 we can resric o càdlàg versions of (A, B 0. Theorem 3.1. Suppose ha (A s,, B s, 0 s saisfies Assumpion 1 and ha (A 0 and (B 0 are chosen o be càdlàg. Then here is a unique Lévy process (X, Y in R m m R m such ha X saisfies (2.7 and such ha ( As, = E (X E (Xs B a.s., 0 s. (3.1 s, E (X (s,] E (Xu dy u The Lévy process (X, Y is given by ( ( X = Log A Y A u d(a u B u, 0, (3.2 8

9 where he inegral is defined as a sochasic inegral wih respec o he naural augmened filraion of (A, B 0, for which (A 0 and (B 0 are semimaringales. Conversely, if (X, Y is a Lévy process in R m m R m such ha X saisfies (2.7, hen (A s,, B s, 0 s defined by he righ hand side of (3.1 saisfies Assumpion 1. Since a mulivariae generalized Ornsein-Uhlenbeck process (V 0 was supposed o saisfy (1.2 wih (A s,, B s, 0 s saisfying Assumpion 1, Theorem 3.1 moivaes he following definiion. Definiion 3.2. Le (X, Y = (X, Y 0 be a Lévy process in R m m R m such ha X saisfies (2.7 and le V 0 be a random variable in R m. Then he R m -valued process (V 0, given by ( V := E (X V 0 + E (Xs dy s, 0, (3.3 will be called mulivariae generalized Ornsein-Uhlenbeck (MGOU process driven by (X, Y 0. The MGOU process will be called causal or non-anicipaive, if V 0 is independen of (X, Y, and sricly non-causal if V is independen of (X s, Y s 0 s< for all 0. I is easy o see ha an MGOU process indeed saisfies (1.2. Remark ha even for m = 1 Definiion 3.2 is generalizing he sandard definiion of a generalized Ornsein-Uhlenbeck process since we do no assume a priori ha V 0 is independen of (X, Y 0 and also he condiion of E(X o be sricly posiive is dropped. Neverheless i seems naural o us o include hese cases in he class of generalized Ornsein-Uhlenbeck processes. Observe ha any MGOU process wih saring random variable V 0 independen of (X, Y is a ime-homogeneous Markov process. Example 3.3. (a If X = Λ for some Λ R m m is a pure drif process hen E (X = E (X = e Λ and he MGOU process V = e (V Λ 0 + e Λs dy s, 0, driven by (X, Y is he usual mulivariae Ornsein Uhlenbeck ype process driven by Y as inroduced in [22]. (b If (X, Y = (diag(x (1,1,..., X (m,m, (Y (1,..., Y (m, i.e. if X is a Lévy process concenraed on he diagonal marices, and X saisfies condiion (2.7, hen E (X = E (X = diag(e(x (1,1,..., E(X (m,m, where E( denoes he usual one-dimensional sochasic exponenial, and he ih componen V (i of he MGOU process (V 0 driven by (X, Y saisfies V (i = E(X (i,i ( V (i 0 + E(X (i,i s dy s (i 0, 0, i = 1,..., m. 9

10 I follows ha V (i is a one-dimensional MGOU process driven by (X (i,i, Y (i. If addiionally X (i,i does no have jumps of size less han or equal o and if V (i 0 is independen of (X (i,i, Y (i, hen V (i is a GOU process. Observe ha in general componens of MGOU processes are no MGOU processes if X is no concenraed on he diagonal marices. An MGOU process can also be characerized by he sochasic differenial equaion i saisfies. Theorem 3.4. (a Le (X, Y be a Lévy process in R m m R m such ha (2.7 holds, and le (V 0 be he MGOU process driven by (X, Y wih saring random variable V 0. Le F = (F 0 be some filraion saisfying he usual hypoheses such ha (X, Y is a semimaringale wih respec o F and V 0 is F 0 -measurable. Then (V 0 solves he SDE dv = du V + dl, (3.4 where (U, L is he Lévy process in R m m R m wih U as defined in (2.9 and L given by L = Y + ( (I + Xs I Y s [X, Y ] c, 0. (3.5 The process U saisfies 0<s de(i + U 0 for all 0. (3.6 (b Conversely, if (U, L is a Lévy process in R m m R m such ha U saisfies (3.6, F = (F 0 is a filraion saisfying he usual hypoheses such ha (U, L is an F-semimaringale and V 0 is an R m -valued F 0 -measurable saring random variable, hen he soluion o (3.4 is an MGOU process driven by (X, Y, where (X, Y is he Lévy process defined by ( ( X Log ( E (U = Y L + [ Log (, 0, (3.7 E (U, L] and X saisfies (2.7. Observe ha under he naural assumpion ha V 0 is independen of (X, Y (i.e. for a causal MGOU process, he smalles filraion F which saisfies he usual hypoheses and is such ha V 0 is F 0 measurable and (X, Y is adaped o F is a filraion such ha X, Y, U and L are semimaringales wih respec o i (cf. Corollary 1 of Theorem VI.11 in [21], as required in he saemen of (a. A similar remark holds for (b if V 0 is independen of (U, L. In he following proposiion we sae some cross-relaions beween (X, Y and (U, L defined by (2.9 and (

11 Proposiion 3.5. Le (X, Y be a Lévy process in R m m R m such ha X saisfies (2.7 and le (U, L be defined by (2.9 and (3.5. Then L = Y + [U, Y ], 0, (3.8 and Y = L + [X, L], 0. (3.9 Finally, we show in he nex example ha he sae vecor of he COGARCH(q, m volailiy process is an m-dimensional MGOU process. Example 3.6. Le m, q N, q m, c 1,..., c m, d 0,..., d m R wih c m 0 and d q 0, d q =... = d m = 0. Denoe d d 1 C = , e = 0., d = d m 2 c m c m c m 2 c 1 1 d m wih C R m m, e, d R m, and le M be a one-dimensional Lévy process wih non-rival Lévy measure. Le β > 0. Then, as defined in [5], he COGARCH(q, m process, driven by M and wih parameers C, β and d has (righ-coninuous volailiy process (S 0 given by S = β + d V, 0, (3.10 where he sae vecor process V = (V 0 is he unique càdlàg soluion of he sochasic differenial equaion dv = CV d + es d[m, M] (d = CV d + e(β + d V d[m, M] (d, 0, (3.11 wih iniial value V 0, independen of (M 0. Here, [M, M] (d = 0<s ( M s 2 denoes he discree par of he quadraic variaion of M. If he process (S 0 is non-negaive almos surely, condiions for which are given in Secion 5 of [5], hen G = (G 0, defined by G 0 = 0, dg = S dm, is called a COGARCH(q, m process wih parameers C, d, β and driving Lévy process M. I follows from [5, Theorem 3.3] and is proof ha he sae vecor process (V 0 saisfies (1.2 wih random funcionals (A s,, B s, which saisfy Assumpion 1, so ha (V 0 is an MGOU process. Using he SDE (3.11 and Theorem 3.4, we ge anoher proof of his, observing ha dv = CV d + βed V d[m, M] (d + βe d[m, M] (d = (C d + βed d[m, M] (d V + βe d[m, M] (d = du V + dl, 11

12 where U = C + β[m, M] (d ed and L = β[m, M] (d e. (3.12 Since he jumps of [M, M] (d are non-negaive, i follows ha U saisfies condiion (3.6 and hence ha V is a causal MGOU process by Theorem MGOU Processes Carried by Affine Subspaces In his secion we will classify MGOU processes which are carried by affine subspaces of R m. To do ha, we inroduce he noion of irreducibiliy which we mainly adop from Bougerol and Picard [4] who sudied generalized auoregressive models in discree ime. The proofs for he resuls of his secion are given in Secion 7. Definiion 4.1. Suppose (X, Y 0 is a Lévy process in R m m R m such ha X saisfies (2.7 and define (A s,, B s, 0 s by (3.1. Then an affine subspace H of R m is called invarian under he auoregressive model (1.2 if A s, H + B s, H, almos surely, holds for all 0 s. If R m is he only invarian affine subspace, he model (1.2 is called irreducible. Obviously, by Assumpion 1(c, i is enough o require he above condiion for s = 0 and all 0. Remark ha he given definiion of invarian subspaces is more resricive han he one in [4], since e.g. seing Y = B = 0 and leing A be a roaion operaor wih angle 2π implies ha in he discree ime model V n = A n,n V n + B n,n, n N, every poin is a zero-dimensional invarian affine subspace, while only he roaion axis is invarian for all 0. Accordingly, irreducibiliy of he coninuous ime model does no direcly imply ha for all h > 0 he discree ime model V nh = A (nh,nh V (nh +B (nh,nh, n N, is irreducible in he sense of [4]. Bu we can show he following proposiion which saes ha a leas here is some h > 0 for which he corresponding discree ime model is irreducible. This will be an imporan ingredien when proving Theorems 5.3, 5.4 and 5.7 below on he exisence of sricly saionary soluions. Proposiion 4.2. Suppose (X, Y 0 is a Lévy process in R m m R m such ha X saisfies (2.7 and define (A s,, B s, 0 s by (3.1. Suppose ha he auoregressive model (1.2 is irreducible. Then here exiss h > 0 for which he discree-ime auoregressive model V nh = A (nh,nh V (nh + B (nh,nh, n N, (4.1 is irreducible in he sense ha here exiss no affine subspace H of R m, H R m, such ha for all n N, A (nh,nh H + B (nh,nh H almos surely. The nex heorem reas MGOU processes where he corresponding auoregressive model admis a d-dimensional invarian affine subspace H. I urns ou ha in his case we can spli up he process carried by H in a consan par and an R m d -valued MGOU process. For convenience we firs assume ha H is parallel o he axes. 12

13 Theorem 4.3. Suppose (V 0 is an MGOU process wih saring random variable V 0, driven by he Lévy process (X, Y 0 in R m m R m, where X fulfills (2.7, and le (A s,, B s, 0 s as defined in (3.1. (a Assume ha H = {(k 1,..., k d, h d+1,..., h m, h d+1,..., h m R} wih 1 d m and consans k 1,..., k d R is an invarian, affine subspace of( R m wih respec o K he model (1.2. Then, given ha V 0 H a.s., i holds V = H a.s. for each 0 wih K = (k 1,..., k d and V R m d, and he Lévy processes X and Y saisfy for all 0 ( X 1 X = 0 X 2 X 3 a.s. where X 1 R d d and (4.2 ( ( Y 1 Y = X 1 = K a.s. where Y 1 R d. (4.3 Y 2 Y 2 The process (V 0 is an MGOU process driven by he Lévy process ( X 3, Y 2 X 2 K 0 (4.4 V in R (m d (m d R m d. Furher, if (U, L is defined as in (2.9 and (3.5, and if V 0 is F 0 -measurable for a filraion F = (F 0 saisfying he usual assumpions such ha U and L are semimaringales wih respec o F (hence (V 0 solves he SDE (3.4 by Theorem 3.4, hen we have a.s. for each 0 ( ( U 1 U = 0 L 1 U 2 U 3 and L = L 2 wih U 1 R d d, L 1 R d, (4.5 where L 1 = U 1 K a.s. and (V 0 solves he SDE dv = du 3 V + d(l 2 + U 2 K, 0. (4.6 (b Conversely, if (4.2 and (4.3 hold for K = (k 1,..., k d R consan, hen he affine subspace H = {(k 1,..., k d, h d+1,..., h m, h d+1,..., h m R} of R m is invarian wih respec o he model (1.2 and for any saring random ( variable V 0 H he K MGOU process defined by (3.3 can be wrien as V = a.s., where (V 0 is an MGOU process driven by he Lévy process (4.4. V Remark 4.4. Observe ha if in he seing of Theorem 4.3 he invarian affine subspace H is no parallel o he axes, hen here exiss an orhogonal ransformaion marix O, such ha OH fulfills he assumpions of Theorem 4.3 for he ransformed MGOU process V = OV. The process (V 0 fulfills he random recurrence equaion V = A s,v s + B s, for 0 s where A s, = OA s, O and B s, = OB s, and hence by Theorem 3.1 i is an MGOU process driven by (OX O, OY 0. Thus he sudy of arbirary invarian affine subspaces reduces o he case reaed in Theorem

14 This observaion and Theorem 4.3 imply he following characerizaion of irreducibiliy of he model (1.2. Corollary 4.5. Suppose (X, Y 0 in R m m R m is a Lévy process such ha X fulfills (2.7. Then he auoregressive model (1.2 wih (A s,, B s, 0 s as defined in (3.1 is irreducible if and only if here exiss no pair (O, K of an orhogonal ransformaion O R m m and a consan K = (k 1,..., k d R d, 1 d m, such ha a.s. ( ( X OX O 1 = 0 X 1 and OY = K where X 1 R d d, 0. (4.7 X 2 X 3 Wih (U, L 0 as defined in (2.9 and (3.5, Equaion (4.7 is furher equivalen o ( ( U OU O 1 = 0 U 1 and OL = K a.s. wih U 1 R d d. (4.8 U 2 U 3 Y 2 L 2 5 Saionary Soluions of MGOU Processes In his secion we invesigae condiions for he exisence of sricly saionary soluions of mulivariae generalized Ornsein-Uhlenbeck processes. The proofs of he resuls are given in Secion 8. Given some exra informaion on he limi behaviour of E (X our firs heorem provides necessary and sufficien condiions for he exisence of saionary soluions of MGOU processes. Before we sae i we give he following lemma on sochasic exponenials which is ineresing in is own righ. Lemma 5.1. Le (X 0 be a Lévy process in R m m. Then for any 0 fixed we have ha E (X d = E (X. In paricular his implies P- lim E (X = 0 P- lim E (X = 0. (5.1 Since E (U = E (X, he condiion P- lim E (U = 0 appearing in Theorem 5.2(a below is equivalen o P- lim E (X = 0. Hence, Theorem 5.2 gives necessary and sufficien condiions for saionariy if eiher P- lim E (X = 0 or P- lim E (X = 0 and hus exends [3, Theorem 2.1]. Theorem 5.2. Suppose (V 0 is an MGOU process driven by he Lévy process (X, Y 0 in R m m R m such ha X saisfies (2.7. Le (U, L 0 be he Lévy process defined in (2.9 and (

15 (a Suppose lim E (U = 0 in probabiliy. Then a finie random variable V 0 can be chosen such ha (V 0 is sricly saionary if and only if he inegral E (Us dl s converges in disribuion for o a finie random variable. In his case, he disribuion of he sricly saionary process (V 0 is uniquely deermined and is obained by choosing V 0 independen of (X, Y 0 wih V 0 d = d- lim E (Us dl s. (b Suppose lim E (X = 0 in probabiliy. Then a finie random variable V 0 can be chosen such ha (V 0 is sricly saionary if and only if he inegral E (Xs dy s converges in probabiliy o a finie random variable as. In his case he sricly saionary soluion is unique and given by V = E (X E (Xs dy s a.s. for all 0. (, Observe ha he soluion obained in Theorem 5.2(a is causal and ha he one in (b is sricly non-causal. By adding he assumpion of irreducibiliy of he underlying model, as characerized in Corollary 4.5, he above heorem can be sharpened as follows. Theorem 5.3. Suppose (X, Y 0 is a Lévy process in R m m R m such ha X saisfies (2.7 and such ha he corresponding auoregressive model (1.2 wih (A s,, B s, 0 s as defined in (3.1 is irreducible. Le (V 0 be he MGOU process driven by (X, Y 0 and le (U, L 0 be he Lévy process defined in (2.9 and (3.5. (a A finie random variable V 0, independen of (X, Y 0, can be chosen such ha (V 0 is sricly saionary if and only if lim E (U = 0 in probabiliy and he inegral E (Us dl s converges in disribuion for o a finie random variable. (b A finie random variable V 0 can be chosen such ha (V 0 is sricly saionary and sricly non-causal if and only if lim E (X = 0 in probabiliy and he inegral E (Xs dy s converges in probabiliy as. Given ha he processes U and L have a finie log-momen Theorem 5.3 can be sharpenend o obain a necessary and sufficien condiion for he exisence of sricly saionary soluions of MGOU processes in erms of he driving Lévy process as saed in Theorem 5.4. To explain is condiions (iv and (v and relae i o he corresponding discree ime resuls, le be a fixed, submuliplicaive marix norm. Recall ha he op Lyapunov exponen of an R m m -valued i.i.d. sequence (C n n N wih E[log + C 1 ] is given by 1 γ := inf n N n E[log C 1 C n ]. (5.2 I is independen of he specific submuliplicaive marix norm used and i holds 1 γ = lim n n log C 1 1 C n = lim n n log C n C 1 a.s., (5.3 15

16 cf. Fursenberg and Kesen [11] and Bougerol and Picard [4]. In [4, Theorem 2.5] i is also shown ha if he discree ime model W n = C n W n + D n, n Z, is irreducible, where (C n, D n n Z is an i.i.d. R m m R m -valued sequence wih E[log + C 1 ] < and E[log + D 1 ] <, hen he discree model admis a sricly saionary causal soluion if and only if he op Lyapunov exponen of he sequence (C n n N is sricly negaive. Now if X and U are as in Theorem 5.4 (v, hen E[log + U 1 ] < implies E[log + E (U ] < for every > 0 as will be shown in Proposiion 8.4 below. Since for each h > 0 he sequence (A (nh,nh = E (U nh E (U (nh n N is i.i.d. and E (U nh = A (nh,nh A 0,h, i follows ha here is h > 0 such ha he op Lyapunov exponen of (A (nh,nh n N is sricly negaive if and only if here is 0 > 0 such ha E[log + E (U 0 ] <, which is equivalen o condiion (iv below by Lemma 5.1. Theorem 5.4. Suppose (X, Y 0 is a Lévy process in R m m R m such ha X saisfies (2.7 and ha he corresponding auoregressive model (1.2 wih (A s,, B s, 0 s as defined in (3.1 is irreducible. Le (V 0 be he MGOU process driven by (X, Y 0 and le (U, L 0 be he Lévy process defined in (2.9 and (3.5. Suppose ha E[log + U 1 ] < and E[log + L 1 ] <. Then he following are equivalen: (i A finie random variable V 0, independen of (X, Y 0, can be chosen such ha (V 0 is sricly saionary. (ii I holds lim E (U = 0 in probabiliy and he inegral E (Us dl s converges in disribuion for o a finie random variable. (iii I holds lim E (U = 0 a.s. and he inegral E (Us dl s converges a.s. for o a finie random variable. (iv There exiss 0 > 0 such ha E[log E (U 0 ] < 0. If addiionally U is a compound Poisson process wih jump heighs (S k k N, hen he above condiions (i o (iv are furher equivalen o (v The op Lyapunov exponen of he sequence (I + S k k N is sricly negaive. Remark 5.5. (a A similar resul as Theorem 5.4 also holds rue for sricly non-causal sricly saionary soluions of MGOU processes in he irreducible case. (b The proof of Theorem 5.4 given in Secion 8 shows ha he implicaions (iv = (iii = (ii = (i and (v = (iii = (ii = (i also hold wihou assuming irreducibiliy of he underlying model. Example 5.6. Consider he sae vecor process (V 0 of he COGARCH(q, m-volailiy process (S 0 as defined in Example 3.6 wih dv = du V + dl and (U, L 0 given by (3.12. Suppose ha m = 2. Then i follows from Corollary 4.5 by a sraighforward bu edious calculaion, using ha c 2 0, ha he corresponding auoregressive 16

17 model (1.2 is irreducible. In paricular, by Theorem 5.3(a, a sricly saionary (causal COGARCH(q, 2-volailiy sae vecor process exiss if and only if lim E (U = 0 in probabiliy and E (Us dl s = β 0<s E (Us e( M s 2 converges in disribuion o a finie random variable as. If in addiion log x ν x >1 M(dx <, where ν M denoes he Lévy measure of M, hen E log + U 1 < and E log + L 1 <, and by Theorem 5.4 he above condiions are equivalen o E log + E (U 0 < 0 for some 0 > 0. Tha he laer condiion is sufficien for a (causal sricly saionary sae vecor o exis was already observed in Remark 3.4(a of [5], bu having he irreducibiliy of he model we now also know ha i is necessary under he finie log-momen assumpion on ν M. Observe however ha he volailiy process (S 0 defined in (3.10 may be sricly saionary even wihou (V 0 being sricly saionary, since i is only a specific linear combinaion of (V 0 plus a consan. We shall no pursue he issue of sric saionariy of (S 0 furher. Also, we have no invesigaed if he auoregressive model (1.2 for he COGARCH(q, m volailiy process wih m 3 is always irreducible. In he case ha he underlying model is no irreducible, P- lim E (U = 0 is no necessary for he exisence of a causal sricly saionary soluion as shown in he following heorem. Theorem 5.7. Suppose (X, Y 0 is a Lévy process in R m m R m such ha X saisfies (2.7 and le (V 0 be he MGOU process driven by (X, Y 0 saisfying he auoregressive model (1.2 wih (A s,, B s, 0 s as defined in (3.1. Define (U, L 0 via (2.9 and (3.5. Then a finie random variable V 0 can be chosen such ha (V 0 is sricly saionary and causal if and only if here exiss a pair (O, K of an orhogonal ransformaion O R m m and a consan K = (k 1,..., k d, 0 d m such ha (4.7 and hence (4.8 hold and such ha P- lim E (U 3 = 0 and E (U 3 s d(l 2 s + U 2 sk converges in disribuion o a finie random variable as. If hese condiions are saisfied a sricly saionary soluion can be obained by choosing V 0 independen of (X, Y 0 wih he same disribuion as he disribuional limi as of ( K O. E (U 3 s d(l 2 s + U 2 sk If d = 0 in he above condiions hen L 2 s + U 2 sk has o be inerpreed as L 2 s, and if d = m hen U 3 is zero-dimensional and he convergence condiions regarding E (U 3 and E (U 3 s d(l 2 s + U 2 sk do no appear. Remark 5.8. Using argumens as in he proof of Theorem 5.3(b a similar resul as Theorem 5.7 for sricly noncausal sricly saionary soluions of MGOU processes can be obained, oo. 17

18 Remark 5.9. The resuls in Secions 3 and 5 remain valid if we rea an MGOU process (V 0 wih V R m l and drop he condiion of l = 1. As he value of l has no influence on he proofs we can simply replace he vecor valued processes (Y 0 and (L 0 by R m l - valued processes. Theorem 4.3 may be applied column-by-column or, alernaively, i is possible o inerpre he MGOU process (V 0 in R m l driven by (X, Y 0, X R m m, Y = (Y 1,..., Y l R m l as an MGOU process in R ml driven by he Lévy process X 0 Y 1...,. in R ml ml R ml. 0 X 6 Proofs for Secion 3 Y l Before proving Theorem 3.1 we give he following proposiion which esablishes in paricular he semimaringale propery of (A 0 and (B 0. Proposiion 6.1. Suppose (A s,, B s, 0 s is a process saisfying Assumpion 1 and such ha (A, B 0 is càdlàg. Le H be he naural augmened filraion of (A, B 0. Then (A 0 and (B 0 are H-semimaringales. Furher, he R m m R m R m m R m - valued process (U, L, X, Y 0 defined by U L X Y is an H-Lévy process. = 0 Log A = da s A s B da s A s B s Log A = A s da s A s d(a s B s, 0, (6.4 Proof. Observe ha (A 0 is a righ H-Lévy process by Assumpion 1 and Lemma 2.1 and hence an H-semimaringale by Proposiion 2.4. I follows ha (U, L, X 0 as given in (6.4 is well defined. By compuaions similar o hose in he proof of [12, Theorem 2.2] one can show ha for 0 s U U s L L s X X s = (s,] d(a s, u A s,u B s, (s,] d(a s, u A s,u B s,u (s,] A s,u d(a s, u, 0 s. (6.5 By Assumpion 1(b,c we observe ha (A s,s+u, B s,s+u d u 0 = (A 0,u, B 0,u u 0 and hus we obain from (6.5 ha (U, L, X has saionary incremens. By Assumpion 1(b, (U U s, L L s, X X s is independen from H s for 0 s, where H = (H 0. We also know ha (U 0, L 0, X 0 = 0 a.s., ha he pahs of (U, L, X are càdlàg since ha held rue for (A, B 0, and ha clearly (U, L, X is adaped o H. Hence (U, L, X 0 is an H-Lévy process. In paricular, L is an H-semimaringale, so ha by (6.4, (B 0 = (L + da s A s B s 0 18

19 is an H-semimaringale, oo. Consequenly Y as given in (6.4 is well defined. For 0 s we hen have from Assumpion 1(a ha Y Y s = A u d(a u B u (s,] = A s,u A s d(a s A s, A s, B s + A s A s, B s, u (s,] = A s,u d(a s, B s, u. (s,] I now follows in complee analogy o he reasoning given above ha (U, L, X, Y is an H-Lévy process. Proof of Theorem 3.1. Le (A s,, B s, 0 s saisfy Assumpion 1. By Proposiion 6.1, (A 0 and (B 0 are semimaringales wih respec o heir naural augmened filraion, and (X, Y defined by (3.2 is a Lévy process. Clearly, X saisfies (2.7, and for 0 s i holds A s, = A A s = (A A s Furhermore, A dy = d(a B from (3.2, so ha B = A A u dy u = E (X giving B s, = B A s, B s = E (X = E (X (s,] E (Xu dy u E (X E (Xu dy u. = E (X E (Xs. E (Xu dy u, E (Xs E (X s This is (3.1. The uniqueness of (X, Y is clear from (3.1. (0,s] E (Xu dy u For he converse, le (X, Y a Lévy process in R m m R m such ha X saisfies (2.7. Le F be he augmened naural filraion of (X, Y 0, hen (A s,, B s, 0 s as given in (3.1 is well defined wih respec o F and we know from Proposiion 2.4 ha E (X is a righ F-Lévy process in GL(R, m whose lef incremens are given by A s,. Thus we have ha for all 0 s u almos surely A s, = A u, A s,u holds. Also i follows direcly from he definiions of A s, and B s, ha B s, = A u, B s,u + B u, a.s. such ha Assumpion 1(a is fulfilled. For he common process (A s,, B s, 0 s observe ha for 0 s we have ( As, = E (X E (Xs B. s, E (X E (Xs (s,] E (X s E (Xu dy u 19

20 Since (A 0 is a righ F-Lévy process he common incremens ( E (X E (Xs, Y Y s s are independen of (X u, Y u 0 u s. Hence i follows ha {(A s,, B s,, a s b} and {(A s,, B s,, c s d} wih b c are independen. Similarly we conclude ha ( As, B s, d = which yields Assumpion 1(c. E (X s E (X0 E (X 0 E (X s E (X0 (0, s] E (Xu dy u = ( A0, s B 0, s The coninuiy in probabiliy a 0 of A = A 0, is clear, while for B = B 0, i follows from ha of A and Y and he coninuiy of he inegral. Proof of Theorem 3.4. (a I is easy o see ha (U, L as consruced in (2.9 and (3.5 is a Lévy process and ha U saisfies (3.6. Define A = A 0, and B = B 0, for 0 by he righ hand side of (3.1. Then V = A V 0 + B. By he definiion of U, we furher have ha da = du A. Hence, denoing L := B da u A u B u as in (6.4, we obain dv = da V 0 + db = du A V 0 + db = du (A V 0 + B + db du B = du V + db da A B = du V + dl. I remains o show ha L = L. Using he inegraion by pars formula (2.3 and (2.2, we obain L = B da s A s B s = E (X E (Xs dy s d( E (X s E (Xu dy u (0,s = E (X E (Xs dy s E (X E (Xs dy s ( ] E (Xu dy u E (Xu dy u + E (X s d (0,s] [ = E (X s E (Xs dy u + [ ] = Y + da s A s, Y (0, ] = Y + [U, Y ]. (0, ] + [ E (X, (0, ] d( E (X s E (X s, Y Tha L = L hen follows from he definiion of U in (2.9 since [ ] ( [U, Y ] = [X, Y ] + [[X, X] c, Y ] + (I + Xs I + X s, Y = [X, Y ] c + 0<s 0<s ( (I + Xs I Y s. ] 20

21 (b Is is clear ha (X, Y as defined in (3.7 is a Lévy process wih X saisfying (2.7. Observe ha he given definiion of (X 0 is equivalen o (2.9 and ha from he definiion of (Y 0 we deduce Y = L + [X, L] = L + ( X s L s + [X, L] c 0<s = L 0<s ( (I + Xs I (I + X s L s + [X, L] c, 0. Hence Y = (I + X L and [X, Y ] c = [X, L] c and we conclude ha Y = L ( (I + Xs I Y s + [X, Y ] c, 0, 0<s ( E ( which is equivalen o (3.5. Thus he MGOU process (X V 0 + E (Xs dy s 0 solves he SDE (3.4 by par (a, giving he claim by he uniqueness of he soluion o (3.4. Proof of Proposiion 3.5. Equaion (3.8 has been esablished when showing ha L = L in he proof of Theorem 3.4(a. Equaion (3.9 follows from he fac ha by (3.5, L = (I + X Y and [X, L] c = [X, Y ] c, so ha by he same calculaion as in he proof of Theorem 3.4(b, L + [X, L] = L ( (I + Xs I Y s + [X, Y ] c, 0. 0<s By (3.5 his implies ( Proofs for Secion 4 In his secion we give he proofs of he resuls of Secion 4. Proof of Proposiion 4.2. For ν N 0 le C ν := A 0,2 ν, D ν := B 0,2 ν and consider he random affine ransformaion f ν : R m R m, x C ν x + D ν, ν N 0. For 0 d < m, a d-dimensional affine subspace H of R m will be called an affine d-fla, and i is f ν -invarian if f ν (H = C ν H + D ν H a.s., which by Assumpion 1(c is equivalen o saying ha H is invarian for he discree ime model V nh = A (nh,nh V (nh + B (nh,nh, n N 0, as defined in (4.1 wih h = 2 ν. Since by (1.3 any subspace which is invarian for he model (4.1 for some h > 0 is also invarian for he model (4.1 for every h := kh wih k N, i is clear ha any f ν invarian affine d-fla is also f ν -invarian. Hence, denoing he se of all f ν -invarian affine d-flas by Hν, d 0 d < m, i follows ha Hν+1 d Hν d 21

22 for all ν N 0, 0 d < m, such ha H d := lim n Hν d = ν=0 Hd ν can be defined for 0 d < m. We furher denoe H ν := m d=0 H d ν for ν N 0 and H := m d=0 H d = H ν. The proof of he proposiion will be given in wo seps: firs i will be shown ha irreducibiliy of he coninuous ime model (1.2 implies H =, and in a second sep ha H = implies he exisence of some ν 0 N 0 such ha H ν0 =, i.e. ha he discree ime model (4.1 is irreducible for h = 2 ν 0. The firs sep will be shown by conradicion, i.e. we assume ha H, i.e. ha here exiss an affine subspace H R m of R m which is invarian under f ν for all ν N. Thus, as argued above, H is also invarian under he model (4.1 for all h = k2 ν, k N, ν N. I hen remains o show ha A 0, H + B 0, H holds for all > 0, i.e. ha H is invarian under (1.2. Bu his follows easily from he fac ha (A 0, 0 and (B 0, 0 have almos surely càdlàg pahs, so ha for every number > 0 we can find a sequence ( n n N of he form n = k n 2 νn converging from he righ o such ha almos surely, A 0, H + B 0, = lim n (A 0,n H + B 0,n H. Hence H is an invarian affine subspace of he coninuous ime model (1.2 giving he desired conradicion. I remains o show ha H = m d=0 Hd = implies he exisence of some ν 0 such ha H ν0 =. Since any affine 0-fla H of he form H = {x} is f ν -invarian if and only if (f ν I(x = 0 a.s., and since f n I is an affine linear mapping, is kernel is an affine linear subspace of R m, S ν say, and we have Hν 0 = {{x}, x S ν }. Since S ν+1 S ν, i follows ha here is ν 1 N such ha S ν1 +n = S ν1 for all n N 0. Hence, H = implies ha here is ν 1 N such ha Hν 0 = for all ν ν 1, and in he following we can concenrae on invarian affine d-flas wih 0 < d < m. Fix a family (O, K = {(O H, K H, H affine d-fla wih 0 < d < m} of pairs of an orhogonal ransformaion O H R m m and a consan K H R m such ha H := O H H K H = {0} m d R d. Then given ν N and some affine d-fla H we obain by easy compuaions ha H is invarian under f ν : x C ν x + D ν if and only if he subspace H is invarian under he mapping g ν : R m R m, x C ν x + D ν wih C ν = O H C ν O H and D ν = O H D ν + (O H C ν O H IK H. Using he special srucure of H his yields ha H is invarian under g ν if and only if i holds almos surely C (i,j ν = 0 and D (i ν ν=0 = 0 for all (i, j J H := {(i, j, 1 i d, d < j m}. This is by definiion of C ν and D ν equivalen o sae ha, almos surely, m m C(k,l ν (O O(j,l C(k,l ν = 0 and (7.1 m k=1 k,l=1 O (i,k H O (i,k H D(k ν + ( m m q=1 k,l=1 H (l,j = O (i,k H k,l=1 O(q,l H C(k,l ν O (i,k H H K (q H K(i H = 0 for all (i, j J H. (7.2 22

23 By inroducing he marices M H,i,j and N H,i in R m m via M (k,l H,i,j := O(i,k H O(j,l H and N (k,l H,i := m q=1 O (i,k H O(q,l H K(q H, k, l = 1,..., m, denoing he ih row of he marix O H by O (i, H and leing vec ( : Rm m R m2 be he vecorizaion operaor which sacks he colums of a given marix below one anoher, (7.1 and (7.2 urn ou o be equivalen o vec (M H,i,j vec (C ν vec (N H,i vec (C ν 0, D ν = 0 = (O (i, H, D ν (i, j J H, 0 K (i H almos surely, where, denoes he sandard scalar produc in R m2 +m+1. Now se R ν := span H H ν \H 0 ν vec (M H,i,j vec (N H,i 0, (O (i, H ; (i, j J H 0 K (i +m+1 Rm2, H where span denoes he linear span. Then by he above we have esablished ha R ν is orhogonal o (vec (C ν, Dν,, a.s., and ha, given an affine d-fla H, 0 < d < m, i is invarian under f ν if and only if all corresponding vecors (vec (M H,i,j, 0, 0 and (vec (N H,i, O (i, H, K(i H for (i, j J H are in R ν. Finally, as R ν is a vecor space and we have ha R ν+1 R ν we observe ha is limi for ν can only be empy (which is equivalen o m d=1 Hd = if here exiss some ν 2 N such ha for all ν ν 2 we have R ν =. Hence i holds H ν = for all ν ν 0 := max{ν 1, ν 2 } so ha he discree ime model (4.1 is irreducible for all ν ν 0 as had o be shown. Proof of Theorem 4.3. (a We sar by verifying (4.2 and (4.3. Since H is an invarian affine subspace we deduce from (1.2 ha for any 0 and all h d+1,..., h m R he equaion A (k 1,..., k d, h d+1,..., h m + B = (k 1,..., k d, g d+1,..., g m a.s. has o admi a soluion g d+1,..., g m R. This is equivalen o d j=1 d j=1 k j A (i,j + k j A (i,j + m j=d+1 m j=d+1 h j A (i,j + b i = k i, i = 1,..., d h j A (i,j + b i = g i, i = d + 1,..., m. Thus we can conclude ha A (i,j = 0 holds a.s. for i d, j > d. Observe by simple algebraic calculaions ha if wo marices M and N in R m m have a d (m d block of 23

24 zero enries in he upper righ corner, hen so do M and MN. More deailed we have for ( ( M1 0 N1 0 M = GL(R, m and N = R m m, M M 2 M 3 N 2 N 1, N 1 R d d, 3 ha M 1 and M 3 are non-singular and i holds ( ( M M 1 0 M = M 3 M 2 M 1 M and MN = 1 N M 2 N 1 + M 3 N 2 M 3 N 3 Now recall ha A = E (X and hus we know ha E (X and E (X a.s. admi a d (m d zero block for all 0. Hence i follows from (2.8 ha also X a.s. has such a zero block which is (4.2. Thus we deduce from (2.4 ha ( E 1 E (X =: 0 E 2 E 3 = ( I + 0 E1 s dx 1 s 0 0 E2 s dx 1 s + 0 E3 s dx 2 s I +, 0, (7.3 0 E3 s dx 3 s and observe in paricular ha E 1 = E (X 1 and E 3 = E (X 3 hold for 0. Insering he previous resuls in (3.3 yields for all 0 a.s. ( ( [( ( ( ] K (E V = = 1 0 K E 1 V (E 3 E 2 (E 1 (E 3 + s 0 Y 1 V 0 E 2 s E 3 d s s Y 2 s ( (E = 1 K (E 3 E 2 (E 1 K + (E 3 (7.4 V 0 ( (E E1 s dy 1 s (E 3 E 2 (E 1 0 E1 s dy 1 s + (E 3 [ 0 E3 s dy 2 s +. 0 E2 s dy 1 s] The equaion for he firs d enries of (7.4 is equivalen o K + E (X 1 s dy 1 s = E (X 1 K = K + E (X 1 s dx 1 sk a.s., 0, from where we deduce (4.3. From he equaion for he second m d enries of (7.4 we derive under use of (4.3, (2.4 and (7.3 ha V E (X 3 V 0 = E (X 3 = E (X 3 = E (X 3 = E (X 3 ( ( ( E 3 s dy 2 s + E 3 s dy 2 s + E 3 s dy 2 s + E 2 s dx 1 sk E 2 (E 1 (I + E 2 s dx 1 sk E 2 ( E 2 s dx 1 sk E (X 3 s d(y 2 s X 2 sk a.s., 0, E (X 1 E (X 1 K E 2 s dx 1 s + 0,] E 1 s dx 1 s K E 3 s dx 2 s K 24

25 such ha (4.4 is shown. Finally le (U, L 0 be he Lévy process defined in (2.9 and (3.5. By Theorem 3.4, (V 0 solves he SDE (3.4 wih respec o F. Observe ha by he same argumenaion as for X or alernaively by (2.9 we deduce ha for all 0 i holds U (i,j = 0 a.s. for i d, j > d. By insering U and L as given in (4.5 in he SDE (3.4 we obain L 1 = U 1 K in he firs and (4.6 in he second componen. This complees he proof. (b Insering (4.2 and (4.3 in (3.3 direcly gives he assumpion by calculaions similar as under (a. Proof of Corollary 4.5. By Remark 4.4, i only remains o show he equivalence of (4.7 and (4.8. For ha, observe ha E (U = E (X and hus E (OUO = E (OXO. Hence, as shown in he proof of Theorem 4.3, OXO has a d (m d block of zero enries in he upper righ corner, if and only if he same is rue for E (OXO, equivalenly for E (OXO, and hence equivalenly for OUO = Log ( E (OXO. I follows ha OX O is of he form as specified in (4.7 if and only if OU O is of he form specified in (4.8, and as seen in he proof of Theorem 4.3, i furher holds ha E (U 1 = E (X 1. To see he equivalence of he relaions regarding OY and OL, suppose firs ha Y saisfies (4.7. Then by (3.8, ( X OL = OY + [OUO 1, OY ] = K + Y 2 [( U 1 0 U 2 U 3, ( ] X 1 K Y 2, and he upper d componens on he righ hand side of his equaion are given by X 1 K + [U 1, X 1 K] = U 1 K, where he las equaion follows from (2.10 since E (U 1 = E (X 1. I follows ha L saisfies (4.8. Conversely, if (4.8 holds, hen i follows from (3.9 ha ( U OY = OL + [OXO 1, OL] = K + L 2 [( X 1 0 X 2 X 3, ( U 1 K L 2 ], and as above i follows from his equaion and (2.10 ha Y saisfies ( Proofs for Secion 5 In his secion we give he proofs for Secion 5 along wih a few resuls on mulivariae sochasic exponenials which will be needed bu are also ineresing in heir own righ. We sar by inroducing an approximaion of he sochasic exponenial which will be a useful ool. Namely, he following resul is due o Emery [8]. 25

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

6. Stochastic calculus with jump processes

6. Stochastic calculus with jump processes A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance. 1 An Inroducion o Backward Sochasic Differenial Equaions (BSDEs) PIMS Summer School 2016 in Mahemaical Finance June 25, 2016 Chrisoph Frei cfrei@ualbera.ca This inroducion is based on Touzi [14], Bouchard

More information

arxiv: v1 [math.pr] 19 Feb 2011

arxiv: v1 [math.pr] 19 Feb 2011 A NOTE ON FELLER SEMIGROUPS AND RESOLVENTS VADIM KOSTRYKIN, JÜRGEN POTTHOFF, AND ROBERT SCHRADER ABSTRACT. Various equivalen condiions for a semigroup or a resolven generaed by a Markov process o be of

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

4 Sequences of measurable functions

4 Sequences of measurable functions 4 Sequences of measurable funcions 1. Le (Ω, A, µ) be a measure space (complee, afer a possible applicaion of he compleion heorem). In his chaper we invesigae relaions beween various (nonequivalen) convergences

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

Generalized Snell envelope and BSDE With Two general Reflecting Barriers

Generalized Snell envelope and BSDE With Two general Reflecting Barriers 1/22 Generalized Snell envelope and BSDE Wih Two general Reflecing Barriers EL HASSAN ESSAKY Cadi ayyad Universiy Poly-disciplinary Faculy Safi Work in progress wih : M. Hassani and Y. Ouknine Iasi, July

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Optimality Conditions for Unconstrained Problems

Optimality Conditions for Unconstrained Problems 62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x

More information

5. Stochastic processes (1)

5. Stochastic processes (1) Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

More information

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,

More information

t 2 B F x,t n dsdt t u x,t dxdt

t 2 B F x,t n dsdt t u x,t dxdt Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.

More information

arxiv: v1 [math.fa] 9 Dec 2018

arxiv: v1 [math.fa] 9 Dec 2018 AN INVERSE FUNCTION THEOREM CONVERSE arxiv:1812.03561v1 [mah.fa] 9 Dec 2018 JIMMIE LAWSON Absrac. We esablish he following converse of he well-known inverse funcion heorem. Le g : U V and f : V U be inverse

More information

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS LECTURE : GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS We will work wih a coninuous ime reversible Markov chain X on a finie conneced sae space, wih generaor Lf(x = y q x,yf(y. (Recall ha q

More information

Quasi-sure Stochastic Analysis through Aggregation

Quasi-sure Stochastic Analysis through Aggregation E l e c r o n i c J o u r n a l o f P r o b a b i l i y Vol. 16 (211), Paper no. 67, pages 1844 1879. Journal URL hp://www.mah.washingon.edu/~ejpecp/ Quasi-sure Sochasic Analysis hrough Aggregaion H. Mee

More information

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

More information

Representation of Stochastic Process by Means of Stochastic Integrals

Representation of Stochastic Process by Means of Stochastic Integrals Inernaional Journal of Mahemaics Research. ISSN 0976-5840 Volume 5, Number 4 (2013), pp. 385-397 Inernaional Research Publicaion House hp://www.irphouse.com Represenaion of Sochasic Process by Means of

More information

Continuous time autoregressive moving average processes with random Lévy coefficients

Continuous time autoregressive moving average processes with random Lévy coefficients ALEA, La Am J Probab Mah Sa 4, 29 244 (207) Coninuous ime auoregressive moving average processes wih random Lévy coefficiens Dirk-Philip Brandes Ulm Universiy Insiue of Mahemaical Finance, Helmholzsr 8,

More information

Cash Flow Valuation Mode Lin Discrete Time

Cash Flow Valuation Mode Lin Discrete Time IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics

More information

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX J Korean Mah Soc 45 008, No, pp 479 49 THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX Gwang-yeon Lee and Seong-Hoon Cho Reprined from he Journal of he

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE Topics MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES 2-6 3. FUNCTION OF A RANDOM VARIABLE 3.2 PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE 3.3 EXPECTATION AND MOMENTS

More information

On Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems

On Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems MATHEMATICS OF OPERATIONS RESEARCH Vol. 38, No. 2, May 2013, pp. 209 227 ISSN 0364-765X (prin) ISSN 1526-5471 (online) hp://dx.doi.org/10.1287/moor.1120.0562 2013 INFORMS On Boundedness of Q-Learning Ieraes

More information

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL: Ann. Func. Anal. 2 2011, no. 2, 34 41 A nnals of F uncional A nalysis ISSN: 2008-8752 elecronic URL: www.emis.de/journals/afa/ CLASSIFICAION OF POSIIVE SOLUIONS OF NONLINEAR SYSEMS OF VOLERRA INEGRAL EQUAIONS

More information

Homework 10 (Stats 620, Winter 2017) Due Tuesday April 18, in class Questions are derived from problems in Stochastic Processes by S. Ross.

Homework 10 (Stats 620, Winter 2017) Due Tuesday April 18, in class Questions are derived from problems in Stochastic Processes by S. Ross. Homework (Sas 6, Winer 7 Due Tuesday April 8, in class Quesions are derived from problems in Sochasic Processes by S. Ross.. A sochasic process {X(, } is said o be saionary if X(,..., X( n has he same

More information

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,

More information

A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS

A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS Theory of Sochasic Processes Vol. 14 3), no. 2, 28, pp. 139 144 UDC 519.21 JOSEP LLUÍS SOLÉ AND FREDERIC UTZET A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS An explici procedure

More information

Monotonic Solutions of a Class of Quadratic Singular Integral Equations of Volterra type

Monotonic Solutions of a Class of Quadratic Singular Integral Equations of Volterra type In. J. Conemp. Mah. Sci., Vol. 2, 27, no. 2, 89-2 Monoonic Soluions of a Class of Quadraic Singular Inegral Equaions of Volerra ype Mahmoud M. El Borai Deparmen of Mahemaics, Faculy of Science, Alexandria

More information

On a Fractional Stochastic Landau-Ginzburg Equation

On a Fractional Stochastic Landau-Ginzburg Equation Applied Mahemaical Sciences, Vol. 4, 1, no. 7, 317-35 On a Fracional Sochasic Landau-Ginzburg Equaion Nguyen Tien Dung Deparmen of Mahemaics, FPT Universiy 15B Pham Hung Sree, Hanoi, Vienam dungn@fp.edu.vn

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

4.6 One Dimensional Kinematics and Integration

4.6 One Dimensional Kinematics and Integration 4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

Richard A. Davis Colorado State University Bojan Basrak Eurandom Thomas Mikosch University of Groningen

Richard A. Davis Colorado State University Bojan Basrak Eurandom Thomas Mikosch University of Groningen Mulivariae Regular Variaion wih Applicaion o Financial Time Series Models Richard A. Davis Colorado Sae Universiy Bojan Basrak Eurandom Thomas Mikosch Universiy of Groningen Ouline + Characerisics of some

More information

NEW EXAMPLES OF CONVOLUTIONS AND NON-COMMUTATIVE CENTRAL LIMIT THEOREMS

NEW EXAMPLES OF CONVOLUTIONS AND NON-COMMUTATIVE CENTRAL LIMIT THEOREMS QUANTUM PROBABILITY BANACH CENTER PUBLICATIONS, VOLUME 43 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 998 NEW EXAMPLES OF CONVOLUTIONS AND NON-COMMUTATIVE CENTRAL LIMIT THEOREMS MAREK

More information

CHAPTER 2 Signals And Spectra

CHAPTER 2 Signals And Spectra CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par

More information

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0.

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0. Advanced Financial Models Example shee 3 - Michaelmas 217 Michael Tehranchi Problem 1. Le f : [, R be a coninuous (non-random funcion and W a Brownian moion, and le σ 2 = f(s 2 ds and assume σ 2

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

On Carlsson type orthogonality and characterization of inner product spaces

On Carlsson type orthogonality and characterization of inner product spaces Filoma 26:4 (212), 859 87 DOI 1.2298/FIL124859K Published by Faculy of Sciences and Mahemaics, Universiy of Niš, Serbia Available a: hp://www.pmf.ni.ac.rs/filoma On Carlsson ype orhogonaliy and characerizaion

More information

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick

More information

Existence of positive solution for a third-order three-point BVP with sign-changing Green s function

Existence of positive solution for a third-order three-point BVP with sign-changing Green s function Elecronic Journal of Qualiaive Theory of Differenial Equaions 13, No. 3, 1-11; hp://www.mah.u-szeged.hu/ejqde/ Exisence of posiive soluion for a hird-order hree-poin BVP wih sign-changing Green s funcion

More information

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n Lecure 3 - Kövari-Sós-Turán Theorem Jacques Versraëe jacques@ucsd.edu We jus finished he Erdős-Sone Theorem, and ex(n, F ) ( /(χ(f ) )) ( n 2). So we have asympoics when χ(f ) 3 bu no when χ(f ) = 2 i.e.

More information

Sample Autocorrelations for Financial Time Series Models. Richard A. Davis Colorado State University Thomas Mikosch University of Copenhagen

Sample Autocorrelations for Financial Time Series Models. Richard A. Davis Colorado State University Thomas Mikosch University of Copenhagen Sample Auocorrelaions for Financial Time Series Models Richard A. Davis Colorado Sae Universiy Thomas Mikosch Universiy of Copenhagen Ouline Characerisics of some financial ime series IBM reurns NZ-USA

More information

Example on p. 157

Example on p. 157 Example 2.5.3. Le where BV [, 1] = Example 2.5.3. on p. 157 { g : [, 1] C g() =, g() = g( + ) [, 1), var (g) = sup g( j+1 ) g( j ) he supremum is aken over all he pariions of [, 1] (1) : = < 1 < < n =

More information

Short Introduction to Fractional Calculus

Short Introduction to Fractional Calculus . Shor Inroducion o Fracional Calculus Mauro Bologna Deparameno de Física, Faculad de Ciencias Universidad de Tarapacá, Arica, Chile email: mbologna@ua.cl Absrac In he pas few years fracional calculus

More information

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu ON EQUATIONS WITH SETS AS UNKNOWNS BY PAUL ERDŐS AND S. ULAM DEPARTMENT OF MATHEMATICS, UNIVERSITY OF COLORADO, BOULDER Communicaed May 27, 1968 We shall presen here a number of resuls in se heory concerning

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

GMM - Generalized Method of Moments

GMM - Generalized Method of Moments GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

The Strong Law of Large Numbers

The Strong Law of Large Numbers Lecure 9 The Srong Law of Large Numbers Reading: Grimme-Sirzaker 7.2; David Williams Probabiliy wih Maringales 7.2 Furher reading: Grimme-Sirzaker 7.1, 7.3-7.5 Wih he Convergence Theorem (Theorem 54) and

More information

Markov Processes and Stochastic Calculus

Markov Processes and Stochastic Calculus Markov Processes and Sochasic Calculus René Caldeney In his noes we revise he basic noions of Brownian moions, coninuous ime Markov processes and sochasic differenial equaions in he Iô sense. 1 Inroducion

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive. LECTURE 3 Linear/Nonnegaive Marix Models x ( = Px ( A= m m marix, x= m vecor Linear sysems of difference equaions arise in several difference conexs: Linear approximaions (linearizaion Perurbaion analysis

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Backward stochastic dynamics on a filtered probability space

Backward stochastic dynamics on a filtered probability space Backward sochasic dynamics on a filered probabiliy space Gechun Liang Oxford-Man Insiue, Universiy of Oxford based on join work wih Terry Lyons and Zhongmin Qian Page 1 of 15 gliang@oxford-man.ox.ac.uk

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

DISCRETE GRONWALL LEMMA AND APPLICATIONS

DISCRETE GRONWALL LEMMA AND APPLICATIONS DISCRETE GRONWALL LEMMA AND APPLICATIONS JOHN M. HOLTE MAA NORTH CENTRAL SECTION MEETING AT UND 24 OCTOBER 29 Gronwall s lemma saes an inequaliy ha is useful in he heory of differenial equaions. Here is

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a

More information

2 Some Property of Exponential Map of Matrix

2 Some Property of Exponential Map of Matrix Soluion Se for Exercise Session No8 Course: Mahemaical Aspecs of Symmeries in Physics, ICFP Maser Program for M 22nd, January 205, a Room 235A Lecure by Amir-Kian Kashani-Poor email: kashani@lpensfr Exercise

More information

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

More information

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

Generalized Chebyshev polynomials

Generalized Chebyshev polynomials Generalized Chebyshev polynomials Clemene Cesarano Faculy of Engineering, Inernaional Telemaic Universiy UNINETTUNO Corso Viorio Emanuele II, 39 86 Roma, Ialy email: c.cesarano@unineunouniversiy.ne ABSTRACT

More information

POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION

POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Novi Sad J. Mah. Vol. 32, No. 2, 2002, 95-108 95 POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Hajnalka Péics 1, János Karsai 2 Absrac. We consider he scalar nonauonomous neural delay differenial

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Undetermined coefficients for local fractional differential equations

Undetermined coefficients for local fractional differential equations Available online a www.isr-publicaions.com/jmcs J. Mah. Compuer Sci. 16 (2016), 140 146 Research Aricle Undeermined coefficiens for local fracional differenial equaions Roshdi Khalil a,, Mohammed Al Horani

More information

CHERNOFF DISTANCE AND AFFINITY FOR TRUNCATED DISTRIBUTIONS *

CHERNOFF DISTANCE AND AFFINITY FOR TRUNCATED DISTRIBUTIONS * haper 5 HERNOFF DISTANE AND AFFINITY FOR TRUNATED DISTRIBUTIONS * 5. Inroducion In he case of disribuions ha saisfy he regulariy condiions, he ramer- Rao inequaliy holds and he maximum likelihood esimaor

More information

SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990),

SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990), SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F Trench SIAM J Marix Anal Appl 11 (1990), 601-611 Absrac Le T n = ( i j ) n i,j=1 (n 3) be a real symmeric

More information

1 Solutions to selected problems

1 Solutions to selected problems 1 Soluions o seleced problems 1. Le A B R n. Show ha in A in B bu in general bd A bd B. Soluion. Le x in A. Then here is ɛ > 0 such ha B ɛ (x) A B. This shows x in B. If A = [0, 1] and B = [0, 2], hen

More information

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

More information

Math 10B: Mock Mid II. April 13, 2016

Math 10B: Mock Mid II. April 13, 2016 Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.

More information

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

More information

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013 Mahemaical Theory and Modeling ISSN -580 (Paper) ISSN 5-05 (Online) Vol, No., 0 www.iise.org The ffec of Inverse Transformaion on he Uni Mean and Consan Variance Assumpions of a Muliplicaive rror Model

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

THE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2.

THE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2. THE BERNOULLI NUMBERS The Bernoulli numbers are defined here by he exponenial generaing funcion ( e The firs one is easy o compue: (2 and (3 B 0 lim 0 e lim, 0 e ( d B lim 0 d e +e e lim 0 (e 2 lim 0 2(e

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

Heat kernel and Harnack inequality on Riemannian manifolds

Heat kernel and Harnack inequality on Riemannian manifolds Hea kernel and Harnack inequaliy on Riemannian manifolds Alexander Grigor yan UHK 11/02/2014 onens 1 Laplace operaor and hea kernel 1 2 Uniform Faber-Krahn inequaliy 3 3 Gaussian upper bounds 4 4 ean-value

More information

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8) I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

6.2 Transforms of Derivatives and Integrals.

6.2 Transforms of Derivatives and Integrals. SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.

More information

SOME MORE APPLICATIONS OF THE HAHN-BANACH THEOREM

SOME MORE APPLICATIONS OF THE HAHN-BANACH THEOREM SOME MORE APPLICATIONS OF THE HAHN-BANACH THEOREM FRANCISCO JAVIER GARCÍA-PACHECO, DANIELE PUGLISI, AND GUSTI VAN ZYL Absrac We give a new proof of he fac ha equivalen norms on subspaces can be exended

More information

RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY

RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY ECO 504 Spring 2006 Chris Sims RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY 1. INTRODUCTION Lagrange muliplier mehods are sandard fare in elemenary calculus courses, and hey play a cenral role in economic

More information

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange

More information