Verification Theorems for Stochastic Optimal Control Problems via a Time Dependent Fukushima - Dirichlet Decomposition

Size: px
Start display at page:

Download "Verification Theorems for Stochastic Optimal Control Problems via a Time Dependent Fukushima - Dirichlet Decomposition"

Transcription

1 Verificaion Theorems for Sochasic Opimal Conrol Problems via a Time Dependen Fukushima - Dirichle Decomposiion Fauso Gozzi, Francesco Russo To cie his version: Fauso Gozzi, Francesco Russo. Verificaion Theorems for Sochasic Opimal Conrol Problems via a Time Dependen Fukushima - Dirichle Decomposiion. Preprin LAGA-Paris To appear: Sochasic Processes and Their Applicaions. 34 pages <hal > HAL Id: hal hps://hal.archives-ouveres.fr/hal Submied on 13 Apr 2006 HAL is a muli-disciplinary open access archive for he deposi and disseminaion of scienific research documens, wheher hey are published or no. The documens may come from eaching and research insiuions in France or abroad, or from public or privae research ceners. L archive ouvere pluridisciplinaire HAL, es desinée au dépô e à la diffusion de documens scienifiques de niveau recherche, publiés ou non, émanan des éablissemens d enseignemen e de recherche français ou érangers, des laboraoires publics ou privés.

2 Verificaion Theorems for Sochasic Opimal Conrol Problems via a Time Dependen Fukushima - Dirichle Decomposiion Fauso GOZZI Diparimeno di Scienze Economiche e Aziendali Facolà di Economia LUISS - Guido Carli Viale Pola 12, I Roma, Ialy Francesco RUSSO Universié Paris 13 Insiu Galilée, Mahémaiques 99, av. JB Clémen, F Villeaneuse, France ccsd , version 1-13 Apr 2006 A.M.S. Subjec Classificaion: Primary: 35R15, 49L10, 60H05, 93B52, 93E20. Secondary: 35J60, 47D03, 47D07, 49L25, 60J35, 93B06, 93B50. Key words: Hamilon-Jacobi-Bellman (HJB) equaions, sochasic calculus via regularizaion, Fukushima-Dirichle decomposiion, sochasic opimal conrol, verificaion heorems. Absrac. This paper is devoed o presen a mehod of proving verificaion heorems for sochasic opimal conrol of finie dimensional diffusion processes wihou conrol in he diffusion erm. The value funcion is assumed o be coninuous in ime and once differeniable in he space variable (C 0,1 ) insead of once differeniable in ime and wice in space (C 1,2 ), like in he classical resuls. The resuls are obained using a ime dependen Fukushima - Dirichle decomposiion proved in a companion paper by he same auhors using sochasic calculus via regularizaion. Applicaions, examples and comparison wih oher similar resuls are also given. 1 Inroducion In his paper we wan o presen a mehod o ge verificaion heorems for sochasic opimal conrol problems of finie dimensional diffusion processes wihou conrol in he diffusion erm. The mehod is based on a generalized Fukushima - Dirichle decomposiion proved in he companion paper [22]. Since his Fukushima - Dirichle decomposiion holds for funcions u : [0, T] R n R ha are C 0 in ime and C 1 in space (C 0,1 in symbols), our verificaion heorem has he advanage of requiring less regulariy of he value funcion V han he classical ones which need C 1 regulariy in ime and C 2 in space of V (C 1,2 in symbols), see e.g. [13, pp. 140, 163, 172]. There are also oher verificaion heorems ha works in cases when he value funcion is nonsmooh: e.g. i is possible o prove a verificaion heorem in he case when V is only coninuous (see [23], [27], [39, Secion 5.2]) in he framework of viscosiy soluions. However all hese resuls applied o our cases are weaker han ours, for reasons ha are clarified in Secion 8. Since he mehod is a bi complex and ariculaed we presen firs in nex Secion 2 he saemen of our verificaion Theorems 2.7 and 2.8 in a model case wih simplified 1

3 assumpions (which subsanially yield nondegeneracy of he diffusion coefficien). In he same secion we also pu Subsecion 2.1 where we give some commens on he heorem, is applicabiliy and is relaionship wih oher similar resuls. Below we pass o he body of he paper giving firs some noaions in Secion 3 and hen presening he general saemens (including also possible degeneracy of he diffusion coefficien) and heir proof in Secion 4. Secion 5 is devoed o necessary condiions and opimal feedbacks. Secions 6 and 7 conain applicaions of our echnique o more specific classes of problems where oher echniques are more difficul o use. The firs is a case of exi ime problem where he HJB equaion is nondegenerae bu C 1,2 regulariy is no known o hold due o he lack of regulariy of he coefficiens; he second is a case where he HJB equaion is degenerae parabolic. Finally in Secion 8 we compare our resul wih oher verificaion echniques. 2 The saemen of he verificaion heorems in a model case To clarify our resuls we describe briefly and informally below he framework and he saemen of he verificaion heorem ha we are going o prove in a model case. The precise saemens and proofs are given in Secion 4; hen in Secions 6, 7 applicaions o more specific (and somehow more difficul) classes of problems are given. We decided his srucure since i is difficul o provide a single general resul: we can say ha we inroduce a echnique, based on he Fukushima-Dirichle decomposiion and on is represenaion given in [22], ha can be adaped wih some work o differen seings each ime wih a differen adapaion. ( ) Firs we ake a given sochasic basis Ω, (F s ) s 0, P ha saisfies he so-called usual condiions, a finie dimensional Hilber space A = R n (he sae space), a finie dimensional Hilber space E = R m (he noise space), a se U R k (he conrol space). We fix hen a erminal ime T [0, + ] (he horizon of he problem, ha can be finie or infinie bu is fixed), an iniial ime and sae (, x) [0, T] A (which will vary as usual in he dynamic programming approach). The sae equaion is (recall ha T = [, T] R) dy(s) = [F 0 (s, y (s)) + F 1 (s, y (s),z (s))] ds + B (s, y (s))dw(s), s T, y () = x, (1) where he following holds (for a marix B by B we mean i,j b ij and, given E, F finie dimensional spaces, by L(E, F) we mean he se of linear operaors from E o F). Hypohesis z : T Ω U (he conrol process) is measurable, locally inegrable in for a.e. ω Ω, adaped o he filraion (F s ) s 0 ; 2. F 0 : T 0 R n R n, F 1 : T 0 R n U R n and B : T 0 R n L(E, R n ) are coninuous. 3. There exiss C > 0 such ha T 0, x 1, x 2 R n, F 0 (, x 1 ) F 0 (, x 2 ),x 1 x 2 + B (, x 1 ) B (, x 2 ) 2 C x 1 x 2 2, F 0 (, x),x + B (, x) 2 C [1 + x 2] and here exiss a consan K such ha, T 0, x R n, z U, F 1 (, x 1, z) F 1 (, x 2, z) K x 1 x 2 F 1 (, x, z),x K (1 + x + z ). 2

4 4. For every x R n 0 [ F 0 (, x) + B (, x) 2] d < +. We call Z ad () he se of admissible conrol sraegies when he iniial ime and sae is (, x) defined as { } z : T Ω U, z measurable, adaped o (F Z ad () = s ) s 0, locally inegrable in s for a.e. ω Ω and call y (s;, x, z) he sae process associaed wih a given z Z ad (); his is he srong soluion of he equaion (1) which exiss and is unique for any given z hanks o Theorem 1.2 in [25, p.2]. Consider now he case T < +. We ry o minimize he funcional { } T J (, x; z) = E l (s, y (s;, x, z),z (s))ds + φ(y (T;, x, z)) (2) over all conrol processes z Z ad (). On he coefficiens we assume he following. Hypohesis 2.2 l : [0, T] R n U R and φ : R n R are coninuous and such ha for each z Z ad () he funcion (s, ω) l (s, y (s, ω),z (s, ω)) is lower semiinegrable in [, T] Ω (recall ha a real valued funcion f is lower semiinegrable in a measure space (M, µ) if M f dµ < +, where f sands for he negaive par of f i.e f = ( f) 0). Remark 2.3 The above requiremen implies ha J (, x; z) is well defined and > for each z Z ad (). I is assumed expliciely o cover a more general class of problems. This is obvious under addiional assumpions, e.g. if l and φ are bounded below. However some ineresing problems arising in economics conain funcions l like lnz and in his case he lower semiinegrabiliy of (s, ω) lnz (s, ω) for each z Z ad () follows by ad hoc argumens. These are problems wih sae consrains ha are no expliciely reaed in his paper bu we wan o se our framework so o be able o rea such class of problems. The value funcion is defined as V (, x) = inf J (, x; z) (3) z Z ad () and a conrol z Z ad () is opimal a (, x) if V (, x) = J (, x; z ). The curren value Hamilonian is defined, for (, x, p, z) [0, T] R n R n U as H CV (, x, p; z) = F 0 (, x),p + F 1 (, x, z),p + l (, x, z) and he (minimum value) Hamilonian as H (, x, p) = inf z U H CV (, x, p; z). Since he firs erm of H CV (, x, p; z) does no depend on he conrol z we will usually define H 0 CV (, x, p; z) = F 1 (, x, z),p + l (, x, z) and H 0 (, x, p) = inf z U H0 CV (, x, p; z), (4) 3

