Fourier series - by Ray Betz
Fourier Series
If
- Dirichlet Conditions are satisfied
then we can write
The above equation is called the complex Fourier Series. Given , we may determine by taking the inner product of with . Let us assume a solution for of the form . Now we take the inner product of with over the interval of one period, .
If then,
If then,
We can simplify the above two conclusions into one equation. (What is the delta function below?)
So, we conclude
Orthogonal Functions
The function and are orthogonal on if and only if .
The set of functions are orthonormal if and only if .
Linear Systems
Let us say we have a linear time invarient system, where is the input and is the output. What outputs do we get as we put different inputs into this system? File:Linear System.JPG
If we put in an impulse response, , then we get out . What would happen if we put a time delayed impulse signal, , into the system? The output response would be a time delayed , or , because the system is time invarient. So, no matter when we put in our signal the response would come out the same (just time delayed).
What if we now multiplied our impulse by a coefficient? Since our system is linear, the proportionality property applies. If we put into our system then we should get out .
By the superposition property(because we have a linear system) we may put into the system the integral of with respect to u and we would get out . This is because What would we get if we put into our system? We could find out by plugging in for in the integral that we just found the output for above. If we do a change of variables (, and ) we get . By pulling out of the integral and calling the remaining integral we get .
INPUT | OUTPUT | REASON |
Given | ||
Time Invarient | ||
Proportionality | ||
Superposition | ||
Superposition | ||
Superposition (from above) |
Fourier Series (indepth)
I would like to take a closer look at in the Fourier Series. Hopefully this will provide a better understanding of .
We will seperate x(t) into three parts; where is negative, zero, and positive.
Now, by substituting into the summation where is negative and substituting into the summation where is positive we get:
Recall that
If is real, then . Let us assume that is real.
Recall that Here is further clarification on this property
So, we may write:
In terms of cosine where is an angle.
Fourier Transform
Fourier transforms emerge because we want to be able to make Fourier expransions of non-periodic functions. We can accomplish this by taking the limit of x(t).
Remember that:
Let's substitute in x(t) for k/T substitute f, for 1/T substitute df, and for the summation substitute the integral.
So,
From the above limit we define and .
By using the above transforms we can now change a function from the frequency domain to the time domain or vise versa. We are not limited to just one domain but can use both of them.
We can take the derivitive of and then put it in terms of the reverse fourier transform.
What happens if we just shift the time of ?
In the same way, if we shift the frequency we get:
What would be the Fourier transform of ?
What would happen if we multiplied our time (time scaling) by a constant in ? We will substitute and . If :
Ok, lets take the fourier transform of the fourier series.
Remember:
CD Player
Below is a diagram of how the information on a CD player is read and processed. As you can see the information on the CD is processed by the D/A converter and then sent through a low pass filter and then to the speaker. If you were recording sound, the sound would be captured by a microphone. Then, it should be sent through a low pass filter. The reason you want a low-pass filter is to keep high frequencies (that you don't intend to record) from being recorded. If a high frequency was recorded at say 30 KHz and the maximum frequency you intended to record was 20KHz, then when you played back the recording you would here a tone at 10KHz. From the filter the signal goes onto the A/D converter and then it is ready to be put on the CD. Recording signals (as just described) is essentially the reverse of the operation pictured below.
In Time Domain:
Let's start with a signal , as shown in the below picture. In this signal there is an infinite amount of information. Obviously, we can't hold it all in a computer, but we could take samples every seconds. Lets do that by multiplying by . Since the magnitude of our delta function is one, we get a series of delta functions that record the value of at intervals of . This gives us a result that looks like:
In Frequency Domain:
In the frequency domain we start with . Now we are in frequency, so we must convolve instead of multiply like we did in the time domain. We would have to convolve with .
Aside:
This result looks it could be a fourier series. We would like to get our result in terms of delta functions. As shown below, the periodic delta functions could be represented as a fourier series with coefficients .
Now we can solve for .
Since the only delta function within the integration limits is the delta function at , we can take out the summation and just leave one delta function. Then, evaluating the integral at we get
Now wer are ready to take the convolution.
Time Domain
In order to output as sound any of the signals that we have we must run them through a D/A converter. This is like convolving the below signal by a step function .
This gives us . This is what the signal looks like as it is output through the D/A converter.
Frequency Domain
To find out what we would multiply by in the frequency domain we just take the inverse fourier transform of and we get .
By multiplying . This is hopefully close to what we started with for a signal.
For 2 times oversampling:
In time, multiply: by . This profides points that are interpolated and makes our output sound better because it looks closer to the original wave.
In frequency, convolve: with . The X(f) that you get is great because there is little distortion near the original frequency plot. This means that you can use a cheaper low-pass filter then you would otherwise have been able to.
Nyquist Frequency
If you are sampling at a frequency of 40 KHz, then the highest frequency that you can reproduce is 20 KHz. The nyquist frequency, would be 20 KHz, the highest frequency that can be reproduced for a given sampling rate.
FIR Filters
A finite impulse response filter (FIR filter) is a digital filter that is applied to data before sending it out a D/A converter. This type of filter allows for compensation of the signal before is it destorted so that it will look as it was originally recorded. Using an FIR filter also allows us to put a cheap low-pass filter on after the D/A converter because the signal has been compensated so it doesn't take an expensive low-pass filter, as it would without the FIR filter.
The coefficients that are sent out to the D/A converter are:
where
Example: Design a FIR low-pass filter to pass between and reject the rest.
Our desired response is: , if |f| is less then or equal to or otherwise.
So,
Note: From the Circular Convolution we get:
Discrete Fourier Transforms (DFTs)
The DFT allows us to take a sample of some signal that is not periodic with time and take the Fourier series of it. There is the DFT and the Inverse DFT listed below.
DFT
IDFT
With the DFT all the negative frequency components are just the complex conjugate of the positive frequency components.
One problem with the DFT is that if the sample taken does not begin and end at zero, (or the same point) then we get what is called leakage. Because the DFT is discrete, if the end of the sample is not at the same place it began then it will make a jump back to the point that it began (leakage). This is because the DFT repeats the recorded section of signal over and over. It is this periodic manner of the DFT that allows us to reproduce a discrete signal that is not periodic. The DFT and IDFT are periodic with period N. This can be easily proved by simplifying .
Adaptive FIR Filters
It should be noted that in the above diagram, . The goal of an adaptive FIR filter is to drive the error, e(n), to zero. If we consider that this is a two coefficient filter and we have a contour plot of then we want to travel in the direction of the negative gradient to minimize the error. Let us say that is the stepping size. So...
What would look like?
How might one find an unknown transfer function? Lets use the example of the tuner upper. The idea here is that we want to remove a sine wave from the signal and leave the original signal(voice) in place.