Fourier Transform Analysis of Signals and Systems
Ideal Filters Filters separate what is desired from what is not desired In the signals and systems context a filter separates signals in one frequency range from signals in another frequency range An ideal filter passes all signal power in its passband without distortion and completely blocks signal power outside its passband 5/10/04 M. J. Roberts - All Rights Reserved 2
Distortion Distortion is construed in signal analysis to mean changing the shape of a signal Multiplication of a signal by a constant (even a negative one) or shifting it in time do not change its shape No Distortion Distortion 5/10/04 M. J. Roberts - All Rights Reserved 3
( ) ()= h t Aδ t t0 or h[ n]= A n n Distortion Since a system can multiply by a constant or shift in time without distortion, a distortionless system would have an impulse response of the form, [ ] δ 0 The corresponding transfer functions are H( f)= Ae j 2π ft 0 or H( F)= Ae j π Fn 2 0 5/10/04 M. J. Roberts - All Rights Reserved 4
Filter Classifications There are four commonly-used classification of filters, lowpass, highpass, bandpass and bandstop. 5/10/04 M. J. Roberts - All Rights Reserved 5
Filter Classifications 5/10/04 M. J. Roberts - All Rights Reserved 6
Bandwidth Bandwidth generally means a range of frequencies This range could be the range of frequencies a filter passes or the range of frequencies present in a signal Bandwidth is traditionally construed to be range of frequencies in positive frequency space 5/10/04 M. J. Roberts - All Rights Reserved 7
Bandwidth Common Bandwidth Definitions 5/10/04 M. J. Roberts - All Rights Reserved 8
Impulse Responses of Ideal Filters 5/10/04 M. J. Roberts - All Rights Reserved 9
Impulse Responses of Ideal Filters 5/10/04 M. J. Roberts - All Rights Reserved 10
Impulse Response and Causality All the impulse responses of ideal filters contain sinc functions, alone or in combinations, which are infinite in extent Therefore all ideal filter impulse responses begin before time, t = 0 This makes ideal filters non-causal Ideal filters cannot be physically realized, but they can be closely approximated 5/10/04 M. J. Roberts - All Rights Reserved 11
Examples of Impulse Responses and Frequency Responses of Real Causal Filters 5/10/04 M. J. Roberts - All Rights Reserved 12
Examples of Impulse Responses and Frequency Responses of Real Causal Filters 5/10/04 M. J. Roberts - All Rights Reserved 13
Examples of Causal Filter Effects on Signals 5/10/04 M. J. Roberts - All Rights Reserved 14
Examples of Causal Filter Effects on Signals 5/10/04 M. J. Roberts - All Rights Reserved 15
Examples of Causal Filter Effects on Signals 5/10/04 M. J. Roberts - All Rights Reserved 16
Examples of Causal Filter Effects on Signals 5/10/04 M. J. Roberts - All Rights Reserved 17
Two-Dimensional Filtering of Images Causal Lowpass Filtering of Rows in an Image Causal Lowpass Filtering of Columns in an Image 5/10/04 M. J. Roberts - All Rights Reserved 18
Two-Dimensional Filtering of Images Non-Causal Lowpass Filtering of Rows in an Image Non-Causal Lowpass Filtering of Columns in an Image 5/10/04 M. J. Roberts - All Rights Reserved 19
Two-Dimensional Filtering of Images Causal Lowpass Filtering of Rows and Columns in an Image Non-Causal Lowpass Filtering of Rows and Columns in an Image 5/10/04 M. J. Roberts - All Rights Reserved 20
The Power Spectrum 5/10/04 M. J. Roberts - All Rights Reserved 21
Noise Removal A very common use of filters is to remove noise from a signal. If the noise bandwidth is much greater than the signal bandwidth a large improvement in signal fidelity is possible. 5/10/04 M. J. Roberts - All Rights Reserved 22
Practical Passive Filters H( jω)= = Z c Vout ( jω) Vin( jω) Zc( jω) ( jω)+ Z ( jω) = 1 jωrc+ R RC Lowpass Filter 1 5/10/04 M. J. Roberts - All Rights Reserved 23
Practical Passive Filters RLC Bandpass Filter H( f )= V V out in 2πf j RC 2 2πf 1 j2πf j RC LC ( f ) ( f ) = ( ) + + 5/10/04 M. J. Roberts - All Rights Reserved 24
Log-Magnitude Frequency- Response Plots Consider the two (different) transfer functions, H ( f )= 1 30 and H ( f )= j2πf + 1 30 4π f + j62πf 1 2 2 2 When plotted on this scale, these magnitude frequency response plots are indistinguishable. 5/10/04 M. J. Roberts - All Rights Reserved 25
Log-Magnitude Frequency- Response Plots When the magnitude frequency responses are plotted on a logarithmic scale the difference is visible. 5/10/04 M. J. Roberts - All Rights Reserved 26
Bode Diagrams A Bode diagram is a plot of a frequency response in decibels versus frequency on a logarithmic scale. The Bel (B) is the common (base 10) logarithm of a power ratio and a decibel (db) is one-tenth of a Bel. The Bel is named in honor of Alexander Graham Bell. A signal ratio, expressed in decibels, is 20 times the common logarithm of the signal ratio because signal power is proportional to the square of the signal. 5/10/04 M. J. Roberts - All Rights Reserved 27
Bode Diagrams H ( f )= 1 30 and H ( f )= j2πf + 1 30 4π f + j62πf 1 2 2 2 5/10/04 M. J. Roberts - All Rights Reserved 28
Bode Diagrams Continuous-time LTI systems are described by equations of the general form, D k a d t b d k k y()= k k x() t k = 0 dt k = 0 dt Fourier transforming, the transfer function is of the general form, N k bk jω Y jω k H( jω)= ( ) ( ) D X( jω) = = 0 k a jω N k = 0 k k ( ) 5/10/04 M. J. Roberts - All Rights Reserved 29
5/10/04 M. J. Roberts - All Rights Reserved 30 Bode Diagrams A transfer function can be written in the form, H j A j z j z j z j p j p j p N D ω ω ω ω ω ω ω ( )= 1 1 1 1 1 1 1 2 1 2 L L where the z s are the values of jω (not ω) at which the transfer function goes to zero and the p s are the values of jω at which the transfer function goes to infinity. These z s and p s are commonly referred to as the zeros and poles of the system.
Bode Diagrams From the factored form of the transfer function a system can be conceived as the cascade of simple systems, each of which has only one numerator factor or one denominator factor. Since the Bode diagram is logarithmic, multiplied transfer functions add when expressed in db. 5/10/04 M. J. Roberts - All Rights Reserved 31
Bode Diagrams System Bode diagrams are formed by adding the Bode diagrams of the simple systems which are in cascade. Each simple-system diagram is called a component diagram. One Real Pole 1 H( jω)= jω 1 p k 5/10/04 M. J. Roberts - All Rights Reserved 32
Bode Diagrams One real zero H( jω)= 1 jω z k 5/10/04 M. J. Roberts - All Rights Reserved 33
Bode Diagrams Integrator (Pole at zero) H( jω)= 1 jω 5/10/04 M. J. Roberts - All Rights Reserved 34
Bode Diagrams Differentiator (Zero at zero) H( jω)= jω 5/10/04 M. J. Roberts - All Rights Reserved 35
Bode Diagrams Frequency-Independent Gain H( jω)= A (This phase plot is for A > 0. If A < 0, the phase would be a constant π or - π radians.) 5/10/04 M. J. Roberts - All Rights Reserved 36
5/10/04 M. J. Roberts - All Rights Reserved 37 Bode Diagrams Complex Pole Pair H Re j j p j p j p p j p ω ω ω ω ω ( )= = ( ) + ( ) 1 1 1 1 1 2 1 2 1 1 2 2 1 2
Complex Zero Pair Bode Diagrams jω jω H( jω)= 1 1 = 1 z z 1 2 2Re( z1) ( jω) jω 2 + 2 z z 1 1 2 5/10/04 M. J. Roberts - All Rights Reserved 38
Practical Active Filters Operational Amplifiers The ideal operational amplifier has infinite input impedance, zero output impedance, infinite gain and infinite bandwidth. H( f )= V V o i ( f ) Z f f ( f ) = ( ) Z ( f ) i H( f )= Z f ( f)+ Zi( f) Z f i ( ) 5/10/04 M. J. Roberts - All Rights Reserved 39
Practical Active Filters Active Integrator V o ( f )= ( ) 1 Vi f RC 123 j2πf integral of V f i ( ) 5/10/04 M. J. Roberts - All Rights Reserved 40
Practical Active Filters Active RC Lowpass Filter V V R ( f ) = R j2πfcr + 1 0 f f 1 i ( ) 5/10/04 M. J. Roberts - All Rights Reserved 41 s f
Practical Active Filters Lowpass Filter An integrator with feedback is a lowpass filter. y ( t)+ y()= t x() t H( jω)= 1 jω + 1 5/10/04 M. J. Roberts - All Rights Reserved 42
Practical Active Filters Highpass Filter 5/10/04 M. J. Roberts - All Rights Reserved 43
Discrete-Time Filters DT Lowpass Filter H( F)= 1 4 5 1 e j 2π F h[ n]= 4 5 n u[ n] 5/10/04 M. J. Roberts - All Rights Reserved 44
Discrete-Time Filters Comparison of a DT lowpass filter impulse response with an RC passive lowpass filter impulse response 5/10/04 M. J. Roberts - All Rights Reserved 45
Discrete-Time Filters DT Lowpass Filter Frequency Response RC Lowpass Filter Frequency Response 5/10/04 M. J. Roberts - All Rights Reserved 46
Discrete-Time Filters Moving-Average Filter jπnf H( F)= e drcl F, N + 1 ( ) h n δ n δ [ n 1 ]+ δ [ n 2 ]+ L + δ n N N + 1 [ ]= [ ]+ Always Stable [ ] 5/10/04 M. J. Roberts - All Rights Reserved 47
Discrete-Time Filters Ideal DT Lowpass Filter Impulse Response Almost-Ideal DT Lowpass Filter Impulse Response Almost-Ideal DT Lowpass Filter Magnitude Frequency Response 5/10/04 M. J. Roberts - All Rights Reserved 48
Discrete-Time Filters Almost-Ideal DT Lowpass Filter Magnitude Frequency Response in db 5/10/04 M. J. Roberts - All Rights Reserved 49
Advantages of Discrete-Time Filters They are almost insensitive to environmental effects CT filters at low frequencies may require very large components, DT filters do not DT filters are often programmable making them easy to modify DT signals can be stored indefinitely on magnetic media, stored CT signals degrade over time DT filters can handle multiple signals by multiplexing them 5/10/04 M. J. Roberts - All Rights Reserved 50
Communication Systems A naive, absurd communication system 5/10/04 M. J. Roberts - All Rights Reserved 51
Communication Systems A better communication system using electromagnetic waves to carry information 5/10/04 M. J. Roberts - All Rights Reserved 52
Communication Systems Problems Antenna inefficiency at audio frequencies All transmissions from all transmitters are in the same bandwidth, thereby interfering with each other Solution Frequency multiplexing using modulation 5/10/04 M. J. Roberts - All Rights Reserved 53
Communication Systems Double-Sideband Suppressed-Carrier (DSBSC) Modulation ()= () ( ) y t x t cos 2πf t c Modulator 5/10/04 M. J. Roberts - All Rights Reserved 54
Communication Systems Double-Sideband Suppressed-Carrier (DSBSC) Modulation 1 Y( f)= X( f) [ ( f fc)+ ( f + fc) ] 2 δ δ Modulator Frequency multiplexing is using a different carrier frequency,, for each transmitter. f c 5/10/04 M. J. Roberts - All Rights Reserved 55
Communication Systems Double-Sideband Suppressed-Carrier (DSBSC) Modulation Typical received signal by an antenna Synchronous Demodulation 5/10/04 M. J. Roberts - All Rights Reserved 56
Communication Systems Double-Sideband Transmitted-Carrier (DSBTC) Modulation [ ] c ( c ) ()= + () y t 1 mx t A cos 2πf t Modulator m = 1 5/10/04 M. J. Roberts - All Rights Reserved 57
Communication Systems Double-Sideband Transmitted-Carrier (DSBTC) Modulation Modulator Carrier Carrier 5/10/04 M. J. Roberts - All Rights Reserved 58
Communication Systems Double-Sideband Transmitted-Carrier (DSBTC) Modulation Envelope Detector 5/10/04 M. J. Roberts - All Rights Reserved 59
Communication Systems Double-Sideband Transmitted-Carrier (DSBTC) Modulation 5/10/04 M. J. Roberts - All Rights Reserved 60
Communication Systems Single-Sideband Suppressed-Carrier (SSBSC) Modulation Modulator 5/10/04 M. J. Roberts - All Rights Reserved 61
Communication Systems Single-Sideband Suppressed-Carrier (SSBSC) Modulation 5/10/04 M. J. Roberts - All Rights Reserved 62
Communication Systems Quadrature Carrier Modulation Modulator Demodulator 5/10/04 M. J. Roberts - All Rights Reserved 63
Phase and Group Delay Through the time- shifting property of the Fourier transform, a linear phase shift as a function of frequency corresponds to simple delay Most real system transfer functions have a nonlinear phase shift as a function of frequency Non-linear phase shift delays some frequency components more than others This leads to the concepts of phase delay and group delay 5/10/04 M. J. Roberts - All Rights Reserved 64
Phase and Group Delay To illustrate phase and group delay let a system be excited by x t Acos ω t cos ω t X( jω)= ()= ( ) ( ) m Modulation c Carrier π δω ( ω ω )+ δω ω + ω 2 + δω ( + ω ω )+ δω+ ω + ω ( ) ( ) A c m c m c m c m an amplitude-modulated carrier. To keep the analysis simple suppose that the system has a transfer function whose magnitude is the constant, 1, over the frequency range, and whose phase is ωc ωm < ω < ωc + ωm φ( ω) 5/10/04 M. J. Roberts - All Rights Reserved 65
Phase and Group Delay The system response is Y( jω)= Aπ δω ( ωc ωm)+ δω ωc + ωm 2 + δω ( + ω ω )+ δω+ ω + ω ( ) ( ) c m c m e ( ) jφ ω After some considerable algebra, the time-domain response can be written as φω ( c + ωm)+ φω ( c ωm) φω ( + ω ) φω ω y()= t Acos ω c t+ cos ω m t+ 2ω 2ω Carrier c c m c m m Modulation ( ) 5/10/04 M. J. Roberts - All Rights Reserved 66
Phase and Group Delay φω ( c + ωm)+ φω ( c ωm) φω ( + ω ) φω ω y()= t Acos ω c t+ cos ω m t+ 2ω 2ω Carrier c c m c m m Modulation ( ) In this expression it is apparent that the carrier is shifted in time by φω ( + ω )+ φω ω 2ω c m c m c ( ) and the modulation is shifted in time by φω ( + ω ) φω ω 2ω c m c m m ( ) 5/10/04 M. J. Roberts - All Rights Reserved 67
Phase and Group Delay If the phase function is a linear function of frequency, φ( ω)= Kω the two delays are the same, -K. If the phase function is the non-linear function, ω φ( ω)= tan 1 2 ω which is typical of a single-pole lowpass filter, with the carrier delay is ω c 1. 107 ω c = 10ω m c and the modulation delay is 04. ω c 5/10/04 M. J. Roberts - All Rights Reserved 68
Phase and Group Delay On this scale the delays are difficult to see. 5/10/04 M. J. Roberts - All Rights Reserved 69
Phase and Group Delay In this magnified view the difference between carrier delay and modulation delay is visible. The delay of the carrier is phase delay and the delay of the modulation is group delay. 5/10/04 M. J. Roberts - All Rights Reserved 70
Phase and Group Delay The expression for modulation delay, approaches φω ( + ω ) φω ω 2ω c m c m d df as the modulation frequency approaches zero. In that same limit the expression for carrier delay, φω ( c + ωm)+ φω ( c ωm) 2ω c approaches φω ( c) ω 5/10/04 M. J. Roberts - All Rights Reserved 71 m ( φω ( )) c ( ) ω= ω c
Group delay is defined as τω ( )= d ( ω φω ( )) d When the modulation time shift is negative, the group delay is positive. Phase and Group Delay Carrier time shift is proportional to phase shift at any frequency and modulation time shift is proportional to the derivative with respect to frequency of the phase shift. 5/10/04 M. J. Roberts - All Rights Reserved 72
Pulse Amplitude Modulation Pulse amplitude modulation is like DSBSC modulation except that the carrier is a rectangular pulse train, Modulator t t p()= t rect comb w 1 T T s s 5/10/04 M. J. Roberts - All Rights Reserved 73
Pulse Amplitude Modulation The response of the pulse modulator is t t y()= t x() t p()= t x() t rect comb w 1 Ts T and its CTFT is s where f s s s s k = ( )= ( ) ( ) Y f wf sinc wkf X f kf = 1 T s 5/10/04 M. J. Roberts - All Rights Reserved 74
Pulse Amplitude Modulation The CTFT of the response is basically multiple replicas of the CTFT of the excitation with different amplitudes, spaced apart by the pulse repetition rate. 5/10/04 M. J. Roberts - All Rights Reserved 75
Discrete-time modulation is analogous to continuous-time modulation. A modulating signal multiplies a carrier. Let the carrier be c[ n]= cos( πfn) If the modulation is x[n], the response is Discrete-Time Modulation 2 0 [ ]= [ ] ( ) y n x n cos πfn 2 0 5/10/04 M. J. Roberts - All Rights Reserved 76
Discrete-Time Modulation ( )= ( ) ( ) Y F X F C F 1 = X( F F )+ X F+ F 2 ( ) [ ] 0 0 5/10/04 M. J. Roberts - All Rights Reserved 77
Spectral Analysis The heart of a swept-frequency type spectrum analyzer is a multiplier, like the one introduced in DSBSC modulation, plus a lowpass filter. Multiplying by the cosine shifts the spectrum of x(t) by f c and the signal power shifted into the passband of the lowpass filter is measured. Then, as the frequency, f c, is slowly swept over a range of frequencies, the spectrum analyzer measures its signal power versus frequency. 5/10/04 M. J. Roberts - All Rights Reserved 78
Spectral Analysis One benefit of spectral analysis is illustrated below. These two signals are different but exactly how they are different is difficult to see by just looking at them. 5/10/04 M. J. Roberts - All Rights Reserved 79
Spectral Analysis The magnitude spectra of the two signals reveal immediately what the difference is. The second signal contains a sinusoid, or something close to a sinusoid, that causes the two large spikes. 5/10/04 M. J. Roberts - All Rights Reserved 80