# Lab S-7: Spectrograms of AM and FM Signals. 2. Study the frequency resolution of the spectrogram for two closely spaced sinusoids.

Save this PDF as:

Size: px
Start display at page:

Download "Lab S-7: Spectrograms of AM and FM Signals. 2. Study the frequency resolution of the spectrogram for two closely spaced sinusoids."

## Transcription

2 1.3 Amplitude Modulation If we add several sinusoids, each with a different frequency.f k /, we cannot use the phasor addition theorem, but we can still express the result as a summation of terms with complex amplitudes via: ( NX N ) X x.t/ D A k cos.2f k t C ' k / D < A k e j' k e j 2f kt (2) kd1 where A k e j' k is the complex amplitude of the k th complex exponential term. The choice of f k will determine the nature of the signal for amplitude modulation or beat signals we pick two or three frequencies that are very close together, see Chapter 3. kd Beat Control GUI To assist you in your experiments with beat notes and AM signals, the MATLAB GUI tool called beatcon has been created. The M-file beatcon.m is part of the DSP-First (or SP-First) Toolbox. This user interface controller will exhibit the basic signal shapes for beat signals and play the signals. A small control panel will appear on the screen with buttons and sliders that vary the different parameters for the beat signals. It can also call a user-written function named beat.m. Experiment with the beatcon control panel and use it to produce a beat signal with two frequency components: one at 690 Hz and the other at 700 Hz. Use a longer duration than the default to hear the beat frequency sound. 1.4 Frequency Modulated Signals In this lab, we will examine signals whose frequency content varies as a function of time. Recall that in a constant-frequency sinusoid (1) the argument of the cosine is.2f 0 t C '/ which is also the exponent of the complex exponential. We define the argument of the cosine as the angle function. In (1), the angle function changes linearly versus time, and its time derivative, 2f 0, equals the constant frequency of the cosine. A generalization is available if we adopt the following notation for the class of signals with time-varying angle functions: x.t/ D Acos..t// D <fae j.t/ g (3) where.t/ is the angle function. 1 The time derivative of the angle function.t/ in (3) gives a frequency that we call the instantaneous (radian) frequency:! i.t/ D d dt.t/ (rad/s) If we prefer units in hertz, then we divide by 2 to define the instantaneous (cyclic) frequency: f i.t/ D 1 2 d dt.t/ (Hz) (4) 1.5 Chirp, or Linearly Swept Frequency A linear-fm chirp signal is a sinusoid whose instantaneous frequency changes linearly from a starting value to an ending one. 2 The formula for such a signal can be defined by creating a complex exponential signal with a quadratic angle function.t/. Mathematically, we define.t/ in (3) as.t/ D 2t 2 C 2f 0 t C ' 1 Look for the demo FM Synthesis on the companion website. 2 Look for the demo Spectrograms & Sounds: Wideband FM on the companion website. 2 McClellan, Schafer and Yoder, Signal Processing First.

3 The derivative of.t/ yields an instantaneous cyclic frequency (4) that changes linearly versus time. f i.t/ D 2t C f 0 (hertz) (5) The slope of f i.t/ is equal to 2 and its t D 0 intercept is f 0. The frequency variation in (5) produced by the time-varying angle function is called frequency modulation, so these signals are called FM signals. Finally, since the linear variation of the frequency (5) can produce an audible sound similar to a siren or a bird chirp, linear-fm signals are also called chirps. If the signal starts at time t D t 1 s with a frequency of f 1 Hz, and ends at time t D t 2 s with a frequency of f 2 Hz, then the slope of the line in (5) will be SLOPE D 2 D f 2 f 1 t 2 t 1 (6) Note that if the signal starts at time t D 0 s, then f 1 D f 0 is also the starting frequency. Otherwise, f 0 D 1.6 MATLAB Synthesis of Chirp Signals In MATLAB signals can only be synthesized by evaluating the signal s defining formula at discrete instants of time. These are called samples of the signal. For the chirp we use the following: x.t n / D Acos.2t 2 n C 2f 0 t n C '/ where t n is the n th time sample. The following MATLAB code will synthesize a linear-fm chirp: fsamp = 8000; dt = 1/fSamp; tstart = 0; tstop = 1.5; tt = tstart:dt:tstop; mu = 600; fzero = 400; phi = 2*pi*rand; %-Number of time samples per second %-- random phase % %% psi =????; <=================== FILL IN THE CODE HERE % cc = real( 7.7*exp(j*psi) ); % soundsc( cc, fsamp ); %-- uncomment to hear the sound plotspec( cc+j*1e-12, fsamp, 256 ), colorbar, grid on %-- with negative frequencies (a) Determine the total duration of this synthesized signal in seconds, and also the length of the tt vector. Use MATLAB s size command to check that the signal vector cc has the expected size. (b) Determine the range of frequencies (in hertz) that will be synthesized by the MATLAB script above, i.e., determine the minimum and maximum frequencies (in Hz) that will be heard. This will require that you relate the parameters, f 0, and ' to the minimum and maximum frequencies. Examine the MATLAB spectrogram to determine the instantaneous (cyclic) frequency f i.t/ versus time. Zoom in to verify the correct starting and ending frequencies. (c) The spectrogram usually shows only the frequency components for f 0, but with the tiny imaginary part trick plotspec will show the negative frequency components. We will called this a twosided spectrogram. Since the chirp signal is real-valued, the spectrum must have conjugate symmetry, so the magnitudes of the negative frequency components are a mirror image of those in the positive frequency region. 3 McClellan, Schafer and Yoder, Signal Processing First.

