语音中提取基音频率matlab程序.docVIP

语音中提取基音频率matlab程序 %%corr.m % correlation for pitch estimation % flag =1 from left to right, otherwise, from right to left % r=corr(s,L,flag) % % %flag==1时 % % % %|--i----| %s:..................................................... % . % . % . % . % . % . %s:..................................................... % |_______d________|__________________L-d_____________| %返回值是一个序列 %注意到在做自相关时，错位值从20开始 %考虑到在用这个自相关序列时可能是用于求周期，这个20时有利的 %语音的基音周期一般在20到140，因此小于20没有周期 %而错开20以上可以避免出现不是所求的峰值避免出错 %求到结果后放到rr序列中时，注意到下标为i时实际上已经是错开 %i+20-1,因此pitch.m在调用这个时需要有一句p=19+ind_pitch; function r=corr(s,d,flag) epsilon=1.0e-10; % to avoid divided by zero L=length(s); NM=120; if flag==1 % From left to right for i=20:min(L-d,NM) temp=0.0; temp1=0.0; temp2=0.0; for j=1:d temp=temp+s(j)*s(j+i); temp1=temp1+s(j+i)*s(j+i); temp2=temp2+s(j)*s(j); end pw1=sqrt(temp1)+epsilon; pw2=sqrt(temp2)+epsilon; r(i-19)=temp/(pw1*pw2); end else % From right to left for i=20:min(L-d,NM) temp=0.0; temp1=0.0; temp2=0.0; for j=L-d+1:L temp=temp+s(j-i)*s(j); temp1=temp1+s(j-i)*s(j-i); temp2=temp2+s(j)*s(j); end pw1=sqrt(temp1)+epsilon; pw2=sqrt(temp2)+epsilon; r(i-19)=temp/(pw1*pw2); end end %%maxx.m % Find all the local maxima of a function which are greater than % t times the global maximum % max_index=maxx(s,t) %输入一个序列和一个标量t %首先求出序列s的最大值max_value %然后如果序列中的一个点比相邻两个点都大 %而且还大于max_value的t倍 %就把这个点的位置存放在序列mi中 %在序列mi最后再补上序列的最后一个点 %考虑到第一个点和最后一个点有可能满足要求 %但这两个点都不能跟旁边的比较 %所以还是把这两个点保留了，放在mi序列的第一个和最后一个 %实际上就是以最大值点的t倍为界找出这个界以上的极值点的位置 function mi=maxx(s,t) max_value=max(s); oldstate=0; L=length(s); j=1; mi=[1]; for i=1:L-1 if s(i)s(i+1) newstate=-1; end if s(i)t*max_value