The Smile in Stochastic Volatility Models

Size: px
Start display at page:

Download "The Smile in Stochastic Volatility Models"

Transcription

1 Global Markes Quaniaive Research Modeling and Managing Financial Risks Conference, Paris January 13h, 2011 Join work wih Lorenzo Bergomi

2 Ouline Moivaion Expansion of he smile a order 2 in vol of vol Firs example: a family of Heson-like models Second example: he Bergomi model wih 2 facors on he variance curve Numerical experimens Shor-erm and long-erm asympoics of he smile Rederiving he link beween skew and skewness of log-reurns Conclusion

3 Moivaion Consider he following general dynamics for a diffusive sochasic volailiy model: X = ln S dx = 1 2 ξ d + p ξ dw 1, X 0 = x (1) dξ u = λ(, u, ξ ) dw, ξ u 0 = y u ξ (ξ u, u): insananeous forward variance curve from onwards. For a given mauriy u, ξ u is a drifless process whose iniial value is read on he marke prices of variance swap conracs: ξ0 u = d 2 `ˆσ du u u, where ˆσ u is he implied variance swap volailiy for mauriy u. λ = (λ 1,..., λ d ): volailiy of forward insananeous variances. W = `W 1,..., W d is a d-dimensional Brownian moion. The firs componen of he Brownian moion, W 1, drives he spo dynamics. No dividend. Zero raes and repos (for he sake of simpliciy)

4 No closed-form formula is available for he price of vanilla opions in Model (1). In a few paricular cases of firs generaion sochasic volailiy models, like he Heson model, some approximaions of he price of vanilla opions have been suggesed in he lieraure. Here, we aim a finding a general approximaion of he smile of implied volailiy which does no depend on a paricular specificaion of he model, i.e., on a paricular choice of λ. We will derive general asympoic expansion of he smile, for small volailiy of volailiy, a second order. We inroduce a scaling facor ω for he volailiy of insananeous forward variance: λ ωλ. Abusing language, we will speak of ω as he vol of vol. X and ξ hen depend on ω: X X ω and ξ ξ ω. Our derivaion relies on he fac ha in Model (1), he volailiy of he asse is an auonomous sochasic process, meaning ha i incorporaes no local volailiy componen, and ha λ does no depend on he asse value.

5 Expansion of he price Expansion of he price of a vanilla opion Consider he vanilla opion delivering g(xt ω ) a ime T. Is price is a funcion P ω of (, X ω, ξ,ω ). We wrie P ω (, x, y): he variable y (y u, u T ) is a curve. P ω solves he PDE ( + L ω ) P ω = 0 wih erminal condiion P ω (T, x, y) = g(x), where L ω = L 0 + ωl 1 + ω 2 L 2 wih L 0 = 1 2 y x y 2 x L 1 = L 2 = 1 2 du µ(, u, y) 2 xy u du du ν(, u, u, y) 2 y u y u µ(, u, y) = p y λ 1(, u, y) = E [dxdξu ξ = y] = d h i dx E dξ ν(, u, u, y) = λ i(, u, y)λ i(, u u dξ u ξ = y, y) = d i=1 h i ds E S dξ u ξ = y d

6 Expansion of he price The perurbaion equaions Assume ha P ω = P 0 + ωp 1 + ω 2 P 2 + ω 3 P = ` + L 0 + ωl 1 + ω 2 L 2 `P0 + ωp 1 + ω 2 P 2 + ω 3 P 3 + = ( + L 0) P 0 + ω (( + L 0) P 1 + L 1P 0) +ω 2 (( + L 0) P 2 + L 1P 1 + L 2P 0) +ω 3 (( + L 0) P 3 + L 1P 2 + L 2P 1) + We need o solve he following equaions: ( + L 0) P 0 = 0, P 0(T, x, y) = g(x) ( + L 0) P 1 + L 1P 0 = 0, P 1(T, x, y) = 0 ( + L 0) P n + L 1P n 1 + L 2P n 2 = 0, P n(t, x, y) = 0, n 2

7 Expansion of he price L 0 is he infiniesimal generaor associaed o X 0, he unperurbed diffusion for which ω = 0. I is he sandard one-dimensional Black-Scholes operaor wih deerminisic volailiy y a ime. Each P n is soluion o he radiional one-dimensional diffusion equaion wih a source erm H n = L 1P n 1 + L 2P n 2: ( + L 0) P n + H n = 0 Feynmann-Kac heorem 0,,x P 0(, x, y) = E ˆg `XT,» P n(, x, y) = E H n(s, Xs 0,,x, y)ds, n 1 where X 0,,x is he unperurbed process where ω = 0, saring a log-spo x a ime : dxs 0,,x = 1 2 ys ds + y s dws 1, X 0,,x = x

8 Expansion of he price The price a order 0 P 0 is jus he Black-Scholes price wih ime-dependen volailiy y :» P 0(, x, y) = E g x + ys dws 1 1 «y s ds = P BS x, 2 where «y s ds» P BS(x, v) = E g x + vg 12 «v, G N (0, 1) (2) v = R T y s ds is he oal variance of X 0 inegraed from o T. P 0(, x, y) depends on he curve y (y s, s T ) only hrough v. P BS is soluion o he PDE vp BS = 1 2 ` 2 x x PBS, P BS(x, 0) = g(x) (3) Links he vega and gamma of a vanilla opion in he unperurbed sae.

