1、 DSP 课程设计报告课程设计报告 姓姓 名名 学学 院院 电气与电子工程学院电气与电子工程学院 专专 业业 电子信息科学与技术电子信息科学与技术 学学 号号 班班 级级 电信科学电信科学 091 班班 指导老师指导老师 2012 年 6 月 18 日 一、一、设计题目设计题目 语音信号处理，进行三种类型的滤波器的设计。 二、二、设计要求设计要求 (1) IIR低通滤波器性能指标 fc=1000 Hz,（通带截止频率） fst=2000 Hz, （阻带 截止频率）Rst=30 dB, Rp5 dB。 (2) IIR带阻,双线性变换法 设计指标 fst1=2000 Hz, fst2=3000 H
2、z, fc1=1000 Hz, fc2=4000 Hz, Rst=30 dB, Rp=5 dB。 (3) 带通滤波器性能指标 fc1=2000 Hz, fc2=3000 Hz, fst1=1000 Hz, fst2=4000 Hz, Rst=30 dB, %Rp5 dB FIR 数字滤波器 根据以上的技术指标设计出相应的滤波器，画出滤波器的频率响应曲线。 三、三、设计原理设计原理 1. 语音信号的采集 利用 Windows 下的录音机，录制一段自己的话音，时间在 1 s 内。然后在 Matlab 软件平台下，利用函数 wavread 对语音信号进行采样，记住采样频率和采 样点数。通过 wavr
3、ead 函数的使用，理解采样频率、采样位数等概念。 y,fs=wavread(d:111.wav,1000 60000); 其中 y为 wav文件的音调数据，长度为 60000-1000+1，fs 为该文件的播放频率。 通过 sound 函数播放该文件的声音： Sound(y,fs); 2. 语音信号的频谱分析 首先画出语音信号的时域波形，然后对语音号进行快速傅里叶变换，得到信 号的频谱特性，其程序如下： Y=fft(y,59001); Subplot(231);plot(y);title(滤波前的信号波形); Subplot(232);plot(abs(Y);title(滤波前的信号频谱);
4、 3. 设计数字滤波器和画出其频率响应 给出各滤波器的性能指标： (1) 低通滤波器性能指标 fc=1000 Hz,（通带截止频率） fst=2000 Hz, （阻带截 止频率）Rst=30 dB, Rp5 dB。 (2) 带阻,双线性变换法 IIR 巴特沃斯数字滤波器 设计指标 fst1=2000 Hz, fst2=3000 Hz, fc1=1000 Hz, fc2=4000 Hz, Rst=30 dB, Rp=5 dB。 (3) 带通滤波器性能指标 fc1=2000 Hz, fc2=3000 Hz, fst1=1000 Hz, fst2=4000 Hz, Rst=30 dB, %Rp5 d
5、B FIR 数字滤波器 根据以上的技术指标设计出相应的滤波器，画出滤波器的频率响应曲线。 4. 用滤波器对信号进行滤波 用自己设计的各滤波器分别对采集的信号利用函数 filter 对信号进行滤波， 并比较滤波前后语音信号的波形及频谱： x=filter(b,a,y); X=fft(x,59001); Subplot(235);plot(x);title(滤波后的信号波形); Subplot(236);plot(abs(X);title(滤波后的信号频谱); 5. 回放语音信号 在 Matlab 中，函数 sound 可以对声音进行回放，其调用格式：sound(x,fs)； 可以感觉滤波前后的声
6、音有变化。 6. 设计系统界面 为了使编制的程序操作方便，有能力的学生设计处理系统的用户界面。在所 设计的系统界面上可以选择滤波器的类型，输入滤波器的参数，显示滤波器的频 率响应，选择信号等。 四、四、源程序清单源程序清单 4.1 冲击响应不变法冲击响应不变法,IIR 低通滤波器性能指标低通滤波器性能指标 fc=1000 Hz,（通带截止频率）（通带截止频率） fst=2000 Hz, （阻带截止频率）（阻带截止频率）Rst=30 dB, Rp5 dB。 y,fs=wavread(d:111.wav,1000 60000); Sound(y,fs); Y=fft(y,59001); %时域信号转换为频域的信号时域信号转换为频域的信号 Subplot(231);plot(y);title(滤波前的信号波形滤波前的信号波形); Subplot(232);plot(abs(Y);title(滤波前的信号频谱滤波前的信号频谱);
