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

Figure 20.26 A 25-tap lowpass FIR filter

g[i] = 0.098460si−12 +0.009482si−11 −0.013681si−10 +0.020420si−9

−0.029197si−8

+0.039309si−7

−0.050479si−6

+0.061500si−5

−0.071781si− 4

+0.080612si−3

−0.087404si−2

+0.091742si−1

+0.906788si

+0.091742si+1

−0.087404si+2

+0.080612si+3

−0.071781si+ 4

+0.061500si+5

−0.050479si+6

+0.039309si+7

−0.029197si+8

+0.020420si+9 −0.013681si+10 +0.009482si+11 +0.098460si+12

no control can be exercised over the width of the transition band. The Kaiser window has a single parameter that controls that width. For a given filter order, if the transition band is made narrower, then stopband attenuation is reduced. The Kaiser window parameter allows the designer to determine this tradeoff.

A windowed sinc filter has much better performance than a truncated sinc, and windowed design is so simple that there is no excuse to use sinc without windowing. In most engineering applications, however, filter performance is best characterized in the frequency domain, and the frequency-domain performance of windowed sinc filters is suboptimal: The performance of an n-tap windowed sinc filter can be bettered by an n-tap filter whose design has been suitably optimized.

Few closed-form methods are known to design optimum digital filters. Design of a high-performance filter usually involves successive approximation, optimizing by trading design parameters back and forth between the time and frequency domains. The classic method was published by J.H. McLellan, T.W. Parks, and L.R. Rabiner (“MPR”), based upon an algorithm developed by the Russian mathematician E.Ya. Remez. In the DSP community, the method is often called the “Remez exchange.”

The coefficients of a high-quality lowpass filter for studio video are shown in Figure 20.26 in the margin.

Reconstruction

Digitization involves sampling and quantization; these operations are performed in an analog-to-digital converter (ADC). Whether the signal is quantized then sampled, or sampled then quantized, is relevant only within the ADC: The order of operations is immaterial outside that subsystem. Modern video ADCs quantize first, then sample.

I have explained that filtering is generally required prior to sampling in order to avoid the introduction of aliases. Avoidance of aliasing in the sampled domain has obvious importance. In order to avoid aliasing, an analog presampling filter needs to operate prior to analog-to-digital conversion. If aliasing is avoided, then the sampled signal can, according to Shannon’s theorem, be reconstructed without aliases.

216

DIGITAL VIDEO AND HD ALGORITHMS AND INTERFACES

PRESAMPLING

SAMPLED

POSTSAMPLING

(ANTIALIASING)

DOMAIN

(RECONSTRUCTION)

FILTER

 

FILTER

 

 

A

D

D

A

Figure 20.27 Sampling and reconstruction

To reconstruct an analog signal, an analog reconstruction filter is necessary following digital-to-analog (D-to-A) conversion. The overall flow is sketched in Figure 20.27 above.

Reconstruction close to 0.5fS

Consider the example in Figure 20.28 below of a sine wave at 0.44fS. This signal meets the sampling criterion, and can be perfectly represented in the digital domain. However, from an intuitive point of view, it is difficult to predict the underlying sinewave from samples 3, 4, 5, and 6 in the lower graph. When reconstructed using

a Gaussian filter, the high-frequency signal vanishes. To be reconstructed accurately, a waveform with a significant amount of power near half the sampling rate must be reconstructed with a high-quality filter.

 

1+ sin0.44

1

 

t

 

 

 

 

2

0.5

 

 

 

 

0

 

 

1

 

 

0.5

Figure 20.28 Reconstruction

 

 

close to 0.5fS

0

2

3

4

5

6

7

CHAPTER 20

FILTERING AND SAMPLING

217

1+ sin0.44 t

2

Figure 20.29 D-to-A conversion with a boxcar waveform is equivalent to a DAC producing an impulse train followed by a boxcar filter with its sinc response. Frequencies close to 0.5fS are attenuated.

1

0.5

0

1

0.5

0

2

3

4

5

6

7

I place “(sin x)/x” in quotes: With the argument properly scaled it is (sin πx)/(πx), but it is almost always pronounced sine-ecks-over-ecks, with argument scaling implicit.

You might consider a DAC’s boxcar waveform to be a “sample-and- hold” operation, but that term is normally used in conjunction with an A-to-D converter, or circuitry that lies in front of an ADC.

“(sin x)/x” correction

I have described how it is necessary for an analog reconstruction filter to follow digital-to-analog conversion. If the DAC produced an impulse “train” where the amplitude of each impulse was modulated by the corresponding code value, a classic lowpass filter would suffice: All would be well if the DAC output resembled my “point” graphs, with power at the sample instants and no power in between. Recall that a waveform comprising just unit impulses has uniform frequency response across the entire spectrum.

Unfortunately for analog reconstruction, a typical DAC does not produce an impulse waveform for each sample. It would be impractical to have a DAC with an impulse response, because signal power is proportional to the integral of the signal, and the amplitude of the impulses would have to be impractically high for the integral of the impulses to achieve adequate signal power. Instead, each converted sample value is held for the entire duration of the sample: A typical DAC produces a boxcar waveform. A boxcar waveform’s frequency response is described by the sinc function.

In Figure 20.29 above, the top graph is a sine wave at 0.44fS; the bottom graph shows the boxcar waveform produced by a conventional DAC. Even with

a high-quality reconstruction filter, whose response extends close to half the sampling rate, it is evident that

218

DIGITAL VIDEO AND HD ALGORITHMS AND INTERFACES

1.0

 

 

 

 

 

 

0.8

 

 

 

 

 

 

0.6

 

 

 

 

 

 

0.4

 

 

 

 

 

 

0.2

 

 

 

 

 

 

0

 

 

 

 

 

 

0

0.2

0.4

0.6

0.8

1.0

1.6

 

 

 

 

 

 

1.4

 

 

 

 

 

 

 

 

 

 

 

 

1.2

 

 

 

 

 

 

1.0

 

 

 

 

 

 

0.8

 

 

 

 

 

 

0.6

 

 

 

 

 

 

0.4

 

 

 

 

 

 

0.2

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0.2

0.4

0.6

0.8

1.0

 

 

Frequency, multiple of fs

 

Figure 20.30 “(sin x)/x” correction is necessary following (or in principle, preceding) digital-to- analog conversion when a DAC with a typical boxcar output waveform is used. The frequency response of a boxcar-waveform DAC is shown in the upper graph. The lower graph shows the response of the “(sin x)/x” correction filter necessary to compensate its high frequency falloff.

reconstruction by a boxcar function reduces the magnitude of high-frequency components of the signal.

The DAC’s holding of each sample value throughout the duration of its sample interval (zero-order hold, ZOH) corresponds to a filtering operation, with a frequency response of sinc. The top graph of Figure 20.30 shows the attenuation due to this phenomenon.

The effect is overcome by “(sin x)/xcorrection: The frequency response of the reconstruction filter is modified to include peaking corresponding to the reciprocal of sinc. In the passband, the filter’s response increases gradually to about 4 dB above its response at DC, to compensate the loss. Above the passband edge frequency, the response of the filter must decrease rapidly to produce a large attenuation near half the sampling frequency, to provide alias-free reconstruction. The bottom graph of Figure 20.30 shows the idealized response of a filter having “(sin x)/x” correction.

This chapter has detailed one-dimensional filtering. In Image digitization and reconstruction, I will introduce twoand three-dimensional sampling and filters.

CHAPTER 20

FILTERING AND SAMPLING

219

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