28.02.2013 Views

Introduction to Acoustics

Introduction to Acoustics

Introduction to Acoustics

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.

516 Part D Hearing and Signal Processing<br />

Part D 14.7<br />

which, for the one-pole filter, is Φ(ω) = tan −1 (−ωτ).<br />

The phase shift is zero at zero frequency, and approaches<br />

90 ◦ at high frequency. This phase behavior is typical<br />

of simple filters in that important phase shifts occur<br />

in frequency regions where the magnitude shows large<br />

attenuation.<br />

14.7.2 Phase Delay and Group Delay<br />

The phase shifts introduced by filters can be interpreted<br />

as delays, whereby the output is delayed in time<br />

compared <strong>to</strong> the input. In general, the delay is different<br />

for different frequencies, and therefore, a complex<br />

signal composed of several frequencies is bent out of<br />

shape by the filtering process. Systems in which the delay<br />

is different for different frequencies are said <strong>to</strong> be<br />

dispersive.<br />

Two kinds of delay are of interest. The phase delay<br />

simply reinterprets the phase shift as a delay. The phase<br />

delay Tϕ is given by Tϕ =−Φ(ω)/ω.Thegroup delay Tg<br />

is given by the derivative Tg =−dΦ(ω)/dω.Phaseand<br />

group delays for a one-pole low-pass filter are shown in<br />

Fig. 14.7 <strong>to</strong>gether with the phase shift.<br />

14.7.3 Resonant Filters<br />

Resonant filters, or tuned systems, have an amplitude response<br />

that has a peak at some frequency where ω = ωo.<br />

Such filters are characterized by the resonant frequency,<br />

ωo, and by the bandwidth, 2∆ω. The bandwidth<br />

is specified by half-power points such that |H(ωo +<br />

∆ω)| 2 ≈|H(ωo)| 2 /2and|H(ωo −∆ω)| 2 ≈|H(ωo)| 2 /2.<br />

The sharpness of a tuned system is often quoted as<br />

a Q value, where Q is a dimensionless number given<br />

Phase shift (rad)<br />

2<br />

1<br />

0<br />

0<br />

Tg<br />

Tj<br />

2 4 6 8 10<br />

ùô<br />

Fig. 14.7 The phase shift Φ for a one-pole low-pass filter<br />

can be read on the left ordinate. The phase and group delays<br />

can be read on the right ordinate<br />

j<br />

2ô<br />

ô<br />

by<br />

Q = ωo/(2∆ω) . (14.105)<br />

As an example, a two-pole low-pass filter with a resonant<br />

peak near the angular frequency ωo is described by the<br />

transfer function<br />

H(ω) =<br />

ω2 o<br />

ω2 o − ω2 . (14.106)<br />

+ jωωo/Q<br />

14.7.4 Impulse Response<br />

Because filtering is described as a product of Fourier<br />

transforms, i. e., in frequency space, the temporal representation<br />

of filtering is a convolution<br />

�<br />

y(t) = dt ′ h(t − t ′ )x(t ′ �<br />

) = dt ′ h(t ′ )x(t − t ′ ) .<br />

(14.107)<br />

The two integrals on the right are equivalent.<br />

Equation (14.107) is a special form of linear processor.<br />

A more general linear processor is described by the<br />

equation<br />

�<br />

y(t) = dt ′ h(t, t ′ )x(t ′ ) , (14.108)<br />

where h(t, t ′ ) permits a perfectly general dependence on<br />

t and t ′ . The special system in which only the difference<br />

in time values is important, i. e. h(t, t ′ ) = h(t − t ′ ), is<br />

a linear time-invariant system. Filters are time invariant.<br />

A system that operates in real time obeys a further<br />

filter condition, namely causality. A system is causal if<br />

the output y(t) depends on the input x(t ′ ) only for t ′ < t.<br />

In words, this says that the present output cannot depend<br />

on the future input. Causality requires that h(t) = 0for<br />

t < 0. For the one-pole corona, low-pass filter of (14.101)<br />

the impulse response is<br />

h(t) = 1<br />

for t > 0 ,<br />

τ e−t/τ<br />

h(t) = 0 fort < 0 ,<br />

h(t) = 1<br />

2τ<br />

for t = 0 . (14.109)<br />

For the two-pole low-pass resonant filter of (14.106),<br />

the impulse response is<br />

ωo<br />

h(t) = �<br />

1 −[1/(2Q)] 2<br />

� �<br />

×sin ωot 1 −[1/(2Q)] 2<br />

�<br />

, t ≥ 0 ,<br />

e− ωo<br />

2Q t<br />

h(t) = 0 , t < 0 . (14.110)

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

Saved successfully!

Ooh no, something went wrong!