求助各位大神，我用麥克風通過采集卡采集到了一段聲音數據，數據在附件里，矩陣中第一列為時間，第二列為采集到的聲音的電壓信號，采樣頻率為1萬Hz，做出時域圖和頻域圖的程序如下所示：

clear;

clc;

load sheng NUM;

fs=10000;

a=NUM(:,1);

b=NUM(:,2);

m=size(b);

N=max(m);

s=fft(b,N-1);

X0=abs(s);

F=(0:N-2)*fs/(N-1);%計算采樣點的對應的頻率;

figure(1);

plot(1:N,b);

title('信號');

xlabel('信號采集序號');

ylabel('信號');

figure(2);

plot(F,X0)%原始頻譜圖;

title('原始頻譜圖');

xlabel('頻率/Hz');

ylabel('幅值');

s0(2:N-1)=X0(2:end)/(0.5*(N-1));

s0(1)=X0(1)/(N-1);

figure(3);

plot(F,s0)%實際頻譜圖;

title('實際頻譜圖');

xlabel('頻率/Hz');

ylabel('幅值');

fshift=round(0.5*(N-1));

figure(4);

plot(F(1:fshift),s0(1:fshift))%由于對稱只取一半頻譜;

title('頻譜圖');

xlabel('頻率/Hz');

ylabel('幅值');

figure(5);

P=phase(s); %計算相位

P=rem(P*180/pi,360); %換算為角度

plot(F,P); %顯示相位圖

title('相位-頻率曲線圖');

xlabel('頻率/Hz');

ylabel('相位/度');

得出來的時域圖和頻域圖如下所示：

頻域圖為何在0處顯示最高呢？求助各位

圖片1.jpg

(20.7 KB, 下載次數: 0)

2015-6-11 10:49 上傳

時域圖

圖片2.jpg

(19.7 KB, 下載次數: 0)

2015-6-11 10:49 上傳

頻域圖

2015-6-11 10:39 上傳

點擊文件名下載附件

215.82 KB, 下載次數: 18

數據文件