9 Expansion of he price The price a order 0 An imporan observaion: Because L 0 incorporaes no local volailiy, L 0 and x commue so ( + L 0) xp p 0 = x p ( + L 0) P 0 = 0. xp p BS X 0, R T y s ds p xp 0(, X 0, y) is a maringale for all ineger p. Equaion (3) hen shows ha for all inegers m, n, v m x n P BS X 0, R T y s ds is a maringale. This is crucial in he compuaions of P 1 and P 2.

10 Expansion of he price The price a order 1 Le us define he inegraed spo-variance covariance funcion C Xξ (y): h i Z ds T E s C Xξ S s dξs u ξ s = y (y) = ds du µ(s, u, y) = ds du s s ds We hen have» P 1(, x, y) = E L 1P 0(s, Xs 0,,x, y)ds» = E» = E ds = ds ds s s = C Xξ (y) 2 xvp BS x, s du µ(s, u, y) y u xp BS Xs 0,,x, du µ(s, u, y) xvp 2 BS Xs 0,,x, s» du µ(s, u, y)e xvp 2 BS Xs 0,,x, «y r dr s s «y r dr «y r dr ««y r dr

11 Expansion of he price The price a order 2 A similar resul holds for he second order correcion: C ξξ (y) = C µ (y) = C ξξ P 2 = P L 2P P L 1P 1 2 P L 2P 0 2 (, x, y) = 1 2 Cξξ (y) 2 vp BS x, P L 1P 1 2 = P L 1P 1 2,0 + P L 1P 1 2,1 P L 1P 1 2,0 (, x, y) = 1 2 CXξ (y) 2 2 x 2 vp BS x, P L 1P 1 2,0 (, x, y) = C µ (y) 2 x vp BS x, ds ds du s s s du ν(s, u, u, y) = du µ(s, u, y) y u C Xξ s (y) ds s du «y r dr «y r dr (y): inegraed variance-variance covariance funcion s du «y r dr h E dξ u s dξ u ds i ξ s = y

12 Expansion of he implied volailiy Expansion of he implied volailiy I AT M T = We wrie C Xξ = C Xξ 0 (y), C ξξ = C ξξ 0 (y) and Cµ = C µ 0 (y). In he general diffusive sochasic volailiy model (1), a second order in he vol of vol ω, he implied volailiy for mauriy T and S T = κ T = srike K is quadraic in L = ln I ω (T, K) = I AT M T Coefficiens are r v T + CXξ 4 vt ω + 1 C Xξ 2v 3/2 T ω v 7/2 T K S 0 : ««K + S T ln + κ T ln 2 + O(ω S 0 KS0 3 ) (4) 32v 5/2 T 12C Xξ2 C ξξ v (v + 4) + 4C µ v (v 4) ω 2 4C µ v 3C Xξ2 ω 2 8v 5/2 T 4C µ v + C ξξ v 6C Xξ2 ω 2 v = R T 0 yr dr = inegraed variance.

13 Expansion of he implied volailiy Commens ATM implied volailiy: r IT AT M v = T + CXξ 4 vt ω v 5/2 12C Xξ2 C ξξ v (v + 4) + 4C µ v (v 4) ω 2 T ATM implied volailiy = variance swap volailiy + spread. A firs order, spread = CXξ 4 vt ω. Typically, on he equiy marke, C Xξ < 0: he ATM implied volailiy lies below he variance swap volailiy. When spo reurns and forward variances are uncorrelaed, C Xξ = C µ = 0 so ha r IT AT M v = T Cξξ (v + 4) 32v 3/2 T ω2 Because C ξξ 0, he ATM implied volailiy lies again below he variance swap volailiy. The higher he volailiy of variances, he smaller he ATM implied volailiy.

14 Expansion of he implied volailiy Commens (coninued) ATM skew: S T = CXξ 2v 3/2 ω + 1 `4C T 8v 5/2 µ v 3C Xξ2 ω 2 T The ATM skew S T is of order ω. I has he sign of C Xξ. S T vanishes when spo reurns and forward variances are uncorrelaed, even a second order. The skew is produced by he spo-variance correlaion. If µ 0, S T 0 a firs and a second order. A firs order in ω, he ATM skew has he same sign as he difference beween he ATM implied volailiy and he variance swap volailiy. ATM convexiy: κ T = `4C µ v + C ξξ v 6C Xξ2 ω 2 1 8v 7/2 T The convexiy κ T is of order ω 2. I is an increasing funcion of C ξξ and of C µ, and a decreasing funcion of he absolue value of he spo-variance covariance CXξ. If spo and variances are uncorrelaed, κ T = Cξξ 8v 5/2 T ω2 0 The convexiy is posiive, and produced by variance-variance correlaion.