4 (d) Use soundsc() to listen to the signal in order to determine whether the signal s frequency content is increasing or decreasing. Notice that soundsc() needs to know two things: the vector containing the signal samples, and the rate at which the signal samples are to be played out. This rate should be the same as the rate at which the signal values were created (fsamp in the code above). For more information do help sound and help soundsc in MATLAB. (e) The test case above generates a chirp sound whose frequency starts low and chirps up. Modify the parameters so that the chirp starts at 3500 Hz and chirps down to 500 Hz. 1.7 Spectrogram of an FM Signal: Sinusoidal Modulation Define an FM signal whose instantaneous frequency is sinusoidal, i.e.,! i.t/ D 2f c C 2 cos.2ˇt C / rad/s (7) where f c is the center frequency, and the parameters, ˇ and control the sinusoidal frequency modulation. (a) Determine the mathematical formula for an FM signal that has the instantaneous frequency in (7). (b) Write a MATLAB function (or script) to create sinusoidal FM signals of the form defined in (7). Modify the code in Sect. 1.6 to use the parameters in (7). If you choose to make a function, the MATLAB function should use the following template: makesinusfmvals( alpha, beta, gamma, fc, fsamp, tstart, tstop ). (c) Create a sinusoidal-fm signal with f c D 100 Hz, D 50, ˇ D 1:5, and D =3. Make the signal duration equal to 3.04 secs, starting at t D 0. Use a sampling rate of 1000 samples/s. The signal amplitude should be one. (d) Create a spectrogram of this chirp signal, and use it to verify that you have the correct instantaneous frequency predicted by (7). The section length should be short enough to track the changing instantaneous frequency. 1.8 Review Topic: Spectrograms The main issue in this lab will be the dependence of the spectrogram on the choice of section length. A spectrogram is formed by taking successive short sections of a signal and performing an FFT analysis of each of those sections to get the spectrum. Since this is done repeatedly, the result is the spectrum versus time, where time is the location of the short sections. For a specific example, assume that the section length is 100, and the signal is a MATLAB vector xx. Then the first short section will be xx(1:100). The sections are usually overlapped and the default in plotspec is 50% overlap, so the second short section is xx(51:150), the third xx(101:200), and so on. The spectrogram image is, in effect, the spectrum versus time, so we need a reference time for each short section. In plotspec this reference time is the midpoint of the section. For the length-100 section, the reference index is 50, which is then converted to a time (in secs) by using the sampling rate.f s /. When the spectrogram is displayed as an image, these reference times are used along the horizontal axis. For more information refer to the write-up in a previous lab. 4 McClellan, Schafer and Yoder, Signal Processing First.

5 2 Lab Exercise For the lab exercise, you will synthesize some AM and FM signals, and then verify that these signals have the correct frequency content by using the spectrogram. The objective is to learn enough to be able to discuss the connection between the time-domain definition of the signal and its frequency-domain content. For the instructor verification, you will have to demonstrate that you understand concepts in a given subsection by answering questions from your lab instructor (or TA). 2.1 Beat Notes and Frequency Resolution In the section on beat notes in Chapter 3 of the text, we discussed signals formed as the product of two sinusoidal signals of slightly different frequencies; i.e., x.t/ D B cos.2f t C ' /cos.2f c t C ' c / (8) where f c is the (high) center frequency, and f is the (low) frequency that modulates the envelope of the signal. An equivalent representation for the beat signal is obtained by rewriting the product as a sum: x.t/ D A 1 cos.2f 1 t C ' 1 / C A 2 cos.2f 2 t C ' 2 / (9) It is relatively easy to derive the relationship between the frequencies ff 1 ; f 2 g and ff c ; f g MATLAB Code for Beat Signals A beat signal is defined by five parameters fb; f c ; f ; ' c ; ' g along with the start and end times and the sampling rate.f s /, as shown in the following template: Amp = 10; %-- B in equation above fc = 1024; %-- center frequency phic = 2*pi*rand; %-- phase of 2nd sinusoid (random) fdelta = 4; %-- modulating frequency phidelta = 2*pi*rand; %-- phase of 1st sinusoid (random) tstart = 0; %-- starting time (secs) tstop = 5; %-- ending time (secs) fsamp = 8000; % tt = tstart:(1/fsamp):tstop; %-- vector of times xx = Amp*cos(2*pi*fc*tt+phic).*cos(2*pi*fDelta*tt+phiDelta) Beat Note Spectrograms Beat notes have a simple time-frequency characteristic in a spectrogram. Even though a beat note signal, when defined as a product in (8), may be viewed as a single frequency signal whose amplitude varies with time, the spectrum requires an additive combination as in (9) which turns out to be the sum of two sinusoids with different constant frequencies. Beat notes provide an interesting way to investigate the time-frequency characteristics of spectrograms. Although some of the mathematical details require further study beyond this course, it is not difficult to appreciate the following issue: there is a fundamental trade-off between knowing which frequencies are present in a signal s spectrum and knowing how those frequencies vary with time. As discussed previously, a spectrogram estimates the frequency content over short sections of the signal; this is the Section Length parameter. 3 If we make the section length very short we can track rapid changes in the signal, usually changes in the frequency content. The tradeoff, however, is that shorter sections may not provide enough 3 The section length is often called the window length; the two terms are used interchangeably in DSP. 5 McClellan, Schafer and Yoder, Signal Processing First.

