matlab2c使用c++实现matlab函数系列教程-fft函数

全栈工程师开发手册 （作者：栾鹏）

matlab2c动态链接库下载
matlab库函数大全
matlab2c基础教程
matlab2c开发全解教程

matlab2c调用方法：

1、下载动态链接库
2、将Matlab2c.dll拷贝到exe同目录下
3、将Matlab2c.h、Matlab2c.lib放到项目头文件目录下
4、在cpp文件中引入下面的代码

#include "Matlab2c.h"#pragma comment(lib,"Matlab2c.lib")  using namespace Matlab2c;

---

matlab中fft函数简介

1、fft函数：

快速傅里叶变换

2、用法说明

xn=[4 3 2 6 7 8 9 0];
Xk=fft(xn) 信号xn可以为实数序列或复数序列

fft的c++源码实现

//实数序列FFT，先将实数序列转化为复数序列，在进行fftCMatrix Matlab2c::fft(Matrix& a){    return fft(CMatrix(a));}//复数序列fft//输入：复数序列a//输出：序列a的fftCMatrix Matlab2c::fft(CMatrix& a){    CMatrix p(a.row,a.column);    int len =a.row*a.column;    double x1,y1;    for(int k=0;k<len;k++)    {        x1=0;y1=0;        for(int i=0;i<len;i++)        {            x1+=a.data[i].real()*std::cos(2*Matlab2c::PI*k*i/len)+a.data[i].imag()*std::sin(2*Matlab2c::PI*k*i/len);            y1+=a.data[i].imag()*std::cos(2*Matlab2c::PI*k*i/len)-std::sin(2*Matlab2c::PI*k*i/len)*a.data[i].real();        }        p.data[k]=Complex(x1,y1);    }    return p;}

fft函数的使用测试

#include "Matlab2c.h"#pragma comment(lib,"Matlab2c.lib")  using namespace Matlab2c;int main(){    double a[]={4,3,2,6,7,8,9,0};    Matrix aa=Matrix(1,8,a);    CMatrix cc=Matlab2c::fft(aa);    cout<<cc.toString()<<endl;    system("pause");    return 0;}