关键词 > EENG31400/EENGM1400

EENG31400 / EENGM1400 DIGITAL FILTERS AND SPECTRAL ANALYSIS 2019

发布时间:2023-01-09

Hello, dear friend, you can consult us at any time if you have any questions, add WeChat: daixieit

January 2019

EENG31400 / EENGM1400

DIGITAL FILTERS AND SPECTRAL ANALYSIS

Q1       Figure  1  shows  the  magnitude  of  the  32-point  DFT  X[k],  k=0,....31  of  discrete  time  signal

x[n]= ej!!n  for the case of window w[n]:

0 ≤ n 7

theTwiSe

X [k]

Figure 1

(a)     Sketch the DTFT X(Ω) ofx[n]

 (2 marks)

(b)     Comment on the shape of the spectrum X[k] provided by the DFT relative to that of X(Ω).

Why are there so many frequency components in the DFT spectrum? Explain.

The DTFT spectrum of signal x[n] = ej!!n  consists of a single impulse at frequency  Ω ! . Due to the windowing process that takes place when we calculate the DFT, the spectrum sampled by the DFT is the DTFT spectrum of x[n] x w[n], i.e. the spectrum of the rectangular window convolved  with  the  spectrum  of x[n]  (impulse  at  frequency    Ω !).  As  a  result  multiple frequency components appear. (8 marks)

(c)     Which normalised radial frequency Ω do the following indices k correspond to:

•    k = 16

    k = 6

K = 16: Ω 1 =16x2π/32= π

K = 6: Ω2 = 6x2π/32 (4 marks)

(d)     Which of the two actions described below would bring the shape of the DFT spectrum X[k] closer to that ofX(Ω)

a)   Increasing the length of the DFT to 256

b)   Increasing the length of window w [n] to 32

Justify your answer

The answer is (b). Increasing the length of the DFT (equivalent to padding the windowed signal with zeroes) will not change the shape of the spectrum. Increasing the length of the DFT means sampling the DTFT spectrum of the windowed signal more densely. The shape of the DTFT spectrum (described above) remains the same, as it depends on the length and type of window used. The longer the window the smaller the main lobe width in the spectrum. (6 marks)

Q2      (a)     Determine whether each of the following statements is true or false (where u is the unit step

function). Justify your answers.

i.      The signal x(t) = u(t+T0) - u(t-T0) will not suffer from aliasing if it is sampled with a sampling period T<2T0

ii.      The signal x(t) with Fourier transform X(  ) = u(  +  0)- u(  -  0) will not suffer from aliasing if it is sampled with a sampling period T <  /  !

iii.      The signal x(t) with Fourier transform X(  ) = u(  )- u(  -  0) will not suffer from aliasing if it is sampled with a sampling period T < 2/  !

 (9 marks)

(b)    Consider  an  ideal  discrete-time band-stop  filter with  impulse response  h[n]  for which  the

frequency response in the interval  ≤ Ω ≤  is :

1, Ω    and Ω   3

0,         elsewheTe

Determine the frequency response of the filter whose impulse response is h[2n]

 (4 marks)

(c)    Consider the system of Figure 2 for changing the sampling rate of signal x[n] by a non-integer factor. Determine the output Xd [n] ifx[n] = cos(3πn / 4) , L=6 and M=7.

 

d[n] =  cos(78(冗)n  )  

(d)

Figure 2

(3 marks)

(b)

Figure 3

Suppose that a continuous time speech signal sc(t) with the Fourier transform spectrum Sc(ω) shown in Figure 3(a) is processed by the system shown in Figure 3(b) to produce the discrete time signal sr [n]. H(Ω) in Figure 3(b) is an ideal discrete-time low pass filter with a cut-off frequency of Ωc  and a gain of L in the passband. The signal sr [n] will be used as input to a speech encoder that operates correctly only on discrete-time speech signals sampled at an 8- kHz rate. Choose values of L, M and Ωc  that produce the correct input signal sr [n] for the speech encoder.

The sampling frequency is equal to 44.1 KHz. We want the output to have an 8KHz rate hence we need  =  =  . We choose L=80, M=441 and Qc=冗/441 (4 marks)

Q3        (a)     Following the windowing method of FIR filter design, use a Hamming window to design a

3-tap FIR lowpass filter with a cut-off frequency of 800 Hz and a sampling rate of 8 kHz.

1.     Q = 2.f / f s   = 2800/ 8000 = / 5

(1     Q 

6.    h,[n] = 0.187   0.016   0.187

7.    Hw [z] = 0. 187(z 1  + z -1 ) + 0.016 z 0 (12 marks)

(b)

 

iv) 

Figure 4

(8 marks)