Optimal Portfolio under Fractional Stochastic Environment

Size: px
Start display at page:

Download "Optimal Portfolio under Fractional Stochastic Environment"

Transcription

1 Opimal Porfolio under Fracional Sochasic Environmen Ruimeng Hu Join work wih Jean-Pierre Fouque Deparmen of Saisics and Applied Probabiliy Universiy of California, Sana Barbara Mahemaical Finance Colloquium Universiy of Souhern California Jan 29, 2018 Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

2 Porfolio Opimizaion: Meron s Problem An invesor manages her porfolio by invesing on a riskless asse B and one risky asse S (single asse for simpliciy) { db = rb d ds = µs d + σs dw π amoun of wealh invesed in he risky asse a ime he wealh process associaed o he sraegy π X π dx π = π ds + Xπ π db (self-financing) S B = (rx π + π (µ r)) d + π σ dw Objecive: M(, x; λ) := sup π A(x,) E [U(X π T ) X π = x where A(x) conains all admissible π and U(x) is a uiliy funcion on R + Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

3 Porfolio Opimizaion: Meron s Problem An invesor manages her porfolio by invesing on a riskless asse B and one risky asse S (single asse for simpliciy) { db = rb d ds = µs d + σs dw π amoun of wealh invesed in he risky asse a ime he wealh process associaed o he sraegy π X π dx π = π ds + Xπ π db (self-financing) S B = (rx π + π (µ r)) d + π σ dw Objecive: M(, x; λ) := sup π A(x,) E [U(X π T ) X π = x where A(x) conains all admissible π and U(x) is a uiliy funcion on R + Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

4 Sochasic Volailiy In Meron s work, µ and σ are consan, complee marke Empirical sudies reveal ha σ exhibis random variaion Implied volailiy skew or smile Sochasic volailiy model: µ(y ), σ(y ) incomplee marke Rough Fracional Sochasic volailiy: Gaheral, Jaisson and Rosenbaum 14 Jaisson, Rosenbaum 16 Omar, Masaaki and Rosenbaum 16 We sudy he Meron problem under slowly varying / fas mean-revering fracional sochasic environmen: Nonlinear + Non-Markovian HJB PDE no available Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

5 Sochasic Volailiy In Meron s work, µ and σ are consan, complee marke Empirical sudies reveal ha σ exhibis random variaion Implied volailiy skew or smile Sochasic volailiy model: µ(y ), σ(y ) incomplee marke Rough Fracional Sochasic volailiy: Gaheral, Jaisson and Rosenbaum 14 Jaisson, Rosenbaum 16 Omar, Masaaki and Rosenbaum 16 We sudy he Meron problem under slowly varying / fas mean-revering fracional sochasic environmen: Nonlinear + Non-Markovian HJB PDE no available Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

6 Relaed Lieraure Opion Pricing + Markovian modeling: Fouque, Papanicolaou, Sircar and Solna 11 (CUP) Porfolio Opimizaion + Markovian modeling: Fouque, Sircar and Zariphopoulou 13 (MF) Fouque and H. 16 (SICON) Opion Pricing + Non-Markovian modeling: Garnier and Solna 15 (SIFIN), 16 (MF) Porfolio Opimizaion + Non-Markovian modeling: Fouque and H. (slow facor, under revision a MF) Fouque and H. (fas facor, under revision a SIFIN) Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

7 A General Non-Markovian Model Dynamics of he risky asse S { ds = S [µ(y ) d + σ(y ) dw, { (W Y : a general sochasic process, G := σ Y ) } -adaped, 0 u wih d W, W Y = ρ d. Dynamics of he wealh process X (assume r = 0 for simpliciy): dx π = π µ(y ) d + π σ(y ) dw Define he value process V by V := ess sup π A E [U(X π T ) F where U(x) is of power ype U(x) = x1 γ 1 γ, γ > 0. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

8 If Y is Markovian... For example, Y is a diffusion process dy = k(y ) d + h(y ) dw Y, V V (, x, y) characerized by a nonlinear HJB PDE Wih a disorion ransformaion 1 V (, x, y) = x1 γ 1 γ Ψ(, y)q, Ψ solves he linear PDE ( 1 Ψ + 2 h2 (y) yy + k(y) y + ρ 1 γ γ λ(y)h(y) y and has he probabilisic represenaion Ψ(, y) = [e Ẽ 1 γ T 2qγ λ 2 (Y s) ds Y = y. 1 Zariphopoulou 99 ) Ψ + 1 γ 2qγ λ2 (y)ψ = 0, Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

9 If Y is Markovian... For example, Y is a diffusion process dy = k(y ) d + h(y ) dw Y, V V (, x, y) characerized by a nonlinear HJB PDE Wih a disorion ransformaion 1 V (, x, y) = x1 γ 1 γ Ψ(, y)q, Ψ solves he linear PDE ( 1 Ψ + 2 h2 (y) yy + k(y) y + ρ 1 γ γ λ(y)h(y) y and has he probabilisic represenaion Ψ(, y) = [e Ẽ 1 γ T 2qγ λ 2 (Y s) ds Y = y. 1 Zariphopoulou 99 ) Ψ + 1 γ 2qγ λ2 (y)ψ = 0, Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

10 If Y is Markovian... For example, Y is a diffusion process dy = k(y ) d + h(y ) dw Y, V V (, x, y) characerized by a nonlinear HJB PDE Wih a disorion ransformaion 1 V (, x, y) = x1 γ 1 γ Ψ(, y)q, Ψ solves he linear PDE ( 1 Ψ + 2 h2 (y) yy + k(y) y + ρ 1 γ γ λ(y)h(y) y and has he probabilisic represenaion Ψ(, y) = [e Ẽ 1 γ T 2qγ λ 2 (Y s) ds Y = y. 1 Zariphopoulou 99 ) Ψ + 1 γ 2qγ λ2 (y)ψ = 0, Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

11 If Y is Markovian + Slowly Varying... For example, Y is a diffusion process dy = δk(y ) d + δh(y ) dw Y, V V (, x, y) characerized by a nonlinear HJB PDE Wih a disorion ransformaion 1 V (, x, y) = x1 γ 1 γ Ψ(, y)q, Ψ solves he linear PDE ( 1 Ψ + 2 δh2 (y) yy + δk(y) y + δρ 1 γ γ λ(y)h(y) y and has he expansion Ψ(, y) = ψ (0) (, y) + δψ (1) (, y) + δψ (2) (, y) +. 1 Zariphopoulou 99 ) Ψ+ 1 γ 2qγ λ2 (y)ψ = 0 Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

12 In General: Maringale Disorion Transformaion 2 The value process V is given by V = X1 γ [Ẽ (e 1 γ T ) 2qγ λ 2 (Y s) ds q µ(y) G, λ(y) = 1 γ σ(y) where under P, W Y := W Y + 0 a s ds is a sandard BM. The opimal sraegy π is [ π λ(y ) = γσ(y ) + ρqξ X γσ(y ) where ξ is given by he maringale represenaion dm = M ξ d W Y and M is M = [e Ẽ 1 γ T 2qγ 0 λ2 (Y s) ds G 2 Tehranchi 04: differen uiliy funcion, proof and assumpions Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

13 Remarks only works for one facor model assumpions: inegrabiliy condiions of ξ, X π and π γ = 1 case of log uiliy, can be reaed separaely degenerae case λ(y) = λ 0, M is a consan maringale, ξ = 0 V = X1 γ 1 γ e 1 γ 2γ λ2 0 (T ), π = λ 0 γσ(y ) X. uncorrelaed case ρ = 0, he problem is linear since q = 1 V = X1 γ [ 1 γ E e 1 γ T 2γ λ 2 (Y s) ds G, π = λ(y ) γσ(y ) X. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

14 Skech of Proof (Verificaion) V is a supermaringale for any admissible conrol π V is a rue maringale following π π is admissible Define α = π /X, hen wih he drif facor D (α ) dv = V D (α ) d + d Maringale D (α ) := α µ γ 2 α2 σ 2 λ2 2γ + q 1 γ a q(q 1) ξ + 2(1 γ) ξ2 + ρqα σξ. α and D (α ) = 0 wih he righ choice of a and q: ( ) 1 γ γ a = ρ λ(y ), q = γ γ + (1 γ)ρ 2. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

15 Skech of Proof (Verificaion) V is a supermaringale for any admissible conrol π V is a rue maringale following π π is admissible Define α = π /X, hen wih he drif facor D (α ) dv = V D (α ) d + d Maringale D (α ) := α µ γ 2 α2 σ 2 λ2 2γ + q 1 γ a q(q 1) ξ + 2(1 γ) ξ2 + ρqα σξ. α and D (α ) = 0 wih he righ choice of a and q: ( ) 1 γ γ a = ρ λ(y ), q = γ γ + (1 γ)ρ 2. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

16 Skech of Proof (Verificaion) V is a supermaringale for any admissible conrol π V is a rue maringale following π π is admissible Define α = π /X, hen wih he drif facor D (α ) dv = V D (α ) d + d Maringale D (α ) := α µ γ 2 α2 σ 2 λ2 2γ + q 1 γ a q(q 1) ξ + 2(1 γ) ξ2 + ρqα σξ. α and D (α ) = 0 wih he righ choice of a and q: ( ) 1 γ γ a = ρ λ(y ), q = γ γ + (1 γ)ρ 2. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

17 Muliple Asses Modeling Consider he following model of S 1, S 2,..., S n ds i = µ i (Y i )S i d + n j=1 σ ij (Y i )S i dw j, i = 1, 2,... n. Each S i is driven by a sochasic facor Y i, bu all facors Y i are adaped o he same single Brownian moion W Y wih he correlaion srucure: d W i, W j = 0, d W i, W Y = ρ d, i, j = 1, 2,..., n, nρ2 < 1. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

18 Maringale Disorion Transformaion wih Muliple Asses Then, he porfolio value V can be expressed as V = X1 γ [Ẽ (e 1 γ T ) 2qγ µ(y s) Σ(Y s) 1 µ(y s) ds q G, 1 γ he consan q is chosen o be: q = The opimal conrol π is given by [ π Σ(Y ) 1 µ(y ) = γ γ γ + (1 γ)ρ 2 n. + ρqξ σ 1 (Y ) 1 n X. γ Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

19 Fracional Processes A fracional Brownian moion W (H), H (0, 1) a coninuous Gaussian process zero [ mean E W (H) W s (H) = σ2 H 2 ( 2H + s 2H s 2H) H < 1/2: shor-range correlaion; H > 1/2: long-range correlaion A fracional Ornsein Uhlenbeck process solves dz H = az H d + dw (H) saionary soluion Z H = e a( s) dw s (H) = K( s) dw s Z Gaussian process wih zero mean and consan variance K is non-negaive, K L 2 Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

20 Fracional Processes A fracional Brownian moion W (H), H (0, 1) a coninuous Gaussian process zero [ mean E W (H) W s (H) = σ2 H 2 ( 2H + s 2H s 2H) H < 1/2: shor-range correlaion; H > 1/2: long-range correlaion A fracional Ornsein Uhlenbeck process solves dz H = az H d + dw (H) saionary soluion Z H = e a( s) dw s (H) = K( s) dw s Z Gaussian process wih zero mean and consan variance K is non-negaive, K L 2 Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

21 Fracional Processes A fracional Brownian moion W (H), H (0, 1) a coninuous Gaussian process zero [ mean E W (H) W s (H) = σ2 H 2 ( 2H + s 2H s 2H) H < 1/2: shor-range correlaion; H > 1/2: long-range correlaion A fracional Ornsein Uhlenbeck process solves dz H = az H d + dw (H) saionary soluion Z H = e a( s) dw s (H) = K( s) dw s Z Gaussian process wih zero mean and consan variance K is non-negaive, K L 2 Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

22 Meron Problem under Slowly Varying Fracional SV Consider a rescaled saionary fou process Z δ,h { [ ds = S µ(z δ,h ) d + σ(z δ,h ) dw, Z δ,h = Kδ ( s) dw Z s, K δ () = δk(δ), d W, W Z = ρ d. Our sudy gives, for all H (0, 1): The value process V δ := ess sup π A δ E [U(X π T ) F The corresponding opimal sraegy π Firs order approximaions o V δ and π A pracical sraegy o generae his approximaed value process Apply he maringale disorion ransformaion wih Y = Z δ,h V δ = X1 γ [Ẽ (e 1 γ T ) 2qγ λ 2 (Zs δ,h ) ds q G, 1 γ Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

23 Meron Problem under Slowly Varying Fracional SV Consider a rescaled saionary fou process Z δ,h { [ ds = S µ(z δ,h ) d + σ(z δ,h ) dw, Z δ,h = Kδ ( s) dw Z s, K δ () = δk(δ), d W, W Z = ρ d. Our sudy gives, for all H (0, 1): The value process V δ := ess sup π A δ E [U(X π T ) F The corresponding opimal sraegy π Firs order approximaions o V δ and π A pracical sraegy o generae his approximaed value process Apply he maringale disorion ransformaion wih Y = Z δ,h V δ = X1 γ [Ẽ (e 1 γ T ) 2qγ λ 2 (Zs δ,h ) ds q G, 1 γ Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

24 Meron Problem under Slowly Varying Fracional SV Consider a rescaled saionary fou process Z δ,h { [ ds = S µ(z δ,h ) d + σ(z δ,h ) dw, Z δ,h = Kδ ( s) dw Z s, K δ () = δk(δ), d W, W Z = ρ d. Our sudy gives, for all H (0, 1): The value process V δ := ess sup π A δ E [U(X π T ) F The corresponding opimal sraegy π Firs order approximaions o V δ and π A pracical sraegy o generae his approximaed value process Apply he maringale disorion ransformaion wih Y = Z δ,h V δ = X1 γ [Ẽ (e 1 γ T ) 2qγ λ 2 (Zs δ,h ) ds q G, 1 γ Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

25 Approximaion o he Value Process Theorem (Fouque-H. 17) For fixed [0, T ), he value process V δ akes he form V δ = X1 γ 1 γ e 1 γ 2γ λ2 (Z δ,h 0 )(T ) + X1 γ + δ H ρ X1 γ 1 γ e 1 γ 2γ λ2 (Z δ,h 0 )(T ) λ 2 (Z δ,h 0 )λ (Z δ,h 0 ) + O(δ 2H ), where φ δ is he random componen of order δ H γ [ T ( φ δ = E Zs δ,h e 1 γ 2γ λ2 (Z δ,h 0 )(T ) λ(z δ,h 0 )λ (Z δ,h 0 )φ δ ( ) 1 γ 2 (T ) H+ 3 2 ) Z δ,h 0 γ ds G. Γ(H ) Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

26 Approximaion o he Opimal Sraegy Recall ha [ π λ(z δ,h ) = γσ(z δ,h ) + ρqξ γσ(z δ,h X ) and ξ is from he maringale rep. of M = [e Ẽ 1 γ T 2qγ 0 λ2 (Zs δ,h ) ds G. Theorem (Fouque-H., 17) The opimal sraegy π is approximaed by [ π λ(z δ,h ) ρ(1 γ) (T ) H+1/2 = γσ(z δ,h + δh ) γ 2 σ(z δ,h ) Γ(H ) λ(z δ,h 0 )λ (Z δ,h 0 ) + O(δ 2H ) := π (0) + δ H π (1) + O(δ 2H ). X Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

27 Approximaion o he Opimal Sraegy Recall ha [ π λ(z δ,h ) = γσ(z δ,h ) + ρqξ γσ(z δ,h X ) and ξ is from he maringale rep. of M = [e Ẽ 1 γ T 2qγ 0 λ2 (Zs δ,h ) ds G. Theorem (Fouque-H., 17) The opimal sraegy π is approximaed by [ π λ(z δ,h ) ρ(1 γ) (T ) H+1/2 = γσ(z δ,h + δh ) γ 2 σ(z δ,h ) Γ(H ) λ(z δ,h 0 )λ (Z δ,h 0 ) + O(δ 2H ) := π (0) + δ H π (1) + O(δ 2H ). X Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

28 How Good is he Approximaion? Corollary In he case of power uiliy U(x) = x1 γ 1 γ, π(0) = λ(zδ,h X γσ(z δ,h generaes he ) approximaion of V δ up o order δ H (leading order + wo correcion erms of order δ H ), hus asympoically opimal in A δ. ) H = 1 2, Zδ,H becomes he Markovian OU process, boh approximaion coincides wih resuls in [Fouque Sircar Zariphopoulou 13. The corollary recovers [Fouque -H. 16. Skech of proofs: Apply Taylor ( expansion ) o λ(z) a he poin Z δ,h 0, and hen conrol he momens Z δ,h Z δ,h 0 δ H. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

29 How Good is he Approximaion? Corollary In he case of power uiliy U(x) = x1 γ 1 γ, π(0) = λ(zδ,h X γσ(z δ,h generaes he ) approximaion of V δ up o order δ H (leading order + wo correcion erms of order δ H ), hus asympoically opimal in A δ. ) H = 1 2, Zδ,H becomes he Markovian OU process, boh approximaion coincides wih resuls in [Fouque Sircar Zariphopoulou 13. The corollary recovers [Fouque -H. 16. Skech of proofs: Apply Taylor ( expansion ) o λ(z) a he poin Z δ,h 0, and hen conrol he momens Z δ,h Z δ,h 0 δ H. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

30 Meron Problem under Fas-Varying Fracional SV Consider a ɛ-scaled saionary fou process Y ɛ,h Y ɛ,h = ɛ H e a( s) ɛ dw s (H) = K ɛ ( s) dws Y, K ɛ () = 1 K( ɛ ɛ ) ogeher wih he risky asse [ ds = S µ(y ɛ,h ) d + σ(y ɛ,h ) dw, d W, W Y = ρ d. For power uiliies, we obain: The value process V ɛ and he corresponding opimal sraegy π Firs order approximaions o V ɛ and π A sraegy π (0) o generae his approximaed value process Using ergodic propery of Y ɛ,h, bu only valid for H ( 1 2, 1) Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

31 Meron Problem under Fas-Varying Fracional SV Consider a ɛ-scaled saionary fou process Y ɛ,h Y ɛ,h = ɛ H e a( s) ɛ dw s (H) = K ɛ ( s) dws Y, K ɛ () = 1 K( ɛ ɛ ) ogeher wih he risky asse [ ds = S µ(y ɛ,h ) d + σ(y ɛ,h ) dw, d W, W Y = ρ d. For power uiliies, we obain: The value process V ɛ and he corresponding opimal sraegy π Firs order approximaions o V ɛ and π A sraegy π (0) o generae his approximaed value process Using ergodic propery of Y ɛ,h, bu only valid for H ( 1 2, 1) Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

32 Meron Problem under Fas-Varying Fracional SV Consider a ɛ-scaled saionary fou process Y ɛ,h Y ɛ,h = ɛ H e a( s) ɛ dw s (H) = K ɛ ( s) dws Y, K ɛ () = 1 K( ɛ ɛ ) ogeher wih he risky asse [ ds = S µ(y ɛ,h ) d + σ(y ɛ,h ) dw, d W, W Y = ρ d. For power uiliies, we obain: The value process V ɛ and he corresponding opimal sraegy π Firs order approximaions o V ɛ and π A sraegy π (0) o generae his approximaed value process Using ergodic propery of Y ɛ,h, bu only valid for H ( 1 2, 1) Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

33 Approximaion o he Value Process V ɛ Theorem (Fouque-H. 17) For fixed [0, T ), he value process V ɛ akes he form V ɛ = X1 γ 1 γ e 1 γ 2γ λ2 (T ) + X1 γ + ɛ 1 H ρ X1 γ 1 γ e 1 γ 2γ λ2 (T ) λ + o(ɛ 1 H ), where φ ɛ is he random componen of order ɛ 1 H e 1 γ 2γ λ2 (T ) φ ɛ γ ( ) 1 γ 2 λλ (T ) H+ 1 2 γ aγ(h ) [ 1 T ( φ ɛ = E λ 2 (Ys ɛ,h ) λ 2) ds 2 G. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

34 Opimal Porfolio Theorem (Fouque-H., 17) The opimal sraegy π is approximaed by [ π λ(y ɛ,h ) ρ(1 γ) λλ (T ) H 1/2 = γσ(y ɛ,h + ɛ1 H ) γ 2 σ(y ɛ,h ) aγ(h ) Corollary + o(ɛ 1 H ) ɛ,h λ(y ) γσ(y ɛ,h In he case of power uiliy, π (0) = X generaes he ) approximaion of V ɛ up o order ɛ 1 H (leading order + wo correcion erms of order ɛ 1 H ), hus asympoically opimal in A ɛ. X Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

35 Ergodiciy of Y ɛ,h For H > 1/2, under appropriae assumpion of λ( ), (λ 2 (Y ɛ,h s (λ(y ɛ,h s (λ(y ɛ,h s are small and of order ɛ 1 H. ) λ 2 ) ds, ) λ) ds, )λ (Ys ɛ,h ) λλ ) ds, Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

36 Comparison wih he Markovian Case The value funcion and he opimal sraegy are derived in [FSZ, 13: [ V ɛ (, X ) = X1 γ 1 γ e 1 γ 2γ λ2 (T ) 1 ( ) 1 γ 2 λθ ɛρ (T ) + O(ɛ) γ 2 [ π (, X, Y ɛ,h λ(y ɛ,h ) ) = γσ(y ɛ,h ) + ρ(1 γ) θ (Y ɛ,h ) ɛ γ 2 σ(y ɛ,h X + O(ɛ) ) 2 Formally le H 1 2 in our resuls: [ V ɛ = X1 γ 1 γ e 1 γ 2γ λ2 (T ) 1+ ( ) 1 γ 2 λ λλ ɛρ (T ) + o( ɛ) γ a [ π λ(y ɛ,h ) = γσ(y ɛ,h ) + ρ(1 γ) λλ ɛ γ 2 σ(y ɛ,h X + o( ɛ) ) a Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

37 Work in preparaion Rough Fas-varying fsv Surprisingly, Y ɛ,h order correcion: V ɛ Y ɛ,h = = X1 γ 1 γ e 1 γ 2γ λ2 (T ) Muliscale fsv ds = S [ Y ɛ,h = µ(y ɛ,h K ɛ ( s) dw Y s, H < 1 2 is no visible o he leading order nor in he firs, Z δ,h K ɛ ( s) dw Y, [ 1 + ( ) 1 γ 2 ɛρ D(T ) + o( ɛ) γ ) d + σ(y ɛ,h Z δ,h =, Z δ,h ) dw, K δ ( s) dw Z. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

38 Work in preparaion Rough Fas-varying fsv Surprisingly, Y ɛ,h order correcion: V ɛ Y ɛ,h = = X1 γ 1 γ e 1 γ 2γ λ2 (T ) Muliscale fsv ds = S [ Y ɛ,h = µ(y ɛ,h K ɛ ( s) dw Y s, H < 1 2 is no visible o he leading order nor in he firs, Z δ,h K ɛ ( s) dw Y, [ 1 + ( ) 1 γ 2 ɛρ D(T ) + o( ɛ) γ ) d + σ(y ɛ,h Z δ,h =, Z δ,h ) dw, K δ ( s) dw Z. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

39 Meron Problem under General Uiliy Maringale Disorion Transformaion is no available Sar wih a given sraegy π (0) A firs order approximaion o V π(0),δ obained by epsilon-maringale decomposiion 34 Opimaliy of π (0) in a smaller class of conrols of feedback form Denoe by v (0) (, x, z) he value funcion a he Sharpe-raio λ(z), we define π (0) by π (0) (, x, z) = λ(z) v x (0) (, x, z) σ(z) v xx (0) (, x, z) and he associaed value process V π(0),δ V π(0),δ 3 Fouque Papanicolaou Sircar 01 4 Garnier Solna 15 := E [ U(X π(0) T ) F. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

40 Meron Problem under General Uiliy Maringale Disorion Transformaion is no available Sar wih a given sraegy π (0) A firs order approximaion o V π(0),δ obained by epsilon-maringale decomposiion 34 Opimaliy of π (0) in a smaller class of conrols of feedback form Denoe by v (0) (, x, z) he value funcion a he Sharpe-raio λ(z), we define π (0) by π (0) (, x, z) = λ(z) v x (0) (, x, z) σ(z) v xx (0) (, x, z) and he associaed value process V π(0),δ V π(0),δ 3 Fouque Papanicolaou Sircar 01 4 Garnier Solna 15 := E [ U(X π(0) T ) F. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

41 Meron Problem under General Uiliy Maringale Disorion Transformaion is no available Sar wih a given sraegy π (0) A firs order approximaion o V π(0),δ obained by epsilon-maringale decomposiion 34 Opimaliy of π (0) in a smaller class of conrols of feedback form Denoe by v (0) (, x, z) he value funcion a he Sharpe-raio λ(z), we define π (0) by π (0) (, x, z) = λ(z) v x (0) (, x, z) σ(z) v xx (0) (, x, z) and he associaed value process V π(0),δ V π(0),δ 3 Fouque Papanicolaou Sircar 01 4 Garnier Solna 15 := E [ U(X π(0) T ) F. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

42 Meron Problem under General Uiliy Maringale Disorion Transformaion is no available Sar wih a given sraegy π (0) A firs order approximaion o V π(0),δ obained by epsilon-maringale decomposiion 34 Opimaliy of π (0) in a smaller class of conrols of feedback form Denoe by v (0) (, x, z) he value funcion a he Sharpe-raio λ(z), we define π (0) by π (0) (, x, z) = λ(z) v x (0) (, x, z) σ(z) v xx (0) (, x, z) and he associaed value process V π(0),δ V π(0),δ 3 Fouque Papanicolaou Sircar 01 4 Garnier Solna 15 := E [ U(X π(0) T ) F. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

43 Epsilon-Maringale Decomposiion Finding Q π(0),δ Then such ha Q π(0),δ T = V π(0),δ T = U(XT π(0) ), Q π(0),δ = M δ + R δ, where M δ is a maringale and R δ is of order δ 2H. V π(0),δ [ = E Q π(0),δ T F = Q π(0),δ R δ + E = M δ + E [RT δ F [R δ T F, and Q π(0),δ is he approximaion o V π(0),δ wih error of order O(δ 2H ) Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

44 Epsilon-Maringale Decomposiion Finding Q π(0),δ Then such ha Q π(0),δ T = V π(0),δ T = U(XT π(0) ), Q π(0),δ = M δ + R δ, where M δ is a maringale and R δ is of order δ 2H. V π(0),δ [ = E Q π(0),δ T F = Q π(0),δ R δ + E = M δ + E [RT δ F [R δ T F, and Q π(0),δ is he approximaion o V π(0),δ wih error of order O(δ 2H ) Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

45 Epsilon-Maringale Decomposiion Finding Q π(0),δ Then such ha Q π(0),δ T = V π(0),δ T = U(XT π(0) ), Q π(0),δ = M δ + R δ, where M δ is a maringale and R δ is of order δ 2H. V π(0),δ [ = E Q π(0),δ T F = Q π(0),δ R δ + E = M δ + E [RT δ F [R δ T F, and Q π(0),δ is he approximaion o V π(0),δ wih error of order O(δ 2H ) Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

46 Firs order approximaion o V π(0),δ Proposiion For fixed [0, T ), he F -measurable value process V π(0),δ is of he form V π(0),δ = Q π(0),δ where Q π(0),δ (x, z) is given by: (X π(0), Z δ,h 0 ) + O(δ 2H ), Q π(0),δ (x, z) =v (0) (, x, z) + λ(z)λ (z)d 1 v (0) (, x, z)φ δ + δ H ρλ 2 (z)λ (z)d1v 2 (0) (T )H+3/2 (, x, z) Γ(H ). For power uiliy, Q π(0),δ coincides wih he approximaion of V δ For he Markovian case H = 1 2, recovers he resuls in [Fouque-H. 16 Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

47 Firs order approximaion o V π(0),δ Proposiion For fixed [0, T ), he F -measurable value process V π(0),δ is of he form V π(0),δ = Q π(0),δ where Q π(0),δ (x, z) is given by: (X π(0), Z δ,h 0 ) + O(δ 2H ), Q π(0),δ (x, z) =v (0) (, x, z) + λ(z)λ (z)d 1 v (0) (, x, z)φ δ + δ H ρλ 2 (z)λ (z)d1v 2 (0) (T )H+3/2 (, x, z) Γ(H ). For power uiliy, Q π(0),δ coincides wih he approximaion of V δ For he Markovian case H = 1 2, recovers he resuls in [Fouque-H. 16 Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

48 Asympoically Opimaliy of π (0) Theorem (Fouque-H. 17) The rading sraegy π (0) (, x, z) = λ(z) σ(z) v x (0) (,x,z) v xx (0) (,x,z) is asympoically opimal in he following class: { Ã δ [ π 0, π 1, α := π = π 0 + δ α π } 1 : π A δ, α > 0, 0 < δ 1. Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

49 Thank you! Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

50 Sochasic Volailiy wih Fas Facor Y S is modeled by: { ds = µ(y )S d + σ(y )S dw, dy = 1 ɛ b(y ) d + 1 ɛ a(y ) dw Y, wih correlaion dw W Y = ρ d. Theorem (Fouque-H., in prep.) Under appropriae assumpions, for fixed (, x, y) and any family of rading sraegies A 0 (, x, y) [ π 0, π 1, α, he following limi exiss and saisfies Ṽ ɛ (, x, y) V π(0),ɛ (, x, y) l := lim 0. ɛ 0 ɛ Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

51 Theorem (Fouque-H., in prep.) The residual funcion E(, x, y) := V π(0),ɛ (, x) v (0) (, x) ɛv (1) (, x) is of order ɛ, where in his case, v (0) solves v (0) 1 2 λ2 and v (1) = 1 2 (T )ρ 1BD 2 1 v(0) (, x). ( v (0) x v (0) xx ) 2 = 0, Ruimeng Hu (UCSB) Opimal Porfolio under Fracional Environmen Jan. 29, / 27

AMartingaleApproachforFractionalBrownian Motions and Related Path Dependent PDEs

AMartingaleApproachforFractionalBrownian Motions and Related Path Dependent PDEs AMaringaleApproachforFracionalBrownian Moions and Relaed Pah Dependen PDEs Jianfeng ZHANG Universiy of Souhern California Join work wih Frederi VIENS Mahemaical Finance, Probabiliy, and PDE Conference

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

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

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

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

More information

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

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

Risk Aversion Asymptotics for Power Utility Maximization

Risk Aversion Asymptotics for Power Utility Maximization Risk Aversion Asympoics for Power Uiliy Maximizaion Marcel Nuz ETH Zurich AnSAp10 Conference Vienna, 12.07.2010 Marcel Nuz (ETH) Risk Aversion Asympoics 1 / 15 Basic Problem Power uiliy funcion U(x) =

More information

Stochastic Modelling in Finance - Solutions to sheet 8

Stochastic Modelling in Finance - Solutions to sheet 8 Sochasic Modelling in Finance - Soluions o shee 8 8.1 The price of a defaulable asse can be modeled as ds S = µ d + σ dw dn where µ, σ are consans, (W ) is a sandard Brownian moion and (N ) is a one jump

More information

Optimal Investment, Consumption and Retirement Decision with Disutility and Borrowing Constraints

Optimal Investment, Consumption and Retirement Decision with Disutility and Borrowing Constraints Opimal Invesmen, Consumpion and Reiremen Decision wih Disuiliy and Borrowing Consrains Yong Hyun Shin Join Work wih Byung Hwa Lim(KAIST) June 29 July 3, 29 Yong Hyun Shin (KIAS) Workshop on Sochasic Analysis

More information

IMA Preprint Series # 2056

IMA Preprint Series # 2056 MULTISCALE ASYMPTOTICS OF PARTIAL HEDGING By Gerard Awanou IMA Preprin Series # 056 ( July 005 ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY OF MINNESOTA 400 Lind Hall 07 Church Sree S.E.

More information

CH Sean Han QF, NTHU, Taiwan BFS2010. (Joint work with T.-Y. Chen and W.-H. Liu)

CH Sean Han QF, NTHU, Taiwan BFS2010. (Joint work with T.-Y. Chen and W.-H. Liu) CH Sean Han QF, NTHU, Taiwan BFS2010 (Join work wih T.-Y. Chen and W.-H. Liu) Risk Managemen in Pracice: Value a Risk (VaR) / Condiional Value a Risk (CVaR) Volailiy Esimaion: Correced Fourier Transform

More information

Optimal Investment Strategy Insurance Company

Optimal Investment Strategy Insurance Company Opimal Invesmen Sraegy for a Non-Life Insurance Company Łukasz Delong Warsaw School of Economics Insiue of Economerics Division of Probabilisic Mehods Probabiliy space Ω I P F I I I he filraion saisfies

More information

Quadratic and Superquadratic BSDEs and Related PDEs

Quadratic and Superquadratic BSDEs and Related PDEs Quadraic and Superquadraic BSDEs and Relaed PDEs Ying Hu IRMAR, Universié Rennes 1, FRANCE hp://perso.univ-rennes1.fr/ying.hu/ ITN Marie Curie Workshop "Sochasic Conrol and Finance" Roscoff, March 21 Ying

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

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility Saisics 441 (Fall 214) November 19, 21, 214 Prof Michael Kozdron Lecure #31, 32: The Ornsein-Uhlenbeck Process as a Model of Volailiy The Ornsein-Uhlenbeck process is a di usion process ha was inroduced

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

Optimal Investment under Dynamic Risk Constraints and Partial Information

Optimal Investment under Dynamic Risk Constraints and Partial Information Opimal Invesmen under Dynamic Risk Consrains and Parial Informaion Wolfgang Puschögl Johann Radon Insiue for Compuaional and Applied Mahemaics (RICAM) Ausrian Academy of Sciences www.ricam.oeaw.ac.a 2

More information

A general continuous auction system in presence of insiders

A general continuous auction system in presence of insiders A general coninuous aucion sysem in presence of insiders José M. Corcuera (based on join work wih G. DiNunno, G. Farkas and B. Oksendal) Faculy of Mahemaics Universiy of Barcelona BCAM, Basque Cener for

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

Portfolio optimization for a large investor under partial information and price impact

Portfolio optimization for a large investor under partial information and price impact Mah Meh Oper Res DOI 1.17/s186-17-589-x Porfolio opimizaion for a large invesor under parial informaion and price impac Zehra Eksi 1 Hyejin Ku Received: 4 Ocober 16 / Acceped: 1 April 17 Springer-Verlag

More information

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

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

More information

Chapter 14 Wiener Processes and Itô s Lemma. Options, Futures, and Other Derivatives, 9th Edition, Copyright John C. Hull

Chapter 14 Wiener Processes and Itô s Lemma. Options, Futures, and Other Derivatives, 9th Edition, Copyright John C. Hull Chaper 14 Wiener Processes and Iô s Lemma Copyrigh John C. Hull 014 1 Sochasic Processes! Describes he way in which a variable such as a sock price, exchange rae or ineres rae changes hrough ime! Incorporaes

More information

and Applications Alexander Steinicke University of Graz Vienna Seminar in Mathematical Finance and Probability,

and Applications Alexander Steinicke University of Graz Vienna Seminar in Mathematical Finance and Probability, Backward Sochasic Differenial Equaions and Applicaions Alexander Seinicke Universiy of Graz Vienna Seminar in Mahemaical Finance and Probabiliy, 6-20-2017 1 / 31 1 Wha is a BSDE? SDEs - he differenial

More information

Existence and uniqueness of solution for multidimensional BSDE with local conditions on the coefficient

Existence and uniqueness of solution for multidimensional BSDE with local conditions on the coefficient 1/34 Exisence and uniqueness of soluion for mulidimensional BSDE wih local condiions on he coefficien EL HASSAN ESSAKY Cadi Ayyad Universiy Mulidisciplinary Faculy Safi, Morocco ITN Roscof, March 18-23,

More information

Some new results on homothetic forward performance processes

Some new results on homothetic forward performance processes Some new resuls on homoheic forward performance processes WCMF, Sana Barbara Sepember 2014 Thaleia Zariphopoulou The Universiy of Texas a Ausin Represenaion of homoheic forward performance processes via

More information

FINM 6900 Finance Theory

FINM 6900 Finance Theory FINM 6900 Finance Theory Universiy of Queensland Lecure Noe 4 The Lucas Model 1. Inroducion In his lecure we consider a simple endowmen economy in which an unspecified number of raional invesors rade asses

More information

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

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

More information

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

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

What Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix

What Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix Wha Ties Reurn Volailiies o Price Valuaions and Fundamenals? On-Line Appendix Alexander David Haskayne School of Business, Universiy of Calgary Piero Veronesi Universiy of Chicago Booh School of Business,

More information

School and Workshop on Market Microstructure: Design, Efficiency and Statistical Regularities March 2011

School and Workshop on Market Microstructure: Design, Efficiency and Statistical Regularities March 2011 2229-12 School and Workshop on Marke Microsrucure: Design, Efficiency and Saisical Regulariies 21-25 March 2011 Some mahemaical properies of order book models Frederic ABERGEL Ecole Cenrale Paris Grande

More information

OBJECTIVES OF TIME SERIES ANALYSIS

OBJECTIVES OF TIME SERIES ANALYSIS OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging

More information

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

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

More information

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy Answer 4 of he following 5 quesions. 1. Consider a pure-exchange economy wih sochasic endowmens. The sae of he economy in period, 0,1,..., is he hisory of evens s ( s0, s1,..., s ). The iniial sae is given.

More information

MATH 5720: Gradient Methods Hung Phan, UMass Lowell October 4, 2018

MATH 5720: Gradient Methods Hung Phan, UMass Lowell October 4, 2018 MATH 5720: Gradien Mehods Hung Phan, UMass Lowell Ocober 4, 208 Descen Direcion Mehods Consider he problem min { f(x) x R n}. The general descen direcions mehod is x k+ = x k + k d k where x k is he curren

More information

Option pricing and implied volatilities in a 2-hypergeometric stochastic volatility model

Option pricing and implied volatilities in a 2-hypergeometric stochastic volatility model Opion pricing and implied volailiies in a 2-hypergeomeric sochasic volailiy model Nicolas Privaul Qihao She Division of Mahemaical Sciences School of Physical and Mahemaical Sciences Nanyang Technological

More information

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

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

More information

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing

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

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

Stochastic Model for Cancer Cell Growth through Single Forward Mutation Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com

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

1 THE MODEL. Monique PONTIER U.M.R. CNRS C 5583, L.S.P. Université Paul Sabatier TOULOUSE cedex 04 FRANCE

1 THE MODEL. Monique PONTIER U.M.R. CNRS C 5583, L.S.P. Université Paul Sabatier TOULOUSE cedex 04 FRANCE Comparison of insider s opimal sraegies, hree differen ypes of side informaion RIMS symposium, he 7h workshop on Sochasic Numerics, June 27-29, 2005 Caroline HILLAIRET CMAPX Ecole Polyechnique 91 128 Palaiseau

More information

Pricing and hedging in stochastic volatility regime switching models.

Pricing and hedging in stochastic volatility regime switching models. Pricing and hedging in sochasic volailiy regime swiching models. Séphane GOUTTE Cenre Naional de la Recherche Scienifique, Laboraoire de Probabiliés e Modèles Aléaoires, Universiés Paris 7 Didero, CNRS,

More information

Local risk minimizing strategy in a market driven by time-changed Lévy noises. Lotti Meijer Master s Thesis, Autumn 2016

Local risk minimizing strategy in a market driven by time-changed Lévy noises. Lotti Meijer Master s Thesis, Autumn 2016 Local risk minimizing sraegy in a marke driven by ime-changed Lévy noises Loi Meijer Maser s Thesis, Auumn 216 Cover design by Marin Helsø The fron page depics a secion of he roo sysem of he excepional

More information

Portfolio Optimization with Nondominated Priors and Unbounded Parameters

Portfolio Optimization with Nondominated Priors and Unbounded Parameters Porfolio Opimizaion wih Nondominaed Priors and Unbounded Parameers Kerem Uğurlu Saurday 14 h July, 218 Deparmen of Applied Mahemaics, Universiy of Washingon, Seale, WA 98195 e-mail:keremu@uw.edu Absrac

More information

Consumption investment optimization with Epstein-Zin utility

Consumption investment optimization with Epstein-Zin utility Consumpion invesmen opimizaion wih Epsein-Zin uiliy Hao Xing London School of Economics Dublin Ciy Universiy, March 6, 2015 1/25 Wha are recursive uiliies? Given a consumpion sream c, V = W(c,m(V +1 )).

More information

Simulation of BSDEs and. Wiener Chaos Expansions

Simulation of BSDEs and. Wiener Chaos Expansions Simulaion of BSDEs and Wiener Chaos Expansions Philippe Briand Céline Labar LAMA UMR 5127, Universié de Savoie, France hp://www.lama.univ-savoie.fr/ Workshop on BSDEs Rennes, May 22-24, 213 Inroducion

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

Examples of Dynamic Programming Problems

Examples of Dynamic Programming Problems M.I.T. 5.450-Fall 00 Sloan School of Managemen Professor Leonid Kogan Examples of Dynamic Programming Problems Problem A given quaniy X of a single resource is o be allocaed opimally among N producion

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

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

Vehicle Arrival Models : Headway

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

More information

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial

More information

arxiv: v2 [q-fin.mf] 7 Nov 2016

arxiv: v2 [q-fin.mf] 7 Nov 2016 Asympoic Opimal Sraegy for Porfolio Opimizaion in a Slowly Varying Sochasic Environmen Jean-Pierre Fouque Ruimeng Hu November 8, 216 arxiv:163.3538v2 [q-fin.mf] 7 Nov 216 Absrac In his paper, we sudy he

More information

Satisfying Convex Risk Limits by Trading

Satisfying Convex Risk Limits by Trading Saisfying Convex Risk Limis by rading Kasper Larsen Deparmen of Accouning and Finance Deparmen of Mahemaics and Compuer Science Universiy of Souhern Denmark DK-53 Odense M kla@sam.sdu.dk raian Pirvu Deparmen

More information

arxiv: v2 [q-fin.mf] 19 Mar 2017

arxiv: v2 [q-fin.mf] 19 Mar 2017 CORRECTION TO BLACK-SCHOLES FORMULA DUE TO FRACTIONAL STOCHASTIC VOLATILITY JOSSELIN GARNIER AND KNUT SØLNA arxiv:59.75v [q-fin.mf 9 Mar 7 Absrac. Empirical sudies show ha he volailiy may exhibi correlaions

More information

1 Solutions to selected problems

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

More information

Lecture Notes 5: Investment

Lecture Notes 5: Investment Lecure Noes 5: Invesmen Zhiwei Xu (xuzhiwei@sju.edu.cn) Invesmen decisions made by rms are one of he mos imporan behaviors in he economy. As he invesmen deermines how he capials accumulae along he ime,

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

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

Dual Representation as Stochastic Differential Games of Backward Stochastic Differential Equations and Dynamic Evaluations

Dual Representation as Stochastic Differential Games of Backward Stochastic Differential Equations and Dynamic Evaluations arxiv:mah/0602323v1 [mah.pr] 15 Feb 2006 Dual Represenaion as Sochasic Differenial Games of Backward Sochasic Differenial Equaions and Dynamic Evaluaions Shanjian Tang Absrac In his Noe, assuming ha he

More information

Optimal Investment and Consumption Decisions Under the Ho-Lee Interest Rate Model

Optimal Investment and Consumption Decisions Under the Ho-Lee Interest Rate Model Opimal Invesmen and Consumpion Decisions Under he Ho-Lee Ineres Rae Model HAO CHANG Tianjin Polyechnic Universiy Deparmen o Mahemaics Binshui Wes Road 399, 300387 Tianjin CHINA ch8683897@126.com XI-MIN

More information

Simulation of BSDEs and. Wiener Chaos Expansions

Simulation of BSDEs and. Wiener Chaos Expansions Simulaion of BSDEs and Wiener Chaos Expansions Philippe Briand Céline Labar LAMA UMR 5127, Universié de Savoie, France hp://www.lama.univ-savoie.fr/ Sochasic Analysis Seminar Oxford, June 1, 213 Inroducion

More information

Forward-backward systems for expected utility maximization

Forward-backward systems for expected utility maximization Forward-backward sysems for expeced uiliy maximizaion Ulrich ors, Ying u, Peer Imkeller, Anhony Réveillac and Jianing Zhang December 23, 213 Absrac In his paper we deal wih he uiliy maximizaion problem

More information

Supplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence

Supplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence Supplemen for Sochasic Convex Opimizaion: Faser Local Growh Implies Faser Global Convergence Yi Xu Qihang Lin ianbao Yang Proof of heorem heorem Suppose Assumpion holds and F (w) obeys he LGC (6) Given

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

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

Loss of martingality in asset price models with lognormal stochastic volatility

Loss of martingality in asset price models with lognormal stochastic volatility Loss of maringaliy in asse price models wih lognormal sochasic volailiy BJourdain July 7, 4 Absrac In his noe, we prove ha in asse price models wih lognormal sochasic volailiy, when he correlaion coefficien

More information

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand Excel-Based Soluion Mehod For The Opimal Policy Of The Hadley And Whiin s Exac Model Wih Arma Demand Kal Nami School of Business and Economics Winson Salem Sae Universiy Winson Salem, NC 27110 Phone: (336)750-2338

More information

(MS, ) Problem 1

(MS, ) Problem 1 MS, 7.6.4) AKTUAREKSAMEN KONTROL I FINANSIERING OG LIVSFORSIKRING ved Københavns Universie Sommer 24 Skriflig prøve den 4. juni 24 kl..-4.. All wrien aids are allowed. The wo problems of oally 3 quesions

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

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion

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

A note on high-order short-time expansions for ATM option prices under the CGMY model

A note on high-order short-time expansions for ATM option prices under the CGMY model A noe on high-order shor-ime expansions for ATM opion prices under he CGMY model José E. Figueroa-López uoing Gong Chrisian Houdré June, 3 Absrac The shor-ime asympoic behavior of opion prices for a variey

More information

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

Introduction to Probability and Statistics Slides 4 Chapter 4 Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random

More information

Optimization problem under change of regime of interest rate

Optimization problem under change of regime of interest rate Opimizaion problem under change of regime of ineres rae Bogdan Ifimie Buchares Universiy of Economic Sudies, and Simion Soilow Insiue of Romanian Academy Bogdan.Ifimie@csie.ase.ro homas Lim Laboraoire

More information

On the Timing Option in a Futures Contract

On the Timing Option in a Futures Contract On he Timing Opion in a Fuures Conrac Francesca Biagini, Mahemaics Insiue Universiy of Munich Theresiensr. 39 D-80333 Munich, Germany phone: +39-051-2094459 Francesca.Biagini@mahemaik.uni-muenchen.de Tomas

More information

Nonlinear expectations and nonlinear pricing

Nonlinear expectations and nonlinear pricing Nonlinear expecaions and nonlinear pricing Zengjing Chen Kun He Deparmen of Mahemaics Shandong Universiy, Jinan 250100 China Reg Kulperger Deparmen of saisical and acuarial Science The Universiy of Wesern

More information

EMS SCM joint meeting. On stochastic partial differential equations of parabolic type

EMS SCM joint meeting. On stochastic partial differential equations of parabolic type EMS SCM join meeing Barcelona, May 28-30, 2015 On sochasic parial differenial equaions of parabolic ype Isván Gyöngy School of Mahemaics and Maxwell Insiue Edinburgh Universiy 1 I. Filering problem II.

More information

Investment performance measurement under asymptotically linear local risk tolerance

Investment performance measurement under asymptotically linear local risk tolerance Invesmen performance measuremen under asympoically linear local risk olerance T. Zariphopoulou and T. Zhou The Universiy of Texas a Ausin Ocober 8, 27 Absrac Using forward opimaliy crieria, we analyze

More information

Stochastic control under progressive enlargement of filtrations and applications to multiple defaults risk management

Stochastic control under progressive enlargement of filtrations and applications to multiple defaults risk management Sochasic conrol under progressive enlargemen of filraions and applicaions o muliple defauls risk managemen Huyên PHAM Laboraoire de Probabiliés e Modèles Aléaoires CNRS, UMR 7599 Universié Paris 7 e-mail:

More information

Exponential utility indifference valuation in two Brownian settings with stochastic correlation

Exponential utility indifference valuation in two Brownian settings with stochastic correlation Exponenial uiliy indifference valuaion in wo Brownian seings wih sochasic correlaion Chrisoph Frei and Marin Schweizer Deparmen of Mahemaics EH Zurich 8092 Zurich Swizerland {chrisoph.frei,marin.schweizer}@mah.ehz.ch

More information

Finance Research Letters. Maximizing utility of consumption subject to a constraint on the probability of lifetime ruin

Finance Research Letters. Maximizing utility of consumption subject to a constraint on the probability of lifetime ruin Finance Research Leers 5 (28 24 212 Conens liss available a ScienceDirec Finance Research Leers www.elsevier.com/locae/frl Maximizing uiliy of consumpion subjec o a consrain on he probabiliy of lifeime

More information

The Structure of General Mean-Variance Hedging Strategies

The Structure of General Mean-Variance Hedging Strategies The Srucure of General Mean-Variance Hedging Sraegies Jan Kallsen TU München (join work wih Aleš Černý, London) Pisburgh, February 27, 2006 1 Quadraic hedging S H (discouned) asse price process (discouned)

More information

PROPERTIES OF MAXIMUM LIKELIHOOD ESTIMATES IN DIFFUSION AND FRACTIONAL-BROWNIAN MODELS

PROPERTIES OF MAXIMUM LIKELIHOOD ESTIMATES IN DIFFUSION AND FRACTIONAL-BROWNIAN MODELS eor Imov r. a Maem. Sais. heor. Probabiliy and Mah. Sais. Vip. 68, 3 S 94-9(4)6-3 Aricle elecronically published on May 4, 4 PROPERIES OF MAXIMUM LIKELIHOOD ESIMAES IN DIFFUSION AND FRACIONAL-BROWNIAN

More information

Stochastic Modelling of Electricity and Related Markets: Chapter 3

Stochastic Modelling of Electricity and Related Markets: Chapter 3 Sochasic Modelling of Elecriciy and Relaed Markes: Chaper 3 Fred Espen Benh, Jūraė Šalyė Benh and Seen Koekebakker Presener: Tony Ware Universiy of Calgary Ocober 14, 2009 The Schwarz model S() = S(0)

More information

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

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

More information

Weyl sequences: Asymptotic distributions of the partition lengths

Weyl sequences: Asymptotic distributions of the partition lengths ACTA ARITHMETICA LXXXVIII.4 (999 Weyl sequences: Asympoic disribuions of he pariion lenghs by Anaoly Zhigljavsky (Cardiff and Iskander Aliev (Warszawa. Inroducion: Saemen of he problem and formulaion of

More information

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS Name SOLUTIONS Financial Economerics Jeffrey R. Russell Miderm Winer 009 SOLUTIONS You have 80 minues o complee he exam. Use can use a calculaor and noes. Try o fi all your work in he space provided. If

More information

Mean-Variance Hedging for General Claims

Mean-Variance Hedging for General Claims Projekbereich B Discussion Paper No. B 167 Mean-Variance Hedging for General Claims by Marin Schweizer ) Ocober 199 ) Financial suppor by Deusche Forschungsgemeinschaf, Sonderforschungsbereich 33 a he

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

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

arxiv: v1 [math.pr] 21 May 2010

arxiv: v1 [math.pr] 21 May 2010 ON SCHRÖDINGER S EQUATION, 3-DIMENSIONAL BESSEL BRIDGES, AND PASSAGE TIME PROBLEMS arxiv:15.498v1 [mah.pr 21 May 21 GERARDO HERNÁNDEZ-DEL-VALLE Absrac. In his work we relae he densiy of he firs-passage

More information

Modern Dynamic Asset Pricing Models

Modern Dynamic Asset Pricing Models Modern Dynamic Asse Pricing Models Teaching Noes 6. Consumpion, Porfolio Allocaion and Equilibrium wih Consrains 1 Piero Veronesi Universiy of Chicago CEPR, NBER 1 These eaching noes draw heavily on Cuoco

More information

Lecture 2 April 04, 2018

Lecture 2 April 04, 2018 Sas 300C: Theory of Saisics Spring 208 Lecure 2 April 04, 208 Prof. Emmanuel Candes Scribe: Paulo Orensein; edied by Sephen Baes, XY Han Ouline Agenda: Global esing. Needle in a Haysack Problem 2. Threshold

More information

Time discretization of quadratic and superquadratic Markovian BSDEs with unbounded terminal conditions

Time discretization of quadratic and superquadratic Markovian BSDEs with unbounded terminal conditions Time discreizaion of quadraic and superquadraic Markovian BSDEs wih unbounded erminal condiions Adrien Richou Universié Bordeaux 1, INRIA équipe ALEA Oxford framework Le (Ω, F, P) be a probabiliy space,

More information

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University THE MYSTERY OF STOCHASTIC MECHANICS Edward Nelson Deparmen of Mahemaics Princeon Universiy 1 Classical Hamilon-Jacobi heory N paricles of various masses on a Euclidean space. Incorporae he masses in he

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

Unit Root Time Series. Univariate random walk

Unit Root Time Series. Univariate random walk Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he

More information