无法让 FFT 在 Octave 中工作

Question

我一直在做一个（可能超级简单的）任务，其中一个步骤是获得傅里叶变换。我按照我大学的指导来转换波声which can be found here（通道：2，采样：17600，采样频率：16KHz）。

看图，好像可行：

[y,fs,wmode,fidx]=readwav('piano.wav','r',-1,0);
left=y(:,1);
amountOfSlices = 6;
samplesPerSlice = fix(length(left) / 6);
frames=enframe(left, samplesPerSlice);
frames=transpose(frames);
fftdata=rfft(frames);
fftdata=fftdata.*conj(fftdata);
plot(fftdata);

接下来，我创建了一个以教程代码为基础的代码文件，并添加了接受参数（分配所需的参数，但为简洁起见省略了）。

samplerate = 512;
% Read the file with raw unscaled audio data from begin to end
[multiData,fs,wmode,fidx]=readwav(filename,'r',-1,0);
disp(sprintf('Number of channels: %d', fidx(5)))
disp(sprintf('Number of samples: %d', fidx(4)))
disp(sprintf('Sample frequency: %d Hz', fs))

% Extract the left channel of the data
leftData = multiData(:, 1);

% Slice the left channel into pieces of the size of 'samplerate'
samplesPerSlice = samplerate

% Splits the leftData vector up into frames of length equal to sample rate
slicedLeftData = enframe(leftData, samplerate)';

% Apply the real data fast fourrier transformation on each data slice
fftdata=rfft(slicedLeftData);
fftdata=fftdata.*conj(fftdata);

plot(fftdata);

你们知道我在这里做错了什么吗？

我的实际问题是： 为什么第二个数据不在 0 到 16.000 Hz 的频域中？我做错了什么？

Answer 1

我猜“这里有什么问题”是指多种颜色？如果是这样，看起来两个图像中的 fftdata plottet 都是矩阵（你想要一个数组，对吗？）。 .* 不是对每个元素都执行操作吗？检查您传入参数的尺寸。

此外，请记住奈奎斯特定理：您永远无法解析任何大于采样频率一半的频率。在采样频率为 16KHz 的情况下，您可能有超过 8kHz 的数据点，但它不会包含任何信息，因此无需将其包含在频域图中。