Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
литература / Digital_Video_and_HD_Second_Edition_Algorithms_and_Interfaces.pdf
Скачиваний:
0
Добавлен:
13.05.2026
Размер:
38.02 Mб
Скачать

Strictly speaking, amplitude is an instantaneous measure that may take a positive or negative value. Magnitude is properly either an absolute value, or a squared or root mean square (RMS) value representative of amplitude over some time interval. The terms are often used interchangeably.

See Linearity on page 37.

Bracewell, Ronald N. (1985),

The Fourier Transform and Its

Applications, Second Edition (New

York: McGraw-Hill).

Magnitude frequency response

To gain a general appreciation of aliasing, it is necessary to understand signals in the frequency domain. The previous section gave an example of inadequate filtering prior to sampling that created an unexpected alias upon sampling. You can determine whether a filter has an unexpected response at any frequency by presenting to the filter a signal that sweeps through all frequencies, from zero, through low frequencies, to some high frequency, plotting the response of the filter as you go. I graphed such a frequency sweep signal at the top of Figure 9.1, on page 98. The middle graph of that figure shows a response waveform typical of a lowpass filter (LPF), which attenuates high frequency signals. The magnitude response of that filter is shown in the bottom graph.

Magnitude response is the RMS average response over all phases of the input signal at each frequency. As you saw in the previous section, a filter’s response can be strongly influenced by the phase of the input signal. To determine response at a particular frequency, you can test all phases at that frequency. Alternatively, provided the filter is linear, you can present just two signals – a cosine wave at the test frequency and a sine wave at the same frequency. The filter’s magnitude response at any frequency is the absolute value of the vector sum of the responses to the sine and the cosine waves.

Analytic and numerical procedures called transforms can be used to determine frequency response. The Laplace transform is appropriate for continuous functions, such as signals in the analog domain. The Fourier transform is appropriate for signals that are sampled periodically, or for signals that are themselves periodic. A variant intended for computation on data that has been sampled is the discrete Fourier transform (DFT). An elegant scheme for numerical computation of the DFT is the fast Fourier transform (FFT). The z-transform is essentially a generalization of the Fourier transform. All of these transforms represent mathematical ways to determine a system’s response to sine waves over a range of frequencies and phases. The result of a transform is an expression or graph in terms of frequency.

196

DIGITAL VIDEO AND HD ALGORITHMS AND INTERFACES

1.0

 

 

 

 

 

 

 

0.8

 

 

 

 

 

 

 

0.6

 

 

 

 

 

 

 

0.4

 

 

 

 

 

 

 

0.2

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

-1

0

1

2

3

4

5

6

1.0

 

 

 

Time, t [multiple of ts]

 

0.8

 

 

 

 

 

 

 

0.6

 

 

 

 

 

 

 

0.4

 

 

 

 

 

 

 

0.2

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

-0.2

0

10π

12π

 

 

 

Frequency, ω [rad·s-1]

 

 

 

 

 

 

Figure 20.7 Frequency response of a boxcar filter. The top graph shows a boxcar weighting function, symmetrical around t =0. Its frequency spectrum is a sinc function, shown underneath. The solid line shows that at certain frequencies, the filter causes phase inversion. Filter response is usually plotted as magnitude; phase inversion in the stopband is reflected as the absolute (magnitude) values shown in dashed lines.

1,

 

 

ω = 0

 

 

 

 

sinc ω = sin πω

ω 0

 

 

 

,

 

 

 

πω

 

Eq 20.1 The sinc function

(pronounced sink) is defined by this equation. Its argument is in radians per second (rad·s-1); here I use the conventional symbol ω for that quantity. The term (sin x)/x (pronounced sine ecks over ecks) is often used synonymously with sinc, without mention of the units of the argument. If applied to frequency in hertz, the function could be written (sin 2πf)/f.

sinc is unrelated to sync (synchronization).

Magnitude frequency response of a boxcar

The top graph of Figure 20.7 shows the weighting function of point sampling, as a function of time (in sample intervals). The Fourier transform of the boxcar function – that is, the magnitude frequency response of a boxcar weighting function – takes the shape of

(sin x)/x. The response is graphed at the bottom of Figure 20.7, with the frequency axis in units of

