jiyin.rar_基音 频率_小波 基音_小波 基音_小波 基音频率_小波基音

x= wavread('声调2.wav'); %去直流分量，如果不去会在短时能量和过零时出错 x=x-mean(x); wavLen=length(x); %幅度归一化到[-1,1] x = double(x); x = x / max(abs(x)); %常数设置 FrameLen = 240; FrameInc =80; tmp1 = voiceframe(x(1:end-1), FrameLen, FrameInc); tmp2 = voiceframe(x(2:end) , FrameLen, FrameInc); framecnt0=size(tmp1); framecnt=framecnt0(1);%记录帧数 %计算短时能量 voice=abs(voiceframe( x, FrameLen, FrameInc)); voice1=voice.*voice; amp = sum(voice1, 2); amp=amp/max(amp);%短时幅度归一化 zcrr=zeros(framecnt,1);%记录过零率 for i=1:framecnt for j=1:FrameLen if(tmp1(i,j)>0&&tmp2(i,j)<0) zcrr(i)=zcrr(i)+1; elseif(tmp1(i,j)<0&&tmp2(i,j)>0)zcrr(i)=zcrr(i)+1; end end end %过零率计算 zcrr=zcrr/framecnt/max(zcrr);%归一化处理 % subplot(311); % plot(x); % subplot(312); % plot(amp); % subplot(313); % plot(zcrr); % % cnt=0;%记录浊音长度 outy=zeros(1,framecnt); Th_zcrr=0.001;%过零阈值线 Th_amp=0.2;%能量阈值线 for i=1:framecnt if((amp(i)>Th_amp)) %if((amp(i)>Th_amp)&&(0.0025<zcrr(i)<0.005)) outy(i)=amp(i)-Th_amp; cnt=cnt+1; end end cnt1=0;%记录浊音段的个数 for i=2:framecnt if((outy(i-1)==0)&&(outy(i)~=0)) cnt1=cnt1+1; end end %记录浊音起始点位置deep0 deep0=zeros(1,cnt1); uu=1; for i=2:framecnt if((outy(i-1)==0)&&(outy(i)~=0)) deep0(uu)=i; uu=uu+1; end end uu=1; %%记录浊音终点位置deep1 deep1=zeros(1,cnt1); for i=1:framecnt-1 if((outy(i+1)==0)&&(outy(i)~=0)) deep1(uu)=i; uu=uu+1; end end disp('浊音起点为：'); disp(deep0); disp('浊音终点为：'); disp(deep1); Secpoint0=deep0*160/8000; Secpoint1=deep1*160/8000; disp('对应时间起点为：'); disp(Secpoint0); disp('对应时间终点为：'); disp(Secpoint1); %基音检测 %test=tmp1(deep0(1)+2,:); %test= double(test); %test = test / max(abs(test)); %用pit记录自相关第一峰值用Fs/pit即得到基音频率 %pit=pitch(test); disp('基音频率为：HZ'); %disp(8000/pit); %用PitLen记录基音频率个数(即浊音总帧数) %用PitFreq记录当前帧的基音频率 deep=[-deep0; deep1]'; PitLen=sum(sum(deep))+cnt1;%每个浊音段实际长多1，故加上浊音段个数 PitFreq=zeros(PitLen,1); Turbid=zeros(PitLen,FrameLen); %将浊音段保存在一个数组中 Atmp=tmp1((deep0(1):deep1(1)),:); for i=2:cnt1 b=tmp1((deep0(i):deep1(i)),:); Atmp=[Atmp;b]; end Turbid=Atmp; %for i=1:deep1(1)-deep0(1)+1 % Turbid(i,:)=tmp1(i+deep0(1)-1,:); %end %for i=deep1(1)-deep0(1)+2:PitLen % Turbid(i,:)=tmp1(i+deep0(2)-(deep1(1)-deep0(1))-2,:); %end %subplot(211); %plot(Turbid(4,:)); %subplot(212); %plot(tmp1(41,:)); %循环计算每帧基音频率 for i=1:PitLen test= double(Turbid(i,:)); test = test / max(abs(test)); PitFreq(i)=8000/pitch(test); end disp(PitFreq); figure(1); plot(PitFreq); % plot(PitFreq(1:112)); % figure(2); % plot(PitFreq(113:224)); % figure(4); % plot(PitFreq(225:336)); % figure(5); % plot(PitFreq(337:448));