[PDF] Introduction A Simple Example PDF introFFT.pdf

Introduction In this tutorial, we will discuss how to use the fft (Fast Fourier Transform) command plot the frequency spectrum using the MATLAB fft command

sous Matlab : fft, ifft, fftshift, ifftshift Karsten Plamann, février 2018 1 Définition et syntaxe de fft et ifft1 Y = fft(X) et X = ifft(Y) calculent la transformation de Fourier

[PDF] FFT Tutorial

Matlab's FFT function is an effective tool for computing the discrete Fourier transform of a signal The following code examples will help you to understand the

[PDF] Introduction A Simple Example

Introduction In this tutorial, we will discuss how to use the fft (Fast Fourier Transform) command plot the frequency spectrum using the MATLAB fft command

[PDF] CS425 Lab: Frequency Domain Processing 1 Discrete Fourier

the fast Fourier transform (FFT) is a fast algorithm for computing the discrete MATLAB has three related functions that compute the inverse DFT: 0 ifft 1 ifft2

[PDF] Analyse de Fourier en pratique

Spectre d'un signal sinusoïdal sous Matlab avec l'application SignalAnalyzer Les signaux sont échantillonnés Les spectres sont calculés par FFT : calcul de

[PDF] Fourier series in MATLAB

Find the Fourier coefficients using your MATLAB function: plot the Fourier coefficients vs frequency 5 The FFT Despite the fact that we presented the discrete

[PDF] Notes on FFT-based differentiation - MIT Mathematics

In Matlab, the 1/N normalization is moved from the DFT (Matlab's fft function) to the IDFT (Matlab's ifft function), which doesn't affect any of the procedures in this

[PDF] IMAGE PROCESSING IN FREQUENCY DOMAIN USING MATLAB

15 sept 2008 · FFT: Fast Fourier Transform IFFT: Inverse Fast Fourier Transform MSQE: Mean Square Quantization Error MATLAB: Matrix Laboratory ROI:

[PDF] RECOVERING SIGNAL FREQUENCIES WITH THE - Kenneth L Ho

The DFT is computed in Matlab by the command fft, for the fast Fourier X = FFTFREQ(X,T,DT) returns the discrete Fourier transform of X sampled over the

[PDF] The Fundamentals of FFT-Based Signal Analysis and Measurement

The Fast Fourier Transform (FFT) and the power spectrum are powerful tools for analyzing and measuring signals from plug-in data acquisition (DAQ) devices

[PDF] fiche de lecture a cp

[PDF] fiche de lecture compréhension cp a imprimer

[PDF] fiche de lecture cp a imprimer pdf

[PDF] fiche de lecture cp gratuite a imprimer

[PDF] fiche de lecture cp pdf

[PDF] fiche de lecture cp son a

[PDF] fiche de lecture pour cp a imprimer

[PDF] fiche de poste comptable

[PDF] fiche de poste en anglais

[PDF] fiche de poste exemple word

[PDF] fiche de poste modèle

[PDF] fiche de vocabulaire français par thème pdf

[PDF] fiche lecture ce1 pdf

[PDF] fiche lecture ce2 pdf

[PDF] fiche pedagogique de l'enseignant primaire au maroc

[PDF] Introduction A Simple Example

Introduction

In this tutorial, we will discuss how to use the fft (Fast Fourier Transform) command within MATLAB. The fft command is in itself pretty simple, but takes a little bit of getting used to in order to be used effectively. When we represent a signal within matlab, we usually use two vectors, one for the x data, and one for the y data. The fft command only operates on the y-data (converting the y-data from the time domain into the frequency domain), so it's up to the user to determine what the x-data in the frequency domain will be! This tutorial will show you how to define your x-axis so that your fft results are meaningful. In addition, it will show you how to obtain a two-sided spectrum as well as a positive frequency spectrum for a given signal.

A Simple Example

1. Let's start off with a simple cosine wave, written in the following manner:

2. Next, let's generate this curve within matlab using the following commands:

3. fo = 4; %frequency of the sine wave

4. Fs = 100; %sampling rate

5. Ts = 1/Fs; %sampling time interval

6. t = 0:Ts:1-Ts; %sampling period

7. n = length(t); %number of samples

8. y = 2sin(2pifot); %the sine curve

10. %plot the cosine curve in the time domain

11. sinePlot = figure;

12. plot(t,y)

13. xlabel('time (seconds)')

14. ylabel('y(t)')

15. title('Sample Sine Wave')

16. grid

Here's what we get:

17. When we take the fft of this curve, we would ideally expect to get the following spectrum in the

frequency domain (based on fourier theory, we expect to see one peak of amplitude 1 at -4 Hz, and another peak of amplitude 1 at +4 Hz):

Using Matlab's FFT Command