15 Expansion of he implied volailiy Anoher derivaion which says a he level of operaors Recall ha he price P ω of he vanilla opion is soluion o ( + L ω ) P ω = 0 wih L ω = L 0, + ωl 1, + ω 2 L 2,, and erminal condiion P ω (T, x, y) = g(x). The price can be expressed in erms of he semigroup (U ω s, 0 s T ) aached o he family of differenial operaors L ω : P ω (, ) = U ω T g. The semigroup is defined by U ω s = lim n (1 δlω 0 ) (1 δl ω 1 ) `1 δl ω n 1, s δ =, i = s+iδ n I saisfies Ur ω = UrsU ω s ω for 0 r s T, hence he noaion R : exp s Lω τ dτ :, where :: denoes ime ordering. We can direcly expand U ω s in powers of ω. This is he usual ime-dependen perurbaion echnique in quanum mechanics. U 0 s is called he free propagaor.

16 Expansion of he implied volailiy Consider he general siuaion where a differenial operaor L is perurbed by anoher operaor H : L ε = L + εh From he definiion of he semigroup, Us ε = U (0) s + εu (1) s + ε 2 U (2) s + wih U (1) s = U (2) s = Z s Z s dτ U 0 sτ H τ U 0 τ dτ 1 Z P ω = P 0 + ωp 1 + ω 2 P 2 +, wih P 1 = P 2 = dτ U 0 τ L 1,τ U 0 τt g τ 1 dτ 2 U 0 sτ 1 H τ1 U 0 τ 1 τ 2 H 2 U 0 τ 2 dτ Uτ 0 L 2,τ UτT 0 g + dτ 1 dτ 2 Uτ 0 1 L 1,τ1 Uτ 0 1 τ 2 L 1,τ2 Uτ 0 2 T g τ 1 We recover he expressions of P 1 and P 2. Quicker.

17 Firs example: a Heson-like model dx ω = 1 2 V ω d + p V ω dw 1, X0 ω = x (5), V0 ω = v dv ω = k (V ω v ) d + ω (V ω ) ϕ ρdw 1 + p 1 ρ 2 dw 2 The insananeous forward variance reads ξ u,ω and has dynamics = E [Vu ω V ω ] = v + (V ω dξ u,ω = ωe k(u ),ω `ξ ϕ ρdw 1 v ) e k(u ) + p 1 ρ 2 dw 2 The iniial erm-srucure of insananeous forward variances is y u ξ u 0 = v + (v v ) e ku Like in all classic firs generaion sochasic volailiy models, his erm-srucure is deermined by he model parameers, and he curren value of he insananeous volailiy.

18 The volailiy λ(, u, y) of insananeous forward variances depends on he insananeous forward variance curve y = (y s, s T ) only hrough he insananeous spo variance y : As a consequence, C Xξ = ρ k 2X C ξξ = C µ = i=1 0 0 ϕ ds (y s ) ϕ+ 1 2 λ 1(, u, y) = ρ `y ϕ e k(u ) λ 2(, u, y) = p 1 ρ 2 `y ϕ e k(u ) k(t 1 e s) «2 ds du λ i(s, u, y) = 1 s k 2 «ρ 2 k 0 ds (y s ) ϕ+ 1 2 s 0 ds (y s ) 2ϕ 1 e k(t s) 2 du (y u ) ϕ 1 2 e k(u s) k(t 1 e u) This coincides wih Equaions (3.7) o (3.10) in Lewis [5], where J (1) = C Xξ, J (3) = 1 2 Cξξ, and J (4) = C µ