5 so we can wrie and H CV (, x, p; z) =: F 0 (, x),p + H 0 CV (, x, p; z) (5) H (, x, p) =: F 0 (, x),p + H 0 (, x, p). (6) When T < + he Hamilon-Jacobi-Bellman (HJB) equaion for he value funcion is a semilinear parabolic PDE ((, x) [0, T] R n ) v (, x) = 1 2 Tr[B (, x) xx v(, x)b (, x)] + F 0 (, x), x v (, x) + H 0 (, x, x v (, x)), wih he final condiion v (T, x) = φ(x), x R n. (8) To ensure finieness and coninuiy of he Hamilonian we need also o add he following assumpion. Hypohesis 2.4 The funcions l, F 1 and he se U are such ha he value funcion is always finie and he Hamilonian H 0 (, x, p) is well defined, finie and coninuous for every (, x, p) [0, T] R n R n. Remark 2.5 The above Hypohesis 2.4 is saisfied e.g. when l, φ are coninuous and U is compac. Anoher possibiliy is o ake U unbounded, F 1 sublinear, l (, x, z) = g (x) + h (z) wih g coninuous and bounded, h coninuous and such ha (7) h(z) / z + as z + ; his case was sudied e.g in [6, 18, 21] even in he infinie dimensional case (see on his he Secions 6, 7). The saemen of he classical verificaion heorem for his model problem when T < + is he following, see Definiion 4.5 for he definiion of sric soluion of equaion (7)-(8). Theorem 2.6 Assume ha Hypoheses 2.1, 2.2, 2.4 hold rue. Le v C 1,2 ([0, T] R n ) be a polynomially growing sric soluion of he HJB equaion (7)-(8) on [0, T] R n. Then he wo following properies hold rue. (i) v V on [0, T] R n ; (ii) fix (, x) [0, T] R n ; if z Z ad () is such ha, calling y (s) = y (s;, x, z), H 0 (s, y (s), x v (s, y (s))) = H 0 CV (s, y (s), x v (s, y (s)),z (s)), P a.s., for a.e. s [, T], hen z is opimal a (, x) and v (, x) = V (, x). This heorem saes a sufficien opimaliy condiion and is proof is based on Iô s formula, see e.g. [39, p.268]. In he classical conex also necessary condiions and exisence of opimal feedback can be proved, see again [39, p.268]. The saemen of our Verificaion Theorem in he model problem described above is very similar. We give i in he case when he following nondegeneracy hypohesis hold B 1 (, x)f 1 (, x, z) is bounded on [0, T] R n U, (9) leaving he discussion for he general case o Secion 4 (srong soluions of he HJB equaion (7)-(8) are defined in Definiion 4.6). 4

6 Theorem 2.7 Assume ha Hypoheses 2.1, 2.2, 2.4 and (9) hold rue. Le v C 0,1 ([0, T] R n ) be a polynomially growing srong soluion of he HJB equaion (7)-(8) on [0, T] R n. Then (i) and (ii) of Theorem 2.6 above hold. In fac we will go furher: proving a resul in he case when (9) does no hold; showing ha in our seing we can obain a necessary condiion and he exisence of opimal feedback conrols on he same line of he classical resuls, see Secion 5. Moreover in a paper in preparaion we will show ha also some cases where he drif is a disribuion can be reaed and also cases where he soluion of HJB enjoys weaker regulariy. When T = +, F 0, F 1 and B do no depend on ime, we consider he problem of minimizing { + } J (, x; z) = E e λ l 1 (y (s;, x, z),z (s))ds, where l (, x, z) = e λ l 1 (x, z), wih λ > 0, l 1 coninuous and bounded. In his case he value funcion V is defined as in (3) bu is dependence on becomes rivial as V (, x) = e λ V (0, x) for every and x. Then we se = 0 and call V 0 (x) = V (0, x). The HJB equaion for V 0 (x) becomes an ellipic PDE (x R n ) λv (x) = 1 2 Tr[B (x) xx v(x)b (x)] + F 0 (, x), x v (, x) + H 1 (x, x v (, x)), (10) where H 1 (x, p) = inf { F 1 (x, z),p + l 1 (x, z)} =: inf z U z U H1 CV (x, p; z). (11) In his ellipic case he saemen of our verificaion heorem is he following (for he definiion of srong soluion of he HJB equaion (10) see Definiion 4.14). Theorem 2.8 Assume ha Hypoheses 2.1 and (9) hold. Assume also ha l 1 is coninuous and bounded and ha Hypohesis 2.4 hold rue wih H 1 in place of H 0. If v is a bounded srong soluion of he HJB equaion (10) and v C 1 (R n ) hen v V 0. Moreover fix x R n ; if z Z ad (0) is such ha, calling y (s) = y (s; 0, x, z), H 1 (y (s), x v (s, y (s))) = H 1 CV (y (s), xv (s, y (s)), z (s)), for a.e. s 0, P a.s., hen z is opimal a (0, x) and v (x) = V 0 (x). 2.1 Oher Verificaion Theorems for sochasic opimal conrol problems We recall ha in he lieraure on sochasic opimal conrol oher verificaion heorems for no C 1,2 funcions were proved. Roughly speaking we can say ha each of hem (and he echnique of proof, oo) is sricly conneced wih he concep of weak soluion of he HJB equaion ha is considered. We have in fac he following. (Srong soluions). In [18, 21] (see also [7, 11, 17]) in a seing similar o ours bu in infinie dimension, a verificaion heorem is proved assuming ha here exiss a srong soluion v of HJB which is C 1 in space. In fac he mehod oulined here goes in he 5

7 same direcion, bu improves and generalizes in he finie dimensional case he ideas conained here. Such improvemen is no a sraighforward one, as we need o use compleely differen ools o ge our resuls and we ge more powerful heorems. We refer o Secions 6 and 7 for examples and o Secion 8 for explanaions. (Viscosiy soluions). In [27, 39, 23] a verificaion echnique for viscosiy soluions is inroduced and sudied. Such echnique adaps o he case when v is only coninuous and so is very general and applicable o cases when he conrol eners in he diffusion coefficiens B. However in he case of our ineres, i.e. when v is C 1 in space, such echnique gives weaker resuls as i requires more assumpions on he coefficiens F 0, F 1, B and on he candidae opimal sraegy; see Secion 8 for explanaions. (Mild soluions). In [14, Theorem 7.2] a verificaion heorem is given when he HJB equaion admis mild soluions which are C 1 in space (see Definiion 6.13 or Remark 7.2) and when he soluion can be represened using he soluion of a suiable backward SDE. The resuls available wih his echnique are limied o infinie dimensional cases where mild soluions exis and Girsanov heorem can be applied in a suiable way. Such resricions preven he use of his echnique (a leas o he presen sae of he ar) e.g. in he cases described in Secions 6 and 7, see Secion 8 for explanaions. Moreover we wan o sress he fac ha he Fukushima-Dirichle decomposiion we go in [22] is in fac sronger han wha we need o prove he verificaion heorems above. In fac for such purpose i would be enough o prove a Dynkin-ype formula (roughly speaking an Iô formula afer expecaion) which is in general easier. We decided o sae and prove such Fukushima-Dirichle decomposiion since i can help o deal wih sochasic conrol problems where he crierion o minimize does no conain expecaion (e.g. pahwise opimaliy and opimaliy in probabiliy). In such conex he HJB equaion becomes a sochasic PDE and o ge a verificaion heorem, Dynkin-ype formulae canno be used. See e.g. [10, 28, 29, 33] on his subjec. To sum up we hink ha he ineres of our verificaion resuls is he fac ha we obain hem using only he fac ha he soluion v of HJB belongs o C 0,1 and no necessarily o C 1,2 obaining beer resuls han oher echniques. hey can be applied also o problems wih pahwise opimaliy and opimaliy in probabiliy. Remark 2.9 As a firs sep of our work, we consider here he case of a sochasic opimal conrol problem in finie dimension wih no sae consrains. We are aware of he fac ha in many cases, he HJB equaion associaed wih he conrol problem admis a C 1,2 soluion, so ha our mehod does no give a real advanage in his case. However in some cases, mainly degenerae ones (see Secions 6 and 7) i is known ha he soluion of HJB is C 0 (or C α wih Hölder exponen α > 0) in ime and C 1 in space bu C 1,2 regulariy is no known a his sage. In paricular we show an explici example where he value funcion is C 1 bu canno be C 2 in space; we hink ha such echnique could be exended also in cases where HJB equaion is inended in a generalized sense; in a paper in preparaion, we invesigae a verificaion heorem relaed o an SDE whose drif is he derivaive of a coninuous funcion herefore a Schwarz disribuion, see for insance [12]. Tha sae equaion models a sochasic paricle moving in a random (irregular) medium; 6

