arxiv: v4 [math.pr] 20 Jul 2016

Transcription

1 Submied o he Aals of Applied Pobabiliy ε-strong SIMULATION FOR MULTIDIMENSIONAL STOCHASTIC DIFFERENTIAL EQUATIONS VIA ROUGH PATH ANALYSIS axiv: v4 [mah.pr] 20 Jul 2016 By Jose Blache, Xiyu Che ad Jig Dog Columbia Uivesiy ad Wuha Uivesiy ad Nohwese Uivesiy Coside a mulidimesioal diffusio pocess X = {X : [0, 1]}. Le ε > 0 be a deemiisic, use defied, oleace eo paamee. Ude sadad egulaiy codiios o he dif ad diffusio coefficies of X, we cosuc a pobabiliy space, suppoig boh X ad a explici, piecewise cosa, fully simulaable pocess X ε such ha sup X ε X < ε 0 1 wih pobabiliy oe. Moeove, he use ca adapively choose ε 0, ε so ha X ε also piecewise cosa ad fully simulaable ca be cosuced codiioal o X ε o esue a eo smalle ha ε wih pobabiliy oe. Ou cosucio equies a deailed sudy of coiuiy esimaes of he Iô map usig Lyos heoy of ough pahs. We appoximae he udelyig Bowia moio, joily wih he Lévy aeas wih a deemiisic ε eo i he udelyig ough pah meic. 1. Ioducio. Coside he Iô Sochasic Diffeeial Equaio SDE 1.1 dx = µxd σxdz, X0 = x0 whee Z is a d -dimesioal Bowia moio, ad µ : R d R d ad σ : R d R d d saisfy suiable egulaiy codiios. We shall assume, i paicula, ha boh µ ad σ ae Lipschiz coiuous so ha a sog soluio o he SDE is guaaeed o exis. Addiioal assumpios o he fis ad secod ode deivaives of µ ad σ, which ae sadad i he heoy of ough pahs, will be discussed i he sequel. Ou coibuio i his pape is he joi cosucio of X = {X : [0, 1]} ad a family of pocesses X ε = {X ε : [0, 1]}, fo each ε 0, 1, suppoed o a pobabiliy space Ω, F, P, ad such ha he followig popeies hold: MSC 2010 subjec classificaios: Pimay 34K50, 65C05, 82B80; secoday 97K60 Keywods ad phases: sochasic diffeeial equaio, Moe Calo mehod, Bowia Moio, Lévy aea, Rough pah 1

2 2 J BLANCHET, X. CHEN AND J. DONG T1 The pocess X ε is piecewise cosa, wih fiiely may discoiuiies i [0, 1]. T2 The pocess X ε ca be simulaed exacly ad, sice i akes oly fiiely may values, is pah ca be fully soed. T3 We have ha wih P -pobabiliy oe 1.2 sup X ε X < ε. [0,1] T4 Fo ay m > 1 ad 0 < ε m <... < ε 1 < 1 we ca simulae X εm codiioal o X ε1,...,x εm 1. We efe o he class of pocedues which achieve he cosucio of such family {X ε : ε 0, 1} as Toleace-Efoced Simulaio TES o ε-sog simulaio mehods. Thoughou he pape we use o deoe he max-om o R d. This pape povides he fis cosucio of a Toleace-Efoced Simulaio pocedue fo mulidimesioal SDEs i subsaial geealiy. All ohe TES o ε-sog simulaio pocedues up o ow ae applicable o oe dimesioal pocesses o mulidimesioal pocesses wih cosa diffusio maix i.e. σ x = σ. Le us discuss some cosideaios ha moivae ou sudy. We fis discuss how his pape elaes o he cue lieaue o ε-sog simulaio of sochasic pocesses, which is a ece aea of eseach. The pape of [6] povides he cosucio of X ε saisfyig oly T1 o T3, i oe dimesio. I paicula, boud 1.2 is saisfied fo a give fixed ε 0 = ε > 0 ad i is o clea how o joily simulae {X εm } m 1 as ε m 0 applyig he echique i [6]. The moivaio of cosucig X ε0 fo [6] came fom he desie o poduce exac samples fom a oe dimesioal diffusio X saisfyig 1.1, ad also assumig σ cosa. The auhos i [6] wee ieesed i exedig he applicabiliy of a algoihm ioduced by Beskos ad Robes, see [2]. The pocedue of Beskos ad Robes, applicable o oe dimesioal diffusios, imposed sog boudedess assumpios o he dif coefficie ad is deivaive. The echique i [6] eabled a exesio which is fee of such boudedess assumpios by usig a localizaio echique ha allowed o apply he ideas behid he algoihm i [2]; see also [3] fo aohe appoach which elimiaes boudedess assumpios. All of hese developmes ae i he oe-dimesioal case. The assumpio of a cosa diffusio coefficie comes a basically o cos i geealiy whe cosideig oe dimesioal diffusios because oe ca always apply Lampei s oe-o-oe asfomaio. Such asfoma-

