Information Premium and Weather Forecasts

Size: px
Start display at page:

Download "Information Premium and Weather Forecasts"

Transcription

1 Informaion Premium and Weaher Forecass Energy Finance Workshop of HUB and Uni DuE Solberg, Harz, Richard Biegler-König Chair for Energy Trading and Finance Universiy of Duisburg-Essen

2 Seie 2/35 Table of conens The Informaion Approach Moivaion The Informaion Premium The Informaion Premium wih Weaher Forecass The Spo Model A simple Temperaure Model Enlargemen of Filraions The Info Premium wih Weaher Forecas Sylised Examples and Muliple Informaion Sylised Examples Muliple Fuure Informaion

3 Seie 3/35 The Informaion Approach The Informaion Approach

4 Seie 4/35 The Informaion Approach Moivaion Spo-Forward Relaionship: Classical heory Wha is he relaionship beween spo and forward prices on elecriciy markes? The classical spo-forward relaionship is where F (, T ) = E Q [S(T ) F ] F is informaion generaed by he spo price up unil ime (i.e. he hisorical filraion) Q is a risk-neural probabiliy The underlying elecriciy is non-sorable Buy-and-hold argumen does no work

5 Seie 4/35 The Informaion Approach Moivaion Spo-Forward Relaionship: Classical heory Wha is he relaionship beween spo and forward prices on elecriciy markes? The classical spo-forward relaionship is where F (, T ) = E Q [S(T ) F ] F is informaion generaed by he spo price up unil ime (i.e. he hisorical filraion) Q is a risk-neural probabiliy The underlying elecriciy is non-sorable Buy-and-hold argumen does no work

6 Seie 5/35 The Informaion Approach Moivaion Sylised and Marke Examples Fuure predicions abou he marke will no affec he curren spo price... bu will affec forward prices Sylized examples: planned ouage of a power plan in one monh weaher forecass Inroducion of CO 2 cerificaes in 2008 wih anicipaed seep increase due o exra coss German Aom Moraorium following he Tōhoku earhquake March 2011 immediae shu-down of 7 nuclear power plans for 3 monhs possible complee shu-down!

7 Seie 5/35 The Informaion Approach Moivaion Sylised and Marke Examples Fuure predicions abou he marke will no affec he curren spo price... bu will affec forward prices Sylized examples: planned ouage of a power plan in one monh weaher forecass Inroducion of CO 2 cerificaes in 2008 wih anicipaed seep increase due o exra coss German Aom Moraorium following he Tōhoku earhquake March 2011 immediae shu-down of 7 nuclear power plans for 3 monhs possible complee shu-down!

8 Seie 6/35 The Informaion Approach Moivaion Example: 2008 CO 2 Emission Coss Inroducion of CO2 fees Spo on Ok. 06 Nov. 06 Dez. 06 Jan. 07 Feb. 07 Mrz. 07 Apr Ok. 07 Nov. 07 Dez. 07 Jan. 08 Feb. 08 Mrz. 08 Apr. 08 Figure: EEX Forward prices observed on 01/10/06 (lef) and 01/10/07 (righ) Typical winer and bank holidays behaviour in boh graphs General upward shif in 2008

9 Seie 7/35 The Informaion Approach Moivaion Example: 2011 German Moraorium - Spo 65,00 EEX Spo (7 days MA) Earhquake ( ) "Moraorium" ( ) Decision o shu down he "7" ( ) End of "Moraorium" ( ) 60,00 55,00 50,00 45,00 40,00 35,00 30, Figure: EEX spo prices 2010/2011 No significan impac on spo prices

10 Seie 8/35 The Informaion Approach Moivaion Example: 2011 German Moraorium - Forward 65 May 2011 Forward Earhquake ( ) "Moraorium" ( ) Decision o shu down he "7" ( ) Figure: EEX May 2011 forward price Permanen increase from 14/03/2011

11 Seie 9/35 The Informaion Approach The Informaion Premium Informaion Approach For elecriciy he hisorical filraion is no sufficien F = σ(s s, s ) Idea: enlarge he filraion!... by addiional fuure informaion a ime T Υ Noe: This echnique has been used o modell insider rading on sock exchanges

12 Seie 9/35 The Informaion Approach The Informaion Premium Informaion Approach For elecriciy he hisorical filraion is no sufficien F = σ(s s, s ) Idea: enlarge he filraion!... by addiional fuure informaion a ime T Υ Noe: This echnique has been used o modell insider rading on sock exchanges