6 data to do an accurate frequency measurement. On the other hand, long sections allow the spectrogram to perform excellent frequency measurements, but fail to track sudden frequency changes. For example, if a signal is the sum of two sinusoids whose frequencies are nearly the same, a very long section length is needed to resolve the two sinusoidal components. This trade-off between the section length (in time) and the frequency resolution is akin to Heisenburg s Uncertainty Principle in physics. We can summarize this discussion by stating the following hypothesis: The frequency resolution of the spectrogram is inversely proportional to the Section Length. In other words, when the true spectrum has two lines (at f 1 and f 2 ) these two lines will be visible as distinct lines in the spectrogram if jf 1 f 2 j C=T SECT where C is a proportionality constant and T SECT is the section duration in secs. Note: When using plotspec(xx,fs,lsect), the section length in samples is an input argument to the spectrogram function. We can use the sampling rate to convert to duration, T SECT D L SECT =f s. We will use beat note signals which consist of two closely spaced spectral lines to confirm this hypothesis. A beat note signal may be viewed as a single frequency signal whose amplitude varies with time, or as the sum of two sinusoidal signals with different constant frequencies. Both views can be used to explain the effect of (window) section length when finding the spectrogram of a beat signal. (a) Use the MATLAB code written in Section to create and plot a beat signal defined via: b.t/ D 10cos.2.f /t C ' /cos /t C ' c /; with a duration of 5 s, and a sampling rate of f s D 8000 samples/s. The frequency f should be set to 4 Hz, but will be varied in later parts. The phases can be random. (b) When f D 4 determine the locations of the two spectrum lines that you expect to see in the spectrogram. In other words, derive (mathematically) the spectrum of the signal defined in part (a). (c) Make the spectrogram of b.t/ using a (window) section length of L SECT D 256 using the commands 4 : plotspec(xx,fsamp,256); colorbar, grid on, zoom on Comment on what you see. Are there two spectral lines, i.e., (horizontal lines across the spectrogram)? If necessary, use the zoom tool (in the MATLAB figure window), or zoom on, to examine the important regions of the spectrogram. (d) It should not be possible to see both spectrum lines with L SECT D 256. In order to get both lines a longer section length is needed, so try doubling the section length. Try L SECT D 512, then L SECT D 1024, and so on until you can discern two spectrum lines. 5 Then reduce the value of L SECT little by little to get the smallest L SECT that will work. Getting a value of L SECT to the nearest 500 is sufficient. As before, use zooming to examine the important regions of the spectrogram. Once you have two spectrum lines, record the value of L SECT and determine whether the frequencies present in the spectrogram are correct. In addition, convert L SECT to the section duration in seconds, T SECT. Instructor Verification (separate page) Inverse Relationship: Section Length vs. Frequency Resolution The shortest section length when you are able to discern the two spectrum lines was determined in the previous section. 4 Use plotspec instead of specgram in order to get a linear amplitude scale rather than logarithmic. 5 Usually the window (section) length is chosen to be a power of two, because a special algorithm called the FFT is used in the computation. The fastest FFT programs are those where the FFT length is a power of 2. 6 McClellan, Schafer and Yoder, Signal Processing First.

7 (a) The shortest (window) section length of L SECT samples has been converted into a (window) section duration in seconds (via the sampling rate). Compare the inverse of this (window) section duration to the frequency separation of the spectrum lines. jf 1 f 2 j?! 1 T SECT If we believe that the inverse relationship between (window) section duration and frequency separation is true, then we can calculate a constant C such that From your first experiment, determine C. jf 1 f 2 j D C T SECT (10) (b) Now change f to 16 Hz and repeat the resolution experiment in Sect (d). That is, find a section length that will resolve the two frequency components which are now farther apart. Use the value of C and (10) to predict the section length that you will need. Verify that this section length will work correctly. Note: The relationship is approximate, so the derived section length is not guaranteed to work. If that happens, a small increase in L SECT should make it work. Instructor Verification (separate page) 2.2 Spectrogram for a Chirp with Negative Instantaneous Frequency Use the code provided in the pre-lab section as a starting point in order to write a MATLAB script or function that will synthesize a chirp signal. Then use that M-file in this section. (a) What happens when we make a signal that chirps down and the instantaneous frequency goes negative? Generate a chirp signal that starts at 2000 Hz when t D 0 s, and chirps down to 1000 Hz, at t D 1:5 s. Use f s D 8000 Hz. Determine the parameters needed in (4). (b) Generate the chirp signal in MATLAB and make a spectrogram with L SECT D 200 to verify that you have the correct starting and ending frequencies. For L SECT D 200, determine the section duration T SECT in secs. (c) Will you hear negative frequency? Use a spectrogram that contains negative frequencies to explain your answer. If possible listen to the signal and describe what you will hear; then explain in terms of the two-sided spectrogram. Instructor Verification (separate page) Section Length in Chirp Spectrogram When we have a signal whose frequency is not constant, we can study how the temporal features of the spectrogram depend on the section length L SECT. As we saw in the previous section, the section length has to be short to capture the temporal changes in the signal. However, let s see what happens with a long section length. (a) Generate the same signal as in Sect. 2.2, but make the spectrogram with L SECT D (b) When L SECT D 1600, determine the section duration T SECT in secs. 7 McClellan, Schafer and Yoder, Signal Processing First.

8 (c) The spectrogram uses 50% overlapping and skipping (see Sect. 1.8). Based on T SECT, determine the time locations where the spectrum is being computed. Relate these time locations to what you see in the spectrogram. Instructor Verification (separate page) (d) Optional: Use the slope of the instantaneous frequency to determine how much the frequency changes during one section. Calculate the frequency change from T SECT and. In this spectrogram, the changing frequency appears as rectangular bars that have a measurable width and height. The width should be 0:5T SECT in secs. Compare the calculated frequency change to the vertical height of the bars are they different or equal? 2.3 Lab-HW: Matching Unknown Spectrograms Now you are given a spectrogram in Fig. 1, and you must synthesize a signal that will match that timefrequency plot. Explain and discuss your work. 1. Define a time signal x.t/ whose spectrogram will match the given spectrogram. This signal definition should be a simple mathematical formula. Note: you might have to iterate with the following two steps to get a good approximation. 2. Generate samples of the signal over the appropriate time interval using f s D 4000 Hz. 3. Choose the section length in plotspec carefully so that your spectrogram is an excellent match Freq (Hz) Time (secs) Figure 1: Spectrogram of unknown signal with f s D 4000 Hz. Section length to be determined. Horizontal axis is time in seconds. 8 McClellan, Schafer and Yoder, Signal Processing First.

9 Lab: Spectrograms of AM and FM Signals INSTRUCTOR VERIFICATION SHEET Turn this page in to your lab grading TA before the end of your scheduled Lab time. Name: LoginUserName: Date: Part Record the value of the section length L SECT (in samples) and T SECT (in secs) when you can discern two separate spectral lines for the beat note signal, using f D 4 Hz,. L SECT D T SECT D Verified: Date/Time: Part Calculate the proportionality constant C for the inverse relationship: jf 1 f 2 j C=T SECT. Then determine a new section length L SECT for f D 16 Hz. Synthesize the signal and make its spectrogram with the new section length L SECT. Then verify that the two spectrum lines are resolved. C D L SECT D Verified: Date/Time: Part 2.2 Write MATLAB code for synthesizing a linear-fm chirp whose instantaneous frequency goes negative. Also, display the two-sided spectrogram that includes the negative frequency region, as well as the one-sided spectrogram that has positive frequency components only. Determine the section duration T SECT in secs. T SECTD Verified: Date/Time: Part Spectrogram of Chirp with long-duration section, L SECT D Explain features in the spectrogram, e.g., location of sections along the time axis and duration of sections. T SECT D Section times (centers) = Verified: Date/Time: 9 McClellan, Schafer and Yoder, Signal Processing First.

10 Lab: Spectrograms of AM and FM Signals LAB HOMEWORK QUESTION Turn this page in to your lab grading TA at the very beginning of your next scheduled Lab time. Name: LoginUserName: Date: Part Match the unknown spectrogram. Give the mathematical formula for the signal 2. Include the MATLAB code for generating the signal. 3. Choose a section length to get the desired spectrogram. Include a plot of the spectrogram when you hand in your Lab-HW. 10 McClellan, Schafer and Yoder, Signal Processing First.

### Lab 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

### DSP First. Laboratory Exercise #4. AM and FM Sinusoidal Signals

DSP First Laboratory Exercise #4 AM and FM Sinusoidal Signals The objective of this lab is to introduce more complicated signals that are related to the basic sinusoid. These are signals which implement

### Lab 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

### 1 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

### Lab 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

### Lecture 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

### Digital Signal Processing Lecture 1 - Introduction

Digital Signal Processing - Electrical Engineering and Computer Science University of Tennessee, Knoxville August 20, 2015 Overview 1 2 3 4 Basic building blocks in DSP Frequency analysis Sampling Filtering

