Affine credit risk models under incomplete information

Size: px
Start display at page:

Download "Affine credit risk models under incomplete information"

Transcription

1 Affine credi risk models under incomplee informaion Rüdiger Frey 1, Cecilia Prosdocimi 2, Wolfgang J. Runggaldier 2 1 Deparmen of Mahemaics, Universiy of Leipzig Leipzig, Germany (ruediger.frey@mah.uni-leipzig.de). 2 Diparimeno di Maemaica Pura ed Aplicaa, Universiá di Padova, Via Triese 63, Padova, Ialy (cecilia.prosdocimi@gmail.com, runggal@mah.unipd.i). We consider he problem of compuing some basic quaniies such as defaulable bond prices and survival probabiliies in a credi risk model according o he inensiy based approach. We le he defaul inensiies depend on an exernal facor process ha we assume is no observable. We use sochasic filering o successively updae is disribuion on he basis of he observed defaul hisory. On one hand his allows us o capure aspecs of defaul conagion (informaion-induced conagion). On he oher hand i allows us o evaluae he above quaniies also in our incomplee informaion conex. We consider in paricular affine credi risk models and show ha in such models he nonlinear filer can be compued via a recursive procedure. This hen leads o an explici expression for he filer ha depends on a finie number of sufficien saisics of he observed inerarrival imes for he defauls provided one chooses an iniial disribuion for he facor process ha is of he Gamma ype. Key words: Credi risk, affine models, incomplee informaion, pricing of credi derivaives, nonlinear filering, finie dimensional filers. 1. Inroducion We consider he problem of evaluaing some basic quaniies such as defaulable bond prices and survival probabiliies in a credi risky environmen under incomplee informaion on he underlying model. We use he reduced-form or inensiy-based approach o credi risk wih defaul imes 1

2 2 modeled as he jump imes of a doubly sochasic Poisson process. In his model class defaul inensiies are driven by a common facor process X ; his is a convenien way for generaing dependence beween defaul evens of differen firms. Typically i is assumed ha X is no direcly observable, and his will also be he main seing here. In ha case he disribuion of X can be updaed on he basis of he observed defaul hisory and his leads o wha may be ermed as informaion-induced conagion (see e.g. 8]). In he nex secion 2 we discuss our underlying model and describe hree examples of basic quaniies in a credi risky environmen (prices of defaulable bonds wih and wihou recovery and survival probabiliies) ha we firs evaluae under he assumpion of full knowledge also of he facor process wih he main purpose of moivaing he filering problem ha arises in he evaluaion of hese same quaniies when only he defauls are observable, bu no he values of X. In secion 3 we discuss he filering problem and presen is general soluion. Alhough in explici form, his general soluion will in mos cases be difficul o compue and so he ineres arises o consider paricular classes of models, for which he soluion can acually be compued. One such class corresponds o he case when X is a finie-sae Markov chain and his is discussed in 4] in a more general conex. In his paper we shall concenrae on he so-called affine case and his is he subjec of secion 4, where we show ha he filer can be compued via a recursive algorihm. The acual compuaion of his recursive algorihm is discussed in he las secion 5, where we also show ha for a suiable choice of he iniial disribuion of he facor process he filer can be compued as an explici funcion of a finie number of sufficien saisics of he observed defaul inerarrival imes. In his paper we assume ha he informaion available o an agen comes only from observing he defaul hisory. More general informaion srucures can be envisaged such as in he case when agens can observe also prices of defaulable bonds. This generalizaion is considered in 4], bu no for he affine case. 2. The model and some basic examples Consider a marke wih m firms ha may defaul and denoe by τ j he defaul ime of firm j {1,, m. The defaul sae of he porfolio can be summarized by he defaul indicaor process (1) Y = (Y,1,, Y,m ) 0 wih Y,j = 1 {{τj Given a filered probabiliy space (Ω, F, G, P), all processes will be G adaped and τ j is an G sopping ime. Our inensiy-based defaul model implies ha here are no common defauls among he firms, so ha

3 3 we may inroduce he ordered defaul imes 0 = T 0 < T 1 < < T m. One may hen also consider wha can be called he defaul-ideniy process ξ n ha denoes he ideniy of he firm defauling a T n. The defaul observaion hisory H G can hen be given he following wo equivalen represenaions (2) H = σ{y s ; s = σ{(t n, ξ n ) ; T n The facor process X R may be any Markov process (a specific such process will be considered below, see (6)). We assume ha defaul imes are condiionally independen, doubly sochasic random imes (see 7], Secion 9.6); he defaul inensiy of firm j a ime is given by λ j (X ) for some funcion λ j : R (0, ). Formally, his means ha defaul imes are independen given F X = σ(x : 0) wih condiional survival probabiliies given by P(τ j > F X ) = exp ( λ j (X s )ds ) The affine case We shall say ha we are in he affine case if X saisfies a diffusion equaion (3) dx = µ(x )d + σ(x )dw wih (4) µ(x ) = α X + β σ 2 (X ) = γ X + δ Furhermore, assuming for sake of generaliy ha also he shor rae is driven by he facor process, i.e. r = r(x ), (5) λ j (X ) = λ j X, λ j > 0 r(x ) = r X, r > 0 In paricular, for he process X we shall consider a Cox-Ingersoll-Ross (CIR)-ype model, i.e. (6) dx = (a bx ) d + σ X dw wih a, b, σ > 0 and a σ2 2 so ha X > 0 a.s. which, by (5), will hen also imply ha λ j (X ) > 0, r(x ) > 0 a.s. For his affine case in wha follows

4 4 we shall be able o derive explici expressions boh in he case of full as well as of parial informaion. To his effec we recall here he following proposiion, which in is general form can e.g. be found in 6], secion by making he following idenificaions : = T, λ = β, µ = λ, ψ(t ) = B(, T), aφ(t ) = A(, T). The paricular case for β = 0 can also be found in 7], secion The derivaion is based on he Kolmogorov equaion for funcionals of Markov processes. Proposiion 1. Le X saisfy (6) and define (7) F(, x) := E,x {e βx exp T λx s ds for a generic β 0 and λ > 0. In he presen affine case his funcion F(, x) admis he represenaion (8) F(, x) = exp A(, T) B(, T)x] where, for given T, he funcions A(, T), B(, T) saisfy he following firs order ordinary differenial equaions in 0, T] (9) B (, T) = b B(, T) σ2 B 2 (, T) λ, B(T, T) = β A (, T) = B(, T), A(T, T) = 0 and hey have as soluions (10) B(, T) = βγ + b + eγ(t ) (γ b)] + 2 λ(e γ(t ) 1) βσ 2 (e γ(t ) 1) + γ b + e γ(t ) (γ + b) A(, T) = 2a σ log 2γe (T )(γ+b) 2 2 βσ 2 (e γ(t ) 1) + γ b + e γ(t ) (γ + b) wih γ := b 2 + 2σ 2 λ. 2.2 Examples Of he following hree examples he firs wo concern pricing under full informaion and so he underlying probabiliy measure P has o be seen as a pricing (maringale) measure. The hird one concerns survival probabiliies and here P represens hen he hisorical/real world probabiliy measure. The basic quaniies in hese hree examples may be considered as building blocks for more imporan credi risky producs. ]

