讲解：ECE 425/525、Digital Signal Pr

ECE 425/525 Digital Signal ProcessingSpring 2019Dept. of Electrical and Computer Engineering Miami University1Homework 11 AssignmentWeights: 60ptsDue Date: May 10, 2019.Turn in Method: Turn in your report with your Matlab? program and other supporting materialslisted at the end of this handout via Canvas by the midnight of May 10, 2019.Project Description:In this project, you will use two bandstop IIR filters to filter out noise.Conduct the following task.1. Download the test wave file, ‘test3.wav’, from the Canvas.2. Use Matlab audioread command (or other Matlab command) to find out the samplingfrequency of the wave file and save the voice to x[k].3. Play x[k] with Matlab?.4. Get the magnitude and phase spectra of x[k] with fft command and plot them. Since thelength of x[k] can be too long, you might want to specify the size of FFT output. Refer toMatlab help menu on how to specify the length of FFT result. (Note: x-axis unit needs tobe Hz and y-axis unit needs to be dB (magnitude spectrum) or radian (phase spectrum))a. Since fft(x) will generate a complex valued array, you need to use Matlabcommands ‘abs’ and ‘angle’ to get magnitude and phase spectrum of yoursignal.b. You only need to plot the first-half of you FFT result. Why?5. The signal (x[k]) is corrupted by two narrowband interferences. Based on the magnitudespectrum you obtain in step 4, design two bandstop digital IIR filters so you can use themsequen代写ECE 425/525作业、代做Digital Signal Processing作业、Matlab编程作业调试、Mtially to remove these interferences.6. Plot the frequency response spectra of the bandstop filters with freqz command (x-axisunit: Hz)7. Apply your two bandstop filters in sequence to filter the x[k] and save your result into y[k].8. Play y[k].9. Get the magnitude and phase spectra of y[k] with fft command and plot them. (Note: xaxisunit needs to be Hz and y-axis unit needs to be dB (magnitude spectrum) or radian(phase spectrum))10. Save y[k] into a wav format file.Submission RequestYou need to turn in the followings:(1) Matlab? codes(2) The filtered voice in wav format(3) The frequency response spectra of two bandstop filtersECE 425/525 Digital Signal ProcessingSpring 2019Dept. of Electrical and Computer Engineering Miami University2(4) The magnitude and phase spectra of x[k] and y[k].Matlab? Commands You Might Needaudioplayer, play, wavread, sound, xcorr, zeros, ones, fft, ifft, abs, angle, unwrap, real,imag, freqz, freqzplot, filter, subplot, buttord, cheb1ord, cheb2ord, ellipord, butter, cheby1,cheby2, ellip, lp2bp, bilinear.Note: These commands are provided for you convenience, you should not assume that they arethe only commands you need for this assignment.Example Matlab Program for filter design (Read Matlab help for ellipord and ellip)Fs=8000;Wp1=[300 3000];Ws1=[500 2000];[N, Wp] = ellipord(Wp1/(Fs/2), Ws1/(Fs/2), 0.5, 40)[B1,A1] = ellip(N,0.5, 40,Wp,stop);freqz(B1,A1,1024,Fs);转自：http://www.7daixie.com/2019051221714067.html