490 likes | 836 Views
Fourier Analysis. Joseph Fourier (1768-1830). “Yesterday was my 21st birthday, and at that age, Newton and Pascal had already acquired many claims to immortality.” – Fourier in a letter to CL Bonard. What is the goal of Fourier Analysis?. To express a function in terms
E N D
Joseph Fourier (1768-1830) “Yesterday was my 21st birthday, and at that age, Newton and Pascal had already acquired many claims to immortality.” – Fourier in a letter to CL Bonard
What is the goal of Fourier Analysis? To express a function in terms of its component frequencies. Time domain Frequency domain
Why should I care about Fourier analysis? • “fourier”: 29,246 • “spectral analysis”: 27,205 • “filter”: 29,025 • “RNAi”: 7,019 • “drosophila”: 57,523 Pub Med Hits
But what is Fourier analysis useful for? • CALCULATING THINGS: Often, the frequency domain is much simpler than the time domain. • For example: filtering (convolving) becomes trivial in the frequency domain • DESCRIBING THINGS: • The cochlea transforms a time domain signal into a frequency signal. • Many brain regions have oscillations of a particular frequency.
Last week: Vector Spaces • Space of vectors • Each vector could be completely described by • Basis set of vectors • Linear combinations of basis vectors [2 4] [0 1] [1 0]
Today: Frequency Space • Space of functions • Each function can be completely described by: • Basis set of sine & cosine waves • Linear combination of basis waves SINE COSINE
Every function can be completely expressed as a sum of sines & cosines of various amplitudes & frequencies. FOURIER THEOREM:
First let’s review the sine wave: 1/ A A = amplitude f = frequency 1/ = period = 1/frequency = phase (offset) t
A cosine wave is just a sine wave shifted in phase by 90o (=90o). degrees
Every function can be completely expressed as a sum of sines & cosines of various amplitudes & frequencies.
Every function can be completely expressed as a sum of sines & cosines of various amplitudes & frequencies. SQUARE WAVE Mr. Square (waves) Mr. Men Mr. Tickle
Mr. Square (waves) Every function can be completely expressed as a sum of sines & cosines of various amplitudes & frequencies. SQUARE WAVE
Every function can be completely expressed as a sum of sines & cosines of various amplitudes & frequencies. SQUARE WAVE Function Component sine waves
But how did I know which sines and cosines should be summed to create the square wave? Fourier Series These things can be calculated. When I want to calculate them, I will find a handy computer program to calculate them for me.
Fourier Series Back to the square wave example ….
So far, we have talked about how to write a function in terms of sums of sines and cosines. Once you know the component sines and cosines, it’s easy to rewrite the function in the “frequency domain”.
TIME DOMAIN FREQUENCY DOMAIN ω1 F Power (or Energy) ω2 ω1 ω2 Time (sec) Frequency (Hz)
Power spectrum: Graph of the power (square of the amplitude) of the signal at each frequency FREQUENCY DOMAIN ω1 Power (or Energy) ω2 ω1 ω2 Frequency (Hz)
0 1 2 3 4 5 FREQUENCY 0 1 2 3 4 5 FREQUENCY SQUARE WAVE Back to the square wave example ….
Using complex exponentials instead of sines and cosines Complex exponentials are sums of sines and cosines. Euler equation: Geometric interpretation: ANY function can be written as a sum of complex exponentials rather than a sum of sines and cosines.
Fourier Transform The extension of the Fourier series for non-periodic functions freq domain time domain: time domain freq domain: DON’T WORRY: these formulas will always be calculated on the computer. Your job is to understand how to interpret the output.
Fourier Transform Imaginary F(ω)= a(ω)+ib(ω) Amplitude or Magnitude: b(ω) The “power” or “energy” is the amplitude squared. a(ω) Real Phase:
A function can be described in the frequency domain in terms of:1) The amplitude of sine and cosinesof various frequenciesOR2) The amplitude and phase of complex exponentials of various frequencies http://falstad.com/fourier/
Any signal can be expressed in the time domain OR the frequency domain • You can go back and forth between these 2 representations. • The fourier transform goes from the time domain to the frequency domain. • The inverse fourier transform goes the other way. • Now, some neuroscience applications …
Parkinson Tremors Time domain Frequencydomain Moore GP, Ding L, Bronte-Stewart HM. Concurrent Parkinson tremors. J Physiol. 2000 Nov 15;529 Pt 1:273-81.
Fourier Transforms in the nervous system • Hair cells vibrate in response to sound and are “tuned” for different frequencies
Tonotopy in the cochlea Sound: change in air pressure over time The filtering properties of the cochlea decompose the incoming signal by frequency, thereby taking a FOURIER TRANSFORM. Low frequency tuning High frequency tuning
Spectrogram =A power spectrum taken at subsequent time bins. frequency “Hey” + music time A good approximation of firing rates in a tonotopic nucleus. http://www-users.cs.york.ac.uk/~alistair/research/dphil/enm/asa/complex.html
Spatial Frequency changes in luminance across space
2-D Power Spectra Spaital grating y x x x F F F Power spectrum ωy ωx ωx ωx
Retinal ganglion cells are tuned to specific spatial frequencies What gives them their spatial frequency tuning?
What we’ve covered so far • Frequency lingo • Amplitude, frequency, period, phase • Fourier’s Theorem • Power spectra • Next week: putting it all together to make predictions about data
Noise (terminology) “white noise” “colored noise” TIME FREQUENCY
Acoustics Pitch fork Trumpet Bassoon Sound: change in air pressure over time
Fourier Imaging Correlation Spectroscopy Margineantu D, Capaldi RA, Marcus AH. Dynamics of the mitochondrial reticulum in live cells using Fourier imaging correlation spectroscopy and digital video microscopy. Biophys J. 2000 Oct;79(4):1833-49.
Why can you create any function by combining sines and cosines? Sines & cosines are orthogonal They form a basis set that spans all functions. ODD: f(x)=-f(x) EVEN: f(x)=f(-x)
Tonotopy in the cochlea AN UNROLLED COCHLEA: Small & stiff: high frequency resonance Wide & floppy: low frequency resonance Tonotopy is maintained through many layers of auditory processing
Power spectrum: Graph of the amount of power (square of the amplitude) at each frequency FREQ DOMAIN TIME DOMAIN