8 we consider problems wihou sae consrains; again we hink ha some sae consrains problems can be reaed wih our approach bu we do no handle hem here: we provide an example wih exi ime which usually presens similar difficulies; we hink ha our mehod can be exended o he infinie dimensional case, where C 1,2 regulariy resuls for he soluion of HJB equaion are much less known while C 1 regulariy can be found in a variey of cases (see e.g. [6, 18, 21, 7]). We do no perform his exension here leaving i for furher work. Remark 2.10 Similar ideas as in his work have been expressed wih a differen formalism in [8]. There he auhors have implemened a generalized Iô formula in he Krylov spiri for proving an exisence and uniqueness resul for a generalized soluion of Bellman equaion for conrolled processes, in a non-degeneracy siuaion. 3 Noaions Throughou his paper we will denoe by (Ω, F, P) a given sochasic basis, where F sands for a given filraion (F s ) s 0 saisfying he usual condiions. Given a finie dimensional real Hilber space E, W will denoe a cylindrical Brownian moion wih values in E and adaped o (F s ) s 0. Given 0 T + and seing T = [, T] R he symbol C F (T Ω; E), will denoe he space of all coninuous processes adaped o he filraion F wih values in E. This is a Fréche space if endowed wih he opology of he uniform convergence in probabiliy (u.c.p. from now on). To be more precise his means ha, given a sequence (X n ) C F (T Ω; E) and X C F (T Ω; E) we have if and only if for every ε > 0, 1 T lim sup n + s [, 1] X n X P ( X n s X s E > ε) = 0. Given a random ime τ and a process (X s ) s T, we denoe by X τ he sopped process defined by Xs τ = X s τ. The space of all processes in [, T], adaped o F and square inegrable wih values in E is denoed by L 2 F (, T; E). Sn will denoe he space of all symmeric marices of dimension n. Le k N. As usual C k (R n ) is he space of all funcions : R n R ha are coninuous ogeher wih heir derivaives up o he order k. This is a Fréche space equipped wih he seminorms sup u (x) R + sup x u (x) R n + sup xx u (x) R n n +... (12) x K x K x K for every compac se K R n. This space will be denoed simply by C k when no confusion may arise. If K ia compac subse of R n hen C k (K) is a Banach space wih he norm (12). The symbol C k b (Rn ) will denoe he Banach space of all funcions from R n o R ha are coninuous and bounded ogeher wih heir derivaives up o he order k. This space is endowed wih he usual sup norm. Passing o parabolic spaces we denoe by C 0 (T R n ) he space of all funcions u : T R n R, (s, x) u (s, x) 7

9 ha are coninuous. This space is a Fréche space equipped wih he seminorms sup u (s, x) R (s,x) [, 1] K for every 1 > 0 and every compac se K R n ). Moreover we will denoe by C 1,2 (T R n ) (respecively C 0,1 (T R n )), he space of all funcions u : T R n R, (s, x) u (s, x) ha are coninuous ogeher wih heir derivaives u, x u, xx u (respecively x u). This space is a Fréche space equipped wih he seminorms (respecively sup u (s, x) R + sup s u (s, x) R n (s,x) [, 1] K (s,x) [, 1] K + sup x u (s, x) R n + sup xx u (s, x) R n n (s,x) [, 1] K (s,x) [, 1] K sup u (s, x) R + sup x u (s, x) R n ) (s,x) [, 1] K (s,x) [, 1] K for every 1 > 0 and every compac se K R n. This space will be denoed simply by C 1,2 (respecively C 0,1 ) when no confusion may arise. Similarly, for α, β [0, 1] one defines C α,1+β (T R n ) (or simply C α,1+β ) as he subspace of C 0,1 (T R n ) of funcions u : T R n R such ha are u (, x) is α Hölder coninuous and x u (s, ) is β Hölder coninuous (wih he agreemen ha 0-Hölder coninuiy means jus coninuiy). If such properies hold jus locally hen such space is denoed by C α,1+β loc (T R n ). Similarly, given an open subse of O of R n one can define C α,1+β (T O) and C α,1+β loc (T O). If K is a compac subse of T R n we define C k (K) or C α,1+β (K) as above. Similarly o Cb k (Rn ) we define he Banach spaces Cb 0 (T R n ) C 1,2 b (T R n ), C α,1+β b (T R n ), C 0,1 b (T R n ). 4 Proof of he verificaion heorems in he model case In his secion we give he precise saemen and he proof of he verificaion Theorems 2.7 and 2.8 for he model problem described in Secion 2. Then in Secions 6 and 7 we will consider wo families of problems o which our echnique apply. Consider firs he parabolic case fixing he horizon T < +. Under Hypoheses 2.1, 2.2 and 2.4 we seek o minimize he cos J (, x; z) given in (2) over all admissible z Z ad () where he sae equaion is given by (1). We define he operaor L 0 : D (L 0 ) C 0 ([0, T] R n ) C 0 ([0, T] R n ), D (L 0 ) = C 1,2 ([0, T] R n ), L 0 v (, x) = v (, x) + F 0 (, x), x v (, x) Tr [B (, x) xx v (, x)b (, x)]. The HJB equaion associaed wih he problem (1) - (2) can hen be wrien as where H 0 is given in (4). L 0 v (, x) + H 0 (, x, x v (, x)) = 0, v (T, x) = φ(x), (13) 8

10 Now we wan o apply he represenaion resul proved in [22] Secion 4, which we recall below for he reader s convenience. Firs we consider he following Cauchy problem for h C 0 ([0, T] R n ). wih he following definiions of soluion. L 0 u (s, x) = h (s, x), u (T, x) = φ(x), (14) Definiion 4.1 We say ha u C 0 ([0, T] R n ) is a sric soluion o he backward Cauchy problem (14) if u D (L 0 ) and (14) holds. Definiion 4.2 We say ha u C 0 ([0, T] R n ) is a srong soluion o he backward Cauchy problem (14) if here exiss a sequence (u n ) D (L 0 ) and wo sequences (φ n ) C 0 (R n ), (h n ) C 0 ([0, T] R n ), such ha 1. For every n N u n is a sric soluion of he problem L 0 u n (, x) = h n (, x), u n (T, x) = φ n (x). 2. The following limis hold u n u in C 0 ([0, T] R n ), h n h in C 0 ([0, T] R n ), φ n φ in C 0 (R n ). The represenaion resul is he following (Corollaries 4.6 and 4.8 of [22]). Theorem 4.3 Le b 1 : [0, T] R n Ω R n, be a coninuous progressively measurable field (coninuous in (s, x)) and b : [0, T] R n R n, σ : [0, T] R n L(R m, R n ), be coninuous funcions. Le u C 0,1 ([0, T] R n ) be a srong soluion of he Cauchy problem (14). Fix [0, T], x R n and le (S s ) be a soluion o he SDE Then u (s, S s ) = u (, S ) + + ds s = b 1 (s, S s )ds + σ (s, S s )dw s, S = x. s 0 s h (r, S r )dr + s x u (r, S r )σ (r, S r )dw r. x u (r, S r ), b 1 (r, S r ) b (r, S r ) dr provided ha: eiher we can choose he approximaing sequence (u n ) of Definiion 4.2 so ha for every 0 s T lim n + s or he funcion x u n (r, S r ) x u (r, S r ),b 1 (r, S r ) b (r, S r ) dr = 0, u.c.p., (15) (, x, ω) σ 1 (, x)[b 1 (r, x, ω) b (r, x)], (where σ 1 sands for he pseudo-inverse of σ), is well defined and bounded on [0, T] R n Ω. 9

11 Remark 4.4 If lim xu n = x u, in C 0 ([0, T] R n ), n + hen Assumpion (15) is verified. This means ha he resul of Theorem 4.3 above applies if we know ha u is a srong soluion in a more resricive sense, i.e. subsiuing he poin 2 of Definiion 4.2 wih u n u in C 0 ([0, T] R n ), x u n x u in C 0 ([0, T] R n ), h n h in C 0 ([0, T] R n ), φ n φ in C 0 (R n ). This is a paricular case of our seing and i is he one used e.g in [18, 21] o ge he verificaion resul. We can say ha in hese papers a resul like Theorem 4.3 is proved under he assumpion ha u is a srong soluion in his more resricive sense. I is worh o noe ha in such simplified seing he proof of Theorem 4.3 follows simply by using sandard convergence argumens. In paricular here one does no need o use he Fukushima-Dirichle decomposiion presened in Secion 3 of [22]. So, from he mehodological poin of view here is a serious difference wih he resul of Theorem 4.3, see Secion 8 for commens. To apply Theorem 4.3 o our case we need firs o adap he noion of srong soluion and o rewrie he main assumpions. Le us give he following definiions. Definiion 4.5 We say ha v C 0 ([0, T] R n ) is a sric soluion of he HJB equaion (13) if v D (L 0 ) and (13) holds on [0, T] R n. Definiion 4.6 A funcion v C 0,1 ([0, T] R n ) is a srong soluion of he HJB equaion (13) if, seing h (, x) = H 0 (, x, x v (, x)), v is a srong soluion of he backward Cauchy problem L 0 v (, x) = h (, x), v (T, x) = φ(x), in he sense of Definiion 4.2. We now presen our verificaion heorem in he exended version, as announced in Secion 2. We need he following assumpion. Hypohesis 4.7 There exiss a funcion v C 0,1 ([0, T] R n ) which is a srong soluion of he HJB equaion (13) in he sense of Definiion 4.6 and is polynomially growing wih is space derivaive in he variable x. Moreover: (i) eiher we can choose he approximaing sequence (v n ) of Definiion 4.6 so ha for every 0 s T and for every admissible conrol z Z ad () s lim n + (ii) or he funcion x v n (r, y (r)) x v (r, y (r)),f 1 (r, y (r), z (r)) dr = 0, u.c.p., (, x, z) B 1 (, x) F 1 (, x, z), where B 1 sands for he pseudo-inverse of B, is well defined and bounded on [0, T] R n U. Remark 4.8 All he resuls below sill hold rue wih suiable modificaions if we assume ha: 10

