基于matlab的频谱分析与信号去噪

《基于matlab的频谱分析与信号去噪》由会员分享，可在线阅读，更多相关《基于matlab的频谱分析与信号去噪（27页珍藏版）》请在金锄头文库上搜索。

1、学生毕业设计报告基于 MATLAB 的频谱分析与信号去噪作者 系（院） 专业 年级 学号 指导教师 日期 学生诚信承诺书本人郑重承诺：所呈交的设计报告是我个人在导师指导下进行的研究工作及取得的研究成果。尽我所知，除了文中特别加以标注和致谢的地方外，报告中不包含其他人已经发表或撰写的研究成果，也不包含为获得安阳师范学院或其他教育机构的学位或证书所使用过的材料。与我一同工作的同志对本研究所做的任何贡献均已在报告中作了明确的说明并表示了谢意。签名： 日期：报告使用授权说明本人完全了解 有关保留、使用学位报告的规定，即：学校有权保留送交报告的复印件，允许报告被查阅和借阅；学校可以公布报告的全部或部分内

2、容，可以采用影印、缩印或其他复制手段保存报告。签名：导师签名：日期：第 1 页基于 MATLAB 的频谱分析与信号去噪摘要：本课题是基于数字信号处理的理论知识对语音信号、图像信号等的频谱分析以及对加噪声语音信号进行时域、频域分析和滤波设计。然后利用 MATLAB 软件进行编程，调试并完善程序，最终在计算机上得以实现。滤波器设计在数字信号处理中占有极其重要的地位，数字滤波器又有 FIR 数字滤波器和 IIR 滤波器两种。利用 MATLAB 信号处理工具箱可以快速有效地设计各种数字滤波器。在设计实现的过程中，使用窗函数法来设计 FIR数字滤波器，用巴特沃斯、切比雪夫和双线性变法设计 IIR 数字滤

3、波器，并利用 MATLAB 作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用 MATLAB 信号处理工具箱可以有效快捷地设计 FIR 和 IIR 数字滤波器，过程简单方便，结果的各项性能指标均符合指定要求。关键词 频域分析， 滤波器 ， MATLAB1 引言1.1 课题背景随着信息时代和数字世界的到来，数字信号处理己成为当今一门极其重要的学科和技术领域，数字信号处理在通信、语音、图像、自动控制、医疗和家用电器等众多领域得到了广泛的应用。任意一个信号都具有时域与频域特性，信号的频谱完全代表了信号，因而研究信号的频谱就等于研究信号本身。通常从频域角度对信

4、号进行分析与处理，容易对信号的特性获得深入的了解。因此，信号的频谱分析是数字信号处理技术中的一种较为重要的工具。数字滤波器, 是数字信号处理中及其重要的一部分。随着数字技术的发展，受到人们越来越多的重视。数字滤波器可以通过数值运算实现滤波，所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。数字滤波器种类很多，根据其实现的网络结构或者其冲激响应函数的时域特性，可分为两种，即有限冲激响应( FIR，Finite Impulse Response)滤波器和无限冲激响应( IIR，Infinite Impulse Response)滤波器。

5、在工程领域中，MATLAB 是一种倍受程序开发人员青睐的语言，对于一些需要做大量数据运算处理的复杂应用以及某些复杂的频谱分析算法 MATLAB 显得游刃有余。1.2 研究意义信号处理几乎涉及到所有的工程技术领域，而频谱分析又是信号处理中一个非常重要的分析手段。一般的频谱分析都依靠传统频谱分析仪来完成，价格昂贵，体积庞大，第 2 页不便于工程技术人员的携带。虚拟频谱分析仪改变了原有频谱分析仪的整体设计思路，用软件代替了硬件，使工程技术人员可以用一部笔记本电脑到现场就可轻松完成信号的采集、处理及频谱分析。 信 号 去 噪 被 用 于 从 一 堆 波 音 资 料 中 提 取 有 用 信 息 ， 去