19 Second example: he Bergomi model dx ω = 1 q 2 ξ,ω d + ξ,ω dw S dξ u,ω = ωξ u,ω α θ (1 θ) e k X (u ) dw X + θe k Y (u ) dw Y = ωλ(, u, ξ,ω ) dw d W S, W X = ρ SXd, d W S, W Y = ρ SY d, d W X, W Y = ρ XY d. The normalizing facor α θ = `(1 θ) 2 + 2ρ XY θ (1 θ) + θ 2 1/2 is such ha he very-shor erm variance ξ,ω has log-normal volailiy ω. We pick k X > k Y, θ is a parameer which mixes he shor-erm facor W X and he long-erm facor W Y.

20 Afer a Cholesky ransform, his can be resaed using independen Brownian moions W 1, W 2 and W 3 as follows: W S = W 1 W X = ρ SXW 1 + q 1 ρ 2 SX W 2 W Y = ρ SY W 1 + χ XY q1 ρ 2 SY W 2 + q (1 χ 2 XY ) (1 ρ2 SY )W 3 where χ XY = ρ XY ρ SX ρ SY 1 ρ 2 SX 1 ρ 2 SY ρ SX, ρ SY and ρ XY define a correlaion marix χ XY [ 1, 1]. The volailiy of variance λ = (λ 1, λ 2, λ 3) reads λ 1(, u, y) = y u α θ (1 θ) ρ SXe k X (u ) + θρ SY e k Y (u ) λ 2(, u, y) = y u α θ (1 θ) q «q1 ρ 2SX e k X (u ) + θχ XY 1 ρ 2 SY e k Y (u ) λ 3(, u, y) = y u α θ θ q (1 χ 2 XY ) (1 ρ2 SY )e k Y (u ) We wrie λ i(, u, y) = y u α θ w ixe k X (u ) + w iy e k Y (u )

21 The covariance funcions read C Xξ = C ξξ = C µ = 0 du Z u 0 d p y λ 1(, u, y) Z u = α θ (1 θ) ρ SX du y u d p y e k X (u ) 0 Z u +α θ θρ SY du y u d p y e k Y (u ) 3X i=1 = α 2 θ 0 3X i=1 0 0 Z u + ds s 0 ds du λ i(s, u, y) s 0 0 «2 ds w ix du y u e k X (u s) + w iy s du y s λ 1(s, u, y u 1 ) s 2 y u dr «y r λ1 (r, u, z) z=yu z u s dλ 1(u,, y ) «2 du y u e k Y (u s)

22 In he case of a fla iniial erm srucure of variance swaps (y 0 ξ), his reads C Xξ = α θ ξ 3/2 (w 1XJ (k X, T ) + w 1Y J (k Y, T )) C ξξ = α 2 θξ 2 (w 0T + w XI(k X, T ) + w Y I(k Y, T ) + w XXI(2k X, T ) +w Y Y I(2k Y, T ) + w XY I(k X + k Y, T )) C µ = «αθξ A1 + A2 wih I(k, u) = J (k, T ) = K(k, T ) = w 0 = Z u = 1 e ku 0 k Z u du = kt 1 + e kt 0 0 k 2 ds (T s)e k(t s) = 0 0 = 3X i=1 w ix k X 3X w iy w Y = 2 i=1 k Y + w iy k Y w ix k X! 2 X 3, wx = 2 w ix + w iy k Y 1 (1 + kt )e kt w ix i=1 k X k X! 3X w, w XX = ix 2 i=1 k X 2 k 2 + w iy k Y! w Y Y = 3X w iy 2 3X w ix w iy i=1 k Y 2, w XY = 2 i=1 k X k Y

23 and A 1 = w2 1X K(k X, T ) w 2 1Y k X k Y K(k Y, T ) + w X J `k X, T + w Y J `k Y, T w 1X w 1Y k X + k Y k X T 1 e k X T k X 2 + e k Y T 1 e k Y T k Y 2 e 2k X T 1 e k Y T + e 2k Y T 1 e k X T 1 A k X k Y A 2 = w X J `k X, T + w Y Y, T + w XX J `2k X, T + w Y Y J `2k Y, T + w XY J `k X + k Y, T wih w X = w 2 1X k X w X = w 2 1X k X + w 1X w 1Y, w Y = w 2 1Y k Y k Y + w 1X w 1Y, w Y = w2 1Y k Y k Y + w 1X w 1Y k X + w 1X w 1Y k X w XX = w2 1X, w Y Y = w2 1Y, w XY = w 1X w 1Y k X k Y k X w 1X w 1Y k Y

24 Firs order Numerical experimens We pick he Bergomi model wih a fla iniial erm srucure of variance swap prices and θ k X k Y ρ SX ρ SY ρ XY χ XY ξ (0.2) 2 ATM implied volaily 20.2% 20.0% 19.8% 19.6% 19.4% 19.2% 19.0% 18.8% 18.6% 18.4% mauriy in years omega = 20%, order 1 omega = 20%, MC omega = 60%, order 1 omega = 60%, MC omega = 200%, order 1 omega = 200%, MC

25 Firs order Firs order ATM skew 4.5% 4.0% 3.5% 3.0% 2.5% 2.0% 1.5% 1.0% omega = 20%, order 1 omega = 20%, MC omega = 60%, order 1 omega = 60%, MC omega = 200%, order 1 omega = 200%, MC 0.5% 0.0% mauriy in years

26 Firs order Firs order smile, omega = 60% 24% 23% 22% 21% 20% 1Y, order 1 1Y, MC 3Y, order 1 3Y, MC 8Y, order 1 8Y, MC 19% 18% 0% 50% 100% 150% 200% 250% srike as percenage of iniial spo

27 Firs order Firs order smile, omega = 200% 40% 35% 30% 25% 20% 1Y, order 1 1Y, MC 3Y, order 1 3Y, MC 8Y, order 1 8Y, MC 15% 10% 0% 50% 100% 150% 200% 250% srike as percenage of iniial spo

28 Firs order Firs order The ATM skew is very sharply esimaed by he firs order expansion, even for large values of he volailiy of variance ω. The ATM volailiy is well capured by he expansion a firs order in ω only for small values of ω (say, up o 60%). True ATM implied volailiies are below heir firs order approximaes he ATM volailiy is a very concave funcion of ω, around ω = 0. In view of he expression for IT AT M, his means ha, for he se of parameers picked, 12C Xξ2 C ξξ v (v + 4) + 4C µ v (v 4) 0 The global shape of he smile is well capured by he firs order expansion: he rue implied volailiy for srike K is indeed approximaely affine in ln(k/s 0), Bu he level of he smile is well capured only for small valus of ω.

29 Second order Second order We firs consider he siuaion when spo reurns and forward variances are uncorrelaed. In his case, he ATM skew vanishes, and so does is expansion a second order in ω. We pick θ k X k Y ρ SX ρ SY ρ XY ξ (0.2) 2 a-he-money implied volailiy 20.50% 20.00% 19.50% 19.00% 18.50% 18.00% omega = 60%, order 2 omega = 60%, MC omega = 120%, order 2 omega = 120%, MC omega = 200%, order 2 omega = 200%, MC omega = 400%, order 2 omega = 400%, MC 17.50% mauriy in years

30 Second order Second order smile, omega = 120% 22.0% 21.5% 21.0% 20.5% 20.0% 1Y, order 2 1Y, MC 3Y, order 2 3Y, MC 8Y, order 2 8Y, MC 19.5% 0% 50% 100% 150% 200% 250% 300% srike as percenage of iniial spo

31 Second order Second order smile, omega = 200% 26% 25% 24% 23% 22% 21% 20% 1Y, order 2 1Y, MC 3Y, order 2 3Y, MC 8Y, order 2 8Y, MC 19% 18% 0% 50% 100% 150% 200% 250% 300% srike as percenage of iniial spo

32 Second order Second order smile, omega = 400% 29% 27% 25% 23% 21% 1Y, order 2 1Y, MC 3Y, order 2 3Y, MC 8Y, order 2 8Y, MC 19% 17% 0% 50% 100% 150% 200% 250% 300% srike as percenage of iniial spo

33 Second order Second order The ATM implied volailiy is very sharply esimaed by he second order expansion, even up o ω = 400% and o long mauriies. For T = 15 years, he esimae is less han 15 bps above he rue ATM volailiy. Looking a he whole smile: he second order expansion of he implied volailiy is excellen around he money, bu becomes oo large for srikes far from he money. No surprising, because no arbirage requires ha for very small and very large srikes, he implied volailiy squared I(T, K) 2 grows a mos linearly wih ln(k/s 0) (see Lee [4]), whereas he second order esimae for I(T, K) 2 grows like ln 4 (K/S 0), see (4). This means ha he remainder O(ω 3 ) = R(ω, T, K) is large for large K, for finie ω. Neverheless, even for ω = 400%, a mauriy of 8 years and a deep ou-he-money srike of 250%, he error is only 1.5 poin of volailiy.

34 Second order Second order We now check numerically he accuracy of he second order expansion of he smile in he general case of correlaed spo reurns and variances. ATM implied volailiy, omega = 120% 20.00% 19.90% 19.80% 19.70% 19.60% 19.50% order 1 order 2 MC 19.40% 19.30% 19.20% mauriy in years

35 Second order Second order 20.2% ATM implied volailiy, omega = 200% 20.0% 19.8% 19.6% 19.4% 19.2% 19.0% order 1 order 2 MC 18.8% 18.6% 18.4% mauriy in years

36 Second order Second order smile, omega = 120% 24% 23% 22% 21% 20% 19% 18% 17% 1Y, order 2 1Y, MC 3Y, order 2 3Y, MC 8Y, order 2 8Y, MC 16% 15% 0% 50% 100% 150% 200% 250% srike as percenage of iniial spo

37 Second order Second order smile, omega = 200% 26% 24% 22% 20% 18% 16% 14% 1Y, order 2 1Y, MC 3Y, order 2 3Y, MC 8Y, order 2 8Y, MC 12% 10% 0% 50% 100% 150% 200% 250% srike as percenage of iniial spo

38 Shor-erm asympoics of implied volailiy Shor-erm asympoics of implied volailiy Assume dξ = d + ω(ξ) ϕ db Le ρ SV be he correlaion beween S and insananeous variance V = ξ Heson: ϕ = 1, ρsv = ρ; 2 Bergomi: ϕ = 1, ρ SV = α θ ((1 θ)ρ SX + θρ SY ) Then for shor mauriies qξ `IAT 00 + ρsv M 2ϕ T T ω + T 4ϕ 3 IT AT M `(5 8ϕ) ρ 2 SV 4 ω 2 + O `ω `IAT S T ρsv M 2ϕ 2 T ω + T 4ϕ 3 IT AT M ρ 2 2 SV ϕ 5 «ω 2 + O `ω 3 12 I AT M T κ T `IAT M 4ϕ 5 T ϕ 7 «ρ 2 SV + 1 «ω 2 + O `ω Shor-erm ATM skew does no depend on shor-erm ATM vol iff ϕ = 1 (observed in equiy markes) Shor-erm ATM convexiy does no depend on shor-erm ATM vol iff ϕ = 5 4. And ( ρsv, κt 0) ϕ 5 4

