DSP Java FFt Jtransform

Question

在我的项目中，我有一个带EEG探针的文件（16338探针，256 Hz采样率）。它看起来像这样：

在我的程序中，我需要从时域切换到频域。我使用Jtransform lib来做到这一点，但我一直坚持下去。 我已经在double数组中加载了带有探针的列，并从Jtransform lib执行双复用FFT。代码如下所示：

int N= chanel.length;
int Fs=256;
double re;
double im;
double fft[]=new double[N*2];
double[] magnitude = new double[N/2];

//System.arraycopy(chanel,0,fft,0,chanel.length);

for (int i=0; i<N-1; i++)
{
  fft[2*i]=chanel[i];
  fft[2*i+1]=0;
}

DoubleFFT_1D fftDo =new DoubleFFT_1D(N);
fftDo.complexForward(fft);

for(int i=0;i<N/2-1;i++)
{
  re=fft[2*i];
  im=fft[2*i+1];

  magnitude[i]=Math.sqrt(re*re+im*im);
}

double max_magnitude =-999999999;
int max_index = -1;
for (int i=0; i<N/2-1;i++)
{
  if (magnitude[i]>max_magnitude)
  {
    max_magnitude=magnitude[i];
    max_index=i;
  }
}

double dominantFreq =  max_index*Fs/N;

这些是该代码不同部分的输出：

我读了很多关于DSP的文章，并试图从我的教授那里得到一些建议，但他只给了我一些我不太了解的数学公式。我需要绘制频域图表，但我仍然不知道如何从那一点获得这些频率（我知道这很愚蠢，因为它非常接近解决方案，但我可以＆＃39经过一周的DSP阅读后，我想出来了。所以我的问题是：

有人可以尝试以最简单的方式解释它或以伪代码或某些示例显示吗？

Answer 1

因此，您开始使用256Hz采样数据，并且您有16388个分档，结果应该是每箱0.015621186 Hz。您可以根据频率计算每个箱子的频率计数0.01562 Hz，最高可达128 Hz（采样率/ 2）。一种方法是将幅度数据放入Excel并将单元格数乘以0.01562以获得频率，然后选择全部并创建图形。

这里我修改了一个循环来输出频率。

double step = Fs/N; // 0.015621186;
double freq = 0;
double freqs[] = double[N/2];
for(int i=0;i<N/2-1;i++)
{
    re=fft[2*i];
    im=fft[2*i+1];

    magnitude[i]=Math.sqrt(re*re+im*im);
    freqs[i] = freq;

    freq = freq + step;
}

否则您可以使用Java图形库。

从概念上讲，这些幅度数中的每一个都是加在一起的箱的频率界限之间的箱中所有正弦波的幅度。

总之，您拥有数据并且只需要绘制图表。您还可以显示在图表上找到的“主导频率”，以查看其外观。