Actuarial Models 1: solutions example sheet 4

Similar documents
EXERCISES FOR SECTION 6.3

ERTH403/HYD503, NM Tech Fall 2006

Introduction to Mechanism Design

Correction for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002

Bogoliubov Transformation in Classical Mechanics

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Adelic Modular Forms

1 Routh Array: 15 points

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281

IEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is

Table of Common Derivatives By David Abraham

Problem Set 8 Solutions

Chapter 4. The Laplace Transform Method

Integration Review. May 11, 2013

CHE302 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang

Math 273 Solutions to Review Problems for Exam 1

Chapter 3 : Transfer Functions Block Diagrams Signal Flow Graphs

Problem set 2: Solutions Math 207B, Winter 2016

TAYLOR POLYNOMIALS FOR NABLA DYNAMIC EQUATIONS ON TIME SCALES

Nuclear and Particle Physics - Lecture 16 Neutral kaon decays and oscillations

Physics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014

TMA4125 Matematikk 4N Spring 2016

TMA4195 Mathematical modelling Autumn 2012

Diagonalization of Matrices Dr. E. Jacobs

Practice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions

MATEMATIK Datum: Tid: eftermiddag. A.Heintz Telefonvakt: Anders Martinsson Tel.:

CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang

Math 342 Partial Differential Equations «Viktor Grigoryan

MA 266 FINAL EXAM INSTRUCTIONS May 2, 2005

Lectures on Markov models in life and sickness insurance.

March 18, 2014 Academic Year 2013/14

Solving Differential Equations by the Laplace Transform and by Numerical Methods

Semilinear obstacle problem with measure data and generalized reflected BSDE

Calculus Class Notes for the Combined Calculus and Physics Course Semester I

RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS

Physics 443 Homework 2 Solutions

Problem 1. Construct a filtered probability space on which a Brownian motion W and an adapted process X are defined and such that

Math 1271 Solutions for Fall 2005 Final Exam

ON ITERATIVE FEEDBACK TUNING AND DISTURBANCE REJECTION USING SIMPLE NOISE MODELS. Bo Wahlberg

Laplace Transformation

You should also review L Hôpital s Rule, section 3.6; follow the homework link above for exercises.

Solution 3.1 Prove the following: γ d. (a) Start with fundamental definitions: V = (b) e = 1 n. wg e S =

Nonlinear Single-Particle Dynamics in High Energy Accelerators

Mathematical Review Problems

Introduction to the Vlasov-Poisson system

Singular perturbation theory

Reading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

These are practice problems for the final exam. You should attempt all of them, but turn in only the even-numbered problems!

Performance Evaluation of Acoustic Scene Classification Using DNN-GMM and Frame-Concatenated Acoustic Features

Section 7.1: Integration by Parts

into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get

LECTURE 12: LAPLACE TRANSFORM

Lecture 15 - Current. A Puzzle... Advanced Section: Image Charge for Spheres. Image Charge for a Grounded Spherical Shell

The machines in the exercise work as follows:

Fermi Distribution Function. n(e) T = 0 T > 0 E F

TRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL

The statistical properties of the primordial fluctuations

SOME RESULTS ON INFINITE POWER TOWERS

Markov Chains in Continuous Time

Lecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004

EconS Advanced Microeconomics II Handout on Bayesian Nash Equilibrium

Compensation of backlash effects in an Electrical Actuator

Solutions to MATH 271 Test #3H

THE BICYCLE RACE ALBERT SCHUELLER

Optimal scheduling in call centers with a callback option

A note on strong convergence for Pettis integrable functions (revision)

ME 375 FINAL EXAM Wednesday, May 6, 2009

Framework Model For Single Proton Conduction through Gramicidin