6、除 干 扰 ， 提 高 波 音 资 料 信 噪比 。 为 了 提 高 信 噪 比 ， 人 们 根 据 信 号 和 噪 声 的 各 种 特 征 差 异 ， 设 计 了 许 多 去 噪 方 法 ，并 在 应 用 中 并 取 得 了 很 好 的 成 果 。 信 号 去 噪 的 很 多 方 法 都 是 利 用 短 时 傅 立 叶 变 换 来 滤波 去 噪 ， 但 是 短 时 傅 立 叶 变 换 不 能 同 时 兼 顾 时 间 分 辨 率 和 频 率 分 辨 率 。 用 不 同 种 滤波 器 滤 波 也 是 一 种 有 效 可 行 的 方 法 。1.3 本文研究内容信号的频谱分析就是利用傅里叶分析的方法

7、，求出与时域描述相对应的频域描述,从中找出信号频谱的变化规律，以达到特征提取的目的。不同信号的傅里叶分析理论与方法，在有关专业书中都有介绍。但实际的待分析信号一般没有解析式，直接利用公式进行傅里叶分析非常困难.。DFT 是一种时域和频域均离散化的傅里叶变换，适合数值计算且有快速算法，是分析信号的有力工具。DFT 及 FFT 是数字信号处理的重要内容。 DFT 是 FFT 的基础，FFT 是 DFT 的快速算法，在 MAT LAB 中可以利用函数 FFT 来计算序列的离散傅里叶变换 DFT。基于此首先介绍了 MAT LAB 软件；其次给出了基于 MAT LAB 软件的 DFT 和 FFT 频谱分

8、析的方法，利用 MAT LAB 软件方法，使得设计方便、快捷，大大减轻了工作量；再而我们提取一段语音信号，加噪，再通过设计的不同滤波器进行滤波后的频谱分析。本文将重点介绍基于 MATLAB 的频谱分析设计，包括：(1) 音频信号频谱分析；(2) 图像信号频谱分析；(3) 离散信号 /序列频谱分析；(4)语音信号提取，分析和加噪；(5)设计 IIR 和 FIR 的各种滤波器；(6)用设计的滤波器进行滤波；(7)分析滤波前后信号特征。2 频谱分析技术及MATLAB简介2.1 时域抽样定理时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件：对于基带信号，信号抽样频率 大于等于 2 倍的信号最高

9、频率 ，即 。时域抽样是把连续信samf mf2samf号 变成适于数字系统处理的离散信号 Xk。对连续信号 以间隔 T 抽样，则可得Xt Xt到的离散序列为 ，如图 2-1 所示。kTtXk第 3 页x(t)xk tk0T2012图 2-1 连续信号抽样的离散序列若 ，则信号 与 Xk的频谱之间存在：kTtXk Xtj sam1(e)j()nT其中， 的频谱为 ，Xk的频谱为 。j(e)(j)je可见，信号时域抽样导致信号频谱的周期化。 (rad/s)为抽样角频率， 为2/samT1samfT抽样频率。数字角频率 与模拟角频率 的关系为：=T。2.2 离散傅立叶变换(DFT) 有限长序列 的

