Topics in Time-Frequency Analysis

Total Page:16

File Type:pdf, Size:1020Kb

Topics in Time-Frequency Analysis Kernels Multitaper basics MT spectrograms Concentration Topics in time-frequency analysis Maria Sandsten Lecture 7 Stationary and non-stationary spectral analysis February 10 2020 1 / 45 Kernels Multitaper basics MT spectrograms Concentration Wigner distribution The Wigner distribution is defined as, Z 1 τ ∗ τ −i2πf τ Wz (t; f ) = z(t + )z (t − )e dτ; −∞ 2 2 with t is time, f is frequency and z(t) is the analytic signal. The ambiguity function is defined as Z 1 τ ∗ τ −i2πνt Az (ν; τ) = z(t + )z (t − )e dt; −∞ 2 2 where ν is the doppler frequency and τ is the lag variable. 2 / 45 Kernels Multitaper basics MT spectrograms Concentration Ambiguity function Facts about the ambiguity function: I The auto-terms are always centered at ν = τ = 0. I The phase of the ambiguity function oscillations are related to centre t0 and f0 of the signal. I jAz1 (ν; τ)j = jAz0 (ν; τ)j when Wz1 (t; f ) = Wz0 (t − t0; f − f0). I The cross-terms are always located away from the centre. I The placement of the cross-terms are related to the interconnection of the auto-terms. 3 / 45 Kernels Multitaper basics MT spectrograms Concentration The quadratic class The quadratic or Cohen's class is usually defined as Z 1 Z 1 Z 1 Q τ ∗ τ i2π(νt−f τ−νu) Wz (t; f ) = z(u+ )z (u− )φ(ν; τ)e du dτ dν: −∞ −∞ −∞ 2 2 The filtered ambiguity function can also be formulated as Q Az (ν; τ) = Az (ν; τ) · φ(ν; τ); and the corresponding smoothed Wigner distribution as Q Wz (t; f ) = Wz (t; f ) ∗ ∗Φ(t; f ); using the corresponding time-frequency kernel. 4 / 45 Kernels Multitaper basics MT spectrograms Concentration Properties of the ambiguity kernel In order to preserve the marginals and the total energy condition of the time-frequency distribution, the ambiguity kernel must fulfill φ(0; τ) = φ(ν; 0) = 1; φ(0; 0) = 1: For a real-valued distribution, φ∗(−ν; −τ) = φ(ν; τ): 5 / 45 Kernels Multitaper basics MT spectrograms Concentration The Choi-Williams distribution The Choi-Williams distribution or the Exponential distribution (ED), is defined from the ambiguity kernel 2 2 − ν τ φED (ν; τ) = e σ ; where σ is a design parameter. The distribution satisfies the marginals and total energy condition, φED (0; τ) = φED (ν; 0) = 1; φED (0; 0) = 1: The Choi-Williams distribution is real-valued, 2 2 ∗ − ν τ φED (−ν; −τ) = e σ = φED (ν; τ): The Choi-Williams distribution works well for suppression of cross-terms found on a diagonal line between the auto-terms. 6 / 45 Kernels Multitaper basics MT spectrograms Concentration Example 7 / 45 Kernels Multitaper basics MT spectrograms Concentration Separable kernels Separable kernels is defined by φ(ν; τ) = G1(ν) · g2(τ); The smoothed Wigner-Ville distribution becomes F Wz (t; f ) = Wz (t; f ) ∗ ∗Φ(t; f ) = g1(t) ∗ Wz (t; f ) ∗ G2(f ): I If G1(ν) = 1 a doppler-independent kernel is given as φ(ν; τ) = g2(τ), resulting in the pseudo-Wigner distribution, with smoothing only in frequency. I If g2(τ) = 1 a lag-independent kernel is φ(ν; τ) = G1(ν) where the time-frequency formulation gives smoothing only in time. 8 / 45 Kernels Multitaper basics MT spectrograms Concentration Doppler-independent kernel 9 / 45 Kernels Multitaper basics MT spectrograms Concentration Lag-independent kernel 10 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation The Rihaczek distribution The energy of a complex-valued signal is Z 1 Z 1 E = jz(t)j2dt = z(t)z∗(t)dt −∞ −∞ Z 1 Z 1 = z(t)Z ∗(f )e−i2πft dfdt: −∞ −∞ The Rihaczek distribution is defined as ∗ −i2πft Rz (t; f ) = z(t)Z (f )e : 11 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation The Rihaczek distribution The marginals are satisfied as Z 1 Z 1 ∗ −i2πft 2 Rz (t; f )df = z(t) Z (f )e df = jz(t)j ; −∞ −∞ and Z 1 Z 1 ∗ −i2πft 2 Rz (t; f )dt = Z (f ) z(t)e dt = jZ(f )j : −∞ −∞ The ambiguity kernel is found as φ(ν; τ) = e−iπντ : (Assignment!) From this expression we also easily verify that the marginals are satisfied as φ(ν; 0) = φ(0; τ) = 1. Is the Rihaczek distribution real-valued? (φ∗(−ν; −τ) = φ(ν; τ)?): 12 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation Example 13 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation Example 14 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation The spectrogram again We can also formulate the quadratic class as Z 1 Z 1 Q −i2πf τ Wz (t; f ) = rz (u; τ)ρ(t − u; τ)e du dτ; −∞ −∞ where the time-lag kernel is defined as Z 1 ρ(t; τ) = φ(ν; τ)ei2πνt dν; −∞ and the IAF as τ τ r (t; τ) = z(t + )z∗(t − ): z 2 2 (Assignment!) 15 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation The spectrogram kernel The spectrogram also belongs to the quadratic class which can be seen from the following derivation: Z 1 ∗ −i2πft1 2 Sz (t; f ) = j h (t1 − t)z(t1)e dt1j ; −∞ Z 1 Z 1 ∗ −i2πft1 ∗ i2πft2 = ( h (t1 − t)z(t1)e dt1)( h(t2 − t)z (t2)e dt2): −∞ −∞ τ τ With t1 = u + 2 and t2 = u − 2 , we get Sz (t; f ) = Z 1 Z 1 τ τ τ τ = z(u+ )z∗(u− )h∗(u+ −t)h(u− −t)e−i2πf τ dτdu: −∞ −∞ 2 2 2 2 16 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation The spectrogram kernel We can identify the time-lag distribution τ τ r (u; τ) = z(u + )z∗(u − ); z 2 2 and the time-lag kernel τ τ ρ (t − u; τ) = h∗(u − t + )h(u − t − ); h 2 2 giving Z 1 Z 1 −i2πf τ Sz (t; f ) = rz (u; τ)ρh(t − u; τ)e dτdu; −∞ −∞ which we recognize as the quadratic class. 17 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation The spectrogram kernel The window function can be expressed as an ambiguity kernel φh(ν; τ) which is the ambiguity function of a window h(t), i.e., Z 1 ∗ τ τ −i2πνt φh(ν; τ) = h (−t + )h(−t − )e dt: −∞ 2 2 The spectrogram does not fulfill the marginals. Why? 18 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation Example Wigner−Ville Spectrogram Choi−Williams 1 1 1 0.8 0.8 0.8 0.6 0.6 0.6 0.4 0.4 0.4 0.2 0.2 0.2 0 0 0 20 20 20 0.1 0.1 0.1 0 0.05 0 0.05 0 0.05 0 0 0 −20 −0.05 −20 −0.05 −20 −0.05 τ −0.1 ν τ −0.1 ν τ −0.1 ν 19 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation Multitaper spectrogram analysis Replacement with u = (t1 + t2)=2 and τ = t1 − t2 in Z 1 Z 1 Q −i2πf τ Wz (t; f ) = rz (u; τ)ρ(t − u; τ)e du dτ; −∞ −∞ give ZZ t1 + t2 t1 + t2 W Q (t; f ) = r ( ; t −t )ρ(t− ; t −t )e−i2πf (t1−t2)dt dt ; z z 2 1 2 2 1 2 1 2 where the IAF t + t t + t t − t t + t t − t r ( 1 2 ; t −t ) = z( 1 2 + 1 2 )z∗( 1 2 − 1 2 ) = z(t )z∗(t ): z 2 1 2 2 2 2 2 1 2 20 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation Multitaper spectrogram analysis Define a rotated time-lag kernel as t + t ρrot (t ; t ) = ρ( 1 2 ; t − t ); 1 2 2 1 2 and we get ZZ Q ∗ rot −i2πft1 i2πft2 Wz (t; f ) = z(t1)z (t2)ρ (t − t1; t2)e e dt1dt2: rot rot ∗ If the kernel ρ (t1; t2) = (ρ (t2; t1)) , then solving the integral Z rot ρ (t1; t2)u(t1)dt1 = λ u(t2); results in eigenvalues λk and eigenfunctions uk (t). 21 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation Multitaper spectrogram analysis The quadratic class for any kernel is now rewritten as a weighted sum of spectrograms, 1 Z Q X −i2πft1 ∗ 2 Wz (t; f ) = λk j z(t1)e uk (t − t1)dt1j : k=1 The quadratic class of time-frequency distributions can always be approximated, K Z Q X −i2πft1 ∗ 2 Wz (t; f ) ≈ λk j z(t1)e uk (t − t1)dt1j ; k=1 where K is chosen as the number including all λk > . If K is reasonable small, the computations are efficient. 22 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation Classification of time-frequency representations M h X T Sz (t; f ) = Φ (t; f )∗∗ Wz (t; f ) ≈ sj ·uj (t)∗ vj (f ) ∗ Wz (t; f ) ; j=1 Johan Brynolfsson and Maria Sandsten, ”Classification of One-Dimensional Non-Stationary Signals Using the Wigner-Ville Distribution in Convolutional Neural Networks, EUSIPCO, Kos island, Greece, 2017. 23 / 45 Kernels Multitaper basics MT spectrograms Concentration Rihaczek dist. Spectrogram kernel Multitaper interpretation Classification of time-frequency representations Eight different classes of multi-component chirp signals (class 1 and 2), sinusoids (class 3-6) and multi-component short Gaussian signals (class 7 and 8).
Recommended publications
  • Configuring Spectrogram Views
    Chapter 5: Spectrographic Analysis Configuring spectrographic views To create a new spectrogram, spectrogram slice, or selection spectrum view, click on the appropriate New View button in the view toolbar (Figure 5.3) or choose a view type from the View > New menu. Figure 5.3. New View buttons Figure 5.3. The New View buttons, in Raven’s view toolbar. 112 Raven Pro 1.4 User’s Manual Chapter 5: Spectrographic Analysis A dialog box appears, containing parameters for configuring the requested type of spectrographic view (Figure 5.4). The dialog boxes for configuring spectrogram and spectrogram slice views are identical, except for their titles. The dialog boxes are identical because both view types calculate a spectrogram of the entire sound; the only difference between spectrogram and spectrogram slice views is in how the data are displayed (see “How the spectrographic views are related” on page 110). The dialog box for configuring a selection spectrum view is the same, except that it lacks the Averaging parameter. The remainder of this section explains each of the parameters in the configuration dialog box. Figure 5.4 Configure Spectrogram dialog Figure 5.4. The Configure New Spectrogram dialog box. Window type Each data record is multiplied by a window function before its spectrum is calculated. Window functions are used to reduce the magnitude of spurious “sidelobe” energy that appears at frequencies flanking each analysis frequency in a spectrum. These sidelobes appear as a result of Raven Pro 1.4 User’s Manual 113 Chapter 5: Spectrographic Analysis analyzing a finite (truncated) portion of a signal.
    [Show full text]
  • Windowing Techniques, the Welch Method for Improvement of Power Spectrum Estimation
    Computers, Materials & Continua Tech Science Press DOI:10.32604/cmc.2021.014752 Article Windowing Techniques, the Welch Method for Improvement of Power Spectrum Estimation Dah-Jing Jwo1, *, Wei-Yeh Chang1 and I-Hua Wu2 1Department of Communications, Navigation and Control Engineering, National Taiwan Ocean University, Keelung, 202-24, Taiwan 2Innovative Navigation Technology Ltd., Kaohsiung, 801, Taiwan *Corresponding Author: Dah-Jing Jwo. Email: [email protected] Received: 01 October 2020; Accepted: 08 November 2020 Abstract: This paper revisits the characteristics of windowing techniques with various window functions involved, and successively investigates spectral leak- age mitigation utilizing the Welch method. The discrete Fourier transform (DFT) is ubiquitous in digital signal processing (DSP) for the spectrum anal- ysis and can be efciently realized by the fast Fourier transform (FFT). The sampling signal will result in distortion and thus may cause unpredictable spectral leakage in discrete spectrum when the DFT is employed. Windowing is implemented by multiplying the input signal with a window function and windowing amplitude modulates the input signal so that the spectral leakage is evened out. Therefore, windowing processing reduces the amplitude of the samples at the beginning and end of the window. In addition to selecting appropriate window functions, a pretreatment method, such as the Welch method, is effective to mitigate the spectral leakage. Due to the noise caused by imperfect, nite data, the noise reduction from Welch’s method is a desired treatment. The nonparametric Welch method is an improvement on the peri- odogram spectrum estimation method where the signal-to-noise ratio (SNR) is high and mitigates noise in the estimated power spectra in exchange for frequency resolution reduction.
    [Show full text]
  • Windowing Design and Performance Assessment for Mitigation of Spectrum Leakage
    E3S Web of Conferences 94, 03001 (2019) https://doi.org/10.1051/e3sconf/20199403001 ISGNSS 2018 Windowing Design and Performance Assessment for Mitigation of Spectrum Leakage Dah-Jing Jwo 1,*, I-Hua Wu 2 and Yi Chang1 1 Department of Communications, Navigation and Control Engineering, National Taiwan Ocean University 2 Pei-Ning Rd., Keelung 202, Taiwan 2 Bison Electronics, Inc., Neihu District, Taipei 114, Taiwan Abstract. This paper investigates the windowing design and performance assessment for mitigation of spectral leakage. A pretreatment method to reduce the spectral leakage is developed. In addition to selecting appropriate window functions, the Welch method is introduced. Windowing is implemented by multiplying the input signal with a windowing function. The periodogram technique based on Welch method is capable of providing good resolution if data length samples are selected optimally. Windowing amplitude modulates the input signal so that the spectral leakage is evened out. Thus, windowing reduces the amplitude of the samples at the beginning and end of the window, altering leakage. The influence of various window functions on the Fourier transform spectrum of the signals was discussed, and the characteristics and functions of various window functions were explained. In addition, we compared the differences in the influence of different data lengths on spectral resolution and noise levels caused by the traditional power spectrum estimation and various window-function-based Welch power spectrum estimations. 1 Introduction knowledge. Due to computer capacity and limitations on processing time, only limited time samples can actually The spectrum analysis [1-4] is a critical concern in many be processed during spectrum analysis and data military and civilian applications, such as performance processing of signals, i.e., the original signals need to be test of electric power systems and communication cut off, which induces leakage errors.
    [Show full text]
  • Discrete-Time Fourier Transform
    .@.$$-@@4&@$QEG'FC1. CHAPTER7 Discrete-Time Fourier Transform In Chapter 3 and Appendix C, we showed that interesting continuous-time waveforms x(t) can be synthesized by summing sinusoids, or complex exponential signals, having different frequencies fk and complex amplitudes ak. We also introduced the concept of the spectrum of a signal as the collection of information about the frequencies and corresponding complex amplitudes {fk,ak} of the complex exponential signals, and found it convenient to display the spectrum as a plot of spectrum lines versus frequency, each labeled with amplitude and phase. This spectrum plot is a frequency-domain representation that tells us at a glance “how much of each frequency is present in the signal.” In Chapter 4, we extended the spectrum concept from continuous-time signals x(t) to discrete-time signals x[n] obtained by sampling x(t). In the discrete-time case, the line spectrum is plotted as a function of normalized frequency ωˆ . In Chapter 6, we developed the frequency response H(ejωˆ ) which is the frequency-domain representation of an FIR filter. Since an FIR filter can also be characterized in the time domain by its impulse response signal h[n], it is not hard to imagine that the frequency response is the frequency-domain representation,orspectrum, of the sequence h[n]. 236 .@.$$-@@4&@$QEG'FC1. 7-1 DTFT: FOURIER TRANSFORM FOR DISCRETE-TIME SIGNALS 237 In this chapter, we take the next step by developing the discrete-time Fouriertransform (DTFT). The DTFT is a frequency-domain representation for a wide range of both finite- and infinite-length discrete-time signals x[n].
    [Show full text]
  • Understanding the Lomb-Scargle Periodogram
    Understanding the Lomb-Scargle Periodogram Jacob T. VanderPlas1 ABSTRACT The Lomb-Scargle periodogram is a well-known algorithm for detecting and char- acterizing periodic signals in unevenly-sampled data. This paper presents a conceptual introduction to the Lomb-Scargle periodogram and important practical considerations for its use. Rather than a rigorous mathematical treatment, the goal of this paper is to build intuition about what assumptions are implicit in the use of the Lomb-Scargle periodogram and related estimators of periodicity, so as to motivate important practical considerations required in its proper application and interpretation. Subject headings: methods: data analysis | methods: statistical 1. Introduction The Lomb-Scargle periodogram (Lomb 1976; Scargle 1982) is a well-known algorithm for de- tecting and characterizing periodicity in unevenly-sampled time-series, and has seen particularly wide use within the astronomy community. As an example of a typical application of this method, consider the data shown in Figure 1: this is an irregularly-sampled timeseries showing a single ob- ject from the LINEAR survey (Sesar et al. 2011; Palaversa et al. 2013), with un-filtered magnitude measured 280 times over the course of five and a half years. By eye, it is clear that the brightness of the object varies in time with a range spanning approximately 0.8 magnitudes, but what is not immediately clear is that this variation is periodic in time. The Lomb-Scargle periodogram is a method that allows efficient computation of a Fourier-like power spectrum estimator from such unevenly-sampled data, resulting in an intuitive means of determining the period of oscillation.
    [Show full text]
  • Airborne Radar Ground Clutter Suppression Using Multitaper Spectrum Estimation Comparison with Traditional Method
    Airborne Radar Ground Clutter Suppression Using Multitaper Spectrum Estimation Comparison with Traditional Method Linus Ekvall Engineering Physics and Electrical Engineering, master's level 2018 Luleå University of Technology Department of Computer Science, Electrical and Space Engineering Airborne Radar Ground Clutter Suppression Using Multitaper Spectrum Estimation & Comparison with Traditional Method Linus C. Ekvall Lule˚aUniversity of Technology Dept. of Computer Science, Electrical and Space Engineering Div. Signals and Systems 20th September 2018 ABSTRACT During processing of data received by an airborne radar one of the issues is that the typical signal echo from the ground produces a large perturbation. Due to this perturbation it can be difficult to detect targets with low velocity or a low signal-to-noise ratio. Therefore, a filtering process is needed to separate the large perturbation from the target signal. The traditional method include a tapered Fourier transform that operates in parallel with a MTI filter to suppress the main spectral peak in order to produce a smoother spectral output. The difference between a typical signal echo produced from an object in the environment and the signal echo from the ground can be of a magnitude corresponding to more than a 60 dB difference. This thesis presents research of how the multitaper approach can be utilized in concurrence with the minimum variance estimation technique, to produce a spectral estimation that strives for a more effective clutter suppression. A simulation model of the ground clutter was constructed and also a number of simulations for the multitaper, minimum variance estimation technique was made. Compared to the traditional method defined in this thesis, there was a slight improve- ment of the improvement factor when using the multitaper approach.
    [Show full text]
  • The Sliding Window Discrete Fourier Transform Lee F
    The Sliding Window Discrete Fourier Transform Lee F. Richardson Department of Statistics and Data Science, Carnegie Mellon University, Pittsburgh, PA The Sliding Window DFT Sliding Window DFT for Local Periodic Signals Leakage The Sliding Window Discrete Fourier Transform (SW- Leakage occurs when local signal frequency is not one DFT) computes a time-frequency representation of a The Sliding Window DFT is a useful tool for analyzing data with local periodic signals (Richardson and Eddy of the Fourier frequencies (or, the number of complete signal, and is useful for analyzing signals with local (2018)). Since SW-DFT coefficients are complex-numbers, we analyze the squared modulus of these coefficients cycles in a length n window is not an integer). In this periodicities. Shortly described, the SW-DFT takes 2 2 2 case, the largest squared modulus coefficients occur at sequential DFTs of a signal multiplied by a rectangu- |ak,p| = Re(ak,p) + Im(ak,p) frequencies closest to the true frequency: lar sliding window function, where the window func- Since the squared modulus SW-DFT coefficients are large when the window is on a periodic part of the signal: Frequency: 2 Frequency: 2.1 Frequency: 2.2 Frequency: 2.3 Frequency: 2.4 20 Freq 2 20 Freq 2 20 Freq 2 20 Freq 2 20 Freq 2 tion is only nonzero for a short duration. Let x = Freq 3 Freq 3 Freq 3 Freq 3 Freq 3 15 15 15 15 15 [x0, x1, . , xN−1] be a length N signal, then for length 10 10 10 10 10 n ≤ N windows, the SW-DFT of x is: 5 5 5 5 5 0 0 0 0 0 0 50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 1 n−1 X −jk Frequency: 2.5 Frequency: 2.6 Frequency: 2.7 Frequency: 2.8 Frequency: 2.9 ak,p = √ xp−n+1+jω 20 Freq 2 20 Freq 2 20 Freq 2 20 Freq 2 20 Freq 2 n n Freq 3 Freq 3 Freq 3 Freq 3 Freq 3 j=0 15 15 15 15 15 k = 0, 1, .
    [Show full text]
  • Fourier Analysis
    FOURIER ANALYSIS Lucas Illing 2008 Contents 1 Fourier Series 2 1.1 General Introduction . 2 1.2 Discontinuous Functions . 5 1.3 Complex Fourier Series . 7 2 Fourier Transform 8 2.1 Definition . 8 2.2 The issue of convention . 11 2.3 Convolution Theorem . 12 2.4 Spectral Leakage . 13 3 Discrete Time 17 3.1 Discrete Time Fourier Transform . 17 3.2 Discrete Fourier Transform (and FFT) . 19 4 Executive Summary 20 1 1. Fourier Series 1 Fourier Series 1.1 General Introduction Consider a function f(τ) that is periodic with period T . f(τ + T ) = f(τ) (1) We may always rescale τ to make the function 2π periodic. To do so, define 2π a new independent variable t = T τ, so that f(t + 2π) = f(t) (2) So let us consider the set of all sufficiently nice functions f(t) of a real variable t that are periodic, with period 2π. Since the function is periodic we only need to consider its behavior on one interval of length 2π, e.g. on the interval (−π; π). The idea is to decompose any such function f(t) into an infinite sum, or series, of simpler functions. Following Joseph Fourier (1768-1830) consider the infinite sum of sine and cosine functions 1 a0 X f(t) = + [a cos(nt) + b sin(nt)] (3) 2 n n n=1 where the constant coefficients an and bn are called the Fourier coefficients of f. The first question one would like to answer is how to find those coefficients.
    [Show full text]
  • Understanding FFT Overlap Processing Fundamentals
    Understanding FFT Overlap Processing Fundamentals Primer Primer Contents Introduction ....................................................................3 Amplitude Effects .........................................................14 The need for seeing faster time-varying signals .................3 Brief pulses .....................................................................14 The Kaiser Window .........................................................17 Expand the Spectrogram View ......................................4 Calculating Amplitude Reduction .....................................17 How Overlap works ........................................................5 Other windows ................................................................18 Overlapping many Time-to-Frequency transforms .............5 Super-short pulses ..........................................................18 Spectrum and Spectrogram comparisons ....................7 Gaps in Time Lose Signals ..........................................20 Non-overlapped transform processing ..............................7 Actual gaps .....................................................................20 Overlapped time transform processing ..............................7 Effective gaps .................................................................20 Stretching time and impact on Spectrum view ..................8 Overlap – with no gaps ...................................................21 Example: a pseudo-random modulated pulse ...................9 Digital Phosphor Spectrum
    [Show full text]
  • Short-Window Spectral Analysis Using AMVAR and Multitaper Methods: a Comparison
    Biol Cybern (2009) 101:71–80 DOI 10.1007/s00422-009-0318-5 ORIGINAL PAPER Short-window spectral analysis using AMVAR and multitaper methods: a comparison Hariharan Nalatore · Govindan Rangarajan Received: 6 October 2008 / Accepted: 6 May 2009 / Published online: 27 May 2009 © Springer-Verlag 2009 Abstract We compare two popular methods for estimat- even over time intervals of the order of one second, one ing the power spectrum from short data windows, namely must examine them using short time windows of 50–100ms the adaptive multivariate autoregressive (AMVAR) method duration or even shorter. The theoretical basis for divid- and the multitaper method. By analyzing a simulated sig- ing the signal into short time windows lies in treating the nal (embedded in a background Ornstein–Uhlenbeck noise time series within such short windows as being generated process) we demonstrate that the AMVAR method performs by an approximately stationary stochastic process. In other better at detecting short bursts of oscillations compared to words, the full time series is viewed as being generated by the multitaper method. However, both methods are immune a non-stationary stochastic process with locally stationary to jitter in the temporal location of the signal. We also show segments. that coherence can still be detected in noisy bivariate time For the 200Hz sampling frequency used in our study, a series data by the AMVAR method even if the individual duration of 50–100ms corresponds to a data string of only power spectra fail to show any peaks. Finally, using data 11–21 points. The conventional non-parametric approach to from two monkeys performing a visuomotor pattern discrim- spectral analysis, based on the discrete Fourier transform, ination task, we demonstrate that the AMVARmethod is bet- is not effective for such short data lengths.
    [Show full text]
  • Some Notes on the Use of the Windowed Fourier Transform for Spectral Analysis of Discretely Sampled Data
    Submitted to Axioms. Pages 1 - 24. OPEN ACCESS axioms ISSN 2075-1680 www.mdpi.com/journal/axioms Article Some Notes on the Use of the Windowed Fourier Transform for Spectral Analysis of Discretely Sampled Data Robert W. Johnson Alphawave Research, 29 Stanebrook Ct., Jonesboro, GA 30238, USA Version May 23, 2013 submitted to Axioms. Typeset by LATEX using class file mdpi.cls Abstract: The properties of the Gabor and Morlet transforms are examined with respect to the Fourier analysis of discretely sampled data. Forward and inverse transform pairs based on a fixed window with uniform sampling of the frequency axis can satisfy numerically the energy and reconstruction theorems; however, transform pairs based on a variable window or nonuniform frequency sampling in general do not. Instead of selecting the shape of the window as some function of the central frequency, we propose constructing a single window with unit energy from an arbitrary set of windows which is applied over the entire frequency axis. By virtue of using a fixed window with uniform frequency sampling, such a transform satisfies the energy and reconstruction theorems. The shape of the window can be tailored to meet the requirements of the investigator in terms of time/frequency resolution. The algorithm extends naturally to the case of nonuniform signal sampling without modification beyond identification of the Nyquist interval. Keywords: Fourier transform, Gabor transform, Morlet transform, multiresolution analysis 1. Introduction arXiv:1303.0306v2 [physics.data-an] 22 May 2013 The primary criticism leveled at the use of the continuous wavelet transform for the spectral analysis of discretely sampled data is that it fails to give quantitatively meaningful results.
    [Show full text]
  • B3. Short Time Fourier Transform (STFT)
    B3. Short Time Fourier Transform (STFT) Objectives: • Understand the concept of a time varying frequency spectrum and the spectrogram • Understand the effect of different windows on the spectrogram; • Understand the effects of the window length on frequency and time resolutions. 1. Introduction VIDEO: Short Time Fourier Transform (19:24) http://faculty.nps.edu/rcristi/eo3404/b-discrete-fourier- transform/videos/b3_01_timeFrequencyAnalysis.mp4 In a lot of applications, signal carry information and information changes with time. Unless we talk about a beacon, or some sort of synchronizing tone, most signals of interest present characteristics which change with time. In particular the frequency composition (ie the frequency spectrum) most of the time is time varying: just think of a piece of music or a sound from a speaker. The very fact that the pitch and the tone changes with time makes the signal interesting and suitable to carry information. It would be very boring to listen to a recording playing the same notes over and over or listening to someone who keeps repeating the same sound. In this section we address the problem of representing the instantaneous spectrum of a signal. This is can be done as a simple extension of the Discrete Fourier Transform (DFT) introduced in the previous section, applied to a window “sliding” on the signal. The end result is the spectrogram, which shows the evolution of frequencies in time. This information is very usefull in the analysis of a signal, since it gives a sort of signature in the time and frequency domain. Like a music score, as shown in the figure below, describes how the musical notes evolve in time, so the spectrogram shows how the various frequency components of a signal evolve with time.
    [Show full text]