matlab解析出声音文件的频率

 1 %1、%[X,Fs,Bit]=wavread(‘d:\\ff.wav‘);
 2 %2、然后进行特征提取
 3 %3、然后分类器识别
 4 %4、输出识别结果
 5 %close all;
 6 %clear all;
 7 [X,Fs]=audioread(‘d:\\ff.wav‘);
 8 N=2048;
 9 n=0:N-1;
10 Xk=fft(X,N); %对信号进行快速Fourier变换
11 nUniquePts = ceil((N+1)/2);
12 Xk = Xk(1:nUniquePts); % 选择前半部，因为后半部是前半部的一个镜像
13 magX=abs(Xk);%求得Fourier变换后的振幅
14 %作完傅氏变换后要加绝对值abs(Xk)；dB值：20*log(abs(Xk)/max(abs(Xk)))。
15 ddd=20*log(magX/max(magX))
16 iPoint = -1 ;
17 for i=1:N
18     if magX(i) == max(magX)
19         iPoint = i;
20         i
21         break ;
22     end
23 end
24 % 创建频率数组begin
25 TXBB_THRESHOLD = 20 ;%误差范围
26 TXBB_CHARACTERS_NUM = 64 ;%假设一共64个频率
27 TXBB_BASEFREQUENCY_H = 17000; %起始频率
28 g_fFrequencies =  TXBB_CHARACTERS_NUM*(0:TXBB_CHARACTERS_NUM-1) + TXBB_BASEFREQUENCY_H;
29 % 创建频率数组end
30               
31 % 查找码表 begin
32 k = 1 ;
33 for i=1:TXBB_CHARACTERS_NUM
34     if g_fFrequencies(i)-20 < Fs/N*iPoint && g_fFrequencies(i)+20 > Fs/N*iPoint
35         k = i ;
36         break ;
37     end
38 end 
39 % 查找码表 end
40 
41 %打印结果
42 g_fFrequencies(k)
43 fprintf(‘k:%d.\n‘,k);
44 disp(‘k:.\n‘);