V = 4 3 πr3. d dt V = d ( 4 dv dt. = 4 3 π d dt r3 dv π 3r2 dv. dt = 4πr 2 dr

Chapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem

The Laplace Transform , Haynes Miller and Jeremy Orloff

ME 375 EXAM #1 Tuesday February 21, 2006

Kinematics of Rotations: A Summary

Chapter 7: The Laplace Transform Part 1

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask

Lecture 7: Testing Distributions

Schrödinger s equation.

2.0 ANALYTICAL MODELS OF THERMAL EXCHANGES IN THE PYRANOMETER

arxiv: v1 [math-ph] 5 May 2014

arxiv: v2 [math.nt] 30 Apr 2015

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations

Final Exam Study Guide and Practice Problems Solutions

More from Lesson 6 The Limit Definition of the Derivative and Rules for Finding Derivatives.

Hybrid Control and Switched Systems. Lecture #6 Reachability

Calculus of Variations

Lecture 10: Forward and Backward equations for SDEs

PHASE-FIELD SIMULATION OF SOLIDIFICATION WITH DENSITY CHANGE

Convex Hulls of Curves Sam Burton

Many problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule

Linear First-Order Equations

Introduction to Laplace Transform Techniques in Circuit Analysis

05 The Continuum Limit and the Wave Equation

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam

Reading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions

New bounds for Morse clusters

Calculus of Variations

DIFFERENTIAL EQUATIONS

Transcription:

Actuarial Moel 1: olution example heet 4 (a) Anwer to Exercie 4.1 Q ( e u ) e σ σ u η η. (b) The forwar equation correponing to the backwar tate e are t p ee(t) σp ee (t) + ηp eu (t) t p eu(t) σp ee (t) ηp eu (t), with bounary conition p ee (0) 1 an p eu (0) 0. (c) Working with the firt equation an replacing p eu (t) by 1 p ee (t) (a p ee (t)+p eu (t) 1) we get t p ee(t) (σ + η)p ee (t) + η. T i a firt orer, linear, non-homogeneou ODE that we nee to olve. There are everal way to o t, but here we follow the metho of firt etermining the general olution of the homogeneou verion of the ODE an then fining a particular olution to the non-homogeneou ODE 1. Let firt olve the homogeneou part of above ODE, i.e. (t) (σ + η)p hom (t). We ee that t phom ee ee t log ( p hom ee (t) ) t phom ee p hom ee (t) (σ + η) (t) an o by integrating on both ie we ee that the general olution to the homogeneou part i given by p hom ee (t) Ce (σ+η)t for ome contant C. Secon, we try to fin a particular olution to the non-homogeneou ODE (note that for uch a particular olution the bounary conition oe not have to be atifie). For t we ue the metho of unetermine coefficient. A the non-homogeneou part of the ODE i η wch i contant in t, we gue that a particular olution might be alo a contant in t. So we try p part ee (t) A. Then by the ODE we mut have 0 (σ + η)a + η, i.e. A η. Hence ppart σ+η ee (t) η i a particular olution to the non-homogeneou ODE. σ+η Hence (by a theorem from the theory of ODE) the general olution to the original ODE i given by p ee (t) p hom ee (t) + p part ee (t) Ce (σ+η)t + η σ + η. 1 Alternative metho are uing an integrating factor or uing the general formula (3). 1

Invoking the bounary conition p ee (0) 1 lea to C σ. Hence in the en we σ+η have p ee (t) σ σ + η e (σ+η)t + η σ + η, p eu (t) 1 p ee (t) σ ( ) 1 e (σ+η)t. σ + η Remark In general for the non-homogeneou, firt orer, linear ODE, f (t) g(t)f(t) + h(t), (1) where g an h are given function (that are (locally) boune an continuou) an with bounary conition f(t 0 ) a, (2) with t 0, a R given, we have that it olution i given by ( t t ) t g(u)u f(t) e 0 h()e t g(u)u 0 + a. (3) t 0 Inee, it i eay to check, via the prouct rule of ifferentation an the funamental theorem of calculu, that (3) i a - an thu the unique - olution to the ODE (1) t with bounary conition (2). Note that Ce g(u)u 0 with C a contant i the general olution of the homogeneou part of the ODE (i.e. f (t) g(t)f(t)), wherea t g(u)u e 0 t 0 h()e t g(u)u 0 i a particular olution to the ODE (1) (but oe not atify the bounary conition (2) unle a 0). Anwer to Exercie 4.2 The Q-matrix, wch epen on age t, i given by (µ(t) + ν(t)) µ(t) ν(t) Q 0 0 0. 0 0 0 So for each 0, the Kolmogorov forwar equation for the tranition probabilitie with backwar tate 1 are given by t p 11(, t) (µ(t) + ν(t))p 11 (, t), t p 12(, t) µ(t)p 11 (, t), t p 13(, t) ν(t)p 11 (, t), where t > with bounary conition p 11 (, ) 1, p 12 (, ) 0 an p 13 (, ) 0. The olution of the firt equation i eaily een to be given by p 11 (, t) e (µ(u)+ν(u))u. (4) 2

In orer to ee t, note that by the firt ODE we have for all u >, u log (p 11(, u)) p u 11(, u) p 11 (, u) (µ(u) + ν(u)) an thu integrating both ie with repect to u from to t an uing the bounary conition, log(p 11 (, t)) log(p 11 (, t)) log(p 11 (, )) u log (p 11(, u)) u (µ(u)+ν(u))u. Hence (4) follow. For p 12 (, t), we have by uing repectively the bounary conition p 12 (, ) 0, the funamental theorem of calculu, the econ ODE an (4), p 12 (, t) p 12 (, t) p 12 (, ) u p 12(, u)u µ(u)p 11 (, u)u µ(u)e u (µ(v)+ν(v))v u. Similarly, we get p 12 (, t) ν(u)e u (µ(v)+ν(v))v u. Anwer to Exercie 4.3 The Q-matrix for t moel i given by (λ + µ) λ µ Q 0 ν ν, 0 0 0 whereby h i the firt tate, i i the econ an the tr. (a) The 3-imenional ytem of Kolmogorov forwar equation correponing to the backwar tate i are: t p ih(t) (λ + µ)p ih (t) t p ii(t) λp ih (t) νp ii (t) t p i(t) µp ih (t) + νp ii (t). We ee that the firt ODE involve only p ih (t). Hence we can olve it by eparation of variable to get p ih (t) Ae (λ+µ)t, t 0. 3

Since p ih (0) 0, we mut have A 0 an o p ih (t) 0 for all t 0. Now the econ ODE become t p ii(t) νp ii (t) an ince p ii (0) 1, we euce p ii (t) e νt. We can fill t in into the tr ODE an olve for p i (t), or intea we ue the fact that p ih (t)+p ii (t)+p i (t) 1 to euce that p i (t) 1 p ii (t) 1 e νt. (b) The 3-imenional ytem of Kolmogorov forwar equation correponing to the backwar tate h are: t p hh(t) (λ + µ)p hh (t) t p (t) λp hh (t) νp (t) t p h(t) µp hh (t) + νp (t). Note that they are baically the ame a in part (a); only the backwar tate change from i to h. However, the olution will be ifferent becaue the bounary conition are not the ame. The firt ODE together with the bounary conition p hh (0) 1 give p hh (t) e (λ+µ)t, t 0. The econ ODE now become t p (t) λe (λ+µ)t νp (t). The olution to t ODE i of the form p (t) p hom (t) + p part (t), where p hom (t) i the general olution to the homogeneou problem: t phom (t) νp hom (t) an p part (t) i a particular olution to the original ODE. We can eaily euce that p hom (t) Ae νt, where A i a contant. Further, a the non-homogeneou part i equal to e (λ+µ)t it i a goo iea to tart looking for a particular olution of the form p part (t) Be (λ+µ)t, with B another contant uch that the original ODE i atifie. We eaily euce that the contant B mut be uch that B(λ + µ) λ νb. T equation only ha a olution if ν λ + µ in wch cae B λ/(λ + µ ν). Hence if ν λ + µ, then p (t) Ae νt 4 λ λ + µ ν e (λ+µ)t.

The contant A can now be euce from p (0) 0. T give A λ. λ+µ ν If ν λ + µ we have to fin a particular olution of a ifferent form. The next bet gue i to try p part (t) Bte νt. We ee that p part (t) atifie the ODE if B λ. In t cae the contant A will then be equal to zero. Hence we have { λ ( λ+µ ν e νt e (λ+µ)t) if ν λ + µ, p (t) λte νt if ν λ + µ. Note that by l Hôpital rule lim (λ+µ) ν λ ( e νt e (λ+µ)t) λte νt λ + µ ν an t coul alo have been ue to gue the olution of the ODE in the cae where ν λ + µ. (Alternatively, one can ue the general formula (3).) Uing p hh (t) + p (t) + p h (t) 1, we euce that { ( 1 1 λ+µ ν λe νt + (µ ν)e (λ+µ)t) if ν λ + µ, p h (t) 1 (1 + λt)e νt if ν λ + µ. Anwer to Exercie 4.4 We nee to prove p ik (, t) p ij(, u)p jk (u, t) for all 0 < u < t an i, k S. We have p ik (, t) Pr(X t k X i) Pr(X t k, X u j X i) Pr(X t k X u j, X i)pr(x u j X i) Pr(X t k X u j)pr(x u j X i) p ij (, u)p jk (u, t), where we have ue the aitivity property of (conitional) probability meaure in the econ equality (note that j S {X t k, X u j} {X t k} an that {X t k, X u j} {X t k, X u i} for i j), the efinition of conitional probability in the tr equality an the Markov property in the fourth equality. 5

We have that Hence Anwer to Exercie 4.5 p il (, t) 1 ( t p il(, t) 1 t j1 p ij (, t) j1 p ik (, t). p ik (, t) t p ik(, t) ) p ij (, t)µ jk (t) p ij (, t)µ jl (t), j1 µ jk (t) where we have ue that ( ) i atifie for k l in the tr equality an the fact that the rowum of Q(t) are 0 (becaue Q(t) i a Q-matrix) in the lat equality. Hence ( ) i atifie for k l a well. 6