Brace-Gatarek-Musiela model

Similar documents
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 15 10/30/2013. Ito integral for simple processes

10.7 Temperature-dependent Viscoelastic Materials

5.1 Angles and Their Measure

AP Physics 1 MC Practice Kinematics 1D

Lecture 4 ( ) Some points of vertical motion: Here we assumed t 0 =0 and the y axis to be vertical.

The Components of Vector B. The Components of Vector B. Vector Components. Component Method of Vector Addition. Vector Components

Physics Courseware Physics I Constant Acceleration

Kinematics Review Outline

i-clicker Question lim Physics 123 Lecture 2 1 Dimensional Motion x 1 x 2 v is not constant in time v = v(t) acceleration lim Review:

The Buck Resonant Converter

Lecture 3: Resistive forces, and Energy

Solutions to Assignment 1

GAMS Handout 2. Utah State University. Ethan Yang

and Sun (14) and Due and Singlen (19) apply he maximum likelihd mehd while Singh (15), and Lngsa and Schwarz (12) respecively emply he hreesage leas s

Ramsey model. Rationale. Basic setup. A g A exogenous as in Solow. n L each period.

independenly fllwing square-r prcesses, he inuiive inerpreain f he sae variables is n clear, and smeimes i seems dicul nd admissible parameers fr whic

Section 5.8 Notes Page Exponential Growth and Decay Models; Newton s Law

Physics 111. Exam #1. September 28, 2018

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

LHS Mathematics Department Honors Pre-Calculus Final Exam 2002 Answers

Differentiation Applications 1: Related Rates

Motion Along a Straight Line

SMKA NAIM LILBANAT KOTA BHARU KELANTAN. SEKOLAH BERPRESTASI TINGGI. Kertas soalan ini mengandungi 7 halaman bercetak.

CHAPTER 7 CHRONOPOTENTIOMETRY. In this technique the current flowing in the cell is instantaneously stepped from

Chapter 2. First Order Scalar Equations

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

Impact Switch Study Modeling & Implications

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

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

This section is primarily focused on tools to aid us in finding roots/zeros/ -intercepts of polynomials. Essentially, our focus turns to solving.

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.

Announcements: Warm-up Exercise:

4.6 One Dimensional Kinematics and Integration

Unit-I (Feedback amplifiers) Features of feedback amplifiers. Presentation by: S.Karthie, Lecturer/ECE SSN College of Engineering

Lecture 13: Markov Chain Monte Carlo. Gibbs sampling

Economics 8105 Macroeconomic Theory Recitation 6

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

8. Basic RL and RC Circuits

Visco-elastic Layers

Math 333 Problem Set #2 Solution 14 February 2003

AP Statistics Notes Unit Two: The Normal Distributions

EECE 301 Signals & Systems Prof. Mark Fowler

PRINCE SULTAN UNIVERSITY Department of Mathematical Sciences Final Examination Second Semester (072) STAT 271.

Direct Current Circuits. February 19, 2014 Physics for Scientists & Engineers 2, Chapter 26 1

Coherent PSK. The functional model of passband data transmission system is. Signal transmission encoder. x Signal. decoder.

System of Linear Differential Equations

Lecture 13 RC/RL Circuits, Time Dependent Op Amp Circuits

13.3 Term structure models

An application of nonlinear optimization method to. sensitivity analysis of numerical model *

Thermodynamics Partial Outline of Topics

SPH3U1 Lesson 06 Kinematics

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

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

Differential Equations

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

Successive ApproxiInations and Osgood's Theorenl

Part a: Writing the nodal equations and solving for v o gives the magnitude and phase response: tan ( 0.25 )

Algebra2/Trig: Trig Unit 2 Packet

EXERCISES FOR SECTION 1.5

Math 116 Second Midterm March 21, 2016

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

ELEG 205 Fall Lecture #10. Mark Mirotznik, Ph.D. Professor The University of Delaware Tel: (302)

Department of Economics, University of California, Davis Ecn 200C Micro Theory Professor Giacomo Bonanno. Insurance Markets

2 Models of Forward LIBOR and Swap Raes 2 Preliminaries Le T > be a horizon dae for our model of economy. In his preliminary secion, we focus on LIBOR

Unit Root Time Series. Univariate random walk

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

Homework #7. True False. d. Given a CFG, G, and a string w, it is decidable whether w ε L(G) True False

15. Vector Valued Functions

x x

DC-DC Switch-Mode Converters

(2) Even if such a value of k was possible, the neutrons multiply

2.4 Cuk converter example

Money in OLG Models. 1. Introduction. Econ604. Spring Lutz Hendricks

Week #13 - Integration by Parts & Numerical Integration Section 7.2

