77范文网 - 专业文章范例文档资料分享平台

语音信号采集与处理(2)

来源:网络收集 时间:2019-06-11 下载这篇文档 手机版
说明:文章内容仅供预览,部分内容可能不全,需要完整文档或者需要复制内容,请下载word后使用。下载word有问题请添加微信号:或QQ: 处理(尽可能给您提供完整文档),感谢您的支持与谅解。点击这里给我发消息

程序: fs=22050;

[ x,fs,Nbits] =wavread('C:\\ 2.wav') ; [ x,fs,Nbits] =wavread('C:\\ 2.wav') ; t=0:1/22050:(length(x)-1)/22050;

sound(x,fs) ; %对加载的语音信号进行回放 figure(1)

plot(x) %做原始语音信号的时域图形 title('原始语音信号'); xlabel('time n'); ylabel('fuzhi n'); figure(2)

freqz(x) %绘制原始语音信号的频率响应图 title('频率响应图')

n=length(x); %求出语音信号的长度 y1=fft(x,n) ; %傅里叶变换 y2=fftshift(y1); %对频谱图进行平移 f=0:fs/n:fs*(n-1)/n; %得出频点 figure(3) subplot(2,1,1);

plot(abs(y2)) %做原始语音信号的FFT频谱图 title('原始语音信号FFT频谱') subplot(2,1,2);

plot(f,abs(y2)); %绘制原始语音信号的频谱图 title('原始语音信号频谱') 时域波形和频谱:

2. 给原始的语音信号加上一个高频余弦噪声,频率为(2Khz)。对加噪后的语音进行分析,并画出其信号时域和频谱图。 程序: %加噪声

k=1:1:n; %定义k 值,噪声与原始语音长度一致 c1=10*sin(2*pi*2000*k); %噪声信号的函数 figure(4) plot(k,c1)

title('噪声信号时域图'); %绘制噪声信号的时域波形图 figure(5)

freqz(c1) %绘制噪声信号的频率响应图 title('噪声信号频率响应图')

s=length(c1); %噪声信号的长度

xc=fft(c1,s); %对噪声信号进行频谱变换 xcc=fftshift(xc); f=0:fs/s:fs*(s-1)/s;

figure(6) %绘制噪声语音的FFT频谱图 plot(f,abs(xcc))

title('噪声信号的频谱图'); %噪音与原始信号叠加 z=x'+c1;

y3=fft(z,s); %对叠加信号进行频谱变换 y4=fftshift(y3); %对频谱图平移 f=0:fs/s:fs*(s-1)/s; %得出频点 figure(7)

plot(f,abs(y4)); %绘出加噪语音频谱图 axis([0 8000 0 400]); title('含噪语音信号的频谱'); figure(8) subplot(2,1,1);

plot(f,abs(y2)) %做原始语音信号的FFT频谱图 title('原始语音信号FFT频谱') subplot(2,1,2);

plot(f,abs(y4)) %做加噪语音信号的FFT频谱图 title('加噪语音信号频谱') 时域波形和频谱:

与原始信号对比,区别:

3. 设计合适的滤波器,滤除高频噪声,绘出滤波后的信号频域和时域波形: 程序:

%用双线性变换法设计巴特沃思低通滤波器 wp=0.25*pi; %通带截止频率 ws=0.3*pi; %阻带截止频率

Rp=1; %通带最大衰减 Rs=15; %阻带最小衰减 Fs=22050; Ts=1/Fs;

wp1=2/Ts*tan(wp/2); %将模拟指标转换成数字指标 ws1=2/Ts*tan(ws/2);

[N,Wn]=buttord(wp1,ws1,Rp,Rs,'s'); %选择滤波器的最小阶数 [Z,P,K]=buttap(N); %创建butterworth模拟滤波器 [Bap,Aap]=zp2tf(Z,P,K); [b,a]=lp2lp(Bap,Aap,Wn);

[bz,az]=bilinear(b,a,Fs); %用双线性变换法实现模拟滤波器到数字滤波器的转换

[H,W]=freqz(bz,az); %绘制频率响应曲线 figure(9)

plot(W*Fs/(2*pi),abs(H)) grid

xlabel('频率/Hz') ylabel('频率响应幅度') title('Butterworth') f1=filter(bz,az,z); figure(10) subplot(2,1,1)

plot(t,z) %画出滤波前的时域图 title('滤波前的时域波形'); subplot(2,1,2)

plot(t,f1); %画出滤波后的时域图 title('滤波后的时域波形');

sound(f1,22050); %播放滤波后的信号 p=length(f1); F0=fft(f1,p);

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,免费范文网,提供经典小说综合文库语音信号采集与处理(2)在线全文阅读。

语音信号采集与处理(2).doc 将本文的Word文档下载到电脑,方便复制、编辑、收藏和打印 下载失败或者文档不完整,请联系客服人员解决!
本文链接:https://www.77cn.com.cn/wenku/zonghe/658218.html(转载请注明文章来源)
Copyright © 2008-2022 免费范文网 版权所有
声明 :本网站尊重并保护知识产权,根据《信息网络传播权保护条例》,如果我们转载的作品侵犯了您的权利,请在一个月内通知我们,我们会及时删除。
客服QQ: 邮箱:tiandhx2@hotmail.com
苏ICP备16052595号-18
× 注册会员免费下载(下载后可以自由复制和排版)
注册会员下载
全站内容免费自由复制
注册会员下载
全站内容免费自由复制
注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: