Simulate AR(1) Poisson Jump Process

조회 수: 5 (최근 30일)
Julian Williams
Julian Williams 2014년 2월 20일
I am being a bit thick, so apologies if this is a trivial question.
I would like to simulate a mixed frequency (intra-day versus daily) discrete time AR(1) jump process.
So the data generating process of interest is:
N(t) is the counting process.
with intensity eta(t)
I would like the intensity to have the following simple dynamics:
eta(t+1) = rho*eta(t) + gamma*noise!
Now the problem I have is the noise! bit.
when writing out the density function I follow the textbook treatment and write:
E(Nt) - eta(t) = sum j*Probability(Nt = j) - eta(t)
Probability(N(t+1) = j) = (1/factorial(j))*exp(eta(t+1))*(eta(t+1))^j;
So I use the law of iterated expectations to recover noise!
Sad thing is that for a Hawks process I know how to generate this, but I just seem to be clueless how to do it for an independent noise process.
If anyone can help, it would be much appreciated, I am sure it is something trivial and I am just missing something practical.
Thanks.

답변 (0개)

카테고리

Help CenterFile Exchange에서 Communications Toolbox에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by