1.4 Fourier Analysis. Random Data 17
|
|
- Nickolas Morrison
- 5 years ago
- Views:
Transcription
1 Random Data Fourier Analysis In the physical world we sample (and acquire) signals in either the temporal domain or the spatial domain. However, we are frequently interested in thefrequencycontentof theses signals. For example consider a tidal gauge sitting in coastalwaterswhichwe sample every 15 minutes or so. We expect that the mean water surface will rise and fall twice a day (semi-diurnal tide) with a period of just over halfaday. Whatarethe driving forces behind the tides? The primary forces are the gravitational forces of the moon and the sun. However, the moon and the sun have slightly different orbital periods (the moon s is 5 minutes or so longer). Furthermore, the sun is much more massive but much further away and it turns out that the moon accounts for something like 7% of tides. How can we extract information on the frequency content of our signal? One method would be to look at the autocorrelation as this will identifydominantperiodic length scales as the signal will correlate with itself strongly over its natural periods. The dominant frequencies are simply the inverse of these periodic time scales. However, asomewhatmoredirectmethodistoconvertoursignaltothefrequency or spectral domain. Infrequencyspacewewillgetlargeamplitudesatthe dominant frequencies in our signal. Consider our sine wave example back in Section If we represent this signal, a pure sine wave, in the frequency domain we get a single non-zero amplitude at frequency f =1/T while all of the other frequencies have zero amplitude. In fact, if we have a complex periodic signal we will get non-zero amplitudes at the frequency components in the signal where the ratios of these amplitudes give us quantitative information about the relative contribution of the various frequencies to the variance of the signal. O.K, but how do we transform our signal into the frequency domain? Fourier transform. Through the The Fourier transform can be written G(f) = g(x)e i2πxf dx (1.4) where i = 1, f is the frequency, andg(f) isthefourier transform of g(x). The c 215 Edwin A. Cowen
2 18 Random Data inverse Fourier transform is defined as g(x) = G(f)e i2πxf df (1.41) Note the only change is that the argument of the exponent term is of opposite sign. Now, as we will be working with discrete data we will work with the discrete form of the Fourier transform (e.g., we will operate on discrete data notcontinuousdata). The discrete form of the Fourier transform, known as the DFT, is given by G k = N 1 n= g n e i2π(n/n)k k =, 1, 2,...,N 1 (1.42) where k is the index of the sample record in the frequency domain and n is the index of the sample record in the spatial/temporal domain. A few remarks on the DFT: As with the discrete autocorrelation function, the DFT is cyclic - e.g., g n is assumed periodic. G()= N 1 n= g n The resolution in the frequency domain is set by the length of the sample record. I.e., if we are working in the time domain with a signal sampled atf s =1Hz(1 samples per second) and the signal is 1 samples long (T =1seconds)thenour t between samples is.1 s. Now, our f (frequency resolution) is f =1/T and our frequency domain is f k = k/t. From the above it should be clear that the maximum frequency inthefrequency domain is essentially the sample rate, f s =1/ t (more precisely it is f s f/2if we take the discrete frequency of each point to be the average frequency across the interval, hence the lowest resolved frequency is f/2). For real data (which is what we will always be interested in) the DFT is symmetric -hencethemaximumfrequencyatwhichwehaveanyindependentinformationis at f N/2 = f s /2 thisisknownasthenyquist frequency. c 215 Edwin A. Cowen
3 Random Data 19 What happens if the frequency content of the sampled signal is higherthanthe Nyquist frequency based on the sample rate? (e.g., the signal containsfrequency components up to 6 Hz but we only sample at 4 Hz) aliasing! We will discuss this further shortly. The cost of calculating the DFT is O(N 2 )(read ordern 2 ). The exact definition of the DFT and FFT (discussed below) can vary slightly the sign of the terms of the exponent may vary, the 2π term may vary, and the location of the 1/N term (in our case in the inverse DFT) may vary. It is critically important that you verify how the DFT is being implemented on any software packages you choose to use. I have presented the form that is consistent with Matlab (note however in Matlab the indices run from k= 1 to N as Matlab does not allow indices but it is generally more common to think of the lowest index as soihaveretainedthatnotationhere)butihavenotverifiedthis with any other software (e.g., Excel, MathCAD, Python) - if you are using canned software other than Matlab you should! The inverse DFT is g n = 1 N N 1 k= G k e i2π(k/n)n n =, 1, 2,...,N 1 (1.43) The penultimate item above is significant - the cost of a DFT is not insignificant, particularly for large N. Fortunately for us there are algorithms to speed up the calculation. The best known of these is the Fast Fourier Transform (FFT), first published in 1965 (Cooley, J.W.; Tukey, J.W.). The algorithm takes advantage of factorizations in the DFT when N is a power of 2 (e.g., log 2 N is an integer). It is not important that we understand the algorithm precisely (See Bracewell, 1986 for avarietyofexplanations) other than to say that it is an exact calculation of the DFT we have defined. Importantly the cost of the calculation is reduced to O(N log N) -afactoroflogn/n less costly than the DFT. Clearly for large N this savings can become significant quickly! For example with N =124,stillrelativelysmall,thecostoftheFFTislessthan1%ofthecostof the DFT. c 215 Edwin A. Cowen
4 2 Random Data You should also be aware that the restriction to multiples of 2isreallynolongerrelevant. AvarietyofalgorithmsforacceleratedDFTcalculationforvectors of almost arbitrary lengths now exist. In particular a group known as FFTW.ORG ( has assembled an incredibly fast and broad set of these tools such that the DFT of vectors of essentially arbitrary length can be calculated in an accelerated manner. The Matlab FFT functions take advantage of the FFTW library so if you are using Matlab don t worry about your vector lengths! Note due to the historical importance of the FFT algorithm and the fact that computers were once a lot slower a lot of legacy algorithms require vector lengths to be powers of 2. This is easy to upgrade if you are working with such packages. One final note, however, that for multi-dimensional data sets on parallel computer architecture there may still be advantages to working with FFTs that have lengths that are powers of Aliasing To get a feel for the spectral domain lets look at some simple sine waves. If we sample three separate sine waves at 1 Hz (Note: the actual frequencies differ slightly from the integer values shown below to ensure zero side-lobe energy, yet another form of aliasing we will talk about later. The actual frequencies are exact multiples of the reciprocal of the record length, which is 124 samples at 1 Hz or 12.4 s, so multiples of 1/12.4 Hz): a =1sin(2.π2.t); b =2sin(6.π2.t); c =4sin(23.π2.t) for 124 samples (12.4 seconds) and we construct the linear superposition of these three as d = a + b + c we can then look at a plot of the magnitude of their Fourier transforms ( G(k) ) versusf yielding the following: c 215 Edwin A. Cowen
5 Random Data t (s) f (Hz) c 215 Edwin A. Cowen
6 22 Random Data If we sample the same signals for only 128 samples (12.8 seconds) we get t (s) f (Hz) c 215 Edwin A. Cowen
7 Random Data 23 Things to notice in general about the spectra: The spectra of both the 12.4 second and 12.8 second sets (sampled at the same rate) have the same maximum frequency, 1 Hz, the sample frequency. The information from half the sampling rate, known as the Nyquist frequency, on is redundant and is the mirror image of the first half of the frequency spectrum. The resolution of the 12.4 second spectrum is finer ( f =1/(12.4 s)=.98 Hz) than that of the 12.8 second spectrum ( f =1/(128 s) =.781 Hz). The amplitude of the spectral peaks vary in a linear manner with the amplitude of the original signal (note this is only true when sampled such that no energy leaks into the side lobes due to the frequency of the signal (f i )beinganon-integer multiple of the frequency resolution (e.g., f i / f integer),however,ingeneral the area under a peak (which may cross several frequencies) is proportionaltothe amplitude of the temporal component. Now, looking at the 12.4 second spectra, A is N a where N is the number of samples (124). 2 For A vs f the only non-zero components are at 2., the frequency of a(t), and 8. Hz which we can think of as 2. Hz. B is N 2 b. For B vs f the only non-zero components are at 4. Hz and 6. Hz (which we think of as 4. Hz)butwehavenotemporalsignalwiththisfrequency,what happened? Aliasing! A 6 Hz signal sampled at 1 Hz looks like a 4Hzsignal. In fact, if a signal being sampled has frequency components greater than the Nyquist frequency (f Nyq )youwillgetaliasing.ifthesignalbeingsampledhasafrequency component (f i )suchthatf Nyq <f i <f s the peak will simply appear at f i with its mirror appearing at f s f i. c 215 Edwin A. Cowen
8 24 Random Data C is N 2 c. For C vs f the only non-zero components are at 3. Hz and 7. Hz (which we think of as 3. Hz).Thisisaliasingagain-butabitmorecomplexversionthan we just encountered. Our original signal (c(t)) has a single frequency component of 23. Hz. Now, since f i >f s the peak folds back into the spectrum. Thus had f i =11Hzwewouldgetapeakat9Hzandhadf i =12Hzwewouldgetapeak at 8 Hz, etc. This works all the way back to Hz which is equivalent to an aliased 2 Hz signal and then the signal folds again meaning a 21 Hz signal would land at 1Hzandour23Hzsignallandsat3Hz! Thus looking at frequency spectrum like that shown for D can be challenging to interpret unless you know that you have sampled such that the signal contains no frequencies greater than the Nyquist frequency. Hence a general rule for spectral analysis is either set your sample frequency at or greater than twice the maximum frequency in the signal or low-pass filter this signal at the Nyquist frequency. If somebody just handed us these spectra, unless we were also given information about the original signals, we would have no choice but to assume that the spectra are telling us that there is energy at frequencies of 2, 3, and 4Hz. Note aliasing is not limited to the frequency domain. As we have indicated there are two methods for determining the autocorrelation function - circular and non-circular. One approach to the non-circular method assumes that all non-specified values are zero (others use data that may exist outside the bounds of the data record length of the calculation, as shown in section 1.3.4). Thus instead of wrapping the signal back onto itself you simply assume the value is zero. If we take our signals a and b (Eqs 1.2 and 1.29) the non-circular autocorrelation functions are: R aa = [ ] R bb = [ ] (1.44) where I have only showed the positive lags (the system is symmetric about the first component). If we compare the circular and non-circular results you see that the circular c 215 Edwin A. Cowen
9 Random Data 25 results are actually an aliased version of the non-circular results where the terms other than R ξξ () in the circular autocorrelation are the sum of two terms in thenon-circular correlation. For example R aa (2) = in the non-circular calculation while R aa (2) = 1.5 inthecircularcalculationwhichisthesumofr aa (2) and R aa (N 2) in the noncircular calculation. In Matlab you can calculate the circular correlations using cxcorr, whichfor1-dcorrela- tions, i.e., vectors, is an order N 2 calculation cost as you must run through an inner and an outer loop of length N to do all the calculations. As we will see in the next section there is a faster way to do the circular correlation. Matlab calculates the non-circular correlation using xcorr, whichhasfouroptionsthatareworthdiscussingastheyshed light on the non-circular correlation. The options each involve what you divide (normalize) each of the calculated R uv (m) inequation1.19with. Aswrittenthenormalizationis N, thiscorrespondstotheoption biased. Thedefaultoptionis none meaning you do not normalize the coefficients - each is just the sum of the N terms at the indicated lag. The most interesting option is unbiased so let s talk about this a bit. The choice of N for the normalization comes from the fact that we have N values being summed in the calculation of each correlation at a specific lag. However, in the noncircular calculation we only truly have N values at a lag of. At a lag of 1 we will have one of the N required products in the summation being simply a multiplication by zero and hence only N 1actualnon-zerocomponentsmakingupthesum.Similarlywith a lag of 2 there will be two multiplications by zero and hence only N 2actualnon-zero components making up the sum; and so on and so on. Thus by normalizing by N we are including values in the summation that we know will be zero and thisunder-biases the correlation value at non-zero lags. This downward bias gets more significant for larger and larger lags. How do we solve it? By normalizing by the actual number of non-zero components in the sum of each autocorrelation at a specific lag, i.e., N 1atalagof1, N 2atalagof2,andsoon. I.e.,theunbiasednon-circularcorrelation is calculated as: c 215 Edwin A. Cowen
10 26 Random Data N 1 1 R uv (m) = u n v n+m (1.45) N m For the positive lags only. for positive and negative lags this becomes: R uv (m) = 1 N m N 1 (N 1) u n v n+m (1.46) The fourth option for xcorr is coeff which is simply the biased non-circular correlation coefficient, that is, none or biased normalized by their zero-lag values Auto and Cross-Correlation via the FFT As mentioned earlier one of the powers of Fourier analysis is that we can compute correlations more inexpensively in the spectral domain then in the temporal or spatial domain. First for simplicity let us define some nomenclature for correlations. We can write autocorrelation in a short-hand form as: 1 R gg (τ) = lim T T T and similarly we can write the cross-correlation as: 1 R gh (τ) = lim T T T g(t)g(t + τ) dt = g g (1.47) g(t)h(t + τ) dt = g h (1.48) Now we can write down the auto and cross-correlation theorems, respectively, as F [g g] = GG = G(s) 2 (1.49) F [g h] = GH (1.5) where F []means Fouriertransformof and() means complex conjugate. Thus we can find the auto and cross-correlations simply by Fourier transforming the sample records, taking a complex conjugate, performing a multiplication in the spectral domain, and inverse Fourier transforming back. In other words R gg = g g= F 1 [ G(s) 2] (1.51) R gh = g h= F 1 [GH ] (1.52) c 215 Edwin A. Cowen
11 Random Data 27 where F 1 []means inversefouriertransformof. Exercise for the student: choose your favorite analysis package, determine how it is performing FFT s and use the package to determine R aa, R bb,andr ab using the FFT. Note, you can determine the non-circular correlation functions in the spectral domain by padding the original sample records with zeros. How many zeros? If you look back at our discussion of aliasing in the previous section clearly the minimum (and hence most efficient) is N -thelengthoftherecord.henceyoumustdoubletherecordlength. Try this out too! Power Spectrum Power spectrum is a common term used by experimentalists which is more specifically known as a power spectral density (PSD) function which intern statisticians will refer to as either autospectral density functions or cross-spectral density functions. What does this all mean? Let s begin with the definitions S uu (s) = S vv (s) = S uv (s) = R uu (τ)e i2πτs dτ = F [R uu ]= U(s) 2 (1.53) R vv (τ)e i2πτs dτ = F [R vv ]= V (s) 2 (1.54) R uv (τ)e i2πτs dτ = F [R uv ]= U(s)V (s) (1.55) and of course the inverse relationships exist, namely R uu (τ) = R vv (τ) = R uv (τ) = S uu (s)e i2πτs ds = F 1 [S uu ] (1.56) S vv (s)e i2πτs ds = F 1 [S vv ] (1.57) S uv (s)e i2πτs ds = F 1 [S uv ] (1.58) Where S ξξ is our nomenclature for spectral density functions. Forrealfunctions(andour sampled signals will be real data) it can be shown that S uu and S vv are real-valued even c 215 Edwin A. Cowen
12 28 Random Data functions of s while S uv is a complex valued function of s. Sincetheautospectraldensity functions are real valued they are frequently written as one-sided spectra, i.e. G uu (s) = 2 G vv (s) = 2 R uu (τ)e i2πτs dτ (1.59) R vv (τ)e i2πτs dτ (1.6) where the factor of two accounts for the symmetric energy in the negative frequencies. Now, based on our definition of the Fourier transform it should be clear that R ξξ () = S ξξ (s) ds = G ξξ (s) = 1 T T ξ 2 (t) dt (1.61) or in other words the integral of the autospectral density function of a signal is the mean squared value of that signal. Further, if we remove the mean of thesignalfromthe original signal (e.g., ζ = ξ ξ where ξ signifies the mean value of ξ) wefindthat R ζζ () = S ζζ (s) ds = G ζζ (s) ds = 1 T T ζ 2 (t) dt = σ 2 ζ (1.62) or in other words the integral of the autospectral density function of a signal is the variance of that signal. O.K., now that we have the definitions, what do spectral density functions tell us? They tell us the portion of the total signal variance that is accounted for at a specific frequency. If we translate our continuous expressions into the discrete domainwefindthatthe spectral density functions tell us the portion of the total signal variance that is accounted for in a specific frequency band, f. If we consider our three sine waves in Section sampled over 12.4 seconds we can determine the autospectral density functions shown at the top of the next page. How did I calculate these in the discrete domain? It is left to the student to verify that S aa = 1 f s N AA = 1 f s N A2 (1.63) f N 1 k= S aa (k) =σ 2 a (1.64) c 215 Edwin A. Cowen
13 Random Data 29 5 S AA 2 S BB 1 5 S CC 5 S DD f (Hz) Note, these are sine waves so the mean of the signals are zero hence the area under the autospectral density function is the variance and the mean squared value. Now we have a method for determining the frequency content of asignal. Realsignals will be considerably more noisy than the above. For example, if we consider just signal a with the addition of uniform noise (known as white noise) withameanamplitudeof 1 we get the signal and spectra shown on the next page. The final spectrum in the series of plots is the average of the spectrum of 16 sample records. Since the noise is random it has infinite frequency content. Thus by averaging the spectrum of sample records we smooth the noise spectrum towards its true spectral shape (a flat line for white noise) which allows our signal of interest to develop an improved signal-to-noise ratio. Ultimately this is how we determine spectra in real signals where noise and/or smooth spectra are of interest. c 215 Edwin A. Cowen
14 3 Random Data a(t) + noise S EE t (s) <S EE > f (Hz) Summarizing our approach to determining spectra (given an ergodic process): Sample data at a rate twice the rate of the maximum relevant frequency (Nyquist criteria). Sample long enough such that the record can be broken up into N data sets to allow individual spectra to be averaged. Keep in mind that the lengthofthese individual records must be long enough to achieve the desired frequencyresolution. Determine the spectrum of each of the N records and average to get a final estimate of S ξξ. This gives us what we need for now to determine spectra. There remains one additional consideration relevant to our work when determining spectra that we will look at later, side-lobe leakage due to the finite length of the sample records. c 215 Edwin A. Cowen
Before we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
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 information12.1 Fourier Transform of Discretely Sampled Data
494 Chapter 12. Fast Fourier Transform now is: The PSD-per-unit-time converges to finite values at all frequencies except those where h(t) has a discrete sine-wave (or cosine-wave) component of finite
More informationLecture 3: Linear Filters
Lecture 3: Linear Filters Professor Fei Fei Li Stanford Vision Lab 1 What we will learn today? Images as functions Linear systems (filters) Convolution and correlation Discrete Fourier Transform (DFT)
More informationV(t) = Total Power = Calculating the Power Spectral Density (PSD) in IDL. Thomas Ferree, Ph.D. August 23, 1999
Calculating the Power Spectral Density (PSD) in IDL Thomas Ferree, Ph.D. August 23, 1999 This note outlines the calculation of power spectra via the fast Fourier transform (FFT) algorithm. There are several
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 informationLecture 3: Linear Filters
Lecture 3: Linear Filters Professor Fei Fei Li Stanford Vision Lab 1 What we will learn today? Images as functions Linear systems (filters) Convolution and correlation Discrete Fourier Transform (DFT)
More informationFourier Transforms For additional information, see the classic book The Fourier Transform and its Applications by Ronald N. Bracewell (which is on the shelves of most radio astronomers) and the Wikipedia
More informationMath 56 Homework 5 Michael Downs
1. (a) Since f(x) = cos(6x) = ei6x 2 + e i6x 2, due to the orthogonality of each e inx, n Z, the only nonzero (complex) fourier coefficients are ˆf 6 and ˆf 6 and they re both 1 2 (which is also seen from
More informationCommunication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi
Communication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture - 41 Pulse Code Modulation (PCM) So, if you remember we have been talking
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 informationA=randn(500,100); mu=mean(a); sigma_a=std(a); std_a=sigma_a/sqrt(500); [std(mu) mean(std_a)] % compare standard deviation of means % vs standard error
UCSD SIOC 221A: (Gille) 1 Reading: Bendat and Piersol, Ch. 5.2.1 Lecture 10: Recap Last time we looked at the sinc function, windowing, and detrending with an eye to reducing edge effects in our spectra.
More informationN H 2 2 NH 3 and 2 NH 3 N H 2
Chemical Equilibrium Notes (Chapter 18) So far, we ve talked about all chemical reactions as if they go only in one direction. However, as with many things in life, chemical reactions can go both in the
More information7 The Waveform Channel
7 The Waveform Channel The waveform transmitted by the digital demodulator will be corrupted by the channel before it reaches the digital demodulator in the receiver. One important part of the channel
More informationUCSD SIO221c: (Gille) 1
UCSD SIO221c: (Gille) 1 Edge effects in spectra: detrending, windowing and pre-whitening One of the big challenges in computing spectra is to figure out how to deal with edge effects. Edge effects matter
More informationSHOCK AND VIBRATION RESPONSE SPECTRA COURSE Unit 7A. Power Spectral Density Function
SHOCK AND VIBRATION RESPONSE SPECTRA COURSE Unit 7A. Power Spectral Density Function By Tom Irvine Introduction A Fourier transform by itself is a poor format for representing random vibration because
More informationSound & Vibration Magazine March, Fundamentals of the Discrete Fourier Transform
Fundamentals of the Discrete Fourier Transform Mark H. Richardson Hewlett Packard Corporation Santa Clara, California The Fourier transform is a mathematical procedure that was discovered by a French mathematician
More informationElectronic structure, plane waves and pseudopotentials
Electronic structure, plane waves and pseudopotentials P.J. Hasnip Spectroscopy Workshop 2009 We want to be able to predict what electrons and nuclei will do from first principles, without needing to know
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 informationComputational Methods for Astrophysics: Fourier Transforms
Computational Methods for Astrophysics: Fourier Transforms John T. Whelan (filling in for Joshua Faber) April 27, 2011 John T. Whelan April 27, 2011 Fourier Transforms 1/13 Fourier Analysis Outline: Fourier
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 informationEach problem is worth 25 points, and you may solve the problems in any order.
EE 120: Signals & Systems Department of Electrical Engineering and Computer Sciences University of California, Berkeley Midterm Exam #2 April 11, 2016, 2:10-4:00pm Instructions: There are four questions
More informationChapter 5 Simplifying Formulas and Solving Equations
Chapter 5 Simplifying Formulas and Solving Equations Look at the geometry formula for Perimeter of a rectangle P = L W L W. Can this formula be written in a simpler way? If it is true, that we can simplify
More informationLES of Turbulent Flows: Lecture 3
LES of Turbulent Flows: Lecture 3 Dr. Jeremy A. Gibbs Department of Mechanical Engineering University of Utah Fall 2016 1 / 53 Overview 1 Website for those auditing 2 Turbulence Scales 3 Fourier transforms
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 informationTutorial 2 - Learning about the Discrete Fourier Transform
Tutorial - Learning about the Discrete Fourier Transform This tutorial will be about the Discrete Fourier Transform basis, or the DFT basis in short. What is a basis? If we google define basis, we get:
More informationThe discrete and fast Fourier transforms
The discrete and fast Fourier transforms Marcel Oliver Revised April 7, 1 1 Fourier series We begin by recalling the familiar definition of the Fourier series. For a periodic function u: [, π] C, we define
More informationFourier Series. (Com S 477/577 Notes) Yan-Bin Jia. Nov 29, 2016
Fourier Series (Com S 477/577 otes) Yan-Bin Jia ov 9, 016 1 Introduction Many functions in nature are periodic, that is, f(x+τ) = f(x), for some fixed τ, which is called the period of f. Though function
More informationDesigning Information Devices and Systems I Fall 2018 Lecture Notes Note Positioning Sytems: Trilateration and Correlation
EECS 6A Designing Information Devices and Systems I Fall 08 Lecture Notes Note. Positioning Sytems: Trilateration and Correlation In this note, we ll introduce two concepts that are critical in our positioning
More informationCorrelation, discrete Fourier transforms and the power spectral density
Correlation, discrete Fourier transforms and the power spectral density visuals to accompany lectures, notes and m-files by Tak Igusa tigusa@jhu.edu Department of Civil Engineering Johns Hopkins University
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 information13.7 Power Spectrum Estimation by the Maximum Entropy (All Poles) Method
3.7 Maximum Entropy (All Poles) Method 565 LP coefficients for each segment of the data. The output is reconstructed by driving these coefficients with initial conditions consisting of all zeros except
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 informationLecture Wigner-Ville Distributions
Introduction to Time-Frequency Analysis and Wavelet Transforms Prof. Arun K. Tangirala Department of Chemical Engineering Indian Institute of Technology, Madras Lecture - 6.1 Wigner-Ville Distributions
More informationSection 2: Friction, Gravity, and Elastic Forces
Chapter 10, Section 2 Friction, Gravity, & Elastic Forces Section 2: Friction, Gravity, and Elastic Forces What factors determine the strength of the friction force between two surfaces? What factors affect
More informationx(t+ δt) - x(t) = slope δt t+δt
Techniques of Physics Worksheet 2 Classical Vibrations and Waves Introduction You will have encountered many different examples of wave phenomena in your courses and should be familiar with most of the
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 informationDiscrete Fourier Transform
Discrete Fourier Transform Valentina Hubeika, Jan Černocký DCGM FIT BUT Brno, {ihubeika,cernocky}@fit.vutbr.cz Diskrete Fourier transform (DFT) We have just one problem with DFS that needs to be solved.
More informationSPEECH ANALYSIS AND SYNTHESIS
16 Chapter 2 SPEECH ANALYSIS AND SYNTHESIS 2.1 INTRODUCTION: Speech signal analysis is used to characterize the spectral information of an input speech signal. Speech signal analysis [52-53] techniques
More informationAppendix C: Recapitulation of Numerical schemes
Appendix C: Recapitulation of Numerical schemes August 31, 2009) SUMMARY: Certain numerical schemes of general use are regrouped here in order to facilitate implementations of simple models C1 The tridiagonal
More informationTutorial 4. Fast Fourier Transforms Phase factors runs around on the unit circle contains Wave forms
Tutorial 4. Fast Fourier Transforms Phase factors There are functions that produce roots-of-one as a function of time (t) or place (x). A good example is a Bloch function φ(x) = exp(i kx) or the phase
More informationA SuperDARN Tutorial. Kile Baker
A SuperDARN Tutorial Kile Baker Introduction n Multi-pulse Technique and ACFs Why ACFs? Bad and missing lags n Basic Theory of FITACF n Basic Fitting Technique Getting the Power and Width Getting the Velocity
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
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 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 ANALYSIS using Python
(version September 5) FOURIER ANALYSIS using Python his practical introduces the following: Fourier analysis of both periodic and non-periodic signals (Fourier series, Fourier transform, discrete Fourier
More informationSHOCK AND VIBRATION RESPONSE SPECTRA COURSE Unit 6B. Notes on the Fourier Transform Magnitude
SHOCK AND VIBRATION RESPONSE SPECTRA COURSE Unit 6B. Notes on the Fourier Transform Magnitude By Tom Irvine Introduction Fourier transforms, which were introduced in Unit 6A, have a number of potential
More informationLecture No. 1 Introduction to Method of Weighted Residuals. Solve the differential equation L (u) = p(x) in V where L is a differential operator
Lecture No. 1 Introduction to Method of Weighted Residuals Solve the differential equation L (u) = p(x) in V where L is a differential operator with boundary conditions S(u) = g(x) on Γ where S is a differential
More informationDesigning Information Devices and Systems I Fall 2018 Lecture Notes Note Positioning Sytems: Trilateration and Correlation
EECS 6A Designing Information Devices and Systems I Fall 08 Lecture Notes Note. Positioning Sytems: Trilateration and Correlation In this note, we ll introduce two concepts that are critical in our positioning
More informationIntroduction to FT-IR Spectroscopy
Introduction to FT-IR Spectroscopy An FT-IR Spectrometer is an instrument which acquires broadband NIR to FIR spectra. Unlike a dispersive instrument, i.e. grating monochromator or spectrograph, an FT-IR
More informationLab 1: Damped, Driven Harmonic Oscillator
1 Introduction Lab 1: Damped, Driven Harmonic Oscillator The purpose of this experiment is to study the resonant properties of a driven, damped harmonic oscillator. This type of motion is characteristic
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 informationOrder out of Tidal Force Complexity 1 Randall D. Peters, PhD Prof. Emeritus of Physics and 2 Paul Hutchinson, Director
Order out of Tidal Force Complexity 1 Randall D. Peters, PhD Prof. Emeritus of Physics and 2 Paul Hutchinson, Director Mercer University, Macon, Georgia USA TPSO, South Australia 1,2 Member of the Seismological
More informationCS 124 Math Review Section January 29, 2018
CS 124 Math Review Section CS 124 is more math intensive than most of the introductory courses in the department. You re going to need to be able to do two things: 1. Perform some clever calculations to
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 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 informationLab 1: damped, driven harmonic oscillator
Lab 1: damped, driven harmonic oscillator 1 Introduction The purpose of this experiment is to study the resonant properties of a driven, damped harmonic oscillator. This type of motion is characteristic
More informationSignal Modeling, Statistical Inference and Data Mining in Astrophysics
ASTRONOMY 6523 Spring 2013 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Course Approach The philosophy of the course reflects that of the instructor, who takes a dualistic view
More informationFundamentals of the Discrete Fourier Transform
Seminar presentation at the Politecnico di Milano, Como, November 12, 2012 Fundamentals of the Discrete Fourier Transform Michael G. Sideris sideris@ucalgary.ca Department of Geomatics Engineering University
More informationLecture 2: Linear regression
Lecture 2: Linear regression Roger Grosse 1 Introduction Let s ump right in and look at our first machine learning algorithm, linear regression. In regression, we are interested in predicting a scalar-valued
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 information4.3 The Discrete Fourier Transform (DFT) and the Fast Fourier Transform (FFT)
CHAPTER. TIME-FREQUECY AALYSIS: FOURIER TRASFORMS AD WAVELETS.3 The Discrete Fourier Transform (DFT and the Fast Fourier Transform (FFT.3.1 Introduction In this section, we discuss some of the mathematics
More informationOSU Physics Department Comprehensive Examination #115
1 OSU Physics Department Comprehensive Examination #115 Monday, January 7 and Tuesday, January 8, 2013 Winter 2013 Comprehensive Examination PART 1, Monday, January 7, 9:00am General Instructions This
More informationOSE801 Engineering System Identification. Lecture 05: Fourier Analysis
OSE81 Engineering System Identification Lecture 5: Fourier Analysis What we will study in this lecture: A short introduction of Fourier analysis Sampling the data Applications Example 1 Fourier Analysis
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 informationReview of Concepts from Fourier & Filtering Theory. Fourier theory for finite sequences. convolution/filtering of infinite sequences filter cascades
Review of Concepts from Fourier & Filtering Theory precise definition of DWT requires a few basic concepts from Fourier analysis and theory of linear filters will start with discussion/review of: basic
More informationFrequency Resolution Effects on FRF Estimation: Cyclic Averaging vs. Large Block Size
Frequency Resolution Effects on FRF Estimation: Cyclic Averaging vs. Large Block Size Allyn W. Phillips, PhD Andrew. Zucker Randall J. Allemang, PhD Research Assistant Professor Research Assistant Professor
More informationFig. 1: Fourier series Examples
FOURIER SERIES AND ITS SOME APPLICATIONS P. Sathyabama Assistant Professor, Department of Mathematics, Bharath collage of Science and Management, Thanjavur, Tamilnadu Abstract: The Fourier series, the
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More information8/19/16. Fourier Analysis. Fourier analysis: the dial tone phone. Fourier analysis: the dial tone phone
Patrice Koehl Department of Biological Sciences National University of Singapore http://www.cs.ucdavis.edu/~koehl/teaching/bl5229 koehl@cs.ucdavis.edu Fourier analysis: the dial tone phone We use Fourier
More informationSpectral Analysis of Flood Waves in Open Watercourse
International Symposium on Water Management and Hydraulic Engineering Ohrid/Macedonia, 1-5 September 2009 Paper: A65 Spectral Analysis of Flood Waves in Open Watercourse Neven Kuspilić, Kristina Novak,
More information(Refer Slide Time: 01:30)
Networks and Systems Prof V.G K.Murti Department of Electrical Engineering Indian Institute of Technology, Madras Lecture - 11 Fourier Series (5) Continuing our discussion of Fourier series today, we will
More informationLecture - 30 Stationary Processes
Probability and Random Variables Prof. M. Chakraborty Department of Electronics and Electrical Communication Engineering Indian Institute of Technology, Kharagpur Lecture - 30 Stationary Processes So,
More informationFigure 1 Different gravity force of the moon on different places on earth
TIDES How the Sun and Moon cause tides. Spring tides, neap tides, and the three patterns of tides on earth caused by the rotation of Earth and local topographical features: Diurnal, Semi-Diurnal and Mixed
More informationA = 6561 times greater. B. 81 times greater. C. equally strong. D. 1/81 as great. E. (1/81) 2 = 1/6561 as great Pearson Education, Inc.
Q13.1 The mass of the Moon is 1/81 of the mass of the Earth. Compared to the gravitational force that the Earth exerts on the Moon, the gravitational force that the Moon exerts on the Earth is A. 81 2
More informationCompensation. June 29, Arizona State University. Edge Detection from Nonuniform Data via. Convolutional Gridding with Density.
Arizona State University June 29, 2012 Overview gridding with density compensation is a common method for function reconstruction nonuniform Fourier data. It is computationally inexpensive compared with
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Basic Concepts Paul Dawkins Table of Contents Preface... Basic Concepts... 1 Introduction... 1 Definitions... Direction Fields... 8 Final Thoughts...19 007 Paul Dawkins i http://tutorial.math.lamar.edu/terms.aspx
More informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationIn this Lecture. Frequency domain analysis
In this Lecture Frequency domain analysis Introduction In most cases we want to know the frequency content of our signal Why? Most popular analysis in frequency domain is based on work of Joseph Fourier
More informationThe Fourier spectral method (Amath Bretherton)
The Fourier spectral method (Amath 585 - Bretherton) 1 Introduction The Fourier spectral method (or more precisely, pseudospectral method) is a very accurate way to solve BVPs with smooth solutions on
More informationThinking in Frequency
09/05/17 Thinking in Frequency Computational Photography University of Illinois Derek Hoiem Administrative Matlab/linear algebra tutorial tomorrow, planned for 6:30pm Probably 1214 DCL (will send confirmation
More informationCS168: The Modern Algorithmic Toolbox Lecture #11: The Fourier Transform and Convolution
CS168: The Modern Algorithmic Toolbox Lecture #11: The Fourier Transform and Convolution Tim Roughgarden & Gregory Valiant May 8, 2015 1 Intro Thus far, we have seen a number of different approaches to
More informationRoberto s Notes on Linear Algebra Chapter 4: Matrix Algebra Section 7. Inverse matrices
Roberto s Notes on Linear Algebra Chapter 4: Matrix Algebra Section 7 Inverse matrices What you need to know already: How to add and multiply matrices. What elementary matrices are. What you can learn
More informationAction Concept Inventory Answers and explanations
Action Concept Inventory Answers and explanations July 21, 2015 Instructions: Answer all of the questions to the best of your knowledge. The quantity action is referred to frequently in this inventory.
More informationNotes on Fourier Analysis
Notes on Fourier Analysis ASTP 611-01: Statistical Methods for Astrophysics Spring Semester 2017 Contents 0 Preliminaries 1 0.1 Outline........................ 1 1 Fourier Analysis 2 1.1 Continuous Fourier
More informationSPS8. STUDENTS WILL DETERMINE RELATIONSHIPS AMONG FORCE, MASS, AND MOTION.
MOTION & FORCES SPS8. STUDENTS WILL DETERMINE RELATIONSHIPS AMONG FORCE, MASS, AND MOTION. A. CALCULATE VELOCITY AND ACCELERATION. B. APPLY NEWTON S THREE LAWS TO EVERYDAY SITUATIONS BY EXPLAINING THE
More informationReliability Theory of Dynamically Loaded Structures (cont.)
Outline of Reliability Theory of Dynamically Loaded Structures (cont.) Probability Density Function of Local Maxima in a Stationary Gaussian Process. Distribution of Extreme Values. Monte Carlo Simulation
More informationName: Earth 110 Exploration of the Solar System Assignment 1: Celestial Motions and Forces Due on Tuesday, Jan. 19, 2016
Name: Earth 110 Exploration of the Solar System Assignment 1: Celestial Motions and Forces Due on Tuesday, Jan. 19, 2016 Why are celestial motions and forces important? They explain the world around us.
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 informationOn the Tides' paradox.
On the Tides' paradox. T. De ees - thierrydemees @ pandora.be Abstract In this paper I analyse the paradox of tides, that claims that the oon and not the Sun is responsible for them, although the Sun's
More informationThe Discrete Fourier Transform
The Discrete Fourier Transform / 23 The Fourier Transform Any discrete time series exponentials: can be built up of a sum of complex = X m e i2πmn/, n = 0,, m=0 = [ + + + + ]. X 0 X e i2πn (/) X 2 e i2πn
More informationDiscrete Fourier Transform
Last lecture I introduced the idea that any function defined on x 0,..., N 1 could be written a sum of sines and cosines. There are two different reasons why this is useful. The first is a general one,
More informationQuadratic Equations Part I
Quadratic Equations Part I Before proceeding with this section we should note that the topic of solving quadratic equations will be covered in two sections. This is done for the benefit of those viewing
More informationAnalog Integrated Circuit Design Prof. Nagendra Krishnapura Department of Electrical Engineering Indian Institute of Technology, Madras
Analog Integrated Circuit Design Prof. Nagendra Krishnapura Department of Electrical Engineering Indian Institute of Technology, Madras Lecture No - 42 Fully Differential Single Stage Opamp Hello and welcome
More informationLab 4: Kirchhoff migration (Matlab version)
Due Date: Oktober 29th, 2012 TA: Mandy Wong (mandyman@sep.stanford.edu) Lab 4: Kirchhoff migration (Matlab version) Robert U. Terwilliger 1 ABSTRACT In this computer exercise you will modify the Kirchhoff
More informationChapter 11. Output Analysis for a Single Model Prof. Dr. Mesut Güneş Ch. 11 Output Analysis for a Single Model
Chapter Output Analysis for a Single Model. Contents Types of Simulation Stochastic Nature of Output Data Measures of Performance Output Analysis for Terminating Simulations Output Analysis for Steady-state
More informationVibratory Motion -- Conceptual Solutions
Vibratory Motion Vibratory Motion -- Conceptual Solutions 1.) An ideal spring attached to a mass m =.3 kg provides a force equal to -kx, where k = 47.33 nt/m is the spring's spring constant and x denotes
More informationBasics of Proofs. 1 The Basics. 2 Proof Strategies. 2.1 Understand What s Going On
Basics of Proofs The Putnam is a proof based exam and will expect you to write proofs in your solutions Similarly, Math 96 will also require you to write proofs in your homework solutions If you ve seen
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationCalculating Accurate Waveforms for LIGO and LISA Data Analysis
Calculating Accurate Waveforms for LIGO and LISA Data Analysis Lee Lindblom Theoretical Astrophysics, Caltech HEPL-KIPAC Seminar, Stanford 17 November 2009 Results from the Caltech/Cornell Numerical Relativity
More information