5 Example 1. Defaulable zero-coupon bond on firm j wih mauriy T and zero recovery Using sandard resuls for pricing defaulable claims in models wih doubly-sochasic defaul imes (see e.g. 7], secion 9.4.3) he price a ime T of a zero recovery bond on firm j can be expressed as (11) where p j (, T) = E (e ) T r(x s )ds (1 Y T,j ) G ( = (1 Y,j )E X e ) T R j (X s ) ds := Π j 1 (X, Y ) (12) R j (X ) := r(x ) + λ j (X ) I is hus a funcion Π j 1 (X, Y ), paramerized by, T ha for simpliciy we drop from he noaion, of he curren values of he facor and he defaul indicaor processes. From he previous secion 2.1 i is easily seen ha in he affine case he funcion Π j 1 (X, Y ) akes he following exponenially affine form (13) Π j 1 (X, Y ) = (1 Y,j ) exp α j (, T) β j (, T) x ] where, for X saisfying he CIR model (6), he coefficiens in (13) are given by he formulae in (10) wih β j (, T) given by he expression for B(, T) here and α j (, T) by ha for A(, T). Furhermore, for he presen case he coefficiens in he righ hand side of (10) have o be chosen as follows (we may consider, T as fixed): a, b, σ come from (6), β = 0, λ = λ j + r; γ = b2 + 2σ 2 λ Example 2. Recovery paymen Denoe by Z j τ j 1 {τj T he recovery paymen a he ime τ j of defaul of he j h firm, where Z j X is an F adaped process. I is well-known ha he value in of he recovery paymen is given by (14) ( (1 Y,j )E e ) τ j r(x s )ds Z j τ j 1 {τj T G = (1 Y,j )E X ( T Zsλ j j (X s )e s R(X u ) du ds ) := Π j 2 (X, Y ) ; see again 7] for a proof. From hese building blocks he price of many credi derivaives is obained in a sraighforward manner. For insance, he price of a zero-coupon

6 6 bond wih recovery is simply given by he sum of he price wihou recovery and he value of he recovery paymen, and also spreads of credi defaul swaps are easily compued. In he affine case also Π j 2 (X, Y ) can be given a more explici form ha is parly of he exponenially affine ype. We do no discuss his in deail here referring he reader o 7], secion or direcly o he original paper 3] Example 3. Survival probabiliies As already menioned, in his example he underlying probabiliy P is he hisorical/real world probabiliy measure. We wan in fac o compue he probabiliy, given our informaion, ha firm j does no defaul prior o a given ime T. A similar argumen as in he derivaion of (11) immediaely gives ] (15) P ( τ j > T G ) = (1 Y,j ) E X {exp T λ j (X s )ds := Π j 3 (X, Y ) Noice ha he expression of Π j 3 (X, Y ) is compleely analogous o ha of Π j 1 (X, Y ) in (11) so ha in he affine case i can be given an expression of he exponenially affine form like Π j 1 (X, Y ) in (13). 3. Incomplee informaion (he filering problem) In he examples of he previous secion we have seen ha, under full informaion also of he facor process X, he values of he basic quaniies of ineres can be expresses as an explici funcion Π(X, Y ) of he curren values of he facor and he defaul indicaor processes. If agens do no have access o he full informaion represened by he filraion G, bu only o ha corresponding o H, i.e. he informaion represened by he defaul hisory hen, based on ieraed condiional expecaions, i appears naural o consider as corresponding values for he basic quaniies he following (16) Π (Y ) := E { Π(X, Y ) H where he expecaion is under he measure P ha is a pricing (maringale) measure in he case of he firs wo examples and he hisorical measure in he hird. For he hird example here is no problem wih he definiion (16), bu in he case of he pricing examples one has o make sure ha (16) leads o arbirage-free prices. To his effec we can sae he following simple lemma Lemma 2. Assume he shor rae is H adaped, in paricular deerminisic. Then, aking as numeraire he money-marke accoun B := exp r 0 sds ], formula (16) leads o arbirage-free prices in he sense ha B 1 Π (Y ) is a (P, H ) maringale.

7 7 Proof : Le s ; hen E { Π (Y ) B H s = E { E { Π(X,Y ) B H Hs = E { Π(X,Y ) B H s = E { E { Π(X,Y ) B G s Hs = E { Π(X s,y s ) B s H s = Π s (Y s ) B s where we have used he fac ha P is a pricing (maringale) measure for he numeraire B in he sense ha B 1 Π(X, Y ) is a P maringale in he full filraion G. 3.1 The filering problem We have seen ha he problem of compuing values of risk-sensiive producs under incomplee informaion abou he facor process amouns o ha of compuing condiional expecaions as in (16). Since Y H, in wha follows we shall for simpliciy drop he dependence on Y so ha he righ hand side of (16) becomes of he form E { f (X ) H where f ( ) is a generic (bounded) funcion of he facor process. Denoing by π (dx) he condiional disribuion of X given H, (16) leads hen o he problem of compuing (17) π ( f ) := E { f (X ) H = f (x)π (dx) which is a nonlinear filering problem of a diffusion process, given poin process observaions. 3.2 General soluion of he filering problem Le us firs inroduce he global defaul inensiy m (18) λ(x, Y ) := (1 Y,j ) λ j (X ) j=1 which is he sum of he defaul inensiies of he sill surviving firms. Noe ha, for T n < T n+1, λ(x, Y ) is he inensiy of T n+1. In deriving our filering resuls we disinguish he cases beween defaul imes and a a defaul ime Filer beween defauls Le T n be he generic n h defaul ime and le T n, T n+1 ). I follows from he general filering equaions (Kushner-Sraonovich equaions) for poin process observaions by he so-called innovaions mehod (see 5], 2]) ha for π ( f ) as defined in (17) one has (19) π ( f ) = π Tn ( f ) + T n πs (L f ) π s ( λ f ) + π s ( λ) π s ( f ) ] ds

8 8 where L is he generaor ha corresponds o he diffusion process X. Furhermore, noicing ha for T n, T n+1 ) one has Y = Y Tn, whenever T n, T n+1 ) we shall consider λ as a funcion of x alone, i.e. (20) λ(x) = λ(x, Y Tn ) wih λ(x, y) as in (18). Proposiion 3. Assume he condiions for he uniqueness of he soluion of (19) as described e.g. in Appendix 2 of 2] hold. The soluion o (19) is hen for T n, T n+1 ) and an inegrable (bounded) f given by (21) π ( f ) = ϱ ( f ) ϱ (1) where ϱ ( f ), he unnormalized condiional expecaion of f, can be obained as (22) ϱ ( f ) = ψ (T n, x)( f ) π Tn (dx) wih (23) ψ (T n, x)( f ) := E Tn,x and λ(x ) according o (20). { f (, X ) e λ(x s )ds Tn Remark 4. This proposiion shows ha, beween defauls, he filer evolves deerminisically and is evoluion is deermined by he Markovian semigroup ψ (T n, x)( f ) in (23). Proof : Alhough he proof can be obained from ha of Proposiion 3.4 in 2], we presen here a direc derivaion. By Io s formula we have f (, X ) e λ(x s Tn )ds f (T n, X Tn ) = e s λ(x u )du Tn T (L f )(s, Xs ) λ(x s ) f (s, X s ) ] ds n + e s λ(x u Tn )du σ f X T s n x (s, X s)dw s Assuming he condiions for applying Fubini s heorem are saisfied, le us ake on he lef and righ hand sides he expecaion condiional on X Tn = x

9 9 hus obaining { f (, X )e λ(x s )ds Tn f (T n, x) E Tn,x = E T Tn,x {e s λ(x u )du Tn (L f )(s, Xs ) λ(x s ) f (s, X s ) ] ds n Using he definiion of ψ in (23) his las relaion can be rewrien as ψ (T n, x)( f ) ψ Tn (T n, x)( f ) = T n ψs (T n, x) (L f ) ψ s (T n, x) ( λ f ) ] ds Inegraing wih respec o π Tn (dx), aking ino accoun (22) and applying once more Fubini, one arrives a ϱ ( f ) ϱ Tn ( f ) = T n ϱs (L f ) ϱ s ( λ f ) ] ds From here one obains hen immediaely ( ) ϱ ( f ) dπ ( f ) = d = π (L f ) d π ( λ f ) d + π ( f ) π ( λ) d ϱ (1) from which he resul follows by he assumed uniqueness of he soluion of (19) Filer a a defaul Consider now a generic defaul ime T n. Again from he general filering equaions of he innovaions approach (5],2]) one has (24) π Tn ( f ) = π T n (λ ξ n f ) π T n (λ ξn ), where we implicily use ha π T n (λ ξn ) > 0 a.s. The expression π T n ( f ) denoes here he lef hand limi of π Tn ( f ) in T n, which exiss by (19). Concluding his secion 3.2 one sees ha he crucial poin o obain a soluion of he filering equaions is he possibiliy of explicily compuing he semigroup ψ in (23) and in he nex Secion 4 we shall address his issue in he case of an affine model. 4. Filering in affine models As menioned a he end of he previous secion, in his secion we shall derive an explici represenaion of he semigroup ψ in (23) for affine models and his will hen lead o an explici soluion of he filering problem.