3 ε-strong SIMULATION FOR SDES 3 io allows o ecas he simulaio poblem o oe ivolvig a diffusio wih cosa σ. Lampei s asfomaio cao be geeally applied i highe dimesios. The pape of [4] exeds he wok of [6] i ha hei algoihms saisfy T1 o T4, bu also i he coex of oe dimesioal pocesses. The pape [11] o oly povides a addiioal exesio which allows o deal wih oe dimesioal SDEs wih jumps, bu also coais a compehesive discussio o exac ad ε-sog simulaio fo SDEs. Popey T4 i he defiiio of TES is desiable because i povides aohe appoach a cosucig ubiased esimaos fo expecaios of he fom Ef X, whee f is, say, a coiuous fucio of he sample pah X. I ode o see his, le us assume fo simpliciy ha f is posiive ad Lipschiz coiuous i he uifom om wih Lipschiz cosa K. The, le T be ay posiive adom vaiable wih a sicly posiive desiy g o [0, ad defie 1.3 Z := I f X > T /g T. Obseve ha [ E[Z] = E[E [Z X]] = E 0 I f X > g ] g d = E[f X], so Z is a ubiased esimao fo Ef X. Theefoe, if Popeies T1 o T4 hold, i is possible o simulae Z by oig ha f X ε > T Kε implies f X > T ad if f X ε < T Kε, he f X T. Sice T4 allows o keep simulaig as ε becomes smalle ad T is idepede of X ε wih a posiive desiy g, he oe eveually is able o simulae Z exacly. The majo obsacle ivolved i developig exac samplig algoihms fo mulidimesioal diffusios is he fac ha σ cao be assumed o be cosa. Moeove, eve i he case of mulidimesioal diffusios wih cosa σ, he oe dimesioal algoihms developed so fa ca oly be exeded o he case i which he dif coefficie µ is he gadie of some fucio, ha is, if µ x = v x fo some v. The easo is ha i his case oe ca epese he likelihood aio L, bewee he soluio o 1.1 ad Bowia moio assumig σ = I fo simpliciy ivolvig a Riema iegal as follows 1.4 L = exp µ X s dx s = exp v X exp v X 0 exp µ X s 2 2 ds 0 λ X s ds,

4 4 J BLANCHET, X. CHEN AND J. DONG fo λ x = v x v x 2 2. The fac ha he sochasic iegal ca be asfomed io a Riema iegal faciliaes he execuio of accepaceejecio because oe ca iepe up o a cosa ad usig localizaio as i [6] he expoeial of he iegal of λ as he pobabiliy ha o aivals occu i a Poisso pocess wih a sochasic iesiy. Such eve i.e. o aivals ca be simulaed by hiig. So, ou moivaio i his pape is o ivesigae a ovel appoach ha allows o sudy ε-sog simulaio fo mulidimesioal diffusios i subsaial geealiy, wihou imposig he assumpio ha σ is cosa o ha a Lampei-ype asfomaio ca be applied. Give he pevious discussio o he coecios bewee exac samplig ad ε-sog simulaio, ad he limiaios of he cue echiques, we believe ha ou esuls hee povide a impoa sep i he developme of exac samplig algoihms fo geeal mulidimesioal diffusios. Fo example, i coas o exisig echiques, which demad L o be expessed i ems of a Riema iegal as idicaed i 1.4, ou esuls hee allow o appoximae diecly L i ems of he sochasic iegal epeseaio ad hus oe does o eed o assume ha µ x = v x. We pla o epo o hese implicaios i fuue papes. Ou esuls aleady allow o obai ubiased esimao of expecaios of sample pah fucioals via 1.3. Howeve, i is oed i [4] ha he expeced umbe of adom vaiables equied o simulae Z is ypically ifiie. The ece pape [11] discusses via umeical examples he pacical limiaios of hese ypes of esimaos. The wok of [12], also poposes ubiased esimaos fo he expecaio of Lipschiz coiuous fucios of X1 usig adomized mulilevel Moe Calo. Neveheless, hei algoihm also exhibis ifiie expeced emiaio ime, excep whe oe ca simulae he Lévy aeas exacly, which cuely ca be doe oly i he coex of wo dimesioal SDEs usig he esuls i [9]. The auhos i [1] also use ough pah aalysis fo Moe Calo esimaio, bu hei focus is o coecios o mulilevel echiques ad o o ε-sog simulaio. I his pape we coceae oly o wha is possible o do i ems of ε-sog simulaio pocedues ad how o eable he use of ough pah heoy fo ε-sog simulaio. We shall sudy efficie implemeaios of he algoihms poposed i a sepaae pape. Ohe eseach aveues ha we pla o ivesigae, ad which leveage off ou developme i his pape, ivolve quaificaio of model uceaiy usig he fac ha ou ε-sog simulaio algoihms i he ed ae uifom fo cases wih a lage class of dif ad diffusio coefficies.

5 ε-strong SIMULATION FOR SDES 5 Fially, we oe ha i ode o build ou Toleace-Efoced Simulaio pocedue we had o obai ew ools fo he aalysis of Lévy aeas ad associaed codiioal lage deviaios esuls fo Lévy aeas give he icemes of Bowia moio. We believe ha hese echical esuls migh be of idepede iees. The es of he pape is ogaized as follows. I Secio 2 we descibe he wo mai esuls of he pape. The fis of hem, Theoem 2.1, povides a eo boud bewee he soluio o he SDE descibed i 1.1 ad a suiable piecewise cosa appoximaio. The secod esul, Theoem 2.2, efes o he pocedues ha ae ivolved i simulaig he bouds, joily wih he piecewise cosa appoximaio, heeby yieldig 1.2. Secio 3 is divided io wo subsecios ad i builds he elemes behid he poof of Theoem 2.2. As i us ou, oe eeds o simulae bouds o he so-called Hölde oms of he udelyig Bowia moio ad he coespodig Lévy aeas. Secio 4 lays ou he deails of he simulaio of he Bowia moio ad a uppe boud of is α-hölde om ad Secio 5 lays ou he deails of he simulaio of he Lévy aeas ad a uppe boud of is 2α-Hölde om. Secio 6 is also divided i seveal pas, coespodig o he elemes of ough pah heoy equied o aalyze he SDE descibed i 1.1 as a coiuous map of Bowia moio ude a suiable meic descibed i Secio 2. While he fial fom of he esimaes i Secio 6 migh be somewha diffee ha hose obaied i he lieaue o ough pah aalysis, he echiques ha we use hee ae ceaily sadad i ha lieaue. We have chose o pese he deails because he echiques migh o be well kow o he Moe Calo simulaio commuiy ad also because ou emphasis is i fidig explici cosas i.e. bouds ha ae ameable o simulaio. 2. Mai Resuls. Ou appoach cosiss i sudyig he pocess X as a asfomaio of he udelyig Bowia moio Z. Such asfomaio is kow as he Iô-Lyos map ad is coiuiy popeies ae sudied i he heoy of ough pahs, pioeeed by T. Lyos, i [10]. A ough pah is a effecive way o summaize a iegula pah ifomaio. The heoy of ough pahs allows o defie he soluio o a SDE such as 1.1 i a pahby-pah basis fee of pobabiliy by imposig cosais o he egulaiy of he ieaed iegals of he udelyig pocess Z. Namely, iegals of he fom 2.1 A s, = s Z i u Z i s dz j u. The heoy esuls i diffee iepeaios of he soluio o 1.1

6 6 J BLANCHET, X. CHEN AND J. DONG depedig o how he ieaed iegals of Z ae iepeed. I his pape, we iepe he iegal i 2.1 i he sese of Iô. I us ou ha he Iô-Lyos map is coiuous ude a suiable α- Hölde meic defied i he space of ough pahs. I paicula, such meic ca be expessed as he maximum of he followig wo quaiies: Z Zs Z α := sup 0 s< 1 s α, A s, A 2α := sup max 1 d s 2α. 0 s< 1 As we shall discuss, coiuiy esimaes of he Iô-Lyos map ca be give explicily i ems of hese wo quaiies. I he case of Bowia moio, as we coside hee, we have ha α 1/3, 1/2. I is show i [7], ha ude suiable egulaiy codiios o µ ad σ, which we shall discuss momeaily, he Eule scheme povides a almos sue appoximaio i uifom om o he soluio o he SDE 1.1. Ou fis esul povides a explici chaaceizaio of all of he pahdepede quaiies ha ae ivolved i he fial eo aalysis such as Z α ad A 2α, he diffeece bewee ou aalysis ad wha has bee doe i pevious developmes is ha ulimaely we mus be able o impleme he Eule scheme joily wih he pah-depede quaiies ha ae ivolved i he eo aalysis. So, i is o sufficie o ague ha hee exiss a pah-depede cosa ha seves as a boud of some so, we acually mus povide a suiable epeseaio ha ca be simulaed i fiie ime. I ode o povide ou fis esul, we ioduce some oaios. Le D deoe he dyadic disceizaio of ode ad deoe he mesh of he disceizaio. Specifically, D := { 0, 1,..., 2 } whee k = k/2 fo k = 0, 1, 2,..., 2 ad = 1/2. Give ˆX 0 = x0, defie { ˆX : D } by he followig ecusio: ˆX i k1 = ˆX i k µ i ˆX k σ ˆX k Z j k1 Z j k 2.4 d d d j=1 l=1 m=1 d j=1 l σ ˆX k σ l,m ˆX k à m,j k, k1, whee à i,i k, k1 = A i,i k, k1 = Z i k1 Z i k 2 /2 /2, ad à k, k1 = 0 fo i j. We le ˆX = ˆX whee = max{ k :

7 ε-strong SIMULATION FOR SDES 7 k } fo [0, 1]. We deoe R l, m := m k=l1 ad fo fixed β 1 α, 2α, wie Γ R := sup sup {A k 1, k à i,i k, k1 }. max 0 s< 1,s, D 1 d R s, s β 2α β We oice ha whe i = j, R i,i l, m = 0; whe i j, R l, m = m k=l1 A k 1, k. We also edefie Z α ad A 2α as. Z α := sup A 2α := sup Z Zs sup 0 s< 1,s, D s α, sup 0 s< 1,s, D A s, max 1 d s 2α. The ew defiiios ae equivale o 2.2 ad 2.3 sice boh Z ad A ae coiuous pocesses. I is well kow ha a soluio o X ca be cosuced pah-by-pah see [7] ad Secio 6. The ex esul chaaceizes a explici boud fo he eo obaied by appoximaig X usig ˆX. Theoem 2.1. Suppose ha hee exiss a cosa M such ha µ M, µ M ad σ i M fo i = 0, 1, 2, 3, whee σ i deoes he i-h deivaive of σ. If Z α K α <, A 2α K 2α <, ad Γ R < K R, we ca compue G explicily i ems of M, K α, K 2α ad K R, such ha sup ˆX X G 2α β. [0,1] Remak: A ecipe ha explais sep-by-sep how o compue G i ems of algebaic expessios ivolvig M, K α, K 2α ad K R is give i Pocedue A i he appedix o his secio. Usig Theoem 2.1, we ca poceed o sae he mai coibuio of his pape. Theoem 2.2. I he coex of Theoem 2.1, hee is a explici Moe Calo pocedue ha allows us o simulae adom vaiables K α, K 2α, ad K R joily wih {Z : D } fo ay 1. Cosequely, give ay

8 8 J BLANCHET, X. CHEN AND J. DONG deemiisic ε > 0 we ca selec ε such ha G 2α β ε X ε = ˆX so ha ε ad he se 2.5 sup X ε X ε, [0,1] wih pobabiliy oe. Remak: A explici descipio of he algoihm ivolved i he Moe Calo pocedue of Theoem 2.2 is give i Algoihm II a he ed of Secio 5.3, ad he discussio ha follows i. Give { } Z : D ε so ha 2.5 holds, he discussio i he emak ha follows Algoihm II explais how o fuhe simulae {Z : D } fo ay > ε. This efieme is useful i ode o saisfy he impoa popey T4 give i he Ioducio. I deail, oce K α, K 2α, ad K R have bee simulaed he G has also bee simulaed ad evaluaed. Cosequely, give ay sequece ε m < ε m 1 <... < ε 1 we jus eed o obai i such ha G 2α β i ε i. The simulae {Z : D i } ad cosuc ˆX i accodig o 2.4. We le X εi = ˆX i ad, owig o Theoem 2.1, we immediaely obai wih pobabiliy oe, as desied. sup X εi X ε i [0,1] 2.1. O Relaxig Boudedess Assumpios. The cosucio of ˆX i ode o saisfy 2.5 assumes ha µ M, µ 1 M ad σ i M fo i = 0, 1, 2, 3. Alhough hese assumpios ae sog, hee we explai how o elax hem. Theoem 2.2 exeds diecly o he case i which µ ad σ ae Lipschiz coiuous, wih µ diffeeiable ad σ hee imes diffeeiable. Sice µ ad σ ae Lipschiz coiuous we kow ha X has a sog soluio which is o-explosive. We ca always cosuc µ M ad σ M so ha µ i x = µ i M x fo x c M ad i = 0, 1, ad σ i x = σ i M x fo x c M fo i = 0, 1, 2, 3. Also we ca cosuc c M, whee c M as M, ad µ M M, µ 1 M M ad σ i M M fo i = 0, 1, 2, 3. Fo M 1 we coside he SDE 1.1 wih µ M ad σ M as dif ad diffusio coefficies, especively, ad le X M be he coespodig soluio o 1.1. We sa by pickig some M 0 1 such ha ε < c M0 ad le M = M 0. The u Algoihm II o poduce { ˆX M : [0, 1]}, which

9 accodig o Theoem 2.2 saisfies, ε-strong SIMULATION FOR SDES 9 sup ˆX M X M ε. [0,1] Noe ha oly Seps 5 o 8 i Algoihm II deped o he SDE 1.1, hough he evaluaio of G, which depeds o M ad so we wie G M := G. If sup [0,1] ˆX M c M ε, he we mus have ha X = X M fo [0, 1] ad we ae doe. Ohewise, we le M 2M ad u agai oly Seps 5 o 8 of Algoihm II. We epea doublig M ad e-uig Seps 5 o 8 updaig G M uil we obai a soluio fo which sup [0,1] ˆX M c M ε. Eveually his mus occu because lim sup X M X = 0 M [0,1] almos suely ad X is o explosive The Evaluaio of G. We ex summaize he way o calculae G i ems of M, K α, K 2α ad K R. We wie d = max{d, d }. Pocedue A. 1. Fid δ ad C i δ > 0 fo i = 1, 2, 3 ha saisfies he followig elaios: C 1 δ C 3 δδ 2α Mδ 1 α dmk α d 3 M 2 K 2α δ α C 2 δ C 3 δδ α d 3 M 2 K 2α 2 C 3 δ α {MC 1δ dmc 1 δ 2 K α d 2 MC 2 δk α 2 d 3 M 2 C 1 δk 2α } Refe o he poof of Lemma 6.1 fo oe paicula mehod o fid such C i δ s. 2. Se C 1 = 2 δ C 1δ, C 2 = 2 δ C 2δ MC 1 dmc 1 K α ad C 3 = α MC 1 dmc 2 1K α d 2 MC 2 K α 2 d 3 M 2 C 1 K 2α 3. Fid δ ad B i δ fo i = 1, 2, 3 ha saisfies he followig elaios: B 1 δ >B 3 δ δ 2α 2Mδ 1 α 2MK α 4M 2 K 2α δ α B 2 δ >B 3 δ δ α 4M 2 K 2α B 3 δ 4 > α {MB 1δ MB 1 δ 2 K α MB 2 δ K α 2M 2 B 1 δ K 2α }

10 10 J BLANCHET, X. CHEN AND J. DONG 4. Se B = 2 δ B 1 δ 5. Se G 1 = 1 BC 3 6. Fid δ ad C 4 δ such ha Bδ α 2 αβ 2 C 4 δ 21 2 Bδ α 2 αβ 1 B d 3 M 2 K R 2 d 3 M 2 C 1 K R 7. Se C 4 = 1 Bδ α C 4 δ 3 M 2 K R 2 d 3 M 2 C 1 K R /δ 8. Se G 2 = C 4 d 3 M 2 K R 9. Se G = G 1 G 2 Lemma 2.1. Give K α, K 2α, K R ad M, Pocedue A ca be execued. Poof. We pove he lemma by povidig oe paicula mehod o fid such δ ad C i δ s, i = 1, 2, 3. The mehod o fid δ, B i δ s, fo i = 1, 2, 3, follows exacly he same aioale. Se C 1 δ = dm Z α 1/2, C 2 δ = d 3 M 2 A 2α 1/2 ad C 3 δ = α MC 1 δ dmc 1 δ 2 Z α d 2 MC 2 δ Z α d 2 M 2 Z α 2 d 3 M 2 C 1 δ A α. The we ca pick δ small eough, such ha C 3 δδ 2α Mδ 1 α d 3 M 2 A 2α δ α < 1/2 ad C 3 δδ α < 1/2. 3. The mai idea of he algoihmic developme. Based o Theoem 2.2, ou mai ask is o calculae/simulae he uppe boud fo Z α, A 2α ad Γ R especively. I his secio, we will ioduce he mai idea of ou algoihmic developme. The developme ca be decomposed io wo asks. The fis oe is o fid a ifiie sum epeseaio of he objecs of iees. The secod oe is o ucae he ifiie sum up o a fiie bu adom level so ha he eo iduced by he emaiig ems i he summaio is suiably coolled. The secod ask calls fo ovel algoihmic cosucios. Simulaig ifiiely may ems is impossible. We eed o fid a efficie way o exac eough ifomaio o he emaiig ems afe he ucaio, so ha we ca obai a almos sue boud o he coibuio of he ems ha ae o simulaed. We ex cay ou he wo asks oe by oe Ifiie sum epeseaio of Bowia moio ad Lévy aea. We sa by ioducig a wavele syhesis of Bowia moio, {Z : 0 1}, called he Lévy-Ciesielski cosucio of Bowia moio Seele [13]. Fis we eed o defie a sep fucio H o [0, 1] by H = I 0 < 1/2 I 1/2 1.

11 We he defie a family of fucios ε-strong SIMULATION FOR SDES 11 H k = 2/2 H2 1 k 1 fo all 0 ad 1 k 2 1. Se H0 0 = 1. The oe obais he followig ifiie sum epeseaio of Bowia moio. Theoem 3.1 Lévy-Ciesielski Cosucio. If {Wk : 1 k 2 1, 0} is a sequece of idepede sadad omal adom vaiables, he he seies defied by 3.1 Z = W H 0 0 s ds 2 1 =1 k=1 W k 0 H k s ds coveges uifomly o [0, 1] wih pobabiliy oe. Moeove, he pocess {Z : [0, 1]} is a sadad Bowia moio o [0, 1]. Figue 1 demosaes he basic idea of he Lévy-Ciesielski Cosucio usig popeies of he Bowia bidge. Specifically, as Z1 N0, 1, we se Z1 = W0 0. Codiioal o he value of Z0 = 0 ad Z1, Z1/2 NZ1/2, 1/4. Thus we se Z1/2 = Z1/2 1/2W1 1. I geeal, codiioal o he value of Z 1 k ad Z 1 k1, fo k = 0, 1,..., 2 1, Z 2k1 N Z 1 k Z 1 k1 /2, 1 Thus we se Z 2k1 = Z 1 k Z 1 k1 /2 1/2 1 W k1. Eveually we will simulae he seies up o a fiie bu adom level N 1 o be discussed lae. By level we mea he ode of dyadic disceizaio. As we ae simulaig he disceizaio levels sequeially, we ofe efe o ime whe discussig levels. We ex aalyze he Lévy aea, A k, k1, fo 1 i, j d, 1, 0 k 2 1. Usig he algebaic popey A k, k1 =A 1 2k, 1 2k1 Z i 1 2k1 A 1 Zi 1 2k 2k1, 1 2k2 Zj 1 2k2 Zj 1 2k1, we have he followig ifiie sum epeseaio of A k, k1.

12 12 J BLANCHET, X. CHEN AND J. DONG Fig 1. Lévy-Ciesielski Cosucio of Bowia Moio o [0, 1] N0,1 N0, 1 4 N0,1 N0, 1 8 N0, 1 N0, N0, Lemma 3.1. Fo 1, 0 k 2 1, A k, k1 = h=1 2 h 1 l=1 { Z i h 2 h k2l 1 Z i h 2 h k2l 2 Z j h 2 h k2l Z j h 2 h k2l 1 }. The ie summaio ems i he expessio fo A k, k1 moivae he defiiio of he followig family of pocesses L k : k = 0, 1,..., 2 1, 1. L 0 := 0 L k := L k 1 Z i 2k 1 Z i 2k 2 Z j 2k Z j 2k 1 fo k = 1, 2,..., 2 1. Usig his defiiio ad Lemma 3.1 we ca succicly wie A k, k1 as 3.2 A k, k1 = h=1 L h 2 h k 1 L h 2 h k The idea of ecod beakes. To ucae he ifiie sum up o a fiie bu adom level, we use a saegy called ecod beakes. Specifically, we fis defie a sequece of ecod beakes. We he fomulae he fuue ifomaio we eed o kow as a sequece of yes o o quesios. Specifically, he yes o o quesio is fomulaed as will hee be a

13 ε-strong SIMULATION FOR SDES 13 ew ecod beake? ad asweig he yes/o quesio is equivale o simulaig a popely defied Beoulli adom vaiable. The defiiio of he ecod beakes eed o saisfy he followig wo codiios: C1. The followig eve happes wih pobabiliy oe: beyod some adom bu fiie ime, hee will be o moe ecod beakes. C2. By kowig ha hee ae o moe ecod beakes, he coibuio of he ems ha we have o simulaed ye ae well ude cool i.e. bouded by a use defied oleace eo. We ex explai how he above saegy is applied o he Bowia moio ad he Lévy aea especively. We have d idepede Bowia moios ad we will use Wi,k fo i {1,..., d } o deoe he, k coefficie i he expasio 3.1 fo he i-h Bowia moio. Fo Z α, we say a ecod is boke a i,, k, fo 1 i d, 0 ad 1 k 2 1, if Wi,k > 4 1. Le N 1 := max{ 1 : W i,k > 4 1 fo some 1 k 2 1, 1 i d }. I is he las ime he ecod beake happes. The followig Lemma shows ha codiio C1 is saisfied E[N 1 ] < implies P N 1 < = 1. Lemma 3.2. Thee exiss a iege valued adom vaiable N 1, wih E[N 1 ] <, such ha fo all > N 1, 1 k 2 1 ad 1 i d Wi,k 4 1. We ex check codiio C2. Defie V = max 1 k 2 1 Wk. We have he followig auxiliay lemma. Lemma 3.3. Z α 2 2α1 = α V.

14 14 J BLANCHET, X. CHEN AND J. DONG Oce we foud N 1, we have N 1 Z α 2 2α1 2 1/2 α V 2 2α3 =0 =0 =N /2 α 1 N 1 2 2α1 2 1/2 α V 2 2α3 C 2 1/2N 111/2 α 1 2 1/21/2 α. whee C = max N1 1{2 /21/2 α 1}. ad Fo he Lévy aea, we fis oice ha whe i = j, sup = sup A i,i s, sup 0 s< 1,s, D s 2α B Bs 2 s 2 sup 0 s< 1,s, D 2 s 2α Z 2 α 1, 2 R i,i l, m = 0. Whe i j, he ecod beake is defied fo he adom walk L s. Specifically, fo L, we say a ecod is boke a i, j,, k, k, fo 1 i, j d, i j, 1, 0 k < k < 2 1, if L k L k > k k β 2α, whee β 1 α, 2α. Le N 2 := max{ 1 : L k L k > k k β 2α fo some 0 k < k 2 1, 1 i, j d, i j}. I is he las ime he ecod beake happes. The followig lemma shows ha codiio C1 is saisfied. Lemma 3.4. Thee exiss a iege valued adom vaiable N 2, wih E[N 2 ] = o 1 2α 2, such ha fo all > N 2 ad all 0 l < m 2 1 we have L m L l m lβ 2α fo α 1/3, 1/2 ad β 1 α, 2α. We ex check codiio C2. The followig coollay follows diecly fom 3.2 ad he defiiio of R.

15 Coollay 3.1. Fo i j, R l, m = ε-strong SIMULATION FOR SDES 15 h=1 L h 2 h m L h 2 h l. The we have he followig bouds fo A 2α ad Γ R based o he N 2. Lemma 3.5. defie The ad Γ L := max Suppose ha N 2 is chose accodig o Lemma 3.4. We { 1, max 1 d,i j max Γ R max <N 2 0 l<m α β 1 2 2α β Γ L { L m L l m l β 2α 2 1 α 2 A 2α Γ R 1 2 2α Z 2 α 1 2 α. }}. I wha follows, we shall explai how o simulae he adom umbes N 1 ad N 2 joily wih he wavele cosucio usig he ecod beake saegy ioduced i he pevious secio. Specifically, we fis fid all he ecod beakes i sequece ad he simulae he es of he pocess codiioal o he ifomaio obaied by kowig he locaio of all he fiiely may ecod beakes. The challege lies i he fac ha he pobabiliy of success of he Beoulli ials, which coespods o he yes/o quesios defied i ems of he ecod beakes, is o kow o us. We sa wih he pocedue o simulae N 1 i Secio 4, which is buil o a sadwichig idea. The codiioal o he value of N 1, we ioduce he pocedue o simulaio N 2 i Secio 5 based o a accepace-ejecio scheme, whee he poposal disibuio is buil o some expoeial ilig. 4. Toleace-Efoced Simulaio of Bouds o α-hölde Noms. We fis oe ha N 1 is o a soppig ime wih espec o he filaio geeaed by {Wi,k : 0 k 2 1, 1 i d : 1}. Fo he simpliciy of demosaio, we shall focus o he 1-dimesioal case. Fo d > 1, we apply he same pocedue fo each Bowia moio. I wha follows i his subsecio, we shall dop he subscipio i. We call a pai, k a ecod-boke-pai if Wk > 4 1. All pais boh ecod-boke-pais ad o ecod-boke-pais ca be oally odeed lexicogaphically, i.e. usig 2 1 k. The disibuio of subseque

16 16 J BLANCHET, X. CHEN AND J. DONG pais a which ecods ae boke is o difficul o compue because of he idepedece of Wk s. So, usig a sequeial accepace / ejecio pocedue we ca simulae all of he ecod-boke-pais. Codiioal o hese pais, he disibuio of he {Wk : 0 k 2 1 : 1} is saighfowad o descibe. Pecisely, if k, is a ecod-boke-pai, he Wk is codiioed o Wk > 4 1, ad hus is saighfowad o simulae. Similaly, if k, is o a ecod-boke-pai, he Wk is codiioed o Wk 4 1, ad also ca be easily simulaed. The simulaio of he ecod-boke-pais has bee sudied i [5]. The idea is o fid all he ecod beakes sequeially uil hee ae o moe ecod beakes. The challege lies i samplig he Beoulli adom vaiable coespodig o he quesio whehe hee will be o moe ecod beakes i he fuue. We ake samplig he fis beake as a example. The pobabiliy ha hee ae o moe ecod beakes beyod 1 is p1 := 2 1 =1 k=0 P W k 4 1, which ivolves evaluaig he poduc of ifiie may ems ad we do o kow is value i closed fom. Howeve, we ca fid a sequece of uppe boud ad lowe bouds of p1, which ae defied as U h 1 = h P =1 whee = 2 1 k ad W k 4 log 2 1 D h 1 = 1 h 1 42 /2 U h especively. The uppe ad lowe bouds saisfy ha D h 1 < D h1 1 < p1 < U h 1 < U h1 1 ad lim h D h 1 = p1 = lim h U h 1. We also have ha U h 1 U h1 1 is equal o he pobabiliy ha he fis ecod beake happes a posiio h. Thus we ca check whehe he Beoulli ial is a success o failue by updaig he uppe ad lowe bouds sequeially. Moeove, if he Beoulli ial is a failue hee ae moe ecod beakes beyod he cue idex, we also kow he idex of he ex ecod beake. We syhesize algoihm 2W i [5] fo ou puposes ex. Algoihm I: Simulae N 1 joily wih he ecod-boke-pais Oupu: A veco S which gives all he idices l = 2 k such ha, k is a boke-ecod-pai.

17 ε-strong SIMULATION FOR SDES 17 Sep 0: Iiialize R = 0 ad S o be a empy aay. Sep 1: Se U = 1, D = 0. Simulae V Uifom0, 1. Sep 2: While U > V > D, se R R 1 ad U P Wk 4 log 2 R 1 U ad D 1 R 1 42 /2 U. Sep 3: If V U, add R o he ed of S, i.e. S = [S, R], ad eu o Sep 1. Sep 4: If V D, N 1 = log 2 maxs. Sep 5: Oupu S. Ed of Algoihm I Remak: Obseve ha fo evey l = 2 1 k S, we ca geeae Wk codiioal o he eve { Wk > 4 1}; fo ohe l i.e. l / S, geeae Wk give { W k 4 1}. Noe ha a he ed of Algoihm 1 ad afe simulaig Wk fo N 1 oe ca compue N 1 K α = 2 2α1 2 1/2 α V 2 2α3 C 2 1/2N 111/2 α 1 2 1/21/2 α, =0 whee C = max N1 1{2 /21/2 α 1}. 5. Toleace-Efoced Simulaio fo Bouds o 2α-Hölde Noms of Lévy Aeas. The simulaio of N 2, is a lo moe complicaed, compaig o N 1, because hee is fai amou of depedece o he sucue of he L k s as oe vaies. Le us povide a geeal idea of ou simulaio pocedue i ode o se he sage fo he defiiios ad esimaes ha mus be sudied fis. Defie F =σ { Wi,k m : 1 k 2m 1 : m }. ad fo he codiioal expecaio give F we wie E [ ] := E[ F ]. Suppose we have simulaed {W i,k : 0 k 2 1, 1 i d : N} fo some N ad defie τ 1 N = if{ N 1 : L m L l > m l β 2α fo some 0 l < m 2 1 }. Because of Lemma 3.4 we have ha he eve {τ 1 N = } has posiive pobabiliy. I wha follows, we will explai how o simulae a Beoulli

18 18 J BLANCHET, X. CHEN AND J. DONG adom vaiable wih pobabiliy of success P τ 1 N = F N. If such Beoulli is a success, he we have ha N 2 = N ad we would have basically cocluded he difficul pa of he simulaio pocedue he es of he pocess ca be simulaed ude a seies of codiioig eves whose pobabiliy iceases o oe as gows. If he Beoulli is a failue i.e. is value is zeo, he we will fid τ 1 N ad simule all he ifomaio up o τ 1 N. We epea he above Beoulli ial wih updaed pobabiliy of success uil we obai a successful Beoulli ial. Now, pa of he poblem is ha Algoihm I has bee aleady execued, so N N 1, i ohe wods, while he adom vaiables {Wi,k : 1 k 2 1 } ae idepede fo fixed > N, hey ae o loge ideically disibued. Isead, Wi,k is sadad Gaussia codiioal o he eve { Wi,k 4 1}. Neveheless, if is lage eough, all of he eves { W i,k 4 1} will occu wih high pobabiliy. So, we shall fis poceed o explai how o simulae a Beoulli adom vaiable wih pobabiliy of success P τ 1 = F assumig is a deemiisic umbe. The pocedue acually will poduce boh he oucome of he Beoulli ial ad if such oucome is a failue i.e. τ 1 <, also he sample pah {W m i,k : 1 k 2m 1, < m τ 1 }. Ou pocedue is based o accepace / ejecio usig a caefully chose poposal disibuio fo he Wi,k s, based o expoeial ilig of L k s, codiioal o F. To his ed, we will eed o compue he codiioal mome geeaig fucio codiioal o F of L k s ad he family of disibuios iduced ove Wi,k s ad W j,k s ude he expoeially ilig. This will be doe i Secio 5.1. The, we eed some lage deviaio esimaes o boud he likelihood aio of a ceai adomizaio pocedue. These bouds ae developed i Secio 5.2. These ae he mai elemes eeded o simulae N 2 ogehe wih he wavele cosucio. We ioduce he acual adomizaio pocedue ad he deails of he algoihm i Secio Codiioal Mome Geeaig Fucios ad Associaed Expoeial Tilig. I his secio we chaaceize he disibuio of {W m i,k : 1 k 2 m 1 : m 1} ude he expoeial ilig codiioal o F. I ode o educe he legh of some of he equaios ha follow, we wie, fo each {1, 2,..., 2 }, 5.1 Λ i := Z i Z i 1. The we have he followig ecusive elaios fo Λ i s.

19 ε-strong SIMULATION FOR SDES 19 Lemma 5.1. Fo k = 1, 2,..., 2 1 Λ i 2k 1 = 1 2 Λ 1 i 1 k 1/2 1 W i,k. Fom Lemma 5.1, we ca see ha Λ i 2k = 1 2 Λ 1 i 1 k 1/2 1 W i,k, F = σ { Z m k Zm k : 0 k < k 2 m 1, m }. Assume ha k < k, we will ieaively compue he codiioal mome geeaig fucio as { E [exp θ 0 L m k }] 5.2 L m k [ [ { = E E 1...E m 1 [exp θ 0 L m k }] ]] L m k.... Recall ha, fo 1 k 2 1, L k = k =1 We shall sa fom he expecaio of exp codiioal o F m 1. Λ i Coollay 5.1. Fo i j, E m 1 [exp θ 0 Λ m i m 2 1 Λ m j = 1 θ /2 m exp θ 1 Λj m 1 exp η 1 Λ m 1 j m 1 whee 2 1 Λ j 2. θ 0 Λ m i ] m 2 m 1 2 η1 Λ m 1 i θ 1 := θ 0 1 θ m1 1 /4, η1 := θ 2 0 Moeove, defie P m, m E m 1 [I = Wi, m W m i, A, W m j, B A, Wj, m B E m 1 [exp θ 0 Λ m i exp m 2 1 m 2 1 Λ m 1 i m 1 2, Λ m j m 1 m 2 1 θ m1 m /8. θ 0 Λ m i Λ m j m 2 1 m 2 ] Λ m j m 2 ],

20 20 J BLANCHET, X. CHEN AND J. DONG he ude P, ad give F m, m m 1, we have ha Wi, m, Wj, m follows a Gaussia disibuio wih covaiace maix Σ m m1 1 = 1 θ0 2 2 m1 1 θ 0 m1 θ 0 m1 1, ad mea veco µ m m = Σ m m θ 0 1/2 θ 0 1/2 m1 Λm 1 j m 1 m1 Λm 1 i m 1 /2 /2. So, fom Coollay 5.1 we coclude ha E m 1 [exp θ 0 k =k1 Λ m i = 1 θ0 2 2 k k/2 m1 exp 5.3 exp η 1 k =k1 Λ m 1 j m 2 1 θ 1 k Λ m j =k1 ] m 2 Λ m 1 j m 1 2 k η1 =k1 m 1 Λ m 1 i Λ m 1 i m 1 2. m 1 If m 2, we ca coiue akig he coespodig codiioal expecaio give F m 2. Due o he ecusive aue of 5.2 ad he liea ad quadaic ems ha aise i 5.3, i is coveie o coside m 1 =1 2 m 1 =1 θ 1 m 1 η 1 m 1 Λ m 1 j Λ m 1 j m 1 m 1 Λ m 1 i m 1 2 Λ m 1 j m 1 2, whee θ 1 m 1 = θ1 I {k 1,..., k }, η 1 m 1 = η1 I {k 1,..., k }. We also ioduce he followig oaios o simply he peseaio of ou ilig paamees. Due o he diffeece i he ecusive elaio fo Λ i

21 η l m l ε-strong SIMULATION FOR SDES 21 bewee odd ad eve s, we ecusively defie fo l = 2,..., m. 5.5 θ l m l = θ l 1 m l1 2 1 θ l 1 m l1 2, θ l m l = θ l 1 m l1 2 1 θ l 1 m l1 2, η l m l = η l 1 m l1 2 1 η l 1 m l1 2, = η l 1 η l 1, ad se ρ l h l m l m l η l = ηl h l m l m l1 2 1 m l2 θ l m l = 1 2 m l2 η l = 1 2 m l2 η l m l 4 4θ l θ l = θl m l 8 m l m l 4 h l 1 4 θl m l m l m l {θ l η l {θ l m l 2 gl m l m l m l 2 4η l ρ l m l η l η l m l1 2, m l2 m l m l m l 1 ρ l 2 m l }, m l 2 ρl m l m l }. 2, Fially, we decompose 5.4 io wo pas he coss em ad he quadaic em by defiig A m l = θ l m l1 Λ m l1 j m l1 2 1 Λ m l1 i m l1 2 1 θ l 1 2 m l1 Λ m l1 j m l1 2 Λ m l1 i m l1 2, B = η l 1 m l1 m l η l m l1 2 Λ m l1 j Λ m l1 j m l1 2 1 m l1 2 2 Λ m l1 j m l Λ m l1 j m l1 2 2,

22 22 J BLANCHET, X. CHEN AND J. DONG ad C m l = 1 2 m l1 η l The 5.4 ca be wie as 2 m 2 =1 A m 2 m l 1 1 ρ l B m 2, ad he followig esul is key i evaluaig 5.2. m l Coollay 5.2. Fo i j, l = 2, 3,..., m ad = 1, 2,..., 2 m l ] E m l [exp A m l B m l =C m l exp θ l m l Λ i m l Λ j m l exp η 2 2 l Λ i Λj. Moeove, defie P m l1, m l1 E m l [I = m l m l m l 2 1/2. Wi, m l1 A, Wj, m l1 B A, Wj, m l1 B exp A m l ] B m l ], E m l [exp A B W m l1 i, m l m l he ude P, ad give F m l1, m l1 m l, we have ha Wi, m l1, Wj, m l1 follows a Gaussia disibuio wih covaiace maix Σ m l1 1 = 1 ρ l m l1 2 m l 1 2 m l1 η l m l 1 g l m l g l m l 1 2 m l1 η l whee g l m l = m l2 θ l m l 1 2 m l2 η l m l 1 m l 2.

23 ε-strong SIMULATION FOR SDES 23 ad mea veco µ m m l1 = 1/2 m l1 Σ m l1 Λi m l Λ j m l m l1 η l m l 1 η l m l 1 2 Λ j 2 Λ i m l θ l m l m l θ l m l. Usig Coollay 5.2 we coclude ha 2 m l exp A E m l 2 m l = =1 C m l =1 exp m l 2 m l 1 =1 B A m l m l 1 B m l 1. Theefoe, combiig Coollay 5.1 ad epeaedly ieaig he pevious expessio we coclude ha [ ] E expθ 0 {L m k L m k} = 1 θ0 2 2 k k/2 m m 5.6 exp 2 l=2 2 m l =1 C θ m Λ i Λ j =1 m l 2 =1 { η m Λ i 2 Λ j 2} Codiioal Lage Deviaios Esimaes fo L k. We wish o esimae, fo 1 i, j d, i j, k > k ad k, k {0, 1,..., 2 m 1 }, P L m k L m k > k k β 2α m exp θ 0 k k β 2α m {E [expθ 0 {L m k L m k}] E [exp θ 0 {L m k L m k}]}. We boow some iuiio fom he poof of Lemma 3.4 ad selec 5.7 θ 0 m, k γ, k := θ 0 = k k 1/2. 2α m

24 24 J BLANCHET, X. CHEN AND J. DONG We will dop he depedece o m, k, k fo beviy. I addiio, we pick γ 1/4 ad α α, 1/2 so ha exp θ 0 k k β 2α m = exp γ k k β 1/2 2α α 2α 1 m [ ] Ou ex ask is o cool he E expθ 0 {L m k L m k}, which is he pupose of he followig esul, poved i he appedix o his secio. Lemma 5.2. Fo i j, suppose ha θ 0 is chose accodig o 5.7, ad is chose such ha 5.8 max 2 { Λ i, Λ j } α ad fo ε 0 0, 1/2 m 5.9 Λ i Λ j ε 0m l β 2α =l1 fo all 0 l < m 2 wih α α, 1/2, he E [expθ 0 {L m k L m k}] 4 exp ε 0 γk k β 1/2. Remak: I is vey impoa o oe ha due o Lemma { 3.2 we ca always coiue simulaig he Wi,k W s maybe codiioal o } i,k 4 1 i case > N 1 o make sue ha 5.8 holds fo some. Similaly, codiio 5.9 ca be simulaeously efoced wih 5.8 because of Lemma 3.4. Acually, Lemma 3.2 ad Lemma 3.4 idicae ha codiios 5.8 ad 5.9 will occu eveually fo all lage ha some adom heshold. Ou simulaio algoihms will ulimaely deec such heshold, bu Lemma 5.2 does o equie ha we kow ha heshold. As a cosequece of Lemma 5.2, usig Cheoff s boud, we obai he followig poposiio. Poposiio 5.1. Fo i j, if is chose such ha 5.8 ad 5.9 hold, he P L m k L m k > k k β 2α 8 exp 1 2 γ k k β 1/2 2α α 2α 1 m m.

25 ε-strong SIMULATION FOR SDES Joi Toleace-Efoced Simulaio fo α-hölde Noms ad Poof of Theoem Defie C m = { L m k L m k > k k β 2α m fo some 0 k < k < 2 m 1, 1 i, j d, i j}, ad pu τ 1 = if{m 1 : C m occus}. We wie C m fo he compleme of C m, so ha P τ 1 < = m=1 P C m m 1 l=1 C l. To faciliae he explaaio, we ex ioduce a few moe oaios. Le ω :m := {W l i,k : 0 k 2 1, 1 i d, < l m}. I addiio, defie v k, k m :=8 exp 1 2 γ k k β 1/2 2α α 2α 1 m ad I 0 k < k 2 m 1 I m 1 b m := v k, k m 0 k<k 2 m 1 q k, k m := v k, k m b m P,k,k,m ω :m = E [I ω :m exp E [exp We also deoe { ψ m, i, j, k, k := log E [exp θ 0 θ 0 {L m L m ] θ 0 {L m k L m k} ]. k L m k} }] k L m k Obseve ha b m = 8 exp 1 2 γ k k β 1/2 2α α 2α 1 m 0 k<k 2 m 1 2 2m3 exp 1 2 γ 2α α 2α 1 m.

26 26 J BLANCHET, X. CHEN AND J. DONG Thus, b m 0 as. The we ca selec ay pobabiliy mass fucio {gm : m 1}, e.g. gm = e 1 /m 1! fo m 1, by assumig ha is sufficiely lage, gm d 2 b m Now coside he followig pocedue, which we called Pocedue Aux, Aux fo auxilia, which is give fo pedagogical puposes, because as we shall see sholy i is o diecly applicable bu useful o udesad he aue of he mehod ha we shall ulimaely use. Pocedue Aux Ipu: We assume ha we have simulaed {Wi,k : 0 k < 2l : l }}. Oupu: A Beoulli F wih paamee P τ 1 <, ad if F = 1, also ω :τ1 = {Wi,k l : 1 k 2l 1, 1 i d, < l τ 1 } codiioal o he eve τ 1 <. Sep 1: Sample M accodig o g m. Sep 2: Give M = m sample I ad J I J uifomly ove he se {1, 2,..., d }.The, sample K, K fom q k, k m. Sep 3: Give M = m, I = i,j = j,k = k, ad K = k, simulae ω :m fom P,k,k,m o Coollay 5.2. Sep 4: Compue. Noe ha simulaio fom P,k,k,m ca be doe accodig Ξ m, i, j, k, k, ω :m = 1, gm d d 1 1 q k, k m exp θ 0 {L m k L m k} ψ m, i, j, k, k ad N m = I 1 d,i j 1 h<h 2 m 1 L m h L m h > h h β 2α Sep 5: Simulae U uifomly disibued o [0, 1] idepede of eveyhig else ad oupu { } L m F =I{U < I k L m k > k k β 2α m m 1 C l=1 l Ξ m, i, j, k, k, ω :m /N m}. m.

