Question

我的目的是更多地了解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();

  }

Answer 1

您的输入是正弦波，并且您只绘制FFT结果的实部，它对应于余弦波相关。请注意，根据频率，正弦波可以与余弦基矢量完全正交，在图中只留下数值噪声。

FFTW输入输出显示DSP？

1 个答案: