|
|
- Horace Wilson
- 5 years ago
- Views:
Transcription
1 SECTION FOR DIGITAL SIGNAL PROCESSING DEPARTMENT OF MATHEMATICAL MODELLING TECHNICAL UNIVERSITY OF DENMARK Course Digital Signal Processing: Solutions to Problems in Proakis and Manolakis, Digital Signal Processing, 3rd Edition Jan Larsen 3rd Edition c September 1998 by Jan Larsen
2
3 Preface This note is a supplement to the textbook used in the DTU Course Digital Signal: J.G. Proakis & D.G. Manolakis: Digital Signal Processing: Principles, Algorithms and Applications, 3rd edition, Upper Saddle River, New Jersey: Prentice-Hall, Inc., Unless anything else is mentioned, all page and equation references are with respect to the course textbook. The note is partially based on material by Simon Boel Pedersen used in the former DTU course 4232 Digital Signal Processing. Jan Larsen Lyngby, September 1998 The manuscript was typeset in 12 points Times Roman and Pandora using L A TE2" i
4
5 SOLUTION TO PROBLEM 9.1 The center frequency is F c = 50 and the bandwidth is B = 20. According to Ch the minimum sampling frequency is given by F s min = 1=T 0 = 2Br 0 =r where r 0 = F c =B +1=2 =50=20 + 1=2 =3and r = br 0 c =3. That is, F s min = 2B = 40. In general, we require that the p'th and the (p + 1)'th replication of the spectrum in the sampled signal does not interfere with the original signal spectrum as shown in Fig j(f )j B- 6 F c j(f ; pf s )j - F pf s ; F c pf s j(f ; (p +1)F s )j 6 - F (p +1)F s ; F c (p +1)F s - F Fig. 9.1 This gives the inequalities: F c ; B=2 > pf s ; F c + B=2 F c + B=2 < (p +1)F s ; F c ; B=2 Here we assume F c >B=2 and p is an integer. Rewriting yields: p < 2F c ; B F s p +1 > 2F c + B F s Thus B F s ; 1 <p; 2F C F s < ; B F s 1
6 which implies the natural condition: F s > 2B The minimum sampling frequency is achieved when p is maximized. natural condition, the rst inequality gives p< B F s 2Fc B ; 1 < F c B ; 1 2 Using the and the second gives m p +1> F c B > B F s 2Fc B +1 p> F c B ; 1 2 : Choosing p as the integer, p max = b Fc ; 1 c, ensures that the replicated parts of the B 2 sampled signal are distributed equally along the frequency axis. Applying p max in the inequalities results in: 2B r0 r <F s < 2B s0 s where r 0 = F c =B +1=2, r = br 0 c, s 0 = F c =B ; 1=2, and s = bs 0 c. SOLUTION TO PROBLEM 9.4 Solution to Question (a): Using the sampling frequency F s = 2B, (f) will be a periodic replication of the original spectrum a (F ) with period F s cf. Eq. (4.2.85) (f) = F s P 1 k=;1 a ((f ; k)f s ). The spectrum of y(n), say Y (f) is obtained by using the fact that multiplication in the time-domain corresponds to convolution in the frequency domain. Finally, y 1 (t) is found by low-pass ltering y(n) using a passband F 2 [0 B]. a (F ) has ideal low-pass lter characteristic, say a (f) =1jF jb, and zero otherwise. The spectrum (f) =F s, 8jfj 1=2, i.e., constant. Thus x(n) =F s (n). y(n) = x 2 (n) = F 2 s (n) with spectrum Y (f) = F 2 s, 8jfj 1=2. Consequently, by Eq. (4.2.87) Y 1 (F )=F s for jf jb, and zero otherwise. That is, y 1 (t) =F s x 1 (t). The spectrum of s a (t) =x 2 a(t) is found using convolution in the frequency domain, i.e., Z 1 ( 2B ;jfj ;2B F 2B S a (F )= a(s) a (F ; S)dS = ;1 0 otherwise The sketch of the spectra is shown in Fig
7 S a (F ) 6 2B ;2B ;B B - F 2B Fig (a) S(f) 6 2BF s f = F=2B ;1 ;1=2 1=2 1 Fig (b) Y 1 (F )=Y 2 (F ) 6 2B ;2B ;B B 2B - F Fig (c) 3
8 Thus y 2 (t) =2Bx(t) =2By 1 (t)=f s = y 1 (t). Solution to Question (b): x a (t) =cos(2f 0 t) with spectrum a (F )=0:5(F F 0 ), where F 0 =20Hz. Using F s =50Hz> 2F 0 the sampling theorem is fullled. (F=F s ) has components of amplitude F s =2 for F = (20 + p50) Hz, p = By using the rule of convolution in the frequency domain Y (F=F s ) has components for F = 10 Hz with amplitude F s =4 and a DC component (F =0)with amplitude F s =2. By ideal D/A we have y 1 (t) = 1=2 +1=2 cos(2(10 Hz)t). s a (t) = x 2 a(t) = cos 2 (2F 0 t) = 1=2 +1=2 cos(22f 0 t), i.e., a (F ) has a DC component with amplitude 1=2 and two components at F = 40 Hz with amplitude 1=4. Sampling with F s = 50 Hz will now cause aliasing. S(F=F s ) will have components at F = 0 with amplitude F s =2 and two components at F = 10 Hz with amplitude F s =4. Ideal D/A gives y 2 (t) =y 1 (t). When sampling with F s = 30 Hz the sampling theorem is not fullled. x(n) will have components at F = 10 Hz with amplitudes F s =2. y(n) consequently have components at F 10 Hz with amplitudes F s =4 at a DC component with amplitude 1=2. Reconstruction gives y 1 (t) =1=2+1=2 cos(2(10 Hz)t). Sampling s a (t) with F s = 30 Hz results in s(n) with a DC component with amplitude F s =2 and two components F = 10 Hz with amplitudes F s =4. Reconstruction gives y 2 (t) =y 1 (t). SOLUTION TO PROBLEM 4.83 Solution to Question (a): From Fig. P4.83 we notice that frequency-inversion scrambling is done by shifting the spectrum by f =1=2 (! = ), thus Y (f) =(f ;1=2). The corresponding time domain operation is found by using the frequency shift property { or modulation theorem Ch : y(n) =(;1) n x(n). Solution to Question (b): The unscrambler is a similar frequency shift of f = 1=2, thus x(n) =(;1) n y(n) =(;1) 2n x(n) =x(n). SOLUTION TO PROBLEM 9.5 A analog TV signal x a (t) is corrupted by an echo signal. The received signal is: s a (t) =x a (t)+x a (t ; ) jj < 1 > 0 By using the time shifting property, the spectrum of s a (t) is given by: S a (F ) = a (F )+ a (F )e ;j2f = a (F ) 1+e ;j2f 4
9 Since the bandwidth of x a (t) isb, clearly the bandwidth of s a (t) isalsob. Choosing the sampling frequency in accordance with the sampling theorem F s > 2B and such that = n 0 =F s = n 0 T, where n 0 is a positive integer, and T is the sampling interval, the digital equivalent becomes: s(n) =x(n)+x(n ; n 0 ) with spectrum S(f) =(f) 1+e ;j2fn 0 In order to get y(n) = x(n) and consequently y a (t) = x a (t) when using an ideal reconstruction, i.e., a BL interpolator with cut-o frequency B, we need a lter H(f) which is the inverse of 1 + e ;j2fn 0, thus: H(f) = 1 1+e ;j2fn 0 The corresponding transfer function is H(z) = 1=(1 + z ;n 0 ), i.e., the lter is an all-pole lter with poles found by solving the binomial equation of degree n 0 : The n 0 poles is therefore given by: z p = 8 >< >: q z n 0 = ; n 0 jje j(2p=n 0) ;1 <<0 p =0 1 n 0 ; 1 q jje j(=n 0 +2p=n 0 ) 0 <1 p =0 1 n 0 ; 1 n 0 The poles are inside the unit circle, i.e., jz p j < 1forjj < 1, thus the lter is (BIBO) stable and causal. If n 0 =1the lter is a low-pass lter for <0cf.Ch For >0 the lter is the corresponding high-pass lter found by using the modulation theorem (frequency shifting) Ch This is easily seen as h(n)(;1) n corresponds to H(f ;1=2). Thus shifting H(f) gives: H(f ; 1=2) = 1 1+e ;j2(f;1=2) = 1 1+e j e ;j2f = 1 1 ; e ;j2f Denote by H n0 (f) the lter using n 0. The lters of n 0 is easily found from H 1 (f). Since H 1 (z) =1=(1 ; z ;1 ), then H n0 (z) =H 1 (z n 0 ), see Ch and That is, the spectrum H n0 (f) =H n0 (z)j jzj=e j2f is found by H n0 (f) =H 1 (n 0 f) That is, H n0 (f) is a comb lter cf. Ch spectra. Fig and show typical 5
10 25 n 0 = α=0.1 α=0.5 α=0.9 H(f) (db) f Fig n 0 = α=0.1 α=0.5 α=0.9 H(f) (db) f Fig
11 SOLUTION TO PROBLEM 9.7 x(n) is a stationary random signal with mean E[x(n)] = 0, variance V [x(n)] = 2 x and autocorrelation function, xx (m). We consider using a delta modulator. Solution to Question (a): Since d(n) = x(n) ; ax(n ; 1) and E[x(n)] = 0 also E[d(n)] = 0, i.e., the variance d 2 of d(n) is found by evaluating 2 d = E[d 2 (n)] = E h (x(n) ; ax(n ; 1)) 2i = E[x 2 (n)] ; 2aE[x(n)x(n ; 1)] + a 2 E[x 2 (n ; 1)] = 2 x ; 2a xx (1) + a 2 2 x = 2 x 1+a 2 ; 2a xx (1) where we used that the normalized autocorrelation function (Ch ) xx (m) = xx = xx (0) = xx = 2 x. Solution to Question (b): An extremum of 2 d is found by solving That is, the optimal value a is given 2 2 =22 x (a ; xx (1)) = 0, a = xx (1) 2 2 x=@a 2 =2 2 x > 0 the extremum is a minimum. By substituting a we get: 2 d = 2 x 1+ 2 xx(1) ; 2 xx (1) xx (1) = 2 x 1 ; 2 xx(1) Solution to Question (c): 2 d < 2 x is obtained when 1 ; 2 xx(1) < 1, i.e., 2 xx(1) > 0 which always is the case except when x(n) is a white noise signal for which xx (m) = 2 x(m). Solution to Question (d): In order to study the second order prediction error d(n) =x(n) ; a 1 x(n ; 1) ; a 2 x(n ; 2) we dene two column vectors: a =(a 1 a 2 ) > x =(x(n ; 1) x(n ; 2)) > Then we express d(n) =x(n);x > a. The variance of d(n) is found by rst evaluating d 2 (n): d 2 (n) = x(n) ; x > 2 a = x 2 (n) ; 2x(n)x > a + a > xx > a 7
12 Performing the expectation E[] where E[d 2 (n)] = E[x 2 (n)] ; 2E[x(n)x > ]a + a > E[xx > ]a E[x(n)x > ] = (E[x(n)x(n ; 1)] E[x(n)x(n ; 2)]) = ( xx (1) xx (2)) " # " E[xx > E[x(n ; 1)x(n ; 1)] E[x(n ; 1)x(n ; 2)] xx (0) ] = = xx (1) E[x(n ; 1)x(n ; 2)] E[x(n ; 2)x(n ; 2)] xx (1) xx (0) Since xx (0) = 2 x we dene P > = E[x(n)x> ] 2 x R = E[xx> ] 2 x = =( xx (1) xx (2)) " 1 xx (1) xx (1) 1 Notice that R is positive denite, i.e., q > Rq > 0 for all q 6= 0 since q Rq > = E[q > xx > q]=e[(q > x) 2 ] > 0 q 6= 0 2 x Now we have the expression 2 d = 2 x 1 ; 2P > a + a > Ra Extremum is found by 2 d=@a = 0, i.e., # # That is the optimal a, say 2 = 2 x (;2P +2Ra) =0 a = R ;1 P Note that R ;1 exists since it is positive denite. Moreover, the optimum is a minimum since the second order derivative 2 2 d=@a@a > = R is positive denite. By substituting the optimal value a into the expression for 2 d we get 2 d = 2 x 1 ; 2P > R ;1 P + P > R ;1 RR ;1 P = 2 x 1 ; P > R ;1 P P > R ;1 P 0 since R is positive denite and equals zero if P = 0, i.e., when xx (m) =0form>1 (white noise). Thus, except for the white noise case, 2 d < 2 x is always fullled. Further note it is possible to show P > R ;1 P 1 since 2 d 0. 8
13 SOLUTION TO PROBLEM 9.11 Solution to Question (a): in Fig Consider the rewritten second order SDM model shown e(n) x(n) - - H(z) - ; 6 ; H(z) - z ; ? d q (n) Fig where H(z) cf. Fig. P9.11 is dened by H(z) =1=(1;z ;1 ). Using the superposition principle we can write: That is, D q (z) = D q (z) =;D q (z)h(z)z ;1 (1 + H(z)) + H 2 (z)z ;1 (z)+e(z) H 2 (z)z ;1 1+H(z)z ;1 (1 + H(z)) (z)+ 1 1+H(z)z ;1 (1 + H(z)) E(Z) By using the expression for H(z), the noise transfer function is H n (z) = and the signal transfer function is H s (z) = 1 1+H(z)z ;1 (1 + H(z)) =(1; z;1 ) 2 H 2 (z)z ;1 1+H(z)z ;1 (1 + H(z)) = z;1 Solution to Question (b): The magnitude frequency response is found by evaluating at jh n (z)j at z = e j!, i.e., jh n (!)j = 2 1 ; e ;j! = 1 ; 2e ;j! + e ;j2! = e ;j! e j! ; 2+e ;j! = j;2 + 2 cos(!)j 9
14 The noise magnitude responses of a rst and a second order SDM is shown in Fig Noise Frequency Responses H n2 (f) H n1 (f) f Fig The 6 db dierence is due to an extra zero for the 2nd order system. Solution to Question (c): Using the fact that the quantization noise is white, it has the power spectral density S e (F )=e=f 2 s where F s is the sampling frequency and e 2 is the noise power. Using Eq. (9.2.19) where B is the bandwidth of the signal x a (t), we get 2 n = Z B ;B jh n (F=F s )j 2 S e (F ) df Z = 4 2 B e (;1 + cos(2f=f s )) 2 df F s ;B Z 2B=Fs = 22 e (;1 + cos(!)) 2 d! ;2B=F s = 2 " 3! 2 e 2 + sin(2!) # 2B=Fs ; 2sin(!) 4 ;2B=F s = 2 " 6B 2 e + sin(4b=f s) ; 4 sin(2b=f s ) F s 2 Taylor series expansion of sin(x) = x ; x 3 =3! + x 5 =5! for x 1 corresponding to F s B gives: 2 n 4 e 2 2B 5 5 F s A doubling of F s gives a reduction of n 2 by a factor of 32 which corresponds to approx. 15 db. # 10
15 SOLUTION TO PROBLEM 9.12 Given x(n) = cos(2n=n), n =0 1 N; 1. The signal is reconstructed through an ideal D/A converter with sampling interval T corresponding to a sampling frequency F s =1=T. Solution to Question (a): In general, cf. Ch. 4.2, a periodic digital signal has normalized frequency components at f = k=n, k =0 1 N;1, and consequently frequency components at F = kf s =N. A cosine has spectral components for k =1 and k = N ; 1only. Thus, F 0 = F s =N. Solution to Question (b): If F s is xed and we only use the given lookup table, the only possibility isto sample x(n). Downsampling a discrete-time signal with a factor of D is equivalent to x sam (n) =x(n) 1 q=;1 (n ; qd) The spectrum is found by (see e.g., Ch. 10 or the textbooks in signal analysis) sam (f) = 1 D D;1 q=0 (f ; q=d) The spectrum of the signal x(n) considered as a aperiodic signal of length N can be interpreted as the corresponding period signal multiplied with a square wave signal of length N. Thus the spectrum is the convolution of the spectrum of the square wave signal and the cosine signal: (f) = sin nfn sin f e;jf(n;1) ((f ; 1=N)=2+(f +1=N)=2) The spectrum is (1=N) = (;1=N) = 1=2 and has zeros at f = k=n k 6= 1. The sampled signals spectrum is a replication of (f) with q=d. Thus if D is a divisor in N, say N = KD, sam (f) =1=D for f =1=N + q=d =(1+qK)=N, q = 0 1 D;1 and zero for all other multiples of 1=N. The innitely replicated signal is forming a periodic signal with frequencies at k =1 1+K 1+2K 1+(D;1)K and k = ;1+K ;1+2K ;1+DK. Changing to the new sampling frequency 1=D gives components at D=N corresponding to F 0 = DF s =N. The number of possible frequencies is thus the number of all possible products of prime factors in N. If e.g., N =12=2 2 3thepossible D 2f g. If D is not a divisor in N then sam (f) willhave non-zero values for most multiples of 1=N thus the replicated signal is not a cosine. This is shown in the following Figures {4. 11
16 0.6 Original Signal, N= (k) k Fig Sampled Signal (zeros incl.), N=14, D= (k) k Fig
17 0.6 Sampled Signal, N=14, D= (k) k Fig Sampled Signal, N=14, D= (k) k Fig
18 SOLUTION TO PROBLEM 5.20 Solution to Question (a): (k) DFT! N x(n) where N is an even number and x(n) fulll the symmetry property: x(n + N=2) = ;x(n) n =0 1 N=2 ; 1 Recall that W N = e ;j2=n. By evaluating (k) we get: (k) = = = = = N;1 N=2;1 N=2;1 N=2;1 N=2;1 x(n)w kn N x(n)w kn N + N;1 n=n=2 N=2;1 x(n)w kn N + N=2;1 x(n)w kn N + x(n)w kn N x(n + N=2)W k(n+n=2) N ;x(n)w kn=2 N W kn N h x(n)+(;1) k+1 x(n) i W kn N Here we used the fact that W kn=2 N = e ;jk = (;1) k. Now for k = 2k 0, k 0 = 0 1 N=2 ; 1 we have: x(n)+(;1) k+1 x(n) 0, thus (2k 0 )=0. Solution to Question (b): Evaluating the odd harmonics we get: (2k 0 +1) = = = N=2;1 N=2;1 N=2;1 h x(n)+(;1) 2k 0 +2 x(n) i W (2k0 +1)n N [2x(n)W n N] W 2k0 N [2x(n)W n N] W k0 N=2 That is, (2k 0 +1) DFT! N=2 y(n) where y(n) =2x(n)W n N, n =0 1 N=2 ; 1. 14
Discrete-Time Signals and Systems. Efficient Computation of the DFT: FFT Algorithms. Analog-to-Digital Conversion. Sampling Process.
iscrete-time Signals and Systems Efficient Computation of the FT: FFT Algorithms r. eepa Kundur University of Toronto Reference: Sections 6.1, 6., 6.4, 6.5 of John G. Proakis and imitris G. Manolakis,
More informationAPPENDIX A. The Fourier integral theorem
APPENDIX A The Fourier integral theorem In equation (1.7) of Section 1.3 we gave a description of a signal defined on an infinite range in the form of a double integral, with no explanation as to how that
More informationExamples. 2-input, 1-output discrete-time systems: 1-input, 1-output discrete-time systems:
Discrete-Time s - I Time-Domain Representation CHAPTER 4 These lecture slides are based on "Digital Signal Processing: A Computer-Based Approach, 4th ed." textbook by S.K. Mitra and its instructor materials.
More informationTopic 3: Fourier Series (FS)
ELEC264: Signals And Systems Topic 3: Fourier Series (FS) o o o o Introduction to frequency analysis of signals CT FS Fourier series of CT periodic signals Signal Symmetry and CT Fourier Series Properties
More informationInterchange of Filtering and Downsampling/Upsampling
Interchange of Filtering and Downsampling/Upsampling Downsampling and upsampling are linear systems, but not LTI systems. They cannot be implemented by difference equations, and so we cannot apply z-transform
More informationChapter 2: Problem Solutions
Chapter 2: Problem Solutions Discrete Time Processing of Continuous Time Signals Sampling à Problem 2.1. Problem: Consider a sinusoidal signal and let us sample it at a frequency F s 2kHz. xt 3cos1000t
More informationWiener Filtering. EE264: Lecture 12
EE264: Lecture 2 Wiener Filtering In this lecture we will take a different view of filtering. Previously, we have depended on frequency-domain specifications to make some sort of LP/ BP/ HP/ BS filter,
More informationUNIVERSITY OF OSLO. Faculty of mathematics and natural sciences. Forslag til fasit, versjon-01: Problem 1 Signals and systems.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 1th, 016 Examination hours: 14:30 18.30 This problem set
More informationEE 521: Instrumentation and Measurements
Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA September 23, 2009 1 / 18 1 Sampling 2 Quantization 3 Digital-to-Analog Converter 4 Analog-to-Digital Converter
More informationEEM 409. Random Signals. Problem Set-2: (Power Spectral Density, LTI Systems with Random Inputs) Problem 1: Problem 2:
EEM 409 Random Signals Problem Set-2: (Power Spectral Density, LTI Systems with Random Inputs) Problem 1: Consider a random process of the form = + Problem 2: X(t) = b cos(2π t + ), where b is a constant,
More information3.2 Complex Sinusoids and Frequency Response of LTI Systems
3. Introduction. A signal can be represented as a weighted superposition of complex sinusoids. x(t) or x[n]. LTI system: LTI System Output = A weighted superposition of the system response to each complex
More informationEE-210. Signals and Systems Homework 7 Solutions
EE-20. Signals and Systems Homework 7 Solutions Spring 200 Exercise Due Date th May. Problems Q Let H be the causal system described by the difference equation w[n] = 7 w[n ] 2 2 w[n 2] + x[n ] x[n 2]
More informationENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 2: Fourier Representations Jie Liang School of Engineering Science Simon Fraser University 1 Outline Chap 2.1 2.5: Signal Classifications Fourier Transform Dirac Delta Function
More informationDigital Signal Processing, Homework 1, Spring 2013, Prof. C.D. Chung
Digital Signal Processing, Homework, Spring 203, Prof. C.D. Chung. (0.5%) Page 99, Problem 2.2 (a) The impulse response h [n] of an LTI system is known to be zero, except in the interval N 0 n N. The input
More information! Introduction. ! Discrete Time Signals & Systems. ! Z-Transform. ! Inverse Z-Transform. ! Sampling of Continuous Time Signals
ESE 531: Digital Signal Processing Lec 25: April 24, 2018 Review Course Content! Introduction! Discrete Time Signals & Systems! Discrete Time Fourier Transform! Z-Transform! Inverse Z-Transform! Sampling
More informationEE538 Final Exam Fall :20 pm -5:20 pm PHYS 223 Dec. 17, Cover Sheet
EE538 Final Exam Fall 005 3:0 pm -5:0 pm PHYS 3 Dec. 17, 005 Cover Sheet Test Duration: 10 minutes. Open Book but Closed Notes. Calculators ARE allowed!! This test contains five problems. Each of the five
More informationMultimedia Signals and Systems - Audio and Video. Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2
Multimedia Signals and Systems - Audio and Video Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2 Kunio Takaya Electrical and Computer Engineering University of Saskatchewan December
More informationRadar Systems Engineering Lecture 3 Review of Signals, Systems and Digital Signal Processing
Radar Systems Engineering Lecture Review of Signals, Systems and Digital Signal Processing Dr. Robert M. O Donnell Guest Lecturer Radar Systems Course Review Signals, Systems & DSP // Block Diagram of
More informationAdvanced Digital Signal Processing -Introduction
Advanced Digital Signal Processing -Introduction LECTURE-2 1 AP9211- ADVANCED DIGITAL SIGNAL PROCESSING UNIT I DISCRETE RANDOM SIGNAL PROCESSING Discrete Random Processes- Ensemble Averages, Stationary
More informationDiscrete-Time Signals and Systems. The z-transform and Its Application. The Direct z-transform. Region of Convergence. Reference: Sections
Discrete-Time Signals and Systems The z-transform and Its Application Dr. Deepa Kundur University of Toronto Reference: Sections 3. - 3.4 of John G. Proakis and Dimitris G. Manolakis, Digital Signal Processing:
More information7.1 Sampling and Reconstruction
Haberlesme Sistemlerine Giris (ELE 361) 6 Agustos 2017 TOBB Ekonomi ve Teknoloji Universitesi, Guz 2017-18 Dr. A. Melda Yuksel Turgut & Tolga Girici Lecture Notes Chapter 7 Analog to Digital Conversion
More informationENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM. Dr. Lim Chee Chin
ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM Dr. Lim Chee Chin Outline Introduction Discrete Fourier Series Properties of Discrete Fourier Series Time domain aliasing due to frequency
More informationImage Acquisition and Sampling Theory
Image Acquisition and Sampling Theory Electromagnetic Spectrum The wavelength required to see an object must be the same size of smaller than the object 2 Image Sensors 3 Sensor Strips 4 Digital Image
More informationGeneralized Sidelobe Canceller and MVDR Power Spectrum Estimation. Bhaskar D Rao University of California, San Diego
Generalized Sidelobe Canceller and MVDR Power Spectrum Estimation Bhaskar D Rao University of California, San Diego Email: brao@ucsd.edu Reference Books 1. Optimum Array Processing, H. L. Van Trees 2.
More informationDiscrete Fourier transform (DFT)
Discrete Fourier transform (DFT) Signal Processing 2008/9 LEA Instituto Superior Técnico Signal Processing LEA (IST) Discrete Fourier transform 1 / 34 Periodic signals Consider a periodic signal x[n] with
More informationINTRODUCTION TO DELTA-SIGMA ADCS
ECE37 Advanced Analog Circuits INTRODUCTION TO DELTA-SIGMA ADCS Richard Schreier richard.schreier@analog.com NLCOTD: Level Translator VDD > VDD2, e.g. 3-V logic? -V logic VDD < VDD2, e.g. -V logic? 3-V
More informationContents. Digital Signal Processing, Part II: Power Spectrum Estimation
Contents Digital Signal Processing, Part II: Power Spectrum Estimation 5. Application of the FFT for 7. Parametric Spectrum Est. Filtering and Spectrum Estimation 7.1 ARMA-Models 5.1 Fast Convolution 7.2
More informationThe Discrete Fourier Transform (DFT) Properties of the DFT DFT-Specic Properties Power spectrum estimate. Alex Sheremet.
4. April 2, 27 -order sequences Measurements produce sequences of numbers Measurement purpose: characterize a stochastic process. Example: Process: water surface elevation as a function of time Parameters:
More information[ ], [ ] [ ] [ ] = [ ] [ ] [ ]{ [ 1] [ 2]
4. he discrete Fourier transform (DF). Application goal We study the discrete Fourier transform (DF) and its applications: spectral analysis and linear operations as convolution and correlation. We use
More informationDISCRETE-TIME SIGNAL PROCESSING
THIRD EDITION DISCRETE-TIME SIGNAL PROCESSING ALAN V. OPPENHEIM MASSACHUSETTS INSTITUTE OF TECHNOLOGY RONALD W. SCHÄFER HEWLETT-PACKARD LABORATORIES Upper Saddle River Boston Columbus San Francisco New
More informationElec4621 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis
Elec461 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis Dr. D. S. Taubman May 3, 011 In this last chapter of your notes, we are interested in the problem of nding the instantaneous
More informationEE 3054: Signals, Systems, and Transforms Spring A causal discrete-time LTI system is described by the equation. y(n) = 1 4.
EE : Signals, Systems, and Transforms Spring 7. A causal discrete-time LTI system is described by the equation Test y(n) = X x(n k) k= No notes, closed book. Show your work. Simplify your answers.. A discrete-time
More informationMEDE2500 Tutorial Nov-7
(updated 2016-Nov-4,7:40pm) MEDE2500 (2016-2017) Tutorial 3 MEDE2500 Tutorial 3 2016-Nov-7 Content 1. The Dirac Delta Function, singularity functions, even and odd functions 2. The sampling process and
More informationDigital Speech Processing Lecture 10. Short-Time Fourier Analysis Methods - Filter Bank Design
Digital Speech Processing Lecture Short-Time Fourier Analysis Methods - Filter Bank Design Review of STFT j j ˆ m ˆ. X e x[ mw ] [ nˆ m] e nˆ function of nˆ looks like a time sequence function of ˆ looks
More informationOversampling Converters
Oversampling Converters David Johns and Ken Martin (johns@eecg.toronto.edu) (martin@eecg.toronto.edu) slide 1 of 56 Motivation Popular approach for medium-to-low speed A/D and D/A applications requiring
More informationSignals, Instruments, and Systems W5. Introduction to Signal Processing Sampling, Reconstruction, and Filters
Signals, Instruments, and Systems W5 Introduction to Signal Processing Sampling, Reconstruction, and Filters Acknowledgments Recapitulation of Key Concepts from the Last Lecture Dirac delta function (
More informationDIGITAL SIGNAL PROCESSING UNIT 1 SIGNALS AND SYSTEMS 1. What is a continuous and discrete time signal? Continuous time signal: A signal x(t) is said to be continuous if it is defined for all time t. Continuous
More informationChapter 7: Filter Design 7.1 Practical Filter Terminology
hapter 7: Filter Design 7. Practical Filter Terminology Analog and digital filters and their designs constitute one of the major emphasis areas in signal processing and communication systems. This is due
More informationDiscrete-Time Signals & Systems
Chapter 2 Discrete-Time Signals & Systems 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 2-1-1 Discrete-Time Signals: Time-Domain Representation (1/10) Signals
More informationDIGITAL SIGNAL PROCESSING LECTURE 1
DIGITAL SIGNAL PROCESSING LECTURE 1 Fall 2010 2K8-5 th Semester Tahir Muhammad tmuhammad_07@yahoo.com Content and Figures are from Discrete-Time Signal Processing, 2e by Oppenheim, Shafer, and Buck, 1999-2000
More informationNORWEGIAN UNIVERSITY OF SCIENCE AND TECHNOLOGY DEPARTMENT OF ELECTRONICS AND TELECOMMUNICATIONS
page 1 of 5 (+ appendix) NORWEGIAN UNIVERSITY OF SCIENCE AND TECHNOLOGY DEPARTMENT OF ELECTRONICS AND TELECOMMUNICATIONS Contact during examination: Name: Magne H. Johnsen Tel.: 73 59 26 78/930 25 534
More informationSolutions. Number of Problems: 10
Final Exam February 9th, 2 Signals & Systems (5-575-) Prof. R. D Andrea Solutions Exam Duration: 5 minutes Number of Problems: Permitted aids: One double-sided A4 sheet. Questions can be answered in English
More informationECE503: Digital Signal Processing Lecture 5
ECE53: Digital Signal Processing Lecture 5 D. Richard Brown III WPI 3-February-22 WPI D. Richard Brown III 3-February-22 / 32 Lecture 5 Topics. Magnitude and phase characterization of transfer functions
More informationE : Lecture 1 Introduction
E85.2607: Lecture 1 Introduction 1 Administrivia 2 DSP review 3 Fun with Matlab E85.2607: Lecture 1 Introduction 2010-01-21 1 / 24 Course overview Advanced Digital Signal Theory Design, analysis, and implementation
More informationButterworth Filter Properties
OpenStax-CNX module: m693 Butterworth Filter Properties C. Sidney Burrus This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 3. This section develops the properties
More informationProbability Space. J. McNames Portland State University ECE 538/638 Stochastic Signals Ver
Stochastic Signals Overview Definitions Second order statistics Stationarity and ergodicity Random signal variability Power spectral density Linear systems with stationary inputs Random signal memory Correlation
More informationDigital Signal Processing
Digital Signal Processing Introduction Moslem Amiri, Václav Přenosil Embedded Systems Laboratory Faculty of Informatics, Masaryk University Brno, Czech Republic amiri@mail.muni.cz prenosil@fi.muni.cz February
More informationDetailed Solutions to Exercises
Detailed Solutions to Exercises Digital Signal Processing Mikael Swartling Nedelko Grbic rev. 205 Department of Electrical and Information Technology Lund University Detailed solution to problem E3.4 A
More informationDiscrete-time Fourier Series (DTFS)
Discrete-time Fourier Series (DTFS) Arun K. Tangirala (IIT Madras) Applied Time-Series Analysis 59 Opening remarks The Fourier series representation for discrete-time signals has some similarities with
More informationCITY UNIVERSITY LONDON. MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746
No: CITY UNIVERSITY LONDON MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746 Date: 19 May 2004 Time: 09:00-11:00 Attempt Three out of FIVE questions, at least One question from PART B PART
More informationLecture 4 - Spectral Estimation
Lecture 4 - Spectral Estimation The Discrete Fourier Transform The Discrete Fourier Transform (DFT) is the equivalent of the continuous Fourier Transform for signals known only at N instants separated
More informationUNIT 1. SIGNALS AND SYSTEM
Page no: 1 UNIT 1. SIGNALS AND SYSTEM INTRODUCTION A SIGNAL is defined as any physical quantity that changes with time, distance, speed, position, pressure, temperature or some other quantity. A SIGNAL
More informationCh.11 The Discrete-Time Fourier Transform (DTFT)
EE2S11 Signals and Systems, part 2 Ch.11 The Discrete-Time Fourier Transform (DTFT Contents definition of the DTFT relation to the -transform, region of convergence, stability frequency plots convolution
More information2A1H Time-Frequency Analysis II
2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 209 For any corrections see the course page DW Murray at www.robots.ox.ac.uk/ dwm/courses/2tf. (a) A signal g(t) with period
More informationUniversity Question Paper Solution
Unit 1: Introduction University Question Paper Solution 1. Determine whether the following systems are: i) Memoryless, ii) Stable iii) Causal iv) Linear and v) Time-invariant. i) y(n)= nx(n) ii) y(t)=
More informationSolutions. Number of Problems: 10
Final Exam February 2nd, 2013 Signals & Systems (151-0575-01) Prof. R. D Andrea Solutions Exam Duration: 150 minutes Number of Problems: 10 Permitted aids: One double-sided A4 sheet. Questions can be answered
More informationFundamentals of the DFT (fft) Algorithms
Fundamentals of the DFT (fft) Algorithms D. Sundararajan November 6, 9 Contents 1 The PM DIF DFT Algorithm 1.1 Half-wave symmetry of periodic waveforms.............. 1. The DFT definition and the half-wave
More informationEE538 Final Exam Fall 2007 Mon, Dec 10, 8-10 am RHPH 127 Dec. 10, Cover Sheet
EE538 Final Exam Fall 2007 Mon, Dec 10, 8-10 am RHPH 127 Dec. 10, 2007 Cover Sheet Test Duration: 120 minutes. Open Book but Closed Notes. Calculators allowed!! This test contains five problems. Each of
More informationECE-700 Review. Phil Schniter. January 5, x c (t)e jωt dt, x[n]z n, Denoting a transform pair by x[n] X(z), some useful properties are
ECE-7 Review Phil Schniter January 5, 7 ransforms Using x c (t) to denote a continuous-time signal at time t R, Laplace ransform: X c (s) x c (t)e st dt, s C Continuous-ime Fourier ransform (CF): ote that:
More informationECE 673-Random signal analysis I Final
ECE 673-Random signal analysis I Final Q ( point) Comment on the following statement "If cov(x ; X 2 ) 0; then the best predictor (without limitation on the type of predictor) of the random variable X
More informationA.1 THE SAMPLED TIME DOMAIN AND THE Z TRANSFORM. 0 δ(t)dt = 1, (A.1) δ(t)dt =
APPENDIX A THE Z TRANSFORM One of the most useful techniques in engineering or scientific analysis is transforming a problem from the time domain to the frequency domain ( 3). Using a Fourier or Laplace
More informationELEG-636: Statistical Signal Processing
ELEG-636: Statistical Signal Processing Gonzalo R. Arce Department of Electrical and Computer Engineering University of Delaware Spring 2010 Gonzalo R. Arce (ECE, Univ. of Delaware) ELEG-636: Statistical
More informationFilter Banks II. Prof. Dr.-Ing. G. Schuller. Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany
Filter Banks II Prof. Dr.-Ing. G. Schuller Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany Page Modulated Filter Banks Extending the DCT The DCT IV transform can be seen as modulated
More informationDown-Sampling (4B) Young Won Lim 11/15/12
Down-Sampling (B) /5/ Copyright (c) 9,, Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version. or any later
More informationOptimization using Calculus. Optimization of Functions of Multiple Variables subject to Equality Constraints
Optimization using Calculus Optimization of Functions of Multiple Variables subject to Equality Constraints 1 Objectives Optimization of functions of multiple variables subjected to equality constraints
More information8 The Discrete Fourier Transform (DFT)
8 The Discrete Fourier Transform (DFT) ² Discrete-Time Fourier Transform and Z-transform are de ned over in niteduration sequence. Both transforms are functions of continuous variables (ω and z). For nite-duration
More informationCHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME
CHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME Shri Mata Vaishno Devi University, (SMVDU), 2013 Page 13 CHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME When characterizing or modeling a random variable, estimates
More informationHST.582J / 6.555J / J Biomedical Signal and Image Processing Spring 2007
MIT OpenCourseWare http://ocw.mit.edu HST.58J / 6.555J / 16.456J Biomedical Signal and Image Processing Spring 007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationAnalog to Digital Converters (ADCs)
Analog to Digital Converters (ADCs) Note: Figures are copyrighted Proakis & Manolakis, Digital Signal Processing, 4 th Edition, Pearson Publishers. Embedded System Design A Unified HW Approach, Vahid/Givargis,
More informationMagnitude F y. F x. Magnitude
Design of Optimum Multi-Dimensional Energy Compaction Filters N. Damera-Venkata J. Tuqan B. L. Evans Imaging Technology Dept. Dept. of ECE Dept. of ECE Hewlett-Packard Labs Univ. of California Univ. of
More informationEDISP (NWL2) (English) Digital Signal Processing Transform, FT, DFT. March 11, 2015
EDISP (NWL2) (English) Digital Signal Processing Transform, FT, DFT March 11, 2015 Transform concept We want to analyze the signal represent it as built of some building blocks (well known signals), possibly
More informationSEL4223 Digital Signal Processing. Inverse Z-Transform. Musa Mohd Mokji
SEL4223 Digital Signal Processing Inverse Z-Transform Musa Mohd Mokji Inverse Z-Transform Transform from z-domain to time-domain x n = 1 2πj c X z z n 1 dz Note that the mathematical operation for the
More informationPractical Spectral Estimation
Digital Signal Processing/F.G. Meyer Lecture 4 Copyright 2015 François G. Meyer. All Rights Reserved. Practical Spectral Estimation 1 Introduction The goal of spectral estimation is to estimate how the
More information/ (2π) X(e jω ) dω. 4. An 8 point sequence is given by x(n) = {2,2,2,2,1,1,1,1}. Compute 8 point DFT of x(n) by
Code No: RR320402 Set No. 1 III B.Tech II Semester Regular Examinations, Apr/May 2006 DIGITAL SIGNAL PROCESSING ( Common to Electronics & Communication Engineering, Electronics & Instrumentation Engineering,
More informationDigital Signal Processing Course. Nathalie Thomas Master SATCOM
Digital Signal Processing Course Nathalie Thomas Master SATCOM 2018 2019 Chapitre 1 Introduction 1.1 Signal digitization : sampling and quantization 1.1.1 Signal sampling Let's call x(t) the original deterministic
More informationTheory and Problems of Signals and Systems
SCHAUM'S OUTLINES OF Theory and Problems of Signals and Systems HWEI P. HSU is Professor of Electrical Engineering at Fairleigh Dickinson University. He received his B.S. from National Taiwan University
More informationLecture 19: Discrete Fourier Series
EE518 Digital Signal Processing University of Washington Autumn 2001 Dept. of Electrical Engineering Lecture 19: Discrete Fourier Series Dec 5, 2001 Prof: J. Bilmes TA: Mingzhou
More informationETSF15 Analog/Digital. Stefan Höst
ETSF15 Analog/Digital Stefan Höst Physical layer Analog vs digital Sampling, quantisation, reconstruction Modulation Represent digital data in a continuous world Disturbances Noise and distortion Synchronization
More informationENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University
ENSC37 Communications Systems : Fourier Representations School o Engineering Science Simon Fraser University Outline Chap..5: Signal Classiications Fourier Transorm Dirac Delta Function Unit Impulse Fourier
More informationExamination with solution suggestions SSY130 Applied Signal Processing
Examination with solution suggestions SSY3 Applied Signal Processing Jan 8, 28 Rules Allowed aids at exam: L. Råde and B. Westergren, Mathematics Handbook (any edition, including the old editions called
More informationChapter 4 Discrete Fourier Transform (DFT) And Signal Spectrum
Chapter 4 Discrete Fourier Transform (DFT) And Signal Spectrum CEN352, DR. Nassim Ammour, King Saud University 1 Fourier Transform History Born 21 March 1768 ( Auxerre ). Died 16 May 1830 ( Paris ) French
More informationECGR4124 Digital Signal Processing Exam 2 Spring 2017
ECGR4124 Digital Signal Processing Exam 2 Spring 2017 Name: LAST 4 NUMBERS of Student Number: Do NOT begin until told to do so Make sure that you have all pages before starting NO TEXTBOOK, NO CALCULATOR,
More informationECE-S Introduction to Digital Signal Processing Lecture 3C Properties of Autocorrelation and Correlation
ECE-S352-701 Introduction to Digital Signal Processing Lecture 3C Properties of Autocorrelation and Correlation Assuming we have two sequences x(n) and y(n) from which we form a linear combination: c(n)
More informationOn the Frequency-Domain Properties of Savitzky-Golay Filters
On the Frequency-Domain Properties of Savitzky-Golay Filters Ronald W Schafer HP Laboratories HPL-2-9 Keyword(s): Savitzky-Golay filter, least-squares polynomial approximation, smoothing Abstract: This
More informationDigital Signal Processing Module 1 Analysis of Discrete time Linear Time - Invariant Systems
Digital Signal Processing Module 1 Analysis of Discrete time Linear Time - Invariant Systems Objective: 1. To understand the representation of Discrete time signals 2. To analyze the causality and stability
More informationEach problem is worth 25 points, and you may solve the problems in any order.
EE 120: Signals & Systems Department of Electrical Engineering and Computer Sciences University of California, Berkeley Midterm Exam #2 April 11, 2016, 2:10-4:00pm Instructions: There are four questions
More informationOther Optimization Techniques. Similar to steepest descent, but slightly different way of choosing direction of next step: x r +1 = x.
Conjugate Gradient Other Optimization Techniques Similar to steepest descent, but slightly different way of choosing direction of next step: x r +1 = x r + r s r s 0 = g 0 s r +1 = g r +1 + r +1 s r {
More informationsummable Necessary and sufficient for BIBO stability of an LTI system. Also see poles.
EECS 206 DSP GLOSSARY c Andrew E. Yagle Fall 2005 absolutely impulse response: h[n] is finite. EX: n=0 ( 3 4 )n = 1 = 4 but 1 3 n=1 1 n. 4 summable Necessary and sufficient for BIBO stability of an LI
More informationHilbert Transforms in Signal Processing
Hilbert Transforms in Signal Processing Stefan L. Hahn Artech House Boston London Contents Preface xiii Introduction 1 Chapter 1 Theory of the One-Dimensional Hilbert Transformation 3 1.1 The Concepts
More informationω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the
he ime-frequency Concept []. Review of Fourier Series Consider the following set of time functions {3A sin t, A sin t}. We can represent these functions in different ways by plotting the amplitude versus
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science. Fall Solutions for Problem Set 2
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Issued: Tuesday, September 5. 6.: Discrete-Time Signal Processing Fall 5 Solutions for Problem Set Problem.
More informationEE 521: Instrumentation and Measurements
Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA November 1, 2009 1 / 27 1 The z-transform 2 Linear Time-Invariant System 3 Filter Design IIR Filters FIR Filters
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 17, 2017 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More informationLecture 7 Discrete Systems
Lecture 7 Discrete Systems EE 52: Instrumentation and Measurements Lecture Notes Update on November, 29 Aly El-Osery, Electrical Engineering Dept., New Mexico Tech 7. Contents The z-transform 2 Linear
More informationHigher-Order Σ Modulators and the Σ Toolbox
ECE37 Advanced Analog Circuits Higher-Order Σ Modulators and the Σ Toolbox Richard Schreier richard.schreier@analog.com NLCOTD: Dynamic Flip-Flop Standard CMOS version D CK Q Q Can the circuit be simplified?
More informationStability Condition in Terms of the Pole Locations
Stability Condition in Terms of the Pole Locations A causal LTI digital filter is BIBO stable if and only if its impulse response h[n] is absolutely summable, i.e., 1 = S h [ n] < n= We now develop a stability
More informationFilter Banks with Variable System Delay. Georgia Institute of Technology. Atlanta, GA Abstract
A General Formulation for Modulated Perfect Reconstruction Filter Banks with Variable System Delay Gerald Schuller and Mark J T Smith Digital Signal Processing Laboratory School of Electrical Engineering
More informationPS403 - Digital Signal processing
PS403 - Digital Signal processing III. DSP - Digital Fourier Series and Transforms Key Text: Digital Signal Processing with Computer Applications (2 nd Ed.) Paul A Lynn and Wolfgang Fuerst, (Publisher:
More informationUNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 9th, 011 Examination hours: 14.30 18.30 This problem set
More informationECE 8440 Unit 13 Sec0on Effects of Round- Off Noise in Digital Filters
ECE 8440 Unit 13 Sec0on 6.9 - Effects of Round- Off Noise in Digital Filters 1 We have already seen that if a wide- sense staonary random signal x(n) is applied as input to a LTI system, the power density
More information