Introduction. If there are no physical guides, the motion is said to be unconstrained. Example 2. - Airplane, rocket

Solutions Problem Set 3 Macro II (14.452)

HOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.

1 Consumption and Risky Assets

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du.

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

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

Lead/Lag Compensator Frequency Domain Properties and Design Methods

Laplace Transforms. Examples. Is this equation differential? y 2 2y + 1 = 0, y 2 2y + 1 = 0, (y ) 2 2y + 1 = cos x,

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

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

Revelation of Soft-Switching Operation for Isolated DC to Single-phase AC Converter with Power Decoupling

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

AP Calculus BC Chapter 10 Part 1 AP Exam Problems

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

Week 7: Dynamic Price Setting

Physics 2B Chapter 23 Notes - Faraday s Law & Inductors Spring 2018

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Distributions, spatial statistics and a Bayesian perspective

Announcements. Formulas Review. Exam format

4 Sequences of measurable functions

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Chapter 7 Response of First-order RL and RC Circuits

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.

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

Transcription:

Chaper 34 Brace-Gaarek-Musiela mdel 34. Review f HJM under risk-neural IP where f ( T Frward rae a ime fr brrwing a ime T df ( T ( T ( T d + ( T dw ( ( T The ineres rae is r( f (. The bnd prices saisfy B( T IE " ( ; ( ; ( u du r(u du F( f ( u du db( T r( B( T d ; ( T {z } B( T dw ( vlailiy f T -mauriy bnd. T implemen HJM, yu specify a funcin A simple chice we wuld like use is ( T T ( T f( T where > is he cnsan vlailiy f he frward rae. This is n pssible because i leads # ( T df ( T 2 f ( T f ( u du f ( u du! d + f( T dw ( 335

336 and Heah, Jarrw and Mrn shw ha sluins his equain lde befre T. The prblem wih he abve equain is ha he d erm grws like he square f he frward rae. T see wha prblem his causes, cnsider he similar deerminisic rdinary differenial equain where f ( c>. We have This sluin ldes a c. f ( f 2 ( ; d d f ( ; f ( + f ( f ( f 2 ( du ; f ( ; f ( f ( c ; c c ; ; c c 34.2 Brace-Gaarek-Musiela mdel New variables Curren ime Time mauriy T ; Frward raes Bnd prices r( f ( + r( f ( r( (2. r( f ( + (2.2 T D( B( + (2.3 (u v ; du dv D( ; ; ; + f ( v dv f ( + u du B( + ;r( D( (2.4 T

CHAPTER 34. Brace-Gaarek-Musiela mdel 337 We will nw wrie ( ( T ; raher han ( T. In his nain, he HJM mdel is df ( T ( ( d + ( dw ( (2.5 db( T r(b( T d ; ( B( T dw ( (2.6 where ( ( u du (2.7 ( ( (2.8 We nw derive he differenials f r( and D(, analgus (2.5 and (2.6 We have Als, dr( df ( + {z } + f ( + d T differenial applies nly firs argumen (2.5,(2.2 ( ( d + ( dw ( + r( d (2.8 i hr( + 2 ( ( 2 d + ( dw ( (2.9 dd( db( + {z } + B( + d T differenial applies nly firs argumen (2.6,(2.4 r( B( + d ; ( B( + dw ( ; r( D( d (2. [r( ; r( ] D( d ; ( D( dw ( (2. 34.3 LIBOR Fix > (say, year. $ D( invesed a ime in a ( + -mauriy bnd grws $ a 4 ime +. L( is defined be he crrespnding rae f simple ineres +L( D( ( + L( ( D( nr L( r( u du ;

338 34.4 Frward LIBOR > is sill fixed. A ime, agree inves $ D( + D( D( + D( a ime +, wih payback f $ a ime + +. Can d his a ime by shring bnds mauring a ime + and ging lng ne bnd mauring a ime + +. The value f his prfli a ime is ; D( + D( +D( + D( The frward LIBOR L( is defined be he simple (frward ineres rae fr his invesmen D( + D( Cnnecin wih frward raes ( + L( +L( D( D( + L( R f; g n R + ; ( + nr + ; (4. s ( + r( + r( f ( + r( lim # L( nr + ( + nr + ; ; > fixed r( is he cninuusly cmpunded rae. L( is he simple rae ver a perid f durain. We cann have a lg-nrmal mdel fr r( because sluins lde as we saw in Secin 34.. Fr fixed psiive, wecan have a lg-nrmal mdel fr L(. (4.2 34.5 The dynamics f L( We wan chse (, appearing in (2.5 s ha dl( ( d + L( ( dw (

CHAPTER 34. Brace-Gaarek-Musiela mdel 339 fr sme (. This is he BGM mdel, and is a subclass f HJM mdels, crrespnding paricular chices f (. Recall (2.9 dr( i hr( u+ 2 u ( ( u 2 d + ( u dw ( Therefre, d +! + d (5. + i + hr( u+ 2 u ( ( u 2 du d + ( u du dw ( i hr( + ; r( + 2 ( ( + 2 ; 2 ( ( 2 d +[ ( + ; ( ] dw ( and dl( (4 d 2 nr 4 + ( + (4., (5. d ; 3 5 + ( + + 2 + d [+L( ] (5.2 [r( + ; r( + 2 ( ( + 2 ; 2 ( ( 2 ] d +[ ( + ; ( ] dw ( + 2 [ ( + ; ( ] 2 d [ + L( ] [r( + ; r( ] d + ( + [ ( + ; ( ] d +[ ( + ; ( ] dw (! 2

34 Bu L( 2 nr 4 + ( + ; 3 5 [r( + ; r( ] [ + L( ][r( + ; r( ] Therefre, dl( L( d + [ + L( ][ ( + ; ( ][ ( + d + dw (] Take ( be given by Then ( L( [+L( ][ ( + ; ( ] (5.3 dl( [ L( +( L( ( + ] d + ( L( dw ( (5.4 Ne ha (5.3 is equivalen ( + ( + Plugging his in (5.4 yields " dl( L( +( L( ( + L2 ( 2 ( +L( L( ( +L( (5.3 # d + ( L( dw ( (5.4 34.6 Implemenain f BGM Obain he iniial frward LIBOR curve L( frm marke daa. Chse a frward LIBOR vlailiy funcin (usually nnrandm (

CHAPTER 34. Brace-Gaarek-Musiela mdel 34 Because LIBOR gives n rae infrmain n ime perids smaller han, we mus als chse a parial bnd vlailiy funcin ( < fr mauriies less han frm he curren ime variable. Wih hese funcins, we can fr each 2 [ slve (5.4 bain L( < Plugging he sluin in (5.3, we bain ( fr <2. We hen slve (5.4 bain L( <2 and we cninue recursively. Remark 34. BGM is a special case f HJM wih HJM s ( generaed recursively by (5.3. In BGM, ( is usually aken be nnrandm; he resuling ( is randm. Remark 34.2 (5.4 (equivalenly, (5.4 is a schasic parial differenial equain because f he L( erm. This is n as errible as i firs appears. Reurning he HJM variables and T, se K( T L( T ; Then and (5.4 and (5.4 becme dk( T dl( T ; ; L( T ; d dk( T ( T ; K( T [ ( T ; + d + dw (] K( T ( T ; ( T ; K( T ( T ; d + d + dw ( +K( T (6. Remark 34.3 Frm (5.3 we have If we le #, hen and s We saw befre (eq. 4.2 ha as #, ( L( [ + L( ] ( + ; ( ( L(! ( + ( ( T ; K( T!( T ; L(!r( f ( +

342 s K( T!f ( T Therefre, he limi as # f (6. is given by equain (2.5 df ( T ( T ; [ ( T ; d + dw (] Remark 34.4 Alhugh he d erm in (6. has he erm 2 ( T ;K 2 ( T +K( T his equain d n lde because 2 ( T ; K 2 ( T +K( T 2 ( T ; K 2 ( T K( T 2 ( T ; K( T invlving K 2, sluins 34.7 Bnd prices Le ( nr r(u du Frm (2.6 we have The sluin B( T ( B( T d ( [;r(b( T d + db( T ] ( ; B( T ( T ; dw ( ( his schasic differenial equain is given by B( T (B( T ; (u T ; u dw (u ; 2 This is a maringale, and we can use i swich he frward measure Girsanv s Therem implies ha IP T (A B( T A A (T dip B(T T dip 8A 2F(T (T B( T ( (u T ; u 2 du W T ( W ( + (u T ; u du T is a Brwnian min under IP T.

CHAPTER 34. Brace-Gaarek-Musiela mdel 343 34.8 Frward LIBOR under mre frward measure Frm (6. we have dk( T ( T ; K( T [ ( T ; + d + dw (] ( T ; K( T dw T + ( s and K( T K( T K(T TK( T K( T ( ( (u T ; u dw T + (u ; 2 (u T ; u dw T + (u ; 2 (u T ; u dw T + (u ; 2 2 (u T ; u du 2 (u T ; u du 2 (u T ; u du (8. We assume ha is nnrandm. Then X( (u T ; u dw T + (u ; 2 (u T ; u du (8.2 2 is nrmal wih variance and mean ; 2 2 (. 2 ( 2 (u T ; u du 34.9 Pricing an ineres rae caple Cnsider a flaing rae ineres paymen seled in arrears. A ime T +, he flaing rae ineres paymen due is L(T K(T T he LIBOR a ime T. A caple precs is wner by requiring him pay nly he cap c if K(T T >c. Thus, he value f he caple a ime T + is (K(T T ; c +. We deermine is value a imes T +. Case I T T +. C T + ( IE ( (T + (K(T T ; c+ F( ( F( (K(T T ; c + IE (T + (K(T T ; c + B( T + (9.

344 Case II T. Recall ha IP T + (A A (T + dip 8A 2F(T + where We have ( C T + ( IE B( T + (B( T + ( (T + (K(T T ; c+ F ( 2 3 (B( T + B( T + IE B(T + T + B( T + 6 (K(T T ; c + F( {z } 4(T + B( T + 7 {z } 5 (T + ( B( T + IE T + (K(T T ; c + F( Frm (8. and (8.2 we have K(T TK( T fx(g R where X( is nrmal under IP T + wih variance 2 T ( Furhermre, X( is independen f F(. 2 (u T ; u du and mean ; 2 2 (. C T + ( B( T + IE T + (K( T fx(g;c + F( Se h g(y IE T + (y +i fx(g;c yn ( lg y c + 2 ( ; cn ( lg y c ; 2 ( Then C T + ( B( T + g(k( T T ; (9.2 In he case f cnsan,wehave ( p T ; and (9.2 is called he Black caple frmula.

CHAPTER 34. Brace-Gaarek-Musiela mdel 345 34. Pricing an ineres rae cap Le T T T 2 2 T n n A cap is a series f paymens (K(T k T k ; c + a ime T k+ k n; The value a ime f he cap is he value f all remaining caples, i.e., C( X kt k C Tk ( 34. Calibrain f BGM The ineres rae caple c n L( T a ime T + has ime-zer value where g (defined in he las secin depends n C T + ( B( T + g(k( T 2 (u T ; u du Le us suppse is a deerminisic funcin f is secnd argumen, i.e., Then g depends n ( ( 2 (T ; u du 2 (v dv If we knw he caple price C T + (, we can back u he squared vlailiy R T 2 (v dv. Ifwe knw caple prices C T +( C T +( C Tn+( where T <T <<T n, we can back u 2 (v dv T 2 (v dv 2 (v dv ; 2 (v dv n T n; 2 (v dv (. In his case, we may assume ha is cnsan n each f he inervals ( T (T T (T n; T n

346 and chse hese cnsans make he abve inegrals have he values implied by he caple prices. If we knw caple prices C T + ( fr all T, we can back u R T 2 (v dv and hen differeniae discver 2 ( and ( p 2 ( fr all. T implemen BGM, we need bh (, and ( < Nw ( is he vlailiy a ime f a zer cupn bnd mauring a ime + (see (2.6. Since is small (say year, and <, i is reasnable se 4 ( < We can nw slve (r simulae ge L( r equivalenly, K( T T using he recursive prcedure ulined a he sar f Secin 34.6. 34.2 Lng raes The lng rae is deermined by lng mauriy bnd prices. Le n be a large fixed psiive ineger, s ha n is 2 r 3 years. Then ( n D( n ny k ny k ( k (k; [ + L( (k ; ] where he las equaliy fllws frm (4.. The lng rae is 34.3 Pricing a swap Le T be given, and se n lg D( n n nx k lg[ +L( (k ; ] T T + T 2 T +2 T n T + n

CHAPTER 34. Brace-Gaarek-Musiela mdel 347 The swap is he series f paymens (L(T k ; c a ime T k+ k n; Fr T, he value f he swap is n; X ( IE (T k+ (L(T k ; c Nw s +L(T k We cmpue L(T k k B(T k T k+ B(T k T ; k+ ( IE (T k+ (L(T k ; c ( IE (T k+ IE 2 6 4 F( F ( B(T k T k+ ; ; c F( {z } B(T k T k+ ( (T k B(T k T IE (Tk k+ (T k+ F(T k ( IE (T k+ The value f he swap a ime is n; X ( IE (T (L(T k ; c k+ k n; X k F( ; ( + cb( T k+ B( T k ; ( + cb( T k+ F( [B( T k ; ( + cb( T k+ ] 3 7 F ( 5 ; ( + cb( T k+ B( T ; ( + cb( T +B( T ; ( + cb( T 2 ++ B( T n; ; ( + cb( T n B( T ; cb( T ; cb( T 2 ; ; cb( T n ; B( T n The frward swap rae w T ( a ime fr mauriy T is he value f c which makes he ime- value f he swap equal zer w T ( B( T ; B( T n [B( T ++ B( T n ] In cnras he cap frmula, which depends n he erm srucure mdel and requires esimain f, he swap frmula is generic.