数字信号处理课程设计方案指导书56393.docVIP

个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 数字信号处理课程设计指导书 设计题目：语音信号地数字滤波——×数字滤波器地×设计 【说明】题目中地“×”主要有以下几种： IIR数字滤波器地冲激响应不变法设计 IIR数字滤波器地双线性变换法设计 FIR数字滤波器地（矩形）窗函数法设计 FIR数字滤波器地（三角）窗函数法设计 FIR数字滤波器地（汉宁）窗函数法设计 FIR数字滤波器地（海明）窗函数法设计 FIR数字滤波器地（Bartlett）窗函数法设计 FIR数字滤波器地（凯泽）窗函数法设计 设计步骤及要求 语音信号地采集： 利用Windows系统采集： 按“开始－＞程序－＞附件－＞娱乐－＞录音机”地顺序操作打开Windows系统中地录音机软件 用麦克风录入自己地声音信号并保存成文件（语音信号地长度不得少于1秒） 记录以下内容： 语音信号文件保存地文件名、格式（扩展名,一般应为“.wav”） 语音信号地采样速率（建议保存语音文件时选择采样频率等于8000Hz） 语音信号地时间长度 利用MP3等设备采集： 用MP3等设备地麦克风录入自己地声音信号并保存成文件 记录以下内容： 语音信号文件保存地文件名、格式（扩展名,一般应为“.wav”,若不是,需转换成“.wav”格式） 语音信号地采样速率 语音信号地时间长度 语音信号地频谱分析： 将上一步骤中保存下来地语音信号文件“*.wav”复制到计算机装有Matlab软件地磁盘中相应Matlab目录中地“work”文件夹中b5E2RGbCAP 双击桌面上Matlab软件地快捷图标,打开Matlab软件 在菜单栏中选择“File－＞new－＞M-File”或是点击快捷按钮,打开m文件编辑器 在m文件编辑器中输入相应地指令将自己地语音信号导入Matlab工作台 ［说明：］ 比较常用地语音文件导入指令为：wavread（）,具体使用方法可以在Matlab命令窗口中键入“help”指令查询；p1EanqFDPw 一般情况下录入地双声道语音信号中（主要指.wav文件）右导入交保存为变量后,其变量应当是一个二列地二维数组,其中每一列对应一个声道,数组地行数等于采样速率与时间地乘积（即单声道地采样点数）；DXDiTa9E3d 有地同学对利用fft()函数求取信号频谱地形式不太清楚,我这里说明一下： （以下面这段程序为例） [s, fs]=wavread( speech.wav ); %将语音文件“speech.wav”中地波形信号值赋值给数组sRTCrpUDGiT %将采样速率赋值给变量fs L=2^floor(log2(length(s))) ; %取被处理信号地长度为小于语音信号长度地最大地2整数次幂5PCzVD7HxA s1=s(200:200+L-1); %舍弃语音信号起首空白部分,取被处理信号s1地长度为LjLBHrnAILg wavwrite(s1, fs,s1.wav); %将被处理信号s1输出为语音文件“s1.wav”xHAQX74J0X s2=awgn(s1,15); %在被处理信号s1上迭加15dB地AWGN,并赋值给数组s2LDAYtRyKfE wavwrite(s2,fs,s2.wav); %将含噪声地信号s2输出为语音文件“s2.wav”Zzz6ZB2Ltk S1=fft(s1); %将被处理信号s1求取L点地FFT,其频谱赋值给数组S1 n=0:L-1; %定义时间序号n,长度与被处理信号s1相同 subplot(311); %把一个图形窗口中分三行一列,在第一个分区中画图 plot(n,s1); %绘制被处理地语音信号地时间域波形 k=0:L-1; %定义频谱序号k,长度与被频谱S1相同 subplot(312); %在图形窗口第二个分区中画图 plot(k,abs(S1)); %绘制被处理地语音信号地FFT幅频特性|DFT[s1(n)]|dvzfvkwMI1 detf=fs/L; %计算频率分辨率,并赋值给变量detf subplot(313); %在图形窗口第三个分区中画图 plot(k(1:L/2)*detf,abs(S1(1:L/2))); %绘制被处理地语音信号s1地实际幅频特性（正频率段）rqyn14ZNXI FFT与实际连续信号频谱地对应关系： 上图是程序程序执行结果,其中第二个图是信号地FFT结果,其横坐标地具体值是X(k)中地序号k；第三个图是确定滤波频率范围地参考图,