12 he srong soluion v belongs o C 0 ([0, T] R n ) C 0,1 ([ε, T] R n ) for every small ε > 0; for some β (0, 1) he map (, x) β x v (, x) belongs o C 0 ([0, T] R n ). The proof of Theorem 4.9 in his case is a sraighforward generalizaion of he one presened here: we do no give i here o avoid echnicaliies since here we deal wih a model problem. In Secion 6 we will rea a case wih such difficuly. See also Remark 4.10 of [22] on his. We give now here he precise saemen of Theorem 2.7. Theorem 4.9 Assume ha Hypoheses 2.1, 2.2, 2.4 and 4.7 hold. Le H 0 CV, H0 be as in (4)-(5). Le v C 0,1 ([0, T] R n ) be a srong soluion of (13) and fix (, x) [0, T] R n which is polynomially growing wih is space derivaive in he variable x. Then (i) v V on [0, T] R n. (ii) If z is an admissible conrol a (, x) ha saisfies (seing y (s) = y (s;, x, z)) H 0 (s, y (s), x v (s, y (s))) = H 0 CV (s, y (s), x v (s, y (s));z (s)), (16) for a.e. s [, T], P almos surely, hen z is opimal a (, x) and v (, x) = V (, x). The proof of his heorem follows by he following fundamenal ideniy ha we sae as a lemma. Lemma 4.10 Assume ha Hypoheses 2.1, 2.2, 2.4 and 4.7 hold. Le v C 0,1 ([0, T] R n ) be a srong soluion of (13). Then, for every (, x) [0, T] R n and z Z ad () such ha J (, x; z) < + he following ideniy holds J (, x; z) = v (, x) [ +E H 0 (s, y(s), x v (s, y (s))) + HCV 0 (s, y(s), xv (s, y (s));z(s)) ] ds (17) where y(s) def = y (s;, x, z) is he soluion of (1) associaed wih he conrol z. Proof. For he sake of compleeness we firs show how he proof goes in he case when v C 1,2 ([0, T] R n ) and is a sric soluion of (13). We use Iô s formula applied o he funcion v (, x) by obaining, for every (, x) [0, T] R n, v (T, y (T)) v (, y ()) = x v (s, y (s)), B (s, y (s))dw (s) F 0 (s, y (s)), x v (s, y (s)) ds (18) F 1 (s, y (s),z (s)), x v (s, y (s)) ds s v (s, y (s))ds Tr [B (s, y (s)) xx v (s, y (s))b (s, y (s))] ds, 11

13 (which is exacly he decomposiion of Proposiion 2.4 of [22]) so, by aking expecaion of boh sides (his is finie since we have he polynomial growh assumpion) we ge he so-called Dynkin formula Ev (T, y (T)) = Ev (, y ())+E Now by (13) we ge, for every s 0, L 0 v (s, y (s))ds+e L 0 v (s, y (s)) = H 0 (s, y (s), x v (s, y (s))), F 1 (s, y (s),z (s)), x v (s, y (s)) ds. which yields Ev (T, y (T)) = Ev (, y ()) [ +E H 0 (s, y (s), x v (s, y (s))) + F 1 (s, y (s),z (s)), x v (s, y (s)) ] ds. Now since we assumed ha J (, x; z) < + and since Hypohesis 2.2 implies J (, x; z) > (his also follows observing ha he las line is no smaller han J (, x; z) and ha he firs line is finie by he polynomial growh of v), we add J (, x; z) o boh erms, use ha v (T, y (T)) = φ(y (T)) and ha y () = x o ge which is he claim (17). J (, x; z) = v (, x) + E H 0 (s, y (s), x v (s, y (s)))ds +E [ F 1 (s, y (s),z (s)), x v (s, y (s)) + l (s, y (s),z (s))] ds We now show how he proof goes when v C 0,1 ([0, T] R n ). The difference is due o he fac ha he erm I := s v (s, y (s))ds Tr [B (s, y (s)) xx v (s, y (s))b (s, y (s))] ds in he hird and fourh line of equaion (18) is no well defined now. However by Hypohesis 4.7 we know ha v is a srong soluion (in he sense of Definiion 4.2) of L 0 v (, x) = h 0 (, x), v (T, x) = φ(x), where we have se h 0 (, x) = H 0 (, x, x v (, x)). Then applying now Theorem 4.3 for he operaor L 0 (seing b = F 0, σ = B and b 1 = F 0 + F 1 ) we ge ha v (T, y (T)) = v (, y ()) + + x v (s, y (s)),b (s, y (s))dw (s) [ H 0 (s, y (s), x v (s, y (s))) + F 1 (s, y (s), z (s)), x v (s, y (s)) ] ds. Now, aking expecaion, adding and subracing J (, x; z) (which is a.s. finie using he same argumen for he smooh case), using ha v (T, y (T)) = φ(y (T)) and ha y () = x we ge [ ] T E l (s, y (s),z (s))ds + φ(y (T)) = v (, x) 12

14 [ +E H 0 (s, y (s), x v (s, y (s))) + HCV 0 (s, y (s), xv (s, y (s));z (s)) ] ds and he claim again follows. Proof of Theorem 4.9. By he definiion of H 0 and H 0 CV, for every (, x) [0, T] Rn, z Z ad () and for a.e. s T, he following inequaliy holds P-almos surely (seing y (s) = y (s;, x, z)) H 0 (s, y (s), x v (s, y (s))) + H 0 CV (s, y (s), xv (s, y (s)) ; z (s)) 0, (19) and hen by he fundamenal ideniy (17) i follows v (, x) J (, x; z) for every admissible z. This gives v V and so par (i) of he Theorem. Now consider an admissible conrol z such ha, for a.e. s [, T], P a.s. H 0 (s, y (s), x v (s, y (s))) = H 0 CV (s, y (s), xv (s, y (s));z (s)). (20) By he fundamenal ideniy (17) we have v (, x) = J (, x; z), which implies opimaliy and v (, x) = V (, x). Remark 4.11 To ge he mehodological novely in he above proof i is crucial o noe he following. Since we know ha u is he limi of classical soluions, he firs idea o approach he above proof (and ha has been used e.g. in he papers [18, 21]), is probably o prove he fundamenal ideniy for he approximaing soluions u n and hen pass o he limi for n. However if one ries o do his one needs o use he uniform convergence of he space derivaives x u n o x u. This is no needed wih our mehod, see also on his Remark 4.4. Le us go now o he proof of Verificaion Theorem 2.8 for he infinie horizon case. Se T = +, assume ha F 0, F 1 and B do no depend on ime and ha l (, x, z) = e λ l 1 (x, z) (λ > 0), φ = 0; hen we can se = 0 (since he dependence of he value funcion V on becomes rivial in his case) and he HJB equaion for V 0 becomes an ellipic PDE. We define he operaor L 0 : D (L 0 ) C 0 (R n ) C 0 (R n ), D (L 0 ) = C 2 (R n ), L 0 u (x) = F 0 (x), x u (x) Tr [B (x) xx u (x) B (x)] and we rewrie he HJB equaion (10) as λu (x) + L 0 u (x) + H 1 (x, x u (x)) = 0, x R n. (21) where H 1 is defined in (11). We now define srong soluions for such equaions. Consider firs, for h C 0 (R n ) he inhomogenous ellipic problem λu (x) + L 0 u (x) + h (x) = 0, x R n. (22) Definiion 4.12 We say ha u is a sric soluion o he ellipic problem (22) if u D (L 0 ) and (22) holds. 13

15 Definiion 4.13 We say ha u is a srong soluion o he ellipic problem (22) if here exiss a sequence (u n ) D (L 0 ) and a sequence (h n ) C 0 (R n ), such ha 1. For every n N u n is a sric soluion of he problem λu n (x) L 0 u n (x) = h n (x), x R n. 2. The following limis hold u n u in C 0 (R n ), h n h in C 0 (R n ). For funcions u C 1 (R n ) ha are srong soluions of he Cauchy problem (22) he resul of Theorem 4.3 sill holds: i is indeed a simpler case. Now we go o srong soluions of he HJB equaion (21). Definiion 4.14 A funcion v C 1 (R n ) is a srong soluion of he HJB equaion (21) if, seing h 0 (x) = H 1 (x, x v (x)), v is a srong soluion of he linear problem in he sense of Definiion Assume he following. λv (x) + L 0 v (x) = h 0 (x), Hypohesis 4.15 There exiss a funcion v C 1 (R n ) which is a srong soluion of he HJB equaion (21) in he sense of Definiion Moreover: (i) eiher we can choose he approximaing sequence (v n ) of Definiion 4.14 so ha for every s 0 and for every admissible conrol z Z ad (0) s lim n + 0 (ii) or he funcion x v n (y (r)) x v (y (r)),f 1 (y (r), z (r)) dr = 0, u.c.p., (x, z) B 1 (x)f 1 (x, z), where B 1 sands for he pseudo-inverse of B, is well defined and bounded on R n U. The precise saemen of he verificaion heorem is he following. Theorem 4.16 Assume ha Hypoheses 2.1, and 4.15 hold. Assume also ha l 1 is coninuous and bounded and ha Hypohesis 2.4 hold rue wih H 1 in place of H 0. If v is a bounded srong soluion of he HJB equaion (21) and v C 1 (R n ) hen v V 0 on R n. Morerover if z is an admissible conrol a (0, x) ha saisfies (seing y (s) = y (s; 0, x, z)) H 1 (y (s), x v (y (s))) = H 1 CV (y (s), x v (y (s));z (s)) for a.e. s [0, + ), P almos surely, hen z is opimal and v (x) = V 0 (x). Proof. The proof goes along he same lines as he finie horizon case. Firs we observe ha he boundedness of he daum l 1 implies ha also V 0 is bounded and ha for every x R n, z Z ad (0) he funcional J (0, x; z) is finie. Then we prove he following fundamenal ideniy J (0, x; z) = v (x) + E + 0 e λs [ H 1 (y(s), x v (y (s))) + H 1 CV (y(s), xv (y (s));z(s)) ] ds. 14