39 Long-erm asympoics of implied volailiy Long-erm asympoics of implied volailiy Assume he erm-srucure of variance swaps volailiies is fla: ξ 0 ξ. Assume ha for large u, µ(, u, y) (u ) α, α > 0. Then a higher order in ω, for long mauriies, S T T α if α < 1 S T T 1 if α > 1 Cf Bergomi, Smile Dynamics 4 [2], for he link wih he skew-sickiness raio. Assume ha for large u and u, ν(, u, u, y) (u ) α (u ) α, α > 0. Also assume ha spos and volailiies are uncorrelaed (µ 0). Then a higher order in ω, for long mauriies, Exponenial decay α > 1. κ T T 2α if α < 1 κ T T 2 if α > 1

40 Remember S T = CXξ 2v 3/2 ω + O `ω 2 T Le us now compue he skewness s T of log-reurns: s T = E ˆX T 3, E [XT 2 XT = XT E [XT ] = ]3/2 We have E ˆX R T 2 T =,ω R T E ˆξ 0 d = 0 ξ 0d and E ˆXT 3 = 3ωC Xξ + O `ω 2 0 q ξ,ω dw 1 A firs order in he vol of vol, he skewness of (he disribuion of) ln (S T /S 0) is hus 3ωC Xξ s T = R 3/2 T 0 ξ 0 d The ATM skew S T simply reads S T = st 6 T + O(ω2 )

