George Mason University ECE 201: Introduction to Signal Analysis
|
|
- Claire Copeland
- 5 years ago
- Views:
Transcription
1 Due Date: Week of May 01, 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 down into two parts. As a result, the project will be worth two labs in terms of your lab grade. This document is Part II of the project. Lab Report Your project report that you will submit by the due date will contain write-ups for what is asked for in both parts of the project. There is an Instructor Verification Sheet that must be filled out and signed by your lab TA or instructor, and turned in to your TA. Honor Code Forgeries and plagiarism are a violation of the honor code and any reasonable suspicion of an honor code violation will be reported. You are allowed to discuss lab exercises with other students, but the submitted work should be original and it should be your own work. 1 Introduction The goal of this project is to introduce a practical application where sinusoidal signals are used to transmit information: a touch-tone dialer. Bandpass FIR filters can be used to extract the information encoded in the waveforms. In Part I of the project, you designed a system that will generate the DTMF signals to dial a telephone number. In this part, you will design and implement the bandpass FIR filters that are required in order to automatically decode the DTMF signals. Your complete system will consist of a touch tone dialer and a DTMF decoder. 2 Background 2.1 Telephone Touch Tone 2 Dialing Telephone touch pads generate dual tone multi frequency (DTMF) signals to dial a telephone. When a key is pressed, the tones of the corresponding column and row (in Fig. 1) are generated and summed, hence dual tone. In Part I of this project, you wrote a function dtmfdial, to implement a DTMF dialer that generates the appropriate dual tone for each key on the pad of a telephone as specified in Fig. 1. The features of your program are: 1. The input to the function is a vector of numbers, each one being between 1 and 12, with 1 9 corresponding to the digits, 10 corresponds to the * key, 11 is the 0 key, and 12 is the # key. 2. The output is a vector containing the DTMF tones, sampled at f s = Hz. The duration of each tone is about 0.2 sec. with a silence of about 0.05 sec. long to separate the DTMF tones. 1 The due date is for a single project report that contains the required information from parts I and II. 2 Touch Tone is a registered trademark
2 freqs 1209 Hz 1336 Hz 1477 Hz 697 Hz Hz Hz Hz * 0 # Figure 1: DTMF encoding table for Touch Tone dialing. When any key is pressed the tones of the corresponding column and row are generated and summed. 3. The frequency information is given as two 4 3 matrices (tone rows and tone cols): one containing the row frequencies, the other the column frequencies. A digit is translated into the correct location in these 4 3 matrices by using quotients and remainders. Your function should create the appropriate tone sequence to dial an arbitrary phone number and, when played through a telephone handset, will be able to dial the phone. 3 Project Part II: DTMF Decoding In this part of the project, you are to design a decoder that will take the dual tone output from youra dialer and determine what key was entered. There are several steps to decoding a DTMF signal: 1. Divide the signal into shorter time segments representing individual key presses. 2. Filter the individual segments to extract the possible frequency components. Bandpass filters can be used to isolate sinusoidal components. 3. Determine which two frequency components are present in each time segment by measuring the size of the output signal from all of the bandpass filters. 4. Determine which key was pressed, 0 9, *, or # by converting frequency pairs back into key names according to Fig. 1. It is possible to decode DTMF signals using a simple FIR filter bank. A filterbank is nothing more than a bank of filters that are connected in parallel. The filter bank in Fig. 2 consists of seven bandpass filters where each filter passes only one of the seven possible DTMF frequencies. The input signal for all the filters is the same DTMF signal. Here is how the system works: 1. When the input to the filter bank is a DTMF signal, the outputs from two of the bandpass filters (BPFs) should be larger than the rest. 2. By determining which two filters produce the largest output, the frequencies may be determined. A good measure of the output levels is the peak value at the filter outputs, because when the BPF is working properly it should pass only one sinusoidal signal and the peak value would be the amplitude of that sinusoid. 3. These frequencies are then used to determine the DTMF code.
3 697 Hz y 1 [n] 770 Hz y 2 [n] x[n] 852 Hz y 3 [n] 1336 Hz y 6 [n] 1477 Hz y 7 [n] Figure 2: Filter bank consisting of bandpass filters which pass frequencies corresponding to the seven DTMF component frequencies listed in Fig Normalized Frequency (x π rad/sample) Figure 3: L-point averaging filters for lengths L = 5, 9, 11, 15, Bandpass Filter Design As we have seen, an L-point averaging filter, h[n] = 1 L 1 δ[n k] L k=0 is a lowpass filter, and the width of the passbands of these filters is inversely proportional to L. For example, shown in Fig. 3 are L-point averaging filters for L = 5, 9, 11, 15, 21. Note that as the length L increases, the width of the passband decreases. To isolate the tones generated by a DTMF dialer, we would like to design a bank of bandpass filters. To do so, recall what we leaned about modulation: multiplying a signal by cosine shifts the spectrum of the signal up in frequency and down in frequency by an amount equal to the frequency of the cosine that is being modulated. Therefore, it should be clear that if we multiply the L-point averaging filter by a cosine of frequency ˆω c, h[n] = 2 L cos (ˆω cn), 0 n < L
4 Normalized Frequency (x π rad/sample) Normalized Frequency (x π rad/sample) Figure 4: The impulse response and frequency response magnitude of a bandpass filter of length L = 21 with a passband centered at ˆω = 0.2π. then the result will be a bandpass filter centered at ˆω = ˆω c. For example, to design a bandpass filter with a center frequency of 0.2π, we would set ˆω c = 0.2π. In this case, the impulse response and frequency response magnitude would be as shown in Fig. 4. As with the lowpass filter, the bandwidth of the bandpass filter is controlled by L; the larger the value of L, the smaller the bandwidth. BANDPASS FILTER (centered at 0.4π) 1 Magnitude STOPBAND 0.2 PASSBAND STOPBAND Frequency (radians) Figure 5: The frequency response of an FIR bandpass filter is shown with its passband andvstopband regions. (a) Design a bandpass filter with L = 40 that will pass a cosine with a frequency ˆω = 0.4π. Make a plot of the frequency response magnitude and phase. (b) Measure the response of the filter (magnitude and phase) at the following frequencies: ˆω = 0.1πk, for k = 10, 9,..., 9, 10. Summarize the values in a table.
5 The passband of the bandpass filter is defined by the region of the frequency response where H(e j ˆω ) is close to its maximum value of one. Typically, the passband width is defined as the length of the frequency region where H(e j ˆω ) is greater than 1/ 2 = Note: you can use MATLAB s find function to locate those frequencies where the magnitude satisfies H(e j ˆω ) (similar to Fig. 5). Similarly, the stopband of the bandpass filter is defined by the region of the frequency response where H(e j ˆω ) is close to zero. Here, the stopband will be defined as the region where H(e j ˆω ) is less than 0.1. (c) Use the plot of the frequency response for the length-40 bandpass filter from part (a), and determine the passband width. Make two other plots of bandpass filters for L = 20 and L = 80, and measure the passband width in both. Then explain how the width of the passband is related to filter length L, i.e., what happens when L is doubled or halved. (d) Comment on the selectivity of the length-80 bandpass filter. Use the frequency response to explain how the filter is able to pass the components at ˆω = ±0.4π, while reducing or rejecting the others at ˆω = 0.1πk for k ±4. (e) If the input signal to the length-40 FIR BPF is x[n] = cos(0.4πn + π/2) + 7 cos(0.5πn + π/4) find (by hand) the output signal (this should be a formula or an equation). (Hint: use the magnitude and phase measurements from part (a).) Comment on the relative amplitudes of the three signal components in the output signal. 3.2 DTMF Decoding A DTMF decoding system needs two pieces: a set of bandpass filters (BPF) to isolate individual frequency components, and a detector to determine whether or not a given component is present. The detector must score each BPF output and determine which two frequencies are most likely to be contained in the DTMF tone. In a practical system where noise and interference are also present, this scoring process is a crucial part of the system design, but we will only work with noise-free signals to understand the basic functionality in the decoding system. To make the whole system work, it will be necessary to write three mfiles: dtrun, dtscore and dtdesign. An additional mfile called dtcut that you will use is available on the course web page, in the lab section. The main mfile should be named dtrun.m. It will call dtdesign.m, dtcut.m, and dtscore.m. The following sections lead you through the process of writing these functions Simple Bandpass Filter Design: dtdesign.m The FIR filters that will be used in the filter bank (Fig. 2) are finite-length cosines of the form: ( ) 2πfb n h[n] = β cos, 0 n L 1 f s where L is the filter length, and f s is the sampling frequency. The constant β gives flexibility for scaling the filter s gain to meet a constraint such as making the maximum value of the frequency response equal to one. The parameter f b defines the frequency location of the passband, e.g., we pick f b = 852 if we want to isolate the 852 Hz component. The bandwidth of the bandpass filter is controlled by L; the larger the value of L, the narrower the bandwidth.
6 function hh = dtdesign(fcent, L, fs) DTDESIGN hh = dtdesign(fcent, L, fs) returns a matrix (L by length(fcent)) where each column is the impulse response of a BPF, one for each frequency in fcent fcent = vector of center frequencies L = length of FIR bandpass filters fs = sampling freq The BPFs must be scaled so that the maximum magnitude of the frequency response is equal to one. Figure 6: Skeleton of the dtdesign.m function. Complete this function with additional lines of code. (a) Devise a strategy for picking the constant β so that the maximum value of the frequency response will be equal to one. Write the one or two lines of MATLAB code that will do this scaling operation in general. There are two approaches here: (i) Mathematical: derive a formula for β from the formula for the frequency response of the BPF. Then use MATLAB to evaluate this closed-form expression for β. (ii) Numerical: let MATLAB measure the peak value of the unscaled frequency response, and then have MATLAB compute β to scale the peak to be one. (b) Complete the mfile dtdesign.m given in Fig. 6. This function should produce all seven bandpass filters needed for the DTMF filter bank system. The function should store the filters in the columns of the matrix hh whose size will be L 7. The rest of this section describes how to design the BPF to satisfy the following design criterion: 3 Filter Requirements The filter length L should be designed so that only one frequency lies within the passband of each BPF and all other DTMF frequencies lie in the stopband. (c) Generate the seven (scaled) bandpass filters with L = 25 and f s = Plot the magnitude of the frequency responses together on one plot (the range 0 ˆω π is sufficient because H(e j ˆω ) is symmetric). Indicate the locations of each of the 7 DTMF frequencies (697, 770, 852, 941, 1209, 1336, and 1477 Hz) on this plot to illustrate whether or not the passbands are narrow enough to separate the DTMF frequency components. Hint: use the hold command and markers as you did in the warm-up. (d) Repeat the previous part with L = 100 and f s = The width of the passband is supposed to vary inversely with the filter length L. Explain whether or not that is true by comparing the length 25 and length 100 cases. (e) Determine the smallest value of L so that the filter design specification is met, and demonstrate that your filters satisfy this design requirement. 3 Recall the definitions of passband and stopbands given in Section 3.1.
7 3.2.2 A Scoring Function: dtscore.m The final objective is decoding a process that requires a binary decision on the presence or absence of the individual tones. In order to make the signal detection an automated process, we need a score function that rates the different possibilities. (a) Complete the dtscore function based on the skeleton given in Fig. 7. The input signal xx to the dtscore function must be a short segment from the DTMF signal. The task of breaking up the signal so that each short segment corresponds to one key is done by the function dtcut prior to calling dtscore. The implementation of the FIR bandpass filter is done with the conv function, but your could also use the filter function. The running time of the convolution function is proportional to the filter length L. Therefore, the filter length L must satisfy two competing constraints: L should be large so that the bandwidth of the BPF is narrow enough to isolate individual frequency components, but making it too large will cause the program to run slowly. Try to make your system work reliably with the smallest possible value for L. function sc = dtscore(xx, hh) DTSCORE usage: sc = dtscore(xx, hh) returns a score based on the max amplitude of the filtered output xx = input DTMF tone hh = impulse response of ONE bandpass filter The signal detection is done by filtering xx with a length-l BPF, hh, and then finding the maximum amplitude of the output. The score is either 1 or 0. sc = 1 if max( y[n] ) is greater than or equal to 0.95 sc = 0 if max( y[n] ) is less than 0.95 xx = xx*(2/max(abs(xx))); ---Scale x[n] to the range [-2,+2] Figure 7: Skeleton of the dtscore.m function. Complete this function with additional lines of code. (b) Use the following rule for scoring: the score equals one when max y i[n] 0.95; otherwise, it is zero. The signal y i [n] is the output of the i-th BPF. Note: Prior to filtering and scoring, make sure that the input signal x[n] is normalized to the range [ 2, +2]. With this scaling the two sinusoids that make up x[n] should each have amplitudes of approximately Therefore the scoring threshold of 0.95 corresponds to a 95 level for detecting the presence of one sinusoid. Be careful to scale only x[n] and not y i [n]. (c) The scoring rule above depends on proper scaling of the frequency response of the bandpass filters. Explain why the maximum value of the magnitude for H(e j ˆω ) must be equal to one. Consider the fact that both sinusoids in the DTMF tone will experience a known gain (or attenuation) through the bandpass filter, so the amplitude of the output can be predicted if we control both the frequency response and the amplitude of the input. 4 The two sinusoids in a DTMF tone have frequencies that are not harmonics. When plotted versus time, the peaks of the two sinusoids will eventually line up. n
8 Note: When debugging your program it might be useful to have a plot command inside the dtscore.m function. If you plot the first points of the filtered output, you should be able to see two cases: either y[n] is a strong sinusoid with an amplitude close to one (when the filter is matched to one of the component frequencies), or y[n] is relatively small when the filter passband and input signal frequency are mismatched. 3.3 DTMF Decode Function: dtrun.m The DTMF decoding function, dtrun must use information from dtscore to determine which key was pressed based on an input DTMF tone. The skeleton of this function in Fig. 8 includes the help comments. function keys = dtrun(xx,l,fs) DTRUN keys = dtrun(xx,l,fs) returns the list of key numbers found in xx. xx = DTMF waveform L = filter length fs = sampling freq freqs = [697,770,852,941,1209,1336,1477]; hh = dtdesign( freqs,l,fs ); hh = L by 7 MATRIX of all the filters. Each column contains the impulse response of one BPF (bandpass filter) [nstart,nstop] = dtcut(xx,fs); <--Find the tone bursts keys = []; for kk=1:length(nstart) x_seg = xx(nstart(kk):nstop(kk)); <--Extract one DTMF tone.. <=========================================FILL IN THE CODE HERE. end Figure 8: Skeleton of dtrun.m. Complete the for loop in this function with additional lines of code. The function dtrun should work as follows: (i) First, dtrun will design the seven bandpass filters that are needed. (ii) Next, the input signal that is created with your DTMF dial function needs to be broken into individual segments. This is done with a call from dtrun to dtcut, a program that has been provided to you, (iii) For each segment, dtrun will make a call to your dtscore function to score the different BPF outputs. The score returned from dtscore must be either a 1 or a 0 for each frequency. (iv) The final step is to determine the key for that segment (You must add the logic to decide which key is present). (v) The final output is the list of decoded keys. The decoding (step iv) works as follows: If exactly one row frequency and one column frequency are scored as 1 s, then an unique key is identified and the decoding is probably successful. In this case, you can determine the key by using the row and column index. It is possible that there might
9 be an error in scoring if too many or too few frequencies are scored as 1 s. In this case, you should return an error indicator (perhaps by setting the key equal to 1). There are several ways to write the dtrun function, but you should avoid excessive use of if statements to test all 12 cases. Hint: use MATLAB s logicals (e.g., help find) to implement the tests in a few statements Testing Once you get your system working there should be no errors with a large value of L, but when you try to reduce the filter length, the error indicator (key equal to 1) would tell you that the filter length is getting too small. Run tests to find the minimum value for L that gives reliable operation. In you lab report, describe how you tested the system to get this minimum value for L. 3.4 Telephone Numbers The functions dtdial.m and dtrun.m can be used to test the entire DTMF system as shown in Fig. 9. You could also use random digits (e.g., ceil(11.9*rand(1,22)+0.09)) in place of 1:12 in dtdial. For the dtrun function to work correctly, all the mfiles must be on the MATLAB path. >>fs = 11025; <--use this sampling rate in all functions >>xx = dtdial( 1:12, fs ); >>soundsc(xx, fs) >>L = 201; <--overkill, this filter length is way too long >>dtrun(xx, L, fs) ans = Figure 9: Testing the DTMF system. It is also essential to have short pauses in between the tone pairs so that dtcut can parse out the individual signal segments. If you are presenting this project in a lab report, demonstrate a working version of your programs by running it on the following phone number: In addition, make a spectrogram of the signal from dtdial to illustrate the presence of the dual tones. 3.5 Demo When you submit your lab report, you must demonstrate your work to your TA. Get your code and files ready for the demo. You should call dtrun for a signal xx provided by your TA. The output should be the decoded telephone number. The evaluation criteria is shown at the end of the verification sheet.
10 Lab Project ECE-201 Spring INSTRUCTOR VERIFICATION PAGE Turn in this sheet to your TA. Name: Date of Lab: Part 3.1 in Part I of the project: Explain the row-column indexing into frequency table: Write the table below, and then write down the playing order of the frequencies. Determine the total length of the soundsc() signal in seconds. Verified: Date/Time: Part 3.3 in Part I of the project: Complete the dialing function dtmfdial.m: Verified: Date/Time: DTMF Decoding Evaluation Does the designed DTMF decoder decodes the telephone numbers correctly for the following values of L? L=30 All Numbers Most None L=53 All Numbers Most None L=75 All Numbers Most None
11 Appendix function [nstart,nstop] = dtcut(xx,fs) DTCUT Find the DTMF tones within x[n] usage: [nstart,nstop] = dtcut(xx,fs) xx = input signal vector fs = sampling frequency length of nstart = M = number of tones found nstart is the set of STARTING indices nstop is the set of ENDING indices Looks for regions of silence at least 10 millisecs long with tones longer than 100 msec setpoint = 0.02; zero everything below 2 zero xx = xx(:) / max(abs(xx)); normalize xx Lx = length(xx); Lz = round(0.01*fs); xx = filter(ones(1,lz) / Lz, 1, abs(xx)); xx = diff(xx > setpoint); jkl = find(xx ~= 0) ; xx(jkl); if xx(jkl(1)) < 0, jkl = [1; jkl]; end if xx(jkl(end)) > 0, jkl = [jkl; Lx]; end jkl ; indx = []; while length(jkl)>1 if jkl(2)>(jkl(1)+10*lz) indx = [indx, jkl(1:2)]; end jkl(1:2) = []; end nstart = indx(1,:); nstop = indx(2,:);
ECE 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 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 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 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 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 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 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 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 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 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 informationSignal Processing First Lab 20: Extracting Frequencies of Musical Tones
Signal Processing First Lab 20: Extracting Frequencies of Musical Tones 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
More informationDSP First. Laboratory Exercise #11. Extracting Frequencies of Musical Tones
DSP First Laboratory Exercise #11 Extracting Frequencies of Musical Tones This lab is built around a single project that involves the implementation of a system for automatically writing a musical score
More informationSignal Processing First Lab 02: Introduction to Complex Exponentials Multipath. x(t) = A cos(ωt + φ) = Re{Ae jφ e jωt }
Signal Processing First Lab 02: Introduction to Complex Exponentials Multipath 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
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 informationGEORGIA INSTITUTE OF TECHNOLOGY. SCHOOL of ELECTRICAL and COMPUTER ENGINEERING
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL and COMPUTER ENGINEERING ECE 2026 Summer 2018 Lab #3: Synthesizing of Sinusoidal Signals: Music and DTMF Synthesis Date: 7 June. 2018 Pre-Lab: You should
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 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 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 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 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 informationFaculty of Engineering Electrical Engineering Department Communication Engineering I Lab (EELE 3170) Eng. Adam M. Hammad
Faculty of Engineering Electrical Engineering Department Communication Engineering I Lab (EELE 3170) Eng. Adam M. Hammad EXPERIMENT #2 UNDERSTANDING TELEPHONE BASICS Telephone components: 1. Handset containing
More informationSignal Processing First Lab 02: Introduction to Complex Exponentials Direction Finding. x(t) = A cos(ωt + φ) = Re{Ae jφ e jωt }
Signal Processing First Lab 02: Introduction to Complex Exponentials Direction Finding Pre-Lab and Warm-Up: You should read at least the Pre-Lab and Warm-up sections of this lab assignment and go over
More informationHere are some of Matlab s complex number operators: conj Complex conjugate abs Magnitude. Angle (or phase) in radians
Lab #2: Complex Exponentials Adding Sinusoids Warm-Up/Pre-Lab (section 2): You may do these warm-up exercises at the start of the lab period, or you may do them in advance before coming to the lab. You
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 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 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 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 information1 Introduction and Overview
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL and COMPUTER ENGINEERING ECE 2026 Summer 2018 Lab #2: Using Complex Exponentials Date: 31 May. 2018 Pre-Lab: You should read the Pre-Lab section of
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 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 informationDSP First Lab 4a: Synthesis of Sinusoidal Signals Speech Synthesis
DSP First Lab 4a: Synthesis of Sinusoidal Signals Speech Synthesis FORMAL Lab Report: You must write a formal lab report that describes your system for speech synthesis (Section 4). This lab report will
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 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 informationDSP First Lab 03: AM and FM Sinusoidal Signals. We have spent a lot of time learning about the properties of sinusoidal waveforms of the form: k=1
DSP First Lab 03: AM and FM Sinusoidal Signals 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
More informationECE 201: Introduction to Signal Analysis
ECE 201: Introduction to Signal Analysis Prof. Paris Last updated: October 9, 2007 Part I Spectrum Representation of Signals Lecture: Sums of Sinusoids (of different frequency) Introduction Sum of Sinusoidal
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 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 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 informationLab S-2: Direction Finding: Time-Difference or Phase Difference
DSP First, 2e Signal Processing First Lab S-2: Direction Finding: Time-Difference or Phase Difference Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab. Verification:
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 information6.02 Practice Problems: Modulation & Demodulation
1 of 12 6.02 Practice Problems: Modulation & Demodulation Problem 1. Here's our "standard" modulation-demodulation system diagram: at the transmitter, signal x[n] is modulated by signal mod[n] and the
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 informationLab P-4: AM and FM Sinusoidal Signals. We have spent a lot of time learning about the properties of sinusoidal waveforms of the form: ) X
DSP First, 2e Signal Processing First Lab P-4: AM and FM Sinusoidal Signals 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
More informationLab S-8: Spectrograms: Harmonic Lines & Chirp Aliasing
DSP First, 2e Signal Processing First Lab S-8: Spectrograms: Harmonic Lines & Chirp Aliasing Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab. Verification:
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 informationECE 3500: Fundamentals of Signals and Systems (Fall 2015) Lab 4: Binary Phase-Shift Keying Modulation and Demodulation
ECE 500: Fundamentals of Signals and Systems (Fall 2015) Lab 4: Binary Phase-Shift Keying Modulation and Demodulation Files necessary to complete this assignment: none Deliverables Due: Before Dec. 18th
More informationLecture 7 Frequency Modulation
Lecture 7 Frequency Modulation Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/3/15 1 Time-Frequency Spectrum We have seen that a wide range of interesting waveforms can be synthesized
More informationFall Music 320A Homework #2 Sinusoids, Complex Sinusoids 145 points Theory and Lab Problems Due Thursday 10/11/2018 before class
Fall 2018 2019 Music 320A Homework #2 Sinusoids, Complex Sinusoids 145 points Theory and Lab Problems Due Thursday 10/11/2018 before class Theory Problems 1. 15 pts) [Sinusoids] Define xt) as xt) = 2sin
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 informationLab P-3: Introduction to Complex Exponentials Direction Finding. zvect( [ 1+j, j, 3-4*j, exp(j*pi), exp(2j*pi/3) ] )
DSP First, 2e Signal Processing First Lab P-3: Introduction to Complex Exponentials Direction Finding Pre-Lab and Warm-Up: You should read at least the Pre-Lab and Warm-up sections of this lab assignment
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 informationTHE CITADEL THE MILITARY COLLEGE OF SOUTH CAROLINA. Department of Electrical and Computer Engineering. ELEC 423 Digital Signal Processing
THE CITADEL THE MILITARY COLLEGE OF SOUTH CAROLINA Department of Electrical and Computer Engineering ELEC 423 Digital Signal Processing Project 2 Due date: November 12 th, 2013 I) Introduction In ELEC
More informationLab P-10: Edge Detection in Images: UPC Decoding. Please read through the information below prior to attending your lab.
DSP First, 2e Signal Processing First Lab P-10: Edge Detection in Images: UPC Decoding Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab. Verification: The
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 informationLab S-7: Spectrograms of AM and FM Signals. 2. Study the frequency resolution of the spectrogram for two closely spaced sinusoids.
DSP First, 2e Signal Processing First Lab S-7: Spectrograms of AM and FM Signals Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab. Verification: The Exercise
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 informationFilter Banks I. Prof. Dr. Gerald Schuller. Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany. Fraunhofer IDMT
Filter Banks I Prof. Dr. Gerald Schuller Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany 1 Structure of perceptual Audio Coders Encoder Decoder 2 Filter Banks essential element of most
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 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 informationLab S-3: Beamforming with Phasors. N r k. is the time shift applied to r k
DSP First, 2e Signal Processing First Lab S-3: Beamforming with Phasors 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 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 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 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 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 informationLaboratory Assignment 5 Amplitude Modulation
Laboratory Assignment 5 Amplitude Modulation PURPOSE In this assignment, you will explore the use of digital computers for the analysis, design, synthesis, and simulation of an amplitude modulation (AM)
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 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 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 information1. page xviii, line 23:... conventional. Part of the reason for this...
DSP First ERRATA. These are mostly typos, double words, misspellings, etc. Underline is not used in the book, so I ve used it to denote changes. JMcClellan, February 22, 2002 1. page xviii, line 23:...
More informationECE 5650/4650 MATLAB Project 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, 2011. To work the project you will need access to
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 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 informationLecture 3 Complex Exponential Signals
Lecture 3 Complex Exponential Signals Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/3/1 1 Review of Complex Numbers Using Euler s famous formula for the complex exponential The
More informationPrinciples of Communications ECS 332
Principles of Communications ECS 332 Asst. Prof. Dr. Prapun Suksompong prapun@siit.tu.ac.th 5. Angle Modulation Office Hours: BKD, 6th floor of Sirindhralai building Wednesday 4:3-5:3 Friday 4:3-5:3 Example
More informationSignals and Systems Lecture 9 Communication Systems Frequency-Division Multiplexing and Frequency Modulation (FM)
Signals and Systems Lecture 9 Communication Systems Frequency-Division Multiplexing and Frequency Modulation (FM) April 11, 2008 Today s Topics 1. Frequency-division multiplexing 2. Frequency modulation
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 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 informationLecture 4 Frequency Response of FIR Systems (II)
EE3054 Signals and Systems Lecture 4 Frequency Response of FIR Systems (II Yao Wang Polytechnic University Most of the slides included are extracted from lecture presentations prepared by McClellan and
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 informationFigure 1: Block diagram of Digital signal processing
Experiment 3. Digital Process of Continuous Time Signal. Introduction Discrete time signal processing algorithms are being used to process naturally occurring analog signals (like speech, music and images).
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 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 informationSuggested Solutions to Examination SSY130 Applied Signal Processing
Suggested Solutions to Examination SSY13 Applied Signal Processing 1:-18:, April 8, 1 Instructions Responsible teacher: Tomas McKelvey, ph 81. Teacher will visit the site of examination at 1:5 and 1:.
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 informationLab S-1: Complex Exponentials Source Localization
DSP First, 2e Signal Processing First Lab S-1: Complex Exponentials Source Localization Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab. Verification: The
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 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 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 informationElectronics Design Laboratory Lecture #4. ECEN 2270 Electronics Design Laboratory
Electronics Design Laboratory Lecture #4 Electronics Design Laboratory 1 Part A Experiment 2 Robot DC Motor Measure DC motor characteristics Develop a Spice circuit model for the DC motor and determine
More informationLaboratory 5: Spread Spectrum Communications
Laboratory 5: Spread Spectrum Communications Cory J. Prust, Ph.D. Electrical Engineering and Computer Science Department Milwaukee School of Engineering Last Update: 19 September 2018 Contents 0 Laboratory
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 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 informationFinal Exam Solutions June 14, 2006
Name or 6-Digit Code: PSU Student ID Number: Final Exam Solutions June 14, 2006 ECE 223: Signals & Systems II Dr. McNames Keep your exam flat during the entire exam. If you have to leave the exam temporarily,
More informationAC CURRENTS, VOLTAGES, FILTERS, and RESONANCE
July 22, 2008 AC Currents, Voltages, Filters, Resonance 1 Name Date Partners AC CURRENTS, VOLTAGES, FILTERS, and RESONANCE V(volts) t(s) OBJECTIVES To understand the meanings of amplitude, frequency, phase,
More informationEECS 452 Midterm Exam (solns) Fall 2012
EECS 452 Midterm Exam (solns) Fall 2012 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 I /40 Section
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 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 informationOutline. Communications Engineering 1
Outline Introduction Signal, random variable, random process and spectra Analog modulation Analog to digital conversion Digital transmission through baseband channels Signal space representation Optimal
More informationEE 400L Communications. Laboratory Exercise #7 Digital Modulation
EE 400L Communications Laboratory Exercise #7 Digital Modulation Department of Electrical and Computer Engineering University of Nevada, at Las Vegas PREPARATION 1- ASK Amplitude shift keying - ASK - in
More informationCharan Langton, Editor
Charan Langton, Editor SIGNAL PROCESSING & SIMULATION NEWSLETTER Baseband, Passband Signals and Amplitude Modulation The most salient feature of information signals is that they are generally low frequency.
More information