### Signal 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

### Lab 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:

### 1 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

### ECE 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

### DSP 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.

### DSP 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

### 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

### Lab P-8: Digital Images: A/D and D/A

DSP First, 2e Signal Processing First Lab P-8: Digital Images: A/D and D/A Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab. Verification: The Warm-up section

### Lab 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

Spectrum Music 7a: Modulation Tamara Smyth, trsmyth@ucsd.edu Department of Music, University of California, San Diego (UCSD) October 3, 7 When sinusoids of different frequencies are added together, the

### DSP First Lab 06: Digital Images: A/D and D/A

DSP First Lab 06: Digital Images: A/D and D/A 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

### Interpolation Error in Waveform Table Lookup

Carnegie Mellon University Research Showcase @ CMU Computer Science Department School of Computer Science 1998 Interpolation Error in Waveform Table Lookup Roger B. Dannenberg Carnegie Mellon University

### Lab 3 FFT based Spectrum Analyzer

ECEn 487 Digital Signal Processing Laboratory Lab 3 FFT based Spectrum Analyzer Due Dates This is a three week lab. All TA check off must be completed prior to the beginning of class on the lab book submission

### DFT: 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

### Reading: Johnson Ch , Ch.5.5 (today); Liljencrants & Lindblom; Stevens (Tues) reminder: no class on Thursday.

L105/205 Phonetics Scarborough Handout 7 10/18/05 Reading: Johnson Ch.2.3.3-2.3.6, Ch.5.5 (today); Liljencrants & Lindblom; Stevens (Tues) reminder: no class on Thursday Spectral Analysis 1. There are

### Armstrong 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,

### Creating Digital Music

Chapter 2 Creating Digital Music Chapter 2 exposes students to some of the most important engineering ideas associated with the creation of digital music. Students learn how basic ideas drawn from the

### Sound Waves and Beats

Sound Waves and Beats Computer 32 Sound waves consist of a series of air pressure variations. A Microphone diaphragm records these variations by moving in response to the pressure changes. The diaphragm

### THE 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

### FFT 1 /n octave analysis wavelet

06/16 For most acoustic examinations, a simple sound level analysis is insufficient, as not only the overall sound pressure level, but also the frequency-dependent distribution of the level has a significant

### Chapter 2. Signals and Spectra

Chapter 2 Signals and Spectra Outline Properties of Signals and Noise Fourier Transform and Spectra Power Spectral Density and Autocorrelation Function Orthogonal Series Representation of Signals and Noise

### Fourier Transform Pairs

CHAPTER Fourier Transform Pairs For every time domain waveform there is a corresponding frequency domain waveform, and vice versa. For example, a rectangular pulse in the time domain coincides with a sinc

### Phase demodulation using the Hilbert transform in the frequency domain

Phase demodulation using the Hilbert transform in the frequency domain Author: Gareth Forbes Created: 3/11/9 Revised: 7/1/1 Revision: 1 The general idea A phase modulated signal is a type of signal which

### Lecture 2: SIGNALS. 1 st semester By: Elham Sunbu

Lecture 2: SIGNALS 1 st semester 1439-2017 1 By: Elham Sunbu OUTLINE Signals and the classification of signals Sine wave Time and frequency domains Composite signals Signal bandwidth Digital signal Signal

### 5.1 Graphing Sine and Cosine Functions.notebook. Chapter 5: Trigonometric Functions and Graphs

Chapter 5: Trigonometric Functions and Graphs 1 Chapter 5 5.1 Graphing Sine and Cosine Functions Pages 222 237 Complete the following table using your calculator. Round answers to the nearest tenth. 2

### Spectrum Analysis: The FFT Display

Spectrum Analysis: The FFT Display Equipment: Capstone, voltage sensor 1 Introduction It is often useful to represent a function by a series expansion, such as a Taylor series. There are other series representations

### Linguistic Phonetics. Spectral Analysis

24.963 Linguistic Phonetics Spectral Analysis 4 4 Frequency (Hz) 1 Reading for next week: Liljencrants & Lindblom 1972. Assignment: Lip-rounding assignment, due 1/15. 2 Spectral analysis techniques There

### Signals. Periodic vs. Aperiodic. Signals

Signals 1 Periodic vs. Aperiodic Signals periodic signal completes a pattern within some measurable time frame, called a period (), and then repeats that pattern over subsequent identical periods R s.

### TRANSFORMS / WAVELETS

RANSFORMS / WAVELES ransform Analysis Signal processing using a transform analysis for calculations is a technique used to simplify or accelerate problem solution. For example, instead of dividing two

### Power Spectral Density (PSD) for TH-UWB signals using PPM is derived in this

C H A P T E R 3 The PSD of TH-UWB Signals Power Spectral Density (PSD) for TH-UWB signals using PPM is derived in this chapter. The adopted approach (Di Benedetto and Vojcic, 3) follows the analog PPM

