数字通信课程设计----音频信号分析及去噪.docVIP

课 程 设 计 报 告 课程名称： 数字通信课程设计 设计名称： 音频信号的谱分析及去噪 姓 名： 学 号: 课 程 设 计 任 务 书 学生班级： 通信0801 学生姓名： 学号： 设计名称： 音频信号的谱分析及去噪 起止日期： 2011.6.21-2011.7.3 指导教师： 胥 磊 设计要求： 基本要求： 录制一段音频（如歌曲，说话声等），采用Matlab工具对此音频信号用FFT作谱分析。 录制一段加入噪声的音频（如在歌声中加入尖锐的口哨声或者其他噪声），采用Matlab工具对此音频信号用FFT作谱分析。 选择合适的指标，设计FIR数字滤波器，将音频中加入的噪声信号减弱或滤除。 扩展要求： 将处理后的音频信号重新生成.wav文件，收听该音频，根据效果调整滤波器指标重新设计滤波器。 课 程 设 计 学 生 日 志 时间 设计内容 6月28日 查阅资料，明确题目内涵 6月29日 录制音频，编写程序 6月30日 检查程序，再次确认设计正确与否。开始写报告 7月4日 复习于答辩 课 程 设 计 考 勤 表 周 星期一 星期二 星期三 星期四 星期五 课 程 设 计 评 语 表 指导教师评语： 成绩： 指导教师： 年 月 日  音频信号的谱分析及去噪 设计目的和意义 1、录制音频，采用matlab工具对此信号做谱分析。 2、录制一段加噪的信号，对加噪信号进行谱分析。 3、选择合适的指标对加噪信号进行滤波。 4、学会使用wavread、wavwrite等函数。 设计原理 对于录制的音频信号的读入可以使用函数wavread。由于音频信号是连续的，故绘图应用plot函数。使用函数fft可以得到音频信号的频谱图。 由于音频信号是连续且长度未知，故可以采用巴特沃斯滤波器。滤掉高频部分的噪音，剩下的就是原信号了。 将去噪后的信号写成wav格式的文件可以使用wavwrite函数。 详细设计步骤 录制一段音频（如歌曲，说话声等），采用Matlab工具对此音频信号用FFT作谱分析。 录制一段音频信号并命名为fl.wav存放在课程设计1的文件夹中。 使用wavread函数读出此信号。 用函数FFT进行傅里叶变换，得到频谱图。 绘制时域图和频域图。 录制一段加入噪声的音频（如在歌声中加入尖锐的口哨声或者其他噪声），采用Matlab工具对此音频信号用FFT作谱分析。 重新录制1中的音频文件，录制时加入噪声信号。 2、使用wavread函数读出此信号。 3、用函数FFT进行傅里叶变换，得到频谱图。 4、绘制时域图和频域图。 选择合适的指标，设计FIR数字滤波器，将音频中加入的噪声信号减弱或滤除。 扩展要求： 根据2得到的频谱图，选择合适的通带截止频率和阻带截止频率。根据上述指标设计出合适的滤波器，将信号通过滤波器，得到去噪后的信号。滤波函数可用fftfilt。 将处理后的音频信号重新生成.wav文件，收听该音频，根据效果调整滤波器指标重新设计滤波器。 使用wavwrite函数将去噪后的函数重新生成wav格式的音频文件。 播放音频文件，根据效果修改滤波器中的参数，直到得到与原信号几乎相同的音频文件。 设计程序见附件。 附程序：附件 %得到无噪声的原始信号和有噪声的加噪信号 clc clear for n = 1 : 999 % a(n) = sin(0.5906*n)+sin(1.0245*n) a(n)=sin(0.5906*n); %无噪声信号 b(n)=sin(1.0245*n); end; a=a/2; sound(a); b=conv(a,b/2); %加噪 sound(b); wavwrite(a,a.wav); %写文件 wavwrite(b,b.wav); %没有噪声的