27 If F = 1, also oupu ω :m. Ed of Pocedue Aux We fis oice ha whe ε-strong SIMULATION FOR SDES 27 L m k L m k > k k β 2α m, gm d d 1 1 q k, k m exp θ 0 {L m k L m k} ψ m, i, j, k, k > 1. Thus Ξ m, i, j, k, k, ω :m < 1. Tha is o say he likelihood aio fucio is bouded ad he Beoulli adom vaiable F is well defied. We claim ha he oupu F is disibued as a Beoulli adom vaiable wih paamee P τ 1 <. Moeove, we claim ha if F = 1, he, ω :M is disibued accodig o P ω:τ1 τ 1 <. We fis veify he claim ha he oucome i Sep 5 follows a Beoulli wih paamee P τ 1 <. I ode o see his, le Q deoe he disibuio iduced by Pocedue Aux. Noe ha Q U < I { L M K L M K > K K β 2α Ξ M, I, J, K, K, ω :M /N m { L =E Q M [I K L M K > K K β 2α Ξ M, I, J, K, K, ω :M /N m] = E Q [I m=1 1 d 1 k<k 2 m 1 dp 1 ω :m N m ] = = dp,m,k,k E m=1 1 d 1 k<k 2 m 1 m=1 P C m m 1 l=1 C l =P τ 1 <. { L m { L I m M k L m M } } M 1 l=1 C l M 1 l=1 C l k > k k β 2α m k L m k > k k β 2α N m m } m 1 C l=1 l } m 1 C l=1 l

28 28 J BLANCHET, X. CHEN AND J. DONG Similaly, fo he secod claim, Q ω :M A U < I = = E ω Q :m A, m=1 C M M 1 l=1 dp I,J,K,K,m dp C l Ξ M, I, J, K, K, ω :M ω :m I P ω :m A, τ 1 = m /P τ 1 < m=1 =P ω:τ1 A τ 1 < C m M 1 l=1 The deficiecy of Pocedue Aux is ha i does o ecogize ha > N 1. Le us ow accou fo his fac ad oe ha codiioal o F N1 we have ha Wi,k s ae i.i.d. N0, 1 bu codiioal o { W i,k 4 1} fo all > N 1. Defie H m = { W h i,k 4 h 1 : 0 k 2 h 1, < h m}. I ode o simulae P N1 τ 1 N 1 < we modify sep 3 of Pocedue Aux. Specifically, we have Pocedue B Ipu: We assume ha we have simulaed {Wi,k l : 0 k < 2l : l }. So, he Wi,k m s ae i.i.d. N0, 1 bu codiioal o { W m i,k < 4 m 1} fo all m >. We also assume ha codiios 5.8 ad 5.9 hold i Lemma 5.2; oe he discussio followig Lemma 5.2 which oes ha his ca be assumed a he expese of simulaig addiioal Wi,k m s wih { W m i,k < 4 m 1} if m > N 1. Oupu: A Beoulli F wih paamee P τ 1 <, H, ad if F = 1, also ω :τ1 = {W l i,k : 1 k 2, 1 i d, < l τ 1 } codiioal o τ 1 < ad o H. Sep 1: Sample M accodig o g m. Sep 2: Give M = m sample I ad J I J uifomly ove he se {1, 2,..., d }.The, sample K, K fom q k, k m. Sep 3: Give M = m, I = i,j = j,k = k, ad K = k, simulae ω :m fom P,m,k,k. Noe ha simulaio fom P,m,k,k ca be doe accodig o Coollay 5.2. Check if Hm occus. If yes, coiue o Sep 4; ohewise, go back o Sep 1. C l /P τ 1 <