### CHAPTER 6 INTRODUCTION TO SYSTEM IDENTIFICATION

CHAPTER 6 INTRODUCTION TO SYSTEM IDENTIFICATION Broadly speaking, system identification is the art and science of using measurements obtained from a system to characterize the system. The characterization

### Sound synthesis with Pure Data

Sound synthesis with Pure Data 1. Start Pure Data from the programs menu in classroom TC307. You should get the following window: The DSP check box switches sound output on and off. Getting sound out First,

### ECE 2111 Signals and Systems Spring 2009, UMD Experiment 3: The Spectrum Analyzer

ECE 2111 Signals and Systems Spring 2009, UMD Experiment 3: The Spectrum Analyzer Objective: Student will gain an understanding of the basic controls and measurement techniques of the Rohde & Schwarz Handheld

### G(f ) = g(t) dt. e i2πft. = cos(2πf t) + i sin(2πf t)

Fourier Transforms Fourier s idea that periodic functions can be represented by an infinite series of sines and cosines with discrete frequencies which are integer multiples of a fundamental frequency

### System Inputs, Physical Modeling, and Time & Frequency Domains

System Inputs, Physical Modeling, and Time & Frequency Domains There are three topics that require more discussion at this point of our study. They are: Classification of System Inputs, Physical Modeling,

### SIGNALS 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

### 1. 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:...

### Sinusoids and Sinusoidal Correlation

Laboratory 3 May 24, 2002, Release v3.0 EECS 206 Laboratory 3 Sinusoids and Sinusoidal Correlation 3.1 Introduction Sinusoids are important signals. Part of their importance comes from their prevalence

### Week 15. Mechanical Waves

Chapter 15 Week 15. Mechanical Waves 15.1 Lecture - Mechanical Waves In this lesson, we will study mechanical waves in the form of a standing wave on a vibrating string. Because it is the last week of

### page 7.51 Chapter 7, sections , pp Angle Modulation No Modulation (t) =2f c t + c Instantaneous Frequency 2 dt dt No Modulation

page 7.51 Chapter 7, sections 7.1-7.14, pp. 322-368 Angle Modulation s(t) =A c cos[(t)] No Modulation (t) =2f c t + c s(t) =A c cos[2f c t + c ] Instantaneous Frequency f i (t) = 1 d(t) 2 dt or w i (t)

### C.8 Comb filters 462 APPENDIX C. LABORATORY EXERCISES

462 APPENDIX C. LABORATORY EXERCISES C.8 Comb filters The purpose of this lab is to use a kind of filter called a comb filter to deeply explore concepts of impulse response and frequency response. The

### Local Oscillator Phase Noise and its effect on Receiver Performance C. John Grebenkemper

Watkins-Johnson Company Tech-notes Copyright 1981 Watkins-Johnson Company Vol. 8 No. 6 November/December 1981 Local Oscillator Phase Noise and its effect on Receiver Performance C. John Grebenkemper All

### Performing the Spectrogram on the DSP Shield

Performing the Spectrogram on the DSP Shield EE264 Digital Signal Processing Final Report Christopher Ling Department of Electrical Engineering Stanford University Stanford, CA, US x24ling@stanford.edu

Adaptive Line Enhancer (ALE) This demonstration illustrates the application of adaptive filters to signal separation using a structure called an adaptive line enhancer (ALE). In adaptive line enhancement,

### Experiment 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

### PART I: The questions in Part I refer to the aliasing portion of the procedure as outlined in the lab manual.

Lab. #1 Signal Processing & Spectral Analysis Name: Date: Section / Group: NOTE: To help you correctly answer many of the following questions, it may be useful to actually run the cases outlined in the

### Timbral Distortion in Inverse FFT Synthesis

Timbral Distortion in Inverse FFT Synthesis Mark Zadel Introduction Inverse FFT synthesis (FFT ) is a computationally efficient technique for performing additive synthesis []. Instead of summing partials

### Sampling 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

### DC and AC Circuits. Objective. Theory. 1. Direct Current (DC) R-C Circuit

[International Campus Lab] Objective Determine the behavior of resistors, capacitors, and inductors in DC and AC circuits. Theory ----------------------------- Reference -------------------------- Young

### Poles and Zeros of H(s), Analog Computers and Active Filters

Poles and Zeros of H(s), Analog Computers and Active Filters Physics116A, Draft10/28/09 D. Pellett LRC Filter Poles and Zeros Pole structure same for all three functions (two poles) HR has two poles and

### Integrators, differentiators, and simple filters

BEE 233 Laboratory-4 Integrators, differentiators, and simple filters 1. Objectives Analyze and measure characteristics of circuits built with opamps. Design and test circuits with opamps. Plot gain vs.

### Chapter 4. Digital Audio Representation CS 3570

Chapter 4. Digital Audio Representation CS 3570 1 Objectives Be able to apply the Nyquist theorem to understand digital audio aliasing. Understand how dithering and noise shaping are done. Understand the

### DIGITAL SIGNAL PROCESSING WITH VHDL

