MATLAB Signal Processing Toolbox. Greg Reese, Ph.D Research Computing Support Group Academic Technology Services Miami University
|
|
- Lenard Clarke
- 6 years ago
- Views:
Transcription
1 MATLAB Signal Processing Toolbox Greg Reese, Ph.D Research Computing Support Group Academic Technology Services Miami University October 2013
2 MATLAB Signal Processing Toolbox 2013 Greg Reese. All rights reserved 2
3 Toolbox Toolbox Collection of code devoted to solving problems in one field of research Can be purchased from MATLAB Can be purchased from third parties Can be obtained for free from third parties 3
4 Toolbox MATLAB Signal Processing Toolbox Code for solving problems in signal processing (!) Sold by MATLAB Part of both Miami s student and faculty license 4
5 Overview MATLAB divides Signal Processing Toolbox as follows Waveforms Pulses, modulated signals, peak-to-peak and RMS amplitude, rise time/fall time, overshoot/undershoot Convolution and Correlation Linear and circular convolution, autocorrelation, autocovariance, cross-correlation, cross-covariance Transforms Fourier transform, chirp z-transform, DCT, Hilbert transform, cepstrum, Walsh-Hadamard transform 5
6 Overview Analog and Digital Filters Analog filter design, frequency transformations, FIR and IIR filters, filter analysis, filter structures Spectral Analysis Nonparametric and parametric spectral estimation, high resolution spectral estimation Parametric Modeling and Linear Prediction Autoregressive (AR) models, linear predictive coding (LPC), Levinson-Durbin recursion Multirate Signal Processing Downsampling, upsampling, resampling, anti-aliasing filter, interpolation, decimation 6
7 Overview Will look very briefly at Analog and Digital Filters Spectral Analysis Parametric Modeling and Linear Prediction Multirate Signal Processing Will look in more depth at Waveforms Convolution and Correlation 7
8 Analog and Digital Filters Toolbox especially good for those serious about their filter design! Analog filters Standard filters Bessel, Butterworth, Chebyshev, Elliptic Filter transforms Low pass to: bandpass, bandstop, or highpass Change cutoff frequency of lowpass Analog to digital filter conversion Bilinear transformation 8
9 Analog and Digital Filters Digital Filter Design with functions Standard filters Butterworth, Chebyshev, Elliptic FIR and IIR design Low pass to: bandpass, bandstop, or highpass Change cutoff frequency of lowpass Objects for specification of filters Arbitrary, lowpass, highpass, bandpass, Hilbert 9
10 Analog and Digital Filters Digital Filter Design interactively (GUI) Filterbuilder specify desired characteristics first, then choose filter type Butterworth, Chebyshev, Elliptic FDATool (Filter Design and Analysis Tool) Quickly design digital FIR or IIR filters Quickly analyze filters, e.g., magnitude/phase response, pole-zero plots 10
11 Analog and Digital Filters SPTool composite of four tools 1. Signal Browser analyze signals 2. FDATool 3. FVTool (Filter Visualization Tool) analyze filter characteristics 4. Spectrum Viewer spectral analysis 11
12 Analog and Digital Filters Digital Filter Analysis Magnitude and phase response, impulse response, group delay, pole-zero plot Digital Filter Implementation Filtering, direct form, lattice, biquad, statespace structures 12
13 Spectral Analysis Nonparametric Methods Periodogram, Welch's overlapped segment averaging, multitaper, cross-spectrum, coherence, spectrogram Parametric Methods Yule-Walker, Burg, covariance, modified covariance Subspace Methods Multiple signal classification (MUSIC), eigenvectorestimator, pseudospectrum Windows Hamming, Blackman, Bartlett, Chebyshev, Taylor, Kaiser 13
14 Parametric Modeling and Linear Prediction Parametric Modeling AR, ARMA, frequency response modeling Linear Predictive Waveforms Linear predictive coefficients (LPC), line spectral frequencies (LSF), reflection coefficients (RC), Levinson-Durbin recursion 14
15 Multirate Signal Processing Multirate signal processing Downsampling, upsampling, resampling, anti-aliasing filter, interpolation, decimation 15
16 Waveforms Waveforms part of toolbox lets you create many commonly used signals, which you can use to study models programmed in MATLAB Uses of waveforms Testing E.g., have simple waveform and can analytically determine model s output. Use toolbox to create that waveform, run it through MATLAB model, confirm result 16
17 Uses of waveforms Waveforms Simulation Most of time can t get analytical output Make waveform of known characteristics and study model s response Modeling of real signals Create waveform that looks like the real signal 17
18 Time vectors Waveforms Digital signals usually sampled from analog at fixed intervals Δt. Want time axis with N points 0 1Δt 2Δt (N-2)Δt (N-1)Δt 18
19 For tstart: starting time tend: ending time N: number of points Waveforms deltat: sampling interval If have starting time, number of points, interval, (tstart, N, deltat): >> deltat = 0.1; >> N = 10; >> t0 = 5; >> t = t0 + deltat * (0:N-1) t =
20 Waveforms If have starting time, ending time, interval (tstart, tend, deltat) >> tstart = 5; >> tend = 5.9; >> deltat = 0.1; >> t = tstart:deltat:tend t = If have starting time, ending time, number of points (tstart, tend, N) >> tstart = 5; >> tend = 5.9; >> N = 10; >> t = linspace( tstart, tend, N ) t =
21 Waveforms In multichannel processing, a number of signals are gathered at the same time Will assume all sampled at same time and same rate Signal processing toolbox, and MATLAB in general, treats each column of a matrix (2D array) as an independent column vector 21
22 Waveforms Example >> M = [ 1:3; 4:6; 7:9; 10:12 ] M = >> mean( M ) ans = Result is mean of each column 22
23 TIP Waveforms Make time vector be a column vector Any vectors created from time vector will also be column vectors and so can be processed more easily >> t = t >> y = abs( t - 3 ) t = 1 y = Column vector begeteth column vector
24 Waveforms TIP repmat (replicate matrix) is general purpose function to make large matrix by replicating small one Trick - quick way to replicate column vector, i.e., to make an m x n matrix T out of an m x 1 column vector v, is T = v(:,ones(1,n)) >> v = (1:5) >> v(:,ones(1,3)) v = 1 ans =
25 TIP Waveforms Trick - quick way to replicate row vector, i.e., to make an m x n matrix T out of an 1 x n row vector v, is T = v(ones(m,1),:) >> v = 1:3 >> T = v( ones(6,1), : ) v = ans =
26 Waveforms TIP Can use preceding two tips to make multichannel signal, e.g., Simulate the multichannel signal sin(2πt), sin(2πt/2), sin(2πt/3), sin(2πt/4) sampled for one second at one-tenth second per sample 26
27 Waveforms TIP >> t = (0:0.1:1)' t =
28 Waveforms TIP >> T = t(:,ones(1,4)) T =
29 TIP >> M = T./ C M = Waveforms
30 Waveforms TIP >> signal = sin( 2*pi*M ) signal =
31 Impulse Waveforms Use to compute impulse response of linear, time-invariant system >> t = (0:0.1:1.9)'; >> impulse = zeros( size(t) ); >> impulse( 1 ) = 1; Impulse
32 Step Waveforms Use to model switch turning on >> t = (-1:0.1:0.9)'; >> step = [ zeros(10,1); ones(10,1) ]; 1 Step
33 Ramp Waveforms Use to model something gradually turning on >> t = (0:0.1:1.9)'; >> ramp = t; Ramp
34 Autocorrelation Autocorrelation the detection of a delayed version of a signal In temporal signal, delay often called lag In spatial signal, delay often called translation or offset Delayed signal may also be scaled Can also think of autocorrelation as similarity of a signal to itself as a function of lag 34
35 Autocorrelation Autocorrelation and cross-correlation Common in both signal processing and statistics Definitions and uses are different When looking for help on these topics, make sure you re looking at a signal-processing source 35
36 Autocorrelation Examples of autocorrelation of digital signals Radar determine distance to object Sonar determine distance to object Music Determine tempo Detect and estimate pitch Astronomy Find rotation frequency of pulsars 36
37 Autocorrelation Examples of spatial autocorrelation Optical Character Recognition (OCR) reading text from images of writing/printing X-ray diffraction helps recover the "Fourier phase information" on atom positions Statistics helps estimate mean value uncertainties when sampling a heterogeneous population Astrophysics used to study and characterize the spatial distribution of galaxies 37
38 Autocorrelation Examples of optical autocorrelation Measurement of optical spectra and of very short-duration light pulses produced by lasers Analysis of dynamic light scattering data to determine particle size distributions The small-angle X-ray scattering intensity of some systems related to the spatial autocorrelation function of the electron density In optics, normalized autocorrelations and cross-correlations give the degree of coherence of an electromagnetic field 38
39 Typical use Blip sent to object Autocorrelation Small blip reflected from object back to sender Use autocorrelation to detect small blip at some lag Know velocity of blip in medium so total distance blip traveled is distance = velocity * lag Distance is round trip, so object distance/2 away 39
40 Autocorrelation Autocorrelation Multiply and sum. Result is autocorrelation at that point Slide over one, multiply and sum x x x x = = 0 40
41 Autocorrelation Autocorrelation Repeat, sliding in both directions until have covered all positions What happens when go past end? x x x ? ? 41
42 Autocorrelation When go past end, have two options Zero padding put zeros on both end of both signals. Can either imagine they are there or actually extend arrays in memory and put in zeros x x x Will discuss second option later = 0 42
43 Autocorrelation Suppose the real discrete-time signal x[n] has L points and the real discrete-time signal y[n] has N points, with L N. The autocorrelation of x and y is L m 1 R xy m = x n + m y[n] n=0 for m = -(N-1), -(N-2),, -1, 0, 1, 2,, L-1 43
44 Autocorrelation Aside For p 0, x[n-p] is x[n] shifted to the right by p For p 0, x[n+p] is x[n] shifted to the left by p Example Unit impulse x n = 1 for n = 0 0 for n x[n] x[n+5] x[n-2] n 44
45 TRY IT Autocorrelation At time n=0 a transmitter sends out a pulse of amplitude ten and duration 3. At time time n=5 it gets back the reflected pulse with the same duration but one tenth the amplitude. What is the autocorrelation? 45
46 Autocorrelation TRY IT Sent Received R xy (0) =? R xy (0)= R xy (1) =? R xy (1)= R xy (2) =? R xy (2)=0 46
47 TRY IT R xy (-1)=0 Autocorrelation R xy (-1) =? R xy (-2)= R xy (-2) =? R xy (-3)= R xy (-3) =? R xy (-4)= R xy (-4) =? R xy (-5)= R xy (-5) =?
48 Autocorrelation TRY IT R xy (-6)=20 R xy (-6) =? R xy (-7)=10 R xy (-7) =? R xy (-8)=0 R xy (-8) =? R xy (-9)=0 R xy (-9) =?
49 Autocorrelation TRY IT Put it together >> m = 2:-1:-9; >> R = [ ]; >> [~,maxindex] = max( R ) maxindex = 8 >> m(maxindex) ans = -5 % max when shifted right by
50 Autocorrelation Autocorrelation(m) TRY IT >> plot( -m, R, o ) Note shape is a triangle, not a rectangle, which is shape of pulse. Autocorrelation detects signal of given shape it does not replicate signal Lag m 50
51 Autocorrelation MATLAB considers what we re doing to be cross-correlation Concept is same as what described here for autocorrelation If one array shorter than another, MATLAB pads shorter one with zeros until both same length 51
52 Autocorrelation To compute cross-correlation of vectors x and y in MATLAB, use where c = xcorr( x, y ) c is a vector with 2N-1 elements N is length of longer of x and y If m is lag as previously defined, c(k) is autocorrelation for lag m = k - N 52
53 Autocorrelation TRY IT Let s do previous graphical autocorrelation with MATLAB >> x = [ ]; >> y = [ ]; >> c = xcorr( x, y ); >> [~,maxcix] = max( c ) maxcix = 5 >> m = maxcix - length( y ) m = -5 % Move 5 to right from element 1 53
54 Autocorrelation TRY IT Example of finding signal buried in noise 1. Make a sine wave with a period of 20 and amplitude of 100 >> wave = 100 * sin( 2*pi*(0:19)/20 ); 2. Reset random number generator (so we all get the same random numbers) >> rng default 3. Make 500 points of noise with randn and variance 75% of wave amplitude noisywave = 75 * randn( 1, 500 ); 54
55 Autocorrelation TRY IT 4. Pick a random spot to place the wave, ensuring that the whole wave fits in ix = randi( [ 1, 481 ] ); 5. Add the wave to the noise noisywave(ix:ix+19)=noisywave(ix:ix+19)+wave; Plot wave in noise. Is wave visible? plot( noisywave ) Compute autocorrelation >> c = xcorr( wave, noisywave );
56 Autocorrelation TRY IT 7. Find max of autocorrelation and calculate lag from that >> [~,maxix] = max( c ) maxix = 269 >> m = maxix 500 m = -231 % shift 231 to right 8. Show random spot where wave added to noise. Match? >> ix ix = 231 Very close match! m should be
57 Autocorrelation TRY IT 9. For grins, plot autocorrelation >> lags = -499:499; >> plot( lags, c ) Why is almost all of right size zero? 8 x Right side 0 corresponds to -2 shifting left and once -4 shift wave more than -6 20, rest of wave -8 is zeros
58 Correlation Questions? 58
59 convolution Uses Convolution Polynomial multiplication LTI response Joint PDF Linear and circular Explain, show when equivalent (padding), good for computing convolution with fft. Do example with tic,toc, time-domain convolution vs. fft,ifft, see cconv 59
60 Convolution Applications of convolution Acoustics reverberation is the convolution of the original sound with echoes from objects surrounding the sound source Computational fluid dynamics the large eddy simulation (LES) turbulence model uses convolution to lower the range of length scales necessary in computation and thereby reducing computational cost Probability probability distribution of the sum of two independent random variables is the convolution of their individual distributions 60
61 Convolution Applications of convolution Spectroscopy line broadening can be due to the Doppler effect and/or collision broadening. The effect due to both is the convolution of the two effects Electronic music imposition of a rhythmic structure on a sound done by convolution Image processing blurring, sharpening, and edge enhancement done by convolution Numerical computation can multiply polynomials quickly with convolution 61
62 Convolution Convolution finds many applications because it is central to linear, timeinvariant systems and many things can be modeled by such systems 62
63 Convolution Linear a linear system obeys two principles Principle of superposition the output to a sum of inputs is equal to the sum of the outputs to the individual inputs Scaling the output to the product of an input and a constant is the product of the constant and the output to the input alone In other words, for a linear system L, a L{ x(t) } + bl{ x(t) } = L{ ax(t) + bx(t) } 63
64 Convolution Suppose you put some input into a system and get some output. If you put in the same input at a later time, if the system is time invariant, the output will be the same as the original output except it will occur at that later time In other words, for a time-invariant system S, If y(t) = S{ x(t) }, then y(t-t 0 ) = S{ x(t-t 0 ) } Time-invariance and linearity are independent. A linear system can be time-invariant or not. A timeinvariant system can be linear or not. 64
65 Example Convolution Change machine at a laundromat. Put in dollar bills, press button, get out quarters Linear? 1. Put $1 in, press button, get 4 quarters out 2. Put $2 in, press button, get 8 quarters out 3. (output from $1) + (output from $2) = 12 quarters 4. Put $3, press button, get 12 quarters out 5. Two outputs equal, so system linear 65
66 Example Time invariant? Convolution 1. Put $1 in, press button, get 4 quarters out 2. Put $2 in, press button, get 8 quarters out An hour later 1. Put $1 in, press button, get 4 quarters out 2. Put $2 in, press button, get 8 quarters out Outputs identical except for same delay as input, so system is time invariant 66
67 Convolution (discrete) impulse: δ n = 1 for n = 0 0 for n > 0 impulse response response h[n] of a system S when the input is an impulse, i.e., h[n] = S { δ n } 67
68 Major fact Convolution The output of a linear, time-invariant (LTI) system to any input is the convolution of that input with the system s impulse response Other words: The impulse response of an LTI system completely characterizes that system The impulse response of an LTI system specifies that system 68
69 1 2 3 Convolution Graphical view of convolving two signals Pick one signal Flip it 180 around left edge Position right element of flipped signal over left element of unflipped signal
70 Graphical view Convolution Multiply corresponding elements and sum x
71 Graphical view Convolution Slide right, multiply, sum x x
72 Graphical view Convolution Repeat until fall off right side x x x
73 Graphical view Convolution x x x
74 Graphical view Convolution x x x
75 Graphical view Convolution x x x (-2)
76 Graphical view Convolution x x x (-2) + 3 (-2)
77 Graphical view Convolution x x x (-2) + 2 (-2)
78 Graphical view Convolution x x x (-2)
79 Graphical view Convolution x x x
80 Graphical view Convolution As with correlation, ignore elements that have fallen off or pad bottom array with zeros x x
81 Graphical view Convolution x
82 So convolution of with gives Convolution
83 Convolution Note that when we convolved with the shorter signal was not completely on top of the longer one for the first two and last two elements of the longer. For these cases, the multiply-and-add convolution computation was missing either one or two terms, so those four calculations are not valid and should be ignored. Bad values at the left and right of a convolution are known as edge effects. 83
84 Convolution In general, if the shorter of two signals in a convolution has M elements, you should ignore the first and last M-1 elements in the result 84
85 Convolution Mathematical definition Suppose we have two signals x[n] has M points and y[n] has N points, M N. The convolution of x[n] and y[n] is w n = n k=0 n k=n M+1 n x n k y k for n = 0,1,2, M 1 x n k y[k] k=n M+1 n = 0, 1, 2, N+M-2 for n = M, M + 1, N 1 x n k y k for n = N, N + 1, N + M 2 85
86 Convolution Compute a convolution in MATLAB with w = conv( u, v ) where u and v are vectors. The output vector w has length length(u)+length(v)-1 86
87 Convolution TRY IT We figured out graphically that convolved with gave Try it in MATLAB >> u = [ ]; >> v = [ ]; >> w = conv( u, v ) w =
88 Convolution Typically convolution involves a data signal and another signal. The second signal, called a kernel, is Also called a filter Usually much shorter than the data signal Designed to produce a desired effect on the data 88
89 Convolution TRY IT To introduce a time lag of m units into the data signal, i.e., to shift it to the right by m, use a kernel of m zeros followed by a one. Example Introduce a lag 2 of 5 two -2 into -2 the 0 signal 0 1 >> kernel = [ ]; >> w = conv( kernel, v ) w =
90 Convolution TRY IT Approximate the derivative by replacing a point with the difference between itself and the previous point Example Approximate the derivative of >> u = [ 1-1 ]; >> w = conv( u, v ) w =
91 Convolution TRY IT Reduce noise in a signal by replacing a point with a weighted average of itself and the previous two points Example >> kernel = (1/9)*[ ]; >> w = conv( u, v ); w*9 =
92 Convolution Questions? 92
93 Signal Processing Toolbox Questions? 93
94 The End 94
DISCRETE-TIME SIGNAL PROCESSING
THIRD EDITION DISCRETE-TIME SIGNAL PROCESSING ALAN V. OPPENHEIM MASSACHUSETTS INSTITUTE OF TECHNOLOGY RONALD W. SCHÄFER HEWLETT-PACKARD LABORATORIES Upper Saddle River Boston Columbus San Francisco New
More informationComputer Engineering 4TL4: Digital Signal Processing
Computer Engineering 4TL4: Digital Signal Processing Day Class Instructor: Dr. I. C. BRUCE Duration of Examination: 3 Hours McMaster University Final Examination December, 2003 This examination paper includes
More informationAdvanced Digital Signal Processing -Introduction
Advanced Digital Signal Processing -Introduction LECTURE-2 1 AP9211- ADVANCED DIGITAL SIGNAL PROCESSING UNIT I DISCRETE RANDOM SIGNAL PROCESSING Discrete Random Processes- Ensemble Averages, Stationary
More informationDFT & Fast Fourier Transform PART-A. 7. Calculate the number of multiplications needed in the calculation of DFT and FFT with 64 point sequence.
SHRI ANGALAMMAN COLLEGE OF ENGINEERING & TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR,TRICHY-621105. DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING UNIT I DFT & Fast Fourier
More informationFilter Analysis and Design
Filter Analysis and Design Butterworth Filters Butterworth filters have a transfer function whose squared magnitude has the form H a ( jω ) 2 = 1 ( ) 2n. 1+ ω / ω c * M. J. Roberts - All Rights Reserved
More informationE : Lecture 1 Introduction
E85.2607: Lecture 1 Introduction 1 Administrivia 2 DSP review 3 Fun with Matlab E85.2607: Lecture 1 Introduction 2010-01-21 1 / 24 Course overview Advanced Digital Signal Theory Design, analysis, and implementation
More informationIT DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A
DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING IT6502 - DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A 1. What is a continuous and discrete time signal? Continuous
More informationAutomatic Autocorrelation and Spectral Analysis
Piet M.T. Broersen Automatic Autocorrelation and Spectral Analysis With 104 Figures Sprin ger 1 Introduction 1 1.1 Time Series Problems 1 2 Basic Concepts 11 2.1 Random Variables 11 2.2 Normal Distribution
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Filtering in the Frequency Domain http://www.ee.unlv.edu/~b1morris/ecg782/ 2 Outline Background
More informationCentre for Mathematical Sciences HT 2017 Mathematical Statistics
Lund University Stationary stochastic processes Centre for Mathematical Sciences HT 2017 Mathematical Statistics Computer exercise 3 in Stationary stochastic processes, HT 17. The purpose of this exercise
More informationFilter structures ELEC-E5410
Filter structures ELEC-E5410 Contents FIR filter basics Ideal impulse responses Polyphase decomposition Fractional delay by polyphase structure Nyquist filters Half-band filters Gibbs phenomenon Discrete-time
More informationLecture 3: Linear Filters
Lecture 3: Linear Filters Professor Fei Fei Li Stanford Vision Lab 1 What we will learn today? Images as functions Linear systems (filters) Convolution and correlation Discrete Fourier Transform (DFT)
More informationDHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A
DHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A Classification of systems : Continuous and Discrete
More informationEE538 Final Exam Fall :20 pm -5:20 pm PHYS 223 Dec. 17, Cover Sheet
EE538 Final Exam Fall 005 3:0 pm -5:0 pm PHYS 3 Dec. 17, 005 Cover Sheet Test Duration: 10 minutes. Open Book but Closed Notes. Calculators ARE allowed!! This test contains five problems. Each of the five
More informationLecture 3: Linear Filters
Lecture 3: Linear Filters Professor Fei Fei Li Stanford Vision Lab 1 What we will learn today? Images as functions Linear systems (filters) Convolution and correlation Discrete Fourier Transform (DFT)
More informationVALLIAMMAI ENGINEERING COLLEGE. SRM Nagar, Kattankulathur DEPARTMENT OF INFORMATION TECHNOLOGY. Academic Year
VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur- 603 203 DEPARTMENT OF INFORMATION TECHNOLOGY Academic Year 2016-2017 QUESTION BANK-ODD SEMESTER NAME OF THE SUBJECT SUBJECT CODE SEMESTER YEAR
More information/ (2π) X(e jω ) dω. 4. An 8 point sequence is given by x(n) = {2,2,2,2,1,1,1,1}. Compute 8 point DFT of x(n) by
Code No: RR320402 Set No. 1 III B.Tech II Semester Regular Examinations, Apr/May 2006 DIGITAL SIGNAL PROCESSING ( Common to Electronics & Communication Engineering, Electronics & Instrumentation Engineering,
More informationChirp Transform for FFT
Chirp Transform for FFT Since the FFT is an implementation of the DFT, it provides a frequency resolution of 2π/N, where N is the length of the input sequence. If this resolution is not sufficient in a
More informationPractical Spectral Estimation
Digital Signal Processing/F.G. Meyer Lecture 4 Copyright 2015 François G. Meyer. All Rights Reserved. Practical Spectral Estimation 1 Introduction The goal of spectral estimation is to estimate how the
More informationFourier Series Representation of
Fourier Series Representation of Periodic Signals Rui Wang, Assistant professor Dept. of Information and Communication Tongji University it Email: ruiwang@tongji.edu.cn Outline The response of LIT system
More informationR13 SET - 1
R13 SET - 1 III B. Tech II Semester Regular Examinations, April - 2016 DIGITAL SIGNAL PROCESSING (Electronics and Communication Engineering) Time: 3 hours Maximum Marks: 70 Note: 1. Question Paper consists
More informationLinear Convolution Using FFT
Linear Convolution Using FFT Another useful property is that we can perform circular convolution and see how many points remain the same as those of linear convolution. When P < L and an L-point circular
More informationContents. Digital Signal Processing, Part II: Power Spectrum Estimation
Contents Digital Signal Processing, Part II: Power Spectrum Estimation 5. Application of the FFT for 7. Parametric Spectrum Est. Filtering and Spectrum Estimation 7.1 ARMA-Models 5.1 Fast Convolution 7.2
More informationParametric Method Based PSD Estimation using Gaussian Window
International Journal of Engineering Trends and Technology (IJETT) Volume 29 Number 1 - November 215 Parametric Method Based PSD Estimation using Gaussian Window Pragati Sheel 1, Dr. Rajesh Mehra 2, Preeti
More informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More informationEE538 Final Exam Fall 2007 Mon, Dec 10, 8-10 am RHPH 127 Dec. 10, Cover Sheet
EE538 Final Exam Fall 2007 Mon, Dec 10, 8-10 am RHPH 127 Dec. 10, 2007 Cover Sheet Test Duration: 120 minutes. Open Book but Closed Notes. Calculators allowed!! This test contains five problems. Each of
More informationModule 3. Convolution. Aim
Module Convolution Digital Signal Processing. Slide 4. Aim How to perform convolution in real-time systems efficiently? Is convolution in time domain equivalent to multiplication of the transformed sequence?
More informationRadar Systems Engineering Lecture 3 Review of Signals, Systems and Digital Signal Processing
Radar Systems Engineering Lecture Review of Signals, Systems and Digital Signal Processing Dr. Robert M. O Donnell Guest Lecturer Radar Systems Course Review Signals, Systems & DSP // Block Diagram of
More information1. Calculation of the DFT
ELE E4810: Digital Signal Processing Topic 10: The Fast Fourier Transform 1. Calculation of the DFT. The Fast Fourier Transform algorithm 3. Short-Time Fourier Transform 1 1. Calculation of the DFT! Filter
More informationExercises in Digital Signal Processing
Exercises in Digital Signal Processing Ivan W. Selesnick September, 5 Contents The Discrete Fourier Transform The Fast Fourier Transform 8 3 Filters and Review 4 Linear-Phase FIR Digital Filters 5 5 Windows
More informationIntroduction to Computer Vision. 2D Linear Systems
Introduction to Computer Vision D Linear Systems Review: Linear Systems We define a system as a unit that converts an input function into an output function Independent variable System operator or Transfer
More informationLab 4: Quantization, Oversampling, and Noise Shaping
Lab 4: Quantization, Oversampling, and Noise Shaping Due Friday 04/21/17 Overview: This assignment should be completed with your assigned lab partner(s). Each group must turn in a report composed using
More informationEEO 401 Digital Signal Processing Prof. Mark Fowler
EEO 401 Digital Signal Processing Prof. Mark Fowler Note Set #21 Using the DFT to Implement FIR Filters Reading Assignment: Sect. 7.3 of Proakis & Manolakis Motivation: DTFT View of Filtering There are
More informationDFT-Based FIR Filtering. See Porat s Book: 4.7, 5.6
DFT-Based FIR Filtering See Porat s Book: 4.7, 5.6 1 Motivation: DTFT View of Filtering There are two views of filtering: * Time Domain * Frequency Domain x[ X f ( θ ) h[ H f ( θ ) Y y[ = h[ * x[ f ( θ
More informationTest 2 Electrical Engineering Bachelor Module 8 Signal Processing and Communications
Test 2 Electrical Engineering Bachelor Module 8 Signal Processing and Communications (201400432) Tuesday May 26, 2015, 14:00-17:00h This test consists of three parts, corresponding to the three courses
More informationCONTENTS NOTATIONAL CONVENTIONS GLOSSARY OF KEY SYMBOLS 1 INTRODUCTION 1
DIGITAL SPECTRAL ANALYSIS WITH APPLICATIONS S.LAWRENCE MARPLE, JR. SUMMARY This new book provides a broad perspective of spectral estimation techniques and their implementation. It concerned with spectral
More informationECSE 512 Digital Signal Processing I Fall 2010 FINAL EXAMINATION
FINAL EXAMINATION 9:00 am 12:00 pm, December 20, 2010 Duration: 180 minutes Examiner: Prof. M. Vu Assoc. Examiner: Prof. B. Champagne There are 6 questions for a total of 120 points. This is a closed book
More informationLab 9a. Linear Predictive Coding for Speech Processing
EE275Lab October 27, 2007 Lab 9a. Linear Predictive Coding for Speech Processing Pitch Period Impulse Train Generator Voiced/Unvoiced Speech Switch Vocal Tract Parameters Time-Varying Digital Filter H(z)
More informationCommunications and Signal Processing Spring 2017 MSE Exam
Communications and Signal Processing Spring 2017 MSE Exam Please obtain your Test ID from the following table. You must write your Test ID and name on each of the pages of this exam. A page with missing
More informationMultidimensional digital signal processing
PSfrag replacements Two-dimensional discrete signals N 1 A 2-D discrete signal (also N called a sequence or array) is a function 2 defined over thex(n set 1 of, n 2 ordered ) pairs of integers: y(nx 1,
More informationFILTER DESIGN FOR SIGNAL PROCESSING USING MATLAB AND MATHEMATICAL
FILTER DESIGN FOR SIGNAL PROCESSING USING MATLAB AND MATHEMATICAL Miroslav D. Lutovac The University of Belgrade Belgrade, Yugoslavia Dejan V. Tosic The University of Belgrade Belgrade, Yugoslavia Brian
More informationConvolution. Define a mathematical operation on discrete-time signals called convolution, represented by *. Given two discrete-time signals x 1, x 2,
Filters Filters So far: Sound signals, connection to Fourier Series, Introduction to Fourier Series and Transforms, Introduction to the FFT Today Filters Filters: Keep part of the signal we are interested
More informationLesson 1. Optimal signalbehandling LTH. September Statistical Digital Signal Processing and Modeling, Hayes, M:
Lesson 1 Optimal Signal Processing Optimal signalbehandling LTH September 2013 Statistical Digital Signal Processing and Modeling, Hayes, M: John Wiley & Sons, 1996. ISBN 0471594318 Nedelko Grbic Mtrl
More informationFilter Banks II. Prof. Dr.-Ing. G. Schuller. Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany
Filter Banks II Prof. Dr.-Ing. G. Schuller Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany Page Modulated Filter Banks Extending the DCT The DCT IV transform can be seen as modulated
More informationResponses of Digital Filters Chapter Intended Learning Outcomes:
Responses of Digital Filters Chapter Intended Learning Outcomes: (i) Understanding the relationships between impulse response, frequency response, difference equation and transfer function in characterizing
More informationOne Dimensional Convolution
Dagon University Research Journal 0, Vol. 4 One Dimensional Convolution Myint Myint Thein * Abstract The development of multi-core computers means that the characteristics of digital filters can be rapidly
More informationContinuous and Discrete Time Signals and Systems
Continuous and Discrete Time Signals and Systems Mrinal Mandal University of Alberta, Edmonton, Canada and Amir Asif York University, Toronto, Canada CAMBRIDGE UNIVERSITY PRESS Contents Preface Parti Introduction
More informationSignals & Systems interaction in the Time Domain. (Systems will be LTI from now on unless otherwise stated)
Signals & Systems interaction in the Time Domain (Systems will be LTI from now on unless otherwise stated) Course Objectives Specific Course Topics: -Basic test signals and their properties -Basic system
More informationHST.582J / 6.555J / J Biomedical Signal and Image Processing Spring 2007
MIT OpenCourseWare http://ocw.mit.edu HST.58J / 6.555J / 16.456J Biomedical Signal and Image Processing Spring 007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationFilter Banks II. Prof. Dr.-Ing Gerald Schuller. Fraunhofer IDMT & Ilmenau Technical University Ilmenau, Germany
Filter Banks II Prof. Dr.-Ing Gerald Schuller Fraunhofer IDMT & Ilmenau Technical University Ilmenau, Germany Prof. Dr.-Ing. G. Schuller, shl@idmt.fraunhofer.de Page Modulated Filter Banks Extending the
More informationSyllabus for IMGS-616 Fourier Methods in Imaging (RIT #11857) Week 1: 8/26, 8/28 Week 2: 9/2, 9/4
IMGS 616-20141 p.1 Syllabus for IMGS-616 Fourier Methods in Imaging (RIT #11857) 3 July 2014 (TENTATIVE and subject to change) Note that I expect to be in Europe twice during the term: in Paris the week
More information! Introduction. ! Discrete Time Signals & Systems. ! Z-Transform. ! Inverse Z-Transform. ! Sampling of Continuous Time Signals
ESE 531: Digital Signal Processing Lec 25: April 24, 2018 Review Course Content! Introduction! Discrete Time Signals & Systems! Discrete Time Fourier Transform! Z-Transform! Inverse Z-Transform! Sampling
More informationBasic Design Approaches
(Classic) IIR filter design: Basic Design Approaches. Convert the digital filter specifications into an analog prototype lowpass filter specifications. Determine the analog lowpass filter transfer function
More informationCourse Name: Digital Signal Processing Course Code: EE 605A Credit: 3
Course Name: Digital Signal Processing Course Code: EE 605A Credit: 3 Prerequisites: Sl. No. Subject Description Level of Study 01 Mathematics Fourier Transform, Laplace Transform 1 st Sem, 2 nd Sem 02
More informationMultirate signal processing
Multirate signal processing Discrete-time systems with different sampling rates at various parts of the system are called multirate systems. The need for such systems arises in many applications, including
More informationBiomedical Signal Processing and Signal Modeling
Biomedical Signal Processing and Signal Modeling Eugene N. Bruce University of Kentucky A Wiley-lnterscience Publication JOHN WILEY & SONS, INC. New York Chichester Weinheim Brisbane Singapore Toronto
More informationA New Twist to Fourier Transforms
Hamish D. Meikle A New Twist to Fourier Transforms WILEY VCH WILEY-VCH Verlag GmbH &, Co. KGaA Table ofcontents 1 The Fourier Transform and the Helix 1 1.1 Fourier Transform Conventions 1 1.1.1 Fourier
More informationQuestion Bank. UNIT 1 Part-A
FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY Senkottai Village, Madurai Sivagangai Main Road, Madurai -625 020 An ISO 9001:2008 Certified Institution Question Bank DEPARTMENT OF ELECTRONICS AND COMMUNICATION
More informationConvolution Spatial Aliasing Frequency domain filtering fundamentals Applications Image smoothing Image sharpening
Frequency Domain Filtering Correspondence between Spatial and Frequency Filtering Fourier Transform Brief Introduction Sampling Theory 2 D Discrete Fourier Transform Convolution Spatial Aliasing Frequency
More informationStability Condition in Terms of the Pole Locations
Stability Condition in Terms of the Pole Locations A causal LTI digital filter is BIBO stable if and only if its impulse response h[n] is absolutely summable, i.e., 1 = S h [ n] < n= We now develop a stability
More informationCCNY. BME I5100: Biomedical Signal Processing. Stochastic Processes. Lucas C. Parra Biomedical Engineering Department City College of New York
BME I5100: Biomedical Signal Processing Stochastic Processes Lucas C. Parra Biomedical Engineering Department CCNY 1 Schedule Week 1: Introduction Linear, stationary, normal - the stuff biology is not
More informationSignal Modeling Techniques in Speech Recognition. Hassan A. Kingravi
Signal Modeling Techniques in Speech Recognition Hassan A. Kingravi Outline Introduction Spectral Shaping Spectral Analysis Parameter Transforms Statistical Modeling Discussion Conclusions 1: Introduction
More informationTheory and Problems of Signals and Systems
SCHAUM'S OUTLINES OF Theory and Problems of Signals and Systems HWEI P. HSU is Professor of Electrical Engineering at Fairleigh Dickinson University. He received his B.S. from National Taiwan University
More informationTopic 7. Convolution, Filters, Correlation, Representation. Bryan Pardo, 2008, Northwestern University EECS 352: Machine Perception of Music and Audio
Topic 7 Convolution, Filters, Correlation, Representation Short time Fourier Transform Break signal into windows Calculate DFT of each window The Spectrogram spectrogram(y,1024,512,1024,fs,'yaxis'); A
More informationLecture 19 IIR Filters
Lecture 19 IIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/5/10 1 General IIR Difference Equation IIR system: infinite-impulse response system The most general class
More informationIT6502 DIGITAL SIGNAL PROCESSING Unit I - SIGNALS AND SYSTEMS Basic elements of DSP concepts of frequency in Analo and Diital Sinals samplin theorem Discrete time sinals, systems Analysis of discrete time
More information1 1.27z z 2. 1 z H 2
E481 Digital Signal Processing Exam Date: Thursday -1-1 16:15 18:45 Final Exam - Solutions Dan Ellis 1. (a) In this direct-form II second-order-section filter, the first stage has
More informationLecture Wigner-Ville Distributions
Introduction to Time-Frequency Analysis and Wavelet Transforms Prof. Arun K. Tangirala Department of Chemical Engineering Indian Institute of Technology, Madras Lecture - 6.1 Wigner-Ville Distributions
More informationE2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)
E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,
More informationEEG- Signal Processing
Fatemeh Hadaeghi EEG- Signal Processing Lecture Notes for BSP, Chapter 5 Master Program Data Engineering 1 5 Introduction The complex patterns of neural activity, both in presence and absence of external
More informationChapter 7: IIR Filter Design Techniques
IUST-EE Chapter 7: IIR Filter Design Techniques Contents Performance Specifications Pole-Zero Placement Method Impulse Invariant Method Bilinear Transformation Classical Analog Filters DSP-Shokouhi Advantages
More informationCONVOLUTION & PRODUCT THEOREMS FOR FRFT CHAPTER INTRODUCTION [45]
CHAPTER 3 CONVOLUTION & PRODUCT THEOREMS FOR FRFT M ANY properties of FRFT were derived, developed or established earlier as described in previous chapter. This includes multiplication property, differentiation
More informationQuestion Paper Code : AEC11T02
Hall Ticket No Question Paper Code : AEC11T02 VARDHAMAN COLLEGE OF ENGINEERING (AUTONOMOUS) Affiliated to JNTUH, Hyderabad Four Year B. Tech III Semester Tutorial Question Bank 2013-14 (Regulations: VCE-R11)
More informationMITOCW watch?v=jtj3v Rx7E
MITOCW watch?v=jtj3v Rx7E The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free. To
More informationEE 225D LECTURE ON DIGITAL FILTERS. University of California Berkeley
University of California Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences Professors : N.Morgan / B.Gold EE225D Digital Filters Spring,1999 Lecture 7 N.MORGAN
More informationCentre for Mathematical Sciences HT 2018 Mathematical Statistics
Lund University Stationary stochastic processes Centre for Mathematical Sciences HT 2018 Mathematical Statistics Computer exercise 1 in Stationary stochastic processes, HT 18. The purpose of this computer
More informationz-transforms Definition of the z-transform Chapter
z-transforms Chapter 7 In the study of discrete-time signal and systems, we have thus far considered the time-domain and the frequency domain. The z- domain gives us a third representation. All three domains
More informationMultirate Digital Signal Processing
Multirate Digital Signal Processing Basic Sampling Rate Alteration Devices Up-sampler - Used to increase the sampling rate by an integer factor Down-sampler - Used to decrease the sampling rate by an integer
More informationDiscrete-Time Signals & Systems
Chapter 2 Discrete-Time Signals & Systems 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 2-1-1 Discrete-Time Signals: Time-Domain Representation (1/10) Signals
More informationDiscrete-Time Signals and Systems
ECE 46 Lec Viewgraph of 35 Discrete-Time Signals and Systems Sequences: x { x[ n] }, < n
More informationProbability and Statistics for Final Year Engineering Students
Probability and Statistics for Final Year Engineering Students By Yoni Nazarathy, Last Updated: May 24, 2011. Lecture 6p: Spectral Density, Passing Random Processes through LTI Systems, Filtering Terms
More informationBME 50500: Image and Signal Processing in Biomedicine. Lecture 5: Correlation and Power-Spectrum CCNY
1 BME 50500: Image and Signal Processing in Biomedicine Lecture 5: Correlation and Power-Spectrum Lucas C. Parra Biomedical Engineering Department CCNY http://bme.ccny.cuny.edu/faculty/parra/teaching/signal-and-image/
More informationTherefore the new Fourier coefficients are. Module 2 : Signals in Frequency Domain Problem Set 2. Problem 1
Module 2 : Signals in Frequency Domain Problem Set 2 Problem 1 Let be a periodic signal with fundamental period T and Fourier series coefficients. Derive the Fourier series coefficients of each of the
More informationLaboratory Project 2: Spectral Analysis and Optimal Filtering
Laboratory Project 2: Spectral Analysis and Optimal Filtering Random signals analysis (MVE136) Mats Viberg and Lennart Svensson Department of Signals and Systems Chalmers University of Technology 412 96
More informationFrequency Domain Speech Analysis
Frequency Domain Speech Analysis Short Time Fourier Analysis Cepstral Analysis Windowed (short time) Fourier Transform Spectrogram of speech signals Filter bank implementation* (Real) cepstrum and complex
More informationComputer Vision Lecture 3
Computer Vision Lecture 3 Linear Filters 03.11.2015 Bastian Leibe RWTH Aachen http://www.vision.rwth-aachen.de leibe@vision.rwth-aachen.de Demo Haribo Classification Code available on the class website...
More informationSummary notes for EQ2300 Digital Signal Processing
Summary notes for EQ3 Digital Signal Processing allowed aid for final exams during 6 Joakim Jaldén, 6-- Prerequisites The DFT and the FFT. The discrete Fourier transform The discrete Fourier transform
More informationDIGITAL SIGNAL PROCESSING UNIT 1 SIGNALS AND SYSTEMS 1. What is a continuous and discrete time signal? Continuous time signal: A signal x(t) is said to be continuous if it is defined for all time t. Continuous
More informationLINEAR-PHASE FIR FILTERS DESIGN
LINEAR-PHASE FIR FILTERS DESIGN Prof. Siripong Potisuk inimum-phase Filters A digital filter is a minimum-phase filter if and only if all of its zeros lie inside or on the unit circle; otherwise, it is
More informationMultimedia Signals and Systems - Audio and Video. Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2
Multimedia Signals and Systems - Audio and Video Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2 Kunio Takaya Electrical and Computer Engineering University of Saskatchewan December
More informationCommunication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi
Communication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture - 3 Brief Review of Signals and Systems My subject for today s discussion
More informationLecture 11 FIR Filters
Lecture 11 FIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/4/12 1 The Unit Impulse Sequence Any sequence can be represented in this way. The equation is true if k ranges
More informationExamples. 2-input, 1-output discrete-time systems: 1-input, 1-output discrete-time systems:
Discrete-Time s - I Time-Domain Representation CHAPTER 4 These lecture slides are based on "Digital Signal Processing: A Computer-Based Approach, 4th ed." textbook by S.K. Mitra and its instructor materials.
More informationLAB 6: FIR Filter Design Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 6: FIR Filter Design Summer 011
More informationIndex. p, lip, 78 8 function, 107 v, 7-8 w, 7-8 i,7-8 sine, 43 Bo,94-96
p, lip, 78 8 function, 107 v, 7-8 w, 7-8 i,7-8 sine, 43 Bo,94-96 B 1,94-96 M,94-96 B oro!' 94-96 BIro!' 94-96 I/r, 79 2D linear system, 56 2D FFT, 119 2D Fourier transform, 1, 12, 18,91 2D sinc, 107, 112
More informationADSP ADSP ADSP ADSP. Advanced Digital Signal Processing (18-792) Spring Fall Semester, Department of Electrical and Computer Engineering
Advanced Digital Signal rocessing (18-792) Spring Fall Semester, 201 2012 Department of Electrical and Computer Engineering ROBLEM SET 8 Issued: 10/26/18 Due: 11/2/18 Note: This problem set is due Friday,
More informationChapter 7: Filter Design 7.1 Practical Filter Terminology
hapter 7: Filter Design 7. Practical Filter Terminology Analog and digital filters and their designs constitute one of the major emphasis areas in signal processing and communication systems. This is due
More informationUNIT 4: DIGITAL SYSTEM MODELS
UNIT 4: DIGITAL SYSTEM MODELS 4.1 Introduction This unit is concerned with the description of digital systems, it introduces the concepts of a linear time-invariant system, convolution, the general system
More informationLecture 17: variance in a band = log(s xx (f)) df (2) If we want to plot something that is more directly representative of variance, we can try this:
UCSD SIOC 221A: (Gille) 1 Lecture 17: Recap We ve now spent some time looking closely at coherence and how to assign uncertainties to coherence. Can we think about coherence in a different way? There are
More informationDiscrete-time signals and systems
Discrete-time signals and systems 1 DISCRETE-TIME DYNAMICAL SYSTEMS x(t) G y(t) Linear system: Output y(n) is a linear function of the inputs sequence: y(n) = k= h(k)x(n k) h(k): impulse response of the
More information