matlab解析出声音文件的频率

因为matlab工具集成了很多算法。记录学习了用matlab来解析声音频率。

ff.wav 是一段频率为17640HZ的已知声音文件

贴上代码：

%1、%[X,Fs,Bit]=wavread('d:\\ff.wav');
%2、然后进行特征提取
%3、然后分类器识别
%4、输出识别结果
%close all;
%clear all;
[X,Fs]=audioread('d:\\ff.wav');
N=2048;
n=0:N-1;
Xk=fft(X,N); %对信号进行快速Fourier变换
nUniquePts = ceil((N+1)/2);
Xk = Xk(1:nUniquePts); % 选择前半部，因为后半部是前半部的一个镜像
magX=abs(Xk);%求得Fourier变换后的振幅
%作完傅氏变换后要加绝对值abs(Xk)；dB值：20*log(abs(Xk)/max(abs(Xk)))。
ddd=20*log(magX/max(magX))
iPoint = -1 ;
for i=1:N
    if magX(i) == max(magX)
        iPoint = i;
        i
        break ;
    end
end
% 创建频率数组begin
TXBB_THRESHOLD = 20 ;%误差范围
TXBB_CHARACTERS_NUM = 64 ;%假设一共64个频率
TXBB_BASEFREQUENCY_H = 17000; %起始频率
g_fFrequencies =  TXBB_CHARACTERS_NUM*(0:TXBB_CHARACTERS_NUM-1) + TXBB_BASEFREQUENCY_H;
% 创建频率数组end
              
% 查找码表 begin
k = 1 ;
for i=1:TXBB_CHARACTERS_NUM
    if g_fFrequencies(i)-20 < Fs/N*iPoint && g_fFrequencies(i)+20 > Fs/N*iPoint
        k = i ;
        break ;
    end
end 
% 查找码表 end

%打印结果
g_fFrequencies(k)
fprintf('k:%d.\n',k);
disp('k:.\n');