29 ε-strong SIMULATION FOR SDES 29 Sep 4: Compue Ξ m, i, j, k, k, ω :m = 1, gm d d 1 1 q k, k m exp θ 0 {L m k L m k} ψ m, i, j, k, k ad N m = I 1 d,i j 1 k<k 2 m 1 L m k L m k > k k β 2α Sep 5: Simulae U uifomly disibued o [0, 1] idepede of eveyhig else ad oupu { } I Hm L m k L m k > k k β 2α m M 1 l=1 C l P H m F = I{U < P H Ξ m, i, j, k, k, ω :m /N m} Noice ha P H m /P H = P Hm ad ca be compued i fiie seps. If F = 1, also oupu ω :m. Ed of Pocedue B Le Q deoe he disibuio iduced by Pocedue B. Followig he same aalysis as ha give fo Pocedue Aux, we ca veify ha { } I Hm L m k L m k > k k β 2α m M 1 l=1 C l P H m Q U < P H Ξ m, i, j, k, k, ω :m /N m = P τ 1 < H. Ad if he Beoulli ial is a success, he, ω :M is disibued accodig o P ω:τ1 τ 1 <, H. Fially, if τ 1 =, we may sill eed o simulae ω :m fo ay m 1, bu ow, codiioal o {τ 1 =, H }. Noe ha P ω :m A τ 1 =, H = P ω :m A, τ 1 =, H P τ 1 =, H = E Iω :m A,τ 1 > m, HmP m τ 1 m =, H m P τ 1 =, H. m.

