0

我一直在做一个（可能超级简单的）任务，其中一个步骤是获得傅里叶变换。我按照我大学的指南转换了可以在此处找到的波声（通道：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);

教程 FFT

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

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);

赋值 FFT

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

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

4

1 回答 1

1

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

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

于 2012-10-26T15:09:44.093 回答