DSP Jtransform

Question

在我的项目中，我有一个文件与脑电图探头(16338探头，256赫兹采样率)。看起来是这样的：

​

​

在我的程序中，我需要从时间域切换到频域。我正在使用Jtransform来完成这个任务，但是我还是坚持了下来。我已经用double数组中的探针加载了我的列，并从Jtransform中执行双工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的文章，并试图从我的教授那里得到一些建议，但他只给了我一些我不太理解的数学公式。我需要绘制频域图，但我仍然不知道如何从那个点得到这些频率(我知道这很愚蠢，因为它离解决方案太近了，但在读了一个星期的DSP之后，我还是找不出)。所以我的问题是：

有人能试着用我最简单的方法来解释它吗?或者用伪代码或一些例子来显示？

Answer 1

所以你从256赫兹采样数据开始，你有16388个垃圾箱，结果应该是每箱0.015621186赫兹。你可以用0.01562赫兹的频率来绘制每个垃圾桶的震级，高达128赫兹(采样率/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图形库。

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

总之，您有数据，只需要绘制它。你也可以显示你在图表上找到的“主导频率”，看看它是什么样子。