安建工数字信号处理课程设计+matlab滤波器的设计.docVIP

数字信号处理综合设计 一、实验目的 1.学会MA TLAB的使用,掌握MA TLAB的程序设计方法; 2.掌握在Windows环境下语音信号采集的方法; 3.掌握数字信号处理的基本概念、基本理论和基本方法; 4.掌握MA TLAB设计FIR和IIR数字滤波器的方法; 5.学会用MATLAB对信号进行分析和处理。 二、实验原理 参考《数字信号处理》教材。 三、主要实验仪器及材料 微型计算机、Matlab6.5教学版、TC编程环境。 四、实验内容 1.语音信号的采集 要求利用windows下的录音机(开始—程序—附件—娱乐—录音机,文件—属性—立即转换—8000KHz,8位,单声道或其他软件,录制一段自己的话音,时间控制在1秒左右。然后在MA TLAB软件平台下,利用函数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.语音信号的频谱分析 要求首先画出语音信号的时域波形;然后对语音信号进行频谱分析,在MA TLAB中,可以利用函数fft对信号进行快速付立叶变换,得到信号的频谱特性;从而加深对频谱特性 的理解。 3.设计数字滤波器和画出频率响应 根据语音信号的特点给出有关滤波器的性能指标:1低通滤波器性能指标,fp=1000Hz, fc=1200 Hz,As=100dB,Ap=1dB;2高通滤波器性能指标,fc=2800 Hz,fp=3000 Hz As=100dB,Ap=1dB;3带通滤波器性能指标,fp1=1200 Hz,fp2=3000 Hz,fc1=1000 Hz, fc2=3200 Hz,As=100dB,Ap=1dB。要求学生首先用窗函数法设计上面要求的三种滤波器,在MA TLAB中,可以利用函数fir1设计FIR滤波器,然后在用双线性变换法设计上面要求的三种滤波器;之后再利用函数butter和cheby1设计上面要求的三种IIR滤波器。最后,利用MA TLAB中的函数freqz画出各滤波器的频率响应。 4.用滤波器对信号进行滤波 比较FIR和IIR两种滤波器的性能,然后用性能好的各滤波器分别对采集的信号进行滤波,在MA TLAB中,FIR滤波器利用函数fftfilt对信号进行滤波,IIR滤波器利用函数filter 对信号进行滤波。 5.比较滤波前后语音信号的波形及频谱 要求在一个窗口同时画出滤波前后的波形及频谱。 6.回放语音信号 在MA TLAB中,函数sound可以对声音进行回放。其调用格式: sound(x,fs,bits; 可以感觉滤波前后的声音有变化。 五、实验思考 1.双线性变换法中Ω和ω之间的关系是非线性的,在实验中你注意到这种非线性关系了吗?从哪几种数字滤波器的幅频特性曲线中可以观察到这种非线性关系? 2.能否利用公式完成脉冲响应不变法的数字滤波器设计?为什么? 六、实验报告要求 1.简述实验原理及目的。 2.按照实验步骤及要求,比较各种情况下的滤波性能。 3.总结实验所得主要结论。 4.简要回答思考题。 语音信号的数字滤波 ——利用双线性变换法实现IIR 数字滤波器的设计 一.课程设计的目的 通过对常用数字滤波器的设计和实现,掌握数字信号处理的工作原理及设计方法;熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法,掌握利用数字滤波器对信号进行滤波的方法,掌握数字滤波器的计算机仿真方法,并能够对设计结果加以分析。 二.课程设计原理 1.用窗函数法设计FIR 滤波器 根据过渡带宽及阻带衰减要求,选择窗函数的类型并估计窗口长度N (或阶数M=N-1,窗函数类型可根据最小阻带衰减As 独立选择,因为窗口长度N 对最小阻带衰减As 没有影响,在确定窗函数类型以后,可根据过渡带宽小于给定指标确定所拟用的窗函数的窗口长度N ,设待求滤波器的过渡带宽为Δw ,它与窗口长度N 近似成反比,窗函数类型确定后,其计算公式也确定了,不过这些公式是近似的,得出的窗口长度还要在计算中逐步修正,原则是在保证阻带衰减满足要求的情况下,尽量选择较小的N ,在N 和窗函数类型确定后,即可调用MATLAB 中的窗函数求出窗函数wd (n 。 根据待求滤波器的理想频率响应求出理想单位脉冲响应hd(n,如果给出待求滤波器频率应为Hd ,则理想的单位脉冲响应可以用下面的傅里