资源简介
我用过了,还可以吧,分享一下,基于谱减法的matlab程序,声音自己去录一个就行了

代码片段和文件信息
clear;
j=sqrt(-1);
%读取一段语音
%[filenamepathname]=uigetfile(‘*.wav‘‘请选择语音文件:‘);
[xfsnbits]=wavread(‘abama.wav‘[1 208210]);
[yfsnbits]=wavread(‘ab.wav‘[1 208210]);
wavin=10*x;%+0.3*y;
wbvin=10*x+0.3*y;
sound(wavinfsnbits);
%[wavinfsnbits]=wavread(z); %提取语音特征信息,wavin保存.wav的抽样值,fs保存该信号的抽样频率
wav_length=length(wavin);
%设置基本信息
%基音周期最大为20ms,为使ifft还原后语音失真尽量小,帧长至少要为基音周期的2倍
%根据fs选择帧长 %定义帧长fram_len和步长step_len 步长为帧长的一半
switch fs
case 8000
frame_len=320;step_len=160;
case 10000
frame_len=400;step_len=200;
case 12000
frame_len=480;step_len=240;
case 16000
frame_len=640;step_len=320;
case 44100
frame_len=1800;step_len=900;
otherwise
frame_len=1800;step_len=900;
end;
n_frame=fix((wav_length-frame_len)/step_len)+1; %确定帧数
%分帧
for i=1:n_frame
n1=(i-1)*step_len+1;
n2=(i-1)*step_len+frame_len;
S(i:)=wavin(n1:n2);
end
win_ham=hamming(frame_len); %定义汉明窗函数
%取前15帧学习噪声
noise_foward15frame=zeros(15frame_len);%定义一个数组放前15帧的噪音短时傅立叶变换结果
for nifrm=1:15
nf=fft(win_ham.*(S(nifrm1:frame_len)‘)frame_len);
noise_foward15frame(nifrm:)=noise_foward15frame(nifrm:)+(abs(nf)‘); %noise_foward15frame保存前15帧的噪音短时傅立叶变换幅度结果
end
am_noise=mean(noise_foward15frame(:1:frame_len));%求出每帧的噪声幅度
%用谱减法进行语音增强
sum_timedomain=zeros(n_frameframe_len);%定义一个数组放时域和的结果
voice_timedomain=zeros(n_frameframe_len);%定义一个数组放去噪后的语音的时域
for ifrm=1:n_frame
sum_timedomain(ifrm1:frame_len)=sum_timedomain(ifrm1:frame_len)+(win_ham‘).*S(ifrm1:frame_len);%sum_timedomain保存整个语音信号的时域和的结果
sf=fft((win_ham‘).*S(ifrm:)frame_len);
phase=angle(sf);%保存这帧语音信号的傅立叶变换的结果的相位
am_signal=abs(sf); %保存这帧语音信号的傅立叶变换的结果的幅度
am_voice=am_signal-am_noise; %谱减 %用信号的幅度减去噪声的幅度得到纯净语音的幅度
for b=1:frame_len
if(am_voice(1b)<0)
am_voice(1b)=0;
end
end
voice=am_voice.*exp(j*phase);%组合相位与幅度得到去噪后的纯净语音信号
sif=real(ifft(voiceframe_len));% 求这帧纯净语音信号的傅立叶反变换的实部
voice_timedomain(ifrm:)=voice_timedomain(ifrm:)+sif;%voice_timedomain保存整个去噪后的语音信号的时域和的结果
end
%求出纯净语音信号的真实幅度
wavout=zeros(1wav_length);
for d=1:n_frame
m1=(d-1)*step_len+1;
m2=(d-1)*step_len+frame_len;
wavout(m1:m2)=wavout(m1:m2)+voice_timedomain(d1:frame_len);
end
%将处理结果输出为‘wav‘文件
wavwrite(wavoutfsnbits[‘去噪后的语音‘]);
sound(wavoutfsnbits);
%将处理前后的结果进行作图比较
subplot(211);
plot(wavin);grid on;
axis([1 wav_length -1 1]);
subplot(212);
plot(wbvin);grid on;
axis([1 wav_length -4 4]);
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 2948 2011-08-11 15:09 新建文件夹 (2)\war5.m
目录 0 2011-08-12 14:59 新建文件夹 (2)
----------- --------- ---------- ----- ----
2948 2
- 上一篇:甘特图-project文件
- 下一篇:相参串脉冲matlab
相关资源
- 高灵敏度GPS接收机MATLAB仿真,附捕获
- 基于MATLAB的质点弹道计算与外弹道优
- 阵列天线的matlab仿真
- MATLAB 经典程序源代码大全
- MATLAB小波软阈值去噪代码33473
- 天线阵的波束形成在MATLAB仿真程序及
- 非线性SVM算法-matlab实现
- 《MATLAB 智能算法超级学习手册》-程序
- 组合导航matlab程序
- 读取txt文件内容matlab代码实现
- Matlab实现基于相关的模板匹配程序
- matlab优化工具箱讲解
- 基于MATLAB的快速傅里叶变换
- 光纤传输中的分布傅立叶算法matlab实
- 基于matlab的图像处理源程序
- matlab 椭圆拟合程序
- 算术编码解码matlab源代码
- optical_flow 光流法 matlab 实现程序
- 引导图像滤波器 Matlab实现
- 分形几何中一些经典图形的Matlab画法
- OFDM系统MATLAB仿真代码
- SVM工具箱(matlab中运行)
- 图像小波变换MatLab源代码
- LU分解的MATLAB实现
- 冈萨雷斯数字图像处理matlab版(第三
- 替代数据法的matlab程序
- 用matlab实现的多站定位系统性能仿真
- 通过不同方法进行粗糙集属性约简m
- k近邻算法matlab实现
- matlab识别系统
评论
共有 条评论