DIGITAL SIGNAL PROCESSING WITH VHDL GET HANDS-ON FROM THEORY TO PRACTICE IN 6 DAYS MODEL WITH SCILAB, BUILD WITH VHDL NUMEROUS MODELLING & SIMULATIONS DIRECTLY DESIGN DSP HARDWARE Brought to you by: Copyright(c)

### Computer Music in Undergraduate Digital Signal Processing

Computer Music in Undergraduate Digital Signal Processing Phillip L. De Leon New Mexico State University Klipsch School of Electrical and Computer Engineering Las Cruces, New Mexico 88003-800 pdeleon@nmsu.edu

### speech signal S(n). This involves a transformation of S(n) into another signal or a set of signals

16 3. SPEECH ANALYSIS 3.1 INTRODUCTION TO SPEECH ANALYSIS Many speech processing [22] applications exploits speech production and perception to accomplish speech analysis. By speech analysis we extract

### Laboratory Experiment #1 Introduction to Spectral Analysis

J.B.Francis College of Engineering Mechanical Engineering Department 22-403 Laboratory Experiment #1 Introduction to Spectral Analysis Introduction The quantification of electrical energy can be accomplished

### Tektronix digital oscilloscope, BK Precision Function Generator, coaxial cables, breadboard, the crystal earpiece from your AM radio kit.

Experiment 0: Review I. References The 174 and 275 Lab Manuals Any standard text on error analysis (for example, Introduction to Error Analysis, J. Taylor, University Science Books, 1997) The manual for

### SGN Audio and Speech Processing

Introduction 1 Course goals Introduction 2 SGN 14006 Audio and Speech Processing Lectures, Fall 2014 Anssi Klapuri Tampere University of Technology! Learn basics of audio signal processing Basic operations

### Worksheet for Exploration 31.1: Amplitude, Frequency and Phase Shift

Worksheet for Exploration 31.1: Amplitude, Frequency and Phase Shift We characterize the voltage (or current) in AC circuits in terms of the amplitude, frequency (period) and phase. The sinusoidal voltage

### The Fundamentals of Mixed Signal Testing

The Fundamentals of Mixed Signal Testing Course Information The Fundamentals of Mixed Signal Testing course is designed to provide the foundation of knowledge that is required for testing modern mixed

### Electronics 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

### SGN Audio and Speech Processing

SGN 14006 Audio and Speech Processing Introduction 1 Course goals Introduction 2! Learn basics of audio signal processing Basic operations and their underlying ideas and principles Give basic skills although

### Lab 6: Building a Function Generator

ECE 212 Spring 2010 Circuit Analysis II Names: Lab 6: Building a Function Generator Objectives In this lab exercise you will build a function generator capable of generating square, triangle, and sine

### Advanced Audiovisual Processing Expected Background

Advanced Audiovisual Processing Expected Background As an advanced module, we will not cover introductory topics in lecture. You are expected to already be proficient with all of the following topics,

### Reducing comb filtering on different musical instruments using time delay estimation

Reducing comb filtering on different musical instruments using time delay estimation Alice Clifford and Josh Reiss Queen Mary, University of London alice.clifford@eecs.qmul.ac.uk Abstract Comb filtering

### Fourier Theory & Practice, Part I: Theory (HP Product Note )

Fourier Theory & Practice, Part I: Theory (HP Product Note 54600-4) By: Robert Witte Hewlett-Packard Co. Introduction: This product note provides a brief review of Fourier theory, especially the unique

### L19: Prosodic modification of speech

L19: Prosodic modification of speech Time-domain pitch synchronous overlap add (TD-PSOLA) Linear-prediction PSOLA Frequency-domain PSOLA Sinusoidal models Harmonic + noise models STRAIGHT This lecture

### Post-processing data with Matlab

Post-processing data with Matlab Best Practice TMR7-31/08/2015 - Valentin Chabaud valentin.chabaud@ntnu.no Cleaning data Filtering data Extracting data s frequency content Introduction A trade-off between

### Class #16: Experiment Matlab and Data Analysis

Class #16: Experiment Matlab and Data Analysis Purpose: The objective of this experiment is to add to our Matlab skill set so that data can be easily plotted and analyzed with simple tools. Background:

### AC : FIR FILTERS FOR TECHNOLOGISTS, SCIENTISTS, AND OTHER NON-PH.D.S

AC 29-125: FIR FILTERS FOR TECHNOLOGISTS, SCIENTISTS, AND OTHER NON-PH.D.S William Blanton, East Tennessee State University Dr. Blanton is an associate professor and coordinator of the Biomedical Engineering

### CMPT 468: Delay Effects

CMPT 468: Delay Effects Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University November 8, 2013 1 FIR/Convolution Since the feedforward coefficient s of the FIR filter are

### An Introduction to Time Waveform Analysis

An Introduction to Time Waveform Analysis Timothy A Dunton, Universal Technologies Inc. Abstract In recent years there has been a resurgence in the use of time waveform analysis techniques. Condition monitoring

### Signal 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).

