资源描述：
《数字信号课程设计-语音信号的处理与滤波》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、广西科技大学（筹）数字信号处理课程设计设计题目：语音信号的处理与滤波设计过程：1.语音信号的采集利用电脑上的声卡和WINDOWS操作系统可以进行数字信号的采集。可以直接用电脑自带的录音工具进行录音。也可以运用MATLAB中的wavrecord()、wavwrite()、wavplay()、wavread()等函数进行语音的录制、播放、存储、读取。具体采集步骤如下：用MATLAB工具编程录制一段时间为5秒的内容（例如“老师，新年快乐”）的语音，存放在C:UsersPKLDesktoptft.wav,并把存放录音文
2、件的文件夹导入MATLAB的work目录中。录制程序如下：程序：fs=16000;%设计抽样频率channel=1;%设计录音通道t=5;fprintf('按任意键开始%d秒录音:',t);pause;%暂停命令fprintf('正在录音');y=wavrecord(t*fs,fs,channel,'double');%录制语音fprintf('录音结束');wavwrite(y,fs,'C:UsersPKLDesktoptft.wav');fprintf('按任意键回放语音');pause;wavpla
3、y(y,fs);2.语音信号的时频域分析利用MATLAB中的“wavread”命令来读入已经录制好的语音信号，将它赋值给任意向量，再对其进行采样。然后可用MATLAB的绘图函数对时域和频域进行绘图分析。运用下面程序进行采样和绘制时频域图，时域图为附件中的图——1，频域图为附件中的图——2。程序：[y,fs,nbits]=wavread('tft.wav');%加载语音信号wavplay(y,fs);%回放语音信号N=length(y);%语音信号长度Y=fft(y,N);%用傅里叶变换求频域subplot(2,1,1
4、);plot(y);title('图——1原始信号时域波形图');subplot(2,1,2);plot(abs(Y));title('图——2原始信号频谱图');2.采用窗函数法和双线性变换法设计滤波器所需设计的滤波器的性能指标为：（1）低通滤波器性能指标：fb=1000Hz,fc=1200Hz,As=100dB,Ap=1dB。（2）高通滤波器性能指标：fb=5000Hz,fc=4800Hz,As=100dB,Ap=1dB。（3）带通滤波器性能指标：fp1=1200Hz,fp2=3000Hz,fsc1=1000Hz
5、,fsc2=3200Hz,As=100dB,Ap=1dB。(一）用双线性变换法设计滤波器用双线性变换法设计滤波器就是先按低通滤波器的性能指标设计模拟低通滤波器，再转换成相应的过渡模拟滤波器。然后采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。这样可以避免产生频谱混叠现象。在双线性变换法中可用butter,cheby1和ellip设计相应IIR滤波器。设计程序如下：低通滤波器设计程序：fp=1000;ft=5500;%取抽样频率fs=1200;wp=2*pi*fp/ft;%通带频率ws=2*pi*
6、fs/ft;%截止频率Omgp=2*ft*tan(wp/2);Omgs=2*ft*tan(ws/2);[n,Omgc]=buttord(Omgp,Omgs,1,100,'s');[z,p,k]=buttap(n);B=k*real(poly(z));A=real(poly(p));[b,a]=lp2lp(B,A,Omgc);[ba,aa]=bilinear(b,a,ft);[Ha,w]=freqz(ba,aa);plot(w/pi,abs(Ha));title('图——3IIR低通滤波器');附件图——3为IIR低通
7、滤波器频谱图。高通滤波器设计程序：fp=5000;ft=25000;%取抽样频率fs=4800;Rp=1;As=100;wp=2*pi*fp/ft;%通带频率ws=2*pi*fs/ft;%截止频率Omgp=2*ft*tan(wp/2);Omgs=2*ft*tan(ws/2);[n,Omgc]=ellipord(Omgp,Omgs,Rp,As,'s');%计算阶数n和截止频率[z,p,k]=ellipap(n,Rp,As);B=k*real(poly(z));A=real(poly(p));[b,a]=lp2hp(B,
8、A,Omgc);[ba,aa]=bilinear(b,a,ft);[Ha,w]=freqz(ba,aa);plot(w/pi,abs(Ha));title('图——4IIR高通滤波器');附件图——4为IIR高通滤波器频谱图。带通滤波器设计程序：fp1=1200;fp2=3000;fs1=1000;fs2=3200;ft=10000;As=10
