C+谱图+； 我试图用C++生成正弦信号的频谱图。要生成光谱图，请执行以下操作： 我把真实的正弦信号分成B块 将汉宁窗口应用于每个块（我假设没有重叠）。这应该给我fft的输入 对每个块应用fft，并根据频率系数u[i][0]和u[i][1]计算幅度，并将其放入v[k][i]中，其中k是块数，i是每个窗口的样本 绘制时间tt[k]vsv[k][i]

这给了我以下的情节：

然而，这一情节看起来并不正确

有什么建议可以让它发挥作用吗？这是我的密码：

#include <stdlib.h>
#include <stdio.h>
#include <time.h>
#include <fftw3.h>
#include <iostream>
#include <cmath>
#include <fstream>
using namespace std;
int main()
{
  int i;
  int N=500;//Number of points acquired inside the window
  double Fs=200;//sampling frequency
  int windowsize=100;//Number of samples in each block
  double  T=1/Fs;//sample time 
  double f=50;//frequency(Hz)
  double *in;
  fftw_complex *out;
  double t[N];//time vector 
  double tt[5];
  double  v [5][249];
  fftw_plan plan_forward;
  in = (double*) fftw_malloc(sizeof(double) * N);
  out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * (N/2 + 1));

  for (int i=0; i<= N;i++)
  {
    t[i]=i*T;

    in[i] =0.7 *sin(2*M_PI*f*t[i]);// generate sine waveform
  }

  plan_forward = fftw_plan_dft_r2c_1d (windowsize, in, out, FFTW_ESTIMATE );

  for (int k=0; k< 5;k++){ 
    for (int i = 0; i<windowsize; i++){
      double multiplier = 0.5 * (1 - cos(2*M_PI*i/(windowsize-1)));//Hanning  Window
      in[i] = multiplier * in[i+k*windowsize];
      fftw_execute ( plan_forward );
      v[k][i]=(20*log10(sqrt(out[i][0]*out[i][0]+ out[i][1]*out[i][1])))/N;//The magnitude in dB
    }
  }

  for (int k=0; k< 5;k++){//Center time for each block

    tt[k]=(2*k+1)*T*(windowsize/2);

  }

  fstream myfile;
  myfile.open("example2.txt",fstream::out);

  myfile << "plot '-' using 1:2" << std::endl;

  for (int k=0; k< 5;k++){ 
    myfile << v[k][i]<< " " << tt[k]<< std::endl;
  }
  myfile.close();
  fftw_destroy_plan ( plan_forward );
  fftw_free ( in );
  fftw_free ( out );
  return 0;
}

#包括
#包括
#包括
#包括
#包括
#包括
#包括
使用名称空间std；
int main（）
{
int i；
int N=500；//在窗口内获取的点数
双Fs=200；//采样频率
int windowsize=100；//每个块中的样本数
双T=1/Fs；//采样时间
双f=50；//频率（Hz）
双*英寸；
fftw_复合体*out；
双t[N]；//时间向量
双tt[5]；
双v[5][249]；
fftw_计划向前；
in=（双*）fftw_malloc（sizeof（双）*N）；
out=（fftw_复合体*）fftw_malloc（尺寸）（fftw_复合体）*（N/2+1））；
对于（int i=0；i

#包括
#包括
#包括
#包括
#包括
#包括
#包括
使用名称空间std；
int main（）
{
int i；
int N=500；
双Fs=200；//采样频率
int windowsize=100；//在窗口内获取的点数
双dF=Fs/N；
双T=1/Fs；//采样时间
双f=50；//频率
双*英寸；
fftw_复合体*out；
双t[N]；//时间向量
双ff[N]；
双tt[5]；
双v[5][249]；
fftw_计划向前；
in=（双*）fftw_malloc（sizeof（双）*N）；
out=（fftw_复合体*）fftw_malloc（尺寸）（fftw_复合体）*（N/2+1））；
对于（int i=0；i你什么意思是“看起来不对劲”？我稍后将尝试发布一个更完整的答案，但最明显的问题是，你只需将
tt
（而不是时间和数据）写入文件，然后尝试将其绘制为x vs.y函数（而实际上，二维数据是时间的函数，因此需要三维绘图）。此外，
fftw_execute
应在每个数据块上调用，而不是针对输入缓冲区中的每个样本。这是一个输入错误。我试图查看tt的结果，但忘记将其放回tt vs v，但该绘图用于我的文件。谢谢。我将修复fftw部分。输入错误放在一边，2D绘图仍然不是显示光谱图的正确绘图类型ram，其中所有时间和频率单元（即3D）的强度均以dB为单位。强度通常由颜色映射表示，因此可以在2D平面上看到。
#include <stdlib.h>
#include <stdio.h>
#include <time.h>
#include <fftw3.h>
#include <iostream>
#include <cmath>
#include <fstream>
using namespace std;

int main()
{
int i;
int N=500;
double Fs=200;//sampling frequency
int windowsize=100;//Number of points acquired inside the window
double dF=Fs/N;
double  T=1/Fs;//sample time 
double f=50;//frequency
double *in;
fftw_complex *out;
double t[N];//time vector 
double ff[N];
double tt[5];
double  v [5][249];
fftw_plan plan_forward;
in = (double*) fftw_malloc(sizeof(double) * N);
out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * (N/2 + 1));

for (int i=0; i<= N;i++)
{ t[i]=i*T;
in[i] =0.7 *sin(2*M_PI*f*t[i]);// generate sine waveform
 }

 plan_forward = fftw_plan_dft_r2c_1d (windowsize, in, out, FFTW_ESTIMATE );

for (int k=0; k< 5;k++){ 
  for (int i = 0; i<windowsize; i++){
  double multiplier = 0.5 * (1 - cos(2*M_PI*i/(windowsize-1)));//Hanning Window
 in[i] = multiplier * in[i+k*windowsize];
                                     }
fftw_execute(plan_forward);
for (int j = 0; j <= ((windowsize/2)-1); j++){
v[k][j] = 2*(out[j][0] * out[j][0] + out[j][1] * out[j][1])/(N*N);
v[k][0]   *= 0.5;
 v[k][N/2] *= 0.5;

for (int j = 0; j <= windowsize/2; j++){
   v[k][j] = 10 * log10(v[k][j] + 1e-5);
                                        }
                  }

for (int j = 0; j <= ((windowsize/2)-1); j++)
    {ff[j]=Fs*j/windowsize;
     }

//Center time for each block

for (int k=0; k< 5;k++){

tt[k]=(2*k+1)*T*(windowsize/2);
                        }

double b[6];
fstream myfile;

myfile.open("data.txt",fstream::out);

myfile << "plot '-' using 1:2" << std::endl;

for (int k=0; k< 6;k++) {
                    b[0]=5;
                    b[k+1]=tt[k];
                     myfile <<b[k] << "\t";
                    }
 myfile<< std::endl;

 for (int j = 0; j <= windowsize/2; j++){  myfile << ff[j]<< "\t";

   for (int k=0; k< 5;k++){ myfile << v[k][j]<< "\t";

                          }
 myfile<< std::endl;
                                    }
 myfile.close();

 fftw_destroy_plan ( plan_forward );
 fftw_free ( in );
 fftw_free ( out );
 return 0;
}