数字信号处理matlab滤波器课程设计.doc

课 程 设 计 任 务 书 题 目 《数字信号处理》课程设计 学院(部) 专 业 班 级 学生姓名 学 号 2010.12.27 至 2011.01.02 共 1 周 指导教师(签字) 数字信号处理综合设计 一、实验目的 1．学会MATLAB的使用，掌握MATLAB的程序设计方法； 2．掌握在Windows环境下语音信号采集的方法； 3．掌握数字信号处理的基本概念、基本理论和基本方法； 4．掌握MATLAB设计FIR和IIR数字滤波器的方法； 5．学会用MATLAB对信号进行分析和处理。 二、实验原理 参考《数字信号处理》教材。 三、主要实验仪器及材料 微型计算机、Matlab6.5教学版、TC编程环境。 四、实验内容 1．语音信号的采集 要求利用windows下的录音机或其他软件，录制一段自己的话音，时间控制在1秒左右。然后在MATLAB软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数。通过wavread函数的使用，要求理解采样频率、采样位数等概念。 wavread函数调用格式： y=wavread(file)，读取file所规定的wav文件，返回采样值放在向量y中。 [y,fs,nbits]=wavread(file)，采样值放在向量y中，fs表示采样频率（Hz），nbits表示采样位数。 y=wavread(file,N)，读取前N点的采样值放在向量y中。 y=wavread(file,[N1,N2])，读取从N1点到N2点的采样值放在向量y中。 2．语音信号的频谱分析 要求首先画出语音信号的时域波形；然后对语音信号进行频谱分析，在MATLAB中，可以利用函数fft对信号进行快速付立叶变换，得到信号的频谱特性；从而加深对频谱特性的理解。 3．设计数字滤波器和画出频率响应 根据语音信号的特点给出有关滤波器的性能指标：1）低通滤波器性能指标，fp=1000Hz，fc=1200 Hz， As=100dB，Ap=1dB；2）高通滤波器性能指标，fc=4800 Hz，fp=5000 Hz As=100dB，Ap=1dB；3）带通滤波器性能指标，fp1=1200 Hz，fp2=3000 Hz，fc1=1000 Hz，fc2=3200 Hz，As=100dB，Ap=1dB。要求学生首先用窗函数法设计上面要求的三种滤波器，在MATLAB中，可以利用函数fir1设计FIR滤波器；然后在用双线性变换法设计上面要求的三种滤波器，在MATLAB中，可以利用函数butte、cheby1和ellip设计IIR滤波器；最后，利用MATLAB中的函数freqz画出各滤波器的频率响应。 4．用滤波器对信号进行滤波 比较两种滤波器的性能，然后用性能好的各滤波器分别对采集的信号进行滤波，在MATLAB中，FIR滤波器利用函数fftfilt对信号进行滤波，IIR滤波器利用函数filter对信号进行滤波。 5．比较滤波前后语音信号的波形及频谱 要求在一个窗口同时画出滤波前后的波形及频谱。 6．回放语音信号 在MATLAB中，函数sound可以对声音进行回放。其调用格式： sound(x,fs,bits)； 可以感觉滤波前后的声音有变化。 五、实验思考 1．双线性变换法中Ω和ω之间的关系是非线性的，在实验中你注意到这种非线性关系了吗？从那几种数字滤波器的幅频特性曲线中可以观察到这种非线性关系？ 2．能否利用公式完成脉冲响应不变法的数字滤波器设计？为什么？ 六、实验报告要求 1.简述实验原理及目的。 2.按照实验步骤及要求，比较各种情况下的滤波性能。 3.总结实验所得主要结论。 4.简要回答思考题。 语音信号的数字滤波 ——利用双线性变换法实现IIR数字滤波器的设计 一．课程设计的目的 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法，掌握利用数字滤波器对信号进行滤波的方法，掌握数字滤波器的计算机仿真方法，并能够对设计结果加以分析。 二．课程设计原理 1．用窗函数法设计FIR滤波器 根据过渡带宽及阻带衰减要求，选择窗函数的类型并估计窗口长度N（或阶数