数字信号语音处理设计.docx

数字信号设计作业一、设计内容音频信号频谱分析及滤波;1、用计算机 开始—所有程序---娱乐—录音机程序，录取本人的“物电学院”音频信号，时间约为2秒。格式为8KHz采样，8位量化，单声道，以自已名字命名的.wav文件。（格式转化在录音机的“文件”下拉菜单的“属性”，选择“立即转换”,再到“属性”里选择相应参数计算机录音一般是采样率为44.1kHz，16位量化；为减小计算量，在录机的文件—属性—立即转换, 将声音数据转换为采样率8kHz，8位量化）。对语音信号逐字进行频谱分析，分析自己语音信号的频谱特征。用wavread()读取声音文件，作图画出声音的时域波形，对其进行频谱分析，画出其频域波形。分析自已音频信号的特点。3、设计一个0Hz----3.4KHz的IIR低通滤波器， , , 。对“物”字和“电”字音频信号逐字滤波。要求画出所设计滤波器的幅频特性曲线，并用该滤波器对音频信号滤波，画出滤波后的音频信号的时域和频域波形，结合波形比较滤波前后的时域和频域信号。同时用sound回放滤波后的声音信号。（本题即为设计一个IIR低通滤波器，通带截止频率为3400Hz,阻带截止频率为3550Hz，阻带衰减为25dB,通带衰减为1dB）4、设计一个0Hz----3.4KHz的FIR低通滤波器， , , 。对“物”字和“电”字音频信号逐字滤波。要求画出所设计滤波器的幅频特性曲线，并用该滤波器对音频信号滤波，画出滤波后的音频信号的时域和频域波形，结合波形比较滤波前后的时域和频域信号。同时用sound回放滤波后的声音信号。（本题即为设计一个FIR低通滤波器，通带截止频率为3400Hz,阻带截止频率为3550Hz，阻带衰减为25dB,通带衰减为1dB）5、设计一个100Hz----8KHz的IIR带通滤波器，, , , , ，对“学”字和“院”字的音频信号逐字滤波。画出所设计滤波器的幅频特性曲线，并用该滤波器对音频信号滤波，画出滤波后的音频信号的时域和频域波形，结合波形比较滤波前后的时域 和频域信号。同时用sound回放滤波后的声音信号。（本题即为设计一个IIR带通滤波器，通带截止频率为100Hz、8000Hz,阻带截止频率为0Hz、8100Hz，阻带衰减为25dB,通带衰减为1dB）6、设计一个100Hz----8KHz的FIR带通滤波器，, , , 对“学”字和“院”字音频信号逐字滤波。画出所设计滤波器的幅频特性曲线，并用该滤波器对音频信号滤波，画出滤波后的音频信号的时域和频域波形，结合波形比较滤波前后的时域和频域信号。同时用sound回放滤波后的声音信号。 （本题即为设计一个FIR带通滤波器，通带截止频率为100Hz、8000Hz阻带截止频率为0Hz、8100Hz，阻带衰减为25dB,通带衰减为1dB）二、设计步骤及结果1、首先画出由wavread函数采样后的信号的时域波形和它的频域响应波形；语音信号的程序如下:[x1,fs,bits]=wavread(d:\yangli.wav); %读取语音信号的数据，赋给变量x1figure(20)plot(x1) %做原始语音信号以8k采样后的时域图形title(原始语音采样后时域信号);xlabel(时间轴 n);ylabel(幅值 A);figure(21)freqz(x1) %绘制原始语音信号采样后的频率响应图 title(原始语音信号采样后频率响应图);图（一）图(二)2、 对语音信号进行FFT变换，得到FFT频谱特性曲线，与原语音信号的频谱特性曲线进行比较。语音信号FFT的程序如下:[x1,fs,bits]=wavread(d:\yangli.wav); %读取语音信号的数据，赋给变量x1y1=fft(x1,1024); %对采样后信号做1024点FFT变换f=fs*(0:511)/1024;figure(19)subplot(2,1,1);plot(f,abs(y1(1:512)));title(原始语音信号频谱)xlabel(频率/Hz);ylabel(幅值)subplot(2,1,2);plot(abs(y1(1:1024))) %采样后信号的FFT频谱图title(原始语音信号FFT频谱)xlabel(点数N);ylabel(幅值);图（三）3、设计数字滤波器的性能指标：低通滤波器性能指标fp=3400Hz，fc＝3550Hz， Rs＝25db，Rp＝1db。带通滤波器性能指标 fs1＝0Hz，fp1＝100Hz，fs2＝8000Hz，fc2＝8100Hz，Rs＝25db，Rp＝1db。4、用Matlab设计数字滤波器(1)用双线性变换法设计上面要求的3种滤波器。在Matlab中，可以利用函数butter，cheby1和ellip设计IIR滤波器；利用M