gta5资产网站正在建设,单位做员工招退工在什么网站,市场营销的十大理论,网站建设用哪种语言该楼层疑似违规已被系统折叠 隐藏此楼查看此楼这是我一个学长给的程序#xff0c;他说里面要算得基本上已经实现了改一下就可以用了。。。但是i本人不是主攻matlab的说白了就是不怎么会#xff0c;#xff0c;想问下该怎么改来实现判断清浊音的功能。%对语音信号采样分帧cle…该楼层疑似违规已被系统折叠 隐藏此楼查看此楼这是我一个学长给的程序他说里面要算得基本上已经实现了改一下就可以用了。。。但是i本人不是主攻matlab的说白了就是不怎么会想问下该怎么改来实现判断清浊音的功能。%对语音信号采样分帧clearclcX wavread(E:\003\luyin.wav);%对信号进行预加重xX;x4filter([1,-0.9375],1,x);figure(1)subplot(2,1,1),plot(x)title(原始语音信号);xlabel(样本序列n);ylabel(幅值);subplot(2,1,2),plot(x4)title(原始语音信号的预加重);xlabel(样本序列n);ylabel(幅值);% 幅度归一化到[11]xdouble(x);xx/max(abs(x));% 常数设置FrameLen240; % 帧长取30ms,8kHz的采样率FrameInc80; % 帧移取10ms,1/3amp13;amp22;zcr110;zcr25;maxsilence3; % 3*10ms30msminlen15; % 15*10ms150msstatus0;count0;silence0;% 短时过零率(矢量法)tmp1enframe(x(1:length(x)-1),FrameLen,FrameInc);tmp2enframe(x(2:length(x)),FrameLen,FrameInc);signs(tmp1.*tmp2) 0;diffs(tmp1-tmp2) 0.02;zcrsum(signs.*diffs,2);figure(2)subplot(2,1,1)plot(zcr);title(短时过零率);ylabel(zcr)%计算短时能量ampsum(abs(enframe(filter([1-0.9375],1,x),FrameLen,FrameInc)),2);inzfind(amp1);ammamp(inz);llmin(amm);figure(2)subplot(2,1,2)plot(amp);title(短时能量);ylabel(amp)%调整能量门限amp1ll(max(amp)-ll)/8;amp2ll(max(amp)-ll)/16;%开始端点检测x10;x20;for n1:length(zcr)goto 0;switch statuscase{0,1} % 0静音1可能开始if amp(n) amp1 % 确信进入语音段x1max(n-count-1,1);status2;silence0;countcount1;elseif amp(n) amp2 | zcr(n) zcr2 % 可能处于语音段status1;countcount1;else % 静音状态status0;count0;endcase 2, % 2语音段if amp(n) amp2 | zcr(n) zcr2 % 保持在语音段countcount1;elsesilencesilence1;if silence maxsilence % 静音还不够长尚未结束countcount1;elseif count minlen % 语音长度太短认为是噪声status0;silence0;count0;else % 语音结束status3;endendcase 3, % 3语音结束break;endendcountcount-silence;x2x1count-1;figure(3)subplot(2,1,1)plot(x)title(语音信号的端点检测);axis([1 length(x) -1 1])ylabel(Speech);line([x1*FrameInc x1*FrameInc],[-1,1],color,red);line([x2*FrameInc x2*FrameInc],[-1,1],color,red);%n1(x1*FrameInc-x2*FrameInc)1;yyx(x1*FrameInc:x2*FrameInc);%x1*FrameInc3760,x28320,%yy的长度是4560figure(3)subplot(2,1,2)plot(yy)axis([1 length(yy) -1 1]) %将此处的横坐标改值就可以取不同的语音段现在是整个语音段title(原始语音信号进行端点检测后得到的有用的语音信号段)fs11.025;%设定采样频率yfft(yy);%进行fft变换magabs(y);%求幅值f(0:length(y)-1)*fs/length(y);%进行对应的频率转换figure(4);plot(f,mag);%做频谱图xlabel(频率(Hz));ylabel(幅值);title(信号波幅频谱图);grid;z0.1*rand(1,length(yy));figure(5);plot(z)fs11.025;%设定采样频率Zfft(z);%进行fft变换magabs(Z);%求幅值f(0:length(Z)-1)*fs/length(Z);%进行对应的频率转换figure(6);plot(f,mag);%做频谱图xlabel(频率(Hz));ylabel(幅值);title(噪声波幅频谱图)grid;myyz;figure(7);subplot(2,1,1);plot(m)fs11.025;%设定采样频率Mfft(m);%进行fft变换magabs(M);%求幅值f(0:length(M)-1)*fs/length(M);%进行对应的频率转换figure(8);plot(f,mag);%做频谱图xlabel(频率(Hz));ylabel(幅值);title(混合信号波幅频谱图)grid;%wavwrite(m,s01)
