.WAV 文件上的 FFT

Question

我正在尝试对来自 .wav 文件的信号的特定部分进行快速傅立叶变换。

.wav 文件是每 0.6 秒重复一次信号。

我试图弄清楚如何重复 10 次并进行比较以表明结果应该相似。

这是我现在拥有的代码：

Fs = 44100;
cj = sqrt(-1);

[test,fs]= wavread('3b healthy2.wav'); % File data name


dt = 1/Fs;
time = 45.6;
N = time/dt;

left=test(:,1);
right=test(:,2);


I = left;
Q = right;

t = 0:dt:(time-dt);

n = length(t);

f = -Fs/2:Fs/n:Fs/2-Fs/n;

s = I+cj.*Q;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Smooth the signal
ss = smooth(s,201);

sf = fftshift(fft(ss(1:N))); % taking fft
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

figure(2)

plot(f,(abs(sf))./max(abs(sf)))

因此，一旦我分离程序，我就可以找出时域，并从那里发现一个重复是在 45.2 到 45.8。

在我运行这个程序后，我得到了一个图表，但据我所见，它不包括一次重复，而是所有重复加起来达到 45.2

编辑

下图是信号在 45 秒时的样子。信号连续重复 3 分钟。现在我需要 45 秒重复的 FFT

Answer 1

您似乎想为要分析的信号设置一个窗口。

如果是这样，在时域中你只需截断你的向量。

left=test(:,1);
right=test(:,2);
time = 45;
interval =1; // per second, for 0.1 second use interval = 0.1;
w_range = time*Fs: (time+interval)*Fs-1;
I = left(w_range) // truncate here 
Q = right(w_range) // and here
n = interval * Fs;
f = -Fs/2:Fs/n:Fs/2-Fs/n;
// continue FFT analysis ...

有意义吗？如果这不是您想要的，请发表评论。