13 Seie 9/35 The Informaion Approach The Informaion Premium Informaion Approach For elecriciy he hisorical filraion is no sufficien F = σ(s s, s ) Idea: enlarge he filraion!... by addiional fuure informaion a ime T Υ Noe: This echnique has been used o modell insider rading on sock exchanges

14 Seie 10/35 The Informaion Approach The Informaion Premium Filraions and Spo-Forward-Relaionship Filraions F - he hisorical filraion H - complee informaion, i.e. H = F σ(s(t Υ )) G - he filraion of all informaion publicly available o he marke Hence, we have he relaion F G H New Spo-Forward-Relaionship: F (, T ) = E Q [S(T ) G ]

15 Seie 10/35 The Informaion Approach The Informaion Premium Filraions and Spo-Forward-Relaionship Filraions F - he hisorical filraion H - complee informaion, i.e. H = F σ(s(t Υ )) G - he filraion of all informaion publicly available o he marke Hence, we have he relaion F G H New Spo-Forward-Relaionship: F (, T ) = E Q [S(T ) G ]

16 Seie 11/35 The Informaion Approach The Informaion Premium The Informaion Premium - Definiion Quanify he influence of fuure informaion using: Informaion Premium The informaion premium is defined o be I(, T ) = E[S T G ] E[S T F ] = F G (, T ) F F (, T ) i.e. he difference beween he prices of he forward under G and F. Noe: no delivery period considered (ye)

17 Seie 12/35 The Informaion Approach The Informaion Premium The Informaion Premium - Propery Lemma The informaion premium is orhogonal o L 2 (F, P). Proof: E[I G (, T ) F ] = E[ E[S(T ) G ] E[S(T ) F ] F ] = 0 In paricular, he resul is valid for all measures equivalen o P Usual mehod o aain he marke price of risk is a measure change This is no possible for he Informaion Premium

18 Seie 12/35 The Informaion Approach The Informaion Premium The Informaion Premium - Propery Lemma The informaion premium is orhogonal o L 2 (F, P). Proof: E[I G (, T ) F ] = E[ E[S(T ) G ] E[S(T ) F ] F ] = 0 In paricular, he resul is valid for all measures equivalen o P Usual mehod o aain he marke price of risk is a measure change This is no possible for he Informaion Premium

19 Seie 13/35 The Informaion Approach The Informaion Premium Empirical Tes for he Informaion Premium Showing he exisence of he info premium is problemaic due o orhogonaliy We developed a es based on Hilber space represenaion... and isolaed a non-measurable premium for he inroducion of cerificaes as well as he Moraorium... ha makes sense in shape and size This is discussed in deail in our paper (o appear)

20 Seie 14/35 The Informaion Premium wih Weaher Forecass The Informaion Premium wih Weaher Forecass

21 Seie 15/35 The Informaion Premium wih Weaher Forecass The Spo Model Spo Price Model We will model he spo using a simple, well known wo-facor model: S() = Λ() + X() + Y () where for a BM W () and Lévy process L() T X(T ) = e α(t ) X() + σ Y (T ) = e β(t ) Y () + T e α(t s) dw (s) e β(t s) dl(s) i.e. a seasonal funcion Λ and wo OU processes Base componen X, spike componen Y

22 Seie 16/35 The Informaion Premium wih Weaher Forecass The Spo Model Forward Price wih Delivery under P and F The forward price wih delivery in [T 1, T 2 ] is hen given by [ ] 1 T2 F (, T 1, T 2 ) = E S(u)du F (App.I) T 2 T 1 and furher calculaions yield T 1 1 F(, T 1, T 2 ) = T 2 T 1 ( T2 Λ(u)du + ᾱ(, T 1, T 2 )X() T 1 + β(, T 1, T 2 )Y () + φ (0) ˆβ(, ) T 1, T 2 ) φ(u) is he log-momen generaing funcion of L ᾱ, β, ˆβ are deerminisic funcions (see App. II)

23 Seie 17/35 The Informaion Premium wih Weaher Forecass A simple Temperaure Model Temperaure Model Given weaher forecass in T Υ > T 2, wha is he informaion premium? We will modell emperaure as a Gaussian OU process where W Z T Z T = e α Z (T ) Z + is a F -BM T + σ Z e α Z (T s) dws Z α Z µ Z (s)e α Z (T s) ds Model fis reasonably well (especially wih seasonal α Z, see various Benh papers)