16 To do i we apply Theorem 4.3 for he operaor L 0 = + L 0 (so b = F 0, σ = B and b 1 = F 0 + F 1 ) bu aking = 0 and replacing v (s, x) by e λs v (x) for s 0. Observe firs ha, if v C 2 (R n ) hen he funcion w (, x) = e λ v (x) solves, for 0, x R n he equaion w (, x) + L 0 w (, x) = e λ H 1 (x, x v (x)). (23) Moreover if v C 1 (R n ) is a srong soluion of (21), also w (, x) is a srong soluion of (23) and saisfies he assumpions of Theorem 4.3. So we ge for every T 1 > 0 e λt1 v (y (T 1 )) = v (y (0)) e λs x v (y (s)),b (y (s))dw (s) 1 + e λs [ H 1 (s, y (s), x v (s, y (s))) + F 1 (s, y (s),z (s)), x v (s, y (s)) ] ds. 0 Now, aking expecaion, adding and subracing E 1 0 e λs l 1 (y (s),z (s))ds (which is a.s. finie by he boundedness of l 1 ) and using ha y (0) = x we ge [ ] T1 E e λs l 1 (y (s),z (s))ds + e λt1 v (y (T 1 )) 0 = v (x) + E 1 0 e λs [ H 1 (y (s), x v (y (s))) + H 1 CV (y (s), x v (y (s));z (s)) ] ds and he claim again follows aking he limi for T 1 + and using he boundedness of v. The res of he proof is exacly he same as in he finie horizon case. Remark 4.17 I mus be noed ha he assumpion on he boundedness of l 1 and v is made for simpliciy of exposiion. Wha is really needed in his infinie horizon case is ha J is well defined for each z Z ad (0) and ha lim T1 + e λt1 v (y (T 1 )) = 0 for every z Z ad (0). This allows o pass o he limi in he laer formula. In some cases his can be checked direcly, in oher cases much weaker condiions can be imposed (e.g. sublineariy of l 1 and esimaes on he soluion y as T 1 + ). 5 Necessary Condiions and Opimal Feedback Conrols Here we wan simply o show ha under addiional assumpions (mainly abou exisence and/or uniqueness of he maximum of he Hamilonian and of he soluion of he closed loop equaion) we can ge necessary condiions and opimal feedback conrols. This is a consequence of he verificaion Theorems 4.9 and 4.16 which has some imporance for applicaions. For he sake of breviy we give he resuls only for he finie horizon case observing ha compleely analogous resuls hold rue for he infinie horizon case. We sar by making a remark abou he so-called weak formulaion of a sochasic conrol problem. Remark 5.1 The seing inroduced in Secion 2 corresponds o he so-called srong formulaion of a sochasic conrol problem. For cerain purposes (namely he exisence of opimal feedbacks) i is convenien o consider he weak formulaion, leing he sochasic basis vary. In such formulaion one considers as he se Z ad () of admissible conrols as he se of 5-uples (Ω, F, P, W, z) such ha (Ω, F, P) is a complee filered probabiliy space wih he filraion F saisfying he usual condiions, 15

17 W is an E valued m-dimensional Brownian moion on [, T], he process z is measurable, (F s )-adaped and a.s. locally inegrable. We will use he noaion (Ω, F, P, W, z) Z ad () and we will call his conrol sraegies weakly admissible (weakly opimal when hey are opimal). When no ambiguiy arises we will leave aside he probabiliy space and he whie noise (regarding i as fixed) and simply consider z Z ad (). We will specify when he weak concep of admissible conrol is needed. The same will be done in Secions 6 and 8. One can see e.g. [39, p.64], [13, pp.141, 160] for commens on hese formulaions. We sar recalling a case when he sufficien condiion of he verificaion Theorem 4.9 becomes also necessary. Proposiion 5.2 Assume ha he hypoheses of Theorem 4.9 hold rue. We also assume ha he value funcion V is a srong soluion of he HJB equaion. Then an admissible conrol z Z ad () is opimal a (, x) [0, T] R n if and only if (16) holds wih V in place of v. Proof. The fundamenal ideniy in his case gives J (, x; z) = V (, x) [ +E H 0 (s, y (s), x V (s, y (s))) HCV 0 (s, y (s), x V (s, y (s));z (s)) ] ds and his immediaely gives he claim. Remark 5.3 One case where he above Proposiion 5.2 can be applied is when: i is known ha V is he unique viscosiy soluion of he HJB equaion; i is known ha srong soluions are also viscosiy soluions. This happens e.g. in he example of Secion 7. We now pass o he exisence of feedbacks. Firs we recall he definiion of (weakly) admissible feedback map. Definiion 5.4 A measurable map G : [0, T] R n U is a (weakly) admissible feedback map if for every (, x) [0, T] R n he closed loop equaion (s [, T]) y(s) = x + s F 1 (r, y (r), G(r, y (r)))dr + s F 0 (r, y (r)) dr + s B (r, y (r))dw(r), (24) admis a (weak) soluion y = y G,,x for every (, x) [0, T] R n and he conrol sraegy z G,,x (s) = G(s, y G,,x (s)) is (weakly) admissible. If we have a (weakly) admissible feedback map G hen he conrol sraegy z G,,x (s) = G(s, y G,,x (s)) is by definiion (weakly) admissible and so, given G, we have, for every (, x) [0, T] R n an admissible (weak) couple (z G,,x, y G,,x ). Our goal is o find an admissible feedback map G such ha, for every (, x) [0, T] R n, he sraegy z G,,x is (weakly) opimal. Such G will be called a (weakly) opimal feedback map. 16

