MATLAB课程设计报告书

精品文档 河南理工大学物理与电子信息学院河南理工大学物理与电子信息学院 MATLABMATLAB 基础与应用课程论文基础与应用课程论文 2016 2017 学年第 二 学期 学生姓名 胥守振 学 号 311508000123 专业班级 电科 1501 2017 年 6 月 28 日 教师评语 评定成绩 精品文档 目 录 摘 要 1 设计原理 1 1 1 设计的目的及要求 2 1 2 课题的研究意义 2 2 设计原理 3 2 1 采样频率 3 2 2 采样位数 3 2 3 采样定理 3 2 4 时域信号的 FFT 分析 3 2 5 数字滤波器设计原理和方法 4 2 6 各种不同类型滤波器的性能比较 4 3 设计内容 4 3 1 语音信号的录入与提取 5 3 2 加噪处理 高频噪音 6 3 3 设计窗函数带阻滤波器 8 3 4 滤波处理 9 3 5 加躁处理 低频噪音 11 3 6 设计椭圆函数高通滤波器滤波处理 13 参考文献 15 摘摘 要要 语音信号处理是研究用数字信号处理技术和语音学知识对语音信号进行处理的新 兴的学科 是目前发展最为迅速的信息科学研究领域的核心技术之一 通过语音传递 精品文档 信息是人类最重要 最有效 最常用和最方便的交换信息形式 Matlab 语言是一种数据分析和处理功能十分强大的计算机应用软件 它可以将声 音文件变换为离散的数据文件 然后利用其强大的矩阵运算能力处理数据 如数字滤 波 傅里叶变换 时域和频域分析 声音回放以及各种图的呈现等 它的信号处理与 分析工具箱为语音信号分析提供了十分丰富的功能函数 利用这些功能函数可以快捷 而又方便地完成语音信号的处理和分析以及信号的可视化 使人机交互更加便捷 信 号处理是 Matlab 重要应用的领域之一 本实验设计用电脑自带的录音机采集了一段语音 对其进行了时域分析 频谱分 析 分析语音信号的特性 并应用 matlab 平台对语音信号加入了不同的噪声 进一步 用窗函数法 椭圆函数法分别设计了一个带阻滤波器和一个高通滤波器 然后对加噪 的语音信号进行滤波处理 最后对比滤波前后的语音信号的时域和频域特性 回放加 噪语音信号和去噪语音信号 对比研究处理前和处理后的声音的不同 关键词 语音信号 频域特性 时域特性 滤波器 1 1 设计目的和要求设计目的和要求 1 11 1 设计目的及要求设计目的及要求 本次课程设计要求利用 MATLAB 对语音信号进行处理和分析 要求学生采集语音信 号后 在 MATLAB 软件平台进行频谱分析 并对速配采集的语音信号加入干扰噪声 对 加入噪声的信号进行频谱分析 设计合适的滤波器滤除噪音 恢复信号 1 21 2 课题的研究意义课题的研究意义 语言是我们人类所特有的功能 它是传承和记载人类几千年文明史 没有语言就 没有我们今天人类的文明 语音是语言最基本的表现形式 是相互传递信息最重要的 精品文档 手段 是人类最重要 最有效 最常用和最方便的交换信息的形式 语音信号处理属于信息科学的一个重要分支 大规模集成技术的高度发展和计算 机技术的飞速前进 推动了这一技术的发展 它是研究用数字信号处理技术对语音信 号进行处理的一门新兴学科 同时又是综合性的多学科领域和涉及面很广的交叉学科 因此我们进行语言信号处理具有时代的意义 2 2 设计原理设计原理 2 12 1 采样频率采样频率 采样频率 也称为采样速度或者采样率 定义了每秒从连续信号中提取并组成离 散信号的采样个数 它用赫兹 Hz 来表示 采样频率的倒数是采样周期或者叫作采 样时间 它是采样之间的时间间隔 通俗的讲采样频率是指计算机每秒钟采集多少个 声音样本 是描述声音文件的音质 音调 衡量声卡 声音文件的质量标准 采样频 率越高 即采样的间隔时间越短 则在单位时间内计算机得到的声音样本数据就越多 对声音波形的表示也越精确 采样频率与声音频率之间有一定的关系 根据奎斯特理 论 只有采样频率高于声音信号最高频率的两倍时 才能把数字信号表示的声音还原 成为原来的声音 这就是说采样频率是衡量声卡采集 记录和还原声音文件的质量标 准 2 22 2 采样位数采样位数 采样位数即采样值或取样值 用来衡量声音波动变化的参数 是指声卡在采集和 播放声音文件时所使用数字声音信号的二进制位数 采样频率是指录音设备在一秒钟 内对声音信号的采样次数 采样频率越高声音的还原就越真实越自然 采样位数和采样率对于音频接口来说是最为重要的两个指标 也是选择音频接 口的两个重要标准 无论采样频率如何 理论上来说采样的位数决定了音频数据最 大的力度范围 2 32 3 采样定理采样定理 在进行模拟 数字信号的转换过程中 当采样频率最大值大于信号中最高频率 fmax 的 2 56 倍时 即 fs max 2 56fmax 则采样之后的数字信号完整地保留了原始信号 中的信息 2 42 4 时域信号的时域信号的 FFTFFT 分析分析 FFT 即为快速傅氏变换 是离散傅氏变换的快速算法 它是根据离散傅氏变换 的奇 偶 虚 实等特性 对离散傅立叶变换的算法进行改进获得的 在 MATLAB 的信号处理工具箱中函数 FFT 和 IFFT 用于快速傅立叶变换和逆变换 精品文档 函数 FFT 用于序列快速傅立叶变换 其调用格式为 y fft x 其中 x 是序列 y 是序 列的 FFT x 可以为一向量或矩阵 若 x 为一向量 y 是 x 的 FFT 且和 x 相同长度 若 x 为一矩阵 则 y 是对矩阵的每一列向量进行 FFT 如果 x 长度是 2 的幂次方 函数 fft 执行高速基 2FFT 算法 否则 fft 执行一种混合基的离散傅立叶变换算法 计算 速度较慢 函数 FFT 的另一种调用格式为 y fft x N 式中 x y 意义同前 N 为正 整数 函数执行 N 点的 FFT 若 x 为向量且长度小于 N 则函数将 x 补零至长度 N 若 向量 x 的长度大于 N 则函数截短 x 使之长度为 N 若 x 为矩阵 按相同方法对 x 进行 处理 2 52 5 数字滤波器的设计原理和方法数字滤波器的设计原理和方法 数字滤波是数字信号分析中最重要的组成部分之一 与模拟滤波相比 它具有精 度和稳定性高 系统函数容易改变 灵活性强 便于大规模集成和可实现多维滤波等 优点 在信号的过滤 检测和参数的估计等方面 经典数字滤波器是使用最广泛的一 种线性系统 数字滤波器的作用是利用离散时间系统的特性对输入信号波形 或频谱 进行加工处理 或者说利用数字方法按预定的要求对信号进行变换 2 62 6 各种滤波器的性能比较各种滤波器的性能比较 IIR 滤波器与 FIR 滤波器的性能比较 FIR Finite Impulse response 有限冲击响应 IIR Infinite Impulse response 无限冲击响应 从性能上来说 IIR 滤波器传输函数的极点可位于单位圆内的任何地方 因此可用 较低的阶数获得高的选择性 所用的存贮单元少 所以经济而效率高 但是这个高效 率是以相位的非线性为代价的 选择性越好 则相位非线性越严重 相反 FIR 滤波器 却可以得到严格的线性相位 然而由于 FIR 滤波器传输函数的极点固定在原点 所以 只能用较高的阶数达到高的选择性 对于同样的滤波器设计指标 FIR 滤波器所要求的 阶数可以比 IIR 滤波器高 5 10 倍 结果 成本较高 信号延时也较大 如果按相同的 选择性和相同的线性要求来说 则 IIR 滤波器就必须加全通网络进行相位较正 同样 要大增加滤波器的节数和复杂性 整体来看 IIR 滤波器达到同样效果阶数少 延迟小 但是有稳定性问题 非线性 相位 FIR 滤波器没有稳定性问题 线性相位 但阶数多 延迟大 3 3 设计内容设计内容 精品文档 3 13 1 语音信号的录入与提取语音信号的录入与提取 用电脑所带的录音工具录制一段录音 时间约为 1s 并将声音保存成扩展名为 wav 格式 并且应注意 录取的声音应为单声道的 否则图形会出不来 对采集的信号进行 FFT 分析 并且画出源信号的时域波形图形和频域图形 fs 16000 x0 wavread E xinhao1 t 0 length x0 1 16000 f fs 0 10239 20480 y1 fft x0 20480 画语音信号时域图 figure 1 subplot 1 1 1 plot t x0 grid on axis tight title 原始语音信号时域 xlabel time s ylabel 幅度 画语音信号频域图 figure 2 subplot 1 1 1 plot f abs y1 1 10240 grid on axis 0 8000 0 20 title 原始语音信号频谱 xlabel Hz ylabel 幅度 程序运行结果如下 精品文档 3 23 2 加噪处理 对信号加高频噪音加噪处理 对信号加高频噪音 在 MATLAB 软件平台下 给原始的语音信号叠加上噪声 绘出加噪声后的语音信号 时域和频谱图 在视觉上与原始语音信号图形对比 也可通过 Windows 播放软件从听 觉上进行对比 分析并体会含噪语音信号频谱和时域波形的改变 程序代码如下 加噪后语音信号变化 d 0 15 cos 2 pi 4000 t 加噪 x1 x0 d y2 fft x1 20480 精品文档 加噪后语音信号的时域图 figure 3 subplot 1 1 1 plot t x0 grid on axis tight title 加噪后的信号时域 xlabel time s ylabel 幅度 加噪后语音信号的频域图 figure 4 subplot 1 1 1 plot f abs y2 1 10240 grid on axis 0 8000 0 20 title 加噪语音信号频谱 xlabel Hz ylabel 幅度 程序运行结果如下 精品文档 加噪声后 听起来明显感觉嘈杂 通过图形比较也可以看出 时域波形图的变化 不大 但频谱图的变化明显 显示出加了干扰 3 33 3 设计窗函数带阻滤波器设计窗函数带阻滤波器 要对加入噪声的语音信号进行滤波 需设计合适的滤波器 由图形可看出 用窗 函数法设计带阻滤波器可对对加噪语音信号进行滤波 设计带阻滤波器 并绘制滤波 器的图形 程序设计如下 精品文档 设计带阻滤波器 wlp 2850 pi 2 fs wls 3400 pi 2 fs wup 4900 pi 2 fs wus 4600 pi 2 fs B wus wls M ceil 12 pi B 1 wp wls wlp 2 pi wus wup 2 pi bz az fir1 M wp stop blackman M 2 hn fir1 M wp stop blackman M 2 带阻滤波器窗函数图像 figure 5 subplot 1 1 1 freqz hn grid on axis tight title 窗函数图像 xlabel Hz ylabel 幅度 程序运行结果如下 3 43 4 对信号进行滤波处理对信号进行滤波处理 精品文档 对加入噪声的语音信号进行滤波 用加窗函数的方法设计了巴带阻滤波器对加噪 语音信号进行滤波 并绘制了加噪语音信号滤波后的时域图和频谱图 程序设计如下 画滤波后的图 figure 6 freqz bz az grid on axis tight y3 filter bz az x1 X fft y3 20480 滤波后时域图 subplot 2 1 1 plot y3 grid on ylim 0 03 0 03 title 滤波后信号时域波形 xlabel time s ylabel 幅度 滤波后频谱 subplot 2 1 2 plot abs X 1 10240 title 滤波后信号频谱 grid on ylim 0 20 xlabel Hz ylabel 幅度 程序运行结果如下 精品文档 经过滤波后得出的语音信号相比较 低通滤波后 声音稍微有些发闷 但是很接 近原来的声音 3 53 5 对信号加低频噪音对信号加低频噪音 程序如下 加低频噪音 d1 0 5 cos 2 pi 200 t 加噪 x4 x0 d1 y4 fft x4 20480 加噪后语音信号的时域图 figure 7 subplot 1 1 1 plot t x4 grid on axis tight title 加噪后的信号时域 xlabel time s ylabel 幅度 加噪后语音信号的频域图 精品文档 figure 8 subplot 1 1 1 plot f abs y4 1 10240 grid on axis 0 8000 0 20 title 加噪语音信号频谱 xlabel Hz ylabel 幅度 运行结果 3 63 6 对加高频噪音的程序进行滤波处理对加高频噪音的程序进行滤波处理 精品文档 用椭圆滤波器对信号的高频噪音进行滤波 程序如下 加高通滤波 wp1 100 pi 2 fs ws1 60 pi 2 fs Rp1 1 As1 100 N1 wso ellipord wp1 ws1 Rp1 As1 B1 A1 ellip N1 Rp1 As1 wso high figure 9 freqz B1 A1 grid on ylim 0 1 y4 filter B1 A1 x4 X1 fft y4 20480 滤波后时域图 subplot 2 1 1 plot y4 grid on axis tight ylim 0 03 0 03 title 滤波后信号时域波形 xlabel time s ylabel 幅度 滤波后频谱 subplot 2 1 2 plot abs X1 1 10240 title 滤波后信号频谱 grid on ylim 0 20 axis tight xlabel Hz 精品文档 ylabel 幅度 程序运行结果 椭圆高通滤波器图像 滤波器图像 figure 10 hn1 ellip N1 Rp1 As1 wso high subplot 1 1 1 freqz hn1 grid on