河南理工大学数字信号处理课程设计报告

数字信号处理课程设计报告题目：语音信号去噪处理专业班级XXXXXXXXX学号XXXXXXXXXXXXX学生姓名XXXXXX指导教师张培玲教师评分2016年1月13日数字信号处理课程设计2语音信号滤波去噪——使用脉冲响应不变法设计的巴特沃斯滤波摘要本课程设计主要运用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出合理的结论。关键词课程设计；滤波去噪；巴特沃斯滤波器；脉冲响应不变法；MATLAB数字信号处理课程设计3目录一、设计目的与要求………………………………………..4二、总体设计方案…………….……………………………52.1实验步骤………………………………………………52.2设计内容…………………………………………….62.3语音信号的频谱分析………………………………72.4滤波器设计……………………….…………………82.5结果分析………………………….…………………13三、设计原理结果与仿真分析…………………………133.1设计原理………………………….…………………133.2仿真结果与分析……………………………………13四、结论……………………………………………………14五、心得体会………………………………………………14参考文献……………………………………………………14附录………………………………………………………….15数字信号处理课程设计4一、设计目的与要求1.1课程设计目的《数字信号处理》课程设计是在学生完成数字信号处理和MATLAB的结合后的基本实验以后开设的。本课程设计的目的是为了让学生综合数字信号处理和MATLAB并实现一个较为完整的小型滤波系统。这一点与验证性的基本实验有本质性的区别。开设课程设计环节的主要目的是通过系统设计、软件仿真、程序安排与调试、写实习报告等步骤，使学生初步掌握工程设计的具体步骤和方法，提高分析问题和解决问题的能力，提高实际应用水平。1.2课程设计的要求（1）滤波器指标必须符合工程设计。（2）设计完后应检查其频率响应曲线是否满足指标。（3）处理结果和分析结论应该一致，而且应符合理论。（4）独立完成课程设计并按要求编写课程设计报告。数字信号处理课程设计5二、总体设计方案2.1设计步骤语音信号滤波去噪——使用脉冲不变响应法设计的巴特沃斯滤波器的设计流程如图3.1.1所示：图2.1.1脉冲响应不变法巴特沃斯滤波器对语音信号去噪流程图开始语音信号的采集（wavread函数），画时域图快速傅里叶变换，并且画频谱图设定滤波器性能指标，通带截止频率fb=1100，阻带截止频率fc=1200，通带波纹Ap=1，阻带波纹As=20脉冲响应不变法法设计巴特沃斯滤波器验证并进行频谱分析设计好的滤波器进行滤波处理比较滤波前后语音信号的波形及频谱回放语音信号结束数字信号处理课程设计62.2设计内容用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析点击windows系统桌面的“开始”按钮，点击开始菜单栏里的“附件”，选择“录音机”选项，点击录音机“文件”选项，进入“声音选定”设置，把属性一栏设置成“8000Hz，8位，单声道，7KB/秒”(见图3.2.1)。点击确定，然后开始语言信号的采集，采集时间为1秒左右为最佳。采集的声音文件以“.wav”格式存储（见图3.2.2）。图2.2.1采集声音的参数设置图2.2.2采集声音数字信号处理课程设计72.3语音信号的频谱分析在MATLAB中编辑m函数，使用wavread函数读取采集的声音文件（.wav)将它赋值给某一向量，再对其进行采样,然后使用plot语句画出相关的频谱图形在figure(1）上。（1）Wavread函数调用格式：[y,Fs,nbits]=wavread(file)功能说明：采样值放在向量y中，Fs表示采样频率（Hz），nbits表示采样位数。（2）快速傅里叶变换算法FFT计算DFT的函数fft,其调用格式如下：Xk=fft(x,n)参数x为被变换的时域序列向量，N是DFT变换区间长度，当n大于x的长度时，fft函数自动在x后面补零。，当n小于xn的长度时，fft函数计算x的前n个元素，忽略其后面的元素。在本次课程设计中，我们利用fft函数对语音信号进行快速傅里叶变换，就可以得到信号的频谱特性。（3）声音采样文件读取的程序（文件名：c.wav)%用麦克风采集一段8000Hz，8k的单声道语音信号，绘制波形并观察其频谱[x,Fs,nbits]=wavread('c.wav')figure(1)n=length(x)x1=fft(x,n)x2=abs(x1)x3=angle(x1)subplot(4,1,1);stem(x);title('原信号时域')subplot(4,1,2);stem(x1);title('原信号fft')subplot(4,1,3);stem(x2);title('fft幅值')subplot(4,1,4);stem(x3);title('fft相角')数字信号处理课程设计8所得语音频谱分析图见图2.3.1图2.3.1声音信号的频谱分析（第一幅是时域波形；采样1是FFT后波形；采样2是频域幅度谱；采样3是频域相位谱）2.4滤波器设计设计指标：通带截止频率为1100Hz，阻带截止频率为1200Hz，通带波纹为1dB，阻带波纹为20dB，用脉冲（冲激）响应不变法设计的一个满足上述指标的巴特沃思IIR滤波器%巴特沃思低通滤波器Fs=8000;%采样频率wp=1100*2/Fs;ws=1200*2/Fs;%根据采样频率将滤波器边界进行转化Rp=1;Rs=20;%通带波纹和阻带波纹Nn=128;[N,Wn]=buttord(wp,ws,Rp,Rs);%求滤波器的最小阶数和3dB固有频率[b,a]=butter(N,Wn)%设计BUTTERWORTH低通滤波器，返回滤波器的系数矩阵数字信号处理课程设计9figure(2);[H,f]=freqz(b,a,Nn,Fs);%用Nn点求得频率特性，H为N个点处的频率响应复值输出向量；F为与第N点处对应的频率值f(Hz)subplot(2,1,1);plot(f,20*log10(abs(H)));title('滤波器特性')%20log10(abs(H))DB为频响幅值曲线幅值xlabel('频率/Hz');ylabel('振幅/dB');gridon;%滤波器幅频特性subplot(2,1,2);plot(f,180/pi*unwrap(angle(H)));%angle(H)为频响相位曲线相位值xlabel('频率/Hz');ylabel('相位/^o');gridon;%滤波器相频特性巴特沃思低通滤波器的相位谱和振幅谱见图2.4.1图2.4.1利用脉冲响应不变法设计的数字巴特沃斯滤波器（w(单位：错误!未找到引用源。)使用巴特沃思低通滤波器对语音信号进行去噪处理，编写程序如下：dt=1/Fs;%采样间隔t=[0;(length(x1)-1)]*dt;%时间序列数字信号处理课程设计10figure(3)subplot(2,1,1);plot(x1);title('输入信号');y=filter(b,a,x1);%对输入信号进行滤波subplot(2,1,2);plot(y);title('输出信号');%绘制输出信号xlabel('时间/s')%比较滤波前后语音信号的波形Am=abs(y);pha=angle(y);figure(4);subplot(4,1,1);plot(Am);title('滤波后幅值')subplot(4,1,2);plot(pha);title('滤波后相位')subplot(4,1,3);plot(x2);title('滤波前幅值')subplot(4,1,4);plot(x3);title('滤波前相位')%比较滤波前后语音信号的频谱滤波前后的波形见图2.4.2图2..4.2滤波前后波形比较图滤波前后波形的频域幅值谱和相位谱的比较图2.4.3数字信号处理课程设计11图2.4.3滤波前后波形的频域幅值谱和相位谱的比较图%用麦克风采集一段8000Hz，8k的单声道语音信号，绘制波形并观察其频谱[x,Fs,nbits]=wavread('c.wav')figure(1)n=length(x)x1=fft(x,2000)x2=abs(x1)x3=angle(x1)subplot(4,1,1);plot(x);title('原信号时域')subplot(4,1,2);plot(x1);title('原信号fft')subplot(4,1,3);plot(x2);title('fft幅值')subplot(4,1,4);plot(x3);title('fft相角')%巴特沃思低通滤波器Fs=8000;%采样频率wp=1100*2/Fs;ws=1200*2/Fs;%根据采样频率将滤波器边界进行转化Rp=1;数字信号处理课程设计12Rs=20;%通带波纹和阻带波纹Nn=128;[N,Wn]=buttord(wp,ws,Rp,Rs);%求滤波器的最小阶数和3dB固有频率[b,a]=butter(N,Wn)%设计BUTTERWORTH低通滤波器，返回滤波器的系数矩阵figure(2);[H,f]=freqz(b,a,Nn,Fs);%用Nn点求得频率特性，H为N个点处的频率响应复值输出向量；F为与第N点处对应的频率值f(Hz)subplot(2,1,1);plot(f,20*log10(abs(H)));title('滤波器特性')%20log10(abs(H))DB为频响幅值曲线幅值xlabel('频率/Hz');ylabel('振幅/dB');gridon;%滤波器幅频特性subplot(2,1,2);plot(f,180/pi*unwrap(angle(H)));%angle(H)为频响相位曲线相位值xlabel('频率/Hz');ylabel('相位/^o');gridon;%滤波器相频特性dt=1/Fs;%采样间隔t=[0;(length(x1)-1)]*dt;%时间序列figure(3)subplot(2,1,1);plot(x1);title('输入信号');y=filter(b,a,x1);%对输入信号进行滤波subplot(2,1,2);plot(y);title('输出信号');%绘制输出信号xlabel('时间/s')%比较滤波前后语音信号的波形Am=abs(y);pha=angle(y);figure(4);subplot(4,1,1);plot(Am);title('滤波后幅值')subplot(4,1,2);plot(pha);title('滤波后相位')subplot(4,1,3);plot(x2);title('滤波前幅值')subplot(4,1,4);plot(x3);title('滤波前相位')%比较滤波前后语音信号的频谱数字信号处理课程设计132.5结果分析根据3.5.1的滤波前后波形的频域幅值谱和相位谱的比较图（图3.5.2）得：从幅值的大小来说，滤波后的幅值变小；从相位的稀疏度来说，滤波后的相位较为稀疏。由此可以得出，通过巴特沃思低通滤波器的滤波，滤除了采样文件中夹杂的噪声，滤波成功。三、设计原理结果与仿真分析3.1设计原理用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析。3.2结果与分析在这次的课程设计中我们设计的滤波器对语音信号的滤噪处理或多或少会因为我们的疏忽或者实践能力不够还是会遇