《数字信号处理课程设计-基于FIR的有噪声语音信号处理》.docVIP

1 课程设计的主要目的和内容及要求 课程设计目的： 1. 通过本次课程设计，综合运用数字信号处理技术课程和其他有关先修课程的理论和生产实际知识去分析和解决具体问题，并使所学知识得到进一步巩固、深化和发展。 2. 初步培养学生对工程设计的独立工作能力，学习设计的一般方法。 3. 通过课程设计树立正确的设计思想，提高学生分析问题、解决问题的能力。 4. 进行设计基本技能的训练，如查阅设计资料和手册、程序的设计、调试等。 本课题设计内容： 1. 利用Windows下的录音机或其他软件，选择Windows系统的“叮······”（ding.wav），并对该信号进行采样； 2. 语音信号的频谱分析，画出采样后语音信号的时域波形和频谱图； 3. 产生噪声信号并加到语音信号中，得到被污染的语音信号，并回放语音信号； 4. 污染信号的频谱分析，画出被污染的语音信号时域波形和频谱； 5. 根据有关的频谱特性，采用Kaiser窗设计FIR数字低通，高通，带通，带阻滤波器并画出4种滤波器的幅频、相频图； 6. 用自己设计的4种滤波器分别对被不同噪声污染的信号进行滤波； 7. 分析得到信号的频谱，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化； 8. 回放语音信号。 设计要求： 1. 熟悉离散信号和系统的时域特性。 2. 掌握数字信号处理的基本概念，基本理论和基本方法。 3. 掌握序列快速傅里叶变换方法。 4. 学会MATLAB的使用，掌握MATLAB的程序设计方法。 5. 掌握利用MATLAB对语音信号进行频谱分析。 6. 掌握MATLAB设计FIR数字滤波器的方法和对信号进行滤波的方法。2 课程设计的总体方案 2.1设计框图 x x x x xxx2.2总体设计概述 本次课程设计所采用的主要软件为Matlab。设计中用到的函数主要有wavread,sound,fft,subplot,plot,fir1,freqz,filter。 设计的主要内容为数字滤波器的设计。从功能上分可把滤波器分为低通，高通，带阻，带通滤波器。同时数字滤波器的设计又可分为无限脉冲响应（IIR）数字滤波器设计和有限脉冲响应（FIR）数字滤波器。IIR数字滤波器设计方法主要可分为脉冲响应不变法和双线性变换法，其中主要有巴特沃斯滤波器，切比雪夫I型滤波器，切比雪夫II型滤波器，椭圆滤波器以及贝塞尔滤波器等五种模拟滤波器做为原型；而对于FIR数字滤波器设计方法主要有窗函数法，频率采样法和等波纹最佳逼近法，其中窗函数法中主要有六种典型的窗函数，它们分别是矩形窗，三角形窗，汉宁（Hanning）窗，哈明（Hamming）窗，布莱克曼（Blackman）窗以及凯塞—贝塞尔窗（Kaiser-Basel Window）。本课程设计采用的是FIR数字滤波器中的窗函数法，选用的是凯塞－贝塞尔窗.主要的设计思路为首先获取原始信号x。其次在画出x的时域图和经采样后的频域图后分别加入四种噪声信号noise1,noise2,noise3,noise4分别获得X1,X2,X3,X4四种带噪声信号。然后分别画出X1,X2,X3,X4的时域图并进行采样画出频域图。再次，利用Matlab中的函数fir1,Kaiser,freqz,filter分别设计出低通，高通，带阻，带通四种数字滤波器并分别画出它们的幅频特性和相频特性。最后分别让X1，X2，X3，X4四种信号分别通过四种滤波器得到四种滤波后的信号y1,y2,y3,y4，分别画出y1,y2,y3,y4四种信号的时域图和经采样后的频域图与X1,X2,X3,X4进行比较。 3 设计的详细原理 设计主要可分为原始信号部分，噪声信号部分，低通滤波部分，高通滤波部分，带阻滤波部分，带通滤波部分等六个部分。 3.1原始信号部分 本次使用的原始信号是windows系统的“叮······”（ding.wav）语音信号，可利用Windows下的录音机（wavread）或其他软件进行获取。利用sound函数进行试听。其中wavread的主要参数有采样频率（Fs），采样信号（x），采样位数（bits），本次设计中Fs＝10000，bits＝16从而可获得时间长为1s的语音信号x。sound函数调用格式为sound(x,Fs,bits) 。对x信号进行采样y=fft(x,m)，其中m=length(x)。最后利用plot函数分别画出x信号的时域图和频域图subplot(211);plot(t,x);title(原始信号波形);xlabel(time(s));subplot(212);plot(f,abs(y));title(原始信号频谱);xlabel(frequency(hz));便可完成原始信号部分的设计。同时通过图像可观察x的特性