30 30 J BLANCHET, X. CHEN AND J. DONG Thus we ca sample ω :m fom P ad accep he pah wih pobabiliy Iτ 1 > m, H mp m τ 1 m =, H m. This clealy ca be doe sice we ca easily simulae Beoulli s wih pobabiliy P m τm =, H m We summaize he algoihm as follows: = P m τ 1 m = H m P m H m. Algoihm II: Simulae N 1 ad N 2 joily wih Wi,k s fo 1 N 0, whee N 0 is chose such ha sup [0,1] ˆX N 0 X ε Ipu: The paamees equied o u Algoihm I, ad Pocedues A ad B. These ae he ilig paamees θ 0 s. Sep 1: Simulae N 1 joily wih Wi,k m s fo 0 m N 1 usig Algoihm I see he emak ha follows afe Algoihm I. Le = N 1. Sep 2: If ay of he codiios 5.8 ad 5.9 fom Lemma 5.2 ae o saisfied keep simulaig Wi,k m s fo m > uil he fis level m > fo which codiios 5.8 ad 5.9 ae saisfied. Redefie o be such fis level m. Sep 3: Ru Pocedue B ad obai as oupu F ad if F = 1 also obai ω :τ. Sep 4: If τ < i.e. F = 1 se τ ad go back o Sep 2. Ohewise, go o Sep 4. Sep 5: Calculae G accodig o Pocedue A ad solve fo N 0 such ha G 2α β N 0 < ε. Sep 6: If N 0 > sample ω :N0 fom P ad sample a Beoulli adom vaiable, I wih pobabiliy of success P N0 τn 0 =, H N 0. Sep 7: If I = 0, go back o Sep 6. Sep 8: Oupu ω 0:N0. Ed of Algoihm II We obai {W l i,k : 0 k < 2l, l N 0, 1 i d} fom Algoihm II. We have fom ecusios i Lemma 5.1 how o obai 5.10 {Z i l Z i l 1 : 1 2 l, 1 l N 0, 1 i d} ad he we ca compue { ˆX N 0 : D N0 } usig equaio 2.4. Remak: Obseve ha afe compleio of Algoihm II, oe ca acually coiue he simulaio of icemes i ode o obai a appoximaio

