java - Java fftpack的基本使用

Question

我在 MATLAB 中编写了一些简单的音频处理算法，将用于 Android 应用程序。我在将 FFT 实现转换为 Java 时遇到了一些麻烦。

我正在尝试使用fftpack的 Java 版本，据报道它很有效，但对我来说也有点神秘，也许是因为我对 FFT 缺乏掌握。

我只是有一个双精度数组，我想从中生成一个对应的复数数组，代表 FFT。我将对 FFT 进行的唯一计算是取其元素的绝对值。

为了清楚起见，这里是我想使用 Java fftpack 复制的基本 MATLAB 代码：

X = fft(myDoubleArray);
abs(X[i]);

我希望这相当简单，但我无法确定 fftpack 的入口点。

解决了：

LutzL 的以下答案有效，但实际上不需要创建 Complex1D 对象，因为 fftpack 支持双数组作为输入：

RealDoubleFFT rdfft = new RealDoubleFFT(myArray.length);
rdfft.ft(myArray);

score 1 · Accepted Answer

您通过初始化一个新的 fft 对象

fft = new ComplexDoubleFFT(N);

然后你使用 fft 对象的正向变换

Complex1D X;
X.x=new double[N];
X.y=new double[N];
for(k=0;k<N;k++) X.x[k]=x[k];

fft.ft(X);

然后可以访问绝对值

x=X.x[k]; y=X.y[k]; a=Math.sqrt(x*x+y*y);

可以选择除以N标准化 FFT。

添加适当的变量声明。

通常的约定适用，如果以采样频率 fs 对信号进行采样，则条目 X.(xy)[k] 表示 k=N/2 到 N-1 的 k*fs/N 处的频率分量。

您可能想寻找一种在纯 java 中完成的、受更好支持的 FFT 实现，实现对二元长度的 FFT 并不复杂。唯一的复杂之处在于实施转换以获得最后的加速因素。

1 回答 1