ECE 5650/4650 MATLAB Project 1
|
|
- Winfred Poole
- 5 years ago
- Views:
Transcription
1 This project is to be treated as a take-home exam, meaning each student is to due his/her own work. The project due date is 4:30 PM Tuesday, October 18, To work the project you will need access to MATLAB and at minimum the signal processing toolbox. Introduction This first MATLAB DSP project will get you acquainted/reacquainted with MATLAB and then move into the exploration of Discrete-time signal generation Sequence convolution and the convolution sum Linear constant coefficient difference equations (LCCDEs) The discrete-time Fourier transform (DTFT) It is the student s responsibility to learn the very basics of MATLAB. There are many resources on the internet. One such example is course notes available in PDF form, from A Quick MATLAB Review Recall that MATLAB stands for matrix laboratory, so it goes without saying that MATLAB is very efficient at doing matrix oriented numerical calculations. MATLAB supports many data types, but the two data types we will work with the most, are scalar numbers (real and complex), and matrices (real and complex). For our signal processing needs the matrices we create will most often be vectors, that is matrices of dimension m 1 or 1 n. The default data type for real and complex numbers is double precision. Complex numbers are handled transparently in MATLAB, so there is no need to worry about these details, for the most part. In fact, when MATLAB starts-up, it has both i and j defined as 1, i.e., at the command window, as shown in Figure 1, enter >> j ans = i Note that if you use i and j for something else, these definitions will be overwritten. The matrix specialization of MATLAB means that the most efficient programming requires vectorizing. This is an aspect of using MATLAB that most students master over time. Getting your code to work in some fashion is first and foremost. Writing tight code is something that requires practice. The Colon Operator Suppose we wish to create a vector of time samples for We enter at the command window >> t = 0:0.01:5; >> size(t) t 05 s, with sample spacing of 0.01s. Introduction 1
2 Figure 1: The MATLAB integrated development environment configured with the function/script editor in the upper right, the command window in the lower right, the current directory, workspace, and figure windows all tabbed in the upper left, and finally the command history in the lower left. ans = >> cls % clear the command window screen >> clf % clear the active figure window >> figure(1) % open a new figure window or address #1 if already open The variable t is now a row vector containing one row and 501 columns. If all we wanted was integer steps, we could simply enter >> n = 0:10 n = Note the semicolon at the end of a statement, suppresses the automatic listing of the return variable. We can create a function of the vector t, say xt = cos 2 2 t, by writing >> x = cos(2*pi*2*t); A plot of x versus t can be made using the plot() command. >> plot(t,x) In Figure 1 we see the results. For the new user, the MATLAB environment needs to be customized by the user to take on the appearance of Figure 1, or whatever IDE set-up you prefer. We started out showing how the colon operator can be used to create a vector. The colon operator can also be used to access portions of an array or matrix, e.g. Introduction 2
3 >> V = rand(6,6) V = >> V(1:2,4:5) ans = >> V(:,3:5) ans = Note that rand() produces uniformly distributed random numbers on the interval [0, 1]. It can also be used to store values into an array or matrix, e.g., >> W = zeros(6,3) W = >> W(2:4,2:3) = rand(3,2) W = Functions One of the key building blocks in constructing a MATLAB simulation, is the ability to write custom user functions. To demonstrate this, we will create a function that generates a sinusoid with a delayed turn-on time, i.e., xt = Acos 2 f 0 t + ut t 0 (1) The signal model above is a continuous-time sinusoid, but the MATLAB produced signal will be discrete-time in nature, because we will be letting t nt = n f s, where f s is the sampling rate. The end result will be a signal vector xn with corresponding integer time axis n and a corresponding time axis t. Introduction 3
4 We require the MATLAB function to take as inputs the frequency f 0, the phase, the sampling rate f s, the delay time t 0, and the time duration for the waveform, Tsim. The amplitude is taken to be one. The function listing is given below: function [x,t,n] = cos_gen(f0,phi,t0,fs,tsim) % [x,t,n] = cos_gen(f0,phi,t0,fs,tsim) % Generate a cos signal/sequence for Tsim s %================ Inputs ========================================= % f0 = frequency in Hz % phi = phase in radians % to = turn-on time delay in s % fs = sampling rate in samples per second % Ts = simulation duration in s %================ Outputs ======================================== % x = signal vector (a row vector) % t = time axis in s corresponding to x % n = integer axis corresponding to t %//////////////////////////////////////////////////////////////// % Mark Wickert September 2008 % Create the time axis >> t = 0:1/fs:Tsim; % Create the corresponding integer time axis >> n = 0:length(t)-1; % Cosine wave starting at t = 0 >> x = cos(2*pi*f0*t+phi); % Zero samples for t < t0 using the find function % The vector idx holds the indices where t < t0 >> idx = find(t < t0); % Replace the elements of x having the indices of idx % with a vector of zeros [x,t,n] = cos_gen(10,0,.2,200,0.5); 1 [x,t,n] = cos_gen(10,pi/2,.2,200,0.5); 1 Amplitude Amplitude Time (s) [x,t,n] = cos_gen(10,pi/2,.1,200,0.5); Time (s) [x,t,n] = cos_gen(10,0,.05,200,0.25); 1 Amplitude Amplitude Time (s) Time (s) Figure 2: Plots obtained from the MATLAB function cos_gen(). Introduction 4
5 >> x(idx) = zeros(size(idx)); To demonstrate the function in action consider T sim = 5 s for various value of the other parameters, as shown in Figure 2. The command line used to create each plot is shown as the title of each plot. Each plot is created as an overlay of a line plot and plot of just the sample points, as follows: >> [x,t,n] = cos_gen(10,0,.2,200,0.5); >> plot(t,x) >> hold Current plot held >> plot(t,x,.r ) >> grid >> xlabel( Time (s) ) >> ylabel( Amplitude ) >> title( [x,t,n] = cos_gen(10,0,.2,200,0.5); ) Note that the array of four plots is crated using the subplot() command, e.g., >> subplot(2,2,1) % a 2 x 2 array, the first element, L to R/T to B, is >> plot() % etc We also have created embedded help for the function cos_gen() by including comment lines immediately following the function declaration, e.g., function [x,t,n] = cos_gen(f0,phi,t0,fs,tsim) % [x,t,n] = cos_gen(f0,phi,t0,fs,tsim) % Generate a cos signal/sequence for Tsim s %================ Inputs ========================================= % f0 = frequency in Hz % phi = phase in radians% to = turn-on time delay in s % fs = sampling rate in samples per second % Ts = simulation duration in s %================ Outputs ======================================== % x = signal vector (a row vector) % t = time axis in s corresponding to x % n = integer axis corresponding to t %//////////////////////////////////////////////////////////////// % Mark Wickert September 2008 The result when we type help function_name at the command prompt is: >> help cos_gen [x,t,n] = cos_gen(f0,phi,t0,fs,tsim)generate a cos signal/sequence for Tsim s ================ Inputs ========================================= f0 = frequency in Hz phi = phase in radians to = turn-on time delay in s fs = sampling rate in samples per second Ts = simulation duration in s ================ Outputs ======================================== x = signal vector (a row vector)t = time axis in s corresponding to x n = integer axis corresponding to t Introduction 5
6 //////////////////////////////////////////////////////////////// Mark Wickert September 2008 Problems Signal Generation In this problem you will write a MATLAB function for generating a pulse train signal. You will then explore some properties of this signal. The signal will be a purely discrete-time signal of the form where here pn N p 1 xn = pn km D 0 n N 1 k = 0 is a rectangular pulse shape of the form (2) pn = 1, 0 n W 1 0, otherwise (3) Write your code so that in the future a different pulse shape may be utilized. 1. Problem statement a) Write the function discrete_pulse_train() to have the following function prototype: function [x,n] = discrete_pulse_train(w,m,np,d,n) % [x,n] = discrete_pulse_train(w,m,np,d,n) % Generate a discrete-time pulse train %================ Inputs ========================================= % W = pulse width in samples % M = period in samples, M > W % Np = Number of periods to create (if 0 full out to N samples) % D = turn-on delay in samples of the first period % N = total number of samples to create %================ Outputs ======================================== % x = signal vector (a row vector) % n = integer time axis corresponding to x %//////////////////////////////////////////////////////////////// An example plot from this function obtained using the command line >> [x,n] = discrete_pulse_train(5,18,4,13,200); is shown in Figure 3. Problems 6
7 [x,n] = discrete_pulse_train(5,18,4,13,200); Amplitude Samples (n) Figure 3: Plots obtained from the MATLAB function discretepulse_train(). b) Test the function and produce output stem() function plots for the function parameter entries given in Table 1. Table 1: Parameter values for test plots to be created with the function discrete_ pulse_train.m. W M Np D N (limited by N ) c) Working from the command window, create a co-sinusoidal pulse train signal of the form yn = N p 1 cos 0 n pn km n 0 0 n N p 1 k = 0 (4) Problems 7
8 Convolution This signal could be used as the basis for a pulsed RADAR simulation waveform. As a specific test case, generate and plot the output vector y for >> [x,n] = discrete_pulse_train(50,500,0,0,10000); and 0 = 10. Here with so many points in the signal, it is OK to use plot() instead of stem(). If say the sampling rate is f s = 1 T = 10 MHz, what is the corresponding pulse width in seconds and pulse repetition frequency in Hz? Here we consider the mechanics of the convolution sum and related properties. We will utilize the pulse train signal generator from Problem 1 and MATLAB s conv() function to validate some familiar concepts from continuous-time convolution. 2. Consider sequences x 1 and x 2 each defined on n 019 as follows: >> x1 = discrete_pulse_train(5,20,1,5,20); >> x2 = discrete_pulse_train(8,20,1,0,20); a) Obtain y 1 = x 1 *x 1 and plot it using the stem() function with a properly scaled sequence axis n. Comment your plot relative to theoretical expectations with regard to nonzero duration (support interval) of the output, the shape, and the peak value. b) Repeat part (a), except now explore y 11 = x 1 *x 1 *x 1. What can you say in general about the support interval of y 11 and the pulse duration? c) Repeat part (a), except now explore y 12 = x 1 *x We will now consider multiple pulses in signal x 1 : >> x3 = discrete_pulse_train(5,20,2,5,60); >> x4 = discrete_pulse_train(5,20,0,5,50); a) Obtain y 3 = x 3 *x 3 and plot it using the stem() function with a properly scaled sequence axis n. Comment your plot relative to theoretical expectations with regard to nonzero duration (support interval) of the output, the shape, and the peak value. b) Repeat part (a) for y 4 = x 4 *x 4. Is the result consistent with what your observed in part (a)? What pattern is developing? Problems 8
9 Difference Equations, Frequency Response, and Filtering To begin with we will consider a simple exponential sequence and how it can be generated via a recursion or difference equation. Then we will consider filtering the pulse train signal of Problem 1 using both FIR and IIR filters. 4. The real exponential pulse signal has both finite and infinite duration forms: y 1 = a n un un L, (5) y 2 = a n un where in both cases we assume that 1 a 1. a) At the command window generate and plot using the stem()function, y 1 for L = 40, a = 0.8, and 0 n 50. Hint: You should be able to use the pulse train function to help generate this signal. Now feed y 1 through an accumulator system which has difference equation z 1 = z 1 n 1 + y 1. (6) You can use the MATLAB function filter() to make this easy. Note the final or steady-state value of z 1. Compare it with the theoretical value using the geometric series sum formula S L 1 a n = = n = 0 1 a L (7) 1 a b) Repeat part (a) using y 2, but this time generate the signal by driving into the system having difference equation yn = ay n 1 + xn. (8) In this case let 0 n 100. Again compare the approximate steady-state value with the theoretical value obtained via the infinite geometric series. 5. In this problem we consider the steady-state response of an IIR notch filter. To begin with we know that when a sinusoidal signal of the form xn = Acos 0 n + n (9) is passed through an LTI system having frequency response the form He j, the system output is of yn AH e j 0 0 n + He j 0 = cos + n (10) In this problem the input will be of the form xn = Acos 0 n + un, so the output will consist of both the transient and steady-state responses. For n large the output is in steady-state, and we should be able to determine the magnitude and phase response at 0 from the waveform. The filter of interest is Problems 9
10 He j = e j e j2 1 2cos rcos 0 e j + r 2 e j2 (11) where controls the center frequency of the notch and r controls the bandwidth. 0 a) Using freqz() plot the magnitude and phase response of the notch for 0 = 2 and r = 0.9. Plot the magnitude response as both a linear plot, He j, and in db, i.e., 20log 10 He j. b) Using filter(), input cos 3 n for 0 n 50. Determine from the output sequence plot approximate values for the filter gain and phase shift at = 3. Note that the filter still has center frequency 0 = 2. Also, plot the transient response as a separate plot, from your knowledge of the theoretical steady-state response. c) Assume that the filter operates within an A/D-He j -D/A system. An interfering signal is present at 120 Hz and the system sampling rate is set at 2000 Hz. Determine 0 so that the filter removes the 120 Hz signal. Simulate this system by plotting the filter output in response to the interference signal input. Determine in ms, how long it takes for the filter output to settle to yn 0.01, assuming the input amplitude is one. Sub-Optimal Filtering of Noisy Speech 6. Consider filtering noisy speech signals using a 65-tap (64th-order) FIR linear phase lowpass filter of the form = b k e jk He j The filter will be designed using the MATLAB function M k = 0 >> b = fir1(64,wn); % Wn = cutoff frequency relative to pi <-> % 1.0, thus wc = pi/2 --> set Wn = 0.5 % If fc = 2 khz & fs = 8000 Hz, wc = 2*fc/fs FYI: This filter design function, part of the Signal Processing Toolbox, uses windowing of a truncated ideal brickwall filter impulse response. The system model is shown in Figure 4. In (12) dn 8 ksps speech wave files from zip package + wn xn Noise corrupted speech signal FIR Lowpass 65 Taps f c yn = dˆ Lowpass filter with variable cutoff frequency to achieve minimum mean-square error Figure 4: Noisy speech signal and a MMSE estimator that uses an FIR lowpass filter. Problems 10
11 practice we only have access to the noisy signal xn. For the initial filter design we have access to dn test vectors to complete the design of a sub-optimum filter for the estimation of dn from xn. a) Plot the magnitude and phase response of this filter for f c = 2 khz, assuming f s = 8 khz. When using freqz() you will want to use the form [H, F] = freqz(b,a,n,fs); % Say N = 1024, Fs = 8000 or 8k You may want to explore the use of the filter view function, fvtool(b,a, fs,8000). Note the passband ripple in db and the stopband attenuation level in the magnitude response. b) The zip package for Set 1p contains two 8 ksps audio speech files (in.wav format) which can be read into the workspace or a MATLAB function via: d1 = wavread('osr_us_000_0030_8k.wav'); % Male voice, column vec d2 = wavread('osr_us_000_0018_8k.wav'); % Female voice, column vec In the block diagram the speech files constitute dn. To listen to these sound files inside MATLAB you can use the sound() function sound(d,8000); % d = d1 or d2, or both together as [d1; d2] We are now going to add white Gaussian noise wn to dn as follows: >> SNRdB = 10.0; % note in later parts the SNR in db will be changed >> var_w = var(d)/10^(snrdb/10); >> w = sqrt(var_w)*randn(size(d)); % Gaussian noise of variance var_w >> x = d + w; % form signal plus noise Pass this signal through the filter of part (a) and plot using subplots both the noisy input x and filtered version of x (y) for just 500 samples. It is OK to use plot() instead of stem() due to the large number of points involved. Note that the signal-to-noise ratio SNR is initially set to 10 db. c) Continuing from part (b), we will now adjust the filter cutoff so as to minimize the mean-square error (MSE) between the noise free input, d, and the filtered noisy output, y. The filter forms the estimate dˆ. To obtain the minimum MSE (MMSE), you will first need to calculate the mean-squared error (MSE) in a meaningful way. Since the long 65-tap FIR filter introduces a mean signal delay of ((65-1)/2 = 32 samples, the input and output vectors need to be trimmed. There is also a transient period present in the filter output. Choose a reasonable value for the trim interval based on the waveform plot of part (b). The MSE can then be found using >> MSE = var(d_trim - y_trim); where d_trim is the trimmed noise-free input and y_trim is the trimmed filter output when noise is present. Note that var() is MATLAB s sample variance estimator. You now need to repeat this experiment for different values of fc in the FIR filter design, to find the approximate MMSE. Collect data using a concatenation of the male and female Problems 11
12 speech vectors plus noise. The intent is obtain a filter design that works for both male and female voices. You should plot MMSE in db (10*log10(MSE)) versus fc in Hz for fc values ranging from [500, 3500] Hz. Create a family of curves using the same range of fc values with SNR fixed at 50, 20, 10, and 5 db. Does it makes sense that the optimum fc value changes with SNR? d) For the case of SNR = 10 db export xn and yn for your optimum fc to a.wav file. ZIP these files up and them to me as part of the assignment requirements. Comment the intelligibility of the pre- and post-filter speech quality. Bibliography/References [1] J. H. McClellan, C.S. Burrus, A.V. Oppenheim, T.W. Parks, R.W. Schafer, and H.W. Schuessler, Computer-Based Exercises for Signal Processing Using MATLAB 5, Prentice Hall, New Jersey, [2] S.K. Mitra, Digital Signal Processing: A Computer Based Approach, third edition, Mc Graw Hill, New York, [3] S.K. Mitra, Digital Signal Processing Laboratory Using MATLAB, McGraw Hill, New York, [4] M.J. Smith and R.M. Mersereau, Introduction to Digital Signal Processing: A Computer Laboratory Textbook, John Wiley, New York, Bibliography/References 12
Basic 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 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 informationGEORGIA INSTITUTE OF TECHNOLOGY. SCHOOL of ELECTRICAL and COMPUTER ENGINEERING. ECE 2026 Summer 2018 Lab #8: Filter Design of FIR Filters
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL and COMPUTER ENGINEERING ECE 2026 Summer 2018 Lab #8: Filter Design of FIR Filters Date: 19. Jul 2018 Pre-Lab: You should read the Pre-Lab section of
More informationElectrical & Computer Engineering Technology
Electrical & Computer Engineering Technology EET 419C Digital Signal Processing Laboratory Experiments by Masood Ejaz Experiment # 1 Quantization of Analog Signals and Calculation of Quantized noise Objective:
More informationLab S-5: DLTI GUI and Nulling Filters. Please read through the information below prior to attending your lab.
DSP First, 2e Signal Processing First Lab S-5: DLTI GUI and Nulling Filters Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab. Verification: The Exercise
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 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 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 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 informationDFT: Discrete Fourier Transform & Linear Signal Processing
DFT: Discrete Fourier Transform & Linear Signal Processing 2 nd Year Electronics Lab IMPERIAL COLLEGE LONDON Table of Contents Equipment... 2 Aims... 2 Objectives... 2 Recommended Textbooks... 3 Recommended
More information1. In the command window, type "help conv" and press [enter]. Read the information displayed.
ECE 317 Experiment 0 The purpose of this experiment is to understand how to represent signals in MATLAB, perform the convolution of signals, and study some simple LTI systems. Please answer all questions
More informationECE Digital Signal Processing
University of Louisville Instructor:Professor Aly A. Farag Department of Electrical and Computer Engineering Spring 2006 ECE 520 - Digital Signal Processing Catalog Data: Office hours: Objectives: ECE
More informationLab 8: Frequency Response and Filtering
Lab 8: Frequency Response and Filtering Pre-Lab and Warm-Up: You should read at least the Pre-Lab and Warm-up sections of this lab assignment and go over all exercises in the Pre-Lab section before going
More informationExperiments #6. Convolution and Linear Time Invariant Systems
Experiments #6 Convolution and Linear Time Invariant Systems 1) Introduction: In this lab we will explain how to use computer programs to perform a convolution operation on continuous time systems and
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 informationECE 5650/4650 Computer Project #3 Adaptive Filter Simulation
ECE 5650/4650 Computer Project #3 Adaptive Filter Simulation This project is to be treated as a take-home exam, meaning each student is to due his/her own work without consulting others. The grading for
More informationSIGNALS AND SYSTEMS LABORATORY 3: Construction of Signals in MATLAB
SIGNALS AND SYSTEMS LABORATORY 3: Construction of Signals in MATLAB INTRODUCTION Signals are functions of time, denoted x(t). For simulation, with computers and digital signal processing hardware, one
More informationELT COMMUNICATION THEORY
ELT 41307 COMMUNICATION THEORY Matlab Exercise #1 Sampling, Fourier transform, Spectral illustrations, and Linear filtering 1 SAMPLING The modeled signals and systems in this course are mostly analog (continuous
More informationBiomedical Signals. Signals and Images in Medicine Dr Nabeel Anwar
Biomedical Signals Signals and Images in Medicine Dr Nabeel Anwar Noise Removal: Time Domain Techniques 1. Synchronized Averaging (covered in lecture 1) 2. Moving Average Filters (today s topic) 3. Derivative
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 informationGeorge Mason University ECE 201: Introduction to Signal Analysis Spring 2017
Assigned: March 7, 017 Due Date: Week of April 10, 017 George Mason University ECE 01: Introduction to Signal Analysis Spring 017 Laboratory Project #7 Due Date Your lab report must be submitted on blackboard
More informationEEO 401 Digital Signal Processing Prof. Mark Fowler
EEO 41 Digital Signal Processing Prof. Mark Fowler Note Set #17.5 MATLAB Examples Reading Assignment: MATLAB Tutorial on Course Webpage 1/24 Folder Navigation Current folder name here Type commands here
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 informationSTANFORD UNIVERSITY. DEPARTMENT of ELECTRICAL ENGINEERING. EE 102B Spring 2013 Lab #05: Generating DTMF Signals
STANFORD UNIVERSITY DEPARTMENT of ELECTRICAL ENGINEERING EE 102B Spring 2013 Lab #05: Generating DTMF Signals Assigned: May 3, 2013 Due Date: May 17, 2013 Remember that you are bound by the Stanford University
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 informationSMS045 - DSP Systems in Practice. Lab 1 - Filter Design and Evaluation in MATLAB Due date: Thursday Nov 13, 2003
SMS045 - DSP Systems in Practice Lab 1 - Filter Design and Evaluation in MATLAB Due date: Thursday Nov 13, 2003 Lab Purpose This lab will introduce MATLAB as a tool for designing and evaluating digital
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 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 informationELEC3104: Digital Signal Processing Session 1, 2013
ELEC3104: Digital Signal Processing Session 1, 2013 The University of New South Wales School of Electrical Engineering and Telecommunications LABORATORY 4: DIGITAL FILTERS INTRODUCTION In this laboratory,
More informationWindow Method. designates the window function. Commonly used window functions in FIR filters. are: 1. Rectangular Window:
Window Method We have seen that in the design of FIR filters, Gibbs oscillations are produced in the passband and stopband, which are not desirable features of the FIR filter. To solve this problem, window
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 informationDSP First. Laboratory Exercise #7. Everyday Sinusoidal Signals
DSP First Laboratory Exercise #7 Everyday Sinusoidal Signals This lab introduces two practical applications where sinusoidal signals are used to transmit information: a touch-tone dialer and amplitude
More informationAdaptive Systems Homework Assignment 3
Signal Processing and Speech Communication Lab Graz University of Technology Adaptive Systems Homework Assignment 3 The analytical part of your homework (your calculation sheets) as well as the MATLAB
More informationSignal Processing. Introduction
Signal Processing 0 Introduction One of the premiere uses of MATLAB is in the analysis of signal processing and control systems. In this chapter we consider signal processing. The final chapter of the
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 informationEEM478-WEEK8 Finite Impulse Response (FIR) Filters
EEM478-WEEK8 Finite Impulse Response (FIR) Filters Learning Objectives Introduction to the theory behind FIR filters: Properties (including aliasing). Coefficient calculation. Structure selection. Implementation
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 informationProblem Set 1 (Solutions are due Mon )
ECEN 242 Wireless Electronics for Communication Spring 212 1-23-12 P. Mathys Problem Set 1 (Solutions are due Mon. 1-3-12) 1 Introduction The goals of this problem set are to use Matlab to generate and
More informationMATLAB SIMULATOR FOR ADAPTIVE FILTERS
MATLAB SIMULATOR FOR ADAPTIVE FILTERS Submitted by: Raja Abid Asghar - BS Electrical Engineering (Blekinge Tekniska Högskola, Sweden) Abu Zar - BS Electrical Engineering (Blekinge Tekniska Högskola, Sweden)
More informationSubtractive Synthesis. Describing a Filter. Filters. CMPT 468: Subtractive Synthesis
Subtractive Synthesis CMPT 468: Subtractive Synthesis Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University November, 23 Additive synthesis involves building the sound by
More informationLab 4 An FPGA Based Digital System Design ReadMeFirst
Lab 4 An FPGA Based Digital System Design ReadMeFirst Lab Summary This Lab introduces a number of Matlab functions used to design and test a lowpass IIR filter. As you have seen in the previous lab, Simulink
More informationTHE HONG KONG POLYTECHNIC UNIVERSITY Department of Electronic and Information Engineering. EIE2106 Signal and System Analysis Lab 2 Fourier series
THE HONG KONG POLYTECHNIC UNIVERSITY Department of Electronic and Information Engineering EIE2106 Signal and System Analysis Lab 2 Fourier series 1. Objective The goal of this laboratory exercise is to
More informationESE531 Spring University of Pennsylvania Department of Electrical and System Engineering Digital Signal Processing
University of Pennsylvania Department of Electrical and System Engineering Digital Signal Processing ESE531, Spring 2017 Final Project: Audio Equalization Wednesday, Apr. 5 Due: Tuesday, April 25th, 11:59pm
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 informationUNIVERSITY OF WARWICK
UNIVERSITY OF WARWICK School of Engineering ES905 MSc Signal Processing Module (2010) AM SIGNALS AND FILTERING EXERCISE Deadline: This is NOT for credit. It is best done before the first assignment. You
More informationProject 2 - Speech Detection with FIR Filters
Project 2 - Speech Detection with FIR Filters ECE505, Fall 2015 EECS, University of Tennessee (Due 10/30) 1 Objective The project introduces a practical application where sinusoidal signals are used to
More informationDSP First. Laboratory Exercise #2. Introduction to Complex Exponentials
DSP First Laboratory Exercise #2 Introduction to Complex Exponentials The goal of this laboratory is gain familiarity with complex numbers and their use in representing sinusoidal signals as complex exponentials.
More informationLABORATORY - FREQUENCY ANALYSIS OF DISCRETE-TIME SIGNALS
LABORATORY - FREQUENCY ANALYSIS OF DISCRETE-TIME SIGNALS INTRODUCTION The objective of this lab is to explore many issues involved in sampling and reconstructing signals, including analysis of the frequency
More informationEEL 4350 Principles of Communication Project 2 Due Tuesday, February 10 at the Beginning of Class
EEL 4350 Principles of Communication Project 2 Due Tuesday, February 10 at the Beginning of Class Description In this project, MATLAB and Simulink are used to construct a system experiment. The experiment
More informationUNIVERSITY OF WARWICK
UNIVERSITY OF WARWICK School of Engineering ES905 MSc Signal Processing Module (2004) ASSIGNMENT 1 In this assignment, you will use the MATLAB package. In Part (A) you will design some FIR filters and
More informationLakehead University. Department of Electrical Engineering
Lakehead University Department of Electrical Engineering Lab Manual Engr. 053 (Digital Signal Processing) Instructor: Dr. M. Nasir Uddin Last updated on January 16, 003 1 Contents: Item Page # Guidelines
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 informationDesign Digital Non-Recursive FIR Filter by Using Exponential Window
International Journal of Emerging Engineering Research and Technology Volume 3, Issue 3, March 2015, PP 51-61 ISSN 2349-4395 (Print) & ISSN 2349-4409 (Online) Design Digital Non-Recursive FIR Filter by
More informationDigital Filter Design using MATLAB
Digital Filter Design using MATLAB Dr. Tony Jacob Department of Electronics and Electrical Engineering Indian Institute of Technology Guwahati April 11, 2015 Dr. Tony Jacob IIT Guwahati April 11, 2015
More informationECE 2713 Homework 7 DUE: 05/1/2018, 11:59 PM
Spring 2018 What to Turn In: ECE 2713 Homework 7 DUE: 05/1/2018, 11:59 PM Dr. Havlicek Submit your solution for this assignment electronically on Canvas by uploading a file to ECE-2713-001 > Assignments
More informationGeorge Mason University Signals and Systems I Spring 2016
George Mason University Signals and Systems I Spring 2016 Laboratory Project #4 Assigned: Week of March 14, 2016 Due Date: Laboratory Section, Week of April 4, 2016 Report Format and Guidelines for Laboratory
More informationLab S-4: Convolution & FIR Filters. Please read through the information below prior to attending your lab.
DSP First, 2e Signal Processing First Lab S-4: Convolution & FIR Filters Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab. Verification: The Exercise section
More informationGEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL and COMPUTER ENGINEERING. ECE 2025 Fall 1999 Lab #7: Frequency Response & Bandpass Filters
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL and COMPUTER ENGINEERING ECE 2025 Fall 1999 Lab #7: Frequency Response & Bandpass Filters Date: 12 18 Oct 1999 This is the official Lab #7 description;
More informationLab 6: Sampling, Convolution, and FIR Filtering
Lab 6: Sampling, Convolution, and FIR Filtering Pre-Lab and Warm-Up: You should read at least the Pre-Lab and Warm-up sections of this lab assignment and go over all exercises in the Pre-Lab section prior
More informationEE 5410 Signal Processing
EE 54 Signal Processing MATLAB Exercise Telephone Touch-Tone Signal Encoding and Decoding Intended Learning Outcomes: On completion of this MATLAB laboratory exercise, you should be able to Generate and
More informationLaboratory Assignment 4. Fourier Sound Synthesis
Laboratory Assignment 4 Fourier Sound Synthesis PURPOSE This lab investigates how to use a computer to evaluate the Fourier series for periodic signals and to synthesize audio signals from Fourier series
More informationECE 4213/5213 Homework 10
Fall 2017 ECE 4213/5213 Homework 10 Dr. Havlicek Work the Projects and Questions in Chapter 7 of the course laboratory manual. For your report, use the file LABEX7.doc from the course web site. Work these
More informationInterpolated Lowpass FIR Filters
24 COMP.DSP Conference; Cannon Falls, MN, July 29-3, 24 Interpolated Lowpass FIR Filters Speaker: Richard Lyons Besser Associates E-mail: r.lyons@ieee.com 1 Prototype h p (k) 2 4 k 6 8 1 Shaping h sh (k)
More informationL A B 3 : G E N E R A T I N G S I N U S O I D S
L A B 3 : G E N E R A T I N G S I N U S O I D S NAME: DATE OF EXPERIMENT: DATE REPORT SUBMITTED: 1/7 1 THEORY DIGITAL SIGNAL PROCESSING LABORATORY 1.1 GENERATION OF DISCRETE TIME SINUSOIDAL SIGNALS IN
More informationThe Formula for Sinusoidal Signals
The Formula for I The general formula for a sinusoidal signal is x(t) =A cos(2pft + f). I A, f, and f are parameters that characterize the sinusoidal sinal. I A - Amplitude: determines the height of the
More informationFINITE IMPULSE RESPONSE (FIR) FILTERS
CHAPTER 5 FINITE IMPULSE RESPONSE (FIR) FILTERS This chapter introduces finite impulse response (FIR) digital filters. Several methods for designing FIR filters are covered. The Filter Design and Analysis
More informationProject I: Phase Tracking and Baud Timing Correction Systems
Project I: Phase Tracking and Baud Timing Correction Systems ECES 631, Prof. John MacLaren Walsh, Ph. D. 1 Purpose In this lab you will encounter the utility of the fundamental Fourier and z-transform
More information1 Introduction and Overview
DSP First, 2e Lab S-0: Complex Exponentials Adding Sinusoids Signal Processing First Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab. Verification: The
More informationASN Filter Designer Professional/Lite Getting Started Guide
ASN Filter Designer Professional/Lite Getting Started Guide December, 2011 ASN11-DOC007, Rev. 2 For public release Legal notices All material presented in this document is protected by copyright under
More informationSampling and Reconstruction of Analog Signals
Sampling and Reconstruction of Analog Signals Chapter Intended Learning Outcomes: (i) Ability to convert an analog signal to a discrete-time sequence via sampling (ii) Ability to construct an analog signal
More informationGeorge Mason University ECE 201: Introduction to Signal Analysis
Due Date: Week of May 01, 2017 1 George Mason University ECE 201: Introduction to Signal Analysis Computer Project Part II Project Description Due to the length and scope of this project, it will be broken
More informationIslamic University of Gaza. Faculty of Engineering Electrical Engineering Department Spring-2011
Islamic University of Gaza Faculty of Engineering Electrical Engineering Department Spring-2011 DSP Laboratory (EELE 4110) Lab#4 Sampling and Quantization OBJECTIVES: When you have completed this assignment,
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 informationSignal Processing Toolbox
Signal Processing Toolbox Perform signal processing, analysis, and algorithm development Signal Processing Toolbox provides industry-standard algorithms for analog and digital signal processing (DSP).
More informationDecoding a Signal in Noise
Department of Electrical & Computer Engineering McGill University ECSE-490 DSP Laboratory Experiment 2 Decoding a Signal in Noise 2.1 Purpose Imagine that you have obtained through some, possibly suspect,
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 informationArmstrong Atlantic State University Engineering Studies MATLAB Marina Sound Processing Primer
Armstrong Atlantic State University Engineering Studies MATLAB Marina Sound Processing Primer Prerequisites The Sound Processing Primer assumes knowledge of the MATLAB IDE, MATLAB help, arithmetic operations,
More informationUniversiti Malaysia Perlis EKT430: DIGITAL SIGNAL PROCESSING LAB ASSIGNMENT 1: DISCRETE TIME SIGNALS IN THE TIME DOMAIN
Universiti Malaysia Perlis EKT430: DIGITAL SIGNAL PROCESSING LAB ASSIGNMENT 1: DISCRETE TIME SIGNALS IN THE TIME DOMAIN Pusat Pengajian Kejuruteraan Komputer Dan Perhubungan Universiti Malaysia Perlis
More informationSystem Identification and CDMA Communication
System Identification and CDMA Communication A (partial) sample report by Nathan A. Goodman Abstract This (sample) report describes theory and simulations associated with a class project on system identification
More informationLecture 3, Multirate Signal Processing
Lecture 3, Multirate Signal Processing Frequency Response If we have coefficients of an Finite Impulse Response (FIR) filter h, or in general the impulse response, its frequency response becomes (using
More informationFilters. Phani Chavali
Filters Phani Chavali Filters Filtering is the most common signal processing procedure. Used as echo cancellers, equalizers, front end processing in RF receivers Used for modifying input signals by passing
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 informationAdditive Synthesis OBJECTIVES BACKGROUND
Additive Synthesis SIGNALS & SYSTEMS IN MUSIC CREATED BY P. MEASE, 2011 OBJECTIVES In this lab, you will construct your very first synthesizer using only pure sinusoids! This will give you firsthand experience
More informationDSP First Lab 08: Frequency Response: Bandpass and Nulling Filters
DSP First Lab 08: Frequency Response: Bandpass and Nulling Filters Pre-Lab and Warm-Up: You should read at least the Pre-Lab and Warm-up sections of this lab assignment and go over all exercises in the
More informationECE 2026 Summer 2016 Lab #08: Detecting DTMF Signals
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL and COMPUTER ENGINEERING ECE 2026 Summer 2016 Lab #08: Detecting DTMF Signals Date: 14 July 2016 Pre-Lab: You should read the Pre-Lab section of the
More informationPart One. Efficient Digital Filters COPYRIGHTED MATERIAL
Part One Efficient Digital Filters COPYRIGHTED MATERIAL Chapter 1 Lost Knowledge Refound: Sharpened FIR Filters Matthew Donadio Night Kitchen Interactive What would you do in the following situation?
More informationELEC3104: Digital Signal Processing Session 1, 2013 LABORATORY 3: IMPULSE RESPONSE, FREQUENCY RESPONSE AND POLES/ZEROS OF SYSTEMS
ELEC3104: Digital Signal Processing Session 1, 2013 The University of New South Wales School of Electrical Engineering and Telecommunications LABORATORY 3: IMPULSE RESPONSE, FREQUENCY RESPONSE AND POLES/ZEROS
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 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 informationUNIVERSITY OF UTAH ELECTRICAL AND COMPUTER ENGINEERING DEPARTMENT
UNIVERSITY OF UTAH ELECTRICAL AND COMPUTER ENGINEERING DEPARTMENT ECE1020 COMPUTING ASSIGNMENT 3 N. E. COTTER MATLAB ARRAYS: RECEIVED SIGNALS PLUS NOISE READING Matlab Student Version: learning Matlab
More informationECE 3500: Fundamentals of Signals and Systems (Fall 2014) Lab 4: Binary Phase-Shift Keying Modulation and Demodulation
ECE 3500: Fundamentals of Signals and Systems (Fall 2014) Lab 4: Binary Phase-Shift Keying Modulation and Demodulation Files necessary to complete this assignment: none Deliverables Due: Before your assigned
More informationCS3291: Digital Signal Processing
CS39 Exam Jan 005 //08 /BMGC University of Manchester Department of Computer Science First Semester Year 3 Examination Paper CS39: Digital Signal Processing Date of Examination: January 005 Answer THREE
More informationSolution Set for Mini-Project #2 on Octave Band Filtering for Audio Signals
EE 313 Linear Signals & Systems (Fall 2018) Solution Set for Mini-Project #2 on Octave Band Filtering for Audio Signals Mr. Houshang Salimian and Prof. Brian L. Evans 1- Introduction (5 points) A finite
More informationSignal Processing. Naureen Ghani. December 9, 2017
Signal Processing Naureen Ghani December 9, 27 Introduction Signal processing is used to enhance signal components in noisy measurements. It is especially important in analyzing time-series data in neuroscience.
More informationExperiment 1 Introduction to MATLAB and Simulink
Experiment 1 Introduction to MATLAB and Simulink INTRODUCTION MATLAB s Simulink is a powerful modeling tool capable of simulating complex digital communications systems under realistic conditions. It includes
More informationSpring 2014 EE 445S Real-Time Digital Signal Processing Laboratory Prof. Evans. Homework #2. Filter Analysis, Simulation, and Design
Spring 2014 EE 445S Real-Time Digital Signal Processing Laboratory Prof. Homework #2 Filter Analysis, Simulation, and Design Assigned on Saturday, February 8, 2014 Due on Monday, February 17, 2014, 11:00am
More informationContents. Introduction 1 1 Suggested Reading 2 2 Equipment and Software Tools 2 3 Experiment 2
ECE363, Experiment 02, 2018 Communications Lab, University of Toronto Experiment 02: Noise Bruno Korst - bkf@comm.utoronto.ca Abstract This experiment will introduce you to some of the characteristics
More informationOctave Functions for Filters. Young Won Lim 2/19/18
Copyright (c) 2016 2018 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published
More informationEE 403: Digital Signal Processing
OKAN UNIVERSITY FACULTY OF ENGINEERING AND ARCHITECTURE 1 EEE 403 DIGITAL SIGNAL PROCESSING (DSP) 01 INTRODUCTION FALL 2012 Yrd. Doç. Dr. Didem Kıvanç Türeli didem.kivanc@okan.edu.tr EE 403: Digital Signal
More informationece 429/529 digital signal processing robin n. strickland ece dept, university of arizona ECE 429/529 RNS
ece 429/529 digital signal processing robin n. strickland ece dept, university of arizona 2007 SPRING 2007 SCHEDULE All dates are tentative. Lesson Day Date Learning outcomes to be Topics Textbook HW/PROJECT
More information