31 ε-strong SIMULATION FOR SDES 31 wih a eo ε < ε. I paicula, his is doe by epeaig Seps 4 o 8. Sa fom Sep 4 wih = N 0. The value of G has bee compued, i does o deped o ε. Howeve, oe eeds o ecompue N 0 := N 0 ε such ha < ε. The we ca impleme Seps 5 o 8 wihou chage. Oe obais a oupu ha, as befoe, ca be asfomed io 5.10 via he ecusios 5.1, yieldig { ˆX N 0ε : D N0 ε } wih a guaaeed eo smalle ha ε i uifom om wih pobabiliy 1. G 2α β N 0 6. Rough Diffeeial Equaios, Eo Aalysis, ad The Poof of Theoem 2.1. The aalysis i his secio follows closely he discussio fom [7] Secio 3 ad Secio 7; see also [8] Chape 10. We made some modificaios o accou fo he dif of he pocess ad also o be able o explicily calculae he cosa G. Le us sa wih he defiiio of a soluio o 1.1 usig he heoy of ough diffeeial equaios. We fis povide a defiiio of he soluio of 1.1 i a pahwise sese, followig [7]. Defiiio 6.1. X is a soluio of 1.1 o [0, 1] if X0 = x0 ad fo almos evey sample pah {Z j : j = 1, 2,..., d} i holds X i X i s µ i Xs s σ XsZ j Z j s d d d j=1 l=1 m=1 d j=1 l σ Xsσ l,m XsA m,j s, = o s fo i = 1, 2,..., d ad 0 s < 1, whee A saisfies 6.1 A, = A, s A s, Z i s Z i Z j Z j s fo 0 < s < 1. The pevious defiiio is moivaed by he followig Taylo-ype developme,

32 32 J BLANCHET, X. CHEN AND J. DONG X i h =X i X i d j=1 h h h µ i Xudu µ i Xudu d j=1 h σ XudZ j u σ X µxu σxzu Z dz j u X i µ i Xh σ XZ j h Z j d d d j=1 l=1 m=1 d j=1 l σ Xσ l,m X h Z m u Z m dz j u. The pevious Taylo developme suggess defiig A s, := s Z iu Z i sdz j u. Depedig o how oe iepes As,, e.g. via Iô o Saoovich iegals, oe obais a soluio X which is iepeed i he coespodig coex. I ode o obai he Iô iepeaio of he soluio o equaio 1.1 via defiiio 6.1 we shall iepe he iegals i he sese of Iô. I addiio, as we shall explai, some echical codiios i addiio o he sadad Lipschiz coiuiy ypically equied o obai a sog soluio mus be imposed i ode o efoce he exisece of a uique soluio o 6.1. Thee ae wo souces of eos whe usig ˆX i equaio 2.4 o appoximae X. Oe is he disceizaio o he dyadic gid, bu assumig ha A k, k1 is kow; his ype of aalysis is he oe ha is mos commo i he lieaue o ough pahs see [7]. The secod souce of eo aises due o he fac ha A k, k1 is o kow fo i j. Thus we divide he poof of Theoem 2.1 io wo seps wo poposiios, each dealig wih oe souce of eo. Simila o ˆX, we defie {X : D } by he followig ecusio:

33 give X 0 = X0, ε-strong SIMULATION FOR SDES 33 Xi k1 =X i k µ ix k σ X k Z j k1 Z j k 6.2 d d d j=1 l=1 m=1 d j=1 l σ X k σ l,mx k A m,j k, k1, ad fo [0, 1], we le X = X, whee i his coex = max{s D : s }. Poposiio 6.1. Ude he codiios of Theoem 2.1, we ca compue a cosa G 1 explicily i ems of M, Z α ad A 2α, such ha fo lage eough X X G 1 3α 1. The poof of Poposiio 6.1 will be give afe ioducig some defiiios ad key auxiliay esuls. We deoe Ii, := Xi Xi µ i X σ X Z j Z j ad Ji, := Ii, d d d j=1 l=1 m=1 d j=1 l σ X σ l,m X A m,j,. The followig lemmas ioduce he mai echical esuls fo he poof of Poposiio 6.1. Lemma 6.1. Ude he codiios of Theoem 2.1, hee exis cosas C 1, C 2 ad C 3 ha deped oly o M, Z α ad A 2α, such ha fo ay lage eough ad, D, X X C 1 α, I, C 2 2α, ad J, C 3 3α.

34 34 J BLANCHET, X. CHEN AND J. DONG Poof. Fo s,, s, D, we have he followig impoa ecusios: ad J i, I i, =I i, s I i s, µ i X s µ i X s σ X s σ X Z j Z j s d j=1 =Ji, s Ji s, µ i X s µ i X s d d [σ X s σ X l σ X Xl s X l j=1 l=1 d d l σ X Il, s l σ X µ i X s ]Z j Z j s l=1 6.3 d d d j=1 l=1 m=1 l=1 [ l σ X sσ l,m X s l σ X σ l,m X ] A m,j s, We ex divide he poof io wo pas. We fis pove ha hee exiss a small eough cosa δ > 0 ad hee lage eough cosas C 1 δ, C 2 δ ad C 3 δ, all idepede of, such ha fo < δ, X X C 1 δ α, I, C 2 δ 2α ad J, C 3 δ 3α. We pove i by iducio. Fis we have J, = 0 ad J, = 0. Suppose he esul hold fo all pais of 0, 0 D wih 0 0 <. We he pick s D as he lages poi bewee ad such ha s /2. The we also have s > /2 ad s < /2. Fo simpliciy of oaio, we deoe d = max{d, d }. As Xi Xi s =Ji s, µ i X s s σ X sz j Z j s d d d j=1 l=1 m=1 d j=1 l σ X sσ l,m X sa m,j s,,

35 ε-strong SIMULATION FOR SDES 35 we have X i X i s C 3 δ s 3α M s dm Z α s α d 3 M 2 A 2α s 2α C 3 δδ 2α Mδ 1 α dm Z α d 3 M 2 A 2α δ α s α C 1 δ s α fo C 1 δ C 3 δδ 2α Mδ 1 α dm Z α d 3 M 2 A 2α δ α. Ad as Ii s, = Ji s, we have d d d j=1 l=1 m=1 l σ X sσ l,m X sa m,j s,, I i s, C 3 δ s 3α d 3 M 2 A 2α s 2α C 3 δδ α d 3 M 2 A 2α s 2α C 2 δ s 2α fo C 2 δ C 3 δδ α d 3 M 2 A 2α. We ow aalyze he ecusio 6.3 em by em. Fis, σ X s σ X ad µ i X s µ i X MC 1 δ s α, d l σ X Xl s X l MC 1δ 2 s 2α, l=1 d l σ X Il, s dmc 2 δ s 2α, l=1 d l σ X µ i X s dm 2 s l=1 l σ X sσ l,m X s l σ X σ l,m X 2M 2 C 1 δ s α. The J i, J i, s J i s, MC 1 δ dmc 1 δ 2 Z α d 2 MC 2 δ Z α d 2 M 2 Z α 2 d 3 M 2 C 1 δ A α 3α