24 Seie 18/35 The Informaion Premium wih Weaher Forecass A simple Temperaure Model Temperaure and Elecriciy We will assume dw dw Z = ρd Or, in oher words where W S dw = ρdw Z is he spo-specific BM + 1 ρ 2 dw S The spo s base componen is correlaed wih emperaure F is now generaed by W S, W Z and L

25 Seie 19/35 The Informaion Premium wih Weaher Forecass Enlargemen of Filraions Enlargemen of Filraion We will model he addiional informaion using he Theory of Enlargemen of Filraion mainly developed in he 1980s by French mahemaicians (Jeulin, Yor, Jacod) bu also by Proer, Imkeller Iniiaing heorem by Iō (1976) (see App. IV) Iō s heorem provides he decomposiion of he corresponding Lévy process under enlargemen of filraion by is own fuure value Sill, we observe prices/emperaures, no heir random shocks We wan o enlarge by emperaure, i.e. a funcional of a BM

26 Seie 19/35 The Informaion Premium wih Weaher Forecass Enlargemen of Filraions Enlargemen of Filraion We will model he addiional informaion using he Theory of Enlargemen of Filraion mainly developed in he 1980s by French mahemaicians (Jeulin, Yor, Jacod) bu also by Proer, Imkeller Iniiaing heorem by Iō (1976) (see App. IV) Iō s heorem provides he decomposiion of he corresponding Lévy process under enlargemen of filraion by is own fuure value Sill, we observe prices/emperaures, no heir random shocks We wan o enlarge by emperaure, i.e. a funcional of a BM

27 Seie 20/35 The Informaion Premium wih Weaher Forecass Enlargemen of Filraions Yor s Crierion Theorem Le F be he Brownian filraion and G a random variable and define H = F σ(g). Le P G (dl) be he law of G and P G (dl) he condiional law. If here exiss p such ha P G (dl) = p (l)p G (dl) hen he H decomposiion of BM W is where ξ is a H BM. d < p. (G), W. > s W = ξ + 0 p s (G) Noe: laws needed in conrary o Iō s heorem!

28 Seie 21/35 The Informaion Premium wih Weaher Forecass Enlargemen of Filraions Imkeller s Mehod To calculae he quadraic covariaion of p and W we normally need a edious applicaion of Iō s lemma Imkeller proposes o use he Malliavin derivaive: Theorem Wih he same seup as before he H decomposiion of BM W is D s Ps G (G) W = ξ + 0 Ps G (G) ds where D denoes he Malliavin derivaive.

29 Seie 22/35 The Informaion Premium wih Weaher Forecass Enlargemen of Filraions P G and P G We have H = F σ(z TΥ ) = F σ( T Υ e α Z (T Υ s) dw Z s ) We define G = m = s 2 = TΥ e α Z (T Υ s) dw Z s e α Z (T Υ s) dw Z s e 2α Z (T Υ s) ds = 1 2α Z (e 2α Z (T Υ ) e 2α Z T Υ ) so ha P G 1 (dl) = exp( 2πs 1 )dl 2TΥ 2 st 2 Υ P G (dl) = 1 exp( 2π(s 1 (l m ) 2 )dl 2TΥ s 2 ) 2 st 2 Υ s 2 l 2

30 Seie 22/35 The Informaion Premium wih Weaher Forecass Enlargemen of Filraions P G and P G We have H = F σ(z TΥ ) = F σ( T Υ e α Z (T Υ s) dw Z s ) We define G = m = s 2 = TΥ e α Z (T Υ s) dw Z s e α Z (T Υ s) dw Z s e 2α Z (T Υ s) ds = 1 2α Z (e 2α Z (T Υ ) e 2α Z T Υ ) so ha P G 1 (dl) = exp( 2πs 1 )dl 2TΥ 2 st 2 Υ P G (dl) = 1 exp( 2π(s 1 (l m ) 2 )dl 2TΥ s 2 ) 2 st 2 Υ s 2 l 2

