参考基于matlabgui信号处理系统设计

目录目录I摘要VABSTRACTVII第1章 前言11.1 课题开发背景和发展状况11.2 研究的意义21.3 本文研究的主要内容3第2章 基本信号分析与处理52.1 信号分析概述52.2 滤波52.2.1 滤波器概述52.2.2 FIR滤波器和IIR滤波器的比较62.3 窗函数72.3.1 窗函数的分类72.3.2 常用的窗函数的性质和特点72.4 时域分析92.4.1 信号的强度表述92.4.2 信号的相关分析102.5 频域分析132.5.1 幅值谱142.5.2 功率谱142.5.3 倒频谱152.6 信号频谱的主要性质162.6.1 时移特性162.6.2 频移特性172.6.3 Parseval定理172.7 信号处理中出现的现象182.7.1 频率混叠现象182.7.2 频谱泄漏现象192.7.3 栅栏效应192.8 系统幅频特性202.8.1 频率响应函数202.8.2 二阶系统的幅频特性212.9 信号调制222.9.1 幅值调制222.9.2 频率调制232.9.3 相位调制242.10 本章小结25第3章 基于MATLAB GUI信号处理实验室的设计及实现273.1 MATLAB GUI273.2 系统总体设计303.3 系统首页的设计313.4 主界面的设计323.5 信号源选择模块的设计333.5.1 模拟信号信息设置模块的设计343.5.2 数据文件读取模块的设计343.6 文件操作模块的设计353.7 信号预处理的设计363.7.1 滤波器设计363.7.2 加窗处理的设计373.7.3 噪声处理的设计383.8 信号时域分析的设计393.8.1 信号参数测量的设计403.8.2 信号相关分析的设计403.9 信号频域分析的设计423.9.1 幅值谱分析的设计423.9.2 功率谱分析的设计433.9.3 倒频谱分析的设计443.10 信号频谱主要性质的设计443.10.1 信号时移特性的设计443.10.2 信号频移特性的设计453.10.3 Parseval定理的设计463.11 信号处理中出现的现象的设计463.11.1 频率混叠现象的设计463.11.2 频谱泄漏现象的设计473.11.3 栅栏效应的设计483.12 系统幅频特性的设计493.13 信号调制的设计503.13.1 幅值调制的设计503.13.2 频率调制的设计513.13.3 相位调制的设计523.14 mcc编辑介绍533.15 本章小结54第4章 结论55第5章 总结与体会57谢辞59参考文献61附录63附录 英文原文63附录 中文翻译71第 77 页 摘要数字信号处理DSP(Digital Signal Processing)被广泛地应用于通信工程、电子信息工程、雷达、遥感、生物工程等许多领域，许多高校的相关专业都把这门课程列为专业必修课。但是，该课程涉及到许多较深的数学理论知识，学生普遍感觉抽象难学。基于上述分析，本文对信号处理中常用的知识进行了讲述，主要讲述了信号的预处理、时域分析、频域分析、信号频谱的主要性质、信号处理中常出现的现象、系统幅频特性以及信号调制的相关知识。并以上述的信号处理知识为基础，利用MATLAB的GUIDE(图形用户界面设计向导)工具箱设计并实现了“基于MATLAB GUI的信号处理实验室”。该系统以模块化的设计思路，并以窗口、弹出框、按键、动态文本框等友好界面形式构成，实现了信号处理中常见知识的描述。即便不懂MATLAB语言编程，也可以非常方便地使用该系统，对信号处理的知识进行比较学习，提高学生变被动为主动的学习态度。关键词：数字信号处理；MATLAB；GUI；实验室ABSTRACTDigital Signal Processing DSP (Digital Signal Processing) is widely applied in communication engineering, electronic and information engineering, radar, remote sensing, biological engineering, and many other areas, many universities are related to the course as professional courses. However, this course involves much deep mathematical theory knowledge; the student feel abstraction is difficult to learn.Based on the above analysis, this paper used in signal processing, the knowledge about the signal pretreatment, the analysis of time and frequency domain analysis, the main character, signal spectrum signal processing in the phenomenon, amplitude frequency characteristics and related knowledge of modulation signal. And on the basis of signal processing knowledge, GUIDE the MATLAB graphical user interface design (GUIDE) was designed and implemented toolbox based on MATLAB GUI signal processing laboratories. This system with modular design train of thought, and with Windows, pop-up, buttons, dynamic text box etc, friendly interface, the form of the common knowledge of signal processing, if do not understand, can also MATLAB programming language is very convenient to use this system, signal processing knowledge to improve students learning, compared to active learning attitude.Keyword: Digital Signal Processing; MATLAB; GUI; Laboratory前言课题开发背景和发展状况数字信号处理DSP(Digital Signal Processing)已经被广泛地应用于通信工程、电子信息工程、雷达、遥感、生物工程等许多领域，许多高校的相关专业都把这门课程列为专业必修课Error! Reference source not found.。数字信号处理的目的是对真实世界的连续模拟信号进行测量或滤波。数字信号处理的核心算法是离散傅立叶变换(DFT)，DFT使信号在时域和频域都实现了离散化，从而可以用通用计算机处理离散信号。而使数字信号处理从理论走向实用的是快速傅立叶变换(FFT)，FFT的出现大大减少了DFT的运算量，使实时的数字信号处理成为可能、极大促进了该学科的发展。信号处理是机械设备故障诊断中的关键技术，起着对监测信号进行加工、变换，从中提取出敏感故障征兆的重要作用。机械系统在工作中会产生大量振动信号，如何对这些信号进行有效的分析是该领域急需解决的关键问题之一。数字信号处理在当今高科技领域有着极为重要的地位和广泛的用途，然而它的基础课程却处于教难、难学的境况中Error! Reference source not found.。主要因为：数字信号处理本身具有抽象性。一方面它是基于“高等数学”、“信号分析”等公式和推导较多的学科，公式、性质的推导和证明繁琐，理论性很强；另一方面数字信号处理是适应高速数字集成电路的，其大量的计算算法适于在计算机上实现，对于人来说则运算量大且繁琐，学习者难以亲手验证，因而经常得不到形象化的结果，使得对理论的理解难以透彻，实际应用中总有一层障碍。因此数字信号处理被公认为大学课程中最难的课程之一，学生对上课内容的理解往往跟不上老师授课的进度。现在很多大学教师在课堂上采用了多媒体教学，多媒体授课具有教学信息量大，信息传递便捷、交互性较强、感官刺激效果好等特点Error! Reference source not found.。但是，由于多媒体教学过程中，学生只能根据老师给出的实例对信号处理的知识进行学习，使得学生学习起来还是比较枯燥，这就需要一个可以更系统的展现信号处理的知识，又能够根据学生的需求，把信号处理知识更直观、形象的展现出来，有方便友好的人机交互界面，消除教与学的隔阂，提高变被动为主动的学习态度。Matlab是Mathworks公司于1982年推出的一套高性能的数值计算和可视化软件，经过20多年的发展和各个领域的专家的共同努力，如今它已集数值分析、矩阵运算、信号处理和图形显示于一体，形成了用于各个领域的工具箱函数，诸如模糊集、神经网络、小波分析、控制系统、图像处理、系统辩识、信号处理等，其强大的扩展功能为各个领域的应用提供了基础，被誉为“站在巨人肩上的工具”。Matlab 7.0引入了Matlab桌面，可用来管理各种窗口，有助于提高工作效率，消除了各种版本中文本编辑器的差异，提高运算速度。GUI制作工具得到了极大的改进和扩展补充，使得图像制作更加简便，图像质量也有很大提高。由于Matlab在数字信号处理领域有着巨大的优势，比如产生波形、抽样、DFT与FFT变换、卷积、滤波器设计等都可以很容易地实现，特别是随着SPTOOL的推出，这种应用变得更为广泛。以MATLAB GUI为开发平台开发一个信号处理实验系统作为辅助教学方法在很多大学得到广泛应用，但是在这些系统中，一般只能实现信号分析与处理的某一项或几项功能(例如，一些系统只能实现信号的滤波处理或加窗处理)；有些系统虽然能实现常见的信号的分析与处理方法，但在信号分析处理过程中遇到的现象则没有进行具体的描述；而且大多数的系统只能对系统本身自带的模拟信号进行分析，而对用户自己得到的信号则无能为力。因此，需要一个能对信号的基本分析方法、常见的现象以及用户得到的信号进行分析显示的系统来满足用户需求。在上述的背景下提出了本课题。研究的意义本课程通过利用MATLAB作为分析数据信号的工具，以MATLAB为开发平台编写GUI程序，实现基本的信号分析与处理。该GUI以窗口、弹出框、按键、动态文本框等友好界面形式构成，并实现了信号处理中常见知识的描述，即便不懂MATLAB语言编程，也可以非常方便地使用该GUI，对信号处理的知识进行比较学习。通过MATLAB GUI开发的信号处理实验室可以对信号进行分析处理，得出信号的基本信息。在机械故障诊断过程中，可以根据所得到信号的基本信息，用这些信息与MATLAB GUI中得到的信息进行比较，可以得出机械故障诊断的一些故障信息。该系统还可以用于学生学习信号处理的一种工具，用该工具可以更直观的体现信号的一些基本信息，这样就可以让学生可以更好的学习这些信号处理的基本信息。因此用MATLAB GUI为平台开发的“信号处理实验室”可以广泛的推广，更多的用于学习、比较不同的信号处理的方法，从而达到让学习信号处理的初学者可以更轻松、更形象的学习信号处理的基本分析方法，这样就可以在以后对机械故障诊断中对分析出来的信号进行正确理解，更好的了解机械的工作情况，从而达到对机械设备进行监测的目的。本文研究的主要内容本文旨在运用MATLAB GUI处理平台，设计开发一个信号处理实验室系统，使得该系统可以对模拟信号、重载信号进行分析、显示以及存储功能。第一章首先介绍了课题开发背景、发展现状和研究意义，分析了信号处理实验室系统在教学中应用的优势，阐述了课题设计的必要性。第二章主要介绍了信号处理实验室的具体设计阶段所需的有关信号分析处理方面的必备知识。第三章是全文的重点部分，详细论述了在MATLAB GUI平台下构建信号处理实验室系统的软件设计过程，详细介绍了程序中的各个模块设计过程。包括信号分析(信号预处理、时域分析、频域分析、频谱分析的主要性质、信号处理中出现的现象、系统幅频特性、信号调制)、信号发生、数据保存和重载等多个功能模块的设计。其中，信号预处理包括滤波处理、加窗处理、添加噪声处理；时域分析包括显示波形，信号电压、频率等参数测量，信号的相关分析；频域分析包括幅值谱、相位谱、功率谱等；频谱分析的主要性质包括时移特性、频移特性以及Parseval定理；信号处理中出现的现象包括频率混叠现象、频谱泄漏现象、栅栏效应；信号调制包括幅值调制、频率调制、相位调制；通过信号发生模块可以产生常用信号(如正弦波、余弦波、方波、三角波)。第四章是结论，对课题的工作进行了总结，提出尚待开展的工作。第五章是总结与体会，总结在完成课题之后的一些学习、生活方面的感触。基本信号分析与处理信号分析概述信号是来自于自然界与工程界的各种物理现象，作为载体，信号包含着反映被测系统的状态或特征的某些有用信息，它是人们认识客观事物内在规律，研究事物之间的相互关系、预测未来发展的依据。工程信号的最初也是最基本的形式是物理量随时间的变化，即时间信号，或称时间历程。信号分析就是拾取信号，对其进行计算、分析和处理，从中提取信息的过程。信息是一个带有主观意识的概念。要了解物理过程中存在的各种内在特征，就需要从不同的角度来提取有关这些特征的信息，也既要用不同的技术和方法来处理信号；突出、分离出相关的信息，衰减、摒弃其他成分。这就决定了信号分析原理，技术以及方法的多样性Error! Reference source not found.。滤波滤波器概述滤波器是信号处理不可缺少的一部分，滤波器顾名思义其作用的对输入信号起到滤波的作用。其时域输入输出关系是： (11)若、的傅立叶变换存在，则输入输出的频域关系是： (12)滤波器的种类很多，分类方法也不同，可以从功能上分，也可以从实现方法上分，或从设计方法上来分类。但总的来说，滤波器可分为两类 ，即经典滤波器和现代滤波器。经典滤波器是假定输入信号 中的有用成分和无用成分(如噪声)各自占有不同的频带，当通过滤波器后可将无用成分有效地去除。如果信号中的有用成分和无用成分的频谱相互重叠，那么经典滤波器将无法滤除信号中的无用成分。经典滤波器可以分为四种：即低通(LP, Low Pass)、高通(HP, High Pass)、带通(BP, Band Pass)和带阻(BS, Band Pass)滤波器，每一种又有模拟滤波器(AF)和数字滤波器(DF)两种形式。数字滤波器按照单位取样响应的时域特征可分为无限脉冲响应(IIR，Infinite Impulse Response)系统和有限脉冲响应(FIR，Finite Impulse Response )系统。如果单位取样响应是时宽无限的，则称为IIR 系统；如果单位取样响应是时宽有限的， ，则称为FIR 系统。现代滤波器理论研究的主要内容是从含有噪声的数据记录(又称为时间序列)中估计出信号的某些特征或信号本身。一旦信号被估计出来，那么估计出的信号与原信号相比会有高的信噪比。现代滤波器把信号和噪声都视为随机信号，然后利用它们的统计特征(如自相关函数、功率谱等)推导出一套最佳的估值算法，然后用硬件或软件来实现Error! Reference source not found.。FIR滤波器和IIR滤波器的比较从性能上说，IIR 滤波器系统函数的极点可以位于单位圆内的任何位置，因此在相同设计指标下(与FIR 滤波器相比)，实现IIR 滤波器的阶次低，即所用的存储单元少，从而经济效率高。但是这个高效率是以相位的非线性为代价的；相反，FIR 却可以得到严格的线性相位，但是FIR 滤波器系统函数的极点固定在原点。所以只能通过较高的滤波器阶次实现相同设计指标(与IIR 滤波器相比)的滤波器设计。对于同样的滤波器设计指标，FIR 滤波器的阶次比IIR 滤波器高510倍，结果成本较高，信号延时也较大。从结构上看，IIR 滤波器必须采用递归结构，极点位置必须在单位圆内，否则系统将不稳定。要这种结构中，由于运算过程中对序列进行舍入处理，这种有限字长效应有时会产生寄生振荡；相反，FIR 滤波器采用非递归结构，不存在稳定的问题，同时运算误差也较小。此外，FIR 滤波器可以采用快速傅里叶变换算法，在相同滤波器阶次的条件下，FIR 滤波器的运算速度比IIR 滤波器快得多。从设计工具看，IIR 滤波器可以借助模拟滤波器的成果，因此一般都有有效的封闭形式的设计公式可供参考，计算工作量较小，而且对计算工具的要求不高；FIR 滤波器一般没有封闭形式的设计公式。窗函数设计FIR 滤波器也仅给出了窗函数的计算公式，但是在计算通带、阻带衰减时无显式表达式。一般FIR 滤波器的设计只有计算程序可循，因此它对计算工具要求较高。另外，IIR 滤波器虽然设计简单，但主要用于设计具有设计片段常数特性的滤波器，如低通、高通、带通、带阻等，往往脱离不了模拟滤波器的格局。而FIR滤波器则要灵活得多，易于实现某些特殊的应用。综上所述，本系统采用IIR 滤波器来设计滤波器，Matlab工具箱提供了几个直接设计数字滤波器的函数，它们使IIR 数字滤波器的设计变得简单。窗函数为了减少频谱能量泄漏，可以采用不同的截取函数进行截断，截断函数成为窗函数，简称为窗。窗函数的分类实际应用的窗函数，可以分为以下主要类型幂窗采用时间变量某种幂次的函数，如矩形、三角形、梯形或其他时间(t)的高次幂；三角函数窗应用三角函数，即正弦或余弦函数等组合成复合函数，例如汉宁窗、海明窗等；指数窗采用指数时间函数，如形式，例如高斯窗等常用的窗函数的性质和特点矩形窗矩形窗属于时间变量的零次幂窗，函数形式为 (13)相应的窗谱为： (14)矩形窗使用最多，习惯上不加窗就是使用了矩形窗。这种窗的优点是主瓣比较集中，缺点就是旁瓣较高，并有负旁瓣，导致变换中带进了高频干扰和泄漏，甚至出现负谱现象。三角窗三角窗又称为费杰(Fejer)窗，是幂窗的一次方形式，其定义为; (15)相应的窗谱为： (16)三角窗与矩形窗比较，主瓣宽约等于矩形窗的两倍，但旁瓣小，而且无负旁瓣，如图所示：汉宁(hanning)窗汉宁窗又称为升余弦窗，其时域表达式为： (17)相应的窗谱为; (18)由此式可以看出，汉宁窗可以看作是3个矩形时间窗的频谱之和，或者说是3个型函数之和，而括号中的两项相对于第一个谱窗向左、右各移动了，从而使旁瓣相互抵消，消去高频干扰和漏能。海明(hamming)窗海明窗也是余弦窗的一种，又称改进的升余弦窗，其时间函数表达式为： (19)相应的窗谱为： (110)海明窗与汉宁窗都是余弦窗，只是加权系数不同。海明窗加权的系数能使旁瓣达到更小。分析表明，海明窗的第一旁瓣衰减为-42dB。海明窗的频谱也是由3个矩形时间窗的频谱合成，使其旁瓣衰减速度为20dB/(10oct)，这比汉宁窗衰减速度慢。海明窗与汉宁窗都是很有用的窗函数。布莱克曼(Blackman)窗布莱克曼窗的时间函数表达式为： (111)布莱克曼窗主瓣宽，旁瓣小，频率识别精度最低，但幅值识别精度最高。时域分析在时域分析模块中，主要有信号时域波形的描述、信号的有效值、峰值、峰峰值的求解以及信号的相关分析。信号的强度表述信号的强度以峰值、绝对均值、有效值和平均功率来表述(见图2-1)。周期信号的强度表示峰值是信号可能出现的最大瞬时值，即 (112)峰峰值是在一个周期中最大瞬时值和最小瞬时值之差。对信号的峰值和峰峰值应有足够的估计，以便确定测试系统的动态范围。一般希望信号的峰峰值在测试系统的线性区域内，使所观测(记录)到的信号正比于被测量的变化状态。如果进入非线性区域，则信号将发生畸变，结果不但不能正比于被测信号的幅值，而且会增生大量谐波。周期信号的均值为 (113)它是信号的常值分量。周期信号全波整流后的均值就是信号的绝对均值，即 (114)有效值是信号的均方根，即 (115)有效值的平方均方值就是信号的平均功率，即 (116)它反映信号的功率大小。信号的相关分析相关函数的定义和物理意义设和是两个均值为零(下同)的随机函数，若二者波形完全相同，则存在如下线性关系： (117)式(2-17)中，为常数，该式表示和线性相关，波形完全相似，相互间无时差。在一般情况下，式(2-17)是不可能满足的。为了判别和的波形相似程度，我们研究它们之间的差函数 (118)的平均功率，亦即误差功率可写为 (119)的大小随值变化，适当选择的值，可以使为极小，即 (120)于是： (121)将式(2-21)带入式(2-19)，得出衡量与差异的指标的极小值 (122)式中 (123)由式(2-22)可知，的值决定了两信号差函数的功率的大小。若的绝对值等于，则误差功率等于零，信号与线性相关，波形相似。如果等于零，那么误差功率将等于信号的功率，两信号完全不相关。通常。再看式(2-23)，等式右边的分母分别是和的均方根值，的大小取决于分母上的积分值。改积分表示信号与线性相关程度，称为与的相关值。称为相关系数值。互相关函数实际上，两个信号之间可能有时差，因而需要研究的是信号与的时延信号的线性相关和波形相似程度。很显然，这种相关程度是时延的函数。的量纲和相同，均为秒。定义两个信号与的互相关函数和互相关系数分别为 (124) (125)前述的相关函数值和相关系数值是在时的相关函数和相关系数的值。因此，我们可以用相关函数或相关系数作为一把尺子来判断两个信号的相似性。当变化时，若出现极大值，就表示在此时延下，两个信号具有一定的线性相关性。相关程度的大小由的值决定，越接近，两信号的线性相关性也越好。相关函数和相关系数随的变化规律完全相同，二者只差一个常数，即式(2-25)的分母表示的与的均方根之和。实际上，作相关分析时，常常计算，因为它具有归一化和可比性；而的表达式比较简单，讨论问题时往往用相关函数，有时也把它们都称为相关函数而不加区分。必须指出，只有当样本长度时，式(2-24)和(2-25)才是随机信号与的互相关函数和互相关系数的定义式。但实际分析处理数据时，样本长度只能是有限的。上述式子是所谓的估计式。对于周期信号，T应取信号的周期。而对瞬变非周期能量信号，它们的互相关函数和互相关系数的定义稍有不同。 (126) (127)互相关函数的性质如下：两信号是同频率的周期信号或包含有同频率周期成分，才有互相关函数，即同频相关，不同频不相关。两个相同周期的信号的互相关函数仍是周期函数，其周期与原信号的周期相同，并保留了两个信号的幅值和相位差信息。两信号在相隔一时间间隔处，可能有最大值，它反映了和之间主传输通道的滞后时间。不是偶函数，即。自相关函数当等于时，自相关函数和系数描述的是信号的一个时刻的值与另一个时刻的值之间的线性相关程度，我们称它们为的自相关函数和自相关系数，用和表示，即 (128) (129)与互相关函数类似，对随机信号应有，对周期信号应取信号的周期，而对瞬变能量信号，自相关函数和系数分别是： (130) (131)频域分析在机械故障诊断中，时域分析所能提供的信息量非常有限。时域分析往往只能粗略地回答机械设备是否有故障，有时也能得到故障严重的程度，但不能回答故障发生部位等信息，因此在故障诊断中还应该结合频域分析来进行诊断。幅值谱在信号的频域描述中，以频率作为自变量，以组成信号的各个频率成分的幅值作为因变量，这样的频率函数称为幅值谱，它表征信号的幅值随频率的分布情况。傅立叶变换是平稳信号分析和处理的一个重要工具，一个时域的问题通过傅立叶变换能被转化成频域的问题来分析研究。设x(t)为t的函数，并且满足狄里赫利条件，其计算式如下： (132)连续傅立叶变换实现了测试信号从时域到频域的转换，在理论分析中具有很大的价值。但是连续傅立叶变换不能直接应用于计算机技术，繁琐的计算限制了它的进一步发展。直到离散傅立叶变换的出现，使得数学方法与计算机技术建立了联系，对工程实际来说它有更重要的价值。设x(n)为一时域数字序列，其计算式如下： (133)对于离散傅立叶变换，求出N点x(k)需要N2次复数乘法，N(N-1)次加法运算。每次复数乘法要做四次实数乘法，两次实数加法运算。这样当N很大时，所需的计算工作量也是相当可观的。为了进一步减少计算量，库利图基于1965年提出了快速傅立叶算法，使用快速傅立叶变换，计算N点离散傅立叶的计算量由N2次复数乘法降为Nlog2N次复数乘法，计算量显著减小。频谱能够分析信号的能量(或功率)的频率分布，通常能够提供比时域波形更加直观的特征信息。该方法理论上比较成熟，适合于平稳信号的分析。功率谱信号的幅值谱描述组成信号的谐波分量的幅值随频率的分布，功率谱则描述信号功率在频率域内的分布。周期信号的功率(也称均方差)可写为 (134)将积分内的一个用它的傅立叶级数的复指数函数展开式带入， (135)求和式(2-35)内的积分等于，最后得 (136)这说明，周期信号的功率等于它的双边幅值谱的平方和。随的变化表示信号功率随频率的分布，故是周期信号的功率谱。根据信号双边谱和单边谱的关系，还可以 (137)直流分量的功率等于它的幅值的平方；正弦分量的有效值是它的幅值的，功率是其有效值的平方，等于幅值平方的二分之一。式(2-37)的物理意义很明确，周期信号的功率是它的各分量之和。倒频谱已知时域信号经过傅立叶变换可变为频域函数或功率谱密度函数。当频谱图上呈现出复杂的周期结构时，如果再进行一次对数的功率谱密度函数傅立叶变换并取平方，则可得到倒频谱函数(Power Sepstrum),其数学表达式为： (138)倒频谱也可以表述为：“对数功率谱的功率谱”。式被定义为信号的功率倒频谱。工程上常用的是式(2-38)的开平方，即： (139)式中，被称为幅值倒频谱，有时简称倒频谱。倒频谱中自变量q称为倒频率，它具有与自相关函数中的自变量相同的时间量纲(因为倒频谱是傅立叶正变换，积分变量是频率而不是时间，故倒频谱的自变量具有时间的量纲)单位为s或ms。q值大者称为高倒频谱，表示谱图上的快速波动。q值小者称为低倒频谱，表示谱图上的缓慢波动。倒频谱是频域函数的傅立叶再变换。对功率谱函数取对数的目的，是使变换后的信号能量格外集中。同时还可解析卷积成分，易于对原信号的识别Error! Reference source not found.。信号频谱的主要性质时移特性信号和相比，在时间轴上超前或滞后了时刻，如图2-2所示。信号的时延那么时移特性可写为 (140)原信号的频谱 (141)时移信号的频谱 (142)比较上述两式可见，时移信号和原信号的幅值频谱相同，但相位频谱则增加了与频率成线性的附加项。如果频率以记，线性项的系数就表示时延量，正系数表示超前，负系数表示滞后。幅值谱不变，则时移信号的各谐波分量的幅值和原信号的相同；与频率成线性的附加相移，保证各谐波分量产生同一时延；从而，合成信号与原信号相比波形相同，并存在一时移。频移特性频移特性可表示为 (143)该式说明，信号乘复指数(常称复调制)后，信号的时域描述已大大改变，从实信号改变为复信号，但其频谱现状不变，只是在频域内移动了，如图2-3所示。信号的频移Parseval定理帕塞瓦尔定理指出，在时域中计算的总能量，等于在频域中计算的信号总能量。 (144)式又叫做能量等式。这个定理可以用傅立叶变换的卷积公式导出。称为能谱，它是沿频率轴的能量分布密度。在整个时间轴上信号平均功率为 (145)因此，自功率谱与幅值谱的关系为 (146)信号处理中出现的现象频率混叠现象如果采样的时间间隔Ts太大，既采样频率fs太低，平均距离1/Ts过小，那么移至各采样脉冲所在处的频谱X(f)就会有一部分交叠，新合成的X(f)*S(f)图形与原X(f)不一致，这种现象称为频谱混叠现象Error! Reference source not found.，如图2-4所示。频率混叠现象在时域采样过程中，若已知信号的最高频率为(即频谱分析范围时)，为了避免在DFT运算过程中发生频率混叠现象，要求抽样频率满足下式：也就是抽样间隔T满足一般应取如果不满足的要求，就产生频率响应的周期延拓分量互相重叠的现象，也就是产生频率响应的混叠失真。频谱泄漏现象数字信号处理的主要数学工具是博里叶变换而傅里叶变换是研究整个时间域和频率域的关系Error! Reference source not found.。不过，当运用计算机实现工程测试信号处理时，不可能对无限长的信号进行测量和运算，而是取其有限的时间片段进行分析。做法是从信号中截取一个时间片段，然后用观察的信号时间片段进行周期延拓处理，得到虚拟的无限长的信号，然后就可以对信号进行傅里叶变换、相关分析等数学处理。无线长的信号被截断以后，其频谱发生了畸变，原来集中在f(0)处的能量被分散到两个较宽的频带中去了(这种现象称之为频谱能量泄漏)。为了减少频谱能量泄漏，可采用不同的截取函数对信号进行截断，截断函数称为窗函数，简称为窗。有关窗的知识在2.3节中有详细介绍。信号截断以后产生的能量泄漏现象是必然的，因为窗函数w(t)是一个频带无限的函数，所以即使原信号x(t)是有限带宽信号，而在截断以后也必然成为无限带宽的函数，即信号在频域的能量与分布被扩展了。又从采样定理可知，无论采样频率多高，只要信号一经截断，就不可避免地引起混叠，因此信号截断必然导致一些误差。泄漏与窗函数频谱的两侧旁瓣有关，如果两侧瓣的高度趋于零，而使能量相对集中在主瓣，就可以较为接近于真实的频谱，为此，在时间域中可采用不同的窗函数来截断信号。栅栏效应因为DFT计算频谱只在限制为基频F0的整数倍处的谱，而不是连续的频率函数，这种只能在离散点的地方才能看到真实景象的现象就被称为“栅栏效应”。减小栅栏效应的一个办法就是使频域抽样更密，即增加频域抽样点数N，在不改变时域数据的情况下，必然是在时域数据末端添加一些零值点，使一个周期内的点数增加，但不改变原有的记录数据。频域抽样为，N增加，必然使样点间距离更近，谱线更密，栅栏效应就减小了。如果数据长度T01，抽样点数N1，补零后的数据长度T02，抽样点数N2，则因为N2N1，故F01F02，所以有人认为补零后，频率分辨力就提高了，这是错误的。原因是补零后，不能增加数据的有效长度，上例中的实际长度仍为T01，因而补零是不能提高频率分辨力的。补零的好处：(1)克服栅栏效应；(2)使N为2的整数幂值，便于FFT计算。减小栅栏效应的另一个办法就是使采样频率更小，在采样点数相等的情况下，减小采样频率可以使频率分辨率更小，从而达到减小栅栏效应的目的Error! Reference source not found.。系统幅频特性频率响应函数正弦信号可以用旋转矢量来表示。一个幅值为X，频率为，出相位为的正弦信号可写为 (147)根据频率保持特性，该信号通过线性系统的输出是 (148)将上述的和的表达式代入线性系统的表达通式 (149)可得 (150)式(2-50)等号右边的分母由系统参数和频率变量组成的复函数，定义为系统的频率响应函数，即 (151)的模称为系统的幅频特性，则 (152)的幅角称为系统的相频特性，则 (153)所以，系统的幅频特性表示系统在正弦输入输出时，输出与输入的幅值之比；相频特性表示输出与输入之间的相位差，它们是频率的函数。幅频、相频特性的组合就是系统的频率响应函数 (154)若输入信号，那么通过线性系统的输出就可以写为 (155)式中，和分别是系统的幅频、相频特性在时的值。二阶系统的幅频特性由二阶系统的定义可以得到二阶系统的传递函数为 (156)二阶系统的频率响应函数为 (157)相应的幅频特性和相频特性分别为 (158) (159)信号调制信号的调制是指利用某种低频信号来控制或改变一个高频振荡信号的某个参数(幅值、频率或相位)的过程。当被控制的量是高频振荡信号的幅值时，称为幅值调制或调幅；当被控制的量是高频振荡信号的频率时，称为频率调制或调频；当被控制的量是高频振荡信号的相位时，称为相位调制或调相。在这里，我们称高频振荡信号为载波，控制高频振荡信号的低频信号为调制信号，调制后的高频振荡信号为已调制信号。幅值调制幅值调制是将一个高频载波信号与被测信号(调制信号)相乘，使高频信号的幅值随被测信号的变化而变化。如下所示：为被测信号，为高频载波信号：，则已调制信号为与的乘积： (160)调幅信号的频域分析由傅里叶变换的性质可知：时域中两个信号相乘对应于频域中这两个信号的傅里叶变换的卷积，即 (161)余弦函数的频域波形是一对脉冲谱线，即 (162)由上面几式可得出 (163)一个函数与单位脉冲函数卷积的结果是将这个函数的波形由坐标原点平移至脉冲函数处。所以，把被测信号和载波信号相乘，其频域特征就是把的频谱由频率坐标原点平移至载波频率处，其幅值减半，如图2-5所示。可以看出所谓的调幅过程相当于频谱“搬移”过程。调幅信号的频谱从图2-5可以看出，载波频率必须高于信号中的最高频率，这样才能使已调幅信号保持原信号的频谱图形而不产生混叠现象。为了减小电路可能引起的失真，信号的频宽对于载波频率应越小越好。在实际应用中，载波频率常常在调制信号上限频率的十倍以上。频率调制频率调制是指利用调制信号控制高频载波信号频率变化的过程。在频率调制中载波幅值保持不变，仅载波的频率随调制信号的幅值成比例变化。记载波，这里角频率为一常量。如果保持振幅A为常数，让载波瞬时角频率随调制信号作线性变化，则有 (164)式中 k比例因子此时调频信号可表示为 (165)图2-6是调制信号为三角波时的调制信号波形。三角波调整信号下的调频波由图可见，在区间，调制信号，调频信号的频率保持原始的中心频率不变；在区间，调频波的瞬时频率随调制信号x(t)的增大而逐渐增高；在区间，调频波随调制信号的减小而逐渐降低；在后，调制信号，调制信号的频率又恢复了原始的中心频率。相位调制相位调制是指利用调制信号控制高频载波信号相位变化的过程。在频率调制中载波幅值保持不变，仅载波的相位随调制信号的幅值成比例变化。记载波，这里初相位为一常量。如果保持振幅A为常数，让载波瞬时相位随调制信号作线性变化，则有 (166)式中 k比例因子此时调频信号可表示为 (167)本章小结本章主要介绍了信号处理实验室的具体设计阶段所需的有关信号分析处理方面的必备知识。主要有信号预处理、信号的时域分析、信号的频谱分析、信号频谱的主要性质、信号处理中出现的现象、系统幅频特性、信号调制等知识。其中信号预处理包括信号的滤波处理、加窗处理、噪声添加处理；信号的的时域分析包括信号的强度描述、信号的相关分析(自相关分析、互相关分析)；信号的频谱分析包括信号的幅值谱分析、信号的功率谱分析、信号的倒频谱分析；信号频谱的主要性质包括信号的时移特性、信号的频移特性、信号的Parseval定理；信号处理中出现的现象包括信号的频率混叠现象、信号的频谱泄漏现象、信号的栅栏效应；信号调制包括信号的幅值调制、信号的频率调制、信号的相位调制。基于MATLAB GUI信号处理实验室的设计及实现MATLAB GUIMATLAB GUI简介 MATLAB是一种交互式、面向对象的程序语言，其结构完整，具有优良的移植性。它主要用于矩阵运算，同时在数据分析、自动控制、数字信号处理、绘图等方面也具有较大的功能。MATLAB采用一些常用的数学符号来表示问题及其解决方案，将计算、可视化和编程等功能集成于一个简单易用的开发环境中。MATLAB是一种基于不限维数组数据类型的内部交互系统，即能够进行矩阵和向量计算，也能够采用特定的方法在标量语言(例如C和Rortran)中编写程序。MATLAB为用户工作平台的管理和输入数据提供了便利的方法，同时还提供MATLAB易用程序M文件(该文件的扩展名为.m)的扩展和管理工具。MATLAB还采用一组被称为工具箱的特殊应用解答集，工具箱本身也是可理解的M文件集。在MATLAB中，图形用户界面是一种包含多种对象的图形窗口。用户必须对每一个对象进行界面布局和编程，从而使用户激活GUI每个对象时都能够执行相应的行为。另外，用户必须保存和发布所创建的GUI，使得GUI能够真正地得到应用Error! Reference source not found.。MATLAB为用户开发图形界面提供了一个方便高效的集成环境：MATLAB图形用户界面开发环境GUIDE(MATLABs Graphical User Interface Development Environment)。上述所有工作都能够使用GUIDE方便地实现。GUIDE主要是一个界面设计工具集。MATLAB将所有GUI支持的用户控件都集成在这个环境中并提供界面外观、属性和行为响应方式的设置方法。GUIDE将用户保持设计好的GUI界面保存在一个FIG资源文件中。同时还能够生成包含GUI初始化和组件界面布局控制代码的M文件。这个M文件为实现回调函数(当用户激活GUI某一组件时执行的函数)提供了一个参考框架。虽然使用用户自己编写的、所包含GUI所有发布命令的M文件也能够实现一个GUI，但是使用GUIDE执行效率更高，使用GUIDE不但能够交互式地进行组件界面布局，而且能够生成两个用来保存和发布GUI的文件。MATLAB具有以下几个特点：强大的数值运算功能：MATLAB有超