18 Proposiion 5.5 Assume ha he Hypoheses of Theorem 4.9 hold rue. Assume also ha: (i) he maximum of he Hamilonian (4) exiss and i is possible o define a measurable map G 0 : [0, T] R n R n U such ha G 0 (, x, p) argmin z H CV (, x, p; z) for every (, x, p) [0, T] R n R n ; (ii) he map G(, x) = G 0 (, x, x v (, x)) is a (weakly) admissible feedback map. Then for every (, x) [0, T] R n, z G,,x is a (weakly) opimal conrol, y G,,x is he corresponding (weakly) opimal sae and v = V. Finally if arg min z U H CV (, x, v x (, x); z) is always a singleon and he closed loop equaion (24) admis a unique (weak) soluion, hen he opimal conrol is (weakly) unique. Proof. This is an obvious consequence of he verificaion Theorem 4.9. Remark 5.6 We observe ha if we have exisence of he (weak) soluion of he closed loop equaion (24) for every srong soluion of he HJB equaion (13), hen we auomaically have ha he srong soluion is unique. Remark 5.7 We finally observe ha Proposiions 5.2 and 5.5 holds also in he infinie horizon case wih obvious changes. We will use he infinie horizon version of Proposiion 5.5 in he sudy of he one dimensional example in Secion 7. 6 Applicaion 1: a class of exi ime conrol problems wih non degenerae diffusion Here we show ha our echnique for proving verificaion heorems works for a family of sochasic opimal conrol problems wih exi ime and nondegenerae diffusion where he soluions of he associaed HJB equaion are no known o be C 1, The problem Le firs R n be he sae space, R n be he space of noises and U (a given subse of a Polish space) be he conrol space. Le T be a fixed finie horizon, (Ω, F, P) be a given sochasic basis (where F sands for a given filraion (F s ) s [0,T] saisfying he usual condiions), W be a cylindrical Brownian moion wih values in R n adaped o (F s ) s [0,T]. Consider a sochasic conrolled sysem in R n wih fixed finie horizon T and iniial ime [0, T) governed by he sae equaion { dy (s) = [F0 (s, y (s)) + F 1 (s, y (s),z (s))] ds + B (y (s))dw(s), s [, T] (25) y () = x. We consider an open bounded domain O R n wih uniformly C 2 boundary (see e.g. [30, pp.2-3] for he definiion). The iniial daum x belongs o O and we call τ O he firs exi ime of he process y from his open se, i.e. τ O (ω) = inf {s > : y (s; ω) O c }. Moreover F 0 : [0, T] Ō Rn, F 1 : [0, T] Ō U Rn, B : [0, T] Ō L(Rn ; R n ), (Ō sands for he closure of he se O) saisfy 17

19 Hypohesis 6.1 F 0 and F 1 are coninuous and bounded. Hypohesis 6.2 B is ime independen, uniformly coninuous and bounded and nondegenerae i.e. here exiss λ 0 > 0 such ha λ 1 0 ξ n [ B (x)b T (x) ] ξ i,j iξ j λ 0 ξ (, x) [0, T] Ō, ξ Rn. i,j=1 Remark 6.3 Firs noe ha he above Hypohesis 6.2 implies ha B is inverible wih bounded inverse, so we fall in he Hypohesis 4.7-(ii). Moreover i is possible o exend our resuls o he case when B is ime dependen, provided suiable regulariy condiions are saisfied. Such condiions are needed o apply he resuls of semigroup heory in he subsecion below and subsanially require ha he domain of he operaor L ()u = 1 2 Tr[B (, x)b (, x) xx u] be consan and he coefficien B be Hölder coninuous in ime, see on his [3, 4] and [30, Secion 3.1.2]. We leave aside hese assumpions o keep a more simplified seing where, anyway, C 1,2 regulariy does no hold in general. Abou he conrol sraegy z : [, T] Ω U we assume ha i belongs o a given se Z ad () of sochasic processes defined on [, T] Ω wih values in a fixed Polish space U. More precisely we will assume he following. Hypohesis 6.4 The conrol space U is a Polish space, while Z ad () is he space of all measurable processes z : [, T] Ω U adaped o he filraion F. Remark 6.5 The above seing corresponds o he so-called srong formulaion of a sochasic conrol problem. In his problem, since he sae equaion admis only weak soluions in general (see nex proposiion) we need o consider he weak formulaion, leing he sochasic basis o vary, see Remark 5.1. To avoid heavy noaion we will keep he same symbols as in he srong formulaion, as we did in Secion 5. The following Proposiion can be proved in he same way as in [37], even if he drif is random. In fac Theorem herein makes use of Girsanov ransformaion o eliminae he drif; he same can be done here. Then we can apply Corollary and Theorem of [37]. Proposiion 6.6 Assume ha Hypoheses 6.1, 6.2 hold. Then, for all z Z ad (), equaion (25) has a weak soluion y( ;, x, z) CF 0 (, T; X). This soluion is unique in he sense of probabiliy law. Remark 6.7 In he case of dimension 1, no coninuiy on he diffusion coefficiens is required, see [37] Exercise a page

20 We now consider he following sochasic opimal conrol problem wih exi ime. Minimize he cos funcional J (, x; z) = [ τo T E l (s, y (s;, x, z),z (s))ds (26) +I {τo<t }ψ (τ O, y (τ O ;, x, z)) + I {τo T }φ(y (T;, x, z)) ], over all conrols z Z ad (). Here y ( ;, x, z) is he soluion of he equaion (25) and we assume ha l, ψ, φ saisfy Hypohesis 6.8 l C 0 ( [0, T] Ō), φ C 0 ( Ō ) ψ C 1+β 2,1+β ([0, T] O) (for some β > 0), ψ (T, x) = φ(x) on O. Remark 6.9 The above hypohesis is needed o apply he resuls of semigroup heory in he subsecion below. In paricular, under Hypohesis 6.8 he operaor Lu = 1 2 Tr[B (x)b (x) xx u] (wih 0 Dirichle boundary condiions) generaes an analyic semigroup and he associaed boundary value problem in O is well posed (see [30, Ch.5]). The value funcion of his problem is defined as V (, x) = inf {J (, x; z) : z Z ad ()} (27) and a conrol z Z ad () and such ha V (, x) = J(, x; z ) is said o be opimal wih respec o he iniial ime and sae (, x). The corresponding HJB equaion is where wih v (, x) Tr [B (x) xxv(, x)b (x)] = F 0 (, x), x v (, x) + H 0 (, x, x v (, x)), [0, T], x O, v (T, x) = φ(x), x O, v (, x) = ψ (, x), [0, T], x O, being he curren value Hamilonian. (28) H 0 (, x, p) = inf z U H0 CV (, x, p; z), (29) H 0 CV (, x, p; z) = F 1 (, x, z),p + l (, x, z), Proposiion 6.10 Under Hypoheses 6.1, 6.2, 6.4, 6.8, he Hamilonian H 0 (, x, p) defined in (29) is coninuous in [0, T] Ō Rn. Moreover here exiss a consan C > 0 such ha H 0 (, x, p) H 0 (, x, q) C p q H 0 (, x, p) C (1 + p ) (30) Proof. I is enough o apply he definiion of he Hamilonian and use he coninuiy and he boundedness of F 1 and l. 19

21 6.2 Srong soluions of he HJB equaion The HJB equaion (28) above is a semilinear parabolic equaion wih coninuous coefficiens. Since he second order erm is nondegenerae (Hypohesis 6.2) one expecs inerior regulariy resuls for he soluion even if he boundary daa are merely coninuous, as i is under our assumpions. In paricular he following resul, aken from [9], Theorems 9.1 and 9.2, apply. Theorem 6.11 Under Hypoheses 6.1, 6.2, 6.4, 6.8, here exiss a viscosiy soluion u of (28) and u belongs o he space C 1+α 2,1+α loc ([0, T) O) C ( 0 [0, T] Ō) for some α (0, 1). We noe also ha i is no known if he soluion is classical in he inerior. This is known, using heorems of [5], if he daa are supposed o be Hölder coninuous in (, x), which we do no assume here. Moreover even in he linear case, if he coefficiens are only coninuous, soluions are only W 2,p, so hey are also C 1,α. An example in his direcion is in [16, Ch.4]. This means ha we are exacly in he case where i makes sense o apply our echnique o prove verificaion heorems: no classical soluion bu soluions (in a generalized sense) wih a leas C 0,1 regulariy. Once his is clear, we need o check if our soluions (ha exiss a leas in he viscosiy sense hanks o he above Theorem 6.11) are also srong soluions in he sense of Definiion 4.6 (suiably modified o ake care of he boundary daum ψ). Such a resul is no available in he lieraure in his form bu i can be easily deduced using he resuls on analyic semigroups conained e.g. in [30]. We explain here below how his can be done in he case when ψ = 0 (he general case can be reaed wih he ideas explained in [30, Secion 5.1.2]). In his case he operaor L : D (L) C 0 ( Ō ) C 0 ( Ō ), { D (L) = η p 1 W 2,p loc (O) : η, Lη ( C0 Ō ) }, η O = 0, (Lη)(x) = 1 2 Tr[B (x) xxη(x)b (x)], generaes an analyic semigroup { e L, 0 }, [30, p.97, Corollary ]. Moreover, given any iniial daum φ C ( 0 Ō ) and any funcion H 0 coninuous in [0, T] Ō Rn and saisfying (30) we can apply a modificaion of he argumen used o prove Proposiion in [30, p.281] o ge exisence and uniqueness of a soluion u of he inegral equaion ( ) u (, x) = e (T )L φ (x) (31) + ( e (T s)l [ F 0 (s, ), x u (s, ) + H 0 (s,, x u (s, )) ]) (x) ds, which can be considered as an inegral form of he PDE (28) when ψ = 0 (in he general case one needs o lif he funcion ψ ino he equaion obaining an exra erm in (31), see on his [30, Remark , p.195]) wrien using he variaion of consans formula. More precisely we have he following definiions and resuls. Definiion 6.12 We say ha a funcion w belongs o he space Σ ( 1,α [0, T] Ō) if w C ( 0 [0, T] Ō), w is Fréche differeniable in x O, x w C ( 0 [0, T ε] Ō) for every ε (0, T), and sup (T ) α x w (, x) < +. (,x) [0,T) Ō ( Definiion 6.13 A funcion u Σ 1, 1 2 [0, T] Ō ) ha saisfies he inegral equaion (31) for every (, x) [0, T] Ō is called a mild soluion of he HJB equaion (28). 20

22 Theorem 6.14 Assume ha B C ( 0 Ō ), F 0 C ( 0 [0, T] Ō), φ C ( 0 Ō ) and ψ = 0. Moreover le H 0 be coninuous ( in [0, T] Ō Rn and saisfies (30). Then here exiss a unique mild soluion u Σ 1, 1 2 [0, T] Ō ) of he HJB equaion (28). Proof. The proof is a sandard applicaion of he conracion mapping principle, see [30, Proposiion 7.3.4] and also [6, 18]. Now we wan o show ha such a mild soluion is a srong soluion in he sense ha i can be seen as he limi of smooh soluions. We rewrie here he definiion of srong soluion for our case since i is slighly weaker han he one used in Secion 4 due o he possible singulariy of he spaial gradien a = 0, see Remark 4.8. Definiion 6.15 Le φ C ( 0 Ō ) ( and ψ C 0 ([0, T] O). A funcion u Σ 1, 1 2 [0, T] Ō ) is a srong soluion of he equaion (28) if i saisfies he boundary and final condiions u (T, x) = φ(x), x O, u (, x) = ψ (, x), [0, T],s O, and if, for every ε > 0, i is a srong soluion (in he sense of Definiion 4.2) of he linear parabolic problem w + Lw = h, [0, T], x O, w (T ε, x) = u (T ε, x), x O, where h = [ F 0 (, ), x u (, ) + H 0 (,, x u (, )) ], i.e. if, for every ε > 0 here exis sequences u ε n, hε n, φε n such ha 1. for every n u ε n is a sric soluion (i.e. i belong o ( C1,2 [0, T ε] Ō) and saisfy he equaliies on [0, T ε] Ō) of he approximaing problem v + Lv = h ε n, v (T ε, x) = φε n (x) ; 2. u ε n converges o u uniformly in [0, T ε] Ō, as n + ; 3. h ε n converges o h uniformly in [0, T ε] Ō, as n + ; 4. φ ε n u (T ε, ) converges o zero uniformly in Ō, as n +. Given he above Definiion 6.15 we can apply direcly Proposiion (see also Theorem ) of [30] o ge he following. Theorem 6.16 Under he same assumpions of Theorem 6.14 he mild soluion of (28) is also srong. Proof. I is enough o apply Proposiion and Theorem of [30] o he nonhomogenous linear parabolic problem v + Lv = h, v (T, x) = u (T ε, x), where we se h (, x) = F 0 (, x), x u (, x) + H 0 (, x, x u (, x)). In fac for every ε > 0 such funcion h belongs o he space C 0 ( [0, T ε] ; C 0 ( Ō )) = C 0 ( [0, T ε] Ō). 21

23 Remark 6.17 In he case when ψ is no 0 he above resuls sill hold rue using he same echniques shown in [30], Theorem , and We do no do i here for simpliciy of exposiion. Remark 6.18 Suppose ha he operaor L can be wrien in divergence form and ha he coefficiens are only Borel measurable and no necessarily coninuous; suppose moreover ha he diffusion coefficiens are lower and upper bounded by a consan. Then, he semigroup has a densiy wih respec o he Lebesgue measure and i fulfills he classical Aronson esimaes, see for insance [1, 36]. Fukushima - Dirichle decomposiion for mild or weak soluions o equaions of ype (28) were reaed by [2, 26, 35]. Using such kind of resuls or possible generalizaions in he spiri of [22], one could esablish verificaion heorems relaed o opimal conrol problems even in ha framework. 6.3 The verificaion heorem Now we prove a verificaion heorem for he opimal conrol problem above (25) - (26). The proof is a modificaion of he proof given in Secion 4. The main differences are due o 1. a singulariy of he firs derivaive and so a differen definiion of srong soluion; 2. he consrain on he se O and so he presence of boundary daa in x (exi ime). The saemen of he verificaion heorem in his case is he following ( Theorem 6.19 Assume ha Hypoheses 6.1, 6.2, 6.4, 6.8, hold. Le v Σ 1, 1 2 [0, T] Ō ) be a srong soluion of (28). Then (i) v V on [0, T] Ō. (ii) If z is an admissible conrol a (, x) [0, T] O ha saisfies (seing y (s) = y (s;, x, z)), P-a.s., H 0 (s, y (s), x v (s, y (s))) = H 0 CV (s, y (s), xv (s, y (s));z (s)), for a.e. s [, T τ O ], hen z is opimal a (, x) and v (, x) = V (, x). The proof of his heorem follows as usual by he following fundamenal ideniy ha we sae as a lemma. ( Lemma 6.20 Assume ha Hypoheses 6.1, 6.2, 6.4, 6.8, hold. Le v Σ 1,1 2 [0, T] Ō ) be a srong soluion of (28). Then, for every (, x) [0, T] Ō and z Z ad () he following ideniy holds J (, x; z) = v (, x) τo T [ +E H 0 (s, y(s), x v (s, y (s))) + HCV 0 (s, y(s), x v (s, y (s));z(s)) ] ds, (32) where y(s) def = y (s;, x, z) is he soluion of (25) associaed wih he conrol z. 22

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

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

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

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

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

arxiv: v1 [math.pr] 28 Nov 2016

arxiv: v1 [math.pr] 28 Nov 2016 Backward Sochasic Differenial Equaions wih Nonmarkovian Singular Terminal Values Ali Devin Sezer, Thomas Kruse, Alexandre Popier Ocober 15, 2018 arxiv:1611.09022v1 mah.pr 28 Nov 2016 Absrac We solve a

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

Backward Stochastic Differential Equations with Nonmarkovian Singular Terminal Values

Backward Stochastic Differential Equations with Nonmarkovian Singular Terminal Values Backward Sochasic Differenial Equaions wih Nonmarkovian Singular Terminal Values Ali Sezer, Thomas Kruse, Alexandre Popier, Ali Sezer To cie his version: Ali Sezer, Thomas Kruse, Alexandre Popier, Ali

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

Semilinear Kolmogorov equations and applications to stochastic optimal control

Semilinear Kolmogorov equations and applications to stochastic optimal control Semilinear Kolmogorov equaions and applicaions o sochasic opimal conrol Federica Masiero 1 Advisor: Prof. Marco Fuhrman 2 1 Diparimeno di Maemaica, Universià degli sudi di Milano, Via Saldini 5, 2133 Milano,

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

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

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

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

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

Convergence of the Neumann series in higher norms

Convergence of the Neumann series in higher norms Convergence of he Neumann series in higher norms Charles L. Epsein Deparmen of Mahemaics, Universiy of Pennsylvania Version 1.0 Augus 1, 003 Absrac Naural condiions on an operaor A are given so ha he Neumann

More information

Homogenization of random Hamilton Jacobi Bellman Equations

Homogenization of random Hamilton Jacobi Bellman Equations Probabiliy, Geomery and Inegrable Sysems MSRI Publicaions Volume 55, 28 Homogenizaion of random Hamilon Jacobi Bellman Equaions S. R. SRINIVASA VARADHAN ABSTRACT. We consider nonlinear parabolic equaions

More information

A remark on the H -calculus

A remark on the H -calculus A remark on he H -calculus Nigel J. Kalon Absrac If A, B are secorial operaors on a Hilber space wih he same domain range, if Ax Bx A 1 x B 1 x, hen i is a resul of Auscher, McInosh Nahmod ha if A has

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

The Optimal Stopping Time for Selling an Asset When It Is Uncertain Whether the Price Process Is Increasing or Decreasing When the Horizon Is Infinite

The Optimal Stopping Time for Selling an Asset When It Is Uncertain Whether the Price Process Is Increasing or Decreasing When the Horizon Is Infinite American Journal of Operaions Research, 08, 8, 8-9 hp://wwwscirporg/journal/ajor ISSN Online: 60-8849 ISSN Prin: 60-8830 The Opimal Sopping Time for Selling an Asse When I Is Uncerain Wheher he Price Process

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

Hamilton Jacobi equations

Hamilton Jacobi equations Hamilon Jacobi equaions Inoducion o PDE The rigorous suff from Evans, mosly. We discuss firs u + H( u = 0, (1 where H(p is convex, and superlinear a infiniy, H(p lim p p = + This by comes by inegraion

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

arxiv: v1 [math.pr] 6 Oct 2008

arxiv: v1 [math.pr] 6 Oct 2008 MEASURIN THE NON-STOPPIN TIMENESS OF ENDS OF PREVISIBLE SETS arxiv:8.59v [mah.pr] 6 Oc 8 JU-YI YEN ),) AND MARC YOR 3),4) Absrac. In his paper, we propose several measuremens of he nonsopping imeness of

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

Expert Advice for Amateurs

Expert Advice for Amateurs Exper Advice for Amaeurs Ernes K. Lai Online Appendix - Exisence of Equilibria The analysis in his secion is performed under more general payoff funcions. Wihou aking an explici form, he payoffs of he

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

Lecture 10: The Poincaré Inequality in Euclidean space

Lecture 10: The Poincaré Inequality in Euclidean space Deparmens of Mahemaics Monana Sae Universiy Fall 215 Prof. Kevin Wildrick n inroducion o non-smooh analysis and geomery Lecure 1: The Poincaré Inequaliy in Euclidean space 1. Wha is he Poincaré inequaliy?

More information

Singular control of SPDEs and backward stochastic partial diffe. reflection

Singular control of SPDEs and backward stochastic partial diffe. reflection Singular conrol of SPDEs and backward sochasic parial differenial equaions wih reflecion Universiy of Mancheser Join work wih Bern Øksendal and Agnès Sulem Singular conrol of SPDEs and backward sochasic

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

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

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

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

Essential Maps and Coincidence Principles for General Classes of Maps

Essential Maps and Coincidence Principles for General Classes of Maps Filoma 31:11 (2017), 3553 3558 hps://doi.org/10.2298/fil1711553o Published by Faculy of Sciences Mahemaics, Universiy of Niš, Serbia Available a: hp://www.pmf.ni.ac.rs/filoma Essenial Maps Coincidence

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

Utility maximization in incomplete markets

Utility maximization in incomplete markets Uiliy maximizaion in incomplee markes Marcel Ladkau 27.1.29 Conens 1 Inroducion and general seings 2 1.1 Marke model....................................... 2 1.2 Trading sraegy.....................................

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

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

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013 IMPLICI AND INVERSE FUNCION HEOREMS PAUL SCHRIMPF 1 OCOBER 25, 213 UNIVERSIY OF BRIISH COLUMBIA ECONOMICS 526 We have exensively sudied how o solve sysems of linear equaions. We know how o check wheher

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

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

Empirical Process Theory

Empirical Process Theory Empirical Process heory 4.384 ime Series Analysis, Fall 27 Reciaion by Paul Schrimpf Supplemenary o lecures given by Anna Mikusheva Ocober 7, 28 Reciaion 7 Empirical Process heory Le x be a real-valued

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

Uniqueness of solutions to quadratic BSDEs. BSDEs with convex generators and unbounded terminal conditions

Uniqueness of solutions to quadratic BSDEs. BSDEs with convex generators and unbounded terminal conditions Recalls and basic resuls on BSDEs Uniqueness resul Links wih PDEs On he uniqueness of soluions o quadraic BSDEs wih convex generaors and unbounded erminal condiions IRMAR, Universié Rennes 1 Châeau de

More information

Algorithmic Trading: Optimal Control PIMS Summer School

Algorithmic Trading: Optimal Control PIMS Summer School Algorihmic Trading: Opimal Conrol PIMS Summer School Sebasian Jaimungal, U. Torono Álvaro Carea,U. Oxford many hanks o José Penalva,(U. Carlos III) Luhui Gan (U. Torono) Ryan Donnelly (Swiss Finance Insiue,

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

A New Perturbative Approach in Nonlinear Singularity Analysis

A New Perturbative Approach in Nonlinear Singularity Analysis Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The

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

A Note on Superlinear Ambrosetti-Prodi Type Problem in a Ball

A Note on Superlinear Ambrosetti-Prodi Type Problem in a Ball A Noe on Superlinear Ambrosei-Prodi Type Problem in a Ball by P. N. Srikanh 1, Sanjiban Sanra 2 Absrac Using a careful analysis of he Morse Indices of he soluions obained by using he Mounain Pass Theorem

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

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

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

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

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

FEEDBACK NULL CONTROLLABILITY OF THE SEMILINEAR HEAT EQUATION

FEEDBACK NULL CONTROLLABILITY OF THE SEMILINEAR HEAT EQUATION Differenial and Inegral Equaions Volume 5, Number, January 2002, Pages 5 28 FEEDBACK NULL CONTROLLABILITY OF THE SEMILINEAR HEAT EQUATION Mihai Sîrbu Deparmen of Mahemaical Sciences, Carnegie Mellon Universiy

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

On Gronwall s Type Integral Inequalities with Singular Kernels

On Gronwall s Type Integral Inequalities with Singular Kernels Filoma 31:4 (217), 141 149 DOI 1.2298/FIL17441A Published by Faculy of Sciences and Mahemaics, Universiy of Niš, Serbia Available a: hp://www.pmf.ni.ac.rs/filoma On Gronwall s Type Inegral Inequaliies

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

Differential Harnack Estimates for Parabolic Equations

Differential Harnack Estimates for Parabolic Equations Differenial Harnack Esimaes for Parabolic Equaions Xiaodong Cao and Zhou Zhang Absrac Le M,g be a soluion o he Ricci flow on a closed Riemannian manifold In his paper, we prove differenial Harnack inequaliies

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

Some New Uniqueness Results of Solutions to Nonlinear Fractional Integro-Differential Equations

Some New Uniqueness Results of Solutions to Nonlinear Fractional Integro-Differential Equations Annals of Pure and Applied Mahemaics Vol. 6, No. 2, 28, 345-352 ISSN: 2279-87X (P), 2279-888(online) Published on 22 February 28 www.researchmahsci.org DOI: hp://dx.doi.org/.22457/apam.v6n2a Annals of

More information

QUANTITATIVE DECAY FOR NONLINEAR WAVE EQUATIONS

QUANTITATIVE DECAY FOR NONLINEAR WAVE EQUATIONS QUANTITATIVE DECAY FOR NONLINEAR WAVE EQUATIONS SPUR FINAL PAPER, SUMMER 08 CALVIN HSU MENTOR: RUOXUAN YANG PROJECT SUGGESTED BY: ANDREW LAWRIE Augus, 08 Absrac. In his paper, we discuss he decay rae for

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

EXISTENCE AND UNIQUENESS THEOREMS ON CERTAIN DIFFERENCE-DIFFERENTIAL EQUATIONS

EXISTENCE AND UNIQUENESS THEOREMS ON CERTAIN DIFFERENCE-DIFFERENTIAL EQUATIONS Elecronic Journal of Differenial Equaions, Vol. 29(29), No. 49, pp. 2. ISSN: 72-669. URL: hp://ejde.mah.xsae.edu or hp://ejde.mah.un.edu fp ejde.mah.xsae.edu EXISTENCE AND UNIQUENESS THEOREMS ON CERTAIN

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

Notes for Lecture 17-18

Notes for Lecture 17-18 U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up

More information

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

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

Existence Theory of Second Order Random Differential Equations

Existence Theory of Second Order Random Differential Equations Global Journal of Mahemaical Sciences: Theory and Pracical. ISSN 974-32 Volume 4, Number 3 (22), pp. 33-3 Inernaional Research Publicaion House hp://www.irphouse.com Exisence Theory of Second Order Random

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

CHARACTERIZATION OF REARRANGEMENT INVARIANT SPACES WITH FIXED POINTS FOR THE HARDY LITTLEWOOD MAXIMAL OPERATOR

CHARACTERIZATION OF REARRANGEMENT INVARIANT SPACES WITH FIXED POINTS FOR THE HARDY LITTLEWOOD MAXIMAL OPERATOR Annales Academiæ Scieniarum Fennicæ Mahemaica Volumen 31, 2006, 39 46 CHARACTERIZATION OF REARRANGEMENT INVARIANT SPACES WITH FIXED POINTS FOR THE HARDY LITTLEWOOD MAXIMAL OPERATOR Joaquim Marín and Javier

More information

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3

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

Asymptotic instability of nonlinear differential equations

Asymptotic instability of nonlinear differential equations Elecronic Journal of Differenial Equaions, Vol. 1997(1997), No. 16, pp. 1 7. ISSN: 172-6691. URL: hp://ejde.mah.sw.edu or hp://ejde.mah.un.edu fp (login: fp) 147.26.13.11 or 129.12.3.113 Asympoic insabiliy

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

arxiv: v1 [math.pr] 23 Jan 2019

arxiv: v1 [math.pr] 23 Jan 2019 Consrucion of Liouville Brownian moion via Dirichle form heory Jiyong Shin arxiv:90.07753v [mah.pr] 23 Jan 209 Absrac. The Liouville Brownian moion which was inroduced in [3] is a naural diffusion process

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

Online Appendix to Solution Methods for Models with Rare Disasters

Online Appendix to Solution Methods for Models with Rare Disasters Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,

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

arxiv: v1 [math.ca] 15 Nov 2016

arxiv: v1 [math.ca] 15 Nov 2016 arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy

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

Some Regularity Properties of Three Dimensional Incompressible Magnetohydrodynamic Flows

Some Regularity Properties of Three Dimensional Incompressible Magnetohydrodynamic Flows Global Journal of Pure and Applied Mahemaics. ISSN 973-78 Volume 3, Number 7 (7), pp. 339-335 Research India Publicaions hp://www.ripublicaion.com Some Regulariy Properies of Three Dimensional Incompressible

More information

THE BELLMAN PRINCIPLE OF OPTIMALITY

THE BELLMAN PRINCIPLE OF OPTIMALITY THE BELLMAN PRINCIPLE OF OPTIMALITY IOANID ROSU As I undersand, here are wo approaches o dynamic opimizaion: he Ponrjagin Hamilonian) approach, and he Bellman approach. I saw several clear discussions

More information

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.

More information

CONTRIBUTION TO IMPULSIVE EQUATIONS

CONTRIBUTION TO IMPULSIVE EQUATIONS European Scienific Journal Sepember 214 /SPECIAL/ ediion Vol.3 ISSN: 1857 7881 (Prin) e - ISSN 1857-7431 CONTRIBUTION TO IMPULSIVE EQUATIONS Berrabah Faima Zohra, MA Universiy of sidi bel abbes/ Algeria

More information

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems Essenial Microeconomics -- 6.5: OPIMAL CONROL Consider he following class of opimizaion problems Max{ U( k, x) + U+ ( k+ ) k+ k F( k, x)}. { x, k+ } = In he language of conrol heory, he vecor k is he vecor

More information

A Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs

A Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs PROC. IEEE CONFERENCE ON DECISION AND CONTROL, 06 A Primal-Dual Type Algorihm wih he O(/) Convergence Rae for Large Scale Consrained Convex Programs Hao Yu and Michael J. Neely Absrac This paper considers

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

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

Rare Event Simulation with Vanishing Error for Small Noise Diffusions

Rare Event Simulation with Vanishing Error for Small Noise Diffusions Rare Even Simulaion wih Vanishing Error for Small Noise Diffusions Eric Vanden-Eijnden, Jonahan Weare Augus 7, 2009 Absrac We consruc an imporance sampling mehod for cerain rare even problems involving

More information

On Oscillation of a Generalized Logistic Equation with Several Delays

On Oscillation of a Generalized Logistic Equation with Several Delays Journal of Mahemaical Analysis and Applicaions 253, 389 45 (21) doi:1.16/jmaa.2.714, available online a hp://www.idealibrary.com on On Oscillaion of a Generalized Logisic Equaion wih Several Delays Leonid

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

Boundedness and Exponential Asymptotic Stability in Dynamical Systems with Applications to Nonlinear Differential Equations with Unbounded Terms

Boundedness and Exponential Asymptotic Stability in Dynamical Systems with Applications to Nonlinear Differential Equations with Unbounded Terms Advances in Dynamical Sysems and Applicaions. ISSN 0973-531 Volume Number 1 007, pp. 107 11 Research India Publicaions hp://www.ripublicaion.com/adsa.hm Boundedness and Exponenial Asympoic Sabiliy in Dynamical

More information

L. Tentarelli* ON THE EXTENSIONS OF THE DE GIORGI APPROACH TO NONLINEAR HYPERBOLIC EQUATIONS

L. Tentarelli* ON THE EXTENSIONS OF THE DE GIORGI APPROACH TO NONLINEAR HYPERBOLIC EQUATIONS Rend. Sem. Ma. Univ. Poliec. Torino Bruxelles-Torino Talks in PDE s Turin, May 5, 06 Vol. 74, (06), 5 60 L. Tenarelli* ON THE EXTENSIONS OF THE DE GIORGI APPROACH TO NONLINEAR HYPERBOLIC EQUATIONS Absrac.

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

Singular perturbation control problems: a BSDE approach

Singular perturbation control problems: a BSDE approach Singular perurbaion conrol problems: a BSDE approach Join work wih Francois Delarue Universié de Nice and Giuseppina Guaeri Poliecnico di Milano Le Mans 8h of Ocober 215 Conference in honour of Vlad Bally

More information