10、离散傅立叶变换(DFT)为)(nx 10,)()(10NnWnxDFTkXNk逆变换为 ,)()()(10kXkInxNnn2.3 快速傅立叶变换(FFT)在各种信号序列中，有限长序列占重要地位。对有限长序列可以利用离散傅立叶变换(DFT) 进行分析。DFT 不但可以很好的反映序列的频谱特性，而且易于用快速算法(FFT)在计算机上进行分析。有限长序列的 DFT 是其 z 变换在单位圆上的等距离采样，或者说是序列傅立叶的等距离采样，因此可以用于序列的谱分析。FFT 是 DFT 的一种快速算法，它是对变换式进行一次次分解，使其成为若干小数据点的组合，从而减少运算量。MATLAB 为计算数据的离散快

11、速傅立叶变换，提供了一系列丰富的数学函数，主要有 fft、ifft、fft2 、ifft2, fftn、ifftn 和 fftshift、ifftshift 等。当所处理的数据的长度为 2 的幂次时，采用基-2 算法进行计算，计算速度会显著增加。所以，要尽可能使所要处理的数据长度为 2 的幂次或者用添零的方式来添补数据使之成为 2 的幂次。fft 函数调用方式：(1)Y=fft(X)(2)Y=fft(X,N)第 4 页(3)Yfft(X,dim)或 Yfft(X,N,dim) 。函数 ifft 的参数应用与函数 fft 完全相同。2.4 频谱分析原理时域分析只能反映信号的幅值随时间的变化情况，

12、除单频率分量的简单波形外，很难明确提示信号的频率组成和各频率分量大小，而频谱分析能很好的解决此问题。由于从频域能获得的主要是频率信息，所以本节主要介绍频率(周期)的估计与频谱图的生成。2.4.1 频率和周期的估计对于 Y(kf)，如果当 kf = f 时，Y(k f)取最大值，则 f 为频率的估计值，由于采样间隔的误差，f 也存在误差，其误差最大为 f / 2。周期 T=1/f。从原理上可以看出，如果在标准信号中混有噪声，用上述方法仍能够精确地估计出原标准信号的频率和周期。2.4.2 频谱图为了直观地表示信号的频率特性，工程上常常将 Fourier 变换的结果用图形的方式表示，即频谱图。以频率

13、 f 为横坐标，|Y( f)|为纵坐标，可以得到幅值谱；以频率 f 为横坐标，Arg Y(f)为纵坐标，可以得到相位谱；以频率 f 为横坐标，Re Y(f)为纵坐标，可以得到实频谱；以频率 f 为横坐标，Im Y(f)为纵坐标，可以得到虚频谱。根据采样定理，只有频率不超过 Fs/2 的信号才能被正确采集，即 Fourier 变换的结果中频率大于 Fs/2 的部分是不正确的部分，故不在频谱图中显示。即横坐标 f 0, Fs/22.5 MATLAB简介2.5.1 MATLAB软件的发展MATLAB 软件是由美国 Mathworks 公司推出的用于数值计算和图形处理的科学计算系统环境。MATLAB

14、是英文 MATrix LABoratory(矩阵实验室)的缩写。它的第一版(DOS 版本 1.0)发行于 1984 年，经过 10 余年的不断改进，现今已推出它的 Windows 98/NT 版本(6.1 版)。新的版本集中了日常数字处理中的各种功能，包括高效的数值计算、矩阵运算、信号处理和图形生成等功能。在 MATLAB 环境下，用户可以集成地进行程序设计、数值计算、图形绘制、输入输出、文件管理等各项操作。MATLAB 提供了一个人机交互的数学系统环境，该系统的基本数据结构是矩阵，在生成矩阵对象时，不要求作明确的维数说明。与利用 C 语言或 FORTRAN 语言作数值计算的程序设计相比，利用

15、 MATLAB 可以节省大量的编程时间。在美国的一些大学里，MATLAB 正在成为对数值计算、算法预设计与验证，以及一些特殊的短阵计算应用，如自动控制理论、统计、数字信号处理( 时间序列分拆)等。MATLAB 系统最初是由 CieveMoler 用 FORTRAN 语言设计的，有关矩阵的算法来自 LINPACK 和 EISPACK 课题的研究成果；现在的 MATLAB 程序是 MathWorks 公司用 C 语言开发的，第一版由 steve Bangert 主持开发编译解释程序，Steve Kleiman 完成图形功能的设计，John Little 和 Cleve Moler 主持开发各类数学分析的子模块，撰写用户指南和大部分 M 文件。自从第 1 版发行以来，已第 5 页有众多的科技工作者加入到 MATLAB 的开发队伍中，并为形成今天的 MATLAB 系统做出了巨大的贡献，MATLAB 以商品形式出现后，仅短短几年，就以其良好的开放性和运行可靠性，使原先控制领域里的封闭式软件包(如英国的 UMIST，瑞典的 LUND，德国的 KEDDC)纷纷淘汰，而改以 MATLAB 为平台加以重建。在时间进入 20 世纪九十年代的时候，MATLAB 已经成为国际控制界公认的标准计算