10 10 Recall ha for an affine model we have posulaed an affine dynamics for X ha we ake here as given by he CIR model (6). Furhermore, as in (5), we shall assume λ j (X ) = λ j X, which by (18) and (20) hen implies ha also λ(x, Y ) = λ(y ) X ; for he shor rae we assume ha i is consan, i.e. r(x ) r. Finally, we shall assume f (, x) o be exponenially affine, analogously o F(, x) in (8), namely of he form (25) f (, x) = exp α(, T) β(, T)x ] and is specific form depends on he specific problem a hand (see he examples in secion 2.2). Noice nex ha wih f (, x) of he form (25) he semigroup ψ becomes (26) ψ (T n, x)( f ) = E Tn,x = e α(,t) E Tn,x {e α(,t) β(,t)x e λ Tn X sds {e β(,t)x e λ Tn X s ds where we have simply wrien λ for λ(y ) since Y = Y Tn. The crucial quaniy becomes herefore he second facor in he righmos expression in (26). We have now he following proposiion, he proof of which follows immediaely from Proposiion 1 Proposiion 5. Under he assumpions of his secion one has (27) E Tn,x {e β(,t)x e λ Tn X s ds = exp A(T n, ; β) B(T n, ; β)x ] β=β(,t) where A(T n, ; β) and B(T n, ; β) are given by (28) A(T n, ; β) = 2a σ log 2γe ( Tn)(γ+b) 2 2 βσ 2 (e γ( Tn) 1) + γ b + e γ( Tn) (γ + b) B(T n, ; β) = βγ + b + eγ( T n) (γ b)] + 2 λ(e γ( T n) 1) βσ 2 (e γ( T n) 1) + γ b + e γ( T n) (γ + b) wih γ := b 2 + 2σ 2 λ 4.1 Filer beween defauls Combining Proposiions 3 and 5 and noicing ha he consan 1 can be expressed as he funcion f (, x) in (25) for α(, T) = β(, T) = 0, one immediaely obains

11 11 Proposiion 6. The filer π ( f ) is, for T n, T n+1 ), given by e B(T n,;β(,t))x π Tn (dx) (29) π ( f ) = η(t n ;, T) e B(T n,;0)x π Tn (dx) where (30) η(t n ;, T) = e α(,t) exp A(T n, ; β(, T)) A(T n, ; 0) ] Remark 7. I follows from (28) ha B(T n, ; β) is, for T n, nonnegaive whenever β β() where 2 λ(e (31) β() γ( Tn) 1) = min γ + b + e γ( Tn) (γ b), γ b + ] eγ( Tn) (γ + b) σ 2 (e γ( Tn) 1) Thus β() is sricly posiive for > T n and i is equal o zero for = T n. From he examples in secion 2.2 i follows ha he value of β(, T) o be used for he parameer β in he numeraor of he righ hand side in (29) (as well as in (34) below) is sricly posiive, while in he denominaor of he same formulae i is equal o zero. Consequenly he corresponding value of B(T n, ; β) is always nonnegaive. In (34) below we shall also need he derivaive of B(T n, ; β) wih respec o β, evaluaed a β = β(t n+1, T) and a β = 0. Since T n+1 > T n and hus β(t n+1 ) < 0, no only β = β(t n+1, T), bu also β = 0 are inerior poins of he domain of posiiviy for B(T n, ; β) and in his domain i is easily seen from (28) ha B(T n, ; β) is differeniable wih respec o β. Recall now also ha X as given by (6) for a σ2 2 has is a-priori suppor in he posiive half line and his implies ha also all condiional disribuions π (dx) have heir suppor in he posiive half line. I follows ha he wo inegrals in he righ hand side of (29) are well defined and finie. From he above Remark 7 we have ha, for all values of ineres for β, B(T n, ; β) is sricly posiive. Given ha X > 0 anyway, his leads us o define he momen generaing funcion of he condiional disribuion π (dx) of X by ( ) (32) χ (φ) := π e φ X φ > 0 From Proposiion 6 one hen obains immediaely Corollary 8. For T n, T n+1 ) and f (, x) as in (25) we have ha he filer value π ( f ) is given by (33) π ( f ) = η(t n ;, T) χ T n (B(T n, ; β(, T)) χ Tn (B(T n, ; 0)) where η( ) is as in (30), χ Tn ( ) as in (32) and B(T n, ; β) is given by (28).

12 Filer a a defaul ime Based on he general formula (24) for he filer a a generic defaul ime T n, we can now show he following (since for he case beween defauls we considered he inerval T n, T n+1 ), here we ake T n+1 o denoe he generic defaul ime) Proposiion 9. Assuming he condiions are fulfilled o differeniae under he inegral sign, a he generic defaul ime T n+1 we have (34) π Tn+1 ( f ) = e α(t β e A(T n,t n+1 ;β) χ Tn (B(T n, T n+1 ; β)) ] β=β(tn+1,t) n+1,t) β e A(T n,t n+1 ;β) χ Tn (B(T n, T n+1 ; β)) ] β=0 where A(T n, T n+1 ; β) and B(T n, T n+1 ; β) are given by (28), β(t n+1, T) corresponds o he exponenially affine represenaion of f (, x) in (25) and χ Tn ( ) is as defined in (32). Proof : Saring from (24) for he defaul ime T n+1 and noicing ha π T n+1 ( f ) is he limi, for T n+1, of he filer π ( f ) when T n, T n+1 ), using Proposiion 3 wih f (, x) as in (25) and wih he semigroup ψ ( ) as in (26) combined wih Proposiion 5, we obain he following sequence of equaliies, where we use differeniaion under he inegral sign in wo insances π Tn+1 ( f ) = π T n+1 {λ ξ n+1 X Tn+1 f (T n+1,x Tn+1 ) π T n+1 {λ ξn+1 X Tn+1 = e α(t n+1,t) E Tn,x { = e α(t n+1,t) = e α(t n+1,t) X Tn+1 e β(t n+1,t) X T λ T n+1 n+1 e Tn ETn,x {X Tn+1 e λ T n+1 Xudu Tn π Tn (dx) { β E Tn,x e β X T λ T n+1 n+1 e Tn { β E Tn,x e β X T λ T n+1 n+1 e Tn Xu du Xu du π Tn (dx) π Tn (dx) β= β(t n+1,t) Xu du π Tn (dx) β=0 β {expa(t n,t n+1 ;β) B(T n,t n+1 ;β) x] β=β(tn+1,t) π Tn (dx) β{expa(t n,t n+1 ;β) B(T n,t n+1 ;β) x] β=0 π Tn (dx) = e α(t n+1,t) βe A(Tn,T n+1 ;β) χ Tn (B(T n,t n+1 ;β))] β=β(tn+1,t) βe A(Tn,T n+1 ;β) χ Tn (B(T n,t n+1 ;β))] β=0 Noice ha, by he consideraions in Remark 7 all he above quaniies are well defined.

13 13 Remark 10. Noice ha in he filer updae a T n+1 only he informaion abou he iming T n+1 of he (n + 1) s defaul is being used and no also ha of he defauling firm ξ n+1 ; in fac, he facor λ ξn+1, which conains his informaion, drops ou due o he normalizaion. This happens however only because of our affine form of he dependence of he defaul inensiy on he facor process and he fac ha he laer is aken o be scalar. In any case, he informaion abou ξ n+1 is no los as i becomes par of Y Tn+1. From Corollary 8 and Proposiion 9 i now follows ha he crucial quaniy for compuing he filer in he affine case is he momen generaing funcion χ (φ) in (32), which we need o compue only for he values of corresponding o he defaul imes and for any φ > 0. Noice also ha χ (φ) is nohing bu he filer π ( f ) for f (, x) when his laer funcion has he exponenially affine form in (25) wih α(, T) = 0 and β(, T) = φ. Combining hese remarks wih he resuls of Proposiion 9 one obains immediaely Corollary 11. Under he assumpions of Proposiion 9 one has β e A(T n,t n+1 ;β) χ Tn (B(T n, T n+1 ; β)) ] β=φ (35) χ Tn+1 (φ) = e A(T n,t n+1 ;β) χ Tn (B(T n, T n+1 ; β)) ] β=0 β On he basis of he previous resuls we can now wrie down an algorihm o compue recursively he filer in he affine case, which we do in he nex subsecion. 4.3 Filer algorihm i) Sar a T 0 = 0 wih a given χ 0 (φ). ii) A he generic T n+1 compue (see Corollary 11) β e A(T n,t n+1 ;β) χ Tn (B(T n, T n+1 ; β)) ] β=φ χ Tn+1 (φ) = e A(T n,t n+1 ;β) χ Tn (B(T n, T n+1 ; β)) ] β=0 β iii) For T n, T n+1 ) and wih f (, x) = expα(, T) β(, T)x] he filer is hen given by (see Corollary 8) π ( f ) = η(t n ;, T) χ T n (B(T n, ; β(, T)) χ Tn (B(T n, ; 0)) where η(t n ;, T) is given in (30) and B(T n, ; β) in (28). iv) For = T n+1 he filer is (see Proposiion 9 and Corollary 11) π Tn+1 ( f ) = e α(t n+1,t) χ Tn+1 (β(t n+1, T))