So now that we know what to expect, let's use MATLAB's built in fft command to try to recreate the frequency spectrum: %plot the frequency spectrum using the MATLAB fft command matlabFFT = figure; %create a new figure YfreqDomain = fft(y); %take the fft of our sin wave, y(t) stem(abs(YfreqDomain)); %use abs command to get the magnitude %similary, we would use angle command to get the phase plot! %we'll discuss phase in another post though! xlabel('Sample Number') ylabel('Amplitude') title('Using the Matlab fft command') grid axis([0,100,0,120]) This doesn't quite look like what we predicted above. If you notice, there are a couple of things that are missing. The x-axis gives us no information on the frequency. How can we tell that the peaks are in the right place?

The amplitude is all the way up to 100

The spectrum is not centered around zero

A Custom Function for fft to Obtain a Two-Sided Spectrum

Here is a helpful function that will simplify the process of plotting a two-sided spectrum. Copy this

code into an m-file and save it. function [X,freq]=centeredFFT(x,Fs) %this is a custom function that helps in plotting the two-sided spectrum %x is the signal that is to be transformed %Fs is the sampling rate

N=length(x);

%this part of the code generates that frequency axis if mod(N,2)==0 k=-N/2:N/2-1; % N even else k=-(N-1)/2:(N-1)/2; % N odd end

T=N/Fs;

freq=k/T; %the frequency axis %takes the fft of the signal, and adjusts the amplitude accordingly

X=fft(x)/N; % normalize the data

X=fftshift(X); %shifts the fft data so that it is centered

This is a relatively simple function to use. The function outputs the correct frequency range and the

transformed signal. It takes in as input the signal to be transformed, and the sampling rate. Let's use the sine wave from above and do a quick example (Remember to set the Matlab directory to the location where you saved the previous m-file). Now, copy and paste these commands into the Matlab command prompt. [YfreqDomain,frequencyRange] = centeredFFT(y,Fs); centeredFFT = figure; %remember to take the abs of YfreqDomain to get the magnitude! stem(frequencyRange,abs(YfreqDomain)); xlabel('Freq (Hz)') ylabel('Amplitude') title('Using the centeredFFT function') grid axis([-6,6,0,1.5])

Here's what you should see:

As you can see, this plot is basically identical to what we would expect! We get peaks at both -4 Hz and +4 Hz, and the amplitude of the peaks are 1.

Redundant Information in the FFT

As you can see from the plots above, the information within the frequency spectrum is entirely symmetric. Thus, we only need one side of the spectrum. In general, the positive side of the spectrum is used, while the negative side is ignored. So let's adjust out function above so that we only get the positive frequencies. A Custom Function for fft to Obtain only the Positive Frequencies The following function is a modification of the above function, and will help you plot only the positive frequencies of the spectrum. function [X,freq]=positiveFFT(x,Fs)

N=length(x); %get the number of points

k=0:N-1; %create a vector from 0 to N-1

T=N/Fs; %get the frequency interval

freq=k/T; %create the frequency range

X=fft(x)/N; % normalize the data

%only want the first half of the FFT, since it is redundant cutOff = ceil(N/2); %take only the first half of the spectrum

X = X(1:cutOff);

freq = freq(1:cutOff); Once again, let's use the same sine wave and put it through this function. Copy and paste the following code into the Matlab command prompt. [YfreqDomain,frequencyRange] = positiveFFT(y,Fs); positiveFFT = figure; stem(frequencyRange,abs(YfreqDomain)); xlabel('Freq (Hz)') ylabel('Amplitude') title('Using the positiveFFT function') grid axis([0,20,0,1.5])

Here's what you should get:

These two functions are very useful, and I still use them all the time!

Download source files and Further Reading

You can download the source files here.

Brush up your Fourier by reading about the theory and background at these links: http://www.dspguide.com/ch8.htm

This is the end of the tutorial.

quotesdbs_dbs2.pdfusesText_4

[PDF] [PDF] Introduction A Simple Example

Introduction

A Simple Example

1. Let's start off with a simple cosine wave, written in the following manner:

2. Next, let's generate this curve within matlab using the following commands:

3. fo = 4; %frequency of the sine wave

4. Fs = 100; %sampling rate

5. Ts = 1/Fs; %sampling time interval

6. t = 0:Ts:1-Ts; %sampling period

7. n = length(t); %number of samples

8. y = 2*sin(2*pi*fo*t); %the sine curve

10. %plot the cosine curve in the time domain

11. sinePlot = figure;

12. plot(t,y)

13. xlabel('time (seconds)')

14. ylabel('y(t)')

15. title('Sample Sine Wave')

16. grid

Here's what we get:

17. When we take the fft of this curve, we would ideally expect to get the following spectrum in the

Using Matlab's FFT Command

The amplitude is all the way up to 100

The spectrum is not centered around zero

N=length(x);

T=N/Fs;

X=fft(x)/N; % normalize the data

Here's what you should see:

Redundant Information in the FFT

N=length(x); %get the number of points

T=N/Fs; %get the frequency interval

X=fft(x)/N; % normalize the data

X = X(1:cutOff);

Here's what you should get:

Download source files and Further Reading

You can download the source files here.

This is the end of the tutorial.

8. y = 2sin(2pifot); %the sine curve