36 36 J BLANCHET, X. CHEN AND J. DONG Likewise, we have The fo J i s, J i s, s J i s, MC 1 δ dmc 1 δ 2 Z α d 2 MC 2 δ Z α d 2 M 2 Z α 2 d 3 M 2 C 1 δ A α s 3α = J i s, MC 1 δ dmc 1 δ 2 Z α d 2 MC 2 δ Z α d 2 M 2 Z α 2 d 3 M 2 C 1 δ A α s 3α. J i, J i, s J i s, 2MC 1 δ dmc 1 δ 2 Z α d 2 MC 2 δ Z α d 2 M 2 Z α 2 d 3 M 2 C 1 δ A α s 3α {2 1 3α C 3 δ 2MC 1 δ dmc 1 δ 2 Z α d 2 MC 2 δ Z α d 2 M 2 Z α 2 d 3 M 2 C 1 δ A α } s 3α C 3 δ s 3α, α C 3 δ 2MC 1 δ dmc 1 δ 2 Z α d 2 MC 2 δ Z α d 2 M 2 Z α 2 d 3 M 2 C 1 δ A α. Theefoe, if we delibeaely choose δ, C 1 δ, C 2 δ ad C 3 δ such ha C 1 δ C 3 δδ 2α Mδ 1 α dm Z α d 3 M 2 A 2α δ α C 2 δ C 3 δδ α d 3 M 2 A 2α 2 C 3 δ α MC 1δ dmc 1 δ 2 Z α d 2 MC 2 δ Z α 6.4 d 2 M 2 Z α 2 d 3 M 2 C 1 δ A α The we have fo < δ, X X C 1 δ α, I, C 2 δ 2α, J, C 3 δ 3α,

37 ε-strong SIMULATION FOR SDES 37 The exisece of δ, C 1 δ, C 2 δ ad C 3 δ, saisfyig he sysem of iequaliies 6.4, follows fom Lemma 2.1. We ow exed he aalysis o he case whe > δ. Fo lage eough < δ/2, if > δ, we ca always fid pois s i D ad = s 0 < s 1 < < s k = such ha max 1 i k s i s i 1 < δ ad mi 1 i k s i s i 1 δ/2. The X i X i k Xi s l Xi s l 1 kc 1 δ α 2 δ C 1δ α l=1 Le C 1 = 2 δ C 1δ ad we ca wie X X C 1 α. Nex, I i, k { Ii s l 1, s l µ i X s l µ i X s 0 s l s l 1 l=1 σ ii X s l σ ij X s 0 Z j s l1 Z j s l } d j=1 k[c 2 δ 2α MC 1 1α dmc 1 Z α 2α ] 2 δ C 2δ MC 1 dmc 1 Z α 2α By seig C 2 = 2 δ C 2δ MC 1 dmc 1 Z α, we have I, C 2 2α. Now followig he same iducio aalysis o Ji s, as we did i he case s < δ, we have J i, 2 2 3α C 3 3α If we choose C 3 = 2MC 1 dmc 2 1 Z α d 2 MC 2 Z α 2 d 3 M 2 C 1 A α 3α α MC 1 dmc 2 1 Z α d 2 MC 2 Z α 2 d 3 M 2 C 1 A α, he J, C 3 s 3α. Lemma 6.2. Le x0 ad x0 R d be wo diffee vecos. We deoe X ad X fo D as he -h dyadic appoximaio defied by 6.2 wih iiial value x0 ad x0 especively. Ude he codiios of

38 38 J BLANCHET, X. CHEN AND J. DONG Theoem 2.1, hee exiss a cosa B, idepede of, such ha fo D, X X X 0 X 0 B α X 0 X 0. Moeove, Poof. Le X X 1 B X 0 X 0. Y i,h = X i X i X h 0 X h 0 We defie 0/0 = 0. The followig he ecusio 6.2, we have Y i k1 =Yi k µ ix k µ i X k X 0 X 0 d j=1 6.5 d σ X k σ X k X 0 X 0 Z j k1 Z j k d d j=1 l=1 m=1 l σ X k σ l,mx k lσ X k σ l.m X k X 0 X 0 A m,j k, k1 The 6.2 ad 6.5 ogehe defie a ecusio o geeae X, X ad Y. Followig Lemma 6.1, hee exiss a cosa B ha depeds oly o M, Z α ad A 2α, such ha Thus, Y Y 0 B α. X X X 0 X 0 B α X 0 X 0, ad X X 1 B X 0 X 0. We ae ow eady o pove Poposiio 6.1.

39 ε-strong SIMULATION FOR SDES 39 Poof of Poposiio 6.1. Fom Lemma 6.1 we have X X C 1 α. By Azela-Ascoli Theoem, hee exis a subsequece of {X } ha coveges uifomly o some coiuous fucio X o [0, 1]. Moeove we have X X C 1 α ad X i X i µ i X σ XZ j Z j d d d j=1 l=1 m=1 d j=1 l σ Xσ l.m XA m,j, < C 2 3α Theefoe, he limi X is a soluio o he SDE. Le X,s ; Xs := X s X 0 = Xs. Specifically, we have X,0 ; X0 = X wih X 0 = X0, ad X, ; X = X. The we ca wie X m X m = m X, k m; X k X, k 1 m; X k 1 k=1 By Lemma 6.2, X, k m ; X k X, k 1 m ; X k 1 1B X k X, k 1 k ; X k 1. We also have X i k X, k 1 i k ; X k 1 d = X i k X i k 1 µ ix k 1 k k 1 σ X k 1 Z j k Z j k 1 d d d j=1 l=1 m=1 C 3 k k 1 3α Thus, X m X m j=1 l σ X k 1 σ l,mx k 1 A m,j k 1, k m X, k m; X k X, k 1 m; X k 1 k=1 m1 BC 3 3α 1 BC 3 3α 1. Nex we u o he aalysis of he eo iduced by appoximaig he Lévy aea.

40 40 J BLANCHET, X. CHEN AND J. DONG Poposiio 6.2. Ude he codiios of Theoem 2.1, we ca compue a cosa G 2 explicily i ems of M, Z α, A 2α ad Γ R, such ha fo lage eough whee β 1 α, 2α. ˆX X G 2 2α β, The poof of Poposiio 6.2 uses a simila echique as he poof of Poposiio 6.1 ad also elies o some auxiliay esuls. Le Ui s, := ˆX i X,s i ; ˆX s d d d l σ ˆX sσ l,m ˆX srm,js,. j=1 l=1 m=1 We fis pove he followig echical esul. Lemma 6.3. Ude he codiios of Theoem 2.1, hee exiss a cosa C 4, ha depeds oly o M, Z α, A 2α ad Γ R, such ha U, C 4 αβ 2α β Poof. Fo 0 < s < 1,, s, D, we have U i, =Ui, s Ui s, [ X,s i d d j=1 l=1 m=1 ; ˆX s X, i d Fom Lemma 6.2, X,s i ; ˆX ˆX i s X, s; ˆX ] l σ ˆX sσ l,m ˆX s l σ ˆX σ l,m ˆX Rm,js, ; ˆX s X, i ; ˆX B s α ˆX s X, s; ˆX i ˆX i s X, i s; ˆX Fom Lemma 6.1, l σ ˆX sσ l,m ˆX s l σ ˆX σ l,m ˆX Rm,js, 2M 2 C 1 s α Γ R s β 2α β 2M 2 C 1 Γ R αβ 2α β

41 ε-strong SIMULATION FOR SDES 41 Theefoe, U, U, s U s, B s α ˆX s X, s; ˆX 2 d 3 M 2 C 1 Γ R αβ 2α β U, s U s, B s α U, s d d d B s α max{ l σ ˆX σ l,m ˆX Rm,j, s } i j=1 l=1 m=1 2 d 3 M 2 C 1 Γ R αβ 2α β 1 B s α U, s U s, 6.6 B d 3 M 2 Γ R 2 d 3 M 2 C 1 Γ R αβ 2α β whee d = max{d, d }. Like he poof of Lemma 6.1, we divide he poof io wo pas. We fis pove ha hee exis a small eough cosa δ > 0 ad a lage eough cosa C 4 δ, boh idepede of, such ha fo < δ, U, C 4 δ αβ 2α β. Ad we pove i by iducio. Fis we have U = 0 k, k ad U = 0. Suppose he boud holds fo all pais 0, k, 0 D wih k1 0 0 <. We pick s D as he lages poi bewee ad such ha s 1/2. The we also have s > 1/2 ad s < 1/2. ad U, 1 B s α U, s U s, B d 3 M 2 Γ R 2 d 3 M 2 C 1 Γ R αβ 2α β U s, 1 B α U s, s U s, B d 3 M 2 Γ R 2 d 3 M 2 C 1 Γ R s αβ 2α β U s, B d 3 M 2 Γ R 2 d 3 M 2 C 1 Γ R αβ 2α β Theefoe, U, 1 Bδ α U, s U s, 2B d 3 M 2 Γ R 2 d 3 M 2 C 1 Γ R αβ 2α β 2 Bδα 2 αβ C 4 δ αβ 2α β 2B d 3 M 2 Γ R 2 d 3 M 2 C 1 Γ R αβ 2α β