14 14 Sep ii) above is a recursive formula o compue he momen generaing funcion χ Tn ha hen leads o an explici filer soluion according o Seps iii) and iv). Alhough recursive in naure, his Sep ii) may become increasingly difficul o compue since he analyic expression for χ Tn (φ) migh become more and more involved wih every sep. A naural quesion hen arises o see wheher, wih a suiable choice of he iniial χ 0 (φ), he recursions remain a a level ha allows for feasible compuaions. This will be he subjec of he nex Secion 5, where we also show ha, alhough a priori he filer is no finie-dimensional in he radiional sense, for a suiable choice of he iniial disribuion i can be paramerized by a finie number of sufficien saisics. Remark 12. We conclude his secion by poining ou ha our filer resuls may have a wider scope han only wha we have described here by showing ha he filer can be compued by compuing he condiional momen generaing funcion and ha his can be done recursively according o (35). In fac, since his condiional momen generaing funcion is relaed o he Laplace ransform of he condiional (filer) disribuion, one can, a leas in heory, inver his condiional momen generaing funcion hereby recovering he filer disribuion iself. This would hen allow o compue condiional expecaions no only of exponenially affine funcions of he facor process bu also of any inegrable funcion hereof. 5. Finie dimensional compuaion of he filer In his secion we exhibi a choice for χ 0 (φ) ha allows Sep ii) in he filer algorihm of he previous subsecion 4.3 o remain compuable a every sep. For his purpose recall ha he recursions in Sep ii) correspond o he recursive formula (35), where he coefficiens A(T n, T n+1 ; β), B(T n, T n+1 ; β) are given in (28). Inroduce he shorhand noaions (36) R n+1 = γ + b + e γ(t n+1 T n ) (γ b) U n+1 = σ ( 2 e γ(t n+1 T n ) 1 ) S n+1 = 2 λ ( e γ(t n+1 T n ) 1 ) V n+1 = γ b + e γ(t n+1 T n ) (γ + b) W n+1 = 2 γ e (T n+1 Tn) (γ+b) 2 ha, since γ = b σ 2 λ, are all posiive quaniies. The coefficiens A(T n, T n+1 ; β) and B(T n, T n+1 ; β) from (28) can hen be given he following

