ECE 201: Introduction to Signal Analysis

Similar documents
ECE 201: Introduction to Signal Analysis. Dr. B.-P. Paris Dept. Electrical and Comp. Engineering George Mason University

ECE 201: Introduction to Signal Analysis

Synthesis: From Frequency to Time-Domain

Limitations of Sum-of-Sinusoid Signals

Digital Signal Processing Lecture 1 - Introduction

DSP 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

CMPT 368: Lecture 4 Amplitude Modulation (AM) Synthesis

DSP First. Laboratory Exercise #2. Introduction to Complex Exponentials

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

Spectrum. Additive Synthesis. Additive Synthesis Caveat. Music 270a: Modulation

GEORGIA INSTITUTE OF TECHNOLOGY. SCHOOL of ELECTRICAL and COMPUTER ENGINEERING

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

Solution to Chapter 4 Problems

Linear Frequency Modulation (FM) Chirp Signal. Chirp Signal cont. CMPT 468: Lecture 7 Frequency Modulation (FM) Synthesis

Lecture 3 Complex Exponential Signals

Music 270a: Modulation

Sinusoids. Lecture #2 Chapter 2. BME 310 Biomedical Computing - J.Schesser

DSP First. Laboratory Exercise #7. Everyday Sinusoidal Signals

Signal Processing First Lab 02: Introduction to Complex Exponentials Multipath. x(t) = A cos(ωt + φ) = Re{Ae jφ e jωt }

Principles of Communications ECS 332

Here are some of Matlab s complex number operators: conj Complex conjugate abs Magnitude. Angle (or phase) in radians

Lecture 15. Signal Transmission Radio Spectrum. Duplexing Channel Sharing or Multiplexing Modulation. Elec 1200

The Formula for Sinusoidal Signals

Signal Processing First Lab 02: Introduction to Complex Exponentials Direction Finding. x(t) = A cos(ωt + φ) = Re{Ae jφ e jωt }

Lecture 7 Frequency Modulation

CMPT 468: Frequency Modulation (FM) Synthesis

Laboratory Assignment 4. Fourier Sound Synthesis

Frequency Division Multiplexing Spring 2011 Lecture #14. Sinusoids and LTI Systems. Periodic Sequences. x[n] = x[n + N]

1 Introduction and Overview

1 Introduction and Overview

Continuous time and Discrete time Signals and Systems

Chapter 2. Signals and Spectra

HW 6 Due: November 3, 10:39 AM (in class)

Chapter-2 SAMPLING PROCESS

2.1 BASIC CONCEPTS Basic Operations on Signals Time Shifting. Figure 2.2 Time shifting of a signal. Time Reversal.

Music 171: Amplitude Modulation

Angle Modulated Systems

Basic Signals and Systems

DFT: Discrete Fourier Transform & Linear Signal Processing

Lab S-8: Spectrograms: Harmonic Lines & Chirp Aliasing

Frequency-Domain Sharing and Fourier Series

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

Speech, music, images, and video are examples of analog signals. Each of these signals is characterized by its bandwidth, dynamic range, and the

4.1 REPRESENTATION OF FM AND PM SIGNALS An angle-modulated signal generally can be written as

Digital Video and Audio Processing. Winter term 2002/ 2003 Computer-based exercises

Project 2 - Speech Detection with FIR Filters

Biomedical Signals. Signals and Images in Medicine Dr Nabeel Anwar

Amplitude Modulation, II

Outline. Communications Engineering 1

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL and COMPUTER ENGINEERING. ECE 2025 Fall 1999 Lab #7: Frequency Response & Bandpass Filters

Introduction to signals and systems

Lecture 17 z-transforms 2

Problem Set 1 (Solutions are due Mon )

Laboratory Assignment 5 Amplitude Modulation

Lecture 6. Angle Modulation and Demodulation

Objectives. Presentation Outline. Digital Modulation Lecture 03

Complex Sounds. Reading: Yost Ch. 4