### Digital Signal Processing Laboratory 1: Discrete Time Signals with MATLAB

Digital Signal Processing Laboratory 1: Discrete Time Signals with MATLAB Thursday, 23 September 2010 No PreLab is Required Objective: In this laboratory you will review the basics of MATLAB as a tool

### Graph of the Sine Function

1 of 6 8/6/2004 6.3 GRAPHS OF THE SINE AND COSINE 6.3 GRAPHS OF THE SINE AND COSINE Periodic Functions Graph of the Sine Function Graph of the Cosine Function Graphing Techniques, Amplitude, and Period

### ECE 203 LAB 2 PRACTICAL FILTER DESIGN & IMPLEMENTATION

Version 1. 1 of 7 ECE 03 LAB PRACTICAL FILTER DESIGN & IMPLEMENTATION BEFORE YOU BEGIN PREREQUISITE LABS ECE 01 Labs ECE 0 Advanced MATLAB ECE 03 MATLAB Signals & Systems EXPECTED KNOWLEDGE Understanding

Lecture Topics Doppler CW Radar System, FM-CW Radar System, Moving Target Indication Radar System, and Pulsed Doppler Radar System 1 Remember that: An EM wave is a function of both space and time e.g.

### Time-Frequency Analysis of Shock and Vibration Measurements Using Wavelet Transforms

Cloud Publications International Journal of Advanced Packaging Technology 2014, Volume 2, Issue 1, pp. 60-69, Article ID Tech-231 ISSN 2349 6665, doi 10.23953/cloud.ijapt.15 Case Study Open Access Time-Frequency

### The Discrete Fourier Transform

CHAPTER The Discrete Fourier Transform Fourier analysis is a family of mathematical techniques, all based on decomposing signals into sinusoids. The discrete Fourier transform (DFT) is the family member

### 6.S02 MRI Lab Acquire MR signals. 2.1 Free Induction decay (FID)

6.S02 MRI Lab 1 2. Acquire MR signals Connecting to the scanner Connect to VMware on the Lab Macs. Download and extract the following zip file in the MRI Lab dropbox folder: https://www.dropbox.com/s/ga8ga4a0sxwe62e/mit_download.zip

### University of Pennsylvania Department of Electrical and Systems Engineering Digital Audio Basics

University of Pennsylvania Department of Electrical and Systems Engineering Digital Audio Basics ESE250 Spring 2013 Lab 4: Time and Frequency Representation Friday, February 1, 2013 For Lab Session: Thursday,

### DISCRETE FOURIER TRANSFORM AND FILTER DESIGN

DISCRETE FOURIER TRANSFORM AND FILTER DESIGN N. C. State University CSC557 Multimedia Computing and Networking Fall 2001 Lecture # 03 Spectrum of a Square Wave 2 Results of Some Filters 3 Notation 4 x[n]

### 6.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

### Operational Amplifiers: Part II

1. Introduction Operational Amplifiers: Part II The name "operational amplifier" comes from this amplifier's ability to perform mathematical operations. Three good examples of this are the summing amplifier,

### Introduction to Wavelet Transform. Chapter 7 Instructor: Hossein Pourghassem

Introduction to Wavelet Transform Chapter 7 Instructor: Hossein Pourghassem Introduction Most of the signals in practice, are TIME-DOMAIN signals in their raw format. It means that measured signal is a

### The Fundamentals of FFT-Based Signal Analysis and Measurement Michael Cerna and Audrey F. Harvey

Application ote 041 The Fundamentals of FFT-Based Signal Analysis and Measurement Michael Cerna and Audrey F. Harvey Introduction The Fast Fourier Transform (FFT) and the power spectrum are powerful tools

### Chapter 3 THE DIFFERENTIATOR AND INTEGRATOR Name: Date

AN INTRODUCTION TO THE EXPERIMENTS The following two experiments are designed to demonstrate the design and operation of the op-amp differentiator and integrator at various frequencies. These two experiments

### Image and Video Processing

Image and Video Processing () Image Representation Dr. Miles Hansard miles.hansard@qmul.ac.uk Segmentation 2 Today s agenda Digital image representation Sampling Quantization Sub-sampling Pixel interpolation

### Experiment # 4. Frequency Modulation

ECE 416 Fall 2002 Experiment # 4 Frequency Modulation 1 Purpose In Experiment # 3, a modulator and demodulator for AM were designed and built. In this experiment, another widely used modulation technique

### AC Circuits. "Look for knowledge not in books but in things themselves." W. Gilbert ( )

AC Circuits "Look for knowledge not in books but in things themselves." W. Gilbert (1540-1603) OBJECTIVES To study some circuit elements and a simple AC circuit. THEORY All useful circuits use varying

### Spectral Estimation & Examples of Signal Analysis

Spectral Estimation & Examples of Signal Analysis Examples from research of Kyoung Hoon Lee, Aaron Hastings, Don Gallant, Shashikant More, Weonchan Sung Herrick Graduate Students Estimation: Bias, Variance