15 15 represenaion (37) A(T n, T n+1 ; β) = 2 a log ( ) W n+1 σ 2 β U n+1 +V n+1 B(T n, T n+1 ; β) = β R n+1+s n+1 β U n+1 +V n+1 Choosing as disribuion for he iniial value X 0 of he facor process a specific Gamma-ype disribuion, we shall now prove he following heorem ha gives an explici compuable represenaion for χ (φ) a he various defaul imes = T n. Theorem 13. Le (38) χ 0 (φ) = ( φ ) 2a σ 2 = (1 + φ) 2a σ 2, φ > 0 which according o (32) corresponds o he momen generaing funcion of a Gamma disribuion for X 0 wih parameers ( 2a σ 2, 1 ). Then (39) χ Tn (φ) = c n (φ H n + K n ) 2a σ 2 n p n (φ) where H n and K n saisfy he recursions H n = R n H n 1 + U n K n 1, H 0 = 1 (40) K n = S n H n 1 + V n K n 1, K 0 = 1 he coefficien c n is given by (41) c n = ] K 2a σ 2 n 1 n p n (0) and p n (φ) is a polynomial of degree n 1 given by (42) 1 for n = 0 and n = 1 ( 2a n + 1 ) H σ 2 n (φu n + V n ) p n (φ) p n (φ) = +U n (φh n + K n ) p n (φ) + (φh n + K n )(φu n + V n ) φ p n(φ) wih ( ) φrn (43) p n (φ) = (φu n + V n ) n 2 + S n p n 1, n 2 φu n + V n for n 2

