8202014 1 Click to edit Master subtitle style Sampling and Reconstruction Chapter 2 Ha Hoang Kha, Ph.D. Ho Chi Minh City University of Technology Email: hhkhahcmut.edu.vn Digital Signal Processing Sampling Content 2 Sampling theorem Sampling and Reconstruction Antialiasing prefilter Spectrum of sampling signals Analog reconstruction Ideal prefilter Practical prefilter Ideal reconstructor Practical reconstructor CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 2 Digital Signal Processing Review of Analog Signals 3 Fourier transform X(Ω) of x(t) is frequency spectrum of the signal Sampling and Reconstruction inverse Fourier transform Response of a linear system to an input signal x(t): X () x(t)e j tdt 2f x(t) X ( )ej t d2 y(t) h(t t)x(t)dt Y() H()X() Digital Signal Processing Review of Analog Signals 4 Frequency response of the system, defined as the Fourier transform of the impulse response h(t) Sampling and Reconstruction The steadystate sinusoidal response of the filter, defined as its response to sinusoidal inputs H() h(t)e j tdt j t x(t) e y(t) H()ejt H()ejt jarg H () CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 3 Digital Signal Processing Review of Analog Signals 5 If the input consists of the sum of two sinusoids of frequencies Sampling and Reconstruction After filtering, the steadystate output will be j t j t x t Ae 1 A e 2 ( ) 1 2 j t j t y(t) A1H(1)e 1 A2H(2)e 2 Digital Signal Processing Review of useful equations 6 Linear system Sampling and Reconstruction 1 sin( )sin( ) cos( ) cos( ) 2 a b a b a b 1 cos( )cos( ) cos( ) cos( ) 2 a b a b a b 0 0 0 1 sin(2 ) ( ) ( ) 2 f t j f f f f FT x t A f t ( ) cos(2 ) 0 x t ( ) Linear system h(t) H(f) y t x t h t ( ) ( ) ( ) X f ( ) Y f X f H f ( ) ( ) ( ) Especially, 0 0 0 1 cos(2 ) ( ) ( ) 2 f t f f f f FT y t A H f f t H f ( ) | ( ) | cos(2 arg( ( ))) 0 0 0 Fourier transform: Trigonometric formulas: 1 sin( )cos( ) sin( ) sin( ) 2 a b a b a b CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 4 Digital Signal Processing A typical signal processing system includes 3 stages: 1. Introduction 7 The digital processor can be programmed to perform signal processing operations such as filtering, spectrum estimation. Digital signal processor can be a general purpose computer, DSP chip or other digital hardware. Sampling and Reconstruction The analog signal is digitalized by an AD converter The digitalized samples are processed by a digital signal processor. The resulting output samples are converted back into analog by a DA converter. Digital Signal Processing 2. Analog to digital conversion 8 Analog to digital (AD) conversion is a threestep process. Sampler Quantizer Coder t=nT x x(t) x(nT)≡x(n) Q(n) 11010 AD converter n x Q(n) 000 001 010 011 100 101 110 111 t x(t) n x(n) Sampling and Reconstruction CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 5 Digital Signal Processing 3. Sampling 9 Sampling is to convert a continuous time signal into a discrete time signal. The analog signal is periodically measured at every T seconds x(n)≡x(nT)=x(t=nT), n=….2, 1, 0, 1, 2, 3…….. T: sampling interval or sampling period (second); fs=1T: sampling rate or sampling frequency (samplessecond or Hz) Sampling and Reconstruction Digital Signal Processing 3. Samplingexample 1 10 The analog signal x(t)=2cos(2πt) with t(s) is sampled at the rate fs=4 Hz. Find the discretetime signal x(n) ? Solution: Sampling and Reconstruction CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 6 Digital Signal Processing 3. Samplingexample 2 11 Consider the two analog sinusoidal signals Solution: 1 7 ( ) 2cos(2 ), 8 x t t 2 1 ( ) 2cos(2 ); ( ) 8 x t t t s These signals are sampled at the sampling frequency fs=1 Hz. Find the discretetime signals ? Sampling and Reconstruction at a sampling rate fs=1T results in a discretetime signal x(n). Digital Signal Processing 3. SamplingAliasing of Sinusoids 12 In general, the sampling of a continuoustime sinusoidal signal Remarks: We can that the frequencies fk=f0+kfs are indistinguishable from the frequency f0 after sampling and hence they are aliases of f0 x t A f t ( ) cos(2 ) 0 The sinusoids is sampled at fs , resulting in a discrete time signal xk(n). x t A f t k k ( ) cos(2 ) If f k=f0+kfs, k=0, ±1, ±2, …., then x(n)=xk(n) . Proof: (in class) Sampling and Reconstruction CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 7 Digital Signal Processing 4. Sampling TheoremSinusoids 13 Consider the analog signal where Ω is the frequency (rads) of the analog signal, and f=Ω2π is the frequency in cycless or Hz. The signal is sampled at the three rate fs =8f, fs=4f, and fs=2f. x t A t A ft ( ) cos( ) cos(2 ) Note that sec sec fs samples samples f cycles cycle To sample a single sinusoid properly, we must require fs 2 samples f cycle Fig: Sinusoid sampled at different rates Sampling and Reconstruction Digital Signal Processing 4. Sampling Theorem 14 For accurate representation of a signal x(t) by its time samples x(nT), two conditions must be met: 1) The signal x(t) must be bandlimitted, i.e., its frequency spectrum must be limited to f max . 2) The sampling rate fs must be chosen at least twice the maximum frequency fmax. f f s 2 max Fig: Typical bandlimited spectrum f s=2fmax is called Nyquist rate; fs2 is called Nyquist frequency; fs2, fs2 is Nyquist interval. Sampling and Reconstruction CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 8 Digital Signal Processing 4. Sampling Theorem 15 The values of f max and fs depend on the application Sampling and Reconstruction Application fmax fs Biomedical 1 KHz 2 KHz Speech 4 KHz 8 KHz Audio 20 KHz 40 KHz Video 4 MHz 8 MHz Digital Signal Processing 4. Sampling TheoremSpectrum Replication 16 Let where ( ) ( ) ( ) ( ) ( ) ( ) n x nT x t x t t nT x t s t ( ) ( ) n s t t nT s(t) is periodic, thus, its Fourier series are given by ( ) n j f nt 2 s n s t S e n 1 1 1 ( ) ( ) j f nt 2 s T T S t e dt t dt T T T 1 2 ( ) j f nt s n s t e T 1 2 ( ) ( ) ( ) ( ) j nf t s n x t x t s t x t e T 1 ( ) ( ) s n X f X f nf T where Thus, which results in Taking the Fourier transform of yields x t ( ) Observation: The spectrum of discretetime signal is a sum of the original spectrum of analog signal and its periodic replication at the interval f s. Sampling and Reconstruction CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 9 4. Sampling TheoremSpectrum Replication Digital Signal Processing 17 Fig: Typical badlimited spectrum f s2 ≥ fmax f s2 < fmax Fig: Aliasing caused by overlapping spectral replicas Fig: Spectrum replication caused by sampling Sampling and Reconstruction Digital Signal Processing 5. Ideal Analog reconstruction 18 Fig: Ideal reconstructor as a lowpass filter An ideal reconstructor acts as a lowpass filter with cutoff frequency equal to the Nyquist frequency fs2. X f X f H f X f a( ) ( ) ( ) ( ) An ideal reconstructor (lowpass filter) 2, 2 ( ) 0 T f f f s s H f otherwise Then Sampling and Reconstruction CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 10 Digital Signal Processing 5. Analog reconstructionExample 1 19 The analog signal x(t)=cos(20πt) is sampled at the sampling frequency fs=40 Hz. a) Plot the spectrum of signal x(t) ? b) Find the discrete time signal x(n) ? c) Plot the spectrum of signal x(n) ? d) The signal x(n) is an input of the ideal reconstructor, find the reconstructed signal xa(t) ? Sampling and Reconstruction Digital Signal Processing 5. Analog reconstructionExample 2 20 The analog signal x(t)=cos(100πt) is sampled at the sampling frequency fs=40 Hz. a) Plot the spectrum of signal x(t) ? b) Find the discrete time signal x(n) ? c) Plot the spectrum of signal x(n) ? d) The signal x(n) is an input of the ideal reconstructor, find the reconstructed signal xa(t) ? Sampling and Reconstruction CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 11 Digital Signal Processing 5. Analog reconstruction 21 Remarks: x a(t) contains only the frequency components that lie in the Nyquist interval (NI) fs2, fs2. x(t), f0 NI > x(n) > xa(t), fa=f0 sampling at fs ideal reconstructor x k(t), fk=f0+kfs> x(n) > xa(t), fa=f0 sampling at fs ideal reconstructor f f f a s mod( ) The frequency fa of reconstructed signal xa(t) is obtained by adding to or substracting from f0 (fk) enough multiples of fs until it lies within the Nyquist interval fs2, fs2.. That is Sampling and Reconstruction Digital Signal Processing 5. Analog reconstructionExample 3 22 The analog signal x(t)=10sin(4πt)+6sin(16πt) is sampled at the rate 20 Hz. Find the reconstructed signal xa(t) ? Sampling and Reconstruction CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 12 Digital Signal Processing 5. Analog reconstructionExample 4 23 Let x(t) be the sum of sinusoidal signals x(t)=4+3cos(πt)+2cos(2πt)+cos(3πt) where t is in milliseconds. Sampling and Reconstruction a) Determine the minimum sampling rate that will not cause any aliasing effects ? b) To observe aliasing effects, suppose this signal is sampled at half its Nyquist rate. Determine the signal xa(t) that would be aliased with x(t) ? Plot the spectrum of signal x(n) for this sampling rate? Digital Signal Processing 6. Ideal antialiasing prefilter 24 The signals in practice may not bandlimitted, thus they must be filtered by a lowpass filter Sampling and Reconstruction Fig: Ideal antialiasing prefilter CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 13 Digital Signal Processing 6. Practical antialiasing prefilter 25 Sampling and Reconstruction Fig: Practical antialiasing lowpass prefilter The Nyquist frequency fs2 is in the middle of transition region. A lowpass filter: fpass, fpass is the frequency range of interest for the application (fmax=fpass) The stopband frequency fstop and the minimum stopband attenuation A stop dB must be chosen appropriately to minimize the aliasing effects. f f f s pass stop Digital Signal Processing 6. Practical antialiasing prefilter 26 Sampling and Reconstruction The attenuation of the filter in decibels is defined as 10 0 ( ) ( ) 20log ( ) ( ) H f A f dB H f where f 0 is a convenient reference frequency, typically taken to be at DC for a lowpass filter. α 10 =A(10f)A(f) (dBdecade): the increase in attenuation when f is changed by a factor of ten. α 2 =A(2f)A(f) (dBoctave): the increase in attenuation when f is changed by a factor of two. Analog filter with order N, |H(f)|~1fN for large f, thus α10 =20N (dBdecade) and α10 =6N (dBoctave) CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 14 Digital Signal Processing 6. Antialiasing prefilterExample 27 Sampling and Reconstruction A sound wave has the form where t is in milliseconds. What is the frequency content of this signal ? Which parts of it are audible and why ? ( ) 2 cos(10 ) 2 cos(30 ) 2 cos(50 ) 2 cos(60 ) 2 cos(90 ) 2 cos(125 ) x t A t B t C t D t E t F t This signal is prefilter by an analog prefilter H(f). Then, the output y(t) of the prefilter is sampled at a rate of 40KHz and immediately reconstructed by an ideal analog reconstructor, resulting into the final analog output ya(t), as shown below: Digital Signal Processing 6. Antialiasing prefilterExample 28 Sampling and Reconstruction Determine the output signal y(t) and ya(t) in the following cases: a)When there is no prefilter, that is, H(f)=1 for all f. b)When H(f) is the ideal prefilter with cutoff fs2=20 KHz. c)When H(f) is a practical prefilter with specifications as shown below: The filter’s phase response is assumed to be ignored in this example. CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 15 Digital Signal Processing 7. Ideal and practical analog reconstructors 29 Sampling and Reconstruction An ideal reconstructor is an ideal lowpass filter with cutoff Nyquist frequency fs2. Digital Signal Processing 7. Ideal and practical analog reconstructors 30 Sampling and Reconstruction The ideal reconstructor has the impulse response: which is not realizable since its impulse response is not casual sin( f t) ( ) s s h t f t It is practical to use a staircase reconstructor CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 16 Digital Signal Processing 7. Ideal and practical analog reconstructors 31 Sampling and Reconstruction Fig: Frequency response of staircase reconstructor Digital Signal Processing 7. Practical reconstructorsantiimage postfilter 32 Sampling and Reconstruction An analog lowpass postfilter whose cutoff is Nyquist frequency fs2 is used to remove the surviving spectral replicas. Fig: Spectrum after postfilter Fig: Analog antiimage postfilter CuuDuongThanCong.com https:fb.comtailieudientucntt8202014 17 Digital Signal Processing 8. Homework 33 Sampling and Reconstruction Problems: provided in class CuuDuongThanCong.com https:fb.comtailieudientucntt
8/20/2014 Chapter Sampling and Reconstruction Ha Hoang Kha, Ph.D.Click to edit Master subtitle style Ho Chi Minh City University of Technology Email: hhkha@hcmut.edu.vn Content Sampling Sampling theorem Spectrum of sampling signals Antialiasing prefilter Ideal prefilter Practical prefilter Analog reconstruction Ideal reconstructor Practical reconstructor Digital Signal Processing CuuDuongThanCong.com Sampling and Reconstruction https://fb.com/tailieudientucntt 8/20/2014 Review of Analog Signals Fourier transform X(Ω) of x(t) is frequency spectrum of the signal X ( ) x(t )e jt dt 2f inverse Fourier transform d 2 Response of a linear system to an input signal x(t): x(t ) X ()e jt y (t ) h(t t ' ) x(t ' )dt ' Y () H () X () Digital Signal Processing Sampling and Reconstruction Review of Analog Signals Frequency response of the system, defined as the Fourier transform of the impulse response h(t) H () h(t )e jt dt The steady-state sinusoidal response of the filter, defined as its response to sinusoidal inputs x(t ) e jt Digital Signal Processing CuuDuongThanCong.com y(t ) H ()e jt H () e jt j arg H () Sampling and Reconstruction https://fb.com/tailieudientucntt 8/20/2014 Review of Analog Signals If the input consists of the sum of two sinusoids of frequencies x(t ) A1e j1t A2e j2t After filtering, the steady-state output will be y(t ) A1H (1 )e j1t A2 H (2 )e j2t Digital Signal Processing Sampling and Reconstruction Review of useful equations Linear system Linear system h(t) H(f) x(t ) X( f ) y(t ) x(t ) h(t ) Y ( f ) X ( f )H ( f ) Especially, x(t ) A cos(2 f0t ) y(t ) A | H ( f0 ) | cos(2 f0t arg( H ( f ))) FT cos(2 f 0t ) [ ( f f ) ( f f )] FT sin(2 f 0t ) j[ ( f f ) ( f f )] Trigonometric formulas: cos(a) cos(b) [cos(a b) cos(a b)] sin(a)sin(b) [cos(a b) cos(a b)] sin(a) cos(b) [sin(a b) sin(a b)] Fourier transform: Digital Signal Processing CuuDuongThanCong.com Sampling and Reconstruction https://fb.com/tailieudientucntt 8/20/2014 Introduction A typical signal processing system includes stages: The analog signal is digitalized by an A/D converter The digitalized samples are processed by a digital signal processor The digital processor can be programmed to perform signal processing operations such as filtering, spectrum estimation Digital signal processor can be a general purpose computer, DSP chip or other digital hardware The resulting output samples are converted back into analog by a D/A converter Digital Signal Processing Sampling and Reconstruction Analog to digital conversion Analog to digital (A/D) conversion is a three-step process x(t) Sampler t=nT x(t) x(nT)≡x(n) Quantizer xQ(n) Coder A/D converter x(n) t Digital Signal Processing CuuDuongThanCong.com 11010 n 111 xQ(n) 110 101 100 011 010 001 000 n Sampling and Reconstruction https://fb.com/tailieudientucntt 8/20/2014 Sampling Sampling is to convert a continuous time signal into a discrete time signal The analog signal is periodically measured at every T seconds x(n)≡x(nT)=x(t=nT), n=….-2, -1, 0, 1, 2, 3…… T: sampling interval or sampling period (second); fs=1/T: sampling rate or sampling frequency (samples/second or Hz) Digital Signal Processing Sampling and Reconstruction Sampling-example The analog signal x(t)=2cos(2πt) with t(s) is sampled at the rate fs=4 Hz Find the discrete-time signal x(n) ? Solution: Digital Signal Processing CuuDuongThanCong.com 10 Sampling and Reconstruction https://fb.com/tailieudientucntt 8/20/2014 Sampling-example Consider the two analog sinusoidal signals x1 (t ) 2cos(2 t ), x2 (t ) 2cos(2 t ); t ( s) 8 These signals are sampled at the sampling frequency fs=1 Hz Find the discrete-time signals ? Solution: Digital Signal Processing 11 Sampling and Reconstruction Sampling-Aliasing of Sinusoids In general, the sampling of a continuous-time sinusoidal signal x(t ) A cos(2 f0t ) at a sampling rate fs=1/T results in a discretetime signal x(n) The sinusoids xk (t ) A cos(2 f k t ) is sampled at fs , resulting in a discrete time signal xk(n) If fk=f0+kfs, k=0, ±1, ±2, …., then x(n)=xk(n) Proof: (in class) Remarks: We can that the frequencies fk=f0+kfs are indistinguishable from the frequency f0 after sampling and hence they are aliases of f0 Digital Signal Processing CuuDuongThanCong.com 12 Sampling and Reconstruction https://fb.com/tailieudientucntt 8/20/2014 Sampling Theorem-Sinusoids Consider the analog signal x(t ) A cos(t ) A cos(2 ft ) where Ω is the frequency (rad/s) of the analog signal, and f=Ω/2π is the frequency in cycles/s or Hz The signal is sampled at the three rate fs=8f, fs=4f, and fs=2f Fig: Sinusoid sampled at different rates Note that f s samples / sec samples f cycles / sec cycle To sample a single sinusoid properly, we must require f s samples f Digital Signal Processing 13 cycle Sampling and Reconstruction Sampling Theorem For accurate representation of a signal x(t) by its time samples x(nT), two conditions must be met: 1) The signal x(t) must be bandlimitted, i.e., its frequency spectrum must be limited to fmax Fig: Typical bandlimited spectrum 2) The sampling rate fs must be chosen at least twice the maximum f s f max frequency fmax fs=2fmax is called Nyquist rate; fs/2 is called Nyquist frequency; [-fs/2, fs/2] is Nyquist interval Digital Signal Processing CuuDuongThanCong.com 14 Sampling and Reconstruction https://fb.com/tailieudientucntt 8/20/2014 Sampling Theorem The values of fmax and fs depend on the application Application fmax fs Biomedical KHz KHz Speech KHz KHz Audio 20 KHz 40 KHz Video MHz MHz Digital Signal Processing 15 Sampling and Reconstruction Sampling Theorem-Spectrum Replication Let x(nT ) x (t ) x(t ) (t nT ) x(t )s(t ) where s(t ) n (t nT ) n s(t) is periodic, thus, its Fourier series are given by s(t ) Se n Thus, s(t ) n j 2 f s nt where Sn j 2 f s nt e T n which results in 1 (t )e j 2 fs nt dt (t )dt TT TT T x (t ) x(t ) s(t ) x(t )e j 2 nfst T n Taking the Fourier transform of x (t ) yields X ( f ) X ( f nf s ) T n Observation: The spectrum of discrete-time signal is a sum of the original spectrum of analog signal and its periodic replication at the interval fs Digital Signal Processing CuuDuongThanCong.com 16 Sampling and Reconstruction https://fb.com/tailieudientucntt 8/20/2014 Sampling Theorem-Spectrum Replication fs/2 ≥ fmax Fig: Spectrum replication caused by sampling Fig: Typical badlimited spectrum fs/2 < fmax Fig: Aliasing caused by overlapping spectral replicas Digital Signal Processing 17 Sampling and Reconstruction Ideal Analog reconstruction Fig: Ideal reconstructor as a lowpass filter An ideal reconstructor acts as a lowpass filter with cutoff frequency equal to the Nyquist frequency fs/2 T 0 An ideal reconstructor (lowpass filter) H ( f ) Then Digital Signal Processing CuuDuongThanCong.com f [ f s / 2, f s / 2] otherwise X a ( f ) X ( f )H ( f ) X ( f ) 18 Sampling and Reconstruction https://fb.com/tailieudientucntt 8/20/2014 Analog reconstruction-Example The analog signal x(t)=cos(20πt) is sampled at the sampling frequency fs=40 Hz a) Plot the spectrum of signal x(t) ? b) Find the discrete time signal x(n) ? c) Plot the spectrum of signal x(n) ? d) The signal x(n) is an input of the ideal reconstructor, find the reconstructed signal xa(t) ? Digital Signal Processing 19 Sampling and Reconstruction Analog reconstruction-Example The analog signal x(t)=cos(100πt) is sampled at the sampling frequency fs=40 Hz a) Plot the spectrum of signal x(t) ? b) Find the discrete time signal x(n) ? c) Plot the spectrum of signal x(n) ? d) The signal x(n) is an input of the ideal reconstructor, find the reconstructed signal xa(t) ? Digital Signal Processing CuuDuongThanCong.com 20 Sampling and Reconstruction https://fb.com/tailieudientucntt 10 8/20/2014 Analog reconstruction Remarks: xa(t) contains only the frequency components that lie in the Nyquist interval (NI) [-fs//2, fs/2] sampling at fs ideal reconstructor x(t), f0 NI > x(n) > xa(t), fa=f0 sampling at fs ideal reconstructor xk(t), fk=f0+kfs > x(n) > xa(t), fa=f0 The frequency fa of reconstructed signal xa(t) is obtained by adding to or substracting from f0 (fk) enough multiples of fs until it lies within the Nyquist interval [-fs//2, fs/2] That is f a f mod( f s ) Digital Signal Processing 21 Sampling and Reconstruction Analog reconstruction-Example The analog signal x(t)=10sin(4πt)+6sin(16πt) is sampled at the rate 20 Hz Find the reconstructed signal xa(t) ? Digital Signal Processing CuuDuongThanCong.com 22 Sampling and Reconstruction https://fb.com/tailieudientucntt 11 8/20/2014 Analog reconstruction-Example Let x(t) be the sum of sinusoidal signals x(t)=4+3cos(πt)+2cos(2πt)+cos(3πt) where t is in milliseconds a) Determine the minimum sampling rate that will not cause any aliasing effects ? b) To observe aliasing effects, suppose this signal is sampled at half its Nyquist rate Determine the signal xa(t) that would be aliased with x(t) ? Plot the spectrum of signal x(n) for this sampling rate? Digital Signal Processing 23 Sampling and Reconstruction Ideal antialiasing prefilter The signals in practice may not bandlimitted, thus they must be filtered by a lowpass filter Fig: Ideal antialiasing prefilter Digital Signal Processing CuuDuongThanCong.com 24 Sampling and Reconstruction https://fb.com/tailieudientucntt 12 8/20/2014 Practical antialiasing prefilter A lowpass filter: [-fpass, fpass] is the frequency range of interest for the application (fmax=fpass) The Nyquist frequency fs/2 is in the middle of transition region The stopband frequency fstop and the minimum stopband attenuation Astop dB must be chosen appropriately to minimize the aliasing effects f s f pass f stop Fig: Practical antialiasing lowpass prefilter Digital Signal Processing 25 Sampling and Reconstruction Practical antialiasing prefilter The attenuation of the filter in decibels is defined as A( f ) 20log10 H( f ) (dB) H ( f0 ) where f0 is a convenient reference frequency, typically taken to be at DC for a lowpass filter α10 =A(10f)-A(f) (dB/decade): the increase in attenuation when f is changed by a factor of ten α2 =A(2f)-A(f) (dB/octave): the increase in attenuation when f is changed by a factor of two Analog filter with order N, |H(f)|~1/fN for large f, thus α10 =20N (dB/decade) and α10 =6N (dB/octave) Digital Signal Processing CuuDuongThanCong.com 26 Sampling and Reconstruction https://fb.com/tailieudientucntt 13 8/20/2014 Antialiasing prefilter-Example A sound wave has the form x(t ) A cos(10 t ) B cos(30 t ) 2C cos(50 t ) D cos(60 t ) E cos(90 t ) F cos(125 t ) where t is in milliseconds What is the frequency content of this signal ? Which parts of it are audible and why ? This signal is prefilter by an analog prefilter H(f) Then, the output y(t) of the prefilter is sampled at a rate of 40KHz and immediately reconstructed by an ideal analog reconstructor, resulting into the final analog output ya(t), as shown below: Digital Signal Processing 27 Sampling and Reconstruction Antialiasing prefilter-Example Determine the output signal y(t) and ya(t) in the following cases: a)When there is no prefilter, that is, H(f)=1 for all f b)When H(f) is the ideal prefilter with cutoff fs/2=20 KHz c)When H(f) is a practical prefilter with specifications as shown below: The filter’s phase response is assumed to be ignored in this example Digital Signal Processing CuuDuongThanCong.com 28 Sampling and Reconstruction https://fb.com/tailieudientucntt 14 8/20/2014 Ideal and practical analog reconstructors An ideal reconstructor is an ideal lowpass filter with cutoff Nyquist frequency fs/2 Digital Signal Processing 29 Sampling and Reconstruction Ideal and practical analog reconstructors sin( f t) s The ideal reconstructor has the impulse response: h(t ) f t which is not realizable since its impulse response is not casual s It is practical to use a staircase reconstructor Digital Signal Processing CuuDuongThanCong.com 30 Sampling and Reconstruction https://fb.com/tailieudientucntt 15 8/20/2014 Ideal and practical analog reconstructors Fig: Frequency response of staircase reconstructor Digital Signal Processing 31 Sampling and Reconstruction Practical reconstructors-antiimage postfilter An analog lowpass postfilter whose cutoff is Nyquist frequency fs/2 is used to remove the surviving spectral replicas Fig: Analog anti-image postfilter Fig: Spectrum after postfilter Digital Signal Processing CuuDuongThanCong.com 32 Sampling and Reconstruction https://fb.com/tailieudientucntt 16 8/20/2014 Homework Problems: provided in class Digital Signal Processing CuuDuongThanCong.com 33 Sampling and Reconstruction https://fb.com/tailieudientucntt 17 ... signal The analog signal is periodically measured at every T seconds x(n)≡x(nT)=x(t=nT), n=… .-2 , -1 , 0, 1, 2, 3…… T: sampling interval or sampling period (second); fs=1/T: sampling rate or... Processing Sampling and Reconstruction Sampling-example The analog signal x(t)=2cos(2πt) with t(s) is sampled at the rate fs=4 Hz Find the discrete-time signal x(n) ? Solution: Digital Signal... Sampling-example Consider the two analog sinusoidal signals x1 (t ) 2cos(2 t ), x2 (t ) 2cos(2 t ); t ( s) 8 These signals are sampled at the sampling frequency fs=1 Hz Find the discrete-time