y(n)= Aa n u(n)+bu(n) b m sin(2πmt)= b 1 sin(2πt)+b 2 sin(4πt)+b 3 sin(6πt)+ m=1 x(t)= x = 2 ( b b b b

Lab10: FM Spectra and VCO

Fall Music 320A Homework #2 Sinusoids, Complex Sinusoids 145 points Theory and Lab Problems Due Thursday 10/11/2018 before class

Lab 4 Fourier Series and the Gibbs Phenomenon

Communication Channels

6.02 Fall 2012 Lecture #12

Outline. EECS 3213 Fall Sebastian Magierowski York University. Review Passband Modulation. Constellations ASK, FSK, PSK.

Problem Set 8 #4 Solution

Subtractive Synthesis. Describing a Filter. Filters. CMPT 468: Subtractive Synthesis

Signals and Systems EE235. Leo Lam

Islamic University of Gaza. Faculty of Engineering Electrical Engineering Department Spring-2011

Communication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi

Theory of Telecommunications Networks

Continuous vs. Discrete signals. Sampling. Analog to Digital Conversion. CMPT 368: Lecture 4 Fundamentals of Digital Audio, Discrete-Time Signals

Discrete Fourier Transform (DFT)

Chapter 2. Fourier Series & Fourier Transform. Updated:2/11/15

Digital Processing of Continuous-Time Signals

Massachusetts Institute of Technology Dept. of Electrical Engineering and Computer Science Fall Semester, Introduction to EECS 2

Practical Application of Wavelet to Power Quality Analysis. Norman Tse

George Mason University Signals and Systems I Spring 2016

Digital Processing of

Final Exam Solutions June 7, 2004

Communication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi

1.Explain the principle and characteristics of a matched filter. Hence derive the expression for its frequency response function.

EEE - 321: Signals and Systems Lab Assignment 3

Final Exam Solutions June 14, 2006

CSC475 Music Information Retrieval

Sampling and Reconstruction of Analog Signals

George Mason University ECE 201: Introduction to Signal Analysis Spring 2017

ELE636 Communication Systems

Multirate Digital Signal Processing

Signals A Preliminary Discussion EE442 Analog & Digital Communication Systems Lecture 2

Communication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi

Angle Modulation, II. Lecture topics. FM bandwidth and Carson s rule. Spectral analysis of FM. Narrowband FM Modulation. Wideband FM Modulation

Fourier Methods of Spectral Estimation

6.02 Practice Problems: Modulation & Demodulation

Department of Electronic Engineering NED University of Engineering & Technology. LABORATORY WORKBOOK For the Course SIGNALS & SYSTEMS (TC-202)

HW 6 Due: November 9, 4 PM

George Mason University ECE 201: Introduction to Signal Analysis

Data Acquisition Systems. Signal DAQ System The Answer?

Midterm 1. Total. Name of Student on Your Left: Name of Student on Your Right: EE 20N: Structure and Interpretation of Signals and Systems

Signals and Systems Lecture 9 Communication Systems Frequency-Division Multiplexing and Frequency Modulation (FM)

Topic 2. Signal Processing Review. (Some slides are adapted from Bryan Pardo s course slides on Machine Perception of Music)

Transcription:

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 Signals Introduction Example: Square Wave To this point we have focused on sinusoids of identical frequency f x(t) = N A i cos(2πft + φ i ). i=1 Note that the frequency f does not have a subscript i! Showed (in phasor addition rule) that the above sum can always be written as a single sinusoid of frequency f.

Introduction Sum of Sinusoidal Signals Introduction Example: Square Wave We will consider sums of sinusoids of different frequencies: x(t) = N A i cos(2πf i t + φ i ). i=1 Note the subscript on the frequencies f i! This apparently minor difference has dramatic consequences.

Sum of Two Sinusoids Introduction Example: Square Wave x(t) = 4 π 4 cos(2πft π/2) + cos(2π3ft π/2) 3π

Sum of 25 Sinusoids Introduction Example: Square Wave x(t) = 25 n=0 4 cos(2π(2n 1)ft π/2) (2n 1)π

Introduction Example: Square Wave Non-sinusoidal Signals as Sums of Sinusoids If we allow infinitely many sinusoids in the sum, then the result is a square wave signal. The example demonstrates that general, non-sinusoidal signals can be represented as a sum of sinusoids. The sinusods in the summation depend on the general signal to be represented. For the square wave signal we need sinusoids of frequencies (2n 1) f, and 4 amplitudes. (2n 1)π (This is not obvious).

Introduction Example: Square Wave Non-sinusoidal Signals as Sums of Sinusoids The ability to express general signals in terms of sinusoids forms the basis for the frequency domain or spectrum representation. Basic idea: list the ingredients of a signal by specifying amplitudes and phases as well as frequencies of the sinusoids in the sum.

Introduction Example: Square Wave The Spectrum of a Sum of Sinusoids Begin with the sum of sinusoids introduced earlier x(t) = A 0 + i=1 N A i cos(2πf i t + φ i ). i=1 where we have broken out a possible constant term. The term A 0 can be thought of as corresponding to a sinusoid of frequency zero. Using the inverse Euler formula, we can replace the sinusoids by complex exponentials N { Xi x(t) = X 0 + 2 exp(j2πf it) + X i } 2 exp( j2πf it). where X 0 = A 0 and X i = A i e jφ i.

Introduction Example: Square Wave The Spectrum of a Sum of Sinusoids (cont d) Starting with x(t) = X 0 + N i=1 { Xi where X 0 = A 0 and X i = A i e jφ i. 2 exp(j2πf it) + X i } 2 exp( j2πf it). The spectrum representation simply lists the complex amplitudes and frequencies in the summation: X(f ) = {(X 0, 0), ( 1 2 X 1, f 1 ), ( 1 2 X 1, f 1),..., ( 1 2 X N, f N ), ( 1 2 X N, f N)}

Example Sum of Sinusoidal Signals Introduction Example: Square Wave Consider the signal x(t) = 3 + 5 cos(20πt π/2) + 7 cos(50πt + π/4). Using the inverse Euler relationship x(t) = 3+ 5 2 e jπ/2 exp(j2π10t) + 5 2 ejπ/2 exp( j2π10t)+ 7 2 ejπ/4 exp(j2π25t) + 7 2 e jπ/4 exp( j2π25t). Hence, X(f ) = {(3, 0), ( 5 2 e jπ/2, 10), ( 5 2 ejπ/2, 10), ( 7 2 ejπ/4, 25), ( 7 2 e jπ/4, 25)}

Exercise Sum of Sinusoidal Signals Introduction Example: Square Wave Find the spectrum of the signal: x(t) = 6 + 4 cos(10πt + π/3) + 5 cos(20πt π/7).

Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Lecture: From Time-Domain to Frequency-Domain and back

Time-domain and Frequency-domain Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Signals are naturally observed in the time-domain. A signal can be illustrated in the time-domain by plotting it as a function of time. The frequency-domain provides an alternative perspective of the signal based on sinusoids: Starting point: arbitrary signals can be expressed as sums of sinusoids (or equivalently complex exponentials). The frequency-domain representation of a signal indicates which complex exponentials must be combined to produce the signal. Since complex exponentials are fully described by amplitude, phase, and frequency it is sufficient to just specify a list of theses parameters. Actually, we list pairs of complex amplitudes (Ae jφ ) and frequencies f and refer to this list as X(f ).

Time-domain and Frequency-domain Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation It is possible (but not necessarily easy) to find X(f ) from x(t): this is called Fourier or spectrum analysis. Similarly, one can construct x(t) from the spectrum X(f ): this is called Fourier synthesis. Notation: x(t) X(f ). Example (from last time): Time-domain: signal x(t) = 3 + 5 cos(20πt π/2) + 7 cos(50πt + π/4). Frequency Domain: spectrum X(f ) = {(3, 0), ( 5 2 e jπ/2, 10), ( 5 2 ejπ/2, 10), ( 7 2 ejπ/4, 25), ( 7 2 e jπ/4, 25)}

Plotting a Spectrum Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation To illustrate the spectrum of a signal, one typically plots the magnitude versus frequency. Sometimes the phase is plotted versus frequency as well.

Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Why Bother with the Frequency-Domain? In many applications, the frequency contents of a signal is very important. For example, in radio communications signals must be limited to occupy only a set of frequencies allocated by the FCC. Hence, understanding and analyzing the spectrum of a signal is crucial from a regulatory perspective. Often, features of a signal are much easier to understand in the frequency domain. (Example on next slides). We will see later in this class, that the frequency-domain interpretation of signals is very useful in connection with linear, time-invariant systems. Example: A low-pass filter retains low frequency components of the spectrum and removes high-frequency components.

Example: Original signal Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation

Example: Corrupted signal Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation

Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Synthesis: From Frequency to Time-Domain Synthesis is a straightforward process; it is a lot like following a recipe. Ingredients are given by the spectrum X(f ) = {(X 0, 0), (X 1, f 1 ), (X 1, f 1),..., (X N, f N ), (X N, f N)} Each pair indicates one complex exponential component by listing its frequency and complex amplitude. Instructions for combining the ingredients and producing the (time-domain) signal: x(t) = N n= N X n exp(j2πf n t). You should simplify the expression you obtain.

Example Sum of Sinusoidal Signals Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Problem: Find the signal x(t) corresponding to X(f ) = {(3, 0), ( 5 2 e jπ/2, 10), ( 5 2 ejπ/2, 10), ( 7 2 ejπ/4, 25), ( 7 2 e jπ/4, 25)} Solution: x(t) = 3 + 5 2 e jπ/2 e j2π10t + 5 2 ejπ/2 e j2π10t + 7 2 ejπ/4 e j2π25t + 7 2 e jπ/4 e j2π25t Which simplifies to: x(t) = 3 + 5 cos(20πt π/2) + 7 cos(50πt + π/4).

Exercise Sum of Sinusoidal Signals Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Find the signal with the spectrum: X(f ) = {(5, 0), (2e jπ/4, 10), (2e jπ/4, 10), ( 5 2 ejπ/4, 15), ( 5 2 e jπ/4, 15)

Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Analysis: From Time to Frequency-Domain The objective of spectrum or Fourier analysis is to find the spectrum of a time-domain signal. We will restrict ourselves to signals x(t) that are sums of sinusoids N x(t) = A 0 + A i cos(2πf i t + φ i ). i=1 We have already shown that such signals have spectrum: X(f ) = {(X 0, 0), ( 1 2 X 1, f 1 ), ( 1 2 X 1, f 1),..., ( 1 2 X N, f N ), ( 1 2 X N, f N)} where X 0 = A 0 and X i = A i e jφ i. We will investigate some interesting signals that can be written as a sum of sinusoids.

Beat Notes Sum of Sinusoidal Signals Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Consider the signal x(t) = 2 cos(2π5t) cos(2π400t). This signal does not have the form of a sum of sinusoids; hence, we can not determine it s spectrum immediately.

MATLAB Code for Beat Notes Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation % BeatNote plot and play a beat note waveform % Parameters f s = 8192; dur = 2; NP = round ( f s / 5 ) ; f1 = 5; f2 = 400; A = 2; % time axis t t =0:1/ f s : dur ; xx = A cos(2 pi f1 t t ). cos(2 pi f2 t t ) ; plot ( t t ( 1 :NP), xx ( 1 :NP) ) xlabel ( Time ( s ) ) soundsc ( xx, f s ) ;

Beat Notes as a Sum of Sinusoids Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Using the inverse Euler relationships, we can write x(t) = 2 cos(2π5t) cos(2π400t) = 2 1 2 (ej2π5t + e j2π5t ) 1 2 (ej2π400t + e j2π400t ). Multiplying out yields: x(t) = 1 2 (ej2π405t + e j2π405t ) + 1 2 (ej2π395t + e j2π395t ). Applying Euler s relationship, lets us write: x(t) = cos(2π405t) + cos(2π395t).

Spectrum of Beat Notes Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation We were able to rewrite the beat notes as a sum of sinusoids x(t) = cos(2π405t) + cos(2π395t). Note that the frequencies in the sum, 395 Hz and 405 Hz, are the sum and difference of the frequencies in the original product, 5 Hz and 400 Hz. It is now straightforward to determine the spectrum of the beat notes signal: X(f ) = {( 1 2, 405), (1 2, 405), (1 2, 395), (1 2, 395)}

Spectrum of Beat Notes Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation

Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Lecture: Amplitude Modulation and Periodic Signals

Amplitude Modulation Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Amplitude Modulation is used in communication systems. The objective of amplitude modulation is to move the spectrum of a signal m(t) from low frequencies to high frequencies. The message signal m(t) may be a piece of music; its spectrum occupies frequencies below 20 KHz. For transmission by an AM radio station this spectrum must be moved to approximately 1 MHz.

Amplitude Modulation Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Conventional amplitude modulation proceeds in two steps: A constant A is added to m(t) such that A + m(t) > 0 for all t. The sum signal A + m(t) is multiplied by a sinusoid cos(2πf c t), where f c is the radio frequency assigned to the station. Consequently, the transmitted signal has the form: x(t) = (A + m(t)) cos(2πf c t).

Amplitude Modulation Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation We are interested in the spectrum of the AM signal. However, we cannot compute X(f ) for arbitrary message signals m(t). For the special case m(t) = cos(2πf m t) we can find the spectrum. To mimic the radio case, f m would be a frequency in the audible range. As before, we will first need to express the AM signal x(t) as a sum of sinusoids.

Amplitude Modulated Signal Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation For m(t) = cos(2πf m t), the AM signal equals This simplifies to x(t) = (A + cos(2πf m t)) cos(2πf c t). x(t) = A cos(2πf c t) + cos(2πf m t) cos(2πf c t). Note that the second term of the sum is a beat notes signal with frequencies f m and f c. We know that beat notes can be written as a sum of sinusoids with frequencies equal to the sum and difference of f m and f c : x(t) = A cos(2πf c t)+ 1 2 cos(2π(f c+f m )t)+ 1 2 cos(2π(f c f m )t).

Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Spectrum of Amplitude Modulated Signal The AM signal is given by x(t) = A cos(2πf c t)+ 1 2 cos(2π(f c+f m )t)+ 1 2 cos(2π(f c f m )t). Thus, its spectrum is X(f ) = { ( A 2, f c), ( A 2, f c), ( 1 4, f c + f m ), ( 1 4, f c f m ), ( 1 4, f c f m ), ( 1 2, f c + f m )}

Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Spectrum of Amplitude Modulated Signal For A = 2, fm = 50, and fc = 400, the spectrum of the AM signal is plotted below.

Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Spectrum of Amplitude Modulated Signal It is interesting to compare the spectrum of the signal before modulation and after multiplication with cos(2πf c t). The signal s(t) = A + m(t) has spectrum S(f ) = {(A, 0), ( 1 2, 50), (1 2, 50)}. The modulated signal x(t) has spectrum X(f ) = { ( A 2, 400), ( A 2, 400), ( 1 4, 450), ( 1 4, 450), ( 1 4, 350), ( 1 2, 350)} Both are plotted on the next page.

Spectrum before and after AM Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation

Spectrum before and after AM Why Bother with the Frequency-Domain? Synthesis: From Frequency to Time-Domain Analysis: From Time to Frequency-Domain Amplitude Modulation Comparison of the two spectra shows that amplitude module indeed moves a spectrum from low frequencies to high frequencies. Note that the shape of the spectrum is precisely preserved. Amplitude modulation can be described concisely by stating: Half of the original spectrum is shifted by f c to the right, and the other half is shifted by f c to the left. Question: How can you get the original signal back so that you can listen to it. This is called demodulation.

What are? Harmonic Frequencies Fourier Series A signal x(t) is called periodic if there is a constant T 0 such that x(t) = x(t + T 0 ) for all t. In other words, a periodic signal repeats itself every T 0 seconds. The interval T 0 is called the fundamental period of the signal. The inverse of T 0 is the fundamental frequency of the signal. Example: A sinusoidal signal of frequency f is periodic with period T 0 = 1/f.

Harmonic Frequencies Harmonic Frequencies Fourier Series Consider a sum of sinusoids: x(t) = A 0 + N A i cos(2πf i t + φ i ). i=1 A special case arises when we constrain all frequencies f i to be integer multiples of some frequency f 0 : f i = i f 0. The frequencies f i are then called harmonic frequencies of f 0. We will show that sums of sinusoids with frequencies that are harmonics are periodic.

Harmonic Signals are Periodic Harmonic Frequencies Fourier Series To establish periodicity, we must show that there is T 0 such x(t) = x(t + T 0 ). Begin with x(t + T 0 ) = A 0 + N i=1 A i cos(2πf i (t + T 0 ) + φ i ) = A 0 + N i=1 A i cos(2πf i t + 2πf i T 0 + φ i ) Now, let f 0 = 1/T 0 and use the fact that frequencies are harmonics: f i = i f 0.

Harmonic Signals are Periodic Harmonic Frequencies Fourier Series Then, f i T 0 = i f 0 T 0 = i and hence x(t + T 0 ) = A 0 + N i=1 A i cos(2πf i t + 2πf i T 0 + φ i ) = A 0 + N i=1 A i cos(2πf i t + 2πi + φ i ) We can drop the 2πi terms and conclude that x(t + T 0 ) = x(t). Conclusion: A signal of the form x(t) = A 0 + N A i cos(2πi f 0 t + φ i ) i=1 is periodic with period T 0 = 1/f 0.

Harmonic Frequencies Fourier Series Finding the Fundamental Frequency Often one is given a set of frequencies f 1, f 2,..., f N and is required to find the fundamental frequency f 0. Specifically, this means one must find a frequency f 0 and integers n 1, n 2,..., n N such that all of the following equations are met: f 1 = n 1 f 0 f 2 = n 2 f 0. f N = n N f 0 Note that there isn t always a solution to the above problem. However, if all frequencies are integers a solution exists. Even if all frequencies are rational a solution exists.

Example Sum of Sinusoidal Signals Harmonic Frequencies Fourier Series Find the fundamental frequency for the set of frequencies f 1 = 12, f 2 = 27, f 3 = 51. Set up the equations: 12 = n 1 f 0 27 = n 2 f 0 51 = n 3 f 0 Try the solution n 1 = 1; this would imply f 0 = 12. This cannot satisfy the other two equations. Try the solution n 1 = 2; this would imply f 0 = 6. This cannot satisfy the other two equations. Try the solution n 1 = 3; this would imply f 0 = 4. This cannot satisfy the other two equations. Try the solution n 1 = 4; this would imply f 0 = 3. This can satisfy the other two equations with n 2 = 9 and n 3 = 17.

Example Sum of Sinusoidal Signals Harmonic Frequencies Fourier Series Note that the three sinusoids complete a cycle at the same time at T 0 = 1/f 0 = 1/3s.

Exercise Sum of Sinusoidal Signals Harmonic Frequencies Fourier Series Find the fundamental frequency for the set of frequencies f 1 = 2, f 2 = 3.5, f 3 = 5.

Fourier Series Sum of Sinusoidal Signals Harmonic Frequencies Fourier Series We have shown that a sum of sinusoids with harmonic frequencies is a periodic signal. One can turn this statement around and arrive at a very important result: Any periodic signal can be expressed as a sum of sinusoids with harmonic frequencies. The resulting sum is called the Fourier Series of the signal. Put differently, a periodic signal can always be written in the form x(t) = A 0 + N i=1 A i cos(2πif 0 t + φ i ) = X 0 + N i=1 X ie j2πif0t + Xi e j2πif 0t with X 0 = A 0 and X i = A i 2 ejφ i.

Fourier Series Sum of Sinusoidal Signals Harmonic Frequencies Fourier Series For a periodic signal the complex amplitudes X i can be computed using a (relatively) simple formula. Specifically, for a periodic signal x(t) with fundamental period T 0 the complex amplitudes X i are given by: X i = 1 T0 x(t) e j2πit/t 0 dt. T 0 0 Note that the integral above can be evaluated over any interval of length T 0.

Example: Square Wave Harmonic Frequencies Fourier Series A square wave signal can be written as x(t) = x(t) = n=0 { 1 0 t < T 0 2 1 T 0 2 t < T 0 4 cos(2π(2n 1)ft π/2) (2n 1)π

Harmonic Frequencies Fourier Series 25-Term Approximation to Square Wave x(t) = 25 n=0 4 cos(2π(2n 1)ft π/2) (2n 1)π

Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals Lecture:

Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals Limitations of Sum-of-Sinusoid Signals So far, we have considered only signals that can be written as a sum of sinusoids. N x(t) = A 0 + A i cos(2πf i t + φ i ). i=1 For such signals, we are able to compute the spectrum. Note, that signals of this form are assumed to last forever, i.e., for < t <, and their spectrum never changes. While such signals are important and useful conceptually, they don t describe real-world signals accurately. Real-world signals are of finite duration, their spectrum changes over time.

Musical Notation Sum of Sinusoidal Signals Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals Musical notation ( sheet music ) provides a way to represent real-world signals: a piece of music. As you know, sheet music places notes on a scale to reflect the frequency of the tone to be played, uses differently shaped note symbols to indicate the duration of each tone, provides the order in which notes are to be played. In summary, musical notation captures how the spectrum of the music-signal changes over time. We cannot write signals whose spectrum changes with time as a sum of sinusoids. A static spectrum is insufficient to describe such signals. Alternative: time-frequency spectrum

Example: Musical Scale Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals Note C D E F G A B C Frequency (Hz) 262 294 330 349 392 440 494 523 Table: Musical Notes and their Frequencies

Example: Musical Scale Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals If we play each of the notes for 250 ms, then the resulting signal can be summarized in the time-frequency spectrum below.

MATLAB Spectrogram Function Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals MATLAB has a function spectrogram that can be used to compute the time-frequency spectrum for a given signal. The resulting plots are similar to the one for the musical scale on the previous slide. Typically, you invoke this function as spectrogram( xx, 256, 128, 256, fs), where xx is the signal to be analyzed and fs is the sampling frequency. The spectrogram for the musical scale is shown on the next slide.

Spectrogram: Musical Scale Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals The color indicates the magnitude of the spectrum at a given time and frequency.

Chirp Signals Sum of Sinusoidal Signals Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals Objective: construct a signal such that its frequency increases with time. Starting Point: A sinusoidal signal has the form: x(t) = A cos(2πf 0 t + φ). We can consider the argument of the cos as a time-varying phase function Ψ(t) = 2πf 0 t + φ. Question: What happens when we allow more general functions for Ψ(t)? For example, let Ψ(t) = 700πt 2 + 440πt + φ.

Spectrogram: cos(ψ(t)) Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals Question: How is he time-frequency spectrum related to Ψ(t)?

Instantaneous Frequency Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals For a regular sinusoid, Ψ(t) = 2πf 0 t + φ and the frequency equals f 0. This suggests as a possible relationship between Ψ(t) and f 0 f 0 = 1 d dt Ψ(t). 2π If the above derivative is not a constant, it is called the instantaneous frequency of the signal, f i (t). Example: For Ψ(t) = 700πt 2 + 440πt + φ we find f i (t) = 1 d 2π dt (700πt2 + 440πt + φ) = 700t + 220. This describes precisely the red line in the spectrogram on the previous slide.

Constructing a Linear Chirp Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals Objective: Construct a signal such that its frequency is initially f 1 and increases linear to f 2 after T seconds. Solution: The above suggests that f i (t) = f 2 f 1 T t + f 1. Consequently, the phase function Ψ(t) must be Ψ(t) = 2π f 2 f 1 2T t2 + 2πf 1 t + φ Note that φ has no influence on the spectrum; it is usually set to 0.

Constructing a Linear Chirp Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals Example: Construct a linear chirp such that the frequency decreases from 1000 Hz to 200 Hz in 2 seconds. The desired signal must be x(t) = cos( 2π200t 2 + 2π1000t).

Exercise Sum of Sinusoidal Signals Limitations of Sum-of-Sinusoid Signals Musical Notation Chirp Signals Construct a linear chirp such that the frequency increases from 50 Hz to 200 Hz in 3 seconds. Sketch the time-frequency spectrum of the following signal x(t) = cos(2π500t + 100 cos(2π2t))