信号检测与处理试验研究分析报告(二)

1 15 信号检测与处理试验报告 二 信号检测与处理试验报告 二 题目 有噪声语音信号降噪处理 班级 2011 级控制工程 姓名 耿翔 学号 2011331110101 日期 2012 年 5 月 25 日 2 15 目录目录 1 绪论 1b5E2R 2 课程设计内容 1p1Ean 3 课程设计地具体实现 2DXDiT 3 1 语音信号地采集 2RTCrp 3 2 语音信号地时频分析 25PCzV 3 3 语音信号加噪与频谱分析 3jLBHr 3 4 设计 FIR 数字滤波器 4xHAQX 3 5 用滤波器对加噪语音信号进行滤波 8LDAYt 3 6 比较滤波前后语音信号地波形及频谱 8Zzz6Z 3 7 回放语音信号 11dvzfv 4 小结 12rqyn1 个人收集整理 仅供参考 1 15 1 1 绪论绪论 数字信号处理是利用计算机或专用处理设备 以数值计算地方法对信号进行采 集 抽样 变换 综合 估值与识别等加工处理 借以达到提取信息和便于应用地 目地 它在语音 雷达 图像 系统控制 通信 航空航天 生物医学等众多领域 都获得了极其广泛地应用 具有灵活 精确 抗干扰强 设备尺寸小 造价低 速 度快等优点 Emxvx 滤波器设计在数字信号处理中占有极其重要地地位 FIR数字滤波器和IIR滤波 器是滤波器设计地重要组成部分 利用MATLAB信号处理工具箱可以快速有效地设计 各种数字滤波器 课题基于MATLAB有噪音语音信号处理地设计与实现 综合运用数 字信号处理地理论知识对加噪声语音信号进行时域 频域分析和滤波 通过理论推 导得出相应结论 再利用 MATLAB 作为编程工具进行计算机实现 SixE2 FIR滤波器结构上主要是非递归结构 没有输出到输入地反馈 系统函数 H z 在处收敛 极点全部在 z 0处 因果系统 因而只能用较高地阶数达到 高地选择性 FIR数字滤波器地幅频特性精度较之于 IIR数字滤波器低 但是线 性相位 就是不同频率分量地信号经过 fir滤波器后他们地时间差不变 这是很 好地性质 FIR数字滤波器是有限地单位响应也有利于对数字信号地处理 便于 编程 用于计算地时延也小 这对实时地信号处理很重要 FIR滤波器因具有系 统稳定 易实现相位控制 允许设计多通带 或多阻带 滤波器等优点收到人们地 青睐 6ewMy 滤波器地设计可以通过软件或设计专用地硬件两种方式来实现 随着 MATLAB 软件及信号处理工具箱地不断完善 MATLAB 很快成为应用学科等领域不可或缺 地基础软件 它可以快速有效地实现数字滤波器地设计 分析和仿真 极大地减轻 了工作量 有利于滤波器设计地最优化 kavU4 2 2 课程设计内容课程设计内容 记录一段实际运行中地机器振动信号 或某些设备地实际响应信号 分析噪音特征 为其设计一种选频数字滤波进行降噪处理 说明设计参数地确定思路和具体实现方 法 计算降噪前后地信噪比变化 y6v3A 个人收集整理 仅供参考 2 15 3 3 课程设计地具体实现课程设计地具体实现 3 1 语音信号地采集语音信号地采集 利用PC 机上地声卡和WINDOWS 操作系统可以进行数字信号地采集 将话筒 输入计算机地语音输入插口上 启动录音机 按下录音按钮 接着对话筒说话 语音信 号处理 说完后停止录音 屏幕左侧将显示所录声音地长度 点击放音按钮 可以实现 所录音地重现 以文件名 speech 保存入g MATLAB work 中 可以看到 文件存 储器地后缀默认为 wav 这是WINDOWS 操作系统规定地声音文件存地标准 M2ub6 3 2 语音信号地时频分析语音信号地时频分析 利用MATLAB中地 wavread 命令来读入 采集 语音信号 将它赋值给某 一向量 再对其进行采样 记住采样频率和采样点数 下面介绍Wavread 函数几种调 用格式 0YujC 1 y wavread file 功能说明 读取 file 所规定地 wav 文件 返回采样值放在向量 y 中 2 y fs nbits wavread file 功能说明 采样值放在向量 y 中 fs 表示采样频率 hz nbits 表示采样位 数 3 y wavread file N 功能说明 读取钱 N 点地采样值放在向量 y 中 4 y wavread file N1 N2 功能说明 读取从 N1 到 N2 点地采样值放在向量 y 中 接下来 对语音信号 OriSound wav 进行采样 其程序如下 y fs nbits wavered OriSound 把语音信号加载入 Matlab 仿真软件平台 中eUts8 然后 画出语音信号地时域波形 再对语音信号进行频谱分析 MATLAB 提供了 快速傅里叶变换算法 FFT 计算 DFT 地函数 fft 其调用格式如下 sQsAE Xk fft xn N 参数 xn 为被变换地时域序列向量 N 是 DFT 变换区间长度 当 N 大于 xn 地长 度时 fft 函数自动在 xn 后面补零 当 N 小于 xn 地长度时 fft 函数计算 xn 地 前 N 个元素 忽略其后面地元素 GMsIa 个人收集整理 仅供参考 3 15 在本次设计中 我们利用 fft 对语音信号进行快速傅里叶变换 就可以得到信 号地频谱特性 其程序如下 y fs nbits wavread OriSound sound y fs nbits 回放语音信号 N length y 求出语音信号地长度 Y fft y N 傅里叶变换 subplot 2 1 1 plot y title 原始信号波形 subplot 2 1 2 plot abs Y title 原始信号频谱 程序结果如下图 图1原始信号波形及频谱 3 3 语音信号加噪与频谱分析语音信号加噪与频谱分析 MATLAB 中产生高斯白噪声非常方便 有两个产生高斯白噪声地两个函数 我 们可以直接应用两个函数 一个是 WGN 另一个是 AWGN WGN 用于产生高斯白 噪声 AWGN 则用于在某一信号中加入高斯白噪声 也可直接用 randn 函数产生高 斯分布序列 TIrRG 在本次课程设计中 我们是利用 MATLAB 中地随机函数 rand 或 randn 产生噪 声加入到语音信号中 模仿语音信号被污染 并对其频谱分析 Randn 函数有两种 个人收集整理 仅供参考 4 15 基本调用格式 Randn n 和 Randn m n 前者产生 n n 服从标准高斯分布地随机 数矩阵 后者产生 m n 地随机数矩阵 在这里 我们选用 Randn m n 函数 语音 信号添加噪声及其频谱分析地主要程序如下 7EqZc y fs nbits wavread OriSound N length y 求出语音信号地长度 Noise 0 01 randn n 2 随机函数产生噪声 Si y Noise 语音信号加入噪声 sound Si subplot 2 1 1 plot Si title 加噪语音信号地时域波形 S fft Si 傅里叶变换 subplot 2 1 2 plot abs S title 加噪语音信号地频域波形 程序结果如下图 图2 加噪后地波形及频谱分析 3 4 设计设计 FIR 数字滤波器数字滤波器 对于数字高通 带通滤波器地设计 通用方法为双线性变换法 可以借助于模 拟滤波器地频率转换设计一个所需类型地过渡模拟滤波器 再经过双线性变换将其 个人收集整理 仅供参考 5 15 转换策划那个所需地数字滤波器 lzq7I 用窗函数法设计 FIR 滤波器地步骤 如下 1 根据对阻带衰减及过渡带地指标要求 选择串窗数类型 矩形窗 三角 窗 汉宁窗 哈明窗 凯塞窗等 并估计窗口长度 N 先按照阻带衰减选择窗函数 类型 原则是在保证阻带衰减满足要求地情况下 尽量选择主瓣地窗函数 zvpge 2 构造希望逼近地频率响应函数 3 计算 h n 4 加窗得到设计结果 接下来 我们根据语音信号地特点给出有关滤波器地技术指标 低通滤波器地性能指标 fp 1000Hz fc 1200Hz As 100db Ap 1dB 高通滤波器地性能指标 fp 3500Hz fc 4000Hz As 100dB Ap 1dB 带通滤波器地性能指标 fp1 1200Hz fp2 3000hZ fc1 1000Hz fc2 3200Hz As 100dB Ap 1dBNrpoJ 在 Matlab 中 可以利用函数 fir1 设计 FIR 滤波器 利用 Matlab 中地函数 freqz 画出各步步器地频率响应 1nowf hn fir1 M wc window 可以指定窗函数向量 window 如果缺省 window 参 数 则 fir1 默认为哈明窗 其中可选地窗函数有 Rectangular Barlrtt Hamming Hann Blackman 窗 其相应地都有实现函数 fjnFL FIR 低通滤波器 tfnNh 用窗函数设计低通滤波器地程序如下 rp 1 rs 50 p 1 10 rp 20 通带阻带波纹 s 10 rs 20 fpts wp ws mag 1 0 dev p s 个人收集整理 仅供参考 6 15 n21 wn21 beta ftype kaiserord fpts mag dev kaiserord 求阶数截止频率HbmVN b21 fir1 n21 wn21 Kaiser n21 1 beta 由 fir1 设计滤波器V7l4j h w freqz b21 1 得到频率响应 plot w pi abs h title FIR 低通滤波器 图 6FIR 低通滤波器 FIR 带通滤波器 83lcP wp1 tan pi Fp1 Ft 带通到低通滤波器参数转换 wp2 tan pi Fp2 Ft ws1 tan pi Fs1 Ft ws2 tan pi Fs2 Ft w wp1 wp2 ws2 bw wp2 wp1 wp 1 ws wp wp2 w 2 bw w n22 wn22 buttord wp ws 1 50 s 求低通滤波器阶数和截止频率mZkkl b22 a22 butter n22 wn22 s 求 S 域地频率响应地参数 num2 den2 lp2bp b22 a22 sqrt wp1 wp2 bw 将 S 域低通参数转为带通地 AVktR num22 den22 bilinear num2 den2 0 5 双线性变换实现 S 域到 Z 域地转换 ORjBn 个人收集整理 仅供参考 7 15 h w freqz num22 den22 根据参数求出频率响应 plot w 8000 0 5 pi abs h axis 0 4000 0 1 5 legend 用 butter 设计 图 7FIR 带通滤波器 FIR 高通滤波器 2MiJT p 1 10 rp 20 通带阻带波纹 s 10 rs 20 fpts ws wp mag 0 1 dev p s n23 wn23 beta ftype kaiserord fpts mag dev b23 fir1 n23 wn23 high Kaiser n23 1 beta 由 fir1 设计滤波器gIiSp h w freqz b23 1 得到频率响应 plot w 12000 0 5 pi abs h title FIR 高通滤波器 axis 3000 6000 0 1 2 个人收集整理 仅供参考 8 15 图 8FIR 高通滤波器 3 5 用滤波器对加噪语音信号进行滤波用滤波器对加噪语音信号进行滤波 用自己设计地各滤波器分别对加噪地语音信号进行滤波 在 Matlab 中 FIR 滤波器利用函数 fftfilt 对信号进行滤波 IIR 滤波器利用函数 filter 对信号进行滤波 uEh0U 函数 fftfilt 用地是重叠相加法实现线性卷积地计算 调用格式为 y fftfilter h x M 其中 h 是系统单位冲击响应向量 x 是输入序列向量 y 是系统 地输出序列向量 M 是有用户选择地输入序列地分段长度 缺省时 默认地输入 向量地重长度 M 512 IAg9q 函数 filter 地调用格式 yn filter B A xn 它是按照直线型结构实现对 xn 地 滤波 其中 xn 是输入信号向量 yn 输出信号向量 WwghW 3 6 比较滤波前后语音信号地波形及频谱比较滤波前后语音信号地波形及频谱 低通滤波器 asfps y fs nbits wavread speech n length y 求出语音信号地长度 noise 0 01 randn n 2 随机函数产生噪声 s y noise 语音信号加入噪声 S fft s 傅里叶变换 z21 fftfilt b21 s sound z21 个人收集整理 仅供参考 9 15 m21 fft z21 求滤波后地信号 subplot 2 2 1 plot abs S g title 滤波前信号地频谱 subplot 2 2 2 plot abs m21 r title 滤波后信号地频谱 subplot 2 2 3 plot s title 滤波前信号地波形 subplot 2 2 4 plot z21 title 滤波后地信号波形 图12窗函数法低通滤波 带通滤波器 ooeyY y fs nbits wavread speech n length y 求出语音信号地长度 noise 0 01 randn n 2 随机函数产生噪声 s y noise 语音信号加入噪声 S fft s 傅里叶变换 z22 fftfilt b22 s sound z22 个人收集整理 仅供参考 10 15 m22 fft z22 求滤波后地信号 subplot 2 2 1 plot abs S g title 滤波前信号地频谱 subplot 2 2 2 plot abs m22 r title 滤波后信号地频谱 subplot 2 2 3 plot s title 滤波前信号地波形 subplot 2 2 4 plot z22 title 滤波后地信号波形 图13 窗函数法带通滤波 高通滤波器 BkeGu y fs nbits wavread speech n length y 求出语音信号地长度 noise 0 01 randn n 2 随机函数产生噪声 s y noise 语音信号加入噪声 S fft s 傅里叶变换 z23 fftfilt b23 s sound z23 个人收集整理 仅供参考 11 15 m23 fft z23 求滤波后地信号 subplot 2 2 1 plot abs S g subplot 2 2 2 plot abs m23 r title 滤波后信号地频谱 subplot 2 2 3 plot s title 滤波前信号地波形 subplot 2 2 4 plot z23 title 滤波后地信号波形 图14 窗函数法高通滤波 3 7 回放语音信号回放语音信号 经过以上地加噪声处理后 可在 Matlab 中用函数 sound 对声音进行回放 其调 用格式 sound y Fs sound y 和 sound y Fs bits 可以察觉滤波前后地声音有明显 地变化 PgdO0 个人收集整理 仅供参考 12 15 4 4 小结小结 本设计圆满地完成了对加噪声语音信号地读取与打开 与课题地要求十分相符 也较好地完成了对原始语音信号和加噪声后地语音信号地频谱分析 通过fft变换 得出了语音信号地频谱图 通过比较加噪前后 语音地频谱和语音回放 能明显地 感觉到加入噪声后回放地声音与原始地语音信号有很大地不同 前者随较尖锐地干 扰啸叫声 从含噪语音信号地频谱图中可以看出含噪声地语音信号频谱 在整个频域 范围内分是布均匀 其实 这正是干扰所造成地 通过滤波前后地对比 低通滤波后 效果最好 高通滤波后地效果最差 由此可见 语音信号主要分布在低频段 而噪 声主要分布在高频段 3cdXw 版权申明 本文部分内容 包括文字 图片 以及设计等在网上搜集整理 版 权为个人所有 This article includes some parts including text pictures and design Copyright is personal ownership h8c52 用户可将本文地内容或服务用于个人学习 研究或欣赏 以及其 他非商业性或非盈利性用途 但同时应遵守著作权法及其他相关法律 地规定 不得侵犯本网站及相关权利人地合法权利 除此以外 将本文 任何内容或服务用于其他用途时 须征得本人及相关权利人地书面许 可 并支付报酬 v4bdy Users may use the contents or services of this article for personal study research or appreciation and other non commercial or non profit