04.01.2013 Views

Springer - Read

Springer - Read

Springer - Read

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

4.4 The Spectral Density of an ARMA Process 133<br />

Since {Zt} has spectral density fZ(λ) σ 2 /(2π), the result now follows from Proposition<br />

4.3.1.<br />

For any specified values of the parameters φ1,...,φp,θ1,...,θq and σ 2 , the<br />

Spectrum>Model option of ITSM can be used to plot the model spectral density.<br />

Example 4.4.1 The spectral density of an AR(2) process<br />

For an AR(2) process (4.4.1) becomes<br />

fX(λ) <br />

<br />

σ 2<br />

2π 1 − φ1e −iλ − φ2e −2iλ 1 − φ1e iλ − φ2e 2iλ<br />

σ 2<br />

2π 1 + φ 2 1 + 2φ2 + φ 2 2 + 2(φ1φ2 − φ1) cos λ − 4φ2 cos 2 λ .<br />

Figure 4.14 shows the spectral density, found from the Spectrum>Model option of<br />

ITSM, for the model (3.2.20) fitted to the mean-corrected sunspot series. Notice the<br />

well-defined peak in the model spectral density. The frequency at which this peak<br />

occurs can be found by differentiating the denominator of the spectral density with<br />

respect to cos λ and setting the derivative equal to zero. This gives<br />

cos λ φ1φ2 − φ1<br />

4φ2<br />

0.849.<br />

The corresponding frequency is λ 0.556 radians per year, or equivalently<br />

c λ/(2π) 0.0885 cycles per year, and the corresponding period is therefore<br />

1/0.0885 11.3 years. The model thus reflects the approximate cyclic behavior of<br />

f<br />

0 200 400 600 800<br />

Figure 4-14<br />

The spectral density<br />

fX (λ), 0 ≤ λ ≤ π of the<br />

AR(2) model (3.2.20) fitted<br />

to the mean-corrected<br />

0.0 0.5 1.0 1.5 2.0 2.5 3.0<br />

sunspot series. Frequency

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!