数字信号语音信号的处理与滤波.doc

广西科技大学（筹） 数字信号处理课程设计 设计题目：语音信号的处理与滤波 设计过程： 1.语音信号的采集 利用电脑上的声卡和WINDOWS操作系统可以进行数字信号的采集。可以直接用电脑自带的录音工具进行录音。也可以运用MATLAB中的wavrecord()、wavwrite()、wavplay()、wavread()等函数进行语音的录制、播放、存储、读取。具体采集步骤如下： 用MATLAB工具编程录制一段时间为5秒的内容（例如“老师，新年快乐” ）的语音，存放在C:\Users\PKL\Desktoptft.wav,并把存放录音文件的文件夹导入MATLAB的work目录中。录制程序如下： 程序： fs=16000;%设计抽样频率 channel=1;%设计录音通道 t=5; fprintf(按任意键开始%d秒录音:,t); pause;%暂停命令 fprintf(正在录音); y=wavrecord(t*fs,fs,channel,double);%录制语音 fprintf(录音结束\n); wavwrite(y,fs,C:\Users\PKL\Desktop\tft.wav); fprintf(按任意键回放语音); pause; wavplay(y,fs); 语音信号的时频域分析 利用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); plot(y); title(图——1 原始信号时域波形图); subplot(2,1,2); plot(abs(Y)); title(图——2 原始信号频谱图); 采用窗函数法和双线性变换法设计滤波器 所需设计的滤波器的性能指标为： 低通滤波器性能指标：fb=1000Hz,fc=1200Hz,As=100dB,Ap=1dB。 高通滤波器性能指标：fb=5000Hz,fc=4800Hz,As=100dB,Ap=1dB。 带通滤波器性能指标：fp1=1200Hz,fp2=3000Hz,fsc1=1000Hz,fsc2=3200Hz,As=100dB,Ap=1dB。 用双线性变换法设计滤波器 用双线性变换法设计滤波器就是先按低通滤波器的性能指标设计模拟低通滤波器，再转换成相应的过渡模拟滤波器。然后采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。这样可以避免产生频谱混叠现象。在双线性变换法中可用butter,cheby1和ellip设计相应IIR滤波器。设计程序如下： 低通滤波器设计程序： fp=1000; ft=5500;%取抽样频率 fs=1200; wp=2*pi*fp/ft;%通带频率 =2*pi*fs/ft;%截止频率 [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(图——3 IIR低通滤波器); fp=5000; ft=25000;%取抽样频率 fs=4800; Rp=1;As=100; wp=2*pi*fp/ft;%通带频率 =2*pi*fs/ft;%截止频率 [n,Omgc]=ellipord(