FFTW输入输出显示DSP

FFTW input and output display DSP?

本文关键字:DSP 显示 输入输出 FFTW      更新时间:2023-10-16

我的目的是更多地了解FFTW。我想使用qt/Qwt在图形中显示FFTW输入和输出数据。我觉得我做错了什么,因为我的图表和我想象的不一样。

请让我知道出了什么问题。下面是我的代码:

#include <iostream>
using namespace std;
#include <qwt_plot_curve.h>
#include <qwt_plot.h>
#include <qapplication.h>
#include <cmath>
#include <fftw3.h>
 const int N=256;

int main (int argc, char **argv)
{
QApplication a(argc,argv);
 double Fs=1000;//sampling frequency
 double  T=1/Fs;//sample time 
 double f=500;//frequency
 double t[N-1];//time vector 
 double signal[N-1];
 for (int i=0; i< N-1;i++)
 {
 t[i]=i*T;
 signal[i]=0.7 *sin(2*M_PI*f*t[i]);// generate sine waveform
 }
 fftw_complex  out[N];
 fftw_plan     p3;
 p3 = fftw_plan_dft_r2c_1d(N, signal, out, FFTW_ESTIMATE);//create plan 
 fftw_execute(p3);// FFT
 double reout[N];
 double imgout[N];
for (int i = 0; i < N; i++) {
   reout[i]=out[i][0];
   imgout[i]=out[i][1];
   cout << imgout[i]<<endl;
   // cout << signal[i]<< endl;
  }
    fftw_destroy_plan(p3);

   QwtPlot myPlot;
   QwtPlotCurve *curve =new QwtPlotCurve();
   curve->setSamples(reout,signal,N/2+1);//plot fft 
   curve->attach(&myPlot);
   myPlot.show();
   return a.exec();
  }

您的输入是正弦波,并且您仅绘制FFT结果的实分量,其对应于余弦波相关性。请注意,根据频率的不同,正弦波可以完全正交于余弦基向量,在您的绘图中只留下数值噪声。