midfilter.m

来自「语音处理的短时时域处理技术中的短时中值滤波」· M 代码 · 共 22 行

function midfilter()
wavesource=wavread('D:\mymusic\22.wav');%读取语音信号
wl=length(wavesource);
ml=5000;    %设定截取语音信号的长度
nl=21;  %滑动窗口的长度
waves=wavesource(fix(wl/2)-ml/2+1:fix(wl/2)+ml/2,2);%截取语音信号
waveanoise=waves;
noise=max(waveanoise);
waveanoise(500,1)=waveanoise(500,1)+noise;%叠加了噪声
waveanoise(1000,1)=waveanoise(500,1)+noise;%叠加了噪声
[wavend]=midfiltering(waveanoise,nl);%中值滤波
subplot(3,1,1)
plot(waves);%画出语音信号的波形
title('画出语音信号的波形');
subplot(3,1,2)
plot(waveanoise);%画出叠加了噪声的语音信号的波形
title('画出叠加了噪声的语音信号的波形');
subplot(3,1,3)
plot(wavend);%画出中值滤波后的语音信号的波形
title('画出中值滤波后的语音信号的波形');