《数字信号处理.doc》由会员分享，可在线阅读，更多相关《数字信号处理.doc（28页珍藏版）》请在三一办公上搜索。
1、Understanding Digital Signal Processing: Periodic Sampling By Richard G. Lyons Dec 3, 2004 Sample Chapter is provided courtesy of Prentice Hall ProfessionalPeriodic sampling, the process of representing a continuous signal with a sequence of discrete data values, pervades the field of digital signal
2、 processing. This chapter covers that process in detail, and how you can use it in your business. Periodic sampling, the process of representing a continuous signal with a sequence of discrete data values, pervades the field of digital signal processing. In practice, sampling is performed by applyin
3、g a continuous signal to an analog-to-digital (A/D) converter whose output is a series of digital values. Because sampling theory plays an important role in determining the accuracy and feasibility of any digital signal processing scheme, we need a solid appreciation for the often misunderstood effe
4、cts of periodic sampling. With regard to sampling, the primary concern is just how fast a given continuous signal must be sampled in order to preserve its information content. We can sample a continuous signal at any sample rate we wish, and well get a series of discrete valuesbut the question is, h
5、ow well do these values represent the original signal? Lets learn the answer to that question and, in doing so, explore the various sampling techniques used in digital signal processing. 2.1 ALIASING: SIGNAL AMBIGUITY IN THE FREQUENCY DOMAINThere is a frequency-domain ambiguity associated with discr
6、ete-time signal samples that does not exist in the continuous signal world, and we can appreciate the effects of this uncertainty by understanding the sampled nature of discrete data. By way of example, suppose you were given the following sequence of values, and were told that they represent instan
7、taneous values of a time-domain sinewave taken at periodic intervals. Next, you were asked to draw that sinewave. Youd start by plotting the sequence of values shown by the dots in Figure 2-1(a). Next, youd be likely to draw the sinewave, illustrated by the solid line in Figure 2-1(b), that passes t
8、hrough the points representing the original sequence. Figure 2-1 Frequency ambiguity: (a) discrete-time sequence of values; (b) two different sinewaves that pass through the points of the discrete sequence.Another person, however, might draw the sinewave shown by the shaded line in Figure 2-1(b). We
9、 see that the original sequence of values could, with equal validity, represent sampled values of both sinewaves. The key issue is that, if the data sequence represented periodic samples of a sinewave, we cannot unambiguously determine the frequency of the sinewave from those sample values alone. Re
10、viewing the mathematical origin of this frequency ambiguity enables us not only to deal with it, but to use it to our advantage. Lets derive an expression for this frequency-domain ambiguity and, then, look at a few specific examples. Consider the continuous time-domain sinusoidal signal defined as
11、Equation 2-1 This x(t) signal is a garden variety sinewave whose frequency is fo Hz. Now lets sample x(t) at a rate of fs samples/s, i.e., at regular periods of ts seconds where ts = 1/fs. If we start sampling at time t = 0, we will obtain samples at times 0ts, 1ts, 2ts, and so on. So, from Eq. (2-1
12、), the first n successive samples have the values Equation 2-2 Equation (2-2) defines the value of the nth sample of our x(n) sequence to be equal to the original sinewave at the time instant nts. Because two values of a sinewave are identical if theyre separated by an integer multiple of 2 radians,
13、 i.e., sin() = sin(+2m) where m is any integer, we can modify Eq. (2-2) as Equation 2-3 If we let m be an integer multiple of n, m = kn, we can replace the m/n ratio in Eq. (2-3) with k so that Equation 2-4 Because fs = 1/ts, we can equate the x(n) sequences in Eqs. (2-2) and (2-4) as Equation 2-5 T
14、he fo and (fo+kfs) factors in Eq. (2-5) are therefore equal. The implication of Eq. (2-5) is critical. It means that an x(n) sequence of digital sample values, representing a sinewave of fo Hz, also exactly represents sinewaves at other frequencies, namely, fo + kfs. This is one of the most importan
15、t relationships in the field of digital signal processing. Its the thread with which all sampling schemes are woven. In words, Eq. (2-5) states that When sampling at a rate of fs samples/s, if k is any positive or negative integer, we cannot distinguish between the sampled values of a sinewave of fo
16、 Hz and a sinewave of (fo+kfs) Hz.Its true. No sequence of values stored in a computer, for example, can unambiguously represent one and only one sinusoid without additional information. This fact applies equally to A/D-converter output samples as well as signal samples generated by computer softwar
17、e routines. The sampled nature of any sequence of discrete values makes that sequence also represent an infinite number of different sinusoids. Equation (2-5) influences all digital signal processing schemes. Its the reason that, although weve only shown it for sinewaves, well see in Chapter 3 that
18、the spectrum of any discrete series of sampled values contains periodic replications of the original continuous spectrum. The period between these replicated spectra in the frequency domain will always be fs, and the spectral replications repeat all the way from DC to daylight in both directions of
19、the frequency spectrum. Thats because k in Eq. (2-5) can be any positive or negative integer. (In Chapters 5 and 6, well learn that Eq. (2-5) is the reason that all digital filter frequency responses are periodic in the frequency domain and is crucial to analyzing and designing a popular type of dig
20、ital filter known as the infinite impulse response filter.) To illustrate the effects of Eq. (2-5), lets build on Figure 2-1 and consider the sampling of a 7-kHz sinewave at a sample rate of 6 kHz. A new sample is determined every 1/6000 seconds, or once every 167 microseconds, and their values are
21、shown as the dots in Figure 2-2(a). Figure 2-2 Frequency ambiguity effects of Eq. (2-5): (a) sampling a 7-kHz sinewave at a sample rate of 6 kHz; (b) sampling a 4-kHz sinewave at a sample rate of 6 kHz; (c) spectral relationships showing aliasing of the 7- and 4-kHz sinewaves.Notice that the sample
22、values would not change at all if, instead, we were sampling a 1-kHz sinewave. In this example fo = 7 kHz, fs = 6 kHz, and k = 1 in Eq. (2-5), such that fo+kfs = 7+(16) = 1 kHz. Our problem is that no processing scheme can determine if the sequence of sampled values, whose amplitudes are represented
23、 by the dots, came from a 7-kHz or a 1-kHz sinusoid. If these amplitude values are applied to a digital process that detects energy at 1 kHz, the detector output would indicate energy at 1 kHz. But we know that there is no 1-kHz tone thereour input is a spectrally pure 7-kHz tone. Equation (2-5) is
24、causing a sinusoid, whose name is 7 kHz, to go by the alias of 1 kHz. Asking someone to determine which sinewave frequency accounts for the sample values in Figure 2-2(a) is like asking them When I add two numbers I get a sum of four. What are the two numbers? The answer is that there is an infinite
25、 number of number pairs that can add up to four. Figure 2-2(b) shows another example of frequency ambiguity, that well call aliasing, where a 4-kHz sinewave could be mistaken for a 2-kHz sinewave. In Figure 2-2(b), fo = 4 kHz, fs = 6 kHz, and k = 1 in Eq. (2-5), so that fo+kfs = 4+(1 6) = 2 kHz. Aga
26、in, if we examine a sequence of numbers representing the dots in Figure 2-2(b), we could not determine if the sampled sinewave was a 4-kHz tone or a 2-kHz tone. (Although the concept of negative frequencies might seem a bit strange, it provides a beautifully consistent methodology for predicting the
27、 spectral effects of sampling. Chapter 8 discusses negative frequencies and how they relate to real and complex signals.) Now, if we restrict our spectral band of interest to the frequency range of fs/2 Hz, the previous two examples take on a special significance. The frequency fs/2 is an important
28、quantity in sampling theory and is referred to by different names in the literature, such as critical Nyquist, half Nyquist, and folding frequency. A graphical depiction of our two frequency aliasing examples is provided in Figure 2-2(c). Were interested in signal components that are aliased into th
29、e frequency band between fs/2 and +fs/2. Notice in Figure 2-2(c) that, within the spectral band of interest (3 kHz, because fs = 6 kHz), there is energy at 2 kHz and +1 kHz, aliased from 4 kHz and 7 kHz, respectively. Note also that the vertical positions of the dots in Figure 2-2(c) have no amplitu
30、de significance but that their horizontal positions indicate which frequencies are related through aliasing. A general illustration of aliasing is provided in the sharks tooth pattern in Figure 2-3(a). Note how the peaks of the pattern are located at integer multiples of fs Hz. The pattern shows how
31、 signals residing at the intersection of a horizontal line and a sloped line will be aliased to all of the intersections of that horizontal line and all other lines with like slopes. For example, the pattern in Figure 2-3(b) shows that our sampling of a 7-kHz sinewave at a sample rate of 6 kHz will
32、provide a discrete sequence of numbers whose spectrum ambiguously represents tones at 1 kHz, 7 kHz, 13 kHz, 19 kHz, etc. Lets pause for a moment and let these very important concepts soak in a bit. Again, discrete sequence representations of a continuous signal have unavoidable ambiguities in their
33、frequency domains. These ambiguities must be taken into account in all practical digital signal processing algorithms. Figure 2-3 Sharks tooth pattern: (a) aliasing at multiples of the sampling frequency; (b) aliasing of the 7-kHz sinewave to 1 kHz, 13 kHz, and 19 kHz.OK, lets review the effects of
34、sampling signals that are more interesting than just simple sinusoids.2.2 SAMPLING LOW-PASS SIGNALSConsider sampling a continuous real signal whose spectrum is shown in Figure 2-4(a). Notice that the spectrum is symmetrical about zero Hz, and the spectral amplitude is zero above +B Hz and below B Hz
35、, i.e., the signal is band-limited. (From a practical standpoint, the term band-limited signal merely implies that any signal energy outside the range of B Hz is below the sensitivity of our system.) Given that the signal is sampled at a rate of fs samples/s, we can see the spectral replication effe
36、cts of sampling in Figure 2-4(b), showing the original spectrum in addition to an infinite number of replications, whose period of replication is fs Hz. (Although we stated in Section 1.1 that frequency-domain representations of discrete-time sequences are themselves discrete, the replicated spectra
37、 in Figure 2-4(b) are shown as continuous lines, instead of discrete dots, merely to keep the figure from looking too cluttered. Well cover the full implications of discrete frequency spectra in Chapter 3.) Figure 2-4 Spectral replications: (a) original continuous signal spectrum; (b) spectral repli
38、cations of the sampled signal when fs/2 B; (c) frequency overlap and aliasing when the sampling rate is too low because fs/2 B.Lets step back a moment and understand Figure 2-4 for all its worth. Figure 2-4(a) is the spectrum of a continuous signal, a signal that can only exist in one of two forms.
39、Either its a continuous signal that can be sampled, through A/D conversion, or it is merely an abstract concept such as a mathematical expression for a signal. It cannot be represented in a digital machine in its current band-limited form. Once the signal is represented by a sequence of discrete sam
40、ple values, its spectrum takes the replicated form of Figure 2-4(b). The replicated spectra are not just figments of the mathematics; they exist and have a profound effect on subsequent digital signal processing. The replications may appear harmless, and its natural to ask, Why care about spectral r
41、eplications? Were only interested in the frequency band within fs/2. Well, if we perform a frequency translation operation or induce a change in sampling rate through decimation or interpolation, the spectral replications will shift up or down right in the middle of the frequency range of interest f
42、s/2 and could cause problems1. Lets see how we can control the locations of those spectral replications. In practical A/D conversion schemes, fs is always greater than 2B to separate spectral replications at the folding frequencies of fs/2. This very important relationship of fs 2B is known as the N
43、yquist criterion. To illustrate why the term folding frequency is used, lets lower our sampling frequency to fs = 1.5B Hz. The spectral result of this undersampling is illustrated in Figure 2-4(c). The spectral replications are now overlapping the original baseband spectrum centered about zero Hz. L
44、imiting our attention to the band fs/2 Hz, we see two very interesting effects. First, the lower edge and upper edge of the spectral replications centered at +fs and fs now lie in our band of interest. This situation is equivalent to the original spectrum folding to the left at +fs/2 and folding to
45、the right at fs/2. Portions of the spectral replications now combine with the original spectrum, and the result is aliasing errors. The discrete sampled values associated with the spectrum of Figure 2-4(c) no longer truly represent the original input signal. The spectral information in the bands of
46、B to B/2 and B/2 to B Hz has been corrupted. We show the amplitude of the aliased regions in Figure 2-4(c) as dashed lines because we dont really know what the amplitudes will be if aliasing occurs. The second effect illustrated by Figure 2-4(c) is that the entire spectral content of the original co
47、ntinuous signal is now residing in the band of interest between fs/2 and +fs/2. This key property was true in Figure 2-4(b) and will always be true, regardless of the original signal or the sample rate. This effect is particularly important when were digitizing (A/D converting) continuous signals. I
48、t warns us that any signal energy located above +B Hz and below B Hz in the original continuous spectrum of Figure 2-4(a) will always end up in the band of interest after sampling, regardless of the sample rate. For this reason, continuous (analog) low-pass filters are necessary in practice. We illu
49、strate this notion by showing a continuous signal of bandwidth B accompanied by noise energy in Figure 2-5(a). Sampling this composite continuous signal at a rate thats greater than 2B prevents replications of the signal of interest from overlapping each other, but all of the noise energy still ends up in the range between fs/2 and +fs/2 of our discrete spectrum shown in Figure 2-5(b). This p
