Geol 491 Spectral Analysis tom wilsonmail wvu edu
- Slides: 25
Geol 491: Spectral Analysis tom. wilson@mail. wvu. edu
Introduction to Fourier series and Fourier transforms Fourier said that any single valued function could be reproduced as a sum of sines and cosines 5*sin (2 4 t) Amplitude = 5 Frequency = 4 Hz seconds
We are usually dealing with sampled data 5*sin(2 4 t) Amplitude = 5 Frequency = 4 Hz Sampling rate = 256 samples/second Sampling duration = 1 seconds
Faithful reproduction of the signal requires adequate sampling If our sample rate isn’t high enough, then the output frequency will be lower than the input,
The Nyquist Frequency • The Nyquist frequency is equal to one-half of the sampling frequency. • The Nyquist frequency is the highest frequency that can be measured in a signal. Where t is the sample rate Frequencies higher than the Nyquist frequencies will be aliased to lower frequency
The Nyquist Frequency Thus if t = 0. 004 seconds, f. Ny = Where t is the sample rate
Fourier series: a weighted sum of sines and cosines • Periodic functions and signals may be expanded into a series of sine and cosine functions
This applet is fun to play with & educational too. Experiment with http: //www. falstad. com/fourier/
Try making sounds by combining several harmonics (multiples of the fundamental frequency) An octave represents a doubling of the frequency. 220 Hz, 440 Hz and 880 Hz played together produce a “pleasant sound” Frequencies in the ratio of 3: 2 represent a fifth and are also considered pleasant to the ear. 220, 660, 1980 etc. Pythagoras (530 BC)
You can also observe how filtering of a broadband waveform will change audible waveform properties. http: //www. falstad. com/dfilter/
Fourier series • The Fourier series can be expressed more compactly using summation notation You’ve seen from the forgoing example that right angle turns, drops, increases in the value of a function can be simulated using the curvaceous sinusoids.
Fourier series • Try the excel file step 2. xls
The Fourier Transform • A transform takes one function (or signal) in time and turns it into another function (or signal) in frequency • This can be done with continuous functions or discrete functions
The Fourier Transform • The general problem is to find the coefficients: a 0, a 1, b 1, etc. Take the integral of f(t) from 0 to T (where T is 1/f). Note =2 /T What do you get? Looks like an average! We’ll work through this on the board.
Getting the other Fourier coefficients To get the other coefficients consider what happens when you multiply the terms in the series by terms like cos(i t) or sin(i t).
Now integrate f(t) cos(i t) This is just the average of i periods of the cosine
Now integrate f(t) cos(i t) Use the identity If i=2 then the a 1 term =
What does this give us? And what about the other terms in the series?
In general to find the coefficients we do the following and The a’s and b’s are considered the amplitudes of the real and imaginary terms (cosine and sine) defining individual frequency components in a signal
Arbitrary period versus 2 Sometimes you’ll see the Fourier coefficients written as integrals from - to and
Exponential notation cos t is considered Re ei t where
The Fourier Transform • A transform takes one function (or signal) and turns it into another function (or signal) • Continuous Fourier Transform:
The Fourier Transform • A transform takes one function (or signal) and turns it into another function (or signal) • The Discrete Fourier Transform:
Some useful links • • • http: //www. falstad. com/fourier/ – Fourier series java applet http: //www. jhu. edu/~signals/ – Collection of demonstrations about digital signal processing http: //www. ni. com/events/tutorials/campus. htm – FFT tutorial from National Instruments http: //www. cf. ac. uk/psych/Culling. J/dictionary. html – Dictionary of DSP terms http: //jchemed. chem. wisc. edu/JCEWWW/Features/Mcad. In. Chem/mcad 008/ FT 4 Free. Ind. Decay. pdf – Mathcad tutorial for exploring Fourier transforms of free-induction decay http: //lcni. uoregon. edu/fft. ppt – This presentation
Meeting times? Other questions?
- Geology 351
- Settelünke geol
- Vernier spectroscopy
- Multitaper spectral analysis
- Ba491
- External analysis in strategic management
- Presidential decree 491
- Cs 491
- Bahri tokmak
- Color 091807
- Ece 491
- Ece 491
- Cse 491
- Ece 491
- Ece 491
- Ee 491
- Ee 491
- Tomtom go 910 update
- Symbolism in the devil and tom walker
- Edu.sharif.edu
- Spectral regrowth
- Spectral regrowth
- Spectral classification
- Profil spectral rigel
- Spectral normalization
- Spectral hashing