31 Seie 22/35 The Informaion Premium wih Weaher Forecass Enlargemen of Filraions P G and P G We have H = F σ(z TΥ ) = F σ( T Υ e α Z (T Υ s) dw Z s ) We define G = m = s 2 = TΥ e α Z (T Υ s) dw Z s e α Z (T Υ s) dw Z s e 2α Z (T Υ s) ds = 1 2α Z (e 2α Z (T Υ ) e 2α Z T Υ ) so ha P G 1 (dl) = exp( 2πs 1 )dl 2TΥ 2 st 2 Υ P G (dl) = 1 exp( 2π(s 1 (l m ) 2 )dl 2TΥ s 2 ) 2 st 2 Υ s 2 l 2

32 Seie 23/35 The Informaion Premium wih Weaher Forecass Enlargemen of Filraions The Informaion Drif We can easily calculae he Malliavin derivaive D P G (l) = P G (l)d ( 1 (l m ) 2 ) 2 st 2 Υ s 2 = P G (l) l m D ( e α Z (TΥ s) dw st 2 Υ s 2 s Z ) 0 = P G (l) l m e α Z (T Υ ) st 2 Υ s 2... by he Malliavin chain rule and properies of Wiener polynomials Thus, using Imkeller s mehod he decomposiion is ξ = W Z = W Z G m e α st 2 s 2 Z (T Υ ) ds Υ 0 0 2α Z e α Z s e 2α Z T Υ e 2α Z s }{{} =a(s) TΥ s e α Z u dw Z u ds

33 Seie 23/35 The Informaion Premium wih Weaher Forecass Enlargemen of Filraions The Informaion Drif We can easily calculae he Malliavin derivaive D P G (l) = P G (l)d ( 1 (l m ) 2 ) 2 st 2 Υ s 2 = P G (l) l m D ( e α Z (TΥ s) dw st 2 Υ s 2 s Z ) 0 = P G (l) l m e α Z (T Υ ) st 2 Υ s 2... by he Malliavin chain rule and properies of Wiener polynomials Thus, using Imkeller s mehod he decomposiion is ξ = W Z = W Z G m e α st 2 s 2 Z (T Υ ) ds Υ 0 0 2α Z e α Z s e 2α Z T Υ e 2α Z s }{{} =a(s) TΥ s e α Z u dw Z u ds

34 Seie 24/35 The Informaion Premium wih Weaher Forecass The Info Premium wih Weaher Forecas Calculaing he Informaion Premium 1/3 The info premium is I G(, T 1, T 2 ; T Υ ) = F G(, T 1, T 2 ; T Υ ) F F(, T 1, T 2 ) Λ(), X, Y erms cancel Lévy erms cancels as well (boh filraions coincide) (16) [ T2 ( u I G(, T 1, T 2 ; T Υ ) = 1 E σρ e α(u s) dw T 2 T 1 s Z T 1 +σ u ) ] 1 ρ 2 e α(u s) dws S du G [ T2 = 1 E σρ T 2 T 1 T 1 u... because expecaion of Iō-inegral is zero e α(u s) dw Z s du G ]

35 Seie 24/35 The Informaion Premium wih Weaher Forecass The Info Premium wih Weaher Forecas Calculaing he Informaion Premium 1/3 The info premium is I G(, T 1, T 2 ; T Υ ) = F G(, T 1, T 2 ; T Υ ) F F(, T 1, T 2 ) Λ(), X, Y erms cancel Lévy erms cancels as well (boh filraions coincide) (16) [ T2 ( u I G(, T 1, T 2 ; T Υ ) = 1 E σρ e α(u s) dw T 2 T 1 s Z T 1 +σ u ) ] 1 ρ 2 e α(u s) dws S du G [ T2 = 1 E σρ T 2 T 1 T 1 u... because expecaion of Iō-inegral is zero e α(u s) dw Z s du G ]

36 Seie 25/35 The Informaion Premium wih Weaher Forecass The Info Premium wih Weaher Forecas Calculaing he Informaion Premium 2/3 We plug in he informaion drif [ T2 E T 1 σρ = E = u e α(u s) dw Z s du G ] [ T2 T2 T 1 T 1 σρ σρ u u ( e α(u s) d ξ s + e α(u s) a(s)e [ TΥ s s 0 TΥ ] a(v) e α Z w dww Z dv )du G v e α Z v dw Z v G ] ds du... an auxiliary resul allows moving he expecaion ou of he double-inegral (i.e. replacing s wih, see App. V) We can hen wrie [ TΥ ] TΥ E e α Z v dwv Z G = E[Z TΥ G ] e α Z (TΥ ) Z α Z µ Z (s)e α Z (TΥ s) ds

37 Seie 25/35 The Informaion Premium wih Weaher Forecass The Info Premium wih Weaher Forecas Calculaing he Informaion Premium 2/3 We plug in he informaion drif [ T2 E T 1 σρ = E = u e α(u s) dw Z s du G ] [ T2 T2 T 1 T 1 σρ σρ u u ( e α(u s) d ξ s + e α(u s) a(s)e [ TΥ s s 0 TΥ ] a(v) e α Z w dww Z dv )du G v e α Z v dw Z v G ] ds du... an auxiliary resul allows moving he expecaion ou of he double-inegral (i.e. replacing s wih, see App. V) We can hen wrie [ TΥ ] TΥ E e α Z v dwv Z G = E[Z TΥ G ] e α Z (TΥ ) Z α Z µ Z (s)e α Z (TΥ s) ds

38 Seie 26/35 The Informaion Premium wih Weaher Forecass The Info Premium wih Weaher Forecas Calculaing he Informaion Premium 3/3 Afer a lenghy calculaion we arrive a ( TΥ ) I G( ) = ρ A( ) T 2 T 1 E[Z TΥ G ] e α Z (TΥ ) Z α Z µ Z (s)e α Z (TΥ s) ds where A( ) is a deerminisic funcion A( ) = 2α Z σe α Z (T Υ+) ( 1 α Z (e α Z (T 2 ) e α Z (T 1 ) ) + 1 α (e α(t 2 ) e α(t 1 ) )) σ Z (α + α Z )(e 2α Z T Υ e 2α Z ) Hence, knowing he curren emperaure and our forecas we can easily calculae he info premium Sign of premium depends on ρ and emperaure erms ρ < 0 Germany, ρ > 0 California, ρ = 0 Qaar

39 Seie 26/35 The Informaion Premium wih Weaher Forecass The Info Premium wih Weaher Forecas Calculaing he Informaion Premium 3/3 Afer a lenghy calculaion we arrive a ( TΥ ) I G( ) = ρ A( ) T 2 T 1 E[Z TΥ G ] e α Z (TΥ ) Z α Z µ Z (s)e α Z (TΥ s) ds where A( ) is a deerminisic funcion A( ) = 2α Z σe α Z (T Υ+) ( 1 α Z (e α Z (T 2 ) e α Z (T 1 ) ) + 1 α (e α(t 2 ) e α(t 1 ) )) σ Z (α + α Z )(e 2α Z T Υ e 2α Z ) Hence, knowing he curren emperaure and our forecas we can easily calculae he info premium Sign of premium depends on ρ and emperaure erms ρ < 0 Germany, ρ > 0 California, ρ = 0 Qaar

40 Seie 27/35 Sylised Examples and Muliple Informaion Sylised Examples and Muliple Informaion

41 Seie 28/35 Sylised Examples and Muliple Informaion Sylised Examples Examples - Seup We find relaed formulae for T Υ T 1 and T 1 < T Υ < T 2 (see App. VI) We wan o check our formula and develop some inuiion Firs, wha is he influence of ime and magniude of weaher forecas? Le s assume we have X = 0 as well as Z = µ Z = 0 Tha means we have zero degrees oday and expec zero degrees Furhermore, le ρ = 0.5, σ = σ Z = and = 0, T 1 = 10 and T 2 = 20 We will le α Z = α and check wo possible values

42 Seie 29/35 Sylised Examples and Muliple Informaion Sylised Examples Info Premium wih α Z = α = Info Premium (Euro) Time T Y Weaher forecas (degrees) Figure: Info Premium for differen values of T Υ and E[Z TΥ G ]

43 Seie 30/35 Sylised Examples and Muliple Informaion Sylised Examples Info Premium wih α Z = α = ,8 0,6 Info Premium (Euro) 0,4 0,2 0 0,2 0,4 04 0,6 0, Time T Y Weaher forecas (degrees) Figure: Info Premium for differen values of T Υ and E[Z TΥ G ]

44 Seie 31/35 Sylised Examples and Muliple Informaion Sylised Examples Realisic Parameers Wha happens if we ry more realisic values? Le s assume again X = 0 as well as Z = µ Z = 0 Again, ρ = 0.5 and = 0, T 1 = 10 and T 2 = 20 Now α = 0.5, σ = and α Z = 0.1 and σ Z = 2.0 Our forecas on he oher hand is E[Z TΥ G ] = 10

45 Seie 32/35 Sylised Examples and Muliple Informaion Sylised Examples Info Premium wih α Z α, σ Z σ 4 Info Premium Delivery Period Value of T Y Figure: Info Premium for differen values of T Υ, E[Z TΥ G ] = 10

46 Seie 33/35 Sylised Examples and Muliple Informaion Muliple Fuure Informaion Muliple Fuure Informaion Normally, we would have a whole series of forecass Especially wih long-delivery conracs his is imporan Playing around wih filraions and spliing inegrals we can incorporae muliple fuure informaion a imes T Υ1, T Υ2,... For his, we have o assume precise (i.e. H = G ) informaion Becomes messy even wih wo pieces

47 Seie 34/35 References and Conac References R. B-K, F. E. Benh, R. Kiesel An empirical sudy of he informaion premium on elecriciy markes, o appear R. B-K, F. E. Benh, R. Kiesel Elecriciy opions and addiional informaion, o appear F. E. Benh and Th. Meyer-Brandis The informaion premium for non-sorable commodiies, Journal of Energy Markes, 2009

48 Seie 35/35 References and Conac Conac phone: +49 (0) Chair for Energy Trading and Finance Universiy Duisburg-Essen web: Thank you for your aenion...

49 Seie 35/35 Appendix I Appendix I Risk-neural valuaion formula yields: [ ] T2 0 = e r E Q e r(u ) (S(u) F(, T 1, T 2 ))du F T 1 If selemens only ake place a he final dae T 2 one ges [ ] T2 0 = e r E Q (S(u) F(, T 1, T 2 ))du F T 1 and finally for he fuures price: [ ] T2 F(, T 1, T 2 ) = E Q 1 S(u)du F T 1 T 2 T 1 17

50 Seie 35/35 Appendix II Appendix II Deerminisic Funcions: { 1 ᾱ(, T 1, T 2 ) = α (e α(t2 ) e α(t1 ) ) T 1 1 α (e α(t2 ) 1) > T 1 { 1 β β(, T 1, T 2 ) = (e β(t2 ) e β(t1 ) ) T 1 1 β (e β(t2 ) 1) > T 1 17 ˆβ(, T 1, T 2 ) = { 1 β (T 2 T β (e β(t 2 ) e β(t 1 ) )) T 1 1 β (T β (e β(t 2 ) 1)) > T 1

51 Seie 35/35 Appendix III Appendix III Remember he Lévy par of he spo For he empirics, we used dy () = βy ()d + dl() N L = where N is Poisson process, inensiy λ, D i are i.i.d jump sizes We used double-exponenially disribued (i.e. he Kou model) wih densiy i=1 D i f D (x) = pη 1 e η 1x 1 x 0 + qη 2 e η 2 x 1 x 0 where p + q = 1 and η 1, η 2 0.

52 Seie 35/35 Appendix IV Appendix IV Theorem Iō s heorem for Lévy processes and addiional incomplee informaion. Le L be a Lévy process and G H = F σ(l TΥ ). Then 1. L is sill a semimaringale wih respec o G 2. if E[ L ] < hen is a G -maringale. TΥ E[L TΥ L s G s ] ξ() = L ds 0 T Υ s 20

53 Seie 35/35 Appendix V Appendix V Theorem Le L() be a Lévy process, F he hisorical filraion. Le G H = F σ(l(t Υ )). If he informaion drif is of he form [ ] TΥ µ G s = g(s)e f (u)dl(u) G s where g and f are con. funcion on [0, T Υ ]. Then [ ] [ TΥ ] TΥ E f (u)dl(u) G = E f (u)dl(u) G s for ime poins s T Υ. 27 s e s f (u)g(u)du

54 Seie 35/35 Appendix VI Appendix VI If T Υ T 1 hen I G(, T 1, T 2 ; T Υ ) = 1 T 2 T 1 ᾱ(t Υ, T 1, T 2 )I G(, T Υ ; T Υ ) where I G (, T Υ ; T Υ ) is he premium wih a delivery ime raher han period If T 1 < T Υ < T 2 hen I G( ) = 1 T 2 T 1 ((T Υ T 1 )I G(, T 1, T Υ ; T Υ ) + (T 2 T Υ )I G(, T Υ, T 2 ; T Υ )) 32

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

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

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

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

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

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

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

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate. Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since

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

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

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

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

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

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

Solutions Problem Set 3 Macro II (14.452)

Solutions Problem Set 3 Macro II (14.452) Soluions Problem Se 3 Macro II (14.452) Francisco A. Gallego 04/27/2005 1 Q heory of invesmen in coninuous ime and no uncerainy Consider he in nie horizon model of a rm facing adjusmen coss o invesmen.

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

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

Section 4.4 Logarithmic Properties

Section 4.4 Logarithmic Properties Secion. Logarihmic Properies 59 Secion. Logarihmic Properies In he previous secion, we derived wo imporan properies of arihms, which allowed us o solve some asic eponenial and arihmic equaions. Properies

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

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

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

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

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

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

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

Section 4.4 Logarithmic Properties

Section 4.4 Logarithmic Properties Secion. Logarihmic Properies 5 Secion. Logarihmic Properies In he previous secion, we derived wo imporan properies of arihms, which allowed us o solve some asic eponenial and arihmic equaions. Properies

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

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

arxiv: v1 [math.pr] 6 Oct 2008

arxiv: v1 [math.pr] 6 Oct 2008 MEASURIN THE NON-STOPPIN TIMENESS OF ENDS OF PREVISIBLE SETS arxiv:8.59v [mah.pr] 6 Oc 8 JU-YI YEN ),) AND MARC YOR 3),4) Absrac. In his paper, we propose several measuremens of he nonsopping imeness of

