II. Nonparametric Spectrum Estimation for Stationary Random Signals - Non-parametric Methods -
|
|
- Ruth Patrick
- 5 years ago
- Views:
Transcription
1 II. onparametric Spectrum Estimation for Stationary Random Signals - on-parametric Methods - - [p. 3] Periodogram - [p. 12] Periodogram properties - [p. 23] Modified periodogram - [p. 25] Bartlett s method - [p. 30] Blackman-Tukey procedure - [p. 33] Welch s procedure - [p. 36] Minimum Variance MV) Method - [p. 42] MultiTaper Method MTM) 04/01/09 EC4440.SpFY09/MPF -Section II 1
2 -Goal: To estimate the Power Spectral Density PSD) of a wss process -Two main types of approaches eist: -on parametric methods: do not rely on data model -based on FT and FFT, or -Parametric methods: -based on different parametric models of the data: -Moving Average: MA, -AutoRegressive: AR, -AutoRegressive & Moving Average: ARMA -based on subspace methods 04/01/09 EC4440.SpFY09/MPF -Section II 2
3 I. Classical methods use FT of the data sequence or its correlation sequence 1) Periodogram Easy to compute Limited ability to produce accurate PS estimate Recall PSD defined as: jω ) ) P e = k = R k e jωk In practice: R k) is approimated with estimated correlation sequence R k ) 04/01/09 EC4440.SpFY09/MPF -Section II 3
4 Recall: Biased correlation estimate of n), n=0, -1, defined as 1 k 1 * Rˆ k) = + k) ) 0 k< 1 k= 0 ˆ j jωk P e = R k e L ω ) ˆ ) L k = L L 10% PSD estimate sometimes known as correlogram Assume L=-1 1 j jωk Pˆ e = R k e ω ) ˆ ) k= + 1 periodogram 04/01/09 EC4440.SpFY09/MPF -Section II 4
5 Recall 1 * Rˆ k) = k )* k) for k ) 0 0 k< 1 Define: DFT ) = ) ) ) n n w n X k ˆ 1 P e = FT k * k 1 jω) jω) * = X e X e jω ) * ) ) Periodogram defined in terms of data directly 1 P e X e ˆ j j ω) ω = ) 2 04/01/09 EC4440.SpFY09/MPF -Section II 5
6 How to Compute the Periodogram Rectangular window of length DFT ) = ) ) ) n n w n X k j πk ) ) ) ˆ =, = 0,... 1 P e X k k Zero padding n) in the time-domain increases the number of frequency PSD points available 04/01/09 EC4440.SpFY09/MPF -Section II 6
7 Eample: n) white noise ) Rk = j ) Pe ω = [Hayes] 04/01/09 EC4440.SpFY09/MPF -Section II 7
8 Average behavior of the periodogram: 1 { ˆ jω )} ˆ jkω E P e = E R k) e + 1 = 1 k= + 1 { ˆ )} E R k e 1 k 1 ER k = E + k { )} * ) ) = 0 jω k k = R ˆ k ) k k = R k) wb k) wb k) = 0 ow 04/01/09 EC4440.SpFY09/MPF -Section II 8
9 1 j jω k { ˆ ω )} ˆ ) k= + 1 { } E P e = E R k e k { )} k ER k = R k) wb k), wb k) = 0 ow 1 { ˆ jω )} ) ) B k= + 1 ) E P e = R k w k e jωk 1 E P e P e W e 2π { ˆ j ω )} j )* j ) = ω ω B 1 sin ω /2) sin ω/ 2) 2 04/01/09 EC4440.SpFY09/MPF -Section II 9
10 Eample: n) 1 tone in white noise ) = sin ω + φ) + ) n A n v n 0 { jω } 1 jω jω EP e ) = Pe ) We B ) 2π 2 2 A j 0) j 0) v W B e ωω + ) W B e ωω = σ + + ) 4 [Hayes] 04/01/09 EC4440.SpFY09/MPF -Section II 10
11 Eample: set of tones in white noise n) n) = sin2π*150*t) + 2*sin2*π*140*t) + 0.1*randnsizet)); f s =1KHz Obtained using the MATLAB function periodogram.m periodogramn,[],'twosided',1024,fs); plots the two sided periodogram note: better to stick to one-sided when dealing with real-data) [Mathworks] 04/01/09 EC4440.SpFY09/MPF -Section II 11
12 Properties of the Periodogram 1) Periodogram is biased when dealing with finite windows 2) Periodogram is asymptotically unbiased i.e., lim j j E ω P e ) ω + = P e ) 3) Periodogram has limited ability to resolve closely spaced narrowband components limiting factor due to the rectangular window, the shorter the data length, the worse off) not good) 4) Values of the periodogram spaced in frequency by integer multiples of 2π/ are approimately uncorrelated. As data length increases, rate of fluctuation in the periodogram increases not good). 5) Variance of the periodogram does not decrease as data length increases not good). In statistical terms, it means the periodogram is not a consistent estimator of the PSD. evertheless, the periodogram can be a useful tool for spectral estimation in situations where SR is high, and especially if the data is long. 04/01/09 EC4440.SpFY09/MPF -Section II 12
13 n = sin2*pi*150*t) + 2*sin2*pi*140*t) + 0.1*randnsizet));; fs=1000; nfft=1024; 100 samples Comments Periodogram has limited ability to resolve closely spaced narrowband components limiting factor due to the rectangular window, the shorter the data length, the worse off) not good) 67 samples [MathWorks] 04/01/09 EC4440.SpFY09/MPF -Section II 13
14 Comments: Variance of the periodogram does not decrease as the data length increases not good). For white noise n) one can show [Therrien, p.590] jω P e ) lim var 0 + jω 4 var P e ) = σ [Hayes] 04/01/09 EC4440.SpFY09/MPF -Section II 14
15 Comments: As data length increases, rate of fluctuation in the periodogram increases not good). [Kay] 04/01/09 EC4440.SpFY09/MPF -Section II 15
16 Connections between Periodogram and filterbank Periodogram can be viewed as obtained by using a filterbank structure. Recall periodogram defined for n), n=0, -1, as: 1 ) ) ) DFT ) j2 πkp/ = = ) p= 0 n n w n X k pe j 2πk ) ) 1 P e X ) k k 2 ˆ =, = 0, /01/09 EC4440.SpFY09/MPF -Section II 16
17 For n) defined over a finite time window {n-+1), n)} 1 ) ) ) DFT ) = = j2 πkp/ ) p= 0 n n w n X k n pe 1 p= 0 h p) ) ) ) DFT n = n w n X k) = n ph ) p) k=0, -1 04/01/09 EC4440.SpFY09/MPF -Section II 17
18 For n) defined over a finite time window {n-+1), n)} 1 ) ) ) DFT n = n w n X k) = n ph ) p) p= 0 j πk ) ) ) 2 Pˆ e = X k 2 1 k=0, -1 n) Filter H z) 1 hp ) = e - j2 π kp / 2 2 π )) ˆ j k P e Periodogram DFT filterbank 04/01/09 EC4440.SpFY09/MPF -Section II 18
19 Filter H z) 1 h0 p) =, k = π 0 )) ˆ j P e n) Filter H z) 1 h p = e k = - j2 π p/ 1 ), π )) ˆ j P e Filter H z) 1 h p = e k = - j2 π p 1)/ 1 ), 1 04/01/09 EC4440.SpFY09/MPF -Section II ) )) ˆ j P e π j πk ) ) ) 2 Pˆ e = X k 2 1
20 n) Filter H z) 1 h0 p) =, k = 0 Filter H z) 1 h p = e k = - j2 π p/ 1 ), 1 Filter H z) 1 h p = e k = 2 2 π 0 )) ˆ j P e - j2 π p 1)/ 1 ), 1 04/01/09 EC4440.SpFY09/MPF -Section II 20 2 Lowpass filter 2 π )) ˆ j P e Modulation term 2 2 1) )) ˆ j P e π j πk ) ) ) 2 Pˆ e = X k 2 1
21 Filter H z) 1 h0 p) =, k = π 0 )) ˆ j P e Lowpass filter n) Filter H z) h p h p e k - j2 π p/ 1 ) = 0 ), = π )) ˆ j P e Filter H z) h p = h p e k = - j2 π p 1)/ 1 ) 0 ), 1 Modulation term 2 2 1) )) ˆ j P e π j πk ) ) ) 2 Pˆ e = X k /01/09 EC4440.SpFY09/MPF -Section II 21
22 Periodogram DFT matri formulation X ) j2 / j2 1)/ X 1) π π 1 e e = 2 j2 π 1)/ j2 π 1) / X 1) 1 e e 1/ n ) 1/ n 1) 1/ n + 1) j πk ) ) ) 2 Pˆ e = X k 2 1 May be viewed as a rectangular window applied to data before FT operation 04/01/09 EC4440.SpFY09/MPF -Section II 22
23 2. Modified Periodogram can apply a general window to n) prior to computing the periodogram periodogram etension) + ω 1 P e ) = n) w R n) e ˆ j jnω n = 2 ˆ 1 E P e = P e * WR e 2π jω ) jω ) jω ) 2 ote: amount of smoothing is determined by the window type Matlab implementation: Hs = spectrum.periodogram'hamming'); create modified periodogram object psdhs,n,'fs',fs,'fft',1024,'spectrumtype', onesided') plot 04/01/09 EC4440.SpFY09/MPF -Section II 23
24 Eample: n = sin2*pi*150*t) + 2*sin2*pi*140*t) + 0.1*randnsizet));; fs=1000; nfft=1024; data length=100. Hamming window Rectangular window ote: lower Hamming window sidelobes, but wider mainlobes with the Hamming window than the rect. window [MathWorks] 04/01/09 EC4440.SpFY09/MPF -Section II 24
25 3. Bartlett s Method average several periodograms rectangular window applied to each) leads to consistent estimate of the PSD Define: K ˆ 1 P ˆ B e = P e K jω ) k ) jω ) k = 1 where K: number of segments in the data. K L 1 jω 1 jnω PB e ) = n+ k 1) L) e KL k= 1 n= 0 2 n) L points L points L points n 1 n) 2 n) K n) n n n 04/01/09 EC4440.SpFY09/MPF -Section II 25
26 var ote: ˆ jω ) 1 ˆ k var ) jω P ) B e P e K 1 jω P ) e K Eamples: white noise [Hayes] two sinusoids in noise 04/01/09 EC4440.SpFY09/MPF -Section II 26
27 Compare: basic periodograms Periodogram of white noise with variance 1. a) Overlay 50 plots, =64 & b) average c) Overlay 50 plots, =128 & d) average e) Overlay 50 plots, =256 & f) average [Hayes, Fig. 8.9] 04/01/09 EC4440.SpFY09/MPF -Section II 27
28 Bartlett estimates Spectrum estimation of white noise with unit variance a) 50-plot overlay, =512, b) ensemble average. c) Overlay of 50 Bartlett estimates with K=4 & L=128, =512, d) ensemble average. e) Overlay of 50 Bartlett estimates with K=8 & L=64, =512, f) ensemble average. K: number of segments L: segment length [Hayes, Fig. 8.14] 04/01/09 EC4440.SpFY09/MPF -Section II 28
29 Spectrum estimation of 2 tones in white noise a) 50-plot overlay, =512, b) ensemble average. c) Overlay of 50 Bartlett estimates with K=4 & L=128, =512, d) ensemble average. e) Overlay of 50 Bartlett estimates with K=8 & L=64, =512, f) ensemble average. ote: Spectral peaks broadening, why? π φ) π φ ) n) = Asin 0.2 n+ + sin 0.25 n+ + vn ) 1 2 [Hayes, Fig. 8.15] 04/01/09 EC4440.SpFY09/MPF -Section II 29
30 4. Blackman-Tukey BT) Procedure apply windowing to correlation samples prior to transformation to reduce variance of the estimator M ˆ jω P ) ) ˆ ) BT e = w k R k e k= M jωk where the window wk) length is much smaller than data length, wk) = 0 for k > M where M << Maimum recommended window length M</5 04/01/09 EC4440.SpFY09/MPF -Section II 30
31 n ) = 10 ep2 π j0.15) n) + 20 ep2 π j0.2) n) + wn ) wn )~ 0,0.1) Data length=20 Periodogram BT PSD estimate Triangular window, length=4 [Kay, Fig.4.9] 04/01/09 EC4440.SpFY09/MPF -Section II 31
32 n ) = 10 ep2 π j0.15) n) + 20 ep2 π j0.2) n) + wn ) Periodogram Data length=100 wn )~ 0,0.1) BT PSD estimate Triangular window, length=20 [Kay, Fig.4.10] 04/01/09 EC4440.SpFY09/MPF -Section II 32
33 5. Welch s Procedure Modification to the Bartlett procedure: Averaging modified periodograms Dataset split into K possibly overlapping segments of length L window applied to each segment all K periodograms are averaged where K ˆ 1 P ˆ w e = P e K jω ) k ) jω ) k = 1 ˆ 1 P w n n e L 1 ) ) ) 2 k k jωn = ) n= 0 window data n) located in K th segment MATLAB Hs = spectrum.welch'rectangular', segmentlength, overlap %); define Welch spectrum object psdhs,n,'fs',fs,'fft',512); plot 04/01/09 EC4440.SpFY09/MPF -Section II 33
34 [Hayes] ote: what do we gain with the Welch s method? Reduction in spectral leakage that takes place thru the sidelobes of the data window 04/01/09 EC4440.SpFY09/MPF -Section II 34
35 Performance comparison variability resolution Figure of merit Periodogram π/) π/) Bartlett 1/K 0.89 K2π/) π/) Welch 50% overlap, triang. window) [Table 8.7, Hayes] 9/8K π/L) π/) BT 2M/ π/M) π/) K: number of data segments : overall data length L: data section length M: BT window length All schemes limited by data length Definitions: { jω var P )} e variability normalized variance) 2 jω E P e ) { } figure of merit variability resolution resolution should be as small as possible. mainlobe width at its 1/2 power 6dB) down level 04/01/09 EC4440.SpFY09/MPF -Section II 35
36 6. Minimum Variance MV) Spectrum Estimate Capon PSD PSD estimated by a filterbank where each filter of length P<= data length) adapts itself to data characteristics, i.e., - has center frequency f i spread between 0 and f s ) - has bandwidth Δ P filters spread over [0,f s ] Δ 1/f s ) - depends on R - designed so components at f i are passed without distortion and rejects maimum amount of out-of-band power ˆ jω ) P PMV e =, H 1 e R e e= e e jω j P 1) ω T [1,,..., ] [Manolakis, p.474] 04/01/09 EC4440.SpFY09/MPF -Section II 36
37 MATLAB implementation ˆ ω) P, [1, ω,..., 1) ω P ] MV e = e= e e H 1 e R e ote: j j j P T H 1 H H 1 H 1 H e R e= e UΛ U ) e= e UΛ U e H 1 H H = ) ) e U Λ e U H jω j P 1) ω eu= [1, e,..., e ] u1 u [ u u ] = FT ),, FT ) 1 P P 04/01/09 EC4440.SpFY09/MPF -Section II 37
38 H 1 H 1 H H = ) Λ ) e R e e U e U [ FT u ),, FT u )] 1 1 [ FT u ),, FT u )] 1 = 1/ λ ) FT u ) + +1/ λ ) FT u ) H = Λ e R e P P P 1 2 = 1/ λk ) FT uk ) k=1 P P H Pˆ MV jω e ) = P k=1 1/ λ ) FT u ) k P k 2 04/01/09 EC4440.SpFY09/MPF -Section II 38
39 Application: n e e wn wn j2 0.1) ) 10 π n j π n = + ) + ), )~ 0,1) Filter order M=P Data length = averaged trials MV: solid line Bartlett method: small dashed line All-pole: large dashed line Δ f = 0.02 Wind. Length ~ 50 [Manolakis, e ] 04/01/09 EC4440.SpFY09/MPF -Section II 39
40 MATLAB Implementation R=covar,P); [v,d]=eig ; UI=diaginvd)+eps); VI=absfftv,nfft)).^2; P=10log10P)-10log10VI*UI); In practice, - pick P < - zeropad FFT to insure accurate measurements 04/01/09 EC4440.SpFY09/MPF -Section II 40
41 7. MultiTaper Method MTM) Spectrum Estimate Thompson) Periodogram can be viewed as obtained by using a filterbank structure, using rectangular windows). Multitaper method MTM) can be viewed the same way, with different filters, i.e., different windows. These optimal FIR filters are derived from a set of sequences known as discrete prolate spheroidal sequences. A set of independent estimates of the power spectrum is computed, by premultiplying the data by orthogonal tapered windows designed to minimize spectral leakage. [Victorin] [Mathworks, statistical signal processing, spectral analysis, non parametric methods] 04/01/09 EC4440.SpFY09/MPF -Section II 41
42 7. MultiTaper Method MTM) cont PSD variance is reduced by averaging over different tapered windows using the full data. Since the data length is not shortened, the overall PSD bias is smaller than that obtained when splitting data in segments. MTM method provides a time-bandwidth parameter, W, with which to balance the variance and resolution. W is directly related to the number of tapers used to compute the spectrum. There are always 2*W-1 tapers used to form the estimate. As W increases, - there are more estimates of the power spectrum, and the variance of the estimate decreases, - each estimate ehibits more spectral leakage i.e., wider peaks) and the overall spectral estimate is more biased. For each data set, there is usually a value for W that allows an optimal trade-off between bias and variance. [Mathworks, statistical signal processing, spectral analysis, non parametric methods] 04/01/09 EC4440.SpFY09/MPF -Section II 42
43 Periodogram fs = 1000Hz; =1000; 2 sines, Amplitude = [1 2]; Sinusoid frequencies f = [150;140]; Hs1 = spectrum.mtm4,'adapt'); psdhs1,n,'fs',fs,'fft',1024) MTM, W=4 MTM, W=3/2 04/01/09 EC4440.SpFY09/MPF -Section II 43
44 References: [Manolakis] [Hayes] Statistical digital signal processing and modeling, M. M. Hayes, 1996, Wiley. [Kay] Modern spectral estimation, S. M. Kay, 1988, Prentice Hall. [Victorin] A. Victorin, Multi-Taper Method for Spectral Analysis and Signal Reconstruction of Solar Wind Data, MS Plasma Physics, Sept. 2007, Royal Institute of Technology, Sweeden, EE_2007_054.pdf 04/01/09 EC4440.SpFY09/MPF -Section II 44
Practical 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 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 informationSignal processing Frequency analysis
Signal processing Frequency analysis Jean-Hugh Thomas (jean-hugh.thomas@univ-lemans.r) Fourier series and Fourier transorm (h30 lecture+h30 practical work) 2 Sampling (h30+h30) 3 Power spectrum estimation
More informationSPECTRUM. Deterministic Signals with Finite Energy (l 2 ) Deterministic Signals with Infinite Energy N 1. n=0. N N X N(f) 2
SPECTRUM Deterministic Signals with Finite Energy (l 2 ) Energy Spectrum: S xx (f) = X(f) 2 = 2 x(n)e j2πfn n= Deterministic Signals with Infinite Energy DTFT of truncated signal: X N (f) = N x(n)e j2πfn
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science : Discrete-Time Signal Processing
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.34: Discrete-Time Signal Processing OpenCourseWare 006 ecture 8 Periodogram Reading: Sections 0.6 and 0.7
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 informationDiscrete Fourier Transform
Discrete Fourier Transform Virtually all practical signals have finite length (e.g., sensor data, audio records, digital images, stock values, etc). Rather than considering such signals to be zero-padded
More informationFourier Analysis of Signals Using the DFT
Fourier Analysis of Signals Using the DFT ECE 535 Lecture April 29, 23 Overview: Motivation Many applications require analyzing the frequency content of signals Speech processing study resonances of vocal
More informationINTRODUCTION TO DELTA-SIGMA ADCS
ECE37 Advanced Analog Circuits INTRODUCTION TO DELTA-SIGMA ADCS Richard Schreier richard.schreier@analog.com NLCOTD: Level Translator VDD > VDD2, e.g. 3-V logic? -V logic VDD < VDD2, e.g. -V logic? 3-V
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 informationComplement on Digital Spectral Analysis and Optimal Filtering: Theory and Exercises
Complement on Digital Spectral Analysis and Optimal Filtering: Theory and Exercises Random Signals Analysis (MVE136) Mats Viberg Department of Signals and Systems Chalmers University of Technology 412
More informationPeriodogram and Correlogram Methods. Lecture 2
Periodogram and Correlogram Methods Lecture 2 Lecture notes to accompany Introduction to Spectral Analysis Slide L2 1 Periodogram Recall 2nd definition of (!): (!) = lim N!1 E 8 >< >: 1 N NX t=1 y(t)e
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 informationECE 636: Systems identification
ECE 636: Systems identification Lectures 7 8 onparametric identification (continued) Important distributions: chi square, t distribution, F distribution Sampling distributions ib i Sample mean If the variance
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 informationComplement on Digital Spectral Analysis and Optimal Filtering: Theory and Exercises
Complement on Digital Spectral Analysis and Optimal Filtering: Theory and Exercises Random Processes With Applications (MVE 135) Mats Viberg Department of Signals and Systems Chalmers University of Technology
More informationSNR Calculation and Spectral Estimation [S&T Appendix A]
SR Calculation and Spectral Estimation [S&T Appendix A] or, How not to make a mess of an FFT Make sure the input is located in an FFT bin 1 Window the data! A Hann window works well. Compute the FFT 3
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 informationL6: Short-time Fourier analysis and synthesis
L6: Short-time Fourier analysis and synthesis Overview Analysis: Fourier-transform view Analysis: filtering view Synthesis: filter bank summation (FBS) method Synthesis: overlap-add (OLA) method STFT magnitude
More informationEEO 401 Digital Signal Processing Prof. Mark Fowler
EEO 4 Digital Signal Processing Pro. Mark Fowler ote Set # Using the DFT or Spectral Analysis o Signals Reading Assignment: Sect. 7.4 o Proakis & Manolakis Ch. 6 o Porat s Book /9 Goal o Practical Spectral
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 informationSystem Identification
System Identification Lecture 4: Transfer function averaging and smoothing Roy Smith 28-- 4. Averaging Multiple estimates Multiple experiments: u r pk, y r pk, r,..., R, and k,..., K. Multiple estimates
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 informationDesign Criteria for the Quadratically Interpolated FFT Method (I): Bias due to Interpolation
CENTER FOR COMPUTER RESEARCH IN MUSIC AND ACOUSTICS DEPARTMENT OF MUSIC, STANFORD UNIVERSITY REPORT NO. STAN-M-4 Design Criteria for the Quadratically Interpolated FFT Method (I): Bias due to Interpolation
More information6.435, System Identification
System Identification 6.435 SET 3 Nonparametric Identification Munther A. Dahleh 1 Nonparametric Methods for System ID Time domain methods Impulse response Step response Correlation analysis / time Frequency
More informationUNIVERSITY OF OSLO. Faculty of mathematics and natural sciences. Forslag til fasit, versjon-01: Problem 1 Signals and systems.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 1th, 016 Examination hours: 14:30 18.30 This problem set
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 informationconvenient means to determine response to a sum of clear evidence of signal properties that are obscured in the original signal
Digital Speech Processing Lecture 9 Short-Time Fourier Analysis Methods- Introduction 1 General Discrete-Time Model of Speech Production Voiced Speech: A V P(z)G(z)V(z)R(z) Unvoiced Speech: A N N(z)V(z)R(z)
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 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 informationBispectral resolution and leakage effect of the indirect bispectrum estimate for different types of 2D window functions
Bispectral resolution and leakage effect of the indirect bispectrum estimate for different types of D window functions Teofil-Cristian OROIAN, Constantin-Iulian VIZITIU, Florin ŞERBAN Communications and
More informationData Processing and Analysis
Data Processing and Analysis Rick Aster and Brian Borchers September 10, 2013 Energy and Power Spectra It is frequently valuable to study the power distribution of a signal in the frequency domain. For
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 information-Digital Signal Processing- FIR Filter Design. Lecture May-16
-Digital Signal Processing- FIR Filter Design Lecture-17 24-May-16 FIR Filter Design! FIR filters can also be designed from a frequency response specification.! The equivalent sampled impulse response
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 informationEEM 409. Random Signals. Problem Set-2: (Power Spectral Density, LTI Systems with Random Inputs) Problem 1: Problem 2:
EEM 409 Random Signals Problem Set-2: (Power Spectral Density, LTI Systems with Random Inputs) Problem 1: Consider a random process of the form = + Problem 2: X(t) = b cos(2π t + ), where b is a constant,
More informationEDISP (NWL3) (English) Digital Signal Processing DFT Windowing, FFT. October 19, 2016
EDISP (NWL3) (English) Digital Signal Processing DFT Windowing, FFT October 19, 2016 DFT resolution 1 N-point DFT frequency sampled at θ k = 2πk N, so the resolution is f s/n If we want more, we use N
More informationIMPROVEMENTS IN MODAL PARAMETER EXTRACTION THROUGH POST-PROCESSING FREQUENCY RESPONSE FUNCTION ESTIMATES
IMPROVEMENTS IN MODAL PARAMETER EXTRACTION THROUGH POST-PROCESSING FREQUENCY RESPONSE FUNCTION ESTIMATES Bere M. Gur Prof. Christopher Niezreci Prof. Peter Avitabile Structural Dynamics and Acoustic Systems
More informationINTRODUCTION TO DELTA-SIGMA ADCS
ECE1371 Advanced Analog Circuits Lecture 1 INTRODUCTION TO DELTA-SIGMA ADCS Richard Schreier richard.schreier@analog.com Trevor Caldwell trevor.caldwell@utoronto.ca Course Goals Deepen understanding of
More information! Spectral Analysis with DFT. ! Windowing. ! Effect of zero-padding. ! Time-dependent Fourier transform. " Aka short-time Fourier transform
Lecture Outline ESE 531: Digital Signal Processing Spectral Analysis with DFT Windowing Lec 24: April 18, 2019 Spectral Analysis Effect of zero-padding Time-dependent Fourier transform " Aka short-time
More informationIII.C - Linear Transformations: Optimal Filtering
1 III.C - Linear Transformations: Optimal Filtering FIR Wiener Filter [p. 3] Mean square signal estimation principles [p. 4] Orthogonality principle [p. 7] FIR Wiener filtering concepts [p. 8] Filter coefficients
More informationParametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d)
Parametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d) Electrical & Computer Engineering North Carolina State University Acknowledgment: ECE792-41 slides
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 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 informationAutoregressive tracking of vortex shedding. 2. Autoregression versus dual phase-locked loop
Autoregressive tracking of vortex shedding Dileepan Joseph, 3 September 2003 Invensys UTC, Oxford 1. Introduction The purpose of this report is to summarize the work I have done in terms of an AR algorithm
More informationSignals and Spectra - Review
Signals and Spectra - Review SIGNALS DETERMINISTIC No uncertainty w.r.t. the value of a signal at any time Modeled by mathematical epressions RANDOM some degree of uncertainty before the signal occurs
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 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 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 informationRandom signals II. ÚPGM FIT VUT Brno,
Random signals II. Jan Černocký ÚPGM FIT VUT Brno, cernocky@fit.vutbr.cz 1 Temporal estimate of autocorrelation coefficients for ergodic discrete-time random process. ˆR[k] = 1 N N 1 n=0 x[n]x[n + k],
More informationCorrelator I. Basics. Chapter Introduction. 8.2 Digitization Sampling. D. Anish Roshi
Chapter 8 Correlator I. Basics D. Anish Roshi 8.1 Introduction A radio interferometer measures the mutual coherence function of the electric field due to a given source brightness distribution in the sky.
More informationIntroduction. Spectral Estimation Overview Periodogram Bias, variance, and distribution Blackman-Tukey Method Welch-Bartlett Method Others
Spectral Estimation Overview Periodogram Bias, variance, and distribution Blackman-Tukey Method Welch-Bartlett Method Others Introduction R x (e jω ) r x (l)e jωl l= Most stationary random processes have
More informationEE 505 Lecture 10. Spectral Characterization. Part 2 of 2
EE 505 Lecture 10 Spectral Characterization Part 2 of 2 Review from last lecture Spectral Analysis If f(t) is periodic f(t) alternately f(t) = = A A ( kω t + ) 0 + Aksin θk k= 1 0 + a ksin t k= 1 k= 1
More informationFGN 1.0 PPL 1.0 FD
FGN PPL FD f f Figure SDFs for FGN, PPL and FD processes (top to bottom rows, respectively) on both linear/log and log/log aes (left- and right-hand columns, respectively) Each SDF S X ( ) is normalized
More informationReview of spectral analysis methods applied to sea level anomaly signals
Review of spectral analysis methods applied to sea level anomaly signals C. Mailhes 1, D. Bonacci 1, O. Besson 1, A. Guillot 2, S. Le Gac 2, N. Steunou 2, C. Cheymol 2, N. Picot 2 1. Telecommunications
More informationPart III Spectrum Estimation
ECE79-4 Part III Part III Spectrum Estimation 3. Parametric Methods for Spectral Estimation Electrical & Computer Engineering North Carolina State University Acnowledgment: ECE79-4 slides were adapted
More informationThe Discrete Fourier Transform
In [ ]: cd matlab pwd The Discrete Fourier Transform Scope and Background Reading This session introduces the z-transform which is used in the analysis of discrete time systems. As for the Fourier and
More informationDistortion Analysis T
EE 435 Lecture 32 Spectral Performance Windowing Spectral Performance of Data Converters - Time Quantization - Amplitude Quantization Quantization Noise . Review from last lecture. Distortion Analysis
More informationParametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes
Parametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes Electrical & Computer Engineering North Carolina State University Acknowledgment: ECE792-41 slides were adapted
More informationTime series models in the Frequency domain. The power spectrum, Spectral analysis
ime series models in the Frequency domain he power spectrum, Spectral analysis Relationship between the periodogram and the autocorrelations = + = ( ) ( ˆ α ˆ ) β I Yt cos t + Yt sin t t= t= ( ( ) ) cosλ
More information2. SPECTRAL ANALYSIS APPLIED TO STOCHASTIC PROCESSES
2. SPECTRAL ANALYSIS APPLIED TO STOCHASTIC PROCESSES 2.0 THEOREM OF WIENER- KHINTCHINE An important technique in the study of deterministic signals consists in using harmonic functions to gain the spectral
More informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Lecture 1 Time-Dependent FT Announcements! Midterm: 2/22/216 Open everything... but cheat sheet recommended instead 1am-12pm How s the lab going? Frequency Analysis with
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 information1. FIR Filter Design
ELEN E4810: Digital Signal Processing Topic 9: Filter Design: FIR 1. Windowed Impulse Response 2. Window Shapes 3. Design by Iterative Optimization 1 1. FIR Filter Design! FIR filters! no poles (just zeros)!
More informationExperimental Fourier Transforms
Chapter 5 Experimental Fourier Transforms 5.1 Sampling and Aliasing Given x(t), we observe only sampled data x s (t) = x(t)s(t; T s ) (Fig. 5.1), where s is called sampling or comb function and can be
More informationcovariance function, 174 probability structure of; Yule-Walker equations, 174 Moving average process, fluctuations, 5-6, 175 probability structure of
Index* The Statistical Analysis of Time Series by T. W. Anderson Copyright 1971 John Wiley & Sons, Inc. Aliasing, 387-388 Autoregressive {continued) Amplitude, 4, 94 case of first-order, 174 Associated
More informationClass of waveform coders can be represented in this manner
Digital Speech Processing Lecture 15 Speech Coding Methods Based on Speech Waveform Representations ti and Speech Models Uniform and Non- Uniform Coding Methods 1 Analog-to-Digital Conversion (Sampling
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 informationChapter 6: Nonparametric Time- and Frequency-Domain Methods. Problems presented by Uwe
System Identification written by L. Ljung, Prentice Hall PTR, 1999 Chapter 6: Nonparametric Time- and Frequency-Domain Methods Problems presented by Uwe System Identification Problems Chapter 6 p. 1/33
More informationDigital Speech Processing Lecture 10. Short-Time Fourier Analysis Methods - Filter Bank Design
Digital Speech Processing Lecture Short-Time Fourier Analysis Methods - Filter Bank Design Review of STFT j j ˆ m ˆ. X e x[ mw ] [ nˆ m] e nˆ function of nˆ looks like a time sequence function of ˆ looks
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 informationELEN E4810: Digital Signal Processing Topic 11: Continuous Signals. 1. Sampling and Reconstruction 2. Quantization
ELEN E4810: Digital Signal Processing Topic 11: Continuous Signals 1. Sampling and Reconstruction 2. Quantization 1 1. Sampling & Reconstruction DSP must interact with an analog world: A to D D to A x(t)
More informationsine wave fit algorithm
TECHNICAL REPORT IR-S3-SB-9 1 Properties of the IEEE-STD-57 four parameter sine wave fit algorithm Peter Händel, Senior Member, IEEE Abstract The IEEE Standard 57 (IEEE-STD-57) provides algorithms for
More informationI. Signals & Sinusoids
I. Signals & Sinusoids [p. 3] Signal definition Sinusoidal signal Plotting a sinusoid [p. 12] Signal operations Time shifting Time scaling Time reversal Combining time shifting & scaling [p. 17] Trigonometric
More informationBasics on 2-D 2 D Random Signal
Basics on -D D Random Signal Spring 06 Instructor: K. J. Ray Liu ECE Department, Univ. of Maryland, College Park Overview Last Time: Fourier Analysis for -D signals Image enhancement via spatial filtering
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 informationDSP. Chapter-3 : Filter Design. Marc Moonen. Dept. E.E./ESAT-STADIUS, KU Leuven
DSP Chapter-3 : Filter Design Marc Moonen Dept. E.E./ESAT-STADIUS, KU Leuven marc.moonen@esat.kuleuven.be www.esat.kuleuven.be/stadius/ Filter Design Process Step-1 : Define filter specs Pass-band, stop-band,
More informationEE 435. Lecture 32. Spectral Performance Windowing
EE 435 Lecture 32 Spectral Performance Windowing . Review from last lecture. Distortion Analysis T 0 T S THEOREM?: If N P is an integer and x(t) is band limited to f MAX, then 2 Am Χ mnp 1 0 m h N and
More informationL29: Fourier analysis
L29: Fourier analysis Introduction The discrete Fourier Transform (DFT) The DFT matrix The Fast Fourier Transform (FFT) The Short-time Fourier Transform (STFT) Fourier Descriptors CSCE 666 Pattern Analysis
More informationIntroduction to Biomedical Engineering
Introduction to Biomedical Engineering Biosignal processing Kung-Bin Sung 6/11/2007 1 Outline Chapter 10: Biosignal processing Characteristics of biosignals Frequency domain representation and analysis
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 informationThe Discrete Fourier Transform (DFT) Properties of the DFT DFT-Specic Properties Power spectrum estimate. Alex Sheremet.
4. April 2, 27 -order sequences Measurements produce sequences of numbers Measurement purpose: characterize a stochastic process. Example: Process: water surface elevation as a function of time Parameters:
More informationELEG 305: Digital Signal Processing
ELEG 305: Digital Signal Processing Lecture : Design of Digital IIR Filters (Part I) Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 008 K. E. Barner (Univ.
More informationRepresentation of Digital Signals
Representation of Digital Signals Representation of Digital Signals Why a special lecture? from NMSOP- Bormann (2002) Representation of Digital Signals Almost every analysis in Geophysics (and meanwhile
More informationrepresentation of speech
Digital Speech Processing Lectures 7-8 Time Domain Methods in Speech Processing 1 General Synthesis Model voiced sound amplitude Log Areas, Reflection Coefficients, Formants, Vocal Tract Polynomial, l
More informationReview of Fundamentals of Digital Signal Processing
Solution Manual for Theory and Applications of Digital Speech Processing by Lawrence Rabiner and Ronald Schafer Click here to Purchase full Solution Manual at http://solutionmanuals.info Link download
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 informationProblem Sheet 1 Examples of Random Processes
RANDOM'PROCESSES'AND'TIME'SERIES'ANALYSIS.'PART'II:'RANDOM'PROCESSES' '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''Problem'Sheets' Problem Sheet 1 Examples of Random Processes 1. Give
More informationLecture 4 - Spectral Estimation
Lecture 4 - Spectral Estimation The Discrete Fourier Transform The Discrete Fourier Transform (DFT) is the equivalent of the continuous Fourier Transform for signals known only at N instants separated
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:305:45 CBC C222 Lecture 8 Frequency Analysis 14/02/18 http://www.ee.unlv.edu/~b1morris/ee482/
More information8 The Discrete Fourier Transform (DFT)
8 The Discrete Fourier Transform (DFT) ² Discrete-Time Fourier Transform and Z-transform are de ned over in niteduration sequence. Both transforms are functions of continuous variables (ω and z). For nite-duration
More informationLABORATORY 3 FINITE IMPULSE RESPONSE FILTERS
LABORATORY 3 FINITE IMPULSE RESPONSE FILTERS 3.. Introduction A digital filter is a discrete system, used with the purpose of changing the amplitude and/or phase spectrum of a signal. The systems (filters)
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 informationSpeech Signal Representations
Speech Signal Representations Berlin Chen 2003 References: 1. X. Huang et. al., Spoken Language Processing, Chapters 5, 6 2. J. R. Deller et. al., Discrete-Time Processing of Speech Signals, Chapters 4-6
More informationThe Selection of Weighting Functions For Linear Arrays Using Different Techniques
The Selection of Weighting Functions For Linear Arrays Using Different Techniques Dr. S. Ravishankar H. V. Kumaraswamy B. D. Satish 17 Apr 2007 Rajarshi Shahu College of Engineering National Conference
More informationFilter Design Problem
Filter Design Problem Design of frequency-selective filters usually starts with a specification of their frequency response function. Practical filters have passband and stopband ripples, while exhibiting
More informationSpectral Analysis. 1.1 Introduction ANDERS I. ERIKSSON
Reprinted from Analysis Methods for Multi-Spacecraft Data Götz Paschmann and Patrick W. Daly (Eds.), ISSI Scientific Report SR-1 (Electronic edition 1.1) c 1998, 2 ISSI/ESA 1 Spectral Analysis ANDERS I.
More informationFrequency estimation by DFT interpolation: A comparison of methods
Frequency estimation by DFT interpolation: A comparison of methods Bernd Bischl, Uwe Ligges, Claus Weihs March 5, 009 Abstract This article comments on a frequency estimator which was proposed by [6] and
More informationDepartment of Electrical and Computer Engineering Digital Speech Processing Homework No. 6 Solutions
Problem 1 Department of Electrical and Computer Engineering Digital Speech Processing Homework No. 6 Solutions The complex cepstrum, ˆx[n], of a sequence x[n] is the inverse Fourier transform of the complex
More informationNonparametric and Parametric Defined This text distinguishes between systems and the sequences (processes) that result when a WN input is applied
Linear Signal Models Overview Introduction Linear nonparametric vs. parametric models Equivalent representations Spectral flatness measure PZ vs. ARMA models Wold decomposition Introduction Many researchers
More information