matlab声音信号时域频域转换,关于处理用采集卡采集到的声音时域信号转化成频域信号........

求助各位大神，我用麦克风通过采集卡采集到了一段声音数据，数据在附件里，矩阵中第一列为时间，第二列为采集到的声音的电压信号，采样频率为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(&＃39;信号&＃39;);

xlabel(&＃39;信号采集序号&＃39;);

ylabel(&＃39;信号&＃39;);

figure(2);

plot(F,X0)%原始频谱图;

title(&＃39;原始频谱图&＃39;);

xlabel(&＃39;频率/Hz&＃39;);

ylabel(&＃39;幅值&＃39;);

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

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

figure(3);

plot(F,s0)%实际频谱图;

title(&＃39;实际频谱图&＃39;);

xlabel(&＃39;频率/Hz&＃39;);

ylabel(&＃39;幅值&＃39;);

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

figure(4);

plot(F(1:fshift),s0(1:fshift))%由于对称只取一半频谱;

title(&＃39;频谱图&＃39;);

xlabel(&＃39;频率/Hz&＃39;);

ylabel(&＃39;幅值&＃39;);

figure(5);

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

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

plot(F,P); %显示相位图

title(&＃39;相位-频率曲线图&＃39;);

xlabel(&＃39;频率/Hz&＃39;);

ylabel(&＃39;相位/度&＃39;);

得出来的时域图和频域图如下所示：

频域图为何在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

数据文件