41 Conclusion We have considered general second generaion sochasic volailiy models and derived an expansion of he smile of implied volailiy a second order in he volailiy of variance. A his order, he smile is quadraic in L = ln(k/s 0). This expansion shows ha he smile is driven by hree model-dependen quaniies: C Xξ, he inegraed spo-variance covariance funcion, C ξξ, he inegraed variance-variance covariance funcion, C µ, which, like C Xξ, depends only on he insananeous spo-variance covariance, bu in a more complex way. C Xξ drives he ATM implied volailiy and ATM skew, a firs order. When spo reurns are uncorrelaed wih variances, he smile is U-shaped, C Xξ = C µ = 0 and C ξξ drives boh he ATM implied volailiy and ATM convexiy. In he general case where spo reurns are correlaed wih variances, he second order correcion for ATM implied volailiy, he second order correcion for he ATM skew, and he ATM convexiy are all driven by a linear combinaion of `C Xξ 2, C ξξ and C µ.

42 Our derivaion relies on he fac ha in Model (1), he volailiy of he asse is an auonomous sochasic process, meaning ha i incorporaes no local volailiy componen, and ha λ does no depend on he asse value. These hree fundamenal covariance funcions are compued for Heson-like radiional sochasic volailiy models (In paricular we give a new derivaion of A. Lewis resuls [5]). They are also compued for he Bergomi model wih wo facors on he variance curve. In his case, numerical experimens show an excellen agreemen beween he esimae and he rue quaniy a firs order for he ATM skew, and a second order for he ATM implied volailiy, hence for he whole smile, up o ypical values of he volailiy of variance (say, 400% in he equiy marke). We have given shor-erm and long-erm behaviour of he smile. We have also rederived ha he ATM skew is jus he skewness of he disribuion of he log-spo, divided by 6 imes he square roo of mauriy, a firs order in he volailiy of variance.

43 Bergomi L., Smile Dynamics 2, Risk Magazine, pages 67-73, Ocober Bergomi L., Smile Dynamics 4, Risk Magazine, December Backus D., Foresi S., Li K. and Wu L., Accouning for Biases in Black-Scholes, unpublished. Lee R., The momen formula for implied volailiy a exreme srikes, Sanford Universiy and Couran Insiue, Lewis A., Opion valuaion under sochasic volailiy, Finance Press, 2000.

Stochastic Volatility s Orderly Smiles

Stochastic Volatility s Orderly Smiles Quan Research Fields Quaniaive Finance Seminar Fields Insiue for Research in Mahemaical Sciences Torono, February 6h, 2013 Join work wih Lorenzo Bergomi (SocGen) jguyon2@bloomberg.ne lorenzo.bergomi@sgcib.com

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

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

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

13.3 Term structure models

13.3 Term structure models 13.3 Term srucure models 13.3.1 Expecaions hypohesis model - Simples "model" a) shor rae b) expecaions o ge oher prices Resul: y () = 1 h +1 δ = φ( δ)+ε +1 f () = E (y +1) (1) =δ + φ( δ) f (3) = E (y +)

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

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

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

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

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

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

Chapter 2. First Order Scalar Equations

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

More information

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS

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

More information

Chapter 7 Response of First-order RL and RC Circuits

Chapter 7 Response of First-order RL and RC Circuits Chaper 7 Response of Firs-order RL and RC Circuis 7.- The Naural Response of RL and RC Circuis 7.3 The Sep Response of RL and RC Circuis 7.4 A General Soluion for Sep and Naural Responses 7.5 Sequenial

More information

Problem Set on Differential Equations

Problem Set on Differential Equations Problem Se on Differenial Equaions 1. Solve he following differenial equaions (a) x () = e x (), x () = 3/ 4. (b) x () = e x (), x (1) =. (c) xe () = + (1 x ()) e, x () =.. (An asse marke model). Le p()

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

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

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims Problem Se 5 Graduae Macro II, Spring 2017 The Universiy of Nore Dame Professor Sims Insrucions: You may consul wih oher members of he class, bu please make sure o urn in your own work. Where applicable,

More information

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1 Chaper 5 Heerocedasic Models Inroducion o ime series (2008) 1 Chaper 5. Conens. 5.1. The ARCH model. 5.2. The GARCH model. 5.3. The exponenial GARCH model. 5.4. The CHARMA model. 5.5. Random coefficien

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

Lecture 4: Processes with independent increments

Lecture 4: Processes with independent increments Lecure 4: Processes wih independen incremens 1. A Wienner process 1.1 Definiion of a Wienner process 1.2 Reflecion principle 1.3 Exponenial Brownian moion 1.4 Exchange of measure (Girsanov heorem) 1.5

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

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

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

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

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

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

More information

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

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

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

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

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

Lecture 9: Advanced DFT concepts: The Exchange-correlation functional and time-dependent DFT

Lecture 9: Advanced DFT concepts: The Exchange-correlation functional and time-dependent DFT Lecure 9: Advanced DFT conceps: The Exchange-correlaion funcional and ime-dependen DFT Marie Curie Tuorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dep. of Chemisry and Couran Insiue

More information

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that ODEs, Homework #4 Soluions. Check ha y ( = is a soluion of he second-order ODE ( cos sin y + y sin y sin = 0 and hen use his fac o find all soluions of he ODE. When y =, we have y = and also y = 0, so

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

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

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

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

Testing for a Single Factor Model in the Multivariate State Space Framework

Testing for a Single Factor Model in the Multivariate State Space Framework esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics

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

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

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

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

More information

Math Final Exam Solutions

Math Final Exam Solutions Mah 246 - Final Exam Soluions Friday, July h, 204 () Find explici soluions and give he inerval of definiion o he following iniial value problems (a) ( + 2 )y + 2y = e, y(0) = 0 Soluion: In normal form,

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

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

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

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

More information

DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND

DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND Asymmery and Leverage in Condiional Volailiy Models Michael McAleer WORKING PAPER

More information

ADVANCED MATHEMATICS FOR ECONOMICS /2013 Sheet 3: Di erential equations

ADVANCED MATHEMATICS FOR ECONOMICS /2013 Sheet 3: Di erential equations ADVANCED MATHEMATICS FOR ECONOMICS - /3 Shee 3: Di erenial equaions Check ha x() =± p ln(c( + )), where C is a posiive consan, is soluion of he ODE x () = Solve he following di erenial equaions: (a) x

More information

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

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

More information

Class Meeting # 10: Introduction to the Wave Equation

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

More information

Couplage du principe des grandes déviations et de l homogénisation dans le cas des EDP paraboliques: (le cas constant)

Couplage du principe des grandes déviations et de l homogénisation dans le cas des EDP paraboliques: (le cas constant) Couplage du principe des grandes déviaions e de l homogénisaion dans le cas des EDP paraboliques: (le cas consan) Alioune COULIBALY U.F.R Sciences e Technologie Universié Assane SECK de Ziguinchor Probabilié

More information

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity ANNALES POLONICI MATHEMATICI LIV.2 99) L p -L q -Time decay esimae for soluion of he Cauchy problem for hyperbolic parial differenial equaions of linear hermoelasiciy by Jerzy Gawinecki Warszawa) Absrac.

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

Generalized Chebyshev polynomials

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

More information

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)

More information

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

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

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

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

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

Chapter 6. Systems of First Order Linear Differential Equations

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

More information

CHAPTER 2 Signals And Spectra

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

More information

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

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

More information

ψ(t) = V x (0)V x (t)

ψ(t) = V x (0)V x (t) .93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in

More information

Regular Variation and Financial Time Series Models

Regular Variation and Financial Time Series Models Regular Variaion and Financial Time Series Models Richard A. Davis Colorado Sae Universiy www.sa.colosae.edu/~rdavis Thomas Mikosch Universiy of Copenhagen Bojan Basrak Eurandom Ouline Characerisics of

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

Math 333 Problem Set #2 Solution 14 February 2003

Math 333 Problem Set #2 Solution 14 February 2003 Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial

More information

Two Coupled Oscillators / Normal Modes

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

More information

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

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

More information

ENGI 9420 Engineering Analysis Assignment 2 Solutions

ENGI 9420 Engineering Analysis Assignment 2 Solutions ENGI 940 Engineering Analysis Assignmen Soluions 0 Fall [Second order ODEs, Laplace ransforms; Secions.0-.09]. Use Laplace ransforms o solve he iniial value problem [0] dy y, y( 0) 4 d + [This was Quesion

More information

Theory of! Partial Differential Equations-I!

Theory of! Partial Differential Equations-I! hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

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

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

Modeling Economic Time Series with Stochastic Linear Difference Equations

Modeling Economic Time Series with Stochastic Linear Difference Equations A. Thiemer, SLDG.mcd, 6..6 FH-Kiel Universiy of Applied Sciences Prof. Dr. Andreas Thiemer e-mail: andreas.hiemer@fh-kiel.de Modeling Economic Time Series wih Sochasic Linear Difference Equaions Summary:

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

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

INDEX. Transient analysis 1 Initial Conditions 1

INDEX. Transient analysis 1 Initial Conditions 1 INDEX Secion Page Transien analysis 1 Iniial Condiions 1 Please inform me of your opinion of he relaive emphasis of he review maerial by simply making commens on his page and sending i o me a: Frank Mera

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

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance Lecure 5 Time series: ECM Bernardina Algieri Deparmen Economics, Saisics and Finance Conens Time Series Modelling Coinegraion Error Correcion Model Two Seps, Engle-Granger procedure Error Correcion Model

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

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

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

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

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

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3. Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies

More information

GMM - Generalized Method of Moments

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

More information

Undetermined coefficients for local fractional differential equations

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

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

EXERCISES FOR SECTION 1.5

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

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

More information

Bernoulli numbers. Francesco Chiatti, Matteo Pintonello. December 5, 2016

Bernoulli numbers. Francesco Chiatti, Matteo Pintonello. December 5, 2016 UNIVERSITÁ DEGLI STUDI DI PADOVA, DIPARTIMENTO DI MATEMATICA TULLIO LEVI-CIVITA Bernoulli numbers Francesco Chiai, Maeo Pinonello December 5, 206 During las lessons we have proved he Las Ferma Theorem

More information

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015 Explaining Toal Facor Produciviy Ulrich Kohli Universiy of Geneva December 2015 Needed: A Theory of Toal Facor Produciviy Edward C. Presco (1998) 2 1. Inroducion Toal Facor Produciviy (TFP) has become

More information

Transform Techniques. Moment Generating Function

Transform Techniques. Moment Generating Function Transform Techniques A convenien way of finding he momens of a random variable is he momen generaing funcion (MGF). Oher ransform echniques are characerisic funcion, z-ransform, and Laplace ransform. Momen

More information

Chapter 3 Boundary Value Problem

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

More information

References are appeared in the last slide. Last update: (1393/08/19)

References are appeared in the last slide. Last update: (1393/08/19) SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be

More information

Optimizing heat exchangers

Optimizing heat exchangers Opimizing hea echangers Jean-Luc Thiffeaul Deparmen of Mahemaics, Universiy of Wisconsin Madison, 48 Lincoln Dr., Madison, WI 5376, USA wih: Florence Marcoe, Charles R. Doering, William R. Young (Daed:

More information

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff Laplace ransfom: -ranslaion rule 8.03, Haynes Miller and Jeremy Orloff Inroducory example Consider he sysem ẋ + 3x = f(, where f is he inpu and x he response. We know is uni impulse response is 0 for

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

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations! hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

HPCFinance research project 8

HPCFinance research project 8 HPCFinance research projec 8 Financial models, volailiy risk, and Bayesian algorihms Hanxue Yang Tampere Universiy of Technology March 14, 2016 Research projec 8 12/2012 11/2015, Tampere Universiy of Technology,

More information

Chapter 15. Time Series: Descriptive Analyses, Models, and Forecasting

Chapter 15. Time Series: Descriptive Analyses, Models, and Forecasting Chaper 15 Time Series: Descripive Analyses, Models, and Forecasing Descripive Analysis: Index Numbers Index Number a number ha measures he change in a variable over ime relaive o he value of he variable

More information