More information

On the Timing Option in a Futures Contract

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

More information

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

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

Avd. Matematisk statistik

Avd. Matematisk statistik Avd Maemaisk saisik TENTAMEN I SF294 SANNOLIKHETSTEORI/EXAM IN SF294 PROBABILITY THE- ORY WEDNESDAY THE 9 h OF JANUARY 23 2 pm 7 pm Examinaor : Timo Koski, el 79 7 34, email: jkoski@khse Tillåna hjälpmedel

More information

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),

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

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

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

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

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

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

Internet Traffic Modeling for Efficient Network Research Management Prof. Zhili Sun, UniS Zhiyong Liu, CATR

Internet Traffic Modeling for Efficient Network Research Management Prof. Zhili Sun, UniS Zhiyong Liu, CATR Inerne Traffic Modeling for Efficien Nework Research Managemen Prof. Zhili Sun, UniS Zhiyong Liu, CATR UK-China Science Bridge Workshop 13-14 December 2011, London Ouline Inroducion Background Classical

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

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

From Complex Fourier Series to Fourier Transforms

From Complex Fourier Series to Fourier Transforms Topic From Complex Fourier Series o Fourier Transforms. Inroducion In he previous lecure you saw ha complex Fourier Series and is coeciens were dened by as f ( = n= C ne in! where C n = T T = T = f (e

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

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

Backward stochastic dynamics on a filtered probability space

Backward stochastic dynamics on a filtered probability space Backward sochasic dynamics on a filered probabiliy space Gechun Liang Oxford-Man Insiue, Universiy of Oxford based on join work wih Terry Lyons and Zhongmin Qian Page 1 of 15 gliang@oxford-man.ox.ac.uk

More information

Polymer Engineering (MM3POE)

Polymer Engineering (MM3POE) Polymer Engineering (MM3POE) VISCOELASTICITY hp://www.noingham.ac.uk/~eazacl/mm3poe Viscoelasiciy 1 Conens Wha is viscoelasiciy? Fundamenals Creep & creep recovery Sress relaxaion Modelling viscoelasic

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

A Dynamic Model of Economic Fluctuations

A Dynamic Model of Economic Fluctuations CHAPTER 15 A Dynamic Model of Economic Flucuaions Modified for ECON 2204 by Bob Murphy 2016 Worh Publishers, all righs reserved IN THIS CHAPTER, OU WILL LEARN: how o incorporae dynamics ino he AD-AS model

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

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

Lecture 3: Exponential Smoothing

Lecture 3: Exponential Smoothing NATCOR: Forecasing & Predicive Analyics Lecure 3: Exponenial Smoohing John Boylan Lancaser Cenre for Forecasing Deparmen of Managemen Science Mehods and Models Forecasing Mehod A (numerical) procedure

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

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070

More information

Time series Decomposition method

Time series Decomposition method Time series Decomposiion mehod A ime series is described using a mulifacor model such as = f (rend, cyclical, seasonal, error) = f (T, C, S, e) Long- Iner-mediaed Seasonal Irregular erm erm effec, effec,

More information

Linear Dynamic Models

Linear Dynamic Models Linear Dnamic Models and Forecasing Reference aricle: Ineracions beween he muliplier analsis and he principle of acceleraion Ouline. The sae space ssem as an approach o working wih ssems of difference

More information

Licenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A

Licenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A Licenciaura de ADE y Licenciaura conjuna Derecho y ADE Hoja de ejercicios PARTE A 1. Consider he following models Δy = 0.8 + ε (1 + 0.8L) Δ 1 y = ε where ε and ε are independen whie noise processes. In

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

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

Physics for Scientists & Engineers 2

Physics for Scientists & Engineers 2 Direc Curren Physics for Scieniss & Engineers 2 Spring Semeser 2005 Lecure 16 This week we will sudy charges in moion Elecric charge moving from one region o anoher is called elecric curren Curren is all

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

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

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

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

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

Constant Acceleration

Constant Acceleration Objecive Consan Acceleraion To deermine he acceleraion of objecs moving along a sraigh line wih consan acceleraion. Inroducion The posiion y of a paricle moving along a sraigh line wih a consan acceleraion

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

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n Lecure 3 - Kövari-Sós-Turán Theorem Jacques Versraëe jacques@ucsd.edu We jus finished he Erdős-Sone Theorem, and ex(n, F ) ( /(χ(f ) )) ( n 2). So we have asympoics when χ(f ) 3 bu no when χ(f ) = 2 i.e.

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

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product 11.1 APPCATON OF AMPEE S AW N SYMMETC MAGNETC FEDS - f one knows ha a magneic field has a symmery, one may calculae he magniude of by use of Ampere s law: The inegral of scalar produc Closed _ pah * d

More information

Technical Appendix. Political Uncertainty and Risk Premia

Technical Appendix. Political Uncertainty and Risk Premia echnical Appendix o accompany Poliical Uncerainy and Risk Premia Ľuboš Pásor Universiy of Chicago, CEPR, andber Piero Veronesi Universiy of Chicago, CEPR, andber Sepember 4, 011 Proof of Lemma 1: he same

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

( ) ( ) 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

Generalized Snell envelope and BSDE With Two general Reflecting Barriers

Generalized Snell envelope and BSDE With Two general Reflecting Barriers 1/22 Generalized Snell envelope and BSDE Wih Two general Reflecing Barriers EL HASSAN ESSAKY Cadi ayyad Universiy Poly-disciplinary Faculy Safi Work in progress wih : M. Hassani and Y. Ouknine Iasi, July

More information

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

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

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

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese

More information

Homework 4 (Stats 620, Winter 2017) Due Tuesday Feb 14, in class Questions are derived from problems in Stochastic Processes by S. Ross.

Homework 4 (Stats 620, Winter 2017) Due Tuesday Feb 14, in class Questions are derived from problems in Stochastic Processes by S. Ross. Homework 4 (Sas 62, Winer 217) Due Tuesday Feb 14, in class Quesions are derived from problems in Sochasic Processes by S. Ross. 1. Le A() and Y () denoe respecively he age and excess a. Find: (a) P{Y

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

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

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

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

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

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

Georey E. Hinton. University oftoronto.   Technical Report CRG-TR February 22, Abstract Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical

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

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

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

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

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

Simulating models with heterogeneous agents

Simulating models with heterogeneous agents Simulaing models wih heerogeneous agens Wouer J. Den Haan London School of Economics c by Wouer J. Den Haan Individual agen Subjec o employmen shocks (ε i, {0, 1}) Incomplee markes only way o save is hrough

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

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

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

More information

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

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

More information

Affine term structure models

Affine term structure models Affine erm srucure models A. Inro o Gaussian affine erm srucure models B. Esimaion by minimum chi square (Hamilon and Wu) C. Esimaion by OLS (Adrian, Moench, and Crump) D. Dynamic Nelson-Siegel model (Chrisensen,

More information