b) discrete-time iv) aperiodic (finite energy)
|
|
- Liliana Hamilton
- 5 years ago
- Views:
Transcription
1 EE 464 Frequency Analysis of Signals and Systems Fall 2018 Read Text, Chapter. Study suggestion: Use Matlab to plot several of the signals and their DTFT in the examples to follow. Some types of signal and Fourier analysis a) continuous-time (analog) i) periodic (infinite energy) ii) aperiodic (finite energy) Fourier series Fourier transform xx(tt) = cc nn ee jj2ππππππ TT nn= XX(ff) = xx(tt)ee jj2ππππππ dddd TT 0 cc nn = 1 TT xx(tt)ee jj2ππππππ TT dddd xx(tt) = XX(ff)ee jj2ππππππ dddd b) discrete-time iii) periodic discrete-time Fourier series iv) aperiodic (finite energy) discrete-time Fourier transform 1
2 Discrete-Time Fourier Transform (DTFT) XX(ff) = xx(nn)ee jj2ππππππtt nn= = XX ee jj2ππππtt Here, ff is frequency in Hertz, and the DTFT is periodic in ff with period 1 TT. Note: In lecture, homework, notes, tests, etc., the notation XX ee jj2ππππtt will usually be used (to clearly distinguish between the Fourier transform of an analog signal and the DTFT of a discrete-time signal). The inverse DTFT is given by 1 2TT xx(nn) = TT XX ee jj2ππππππ ee jj2ππππππππ dddd 1 2TT (Set T = 1 in the expression if using normalized frequency.) 2
3 Four expressions for DTFT all equivalent, but using different frequency variables. XX ee jjjj = xx(nn)ee jjjjjj, (1) nn= ωω is in units of normalized frequency, radians/sample, and XX ee jjjj is periodic 2ππ in ωω. XX ee jjjjjj = xx(nn)ee jjjjjjjj, (2) nn= ωω in units of frequency, radians/sec, and XX ee jjjjjj is periodic 2ππ TT in ωω. XX ee jj2ππππ = xx(nn)ee jj2ππππππ, (3) nn= ff is in units of normalized frequency, cycles/sample, and XX ee jj2ππππ is periodic 1 in ff. XX ee jj2ππππππ = xx(nn)ee jj2ππππππππ, nn= (4) ff is in units of frequency, cycles/sec (Hertz), and XX ee jj2ππππππ is periodic 1 in ff. TT 3
4 Example 1. Find the DTFT of xx(nn) = αα nn uu(nn), αα < 1. Plot XX ee jj2ππππtt for real-valued 0 < αα < 1. Modify to units of rad/s. Relate to the Fourier transform of an analog signal xx aa (tt) = ee ββββ uu(tt), sampled at the rate 1 samples/sec. Plot XX TT aa(ff) and relate to XX ee jj2ππππππ. 4
5 Specific case: xx aa (tt) = ee 2tt uu(tt) FF XX aa (ff) = 1 2+jj2ππππ. Select the sampling rate as 1 TT = 8 Hz. Then ee 2TT = αα 0.78 The DTFT of xx(nn) = xx aa (nnnn) = ee 2nnnn uu(nnnn) = αα nn uu(nn) is XX ee jj2ππππππ = XX ee jj2ππππππ = 1 ααee 1 jj2ππππππ with 1 1 2αα cos(2ππππππ)+αα 2. Note that XX ee jj2ππππππ is periodic with period 1 = 8 Hz. The TT difference between the two curves is due to aliasing. Note also the scaling of XX aa (ff) by 1 in the figure. TT 5
6 Low Frequency and High Frequency Since the DTFT is periodic in frequency, normally low frequency is zero, and high frequency is one-half the period (that is, one-half the sampling frequency). Example. Assume a normalized frequency (cycles/sample, or equivalently, TT = 1). For xx(nn) = αα nn uu(nn), αα < 1, the DTFT was found to be XX ee jj2ππππ = 1 ααee 1 jj2ππππ for the normalized frequency variable ff in units of cycles/sample. Then high frequency is 1 2 cycles/sample. Alternatively, using the frequency variable ff in units of Hz, then XX ee jj2ππππππ 1 1 = 1 ααee jj2ππππππ, and high frequency is Hz. 2TT 6
7 Existence of the DTFT. Since the DTFT is a special case of the ZZ-transform, where zz is on the unit circle, the DTFT exists whenever the ZZtransform ROC contains the unit circle. All BIBO stable LTI systems have ROC containing the unit circle, implying the DTFT of the impulse response must exist. HH ee jj2ππππππ = DTFT{h(nn)} is the frequency response of a stable LTI system. 7
8 Example. A LTI system has impulse response h(nn) = 1 2, 1 2. Determine the transfer function, HH(zz), and the DTFT HH ee jj2ππππtt. Plot HH ee jj2ππππtt. Relate the dc gain and the high-frequency gain to the impulse response and to HH(zz). Repeat for a LTI system with impulse response gg(nn) = 1 2,
9 Useful Facts and Properties Let xx(nn) have DTFT XX ee jj2ππππππ = Write in magnitude/phase form: XX ee jj2ππππππ = XX ee jj2ππππππ ee jjjj(ff) where θθ(ff) = tan 1 IIII XX eejj2ππππππ RRRR XX ee jj2ππππππ. nn= xx(nn)ee jj2ππππππππ. 1. The only causal signal with θθ(ff) = 0 for all ff (that is, zero phase) is xx(nn) = AAAA(nn). Hence, the only zerophase causal filter is the simple scaling filter. That is, h(nn) = AAAA(nn), with a system input-output description is yy(nn) = AAAA(nn). This is also a (trivial) all-pass filter. 9
10 2. A (non-causal) real-valued signal can have zero phase if xx(nn) = xx( nn) for all nn. Proof. Write out XX ee jj2ππππππ : XX ee jj2ππππππ = xx(nn)ee jj2ππππππππ nn= = xx(0) + xx(1)ee jj2ππππ1tt + xx( 1)ee jj2ππππ( 1)TT + xx(2)ee jj2ππππ2tt + xx( 2)ee jj2ππππ( 2)TT + XX ee jj2ππππππ XX ee jj2ππππππ Hence, the phase of XX ee jj2ππππππ is = xx(0) + xx(1) ee jj2ππππππ + ee jj2ππππππ + xx(2) ee jj2ππππ2tt + ee jj2ππππ2tt + = xx(0) + xx(1)[2 cos(2ππππππ)] + xx(2)[2 cos(4ππππππ)] + = real valued function of ff θθ(ff) = tan 1 Im XX eejj2ππππππ Re{XX(ee jj2ππππππ )} = tan 1 0 Re{XX(ee jj2ππππππ = 0, )} 10
11 provided Re XX ee jj2ππππππ > 0. So, a symmetric signal, xx(nn) = xx( nn) implies zero phase. Example. The signal h(nn) =,0, 1,2, 6, 2, 1,0, has zero phase. The DTFT is evaluated as HH ee jj2ππππππ = cos(2ππππππ) 2 cos(4ππππππ) Setting TT = 1 so frequency is normalized to cycles/sample, the DTFT is plotted below. 11
12 Note: The period is 1, and high frequency is ff = 1 2 cycles/sample. The phase is zero for all ff. 3. DTFT Shifting Property: If xx(nn) DDDDDDDD XX ee jj2ππππππ, then xx(nn nn 0 ) DDDDDDDD ee jj2ππππππnn 0XX ee jj2ππππππ. Proof: Use the shifting property of the Z-transform with zz = ee jj2ππππππ. Linear phase is of the form θθ(ff) = AAff, for some constant AA, so a linear phase DTFT has the form HH ee jj2ππππππ = HH ee jj2ππππππ ee jjjjff. 4. A causal FIR filter with impulse response h(0),, h(nn 1) has linear phase if h(nn) = h(nn 1 nn), nn = 0,1,, NN 1. Proof. The DTFT is NN 1 HH ee jj2ππππππ = h(nn)ee jj2ππππππππ nn=0 = h(0) + h(1)ee jj2ππππ1tt + h(2)ee jj2ππππ2tt + + h(nn 2)ee jj2ππππ(nn 2)TT + h(nn 1)ee jj2ππππ(nn 1)TT 12
13 HH ee jj2ππππππ = h(0) 1 + ee jj2ππππ(nn 1)TT + h(1) ee jj2ππππ1tt + ee jj2ππππ(nn 2)TT + Now, factor out the same factor from each complex exponential, namely ee jj2ππππ(nn 1)TT 2. Then HH ee jj2ππππππ HH ee jj2ππππππ = ee jj2ππππ(nn 1)TT 2 h(0) ee jj2ππππ(nn 1)TT 2 + ee jj2ππππ(nn 1)TT 2 + h(1) ee jj2ππππ(nn 3)TT 2 + ee jj2ππππ(nn 3)TT 2 + = ee jj2ππππ(nn 1)TT 2ππππ(NN 1)TT 2 h(0) cos 2 2ππππ(NN 3)TT + h(1) cos + 2 Hence, the phase of HH ee jj2ππππππ is 13 θθ(ff) = 2ππππ(NN 1)TT. 2
14 Example. A LTI system has impulse response h(nn) = 1 1, 2,6,2, 1. The length is NN = 5. By 8 inspection, h(nn) = h(nn 1 nn), nn = 0,1,2,3,4, so this is a linear-phase filter. What is the dc gain? What is the high-frequency gain? What type of filter is this? Example. A LTI system has impulse response h(nn) = 1 4 1, 2,1. Find the DTFT HH eejj2ππππππ. Use convolution to find the response to input xx(nn) = 0, 1, 2,3,4,3,2,1,0. Relate the filtering time delay to the DTFT. 14
15 5. Convolution property. yy(nn) = xx(nn) h(nn) DDDDDDDD XX ee jj2ππππππ HH ee jj2ππππππ (1) and yy(nn) = xx(nn)h(nn) DDDDDDDD XX ee jj2ππππππ HH ee jj2ππππππ (2) where convolution in the frequency domain is periodic convolution defined as 1 2TT XX ee jj2ππππππ HH ee jj2ππππππ = TT XX ee jj2ππ(ff λλ)tt HH ee jj2ππλλtt dddd 1 2TT Proof of (1) follows directly from the convolution property of the ZZ-transform, with zz = ee jj2ππππππ. Proof of (2). Begin with the right-hand side of (2) and take the inverse DTFT. 15
16 6. Parseval s Relation Proof. EE xx = 1 2TT xx(nn) 2 = TT XX ee jj2ππππππ 2 ddff nn= 1 2TT 16
17 7. Let xx(nn) DDDDDDDD XX ee jj2ππππππ. In general, xx(nn) and XX ee jj2ππππππ are complex-valued. Write the real and imaginary parts as the signal name with subscript RR or II. So xx(nn) = xx RR (nn) + jjxx II (nn) XX ee jj2ππππππ = XX RR ee jj2ππππππ + jjxx II ee jj2ππππππ a. If xx(nn) real-valued then XX RR ee jj2ππππππ is even in frequency ff and XX II ee jj2ππππππ is odd in frequency, so the magnitude, XX ee jj2ππππππ is even, and the phase, XX ee jj2ππππππ is odd in frequency. b. If xx(nn) real-valued and even then XX ee jj2ππππππ is real and even. c. If xx(nn) real-valued and odd then XX ee jj2ππππππ is imaginary and XX II ee jj2ππππππ is odd. d. If xx(nn) imaginary then XX RR ee jj2ππππππ is odd and XX II ee jj2ππππππ is even. 8. Time reversal. If xx(nn) DDDDDDDD XX ee jj2ππππππ then xx( nn) DDDDDDDD XX ee jj2ππππππ Proof. This property follows directly from the ZZ-transform property. 17
18 9. Energy Density Spectrum and Autocorrelation Let xx(nn) be an energy signal with xx(nn) DDDDDDDD XX ee jj2ππππππ. The magnitude squared of the DTFT is referred to as the energy density spectrum, SS xxxx ee jj2ππππππ = XX ee jj2ππππππ 2. From Parseval s relation, the signal energy is 1 2TT EE xx = TT SS xxxx ee jj2ππππππ ddff. 1 2TT Property: The energy density spectrum is related to the autocorrelation as rr xxxx (kk) DDDDDDDD SS xxxx ee jj2ππππππ, where the autocorrelation function is rr xxxx (kk) = xx(nn + kk)xx (nn) Proof. Directly compute the DTFT of rr xxxx (kk). nn=. 18
19 Relate LTI System Transfer Function Poles and Zeros to Filter Frequency Response Suppose we have a causal LTI system described by NN yy(nn) = aa kk yy(nn kk) kk=1 MM + bb kk xx(nn kk). kk=0 with transfer function HH(zz) = BB(zz) MM = kk=0 bb kkzz kk AA(zz) 1 NN aa kk zz kk kk=1 system has poles and zeros, and the numerator and denominator can be factored into the form HH(zz) = BB(zz) AA(zz) = KK ii=1 (1 rr iizz 1 ) NN (1 pp ii zz 1 ) MM ii=1. The where we assume MM zeros at locations zz = rr ii and NN poles at locations zz = pp ii and no pole-zero cancellations. Simple manipulation yields HH(zz) = KKzzNN MM MM ii=1(zz rr ii ) NN. (zz pp ii ) ii=1 The magnitude frequency response is obtained by setting zz = ee jj2ππππππ and evaluating the magnitude of HH ee jj2ππππππ as HH ee jj2ππffff = KK eejj2ππffff(nn MM) ii=1 ee jj2ππffff rr ii NN (ee jj2ππffff pp ii ) ii=1 MM 19
20 In the zz-plane, this is HH ee jj2ππffff = KK NN ii=1 MM ii=1 dddddddd ffrrrrrr eejj2ππffff tttt zzeeeeee rr ii dddddddd ffrrrrrr ee jj2ππffff tttt pppppppp pp ii Now, what about phase? The frequency response can be written as (assuming KK is positive) MM HH ee jj2ππffff = KK eejj2ππππππ rr ii ee jjφφ ii ii=1 NN ee jj2ππππππ pp ii ii=1 ee jjθθ ii where φφ ii = ee jj2ππππππ rr ii is the angle to the zero and θθ ii = ee jj2ππππππ pp ii is the angle to the pole. Combining terms MM HH ee jj2ππffff = KK eejj2ππππππ rr ii ee jj φφ ii ii=1 NN ee jj2ππππππ pp ii ii=1 MM NN ii=1 ii=1 θθ ii So, the phase of HH ee jj2ππffff equals the sum of the phases from the zeros, minus the sum of the phases from the poles.. 20
21 Example. Let HH(zz) = zz 1 = zz zz 0.9. Sketch the pole-zero locations in the z-plane, and then the magnitude and phase frequency response. 21
22 22
23 Example 2. Let HH(zz) = 2 1 zz Sketch the magnitude (1+0.81zz 2 ) frequency response of the filter. First, sketch the pole and zero locations in the zz-plane. 23
24 Comments about frequency response Begin at ff = 0 and traverse the unit circle as ee jj2ππππππ. As ee jj2ππππππ approaches a zero, one term in the numerator of the frequency response decreases rapidly, so HH ee jj2ππffff decreases. As ee jj2ππππππ moves away from a zero, then HH ee jj2ππffff increases. As ee jj2ππππππ approaches a pole, one term in the denominator of the frequency response decreases rapidly, so HH ee jj2ππffff increases. A zero exactly on the unit circle causes the gain to go to zero. For poles or zeros near the unit circle, the phase suddenly changes as ee jj2ππππππ sweeps past. Poles or zeros at the origin do not affect the magnitude frequency response. 24
25 10. Time shifting and frequency shifting. If xx(nn) DDDDDDDD XX ee jj2ππππππ, then a. xx(nn nn 0 ) DDDDDDDD ee jj2ππππππnn 0XX ee jj2ππππππ b. ee jj2ππff 0nnnn xx(nn) DDDDDDDD XX ee jj(2ππππ 2ππff 0 )TT Proof. Direct computation of the DTFT and inverse DTFT. Example. We already know that XX ee jj2ππππππ is periodic in ff with period 1. So, if ff TT 0 = kk for any integer kk, then TT XX ee jj(2ππππ 2ππff 0 )TT = XX ee jj2ππππππ. Example. If yy(nn) = ( 1) nn xx(nn), find YY ee jj2ππππππ in terms of XX ee jj2ππππππ. Application: Suppose we have a good lowpass FIR filter with transfer function HH(zz) = NN 1 nn=0 h(nn)zz nn, where h(nn) is the length-n impulse response. How can we easily design a length-n FIR high-pass filter? Use frequency shifting. 25
26 Example. For the causal, linear phase lowpass FIR filter with impulse response h(nn) = 1 1, 2,6,2, 1, the DTFT 8 is evaluated as HH ee jj2ππππππ = ee jj2ππππ2tt [6 + 4 cos(2ππππππ) 2 cos(4ππππππ)] 8 and is plotted below for TT = 1. The delay is NN 1 2 = 2 samples. 26
27 Now, to obtain a high-pass filter via frequency translation, we wish to shift by half a period (that is, by 1 Hz) in frequency. Let GG ee jj2ππππππ = HH ee jj 2ππ(ff 1 2TT ) TT = HH ( 1)ee jj2ππππππ, so (replacing ee jj2ππππππ = zz) it follows that GG(zz) = HH( zz), and NN 1 GG(zz) = HH( zz) = h(nn)( zz) nn = nn=0 NN 1 2TT h(nn)( 1) nn zz nn. nn=0 Hence, gg(nn) = ( 1) nn h(nn) is the mapping to convert the lowpass filter to a highpass filter. For h(nn) = 1 1, 2,6,2, 1, this implies that 8 gg(nn) = 1 1, 2,6, 2, 1. 8 Note the dc gain, and the high-frequency gain, of each filter. 27
28 28
29 Comparing the two frequency responses, the highpass filter derived by frequency translation, GG ee jj2ππππππ = HH ee jj 2ππ(ff 1 2TT ) TT In summary, to convert an FIR lowpass filter with impulse response h(nn) to an FIR highpass filter with impulse response gg(nn), using frequency shifting, simply let gg(nn) = ( 1) nn h(nn). 29
30 Example. A few causal lowpass and highpass FIR filters are shown below. The highpass filter impulse response, gg(nn), is formed from the lowpass filter impulse response, h(nn), using gg(nn) = ( 1) nn h(nn). Lowpass filter impulse and frequency responses (N=2,3,4,5). Highpass filter impulse and frequency response (N=2,3,4,5). 30
31 Problem: Suppose we have an IIR lowpass filter given by the difference equation NN yy(nn) = aa kk yy(nn kk) + bb kk xx(nn kk). kk=1 MM kk=0 How can we use frequency translation to modify the lowpass filter into a highpass filter? 31
32 Design Problem. Design a non-causal lowpass filter with sampling rate 1 TT Hz and cutoff frequency ff cc Hz. Solution. The ideal lowpass filter frequency response is 1, 0 ff ff cc HH ee jj2ππππππ = 0, ff cc < ff 1 2TT Use the definition of inverse DTFT to compute h(nn). h iiiiiiiiii (nn) 1 2TT = TT HH ee jj2ππππππ ee jj2ππππππππ dddd 1 2TT ff cc = TT ee jj2ππππππππ ddff = 2ff cc TTsinc(2ff cc TTTT) ff cc A truncated non-causal impulse response is defined as the 2NN + 1 samples of h(nn) over the range nn = NN,,0, NN, h NN (nn) = h iiiiiiiiii(nn), NN nn NN 0, otherwise 32
33 The figures below show h NN (nn) and HH NN ee jj2ππππππ for TT = 1, ff cc = 0.2, and NN = 4, 20, and 50, corresponding to impulse response lengths of 9, 41, and
34 Note: The filters above are non-causal. They satisfy the symmetry condition h NN (nn) = h NN ( nn), so they have zero phase and HH NN ee jj2ππππππ is real-valued. Causal filters with the same frequency response magnitude can be constructed simply by using the DTFT delay property and forming the causal impulse response as h(nn) = h NN (nn NN), nn = 0,,2NN, yielding the length MM = 2NN + 1 FIR filter with frequency response HH ee jj2ππππππ = ee jj2ππππππππ HH NN ee jj2ππππππ. From the construction, the causal filter impulse response satisfies the condition h(nn) = h(mm 1 nn), nn = 0,, MM 1 so the FIR filter has linear phase given by 34
35 θθ(ff) = 2ππππππ(MM 1) 2 = 2ππππππππ. 35
36 This example then provides a design strategy for designing causal, linear phase, FIR filters. 1. Design a zero phase non-causal (generally infinite length) filter impulse response using the inverse DTFT. 2. Truncate the infinite-length impulse response to the sample range nn = NN,,0, NN. Provided h(nn) = h( nn), the truncated filter is still zero phase. 3. Delay the impulse response values by N, yielding a length MM = 2NN 1 causal FIR filter. Note: This general method is known as the window method of linear phase FIR filter design. In general, the truncation window need not be rectangular, but does need to preserve the required symmetry condition, h(nn) = h(mm 1 nn), for nn = 0,, MM 1, with causal impulse response length MM. 36
37 Related Design Problem: Design a length-m linear-phase, causal FIR lowpass filter that is optimum in the MSE sense. That is, if the designed filter has frequency response HH ee jj2ππππππ, and the ideal (causal, linear-phase) lowpass filter has the frequency response, jj2ππππππ(mm 1) HH iiiiiiiiii ee jj2ππππππ = 1ee 2, 0 ff ff cc 0, ff cc < ff 1 2TT Then design HH ee jj2ππππππ to minimize the integral squared error (mean-square error) 1 2TT HH iiiiiiiiii ee jj2ππππππ HH ee jj2ππππππ 2 dddd. 1 2TT Solution. Recall Parseval s relation: 1 2TT xx(nn) 2 = TT XX ee jj2ππππππ 2 dddd. nn= 1 2TT 37
38 Digital Frequency Suppose an analog signal, xx aa (tt), is bandlimited to BB Hz, meaning that XX aa (ff) = 0 for ff BB Hz. The Nyquist sampling rate is 1 TT = 2BB samples/sec. Let xx aa (tt) = cos(2ππff 0 tt) so that the samples are xx(nn) = cos(2ππff 0 nnnn). a) If ff 0 = 0, then xx(nn) =,1, 1, 1,1, (a dc signal) b) If ff 0 = ¼TT, then xx(nn) = cos (2ππ 1 nnnn) = cos 4TT (ππ nn) 2 = (,1,0, 1,0, 1,0, 1,0,1, Hence, there is a sign change every other sample. c) If ff 0 = ½TT, then xx(nn) = cos (2ππ 1 nnnn) = cos ( ππππ) 2TT = (,1, 1, 1, 1,1, There is a sign change every sample. Note the locations of ff = 0, 1, 1 on the unit circle in the 4TT 2TT z-plane. 38
39 Convolution and Steady-State Response Suppose we are given a causal, stable, LTI system with impulse response h(nn). What is the steady-state response to a sinusoidal input, e.g., to input xx(nn) = AA sin(2ππff 0 nn)? Derivation. Assume a normalized sampling rate, so 1 TT = 1. Let the system input be xx(nn) = AAeejj2ππff 0nn for all nn. The system output is given by convolution, Evaluating, yy(nn) = h(nn) xx(nn). yy(nn) = h(kk)xx(nn kk) = h(kk)aaee jj2ππff 0 (nn kk) kk= kk= = AAee jj2ππff 0nn h(kk)ee jj2ππff 0kk kk= = HH ee jj2ππff 0 AAee jj2ππff 0nn Hence, the output is just a scaled version of the input, where the scale factor is the transfer function evaluated at the input frequency. This result also holds for the steady-state 39
40 response if the input begins at some finite time index, and the output is observed after a sufficiently large delay. Example. Let h(nn) = 1 1, 2,1. Determine the steady state 4 response to the following signals (assume TT = 1). i) xx(nn) = 12 uu(nn) ii) xx(nn) = 60 sin ππ nn uu(nn) 2 iii) xx(nn) = 8 cos(ππππ) uu(nn) Solution. HH ee jj2ππππ = ee jj2ππππ + ee jj4ππππ = 1 2 ee jj2ππππ (1 + cos(2ππππ)). xx(nn) = 12 uu(nn) yy ssssssssssssssssssssss (nn) = HH ee jj0 12 = 12. xx(nn) = 60 sin ππ 2 nn uu(nn) yy ssssssssssyyssssssssss(nn) = HH ee jj0.5ππ 60 sin ππ 2 nn + θθ, where θθ = phase HH ee jj0.5ππ. Evaluating, HH ee jj0.5ππ = 1 2 ee jj0.5ππ 1 + cos ππ jj =, so
41 = HH ee jj0.5ππ = 1 ππ and θθ =. Hence, the steady-state 2 2 response to the input in ii) is yy ssssssssssssssssssssss (nn) = sin ππ 2 ππ nn + = 30 sin 2 ππ (nn 1). 2 Note the delay of one sample, matching the linear phase term in the expression for HH ee jj2ππππ. 41
42 Inverse Filters, Minimum Phase, and All-Pass Filters Definition: A LTI system is called minimum phase if all zeros are inside the unit circle. (Hence, a stable minimum phase LTI system is invertible as a causal filter that is also stable.) An interesting filter: HH(zz) = 1 1 αα zz 1 1 ααzz 1. The pole is at zz = αα and we assume that αα < 1 so the filter is stable as a causal filter. The zero is at zz = 1 αα, and since αα < 1, the zero is outside the unit circle. Hence, this filter is non-minimum phase. Assume that αα is real-valued and examine the frequency response. HH ee jj2ππππππ = 1 1 αα ee jj2ππππππ 1 ααee ee jj2ππππππ αα jj2ππππππ = ααee jj2ππππππ 1 1 ααee jj2ππππππ Noting that the two terms in parentheses are related as ααee jj2ππππππ 1 = 1 ααee jj2ππππππ, then HH ee jj2ππππππ = 1, since αα is real-valued. αα. 42
43 Hence, HH(zz) = 1 1 αα zz 1 1 ααzz 1 is an all-pass filter, with magnitude frequency response HH ee jj2ππππππ = 1, i.e., constant for all ff. αα Now, what kind of pole-zero configurations are possible in an all-pass filter? For a difference equation LTI system with real-valued coefficients, all complex-valued poles of zeros must appear in complex-conjugate pairs. So, if zz = αα is a pole, so must be zz = αα. For an all-pass filter, if zz = αα is a pole, then zz = 1 αα must be a zero. So, if αα is complex-valued, then zz = αα and zz = αα are poles, implying that zz = 1 and zz = 1 αα αα must be zeros. Hence, in general, poles and zeros must appear in quads in all-pass filters. If a difference equation all-pass filter is causal and BIBO stable, is the inverse filter BIBO stable as a causal filter? 43
44 Notch Filters A notch filter can be designed to eliminate a specific frequency component, say at frequency ff 0 Hz, by placing a zero of the transfer function at zz = ee jj2ππff 0TT. Assuming this zero location is complex-valued, then a zero need also be placed at zz = ee jj2ππff 0TT so that the filter coefficients are realvalued. Placing poles near the zeros, but inside the unit circle, results in a transfer function of the form HH(zz) = 1 eejj2ππff 0TT zz 1 1 ee jj2ππff 0TT zz 1 (1 αααα jj2ππff 0TT zz 1 )(1 αααα jj2ππff 0TT zz 1 ), where 0 αα < 1 controls how close the pole is to the zero. Simplifying, this becomes HH(zz) = 1 2 cos(2ππff 0TT) zz 1 + zz 2 1 2α cos(2ππff 0 TT) zz 1 + αα 2 zz 2. Example. Design a notch filter to remove 180 Hz, assuming a sampling rate of 1,000 Hz. Solution. Use ff 0 = 180 and TT = 10 3 in the equation above, yielding a transfer function HH(zz) = 1 2 cos(0.36ππ) zz 1 + zz 2 1 2α cos(0.36ππ) zz 1 + αα 2 zz 2. 44
45 The figure below shows the frequency response magnitude for αα = 0, αα = 0.9, and αα = Note how the width of the notch depends on α, and that without poles near the zeros the frequency response is not flat for frequencies significantly distant from ff 0 = 180 Hz. 45
46 Comb Filters These filters can be used to eliminate a periodically spaced set of frequencies, or simply to have a periodic shaped frequency response. We can generate a comb filter from a prototype filter, say HH(zz). A comb filter is formed as HH LL (zz) = HH(zz LL ), for any integer LL > 1. The comb filter frequency response is then HH LL ee jj2ππππππ = HH ee jj2ππππππππ. Since HH ee jj2ππππππ is periodic in frequency with period 1 TT, then as ff varies over the range 0, 1 TT Hz, HH eejj2ππππππππ is varying over LL periods. Hence, HH LL ee jj2ππππππ looks like L periods of HH ee jj2ππππππ, but compressed into the frequency range 0, 1 TT Hz. Example 1. Select LL = 3 and the prototype filter that has impulse response h 0 (nn) = 0.5, 0.5. Then HH(zz) = 1 2 (1 + zz 1 ). The comb filter has transfer function HH LL (zz) = 1 2 (1 + zz LL ), so that with LL = 3, HH 3 (zz) = 1 2 (1 + zz 3 ), implying an impulse 46
47 response h 0 (nn) = 0.5, 0, 0, 0.5. The frequency responses are shown in the figure below. Example 2. A simple IIR filter has transfer function 1 HH(zz) = 1 0.9zz 1. If LL = 4, the comb filter has transfer function HH LL (zz) = zz LL. The frequency responses are shown below. 47
48 Question: Suppose we have a difference equation LTI digital filter. How can we implement the comb filter? Answer: Just replace every delay, zz 1, with zz LL in the realization. 48
ECE 421 Introduction to Signal Processing Project 1 - Solutions
1. (10 credits) Given, xx oooo (tt) = cos (2ππFF cc tt) xx iiii (tt) = cos (2ππππππ) ECE 421 Introduction to Signal Processing Project 1 - Solutions Dror Baron, Spring 2017 The AM modulated output satisfies,
More informationEE 464 Short-Time Fourier Transform Fall and Spectrogram. Many signals of importance have spectral content that
EE 464 Short-Time Fourier Transform Fall 2018 Read Text, Chapter 4.9. and Spectrogram Many signals of importance have spectral content that changes with time. Let xx(nn), nn = 0, 1,, NN 1 1 be a discrete-time
More informationECE 421 Introduction to Signal Processing
ECE 421 Introduction to Signal Processing Dror Baron Assistant Professor Dept. of Electrical and Computer Engr. North Carolina State University, NC, USA Digital Filter Design [Reading material: Chapter
More informationThe University of Texas at Austin Dept. of Electrical and Computer Engineering Final Exam
The University of Texas at Austin Dept. of Electrical and Computer Engineering Final Exam Date: December 18, 2017 Course: EE 313 Evans Name: Last, First The exam is scheduled to last three hours. Open
More informationLecture 17 z-transforms 2
Lecture 17 z-transforms 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/5/3 1 Factoring z-polynomials We can also factor z-transform polynomials to break down a large system into
More informationEE 470 Signals and Systems
EE 470 Signals and Systems 9. Introduction to the Design of Discrete Filters Prof. Yasser Mostafa Kadah Textbook Luis Chapparo, Signals and Systems Using Matlab, 2 nd ed., Academic Press, 2015. Filters
More informationDIGITAL FILTERS. !! Finite Impulse Response (FIR) !! Infinite Impulse Response (IIR) !! Background. !! Matlab functions AGC DSP AGC DSP
DIGITAL FILTERS!! Finite Impulse Response (FIR)!! Infinite Impulse Response (IIR)!! Background!! Matlab functions 1!! Only the magnitude approximation problem!! Four basic types of ideal filters with magnitude
More informationECE 429 / 529 Digital Signal Processing
ECE 429 / 529 Course Policy & Syllabus R. N. Strickland SYLLABUS ECE 429 / 529 Digital Signal Processing SPRING 2009 I. Introduction DSP is concerned with the digital representation of signals and the
More informationSignals. Continuous valued or discrete valued Can the signal take any value or only discrete values?
Signals Continuous time or discrete time Is the signal continuous or sampled in time? Continuous valued or discrete valued Can the signal take any value or only discrete values? Deterministic versus random
More informationELEC-C5230 Digitaalisen signaalinkäsittelyn perusteet
ELEC-C5230 Digitaalisen signaalinkäsittelyn perusteet Lecture 10: Summary Taneli Riihonen 16.05.2016 Lecture 10 in Course Book Sanjit K. Mitra, Digital Signal Processing: A Computer-Based Approach, 4th
More information16.30 Learning Objectives and Practice Problems - - Lectures 16 through 20
16.30 Learning Objectives and Practice Problems - - Lectures 16 through 20 IV. Lectures 16-20 IVA : Sampling, Aliasing, and Reconstruction JVV 9.5, Lecture Notes on Shannon - Understand the mathematical
More informationUNIT IV FIR FILTER DESIGN 1. How phase distortion and delay distortion are introduced? The phase distortion is introduced when the phase characteristics of a filter is nonlinear within the desired frequency
More informationCDS 101/110: Lecture 9.1 Frequency DomainLoop Shaping
CDS /: Lecture 9. Frequency DomainLoop Shaping November 3, 6 Goals: Review Basic Loop Shaping Concepts Work through example(s) Reading: Åström and Murray, Feedback Systems -e, Section.,.-.4,.6 I.e., we
More informationCDS 101/110: Lecture 8.2 PID Control
CDS 11/11: Lecture 8.2 PID Control November 16, 216 Goals: Nyquist Example Introduce and review PID control. Show how to use loop shaping using PID to achieve a performance specification Discuss the use
More informationSignals and Systems Using MATLAB
Signals and Systems Using MATLAB Second Edition Luis F. Chaparro Department of Electrical and Computer Engineering University of Pittsburgh Pittsburgh, PA, USA AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK
More informationDigital Filters IIR (& Their Corresponding Analog Filters) Week Date Lecture Title
http://elec3004.com Digital Filters IIR (& Their Corresponding Analog Filters) 2017 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: Week Date
More informationIIR Filter Design Chapter Intended Learning Outcomes: (i) Ability to design analog Butterworth filters
IIR Filter Design Chapter Intended Learning Outcomes: (i) Ability to design analog Butterworth filters (ii) Ability to design lowpass IIR filters according to predefined specifications based on analog
More informationThe University of Texas at Austin Dept. of Electrical and Computer Engineering Midterm #1
The University of Texas at Austin Dept. of Electrical and Computer Engineering Midterm #1 Date: October 18, 2013 Course: EE 445S Evans Name: Last, First The exam is scheduled to last 50 minutes. Open books
More informationPROBLEM SET 6. Note: This version is preliminary in that it does not yet have instructions for uploading the MATLAB problems.
PROBLEM SET 6 Issued: 2/32/19 Due: 3/1/19 Reading: During the past week we discussed change of discrete-time sampling rate, introducing the techniques of decimation and interpolation, which is covered
More information(i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters
FIR Filter Design Chapter Intended Learning Outcomes: (i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters (ii) Ability to design linear-phase FIR filters according
More informationPart B. Simple Digital Filters. 1. Simple FIR Digital Filters
Simple Digital Filters Chapter 7B Part B Simple FIR Digital Filters LTI Discrete-Time Systems in the Transform-Domain Simple Digital Filters Simple IIR Digital Filters Comb Filters 3. Simple FIR Digital
More information1.5 The voltage V is given as V=RI, where R and I are resistance matrix and I current vector. Evaluate V given that
Sheet (1) 1.1 The voltage across a discharging capacitor is v(t)=10(1 e 0.2t ) Generate a table of voltage, v(t), versus time, t, for t = 0 to 50 seconds with increment of 5 s. 1.2 Use MATLAB to evaluate
More information(i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters
FIR Filter Design Chapter Intended Learning Outcomes: (i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters (ii) Ability to design linear-phase FIR filters according
More information2.1 BASIC CONCEPTS Basic Operations on Signals Time Shifting. Figure 2.2 Time shifting of a signal. Time Reversal.
1 2.1 BASIC CONCEPTS 2.1.1 Basic Operations on Signals Time Shifting. Figure 2.2 Time shifting of a signal. Time Reversal. 2 Time Scaling. Figure 2.4 Time scaling of a signal. 2.1.2 Classification of Signals
More informationMicrophonics. T. Powers
Microphonics T. Powers What is microphonics? Microphonics is the time domain variation in cavity frequency driven by external vibrational sources. A 1.5 GHz structure 0.5 m long will change in frequency
More informationTeam proposals are due tomorrow at 6PM Homework 4 is due next thur. Proposal presentations are next mon in 1311EECS.
Lecture 8 Today: Announcements: References: FIR filter design IIR filter design Filter roundoff and overflow sensitivity Team proposals are due tomorrow at 6PM Homework 4 is due next thur. Proposal presentations
More informationConcordia University. Discrete-Time Signal Processing. Lab Manual (ELEC442) Dr. Wei-Ping Zhu
Concordia University Discrete-Time Signal Processing Lab Manual (ELEC442) Course Instructor: Dr. Wei-Ping Zhu Fall 2012 Lab 1: Linear Constant Coefficient Difference Equations (LCCDE) Objective In this
More informationProblem Point Value Your score Topic 1 28 Discrete-Time Filter Analysis 2 24 Upconversion 3 30 Filter Design 4 18 Potpourri Total 100
The University of Texas at Austin Dept. of Electrical and Computer Engineering Midterm #1 Date: October 17, 2014 Course: EE 445S Evans Name: Last, First The exam is scheduled to last 50 minutes. Open books
More informationy(n)= Aa n u(n)+bu(n) b m sin(2πmt)= b 1 sin(2πt)+b 2 sin(4πt)+b 3 sin(6πt)+ m=1 x(t)= x = 2 ( b b b b
Exam 1 February 3, 006 Each subquestion is worth 10 points. 1. Consider a periodic sawtooth waveform x(t) with period T 0 = 1 sec shown below: (c) x(n)= u(n). In this case, show that the output has the
More informationNH 67, Karur Trichy Highways, Puliyur C.F, Karur District DEPARTMENT OF INFORMATION TECHNOLOGY DIGITAL SIGNAL PROCESSING UNIT 3
NH 67, Karur Trichy Highways, Puliyur C.F, 639 114 Karur District DEPARTMENT OF INFORMATION TECHNOLOGY DIGITAL SIGNAL PROCESSING UNIT 3 IIR FILTER DESIGN Structure of IIR System design of Discrete time
More informationIntroduction to Signals and Systems Lecture #9 - Frequency Response. Guillaume Drion Academic year
Introduction to Signals and Systems Lecture #9 - Frequency Response Guillaume Drion Academic year 2017-2018 1 Transmission of complex exponentials through LTI systems Continuous case: LTI system where
More informationGATE 2014 QUESTION PAPERS EC: Electronics & Communications Engineering - 16th February 2014 (Morning)
ENGISTAN.COM GATE 2014 QUESTION PAPERS EC: Electronics & Communications Engineering - 16th February 2014 (Morning) All GATE 2014 Question Papers & Answers keys are Original Papers & Keys as provided by
More information1 PeZ: Introduction. 1.1 Controls for PeZ using pezdemo. Lab 15b: FIR Filter Design and PeZ: The z, n, and O! Domains
DSP First, 2e Signal Processing First Lab 5b: FIR Filter Design and PeZ: The z, n, and O! Domains The lab report/verification will be done by filling in the last page of this handout which addresses a
More informationEEM478-DSPHARDWARE. WEEK12:FIR & IIR Filter Design
EEM478-DSPHARDWARE WEEK12:FIR & IIR Filter Design PART-I : Filter Design/Realization Step-1 : define filter specs (pass-band, stop-band, optimization criterion, ) Step-2 : derive optimal transfer function
More informationECE438 - Laboratory 7a: Digital Filter Design (Week 1) By Prof. Charles Bouman and Prof. Mireille Boutin Fall 2015
Purdue University: ECE438 - Digital Signal Processing with Applications 1 ECE438 - Laboratory 7a: Digital Filter Design (Week 1) By Prof. Charles Bouman and Prof. Mireille Boutin Fall 2015 1 Introduction
More informationECE503: Digital Filter Design Lecture 9
ECE503: Digital Filter Design Lecture 9 D. Richard Brown III WPI 26-March-2012 WPI D. Richard Brown III 26-March-2012 1 / 33 Lecture 9 Topics Within the broad topic of digital filter design, we are going
More informationHomework Assignment Consider the circuit shown. Assume ideal op-amp behavior. Which statement below is true?
Question 1 (2 points each unless noted otherwise) Homework Assignment 03 1. Consider the circuit shown. Assume ideal op-amp behavior. Which statement below is true? (a) V = VV + = 5 V (op-amp operation)
More informationDigital Filtering: Realization
Digital Filtering: Realization Digital Filtering: Matlab Implementation: 3-tap (2 nd order) IIR filter 1 Transfer Function Differential Equation: z- Transform: Transfer Function: 2 Example: Transfer Function
More informationDiscrete-Time Signal Processing (DTSP) v14
EE 392 Laboratory 5-1 Discrete-Time Signal Processing (DTSP) v14 Safety - Voltages used here are less than 15 V and normally do not present a risk of shock. Objective: To study impulse response and the
More information4. Design of Discrete-Time Filters
4. Design of Discrete-Time Filters 4.1. Introduction (7.0) 4.2. Frame of Design of IIR Filters (7.1) 4.3. Design of IIR Filters by Impulse Invariance (7.1) 4.4. Design of IIR Filters by Bilinear Transformation
More informationUNIT II IIR FILTER DESIGN
UNIT II IIR FILTER DESIGN Structures of IIR Analog filter design Discrete time IIR filter from analog filter IIR filter design by Impulse Invariance, Bilinear transformation Approximation of derivatives
More informationDesign of FIR Filters
Design of FIR Filters Elena Punskaya www-sigproc.eng.cam.ac.uk/~op205 Some material adapted from courses by Prof. Simon Godsill, Dr. Arnaud Doucet, Dr. Malcolm Macleod and Prof. Peter Rayner 1 FIR as a
More informationContinuous-Time Analog Filters
ENGR 4333/5333: Digital Signal Processing Continuous-Time Analog Filters Chapter 2 Dr. Mohamed Bingabr University of Central Oklahoma Outline Frequency Response of an LTIC System Signal Transmission through
More informationAdvanced Digital Signal Processing Part 5: Digital Filters
Advanced Digital Signal Processing Part 5: Digital Filters Gerhard Schmidt Christian-Albrechts-Universität zu Kiel Faculty of Engineering Institute of Electrical and Information Engineering Digital Signal
More informationSECTION 4: TRANSMISSION LINES. ESE 470 Energy Distribution Systems
SECTION 4: TRANSMISSION LINES ESE 470 Energy Distribution Systems 2 Introduction Transmission Lines 3 Transmission and distribution of electrical power occurs over metal cables Overhead AC or DC Underground
More informationDigital Processing of Continuous-Time Signals
Chapter 4 Digital Processing of Continuous-Time Signals 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 4-1-1 Digital Processing of Continuous-Time Signals Digital
More informationAutomatic Control Motion control Advanced control techniques
Automatic Control Motion control Advanced control techniques (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations (I) 2 Besides the classical
More informationGATE 2014: General Instructions during Examination
GATE 2014: General Instructions during Examination 1. Total duration of the GATE examination is 180 minutes. 2. The clock will be set at the server. The countdown timer at the top right corner of screen
More informationFinal Exam Practice Questions for Music 421, with Solutions
Final Exam Practice Questions for Music 4, with Solutions Elementary Fourier Relationships. For the window w = [/,,/ ], what is (a) the dc magnitude of the window transform? + (b) the magnitude at half
More informationDigital Signal Processing
Digital Signal Processing Fourth Edition John G. Proakis Department of Electrical and Computer Engineering Northeastern University Boston, Massachusetts Dimitris G. Manolakis MIT Lincoln Laboratory Lexington,
More informationAdaptive Filters Application of Linear Prediction
Adaptive Filters Application of Linear Prediction Gerhard Schmidt Christian-Albrechts-Universität zu Kiel Faculty of Engineering Electrical Engineering and Information Technology Digital Signal Processing
More informationThe University of Texas at Austin Dept. of Electrical and Computer Engineering Midterm #2
The University of Texas at Austin Dept. of Electrical and Computer Engineering Midterm #2 Date: November 18, 2010 Course: EE 313 Evans Name: Last, First The exam is scheduled to last 75 minutes. Open books
More informationPower Spectral Density and Hilbert Transform
ARL-TR-7920 DEC 2016 US Army Research Laboratory Power Spectral Density and Hilbert Transform by Patrick Jungwirth NOTICES Disclaimers The findings in this report are not to be construed as an official
More informationFinal Exam Solutions June 7, 2004
Name: Final Exam Solutions June 7, 24 ECE 223: Signals & Systems II Dr. McNames Write your name above. Keep your exam flat during the entire exam period. If you have to leave the exam temporarily, close
More informationDigital Processing of
Chapter 4 Digital Processing of Continuous-Time Signals 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 4-1-1 Digital Processing of Continuous-Time Signals Digital
More informationExperiment 4- Finite Impulse Response Filters
Experiment 4- Finite Impulse Response Filters 18 February 2009 Abstract In this experiment we design different Finite Impulse Response filters and study their characteristics. 1 Introduction The transfer
More informationSystem analysis and signal processing
System analysis and signal processing with emphasis on the use of MATLAB PHILIP DENBIGH University of Sussex ADDISON-WESLEY Harlow, England Reading, Massachusetts Menlow Park, California New York Don Mills,
More informationDSP Laboratory (EELE 4110) Lab#10 Finite Impulse Response (FIR) Filters
Islamic University of Gaza OBJECTIVES: Faculty of Engineering Electrical Engineering Department Spring-2011 DSP Laboratory (EELE 4110) Lab#10 Finite Impulse Response (FIR) Filters To demonstrate the concept
More informationEE 422G - Signals and Systems Laboratory
EE 422G - Signals and Systems Laboratory Lab 3 FIR Filters Written by Kevin D. Donohue Department of Electrical and Computer Engineering University of Kentucky Lexington, KY 40506 September 19, 2015 Objectives:
More informationThe performance of AM and FM receivers. Editor: Xuanfeng Li Teacher: Prof. Xiliang Luo
The performance of AM and FM receivers Editor: Xuanfeng Li Teacher: Prof. Xiliang Luo The performance of AM receivers using Envelop Detection In a full AM signal, both sidebands and the carrier wave are
More informationFUNDAMENTALS OF SIGNALS AND SYSTEMS
FUNDAMENTALS OF SIGNALS AND SYSTEMS LIMITED WARRANTY AND DISCLAIMER OF LIABILITY THE CD-ROM THAT ACCOMPANIES THE BOOK MAY BE USED ON A SINGLE PC ONLY. THE LICENSE DOES NOT PERMIT THE USE ON A NETWORK (OF
More informationLab #2: Electrical Measurements II AC Circuits and Capacitors, Inductors, Oscillators and Filters
Lab #2: Electrical Measurements II AC Circuits and Capacitors, Inductors, Oscillators and Filters Goal: In circuits with a time-varying voltage, the relationship between current and voltage is more complicated
More informationELECTRONOTES APPLICATION NOTE NO Hanshaw Road Ithaca, NY August 3, 2017
ELECTRONOTES APPLICATION NOTE NO. 432 1016 Hanshaw Road Ithaca, NY 14850 August 3, 2017 SIMPLIFIED DIGITAL NOTCH FILTER DESIGN Recently [1] we have been involved with an issue of a so-called Worldwide
More informationBibliography. Practical Signal Processing and Its Applications Downloaded from
Bibliography Practical Signal Processing and Its Applications Downloaded from www.worldscientific.com Abramowitz, Milton, and Irene A. Stegun. Handbook of mathematical functions: with formulas, graphs,
More informationGATE 2014: General Instructions during Examination
GATE 2014: General Instructions during Examination 1. Total duration of the GATE examination is 180 minutes. 2. The clock will be set at the server. The countdown timer at the top right corner of screen
More informationDigital Video and Audio Processing. Winter term 2002/ 2003 Computer-based exercises
Digital Video and Audio Processing Winter term 2002/ 2003 Computer-based exercises Rudolf Mester Institut für Angewandte Physik Johann Wolfgang Goethe-Universität Frankfurt am Main 6th November 2002 Chapter
More information(i) Understanding the basic concepts of signal modeling, correlation, maximum likelihood estimation, least squares and iterative numerical methods
Tools and Applications Chapter Intended Learning Outcomes: (i) Understanding the basic concepts of signal modeling, correlation, maximum likelihood estimation, least squares and iterative numerical methods
More informationF I R Filter (Finite Impulse Response)
F I R Filter (Finite Impulse Response) Ir. Dadang Gunawan, Ph.D Electrical Engineering University of Indonesia The Outline 7.1 State-of-the-art 7.2 Type of Linear Phase Filter 7.3 Summary of 4 Types FIR
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #25 Wednesday, November 5, 23 Aliasing in the impulse invariance method: The impulse invariance method is only suitable for filters with a bandlimited
More informationijdsp Workshop: Exercise 2012 DSP Exercise Objectives
Objectives DSP Exercise The objective of this exercise is to provide hands-on experiences on ijdsp. It consists of three parts covering frequency response of LTI systems, pole/zero locations with the frequency
More informationDiode Circuits Recent GATE Problems
Diode Circuits Recent GATE Problems 1. The diodes and capacitors in the circuit shown are ideal. The voltage v(t) across the diode DD 1 is CC 1 DD 2 cos(ωωωω) AC DD 1 CC 1 (a) cos(ωωωω) 1 (b) sin(ωωωω)
More informationDigital Filters FIR and IIR Systems
Digital Filters FIR and IIR Systems ELEC 3004: Systems: Signals & Controls Dr. Surya Singh (Some material adapted from courses by Russ Tedrake and Elena Punskaya) Lecture 16 elec3004@itee.uq.edu.au http://robotics.itee.uq.edu.au/~elec3004/
More informationBasic Signals and Systems
Chapter 2 Basic Signals and Systems A large part of this chapter is taken from: C.S. Burrus, J.H. McClellan, A.V. Oppenheim, T.W. Parks, R.W. Schafer, and H. W. Schüssler: Computer-based exercises for
More informationHomework Assignment 09
Question 1 (2 points each unless noted otherwise) Homework Assignment 09 1. For SPICE, Explain very briefly the difference between the multiplier M and Meg, as in a resistor has value 2M versus a resistor
More informationECE 5650/4650 Exam II November 20, 2018 Name:
ECE 5650/4650 Exam II November 0, 08 Name: Take-Home Exam Honor Code This being a take-home exam a strict honor code is assumed. Each person is to do his/her own work. Bring any questions you have about
More informationCopyright S. K. Mitra
1 In many applications, a discrete-time signal x[n] is split into a number of subband signals by means of an analysis filter bank The subband signals are then processed Finally, the processed subband signals
More informationBrief Introduction to Signals & Systems. Phani Chavali
Brief Introduction to Signals & Systems Phani Chavali Outline Signals & Systems Continuous and discrete time signals Properties of Systems Input- Output relation : Convolution Frequency domain representation
More informationTABLE OF CONTENTS TOPIC NUMBER NAME OF THE TOPIC 1. OVERVIEW OF SIGNALS & SYSTEMS 2. ANALYSIS OF LTI SYSTEMS- Z TRANSFORM 3. ANALYSIS OF FT, DFT AND FFT SIGNALS 4. DIGITAL FILTERS CONCEPTS & DESIGN 5.
More informationPATCH ANTENNA DESIGN
FACULTY OF ELECTRONICS, TELECOMMUNICATIONS, AND INFORMATION TECHNOLOGY Department of Telecommunications and Information Technology Laboratory of Antennas and Propagation Prof. dr. ing. Ion BOGDAN PATCH
More informationANGLE MODULATION. U1. PHASE AND FREQUENCY MODULATION For angle modulation, the modulated carrier is represented by
[4.1] ANGLE MODULATION U1. PHASE AND FREQUENCY MODULATION For angle modulation, the modulated carrier is represented by xx cc (tt) = AA cccccc[ωω cc tt + φφ(tt)] (1.1) Where A ω c are constants the phase
More informationLecture 3 Review of Signals and Systems: Part 2. EE4900/EE6720 Digital Communications
EE4900/EE6720: Digital Communications 1 Lecture 3 Review of Signals and Systems: Part 2 Block Diagrams of Communication System Digital Communication System 2 Informatio n (sound, video, text, data, ) Transducer
More informationDigital Signal Processing
Digital Signal Processing System Analysis and Design Paulo S. R. Diniz Eduardo A. B. da Silva and Sergio L. Netto Federal University of Rio de Janeiro CAMBRIDGE UNIVERSITY PRESS Preface page xv Introduction
More informationProblem Point Value Your score Topic 1 28 Discrete-Time Filter Analysis 2 24 Improving Signal Quality 3 24 Filter Bank Design 4 24 Potpourri Total 100
The University of Texas at Austin Dept. of Electrical and Computer Engineering Midterm #1 Date: March 7, 2014 Course: EE 445S Evans Name: Last, First The exam is scheduled to last 50 minutes. Open books
More informationUnderstanding Digital Signal Processing
Understanding Digital Signal Processing Richard G. Lyons PRENTICE HALL PTR PRENTICE HALL Professional Technical Reference Upper Saddle River, New Jersey 07458 www.photr,com Contents Preface xi 1 DISCRETE
More information8: IIR Filter Transformations
DSP and Digital (5-677) IIR : 8 / Classical continuous-time filters optimize tradeoff: passband ripple v stopband ripple v transition width There are explicit formulae for pole/zero positions. Butterworth:
More informationChapter 2. Fourier Series & Fourier Transform. Updated:2/11/15
Chapter 2 Fourier Series & Fourier Transform Updated:2/11/15 Outline Systems and frequency domain representation Fourier Series and different representation of FS Fourier Transform and Spectra Power Spectral
More informationStudying Noise Contributions in Nonlinear Vector Network Analyzer (NVNA) Measurements
FACULTY OF ENGINEERING AND SUSTAINABLE DEVELOPMENT Studying Noise Contributions in Nonlinear Vector Network Analyzer (NVNA) Measurements Feng Tianyang September 2012 Master s Thesis in Telecommunications
More informationOutline. Discrete time signals. Impulse sampling z-transform Frequency response Stability INF4420. Jørgen Andreas Michaelsen Spring / 37 2 / 37
INF4420 Discrete time signals Jørgen Andreas Michaelsen Spring 2013 1 / 37 Outline Impulse sampling z-transform Frequency response Stability Spring 2013 Discrete time signals 2 2 / 37 Introduction More
More informationEECS 452 Practice Midterm Exam Solutions Fall 2014
EECS 452 Practice Midterm Exam Solutions Fall 2014 Name: unique name: Sign the honor code: I have neither given nor received aid on this exam nor observed anyone else doing so. Scores: # Points Section
More informationMultirate Digital Signal Processing
Multirate Digital Signal Processing Basic Sampling Rate Alteration Devices Up-sampler - Used to increase the sampling rate by an integer factor Down-sampler - Used to increase the sampling rate by an integer
More informationModule 9 AUDIO CODING. Version 2 ECE IIT, Kharagpur
Module 9 AUDIO CODING Lesson 30 Polyphase filter implementation Instructional Objectives At the end of this lesson, the students should be able to : 1. Show how a bank of bandpass filters can be realized
More informationB.Tech III Year II Semester (R13) Regular & Supplementary Examinations May/June 2017 DIGITAL SIGNAL PROCESSING (Common to ECE and EIE)
Code: 13A04602 R13 B.Tech III Year II Semester (R13) Regular & Supplementary Examinations May/June 2017 (Common to ECE and EIE) PART A (Compulsory Question) 1 Answer the following: (10 X 02 = 20 Marks)
More informationLab S-9: Interference Removal from Electro-Cardiogram (ECG) Signals
DSP First, 2e Signal Processing First Lab S-9: Interference Removal from Electro-Cardiogram (ECG) Signals Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab.
More informationEE3079 Experiment: Chaos in nonlinear systems
EE3079 Experiment: Chaos in nonlinear systems Background: November 2, 2016 Revision The theory of nonlinear dynamical systems and Chaos is an intriguing area of mathematics that has received considerable
More informationSummer Homework. Trace each number. Count to 10. Complete the picture. Tell a story about your picture..
Week 1 {June 4} Read every day! Parent Initial Week 2 {June 11} Read every day! Parent Initial June Summer Homework Monday Tuesday Wednesday Thursday Trace the letters. Color each of the pictures. Match
More informationInstruction Manual for Concept Simulators. Signals and Systems. M. J. Roberts
Instruction Manual for Concept Simulators that accompany the book Signals and Systems by M. J. Roberts March 2004 - All Rights Reserved Table of Contents I. Loading and Running the Simulators II. Continuous-Time
More informationE Final Exam Solutions page 1/ gain / db Imaginary Part
E48 Digital Signal Processing Exam date: Tuesday 242 Final Exam Solutions Dan Ellis . The only twist here is to notice that the elliptical filter is actually high-pass, since it has
More information8.0 Overview 8.1 The frequency response for a stable system
8. Frequency Response Analysis 8.0 Overview 8.1 The frequency response for a stable system 8.1.1 Simple system elements 8.1.2 First-order system 8.1.3 Second-order systems 8.1.4 Systems of arbitrary order
More informationExperiment 2 Effects of Filtering
Experiment 2 Effects of Filtering INTRODUCTION This experiment demonstrates the relationship between the time and frequency domains. A basic rule of thumb is that the wider the bandwidth allowed for the
More informationIf You Give A Mouse A Letter He ll Want The Whole Alphabet By
If You Give A Mouse A Letter He ll Want The Whole Alphabet By a TeachWithMe.com A If you give a mouse an A he will want a Bb. AAAA aaaa B If you give a mouse an B he will want a Cc. BBBB bbbb C If you
More information