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，先将实数序列转化为复数序列，在进行fft
CMatrix Matlab2c::fft(Matrix& a)
{
	return fft(CMatrix(a));
}

//复数序列fft
//输入：复数序列a
//输出：序列a的fft
CMatrix 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;
}