快速傅里叶变换的基2FFT算法的C++实现

最新推荐文章于 2025-01-31 00:15:00 发布

dragon_th

最新推荐文章于 2025-01-31 00:15:00 发布

阅读量2.9k

点赞数

分类专栏：图像处理文章标签： fft 算法 c++ delete 工具 ini

本文链接：https://blog.youkuaiyun.com/dragon_th/article/details/6163247

版权

本文介绍了快速傅里叶变换(FFT)中基2算法的C++实现，详细阐述了相关步骤和代码细节。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

[转自 http://hi.baidu.com/kxw102/blog/item/8672eaeddf54532e279791a4.html]

快速傅里叶变换的基2FFT算法的C++实现

2011-01-19 05:26

快速傅里叶变换的基本原理由于公式不好显示请读者参考其它文章或书籍，本文重点给出了时域抽取法FFT的C++实现。

下面是算法的流程图

倒序的流程图

C++实现代码：

1. FFT.h

#pragma once
#ifndef FFT_H
#define FFT_H
#include <cstring>
#include <cmath>
using namespace std;
#define pi 3.141592
class FFT
{
private:
void changeOrder(double* xr,double* xi,int n);
public:
FFT(void);
void FFT_1D(double* ctxr,double* ctxi,double* cfxr,double* cfxi,int len);
void rFFT_1D(double* ctxr,double* cfxr,double* cfxi,int len);
void IFFT_1D(double* cfxr,double* cfxi,double* ctxr,double* ctxi,int len);
void rIFFT_1D(double* cfxr,double* cfxi,double* ctxr,int len);
~FFT(void);
};
#endif

2.FFT.cpp

#include "./fft.h"
FFT::FFT()
{
}

///
//倒序实现
//xr实部，xi虚部，n为2的幂
///
void FFT::changeOrder(double *xr,double *xi,int n)
{
double temp;
int k;
int lh=n/2;
int j=lh;
int n1=n-2;
for(int i=1;i<=n1;i++)
{
if(i<j)
{
temp=xr[i];
xr[i]=xr[j];
xr[j]=temp;