16 16 Proof : The saemen is clearly rue for n = 0. We show i firs for n = 1 and hen inducively for all n 2. i) he case n = 1 : by (35),(36), (37) and he recursions in (40) we have (44) χ T1 (φ) = = ( ) 2a ( ) W 1 σ 2 βr1 +S β βu 1 +V χ 1 T0 1 βu 1 +V 1 β=φ ( ) 2a ( ) W 1 σ 2 βr1 +S β βu 1 +V χ 1 T0 1 βu 1 +V 1 β=0 ( ) 2a ( ) W 1 σ 2 2a βh 1 +K 1 σ 2 β βu 1 +V 1 βu 1 +V 1 β=φ ( ) 2a ( ) W 1 σ 2 2a βh 1 +K 1 σ 2 β βu 1 +V 1 βu 1 +V 1 β=0 = (φh 1+K 1 ) 2a σ 2 1 K 2a σ which indeed corresponds o (39) wih (41) and (42). ii) he general case n 2 : assume (39) holds for n 1. Then, always by (35), (36), (37) and (40) we obain (45) χ Tn (φ) = = { ( ) 2a ( Wn σ 2 βrn+sn ) β χ βun+vn Tn 1 βun+vn β { ( ) 2a Wn βun+vn σ 2 χ Tn 1 ( βrn+sn βun+vn { ( ) 2a ( 1 σ 2 βhn+kn β βun+vn βun+vn β { ( ) 2a 1 βun+vn σ 2 ( βhn+kn βun+vn β=φ ) β=0 ) 2a σ 2 n+1 ( βrn+sn ) p n 1 βun+vn β=φ ) 2a σ 2 n+1 ( βrn+sn ) p n 1 βun+vn Taking ino accoun ha, by (43), ( ) βrn + S n (46) p n 1 = (βu n + V n ) 2 n p n (β) βu n + V n he numeraor in he righmos expression of (45) becomes (47) { β (βun + V n )(βh n + K n ) 2a n+1 p σ 2 n (β) β=φ = (βh n + K n ) 2a σ 2 n{ U n (βh n + K n ) p n (β) β=0 + ( 2a n + 1 ) H σ 2 n (βu n + V n ) p n (β) + (βu n + V n )(βh n + K n ) β p n(β) β=φ = (φ H n + K n ) 2a σ 2 n p n (φ)

17 17 where we have used he definiion of p n (φ) in (42). Reurning o (45) and recalling ha he denominaor in he righmos expression of (45) is he same as he numeraor excep for puing β = 0, one finally obains (48) χ Tn (φ) = (φ H n + K n ) 2a σ 2 n p n (φ) K 2a σ 2 n n p n (0) = c n (φ H n + K n ) 2a σ 2 n p n (φ) Remark 14. I follows from Theorem 13 ha, for a choice of he iniial disribuion corresponding o χ 0 (φ) in (38), he sequence χ Tn (φ) and herefore (see Seps iii) and iv) in Secion 4.3) he enire filer is parameerized by a same finie number of sufficien saisics, namely he pairs (H n, K n ) and he polynomial funcions p n (φ), all of which can be compued recursively on he basis of he funcions R n, S n, U n, V n of he inerarrival imes of he defauls. References 1. Bielecki, T. R., M. Jeanblanc, & M. Rukowski (2004). Modeling and Valuaion of Credi Risk. In : Sochasic Mehods in Finance (M. Frielli and W. Runggaldier, eds), Lecure Noes in Mahemaics. Vol 1856, Springer Verlag. 2. Ceci,C., & A. Gerardi (2006). A Model for High Frequency Daa under Parial Informaion : a Filering Approach. Inernaional Journal of Theoreical and Applied Finance (o appear). 3. Duffie, D., & N. Gârleanu (2001). Risk and Valuaion of Collaeralized Deb Obligaions. Financial Analyss Journal. 57, Frey, R., & W. J. Runggaldier ( 2006). Credi Risk and Incomplee Informaion : a Nonlinear Filering Approach. Preprin. 5. Kliemann,W. H., G. Koch, & F. Marchei (1990). On he Unnormalized Soluion of he Filering Problem wih Couning Process Observaions. IEEE Transacions on Informaion Theory. 36, Lamberon,D., & B. Lapeyre (1995). An Inroducion o Sochasic Calculus Applied o Finance. Chapman and Hall. 7. McNeil A. J., R. Frey, & P. Embrechs (2005). Quaniaive Risk Managemen. Princeon Universiy Press. 8. Schönbucher, P. (2004) Informaion-driven Defaul Conagion. Preprin, Deparmen of Mahemaics, ETH Zürich.

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

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

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

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

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

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

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

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

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

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

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

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

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

Math 10B: Mock Mid II. April 13, 2016 Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.

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

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

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

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

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

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

Notes for Lecture 17-18

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

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

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

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

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

arxiv: v1 [math.pr] 19 Feb 2011

arxiv: v1 [math.pr] 19 Feb 2011 A NOTE ON FELLER SEMIGROUPS AND RESOLVENTS VADIM KOSTRYKIN, JÜRGEN POTTHOFF, AND ROBERT SCHRADER ABSTRACT. Various equivalen condiions for a semigroup or a resolven generaed by a Markov process o be of

More information

6. Stochastic calculus with jump processes

6. Stochastic calculus with jump processes A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio

More information

t 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

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

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

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

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

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,

More information

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

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

5. Stochastic processes (1)

5. Stochastic processes (1) Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

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

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients mahemaics Aricle A Noe on he Equivalence of Fracional Relaxaion Equaions o Differenial Equaions wih Varying Coefficiens Francesco Mainardi Deparmen of Physics and Asronomy, Universiy of Bologna, and he

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

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

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE Topics MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES 2-6 3. FUNCTION OF A RANDOM VARIABLE 3.2 PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE 3.3 EXPECTATION AND MOMENTS

More information

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H. ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple

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

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

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

Convergence of the Neumann series in higher norms

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

More information

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

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

Stochastic models and their distributions

Stochastic models and their distributions Sochasic models and heir disribuions Couning cusomers Suppose ha n cusomers arrive a a grocery a imes, say T 1,, T n, each of which akes any real number in he inerval (, ) equally likely The values T 1,,

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

THE TERM STRUCTURE OF INTEREST RATES IN A MARKOV SETTING

THE TERM STRUCTURE OF INTEREST RATES IN A MARKOV SETTING THE TERM STRUCTURE OF INTEREST RATES IN A MARKOV SETTING Rober J. Ellio Haskayne School of Business Universiy of Calgary Calgary, Albera, Canada rellio@ucalgary.ca Craig A. Wilson College of Commerce Universiy

More information

GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT

GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT Inerna J Mah & Mah Sci Vol 4, No 7 000) 48 49 S0670000970 Hindawi Publishing Corp GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT RUMEN L MISHKOV Received

More information

Lecture 10: The Poincaré Inequality in Euclidean space

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

More information

The Strong Law of Large Numbers

The Strong Law of Large Numbers Lecure 9 The Srong Law of Large Numbers Reading: Grimme-Sirzaker 7.2; David Williams Probabiliy wih Maringales 7.2 Furher reading: Grimme-Sirzaker 7.1, 7.3-7.5 Wih he Convergence Theorem (Theorem 54) 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

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

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

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

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

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

Inequality measures for intersecting Lorenz curves: an alternative weak ordering

Inequality measures for intersecting Lorenz curves: an alternative weak ordering h Inernaional Scienific Conference Financial managemen of Firms and Financial Insiuions Osrava VŠB-TU of Osrava, Faculy of Economics, Deparmen of Finance 7 h 8 h Sepember 25 Absrac Inequaliy measures for

More information

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

More information

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem.

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem. Noes, M. Krause.. Problem Se 9: Exercise on FTPL Same model as in paper and lecure, only ha one-period govenmen bonds are replaced by consols, which are bonds ha pay one dollar forever. I has curren marke

More information

On Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems

On Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems MATHEMATICS OF OPERATIONS RESEARCH Vol. 38, No. 2, May 2013, pp. 209 227 ISSN 0364-765X (prin) ISSN 1526-5471 (online) hp://dx.doi.org/10.1287/moor.1120.0562 2013 INFORMS On Boundedness of Q-Learning Ieraes

More information

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

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering

More information

A New Perturbative Approach in Nonlinear Singularity Analysis

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

More information

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

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

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

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,

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

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

arxiv: v1 [math.ca] 15 Nov 2016

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

More information

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

More information

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

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

More information

III. Module 3. Empirical and Theoretical Techniques

III. Module 3. Empirical and Theoretical Techniques III. Module 3. Empirical and Theoreical Techniques Applied Saisical Techniques 3. Auocorrelaion Correcions Persisence affecs sandard errors. The radiional response is o rea he auocorrelaion as a echnical

More information

A proof of Ito's formula using a di Title formula. Author(s) Fujita, Takahiko; Kawanishi, Yasuhi. Studia scientiarum mathematicarum H Citation

A proof of Ito's formula using a di Title formula. Author(s) Fujita, Takahiko; Kawanishi, Yasuhi. Studia scientiarum mathematicarum H Citation A proof of Io's formula using a di Tile formula Auhor(s) Fujia, Takahiko; Kawanishi, Yasuhi Sudia scieniarum mahemaicarum H Ciaion 15-134 Issue 8-3 Dae Type Journal Aricle Tex Version auhor URL hp://hdl.handle.ne/186/15878

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models Journal of Saisical and Economeric Mehods, vol.1, no.2, 2012, 65-70 ISSN: 2241-0384 (prin), 2241-0376 (online) Scienpress Ld, 2012 A Specificaion Tes for Linear Dynamic Sochasic General Equilibrium Models

More information

A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS

A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS Theory of Sochasic Processes Vol. 14 3), no. 2, 28, pp. 139 144 UDC 519.21 JOSEP LLUÍS SOLÉ AND FREDERIC UTZET A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS An explici procedure

More information

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

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15. SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a

More information

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The

More information

RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY

RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY ECO 504 Spring 2006 Chris Sims RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY 1. INTRODUCTION Lagrange muliplier mehods are sandard fare in elemenary calculus courses, and hey play a cenral role in economic

More information

Representation of Stochastic Process by Means of Stochastic Integrals

Representation of Stochastic Process by Means of Stochastic Integrals Inernaional Journal of Mahemaics Research. ISSN 0976-5840 Volume 5, Number 4 (2013), pp. 385-397 Inernaional Research Publicaion House hp://www.irphouse.com Represenaion of Sochasic Process by Means of

More information

A Bayesian Approach to Spectral Analysis

A Bayesian Approach to Spectral Analysis Chirped Signals A Bayesian Approach o Specral Analysis Chirped signals are oscillaing signals wih ime variable frequencies, usually wih a linear variaion of frequency wih ime. E.g. f() = A cos(ω + α 2

More information

Existence of positive solution for a third-order three-point BVP with sign-changing Green s function

Existence of positive solution for a third-order three-point BVP with sign-changing Green s function Elecronic Journal of Qualiaive Theory of Differenial Equaions 13, No. 3, 1-11; hp://www.mah.u-szeged.hu/ejqde/ Exisence of posiive soluion for a hird-order hree-poin BVP wih sign-changing Green s funcion

More information

RISK-NEUTRAL VALUATION UNDER FUNDING COSTS AND COLLATERALIZATION

RISK-NEUTRAL VALUATION UNDER FUNDING COSTS AND COLLATERALIZATION RISK-NEUTRAL VALUATION UNDER FUNDING COSTS AND COLLATERALIZATION Damiano Brigo Dep. of Mahemaics Imperial College London Andrea Pallavicini Dep. of Mahemaics Imperial College London Crisin Buescu Dep.

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

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

6.2 Transforms of Derivatives and Integrals.

6.2 Transforms of Derivatives and Integrals. SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

18 Biological models with discrete time

18 Biological models with discrete time 8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so

More information

8. Basic RL and RC Circuits

8. Basic RL and RC Circuits 8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics

More information

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017 Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =

More information

Echocardiography Project and Finite Fourier Series

Echocardiography Project and Finite Fourier Series Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every

More information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

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

Estimation of Poses with Particle Filters

Estimation of Poses with Particle Filters Esimaion of Poses wih Paricle Filers Dr.-Ing. Bernd Ludwig Chair for Arificial Inelligence Deparmen of Compuer Science Friedrich-Alexander-Universiä Erlangen-Nürnberg 12/05/2008 Dr.-Ing. Bernd Ludwig (FAU

More information

Empirical Process Theory

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

More information

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

Martingales Stopping Time Processes

Martingales Stopping Time Processes IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765. Volume 11, Issue 1 Ver. II (Jan - Feb. 2015), PP 59-64 www.iosrjournals.org Maringales Sopping Time Processes I. Fulaan Deparmen

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information