FFT
Description
These two Functions will do the 1 dimension Fast Fourier Transform. Data matrix should be of type double. If user have the data matrix in integer form, user should first transform it to double using the member function of matrixbase "CastToDouble".
This function performs 1-dimensional Fast Fourier Transform on each row of data in a matrix. Data matrix should be of type complex. If user have the data matrix in integer form, user should first transform it to complex.
Syntax
int FFT( matrix & matSource, matrix<complex> & matFFTResult, int nColProcess = -1 )
int FFT( matrix<complex> & matSource, matrix<complex> & matFFTResult, int nColProcess = -1 )
int FFT( vector<complex> & vecSource, vector<complex> & vecFFTResult )
Parameters
- matSource
- [input] The data source matrix
- matFFTResult
- [output] The matrix containing the 1 dimension FFT results
- nColProcess
- [input] The number of columns user provided. If it is -1, function will not do any padding or truncating. But if it is larger than the size of the matrix, it will pad zero to meet the size if less than, it will truncate the data matrix.
- matSource
- [input] The data source matrix
- matFFTResult
- [output] The matrix containing the 1 dimension FFT results
- nColProcess
- [input] The number of columns user provided. If it is -1, function will not do any padding or truncating. But if it is larger than the size of the matrix, it will pad zero to meet the size if less than, it will truncate the data matrix.
- vecSource
- [input] The data source vector
- vecFFTResult
- [output] The vector containing the 1-dimensional FFT result
Return
If succeed, it will return 0; Otherwise it returns error indications.
Error:
-1: Provided column number or row number. Must be -1 or positive integer
-2: The data source matrix is empty
-3: When calling member function, it fails.
-4: When calling GetReal or GetImage member function of matrix class, it failed.
If succeed, it will return 0; Otherwise it returns error indications.
Error:
-1: Provided column number or row number. Must be -1 or positive integer
-2: The data source matrix is empty
-3: When calling member function, it fails.
-4: When calling GetReal or GetImage member function of matrix class, it failed.
If succeed, it will return 0; Otherwise it returns error indications.
Error:
-5: source vector is of zero length
n: other error codes returned by matrix-based FFT function
Examples
Prior to compilation, load fft_utils.c to the workspace by executing the following LabTalk command:
Run.LoadOC("Originlab\fft_utils.c", 16);
To retain fft_utils.c in the workspace for successive sessions, drag and drop the file from the Temporary folder to the System folder.
EX1
//This example computes the 1 dimension fast fourier transform on each row of the matrix mRe.
#include <fft_utils.h>
void FFT_ex1()
{
matrix <double> mRe = {{1,2.5,3.5,2,1.0,2.3},{3.2,4.1,2.2,1.0,4,2}};
int ii , jj;
matrix <complex> mResult;
//test 1d fft function for arbitrary source matrix.
int Err = FFT(mRe, mResult, -1);
if ( 0 != Err )
{
ASSERT(false);
printf("one-dimension FFT fails.\n");
}
else
{
for(ii = 0; ii < mResult.GetNumRows(); ii++)
{
for(jj = 0; jj < mResult.GetNumCols(); jj++)
printf("(%f, %fi) ",mResult[ii][jj].m_re,mResult[ii][jj].m_im);
printf("\n");
}
}
}
/**
The output is as following:
(5.021454, 0.000000i) (-0.347011, -0.954594i) (-0.673610, 0.813173i) (-0.530723, 0.000000i) (-0.673610, -0.813173i) (-0.347011, 0.954594i)
(6.736097, 0.000000i) (0.877734, -0.106066i) (-0.796084, -1.378858i) (0.938971, 0.000000i) (-0.796084, 1.378858i) (0.877734, 0.106066i)
**/
EX2
//This example computes the 1 dimension fast fourier transform on each row of the matrix mSource.
#include <fft_utils.h>
void FFT_ex2()
{
matrix <double> mRe = {{1,2.5,3.5,2,1.0,2.3},{3.2,4.1,2.2,1.0,4,2}};
matrix <double> mIm = {{1,0.5,1.0,0.8,3.1,0.3},{2.2,3.2,3.0,1.5,2,2.1}};
int ii , jj;
matrix <complex> mSource, mResult;
mSource.MakeComplex(mRe,mIm);
//test 1d fft function for complex source matrix.
int Err = FFT(mSource, mResult);
if ( 0 != Err )
{
ASSERT(false);
printf("one-dimension FFT fails\n");
}
else
{
for(ii = 0; ii < mResult.GetNumRows(); ii++)
{
for(jj = 0; jj < mResult.GetNumCols(); jj++)
printf("(%f, %fi) ",mResult[ii][jj].m_re,mResult[ii][jj].m_im);
printf("\n");
}
}
}
/**
The output is as following:
(5.021454, 2.735264i) (-1.018762, -1.546554i) (0.139563, 0.547811i) (-0.530723, 1.428869i) (-1.486782, -1.078534i) (0.324740, 0.362634i)
(6.736097, 5.715476i) (1.620196, 0.240945i) (-0.760729, -1.970818i) (0.938971, 0.163299i) (-0.831440, 0.786898i) (0.135272, 0.453077i)
**/
EX3
//This example computes the 1-dimensional Fast Fourier Transform on a complex vector.
#include <fft_utils.h>
void FFT_ex3()
{
vector vR = { 1,2.5,3.5,2,1.0,2.3 };
vector vI = { 1,0.5,1.0,0.8,3.1,0.3 };
vector <complex> vC(vR, vI), vResult(vC.GetSize());
//tedt 1d fft function on data in a complex vector
int Err = FFT(vC, vResult);
if ( 0 != Err )
{
ASSERT(false);
printf("one-dimension FFT fails\n");
}
else
{
for(int ii = 0; ii < vResult.GetSize(); ii++)
{
printf("(%f, %fi) ",vResult[ii].m_re,vResult[ii].m_im);
printf("\n");
}
}
}
/**
The output is as following:
(5.021454, 2.735264i)
(-1.018762, -1.546554i)
(0.139563, 0.547811i)
(-0.530723, 1.428869i)
(-1.486782, -1.078534i)
(0.324740, 0.362634i)
**/
Remark
This function performs 1-dimensional Fast Fourier Transform on data in a complex vector.
Data vector should be of type complex. If user have the data vector in integer or double form, user should first transform it or copy it to a complex vector.
This function in turn calls the matrix-based FFT function.
See Also
header to Included
fft_utils.h
Reference
|