ω =2πfS. Equation 20.1 in the margin defines the function. This function is so important that it has been given the special symbol sinc, introduced by Phillip M. Woodward in 1952 as a contraction of sinus cardinalis.

A presampling filter should have fairly uniform response below half the sample rate, to provide good sharpness, and needs to severely attenuate frequencies at and above half the sample rate, to achieve low aliasing. The bottom graph of Figure 20.7 shows that this requirement is not met by a boxcar weighting function. The graph of sinc predicts frequencies where aliasing can be introduced. Figure 20.6 showed an example of a sine wave at 0.75fS; reading the value of

CHAPTER 20

FILTERING AND SAMPLING

197

A near-ideal filter in analog video is sometimes called a brick wall filter, though there is no precise definition of this term.

sinc at 1.5π from Figure 20.7 shows that aliasing is expected.

You can gain an intuitive understanding of the boxcar weighting function by considering that when the input frequency is such that an integer number of cycles lies under the boxcar, the response will be null. But when an integer number of cycles, plus a half-cycle, lies under the weighting function, the response will exhibit a local maximum that can admit an alias.

To obtain a presampling filter that rejects potential aliases, we need to pass low frequencies, up to almost half the sample rate, and reject frequencies above it.

We need a frequency response that is constant at unity up to just below 0.5fS, whereupon it drops to zero. We need a filter function whose frequency response – not time response – resembles a boxcar.

The sinc weighting function

Remarkably, the Fourier transform possesses the mathematical property of being its own inverse (within a scale factor). In Figure 20.7, the Fourier transform of a boxcar weighting function produced a sinc-shaped frequency response. Figure 20.8 opposite shows a sinc-shaped weighting function; it produces a boxcar-shaped frequency response. So, sinc weighting gives the ideal lowpass filter (ILPF), and sinc is the ideal temporal weighting function for use in a presampling filter. However, there are several theoretical and practical difficulties in using sinc. In practice, we approximate it.

An analog filter’s response is a function of frequency on the positive real axis. In analog signal theory, there is no upper bound on frequency. But in a digital filter the response to a test frequency fT is identical to the response at fT offset by any integer multiple of the sampling frequency: The frequency axis “wraps” at multiples of the sampling rate. Sampling theory also dictates “folding” around half the sample rate. Signal components having frequencies at or above the Nyquist rate cannot accurately be represented.

The temporal weighting functions used in video are usually symmetrical; nonetheless, they are usually graphed in a two-sided fashion. The frequency response of a filter suitable for real signals is symmetrical about

198

DIGITAL VIDEO AND HD ALGORITHMS AND INTERFACES

1.0

0.8

0.6

0.4

0.2

-5

-4

-3

-2

-1 -0.2

1

2

3

4

5

1.0Time, t [multiple of ts]

0.8

0.6

0.4

0.2

-5

-4

-3

-2

-1 -0.2

10π

 

 

 

 

 

 

Frequency, ω [rad·s-1]

 

Figure 20.8 The sinc weighting function is shown in the top graph. Its frequency spectrum, shown underneath, has

a boxcar shape: sinc weighting exhibits the ideal properties for a presampling filter. However, its infinite extent makes it physically unrealizable; also, its negative lobes make it unrealizable for transducers of light such as cameras, scanners, and displays. Many practical digital lowpass filters have coefficients that approximate samples of sinc.

zero; conventionally, frequency response is graphed in one-sided fashion starting at zero frequency (“DC”). Sometimes it is useful to consider or graph frequency response in two-sided style.

Frequency response of point sampling

The Fourier transform provides an analytical tool to examine frequency response: We can reexamine point sampling. Taking an instantaneous sample of a waveform is mathematically equivalent to using a weighting function that is unity at the sample instant, and zero everywhere else – the weighting function is an impulse. The Fourier transform of an impulse function is constant, unity, at all frequencies. A set of equally spaced impulses is an impulse train; its transform is also unity everywhere. The sampling operation is represented as multiplication by an impulse train. An unfiltered signal sampled by a set of impulses will admit aliases equally from all input frequencies.

CHAPTER 20

FILTERING AND SAMPLING

199

Соседние файлы в папке литература