Second-order ARMA filter. Using Example 16.41 as a guide, compose Matlab code to implement...

Second-order ARMA filter. Using Example 16.41 as a guide, compose Matlab code to implementan ARMA sys-tem that has

hi = sin(ωoi) for 0 ≤ i ≤ 64

for ωo = π/8. Hint: the difference equation is

yi = 2cos(ωo)yi−1 − yi−2 + sin(ωo)xi−1

Using subplot(3,1,n),stem(.), plot hi computed directly from the sine function, yi, and the error hi − yi.

Example 16.41

ARMA digital filter

Figure 16.39 shows the Matlab function that implements an ARMA digital filter by cascading an MA filter followed by an AR filter. The input xi to the MA filter produces output y1i , which then acts as the input to the AR filter to produce output yi.

Figure 16.39

Matlab script that implements an ARMA digital filter.