Spectral Analysis - Introduction
|
|
- Kevin Strickland
- 5 years ago
- Views:
Transcription
1 Spectral Analysis - Introduction Initial Code Introduction Summer Semester 011/01 Lukas Vacha In previous lectures we studied stationary time series in terms of quantities that are functions of time. The value of a variable X t at time t was described as a sequence of innovations. The main focus was on the covariance structure between X t and X t-d at a distinct dates t and t - d. This approach is known as analyzing the properties of a process X t in the time domain. Another approach is to describe the value of a variable (signal) X t as a weighted sum of periodic functions of the form coshwl and sinhwl, where w denotes an angular frequency: X t = m + Ÿ 0 p ahwl coshw tl w + Ÿ 0 p dhwl sinhw tl w The goal is to determine how important are cycles of different frequencies in the variable X t. This is known as frequency domain analysis. The spectrum of a stationary time series hhwl in frequency domain, is the counterpart of a covariance function in the time domain. Theoretically it provides a different representation of the same information, this approach can yield both powerful numerical methods of analysis and new insights. Frequency, Period f is frequency of cycles per second, c.p.s., Hz w = p f - angular frequency (radinas per second) T = 1êf - period The complex exponential A sine wave e i x = coshxl + i sinhxl Nice example of sin function and the angular frequency:
2 QF_II_Lecture 3 Spectra1.cdf Plot@Sin@xD,8x,0,Pi<D Plot@Sin@4xD,8x,0,Pi<D Plot@Sin@1xD,8x,0,Pi<D
3 QF_II_Lecture 3 Spectra1.cdf 3 Examples Plot@Sin@xD+Sin@8xD,8x,0,Pi<D Plot@Sin@xD+Sin@4xD+Sin@1xD+Sin@0xD,8x,0,4Pi<D Sampling The Nyquist frequency is equal to one-half of the sampling frequency. w N Y Q U I S T = T-1, T denotes the sample length. The Nyquist frequency is the highest frequency that can be measured in a signal. An undersampled signal Expel of undersampled signal:
4 4 QF_II_Lecture 3 Spectra1.cdf t/nyquist_limit_java_plugin.html The Fourier transform A transform takes one function (or signal) and turns it into another function (or signal). Any infinite sequence xhtl may be viewed as a combination of periodic functions. The Fourier coefficients of a time series (signal) x(t) are given by XHfL = Ÿ - xhtl e - p i f t d t XHfL determines the relative weight of each complex sinusoids. inverse Fourier transform: xhtl = 1 p Ÿ p -p XHfL e p i f t d f The original sequence is a linear combination of complex sinusoids. Discrete Fourier Transform The Discrete Fourier Transform (DFT) transforms one function into another, which is called the frequency domain representation of the original function (which is often a function in the time domain). The input to the DFT is a finite sequence of real or complex numbers.the only requirements of these conventions are that the DFT and Inverse Discrete Fourier Transform (IDFT) have opposite-sign exponents and that the product of their normalization factors be 1êN. A normalization of 1í N for both the DFT and IDFT makes the transforms unitary, which has some theoretical advantages, but it is often more practical in numerical computation to perform the scaling all at once (and a unit scaling can be convenient in other ways). DFT: X k = N-1 t=0 x t e - p i f k t, k = 0, 1,... N - 1 IDFT:
5 QF_II_Lecture 3 Spectra1.cdf 5 IDFT: x t = 1 N-1 N k=0 X k e p i f k t, t = 0, 1,... N - 1 How Does FT Work FT uses complex exponentials (sinusoid) as building blocks. e j w t = coshw tl + j sinhw tl For each frequency of complex exponential, the sinusoid at that frequency is compared to the signal. If the signal consists of that frequency, the correlation is high ö large FT coefficients. XHfL = Ÿ - xhtl e - p i f t dt If the signal does not have any spectral component at a frequency, the correlation at that frequency is low / zero, ö small / zero FT coefficient. Example FT At Work Remove comments (*...*) for modification of the example signal=table@sin@4xdh*+sin@8xd*l,8x,0,pi,.04<d; ListPlot@signal,PlotLabelØ"signal",FillingØAxis,PlotStyleØ8PointSize@MediumD<D ListPlot@Abs@Fourier@signalDD@@1;;Floor@HLength@signalD-1LêDDD^,PlotLabelØ,Pl signal
6 6 QF_II_Lecture 3 Spectra1.cdf Example with noise signal=table@sin@4xd+1.5*randomreal@normaldistribution@0,1ddh*+sin@8xd*l,8x,0,pi,.04<d ListPlot@signal,PlotLabelØ"signal",FillingØAxis,PlotStyleØ8PointSize@MediumD<D ListPlot@Abs@Fourier@signalDD@@1;;Floor@HLength@signalD-1LêDDD^,PlotLabelØ,Pl signal
7 QF_II_Lecture 3 Spectra1.cdf 7 Example: Complete time series reconstruction I Example shows that the original sequence is a linear combination of complex sinusoids. signal=table@sin@4xd+sin@8xd+sin@4xdh*+randomreal@normaldistribution@0,1dd*l,8x,0,pi ListPlot@signal,PlotLabelØ"Original signal",fillingøaxis,plotstyleø8pointsize@mediumd<d ListPlot@Abs@Fourier@signalDD@@1;;Floor@HLength@signalD-1LêDDD^,PlotLabelØ,Pl InverseFourier@Fourier@signalDD; ListPlot@InverseFourier@Fourier@signalDD,PlotLabelØ"Reconstructed signal",fillingøaxis, 4 Original signal Reconstructed signal Nice wiki example:
8 8 QF_II_Lecture 3 Spectra1.cdf Example: Complete time series reconstruction II (Stationarity) FT identifies all spectral components present in the signal, however it does not provide any information regarding the temporal (time) localization of these components. Why? signal1=table@sin@xd,8x,0,5pi,.1<d; signal=table@sin@0xd,8x,0,5pi,.1<d; signal=join@signal1,signald; ListPlot@signal,PlotLabelØ"Original signal",fillingøaxis,plotstyleø8pointsize@mediumd<d ListPlot@Abs@Fourier@signalDD@@1;;Floor@HLength@signalD-1LêDDD^,PlotLabelØ,Pl H*InverseFourier@Fourier@signalDD;*L ListPlot@InverseFourier@Fourier@signalDD,PlotLabelØ"Reconstructed signal",fillingøaxis, Original signal Reconstructed signal What happens when we use for the reconstruction only the most significant frequencies?
9 QF_II_Lecture 3 Spectra1.cdf 9 Example: time series reconstruction - a real time series There is a problem with the trend!. Why? signal=log@financialdata@"^dji","nov. 1, 009"D@@All,DDD; ListLinePlot@signal,PlotLabelØ"Original signal"d ListPlot@Abs@Fourier@signalDD@@1;;Floor@HLength@signalD-1LêDDD^,PlotLabelØ,Pl InverseFourier@Fourier@signalDD; ListLinePlot@InverseFourier@Fourier@signalDDH*@@1;;10DD*L,PlotLabelØ"Reconstructed sign Original signal Reconstructed signal Parseval's theorem The energy spectral density describes how the energy (variance) of a signal (time series) is distributed with frequency. Total energy of the signal may be obtained by integrating the energy per unit frequency 1ê p XHfL over the p interval. t=- x t = 1 p Ÿ p -p XHfL df The squared magnitude of the Fourier transform X HfL is called the power spectrum, denoted hhwl, (or
10 10 QF_II_Lecture 3 Spectra1.cdf the p interval. t=- x t = 1 p Ÿ p -p XHfL df The squared magnitude of the Fourier transform X HfL is called the power spectrum, denoted hhwl, (or energy-density spectrum, or spectrum) of the signal x t. (Periodogram is another name for essentially the same quantity) Parseval's theorem for the finite sequence: N-1 t=0 x t = 1 N-1 N k=0 X k Estimation of the The spectrum of a stationary time series hhwl is the counterpart of a covariance function in frequency domain. That is, it is the Fourier transform of the covariance function ghkl and vice versa. hhwl = 1 p k=- ghkl e -i k w. p ghkl =Ÿ -p hhwl e i w k w In practice, h`hwl is computed only at the Fourier frequencies: w j = p jên for j œ I1,,..., n è M where n è = dn - 1êt h`iw j M = 1 p n n t=1 xhtl e iw j t = 1 p XIw j M Examples: Random Numbers, ARFIMA ListLinePlot@Abs@Fourier@RandomReal@NormalDistribution@0,1D,300DDD@@1;;Floor@H300-1LêDD
11 QF_II_Lecture 3 Spectra1.cdf 11 ListLinePlot@Abs@Fourier@ARFIMA@0.6,0.4,-0.0,100DDD@@1;;Floor@H100-1LêDDD^,PlotLabelØ Appendix Sample periodogram The periodogram S(w) is a basic estimator for the spectrum (spectral density) of stationary stochastic processes. The periodogram S(w) is often calculated only for N equidistant points in the frequency domain from zero until (N-1)/N Hz. w j = p jêt for j = 1,,..., HT - 1Lê [ a j = T T t=1 xhtl cos Aw j Ht - 1LE [ b j = T T t=1 xhtl sin Aw j Ht - 1LE [ Sx Hw j L = :a j + [ b j >
12 1 QF_II_Lecture 3 Spectra1.cdf Examples: spectral densities of AR and MA processes ü 1) AR signal=arfima@0.9,0,0.0,300d; ListLinePlot@signal,PlotLabelØSignalD ListLinePlot@Abs@Fourier@signalDD@@1;;Floor@H300-1LêDDD^,PlotLabelØ,PlotRange 6 Signal
13 QF_II_Lecture 3 Spectra1.cdf 13 signal=arfima@-0.8,0,0.0,300d; ListLinePlot@signal,PlotLabelØSignalD ListLinePlot@Abs@Fourier@signalDD@@1;;Floor@H300-1LêDDD^,PlotLabelØ,PlotRange 4 Signal
14 14 QF_II_Lecture 3 Spectra1.cdf ü ) MA signal=arfima@0.0,0,0.8,300d; ListLinePlot@signal,PlotLabelØSignalD ListLinePlot@Abs@Fourier@signalDD@@1;;Floor@H300-1LêDDD^,PlotLabelØ,PlotRange 3 Signal
15 QF_II_Lecture 3 Spectra1.cdf 15 signal=arfima@0.0,0,-0.8,300d; ListLinePlot@signal,PlotLabelØSignalD ListLinePlot@Abs@Fourier@signalDD@@1;;Floor@H300-1LêDDD^,PlotLabelØ,PlotRange 3 Signal
16 16 QF_II_Lecture 3 Spectra1.cdf ü MA(q) Needs@"TimeSeries`TimeSeries`"D tsma=timeseries@mamodel@8-0.5, 0.9<, 1D,300D; ListLinePlot@Abs@Fourier@tsMADD@@1;;Floor@HLength@tsMAD-1LêDDD^,PlotLabelØ, tsma=timeseries@mamodel@80.3, 0.4, 0.5, 0.4, 0.3<, 1D,300D; ListLinePlot@Abs@Fourier@tsMADD@@1;;Floor@HLength@tsMAD-1LêDDD^,PlotLabelØ, References: Priestley (1981) Hamilton (1994) Homework #4 Deadline: Monday.4.01, 16:00 1 ) Estimate power spectrum (Periodogram) of simulated AR, MA and ARMA process. ) Estimate power spectrum (Periodogram) of a real financial market time series. 3) to be announced... Homework may be returned in class, or sent via to vachal@utia.cas.cz
SEISMIC WAVE PROPAGATION. Lecture 2: Fourier Analysis
SEISMIC WAVE PROPAGATION Lecture 2: Fourier Analysis Fourier Series & Fourier Transforms Fourier Series Review of trigonometric identities Analysing the square wave Fourier Transform Transforms of some
More informationLecture 10. Applied Econometrics. Time Series Filters. Jozef Barunik ( utia. cas. cz
Applied Econometrics Lecture 10 Time Series Filters Please note that for interactive manipulation you need Mathematica 6 version of this.pdf. Mathematica 6 is available at all Lab's Computers at IE http://staff.utia.cas.cz/baruni
More informationPure Random process Pure Random Process or White Noise Process: is a random process {X t, t 0} which has: { σ 2 if k = 0 0 if k 0
MODULE 9: STATIONARY PROCESSES 7 Lecture 2 Autoregressive Processes 1 Moving Average Process Pure Random process Pure Random Process or White Noise Process: is a random process X t, t 0} which has: E[X
More informationQuantitative Finance I
Quantitative Finance I Linear AR and MA Models (Lecture 4) Winter Semester 01/013 by Lukas Vacha * If viewed in.pdf format - for full functionality use Mathematica 7 (or higher) notebook (.nb) version
More informationSystem Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang
System Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang 1-1 Course Description Emphases Delivering concepts and Practice Programming Identification Methods using Matlab Class
More informationFourier transform. Stefano Ferrari. Università degli Studi di Milano Methods for Image Processing. academic year
Fourier transform Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Methods for Image Processing academic year 27 28 Function transforms Sometimes, operating on a class of functions
More informationEE 574 Detection and Estimation Theory Lecture Presentation 8
Lecture Presentation 8 Aykut HOCANIN Dept. of Electrical and Electronic Engineering 1/14 Chapter 3: Representation of Random Processes 3.2 Deterministic Functions:Orthogonal Representations For a finite-energy
More informationEmpirical Market Microstructure Analysis (EMMA)
Empirical Market Microstructure Analysis (EMMA) Lecture 3: Statistical Building Blocks and Econometric Basics Prof. Dr. Michael Stein michael.stein@vwl.uni-freiburg.de Albert-Ludwigs-University of Freiburg
More informationTutorial Sheet #2 discrete vs. continuous functions, periodicity, sampling
2.39 utorial Sheet #2 discrete vs. continuous functions, periodicity, sampling We will encounter two classes of signals in this class, continuous-signals and discrete-signals. he distinct mathematical
More informationELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process
Department of Electrical Engineering University of Arkansas ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Definition of stochastic process (random
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 informationTSKS01 Digital Communication Lecture 1
TSKS01 Digital Communication Lecture 1 Introduction, Repetition, and Noise Modeling Emil Björnson Department of Electrical Engineering (ISY) Division of Communication Systems Emil Björnson Course Director
More informationContinuous-time Fourier Methods
ELEC 321-001 SIGNALS and SYSTEMS Continuous-time Fourier Methods Chapter 6 1 Representing a Signal The convolution method for finding the response of a system to an excitation takes advantage of the linearity
More informationBCT Lecture 3. Lukas Vacha.
BCT Lecture 3 Lukas Vacha vachal@utia.cas.cz Stationarity and Unit Root Testing Why do we need to test for Non-Stationarity? The stationarity or otherwise of a series can strongly influence its behaviour
More informationIntroduction to Signal Processing
to Signal Processing Davide Bacciu Dipartimento di Informatica Università di Pisa bacciu@di.unipi.it Intelligent Systems for Pattern Recognition Signals = Time series Definitions Motivations A sequence
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 informationTIME SERIES ANALYSIS
2 WE ARE DEALING WITH THE TOUGHEST CASES: TIME SERIES OF UNEQUALLY SPACED AND GAPPED ASTRONOMICAL DATA 3 A PERIODIC SIGNAL Dots: periodic signal with frequency f = 0.123456789 d -1. Dotted line: fit for
More informationE 4101/5101 Lecture 6: Spectral analysis
E 4101/5101 Lecture 6: Spectral analysis Ragnar Nymoen 3 March 2011 References to this lecture Hamilton Ch 6 Lecture note (on web page) For stationary variables/processes there is a close correspondence
More informationStochastic Processes: I. consider bowl of worms model for oscilloscope experiment:
Stochastic Processes: I consider bowl of worms model for oscilloscope experiment: SAPAscope 2.0 / 0 1 RESET SAPA2e 22, 23 II 1 stochastic process is: Stochastic Processes: II informally: bowl + drawing
More informationSTAD57 Time Series Analysis. Lecture 23
STAD57 Time Series Analysis Lecture 23 1 Spectral Representation Spectral representation of stationary {X t } is: 12 i2t Xt e du 12 1/2 1/2 for U( ) a stochastic process with independent increments du(ω)=
More informationTime Series: Theory and Methods
Peter J. Brockwell Richard A. Davis Time Series: Theory and Methods Second Edition With 124 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vn ix CHAPTER 1 Stationary
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 informationSignal Processing COS 323
Signal Processing COS 323 Digital Signals D: functions of space or time e.g., sound 2D: often functions of 2 spatial dimensions e.g. images 3D: functions of 3 spatial dimensions CAT, MRI scans or 2 space,
More informationEC402: Serial Correlation. Danny Quah Economics Department, LSE Lent 2015
EC402: Serial Correlation Danny Quah Economics Department, LSE Lent 2015 OUTLINE 1. Stationarity 1.1 Covariance stationarity 1.2 Explicit Models. Special cases: ARMA processes 2. Some complex numbers.
More information8.2 Harmonic Regression and the Periodogram
Chapter 8 Spectral Methods 8.1 Introduction Spectral methods are based on thining of a time series as a superposition of sinusoidal fluctuations of various frequencies the analogue for a random process
More informationTime Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY
Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY & Contents PREFACE xiii 1 1.1. 1.2. Difference Equations First-Order Difference Equations 1 /?th-order Difference
More informationLecture 2: ARMA(p,q) models (part 2)
Lecture 2: ARMA(p,q) models (part 2) Florian Pelgrin University of Lausanne, École des HEC Department of mathematics (IMEA-Nice) Sept. 2011 - Jan. 2012 Florian Pelgrin (HEC) Univariate time series Sept.
More informationThe Hilbert Transform
The Hilbert Transform David Hilbert 1 ABSTRACT: In this presentation, the basic theoretical background of the Hilbert Transform is introduced. Using this transform, normal real-valued time domain functions
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 informationSystem Identification & Parameter Estimation
System Identification & Parameter Estimation Wb3: SIPE lecture Correlation functions in time & frequency domain Alfred C. Schouten, Dept. of Biomechanical Engineering (BMechE), Fac. 3mE // Delft University
More informationCentre for Mathematical Sciences HT 2017 Mathematical Statistics. Study chapters 6.1, 6.2 and in the course book.
Lund University Stationary stochastic processes Centre for Mathematical Sciences HT 2017 Mathematical Statistics Computer exercise 2 in Stationary stochastic processes, HT 17. The purpose with this computer
More informationSHOCK AND VIBRATION RESPONSE SPECTRA COURSE Unit 6A. The Fourier Transform. By Tom Irvine
SHOCK ND VIBRTION RESPONSE SPECTR COURSE Unit 6. The Fourier Transform By Tom Irvine Introduction Stationary vibration signals can be placed along a continuum in terms of the their qualitative characteristics.
More information2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit dwm/courses/2tf
Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic
More informationAdvanced Econometrics
Advanced Econometrics Marco Sunder Nov 04 2010 Marco Sunder Advanced Econometrics 1/ 25 Contents 1 2 3 Marco Sunder Advanced Econometrics 2/ 25 Music Marco Sunder Advanced Econometrics 3/ 25 Music Marco
More informationFourier Series and Fourier Transforms
Fourier Series and Fourier Transforms EECS2 (6.082), MIT Fall 2006 Lectures 2 and 3 Fourier Series From your differential equations course, 18.03, you know Fourier s expression representing a T -periodic
More information2A1H Time-Frequency Analysis II
2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 209 For any corrections see the course page DW Murray at www.robots.ox.ac.uk/ dwm/courses/2tf. (a) A signal g(t) with period
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 informationω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the
he ime-frequency Concept []. Review of Fourier Series Consider the following set of time functions {3A sin t, A sin t}. We can represent these functions in different ways by plotting the amplitude versus
More informationA523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011
A523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011 Lecture 6 PDFs for Lecture 1-5 are on the web page Problem set 2 is on the web page Article on web page A Guided
More informationEA2.3 - Electronics 2 1
In the previous lecture, I talked about the idea of complex frequency s, where s = σ + jω. Using such concept of complex frequency allows us to analyse signals and systems with better generality. In this
More informationSpectral Analysis. Jesús Fernández-Villaverde University of Pennsylvania
Spectral Analysis Jesús Fernández-Villaverde University of Pennsylvania 1 Why Spectral Analysis? We want to develop a theory to obtain the business cycle properties of the data. Burns and Mitchell (1946).
More informationIntroduction to Spectral and Time-Spectral Analysis with some Applications
Introduction to Spectral and Time-Spectral Analysis with some Applications A.INTRODUCTION Consider the following process X t = U cos[(=3)t] + V sin[(=3)t] Where *E[U] = E[V ] = 0 *E[UV ] = 0 *V ar(u) =
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 informationTime Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY
Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY PREFACE xiii 1 Difference Equations 1.1. First-Order Difference Equations 1 1.2. pth-order Difference Equations 7
More informationTime Series 3. Robert Almgren. Sept. 28, 2009
Time Series 3 Robert Almgren Sept. 28, 2009 Last time we discussed two main categories of linear models, and their combination. Here w t denotes a white noise: a stationary process with E w t ) = 0, E
More informationAdditional exercises in Stationary Stochastic Processes
Mathematical Statistics, Centre or Mathematical Sciences Lund University Additional exercises 8 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
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 informationLecture 3: Autoregressive Moving Average (ARMA) Models and their Practical Applications
Lecture 3: Autoregressive Moving Average (ARMA) Models and their Practical Applications Prof. Massimo Guidolin 20192 Financial Econometrics Winter/Spring 2018 Overview Moving average processes Autoregressive
More informationSpectral Analysis of Random Processes
Spectral Analysis of Random Processes Spectral Analysis of Random Processes Generally, all properties of a random process should be defined by averaging over the ensemble of realizations. Generally, all
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 informationDiscrete-time Fourier Series (DTFS)
Discrete-time Fourier Series (DTFS) Arun K. Tangirala (IIT Madras) Applied Time-Series Analysis 59 Opening remarks The Fourier series representation for discrete-time signals has some similarities with
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 informationat least 50 and preferably 100 observations should be available to build a proper model
III Box-Jenkins Methods 1. Pros and Cons of ARIMA Forecasting a) need for data at least 50 and preferably 100 observations should be available to build a proper model used most frequently for hourly or
More informationPART 1. Review of DSP. f (t)e iωt dt. F(ω) = f (t) = 1 2π. F(ω)e iωt dω. f (t) F (ω) The Fourier Transform. Fourier Transform.
PART 1 Review of DSP Mauricio Sacchi University of Alberta, Edmonton, AB, Canada The Fourier Transform F() = f (t) = 1 2π f (t)e it dt F()e it d Fourier Transform Inverse Transform f (t) F () Part 1 Review
More informationFourier Analysis of Stationary and Non-Stationary Time Series
Fourier Analysis of Stationary and Non-Stationary Time Series September 6, 2012 A time series is a stochastic process indexed at discrete points in time i.e X t for t = 0, 1, 2, 3,... The mean is defined
More informationA time series is called strictly stationary if the joint distribution of every collection (Y t
5 Time series A time series is a set of observations recorded over time. You can think for example at the GDP of a country over the years (or quarters) or the hourly measurements of temperature over a
More informationFourier Analysis Linear transformations and lters. 3. Fourier Analysis. Alex Sheremet. April 11, 2007
Stochastic processes review 3. Data Analysis Techniques in Oceanography OCP668 April, 27 Stochastic processes review Denition Fixed ζ = ζ : Function X (t) = X (t, ζ). Fixed t = t: Random Variable X (ζ)
More informationElements of Multivariate Time Series Analysis
Gregory C. Reinsel Elements of Multivariate Time Series Analysis Second Edition With 14 Figures Springer Contents Preface to the Second Edition Preface to the First Edition vii ix 1. Vector Time Series
More informationClassic Time Series Analysis
Classic Time Series Analysis Concepts and Definitions Let Y be a random number with PDF f Y t ~f,t Define t =E[Y t ] m(t) is known as the trend Define the autocovariance t, s =COV [Y t,y s ] =E[ Y t t
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 informationStochastic Processes. A stochastic process is a function of two variables:
Stochastic Processes Stochastic: from Greek stochastikos, proceeding by guesswork, literally, skillful in aiming. A stochastic process is simply a collection of random variables labelled by some parameter:
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.6 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts
More informationFourier Series and Transform KEEE343 Communication Theory Lecture #7, March 24, Prof. Young-Chai Ko
Fourier Series and Transform KEEE343 Communication Theory Lecture #7, March 24, 20 Prof. Young-Chai Ko koyc@korea.ac.kr Summary Fourier transform Properties Fourier transform of special function Fourier
More informationEmpirical Macroeconomics
Empirical Macroeconomics Francesco Franco Nova SBE April 5, 2016 Francesco Franco Empirical Macroeconomics 1/39 Growth and Fluctuations Supply and Demand Figure : US dynamics Francesco Franco Empirical
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 informationEEL3135: Homework #3 Solutions
EEL335: Homework #3 Solutions Problem : (a) Compute the CTFT for the following signal: xt () cos( πt) cos( 3t) + cos( 4πt). First, we use the trigonometric identity (easy to show by using the inverse Euler
More informationStatistics 349(02) Review Questions
Statistics 349(0) Review Questions I. Suppose that for N = 80 observations on the time series { : t T} the following statistics were calculated: _ x = 10.54 C(0) = 4.99 In addition the sample autocorrelation
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 informationTIME SERIES ANALYSIS. Forecasting and Control. Wiley. Fifth Edition GWILYM M. JENKINS GEORGE E. P. BOX GREGORY C. REINSEL GRETA M.
TIME SERIES ANALYSIS Forecasting and Control Fifth Edition GEORGE E. P. BOX GWILYM M. JENKINS GREGORY C. REINSEL GRETA M. LJUNG Wiley CONTENTS PREFACE TO THE FIFTH EDITION PREFACE TO THE FOURTH EDITION
More informationContinuous Fourier transform of a Gaussian Function
Continuous Fourier transform of a Gaussian Function Gaussian function: e t2 /(2σ 2 ) The CFT of a Gaussian function is also a Gaussian function (i.e., time domain is Gaussian, then the frequency domain
More informationStatistics of Stochastic Processes
Prof. Dr. J. Franke All of Statistics 4.1 Statistics of Stochastic Processes discrete time: sequence of r.v...., X 1, X 0, X 1, X 2,... X t R d in general. Here: d = 1. continuous time: random function
More information13. Power Spectrum. For a deterministic signal x(t), the spectrum is well defined: If represents its Fourier transform, i.e., if.
For a deterministic signal x(t), the spectrum is well defined: If represents its Fourier transform, i.e., if jt X ( ) = xte ( ) dt, (3-) then X ( ) represents its energy spectrum. his follows from Parseval
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 informationLecture 7 ELE 301: Signals and Systems
Lecture 7 ELE 30: Signals and Systems Prof. Paul Cuff Princeton University Fall 20-2 Cuff (Lecture 7) ELE 30: Signals and Systems Fall 20-2 / 22 Introduction to Fourier Transforms Fourier transform as
More informationTime Series Examples Sheet
Lent Term 2001 Richard Weber Time Series Examples Sheet This is the examples sheet for the M. Phil. course in Time Series. A copy can be found at: http://www.statslab.cam.ac.uk/~rrw1/timeseries/ Throughout,
More informationConvergence of Fourier Series
MATH 454: Analysis Two James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University April, 8 MATH 454: Analysis Two Outline The Cos Family MATH 454: Analysis
More informationDeterministic. Deterministic data are those can be described by an explicit mathematical relationship
Random data Deterministic Deterministic data are those can be described by an explicit mathematical relationship Deterministic x(t) =X cos r! k m t Non deterministic There is no way to predict an exact
More informationNon-parametric identification
Non-parametric Non-parametric Transient Step-response using Spectral Transient Correlation Frequency function estimate Spectral System Identification, SSY230 Non-parametric 1 Non-parametric Transient Step-response
More informationChapter 3 - Temporal processes
STK4150 - Intro 1 Chapter 3 - Temporal processes Odd Kolbjørnsen and Geir Storvik January 23 2017 STK4150 - Intro 2 Temporal processes Data collected over time Past, present, future, change Temporal aspect
More informationIntroduction to Fourier Transforms. Lecture 7 ELE 301: Signals and Systems. Fourier Series. Rect Example
Introduction to Fourier ransforms Lecture 7 ELE 3: Signals and Systems Fourier transform as a limit of the Fourier series Inverse Fourier transform: he Fourier integral theorem Prof. Paul Cuff Princeton
More informationIDENTIFICATION OF ARMA MODELS
IDENTIFICATION OF ARMA MODELS A stationary stochastic process can be characterised, equivalently, by its autocovariance function or its partial autocovariance function. It can also be characterised by
More informationRead Section 1.1, Examples of time series, on pages 1-8. These example introduce the book; you are not tested on them.
TS Module 1 Time series overview (The attached PDF file has better formatting.)! Model building! Time series plots Read Section 1.1, Examples of time series, on pages 1-8. These example introduce the book;
More informationAtmospheric Flight Dynamics Example Exam 1 Solutions
Atmospheric Flight Dynamics Example Exam 1 Solutions 1 Question Figure 1: Product function Rūū (τ) In figure 1 the product function Rūū (τ) of the stationary stochastic process ū is given. What can be
More informationChapter 4: Models for Stationary Time Series
Chapter 4: Models for Stationary Time Series Now we will introduce some useful parametric models for time series that are stationary processes. We begin by defining the General Linear Process. Let {Y t
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 informationLecture 1: Fundamental concepts in Time Series Analysis (part 2)
Lecture 1: Fundamental concepts in Time Series Analysis (part 2) Florian Pelgrin University of Lausanne, École des HEC Department of mathematics (IMEA-Nice) Sept. 2011 - Jan. 2012 Florian Pelgrin (HEC)
More informationInformation and Communications Security: Encryption and Information Hiding
Short Course on Information and Communications Security: Encryption and Information Hiding Tuesday, 10 March Friday, 13 March, 2015 Lecture 5: Signal Analysis Contents The complex exponential The complex
More informationLecture # 06. Image Processing in Frequency Domain
Digital Image Processing CP-7008 Lecture # 06 Image Processing in Frequency Domain Fall 2011 Outline Fourier Transform Relationship with Image Processing CP-7008: Digital Image Processing Lecture # 6 2
More informationEmpirical Macroeconomics
Empirical Macroeconomics Francesco Franco Nova SBE April 21, 2015 Francesco Franco Empirical Macroeconomics 1/33 Growth and Fluctuations Supply and Demand Figure : US dynamics Francesco Franco Empirical
More informationTime Series Analysis. Solutions to problems in Chapter 5 IMM
Time Series Analysis Solutions to problems in Chapter 5 IMM Solution 5.1 Question 1. [ ] V [X t ] = V [ǫ t + c(ǫ t 1 + ǫ t + )] = 1 + c 1 σǫ = The variance of {X t } is not limited and therefore {X t }
More informationECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications
ECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications Yongmiao Hong Department of Economics & Department of Statistical Sciences Cornell University Spring 2019 Time and uncertainty
More informationBME 50500: Image and Signal Processing in Biomedicine. Lecture 2: Discrete Fourier Transform CCNY
1 Lucas Parra, CCNY BME 50500: Image and Signal Processing in Biomedicine Lecture 2: Discrete Fourier Transform Lucas C. Parra Biomedical Engineering Department CCNY http://bme.ccny.cuny.edu/faculty/parra/teaching/signal-and-image/
More informationElec4621 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis
Elec461 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis Dr. D. S. Taubman May 3, 011 In this last chapter of your notes, we are interested in the problem of nding the instantaneous
More informationSignals and Systems Lecture (S2) Orthogonal Functions and Fourier Series March 17, 2008
Signals and Systems Lecture (S) Orthogonal Functions and Fourier Series March 17, 008 Today s Topics 1. Analogy between functions of time and vectors. Fourier series Take Away Periodic complex exponentials
More informationLOPE3202: Communication Systems 10/18/2017 2
By Lecturer Ahmed Wael Academic Year 2017-2018 LOPE3202: Communication Systems 10/18/2017 We need tools to build any communication system. Mathematics is our premium tool to do work with signals and systems.
More informationIntroduction to time-frequency analysis Centre for Doctoral Training in Healthcare Innovation
Introduction to time-frequency analysis Centre for Doctoral Training in Healthcare Innovation Dr. Gari D. Clifford, University Lecturer & Director, Centre for Doctoral Training in Healthcare Innovation,
More informationDiscrete Fourier transform
Discrete Fourier transform Alejandro Ribeiro Dept. of Electrical and Systems Engineering University of Pennsylvania aribeiro@seas.upenn.edu http://www.seas.upenn.edu/users/~aribeiro/ January 2, 216 Signal
More informationAutomatic Control III (Reglerteknik III) fall Nonlinear systems, Part 3
Automatic Control III (Reglerteknik III) fall 20 4. Nonlinear systems, Part 3 (Chapter 4) Hans Norlander Systems and Control Department of Information Technology Uppsala University OSCILLATIONS AND DESCRIBING
More informationEE 435. Lecture 30. Data Converters. Spectral Performance
EE 435 Lecture 30 Data Converters Spectral Performance . Review from last lecture. INL Often Not a Good Measure of Linearity Four identical INL with dramatically different linearity X OUT X OUT X REF X
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 information