第六章高教版

PAGEPAGE297第6章随机信号分析实验随机信号分析是一门理论性较强的专业基础课程，在理论课上主要是通过数学方法讲述随机信号的分析及变换。因此，为了使初学者对随机信号的产生方法、统计特性、数字特征及其应用有较全面地了解；为了使初学者能够有效的将理论与实践相结合，特设置随机信号分析实验。随机信号分析实验是设计性、开放性实验。根据随机信号分析理论课程的内容，有目的性地选择了一些实验题目，以培养初学者的综合能力和创新能力，使初学者能够把所学的基本理论用于当前本专业研究的热点中，解决工程实践中的具体问题，同时提高初学者的学习兴趣，锻炼初学者的独立分析问题与解决问题的能力。本章实验一~实验四为随机信号分析理论课程的基础应用性实验；实验五和实验六为综合性实验，初学者可根据自己的兴趣选择实验内容。由于在进行本实验之前，初学者已有多门相关课程的基础理论知识和基本实验技能。因此，为了减少重复，本章只针对实验一的内容进行全面详细地描述，而其它实验，仅对实验目的、实验内容、实验任务与要求做了叙述。针对于其它五个实验，初学者可参考实验一的实验过程、实验方法及实验步骤完成实验内容并撰写实验报告。6.1实验一随机噪声特性分析6.1.1实验目的1、了解随机噪声的均值、均方值、方差、相关函数、概率密度、频谱及功率谱的特性。2、掌握白噪声、色噪声信号的产生方法。3、掌握随机信号的分析方法。4、熟悉常用的信号处理仿真软件平台：MATLAB、C/C++。6.1.2实验仪器与软件平台1、微计算机；2、双路输入示波器、信号源；3、MATLAB或C/C++仿真平台。6.1.3实验步骤1、根据本实验内容和要求查阅有关资料，设计并编写出相关的程序流程；2、自选MATLAB或C仿真软件；3、测试程序是否达到设计要求；4、分析实验结果是否与理论概念相符。6.1.4、实验内容1白噪声特性分析一、实验原理白噪声是指它的概率统计特性服从某种分布而它的功率谱密度又是均匀的。确切的说，白噪声只是一种理想化的模型，因为实际的噪声功率谱密度不可能具有无限宽的带宽，否则它的平均功率将是无限大，是物理上不可实现的。然而白噪声在数学处理上比较方便，所以它在通信系统的分析中有十分重要的作用。一般地说，只要噪声的功率谱密度的宽度远大于它所作用的系统的带宽，并且在系统的带内，它的功率谱密度基本上是常数，就可以作为白噪声处理了。白噪声的功率谱密度为：其中就是白噪声的均方值。白噪声的自相关函数为：白噪声的自相关函数是位于处、强度为的冲激函数。这表明白噪声在任何两个不同的瞬间的取值是不相关的。同时也意味着白噪声能随时间无限快的变化，因为它的带宽是无限宽的。通过几种常见的、概率分布为均匀分布、正态分布、指数分布、瑞利分布、方分布的随机噪声的统计特性分析，判断在这些随机噪声中哪些是白噪声。二、实验任务与要求(1)通过本实验要求掌握几种分布的随机噪声的共同点和不同点，重点在于系统测试与分析。实验系统框图如图6-1-1、图6-1-2所示图6-1-1图6-1-2(2)自选MATLAB或C/C++仿真软件之一产生几种概率分布的随机信号：随机序列的长度N=1024，随机序列包括均匀分布、正态分布、指数分布、瑞利分布、方分布。并计算这些随机数的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱密度并绘图。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱密度的物理意义。例：均值除了表示信号的平均值，它的物理含义是什么。相关函数当时为什么会有一个冲激，有什么物理意义。信号的时域波形有哪些特征，频域又有哪些特征，频谱及功率谱密度有什么差异。哪些噪声是白噪声，这些噪声是否符合白噪声的定义。在进行系统仿真和用波形图表示信号时如涉及到采样频率要根据输入信号来确定。(3)验证当数据长度增大时，白噪声的功率谱密度逼近。设产生N=20480长度的均值为2方差为3的正态随机数，从中取1024、10240、20480个点的功率谱密度作比较，观察这些随机数的功率谱密度随数据长度的变化。实际的白噪声功率谱密度不为常数，对白噪声功率谱密度求均值等于，初学者可以自己验证。(4)根据白噪声的特性，确定哪些随机信号属于白噪声范畴。根据分析确定白噪声与概率分布有关系吗？(5)试将五个均匀分布的随机噪声相叠加，其结果是否为高斯分布并分析其原因；试将五个指数分布的随机噪声相叠加，其结果是否为高斯分布并分析其原因。分析在这两种情况下的数据长度对叠加结果的影响。2白化滤波器的设计与分析一、实验原理在统计信号处理中，往往遇到等待处理的随机信号是非白色的，例如云雨、海浪、地物反射的杂乱回波等，它们的功率谱即使在信号通带内也非均匀分布。这样会给问题的解决带来困难。克服这一困难的措施之一是对色噪声进行白化处理。主要内容是设计一个稳定的线性滤波器或者一种白化变换方法，将输入的有色噪声变成输出的白噪声。下面探讨两种方法来实现白化问题。(1)白化滤波器将任意随机信号输入一个线性时不变滤波器，滤波器将白化为白噪声，这个滤波器就叫做白化滤波器。可以使用频域技术白化这个信号，用输入信号的功率谱密度，选择最小相位得到极点和零点都位于面左侧，这样就可以用以下关系构造白化滤波器：，选择最小相位滤波器保证逆滤波器稳定，必须保证在所有上都严格为正，这样就不会有奇点。白化噪声的功率谱为：，白化噪声的功率谱为常数，可见随机噪声已白化。(2)白化滤波器的设计方法首先计算色噪声自相关函数，根据色噪声的自相关函数，计算出色噪声的功率谱，然后根据公式，注意求倒数时不能为零。计算出白化滤波器的频谱。白化变换就是要构造一个白化矩阵，使色噪声与白化矩阵相乘后为白噪声。线性变换，使得的协方差矩阵为单位矩阵(即)。这里称作白化矩阵，它可以通过对色噪声矩阵的协方差矩阵的对角化求解来获得：，式中矩阵由的特征向量组成，为的特征值组成的对角矩阵。经过白化处理后，色噪声信号变换为具有单位方差的信号，且中各信号分量相互正交。白化变换方法总结：将生成的色噪声由一行变为的矩阵。计算色噪声的协方差矩阵。计算协方差矩阵的特征值以及特征向量。求白化变换矩阵。白噪声等于色噪声乘白化矩阵。再将生成的色噪声由的矩阵变为一行。二、实验任务与要求(1)实验系统框图如图6-1-3所示，其中、均为白噪声样本函数，其时域、频域波形图如图6-1-4所示。图6-1-3实验系统框图（a）白噪声信号时域波形图（b）白噪声信号频谱波形图图6-1-4高斯白噪声的时域、频域波形图要求测试的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱密度并用波形图表示。(2)产生色噪声：设计低通型衰减滤波器，实现色噪声。色噪声波形和功率谱密度如图6-1-5和6-1-6所示。图6-1-5色噪声信号图6-1-6色噪声功率谱要求测试的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱密度并用波形图表示。（3）白化噪声设计：使用白化变换方法或白化滤波方法产生白化信号。白化信号的时域、频域波形图6-1-7所示：（a）白化后噪声信号（b）白化噪声的功率谱密度图6-1-7白化信号的时域波形和功率谱密度波形计算白化信号的均值、均方值、方差、自相关函数、概率密度函数、频谱及功率谱密度。（4）注意观察色噪声信号、白化后的信号相关特性，并作结果分析。（5）思考在这里为什么要进行频谱和功率谱密度的计算与分析，理想白噪声能否进行傅里叶变换？3理想白噪声、带限白噪声特性比较一、实验原理若一个具有零均值的平稳随机过程，其功率谱密度在某一个有限频率范围内均匀分布，而在此范围外为零，则称这个过程为带限白噪声。带限白噪声可分为低通和带通型。二、实验任务与要求(1)通过实验掌握白噪声的特性以及带限白噪声的意义，重点在于系统测试与分析。选用MATLAB或C/C++仿真软件之一编写程序和仿真，系统框图如图6-1-8所示。图6-1-8带通滤波器系统框图(2)输入信号分别为高斯白噪声和均匀白噪声信号要求测试白噪声的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱，并绘图。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱的物理意义。(3)设计一个低通滤波器和一个带通滤波器。要求白噪声分别通过低通滤波器和带通滤波器后的信号能够表现出带限白噪声的特点。测试低通滤波器和一个带通滤波器的时域特性和频域特性以验证其正确性。(4)分别计算高斯白噪声、均匀白噪声经过低通滤波、带通滤波器后的均值、均方值、方差、概率密度、自相关函数、频谱及功率谱密度，并加以分析。(5)所有结果均用波形图来表示。(6)白噪声在什么情况下为带限白噪声？4色噪声的产生与特性分析一、实验原理把除了白噪声之外的所有噪声都称为有色噪声。就像白光一样，除了白光就是有色光。色噪声中的几个典型实例：(1)粉红噪声。粉红噪音是自然界最常见的噪音，简单说来，粉红噪音的频率分量功率主要分布在中低频段。从波形角度看，粉红噪音是分形的，在一定的范围内音频数据具有相同或类似的能量。从功率（能量）的角度来看，粉红噪音的能量从低频向高频不断衰减。粉红噪声的能量分布在任意的同比例带宽中是相等。在给定频率范围内(不包含直流成分)，随着频率的增加，其功率密度每倍频程下降3dB(密度与频率成反比)。每倍频的功率相同，但要产生每倍频程3dB的衰减非常困难。因此，没有纹波的粉红噪声在现实中很难找到。粉红噪声的能量在低频段频率接近0Hz(不包括0Hz)而在高频段频率接近二十几千赫，而且它在等比例带宽内的能量是相等的(误差只不过0.1dB左右)。粉红噪声的功率谱图如图6-1-9所示。图6-1-9粉红噪声的功率谱密度(2)红噪声(海洋学概念)。这是有关海洋环境的一种噪声，由于它是有选择地吸收较高的频率，因此称之为红噪声。(3)橙色噪声。该类噪声是准静态噪声，在整个连续频谱范围内，功率谱有限，零功率窄带信号数量也有限。这些零功率的窄带信号集中于任意相关音符系统的音符频率中心上。由于消除了所有的合音，这些剩余频谱就称为“橙色”音符。(4)蓝噪声。在有限频率范围内，功率密度随频率的增加每倍频增长3dB(密度正比于频率)。对于高频信号来说，它属于良性噪声。(5)紫噪声。在有限频率范围内，功率密度随频率的增加每倍频增长6dB(密度正比于频率的平方值)。(6)灰色噪声。该噪声在给定频率范围内，类似于心理声学上的等响度曲线，因此在所有频率点的噪声电平相同。(7)棕色噪声。在不包含直流成分的有限频率范围内，功率密度随频率的增加每倍频下降6dB(密度与频率的平方成反比)。该噪声实际上是布朗运动产生的噪声，它也称为随机飘移噪声或醉鬼噪声。(8)黑噪声(静止噪声)。有源噪声控制系统在消除了一个现有噪声后的输出信号。在20kHz以上的有限频率范围内，功率密度为常数的噪声，一定程度上它类似于超声波白噪声。这种黑噪声就像“黑光”一样，由于频率太高而使人们无法感知，但它对你及其周围的环境仍有影响。二、实验任务与要求(1)通过实验掌握白噪声、色噪声特点，以及如何产生色噪声，重点在于系统测试与分析。选用MATLAB或C/C++仿真软件之一编写程序和仿真。(2)产生高斯白噪声。要求测试白噪声的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱并用波形图表示。分析实验结果，搞清楚均值、均方值、方差，自相关函数、频谱及功率谱的物理意义。(3)产生高斯色噪声：要求自己设计一个系统，使高斯白噪声通过该系统后变为高斯色噪声。有一种方法可供参考：让高斯白噪声通过低通型、带通型、高通型衰减滤波器中的任意一个就可以产生高斯色噪声。高斯色噪声及功率谱密度如图6-1-10所示。（a）产生的高斯色噪声时域波形图（b）高斯色噪声的功率谱密度图6-1-10高斯色噪声时域及频域波形图仅做参考的低通衰减滤波器参数如下：通带截至频率1Hz阻带截至频率10KHz通带衰减1dB阻带衰减20dB(3)产生粉红色噪声：要求自己设计一个系统，使高斯白噪声通过该系统后变为粉红色噪声。粉红色噪声的频谱、功率谱必须用dB来表示，以方便观察其衰减程度。(4)计算高斯色噪声、粉红色噪声的均值、均方值、方差、相关函数、概率密度、频谱及功率谱密度。所有结果均用波形图来表示，并能读出具体数值。5随机信号平稳特性分析一、实验原理平稳随机过程是在时间平移下概率性质不变的随机过程。其统计特性是，任意有限维分布函数不随时间的推移而改变；当过程随时间的变化而产生随机波动时，其前后状态相互联系，即不但它的当时情况，而且它的过去情况对未来都有不可忽视的影响。按照描述平稳随机过程的统计特性的不同，平稳随机过程分为严平稳随机过程和宽平稳随机过程。二、实验任务与要求(1)通过实验掌握随机信号的方法以及白噪声信号的特点。自选MATLAB或C/C++仿真软件之一编程并仿真。(2)生成满足几种概率分布的随机信号，具体要求：产生的10000个泊松分布随机数，计算它们的均值、均方值、方差、概率密度、频谱、功率谱密度，自相关函数，用波形图来表示。(3)产生高斯分布随机数与。要求测试这两种噪声的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱，用波形图来表示。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱的物理意义。统计分析：二维正态分布，的联合概率密度函数为其中，求二维正态分布的边沿分布、，并用波形图来表示。(4)对正态分布随机数取几个不同的样本值，计算它们的数字特征，分析是否满足平稳性和遍历性。6.1.5实验报告简述主要实验内容、基本原理和结论，重点叙述实验过程中遇到的问题及解决方法，讨论该实验内容和方法在相关领域的应用和发展，试举例说明。6.2实验二平稳随机过程的抽取、插值方法的探讨6.2.1实验目的1、了解确定信号的采样与平稳随机信号的采样之间的关系，掌握信号的采样定理及其应用；2、掌握随机信号的均值、方差、自相关函数、概率密度、频谱及功率谱密度的特性；3、掌握随机信号的分析方法；4、熟悉常用的信号处理仿真软件平台：MATLAB或C/C++。6.2.2实验内容一、实验原理确定信号的采样符合香农定理，那么随机信号的采样呢？答案是肯定的。若X(t)为平稳随机过程，且具有零均值，它的功率谱密度限于之间。当满足条件时，便可将按它的振幅样本展开为：上式为平稳随机过程的采样定理，式中T为采样周期。二、实验任务与要求(1)通过实验掌握随机信号的抽取与插值的基本原理方法及方法在实际应用中的意义。计算信号的数字特征，理解它们的物理概念。用MATLAB或C/C++仿真软件之一编写程序和仿真。系统框图如图6-2-1图6-2-1抽取、插值系统框图(2)输入信号：是在一个波形周期有80个取样点的正弦信号+。信号的真实频率取决于采样频率。假设采样频率为8000Hz，则输入信号频率就是100Hz。为高斯白噪声。输入信号的时域、频域波形图如图6-2-2所示，高斯白噪声的时域、频域波形图如图6-2-3所示。（a）输入信号时域波形图（b）输入信号频域波形图图6-2-2输入信号的时域、频域波形图（a）白噪声信号时域波形图（b）白噪声信号频域波形图图6-2-3高斯白噪声的时域、频域波形图要求测试白噪声的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱并用波形图表示。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱的物理意义。(3)滤波器设计：要求信号经滤波器后保留有用信号。滤波器的结构及参数和所涉及到的采样频率取值根据滤波器在这里所起的作用、输入信号本身的特点所确定，因此这里不作规定。滤波器设计好之后，要求测试它的频率特性并将频率特性曲线用波形图来表示，以验证是否符合要求。(4)信号抽取：所谓信号抽取也就是信号抽样率的降低。将低通滤波器输出信号作为原始信号，用抽取的方法分别为以、对原始信号进行抽取。(5)信号插值：将、抽取信号分别作为原始信号，用插值的方法，分别还原原始信号（滤波后的信号）。(6)计算、、、的均值、均方值、方差、自相关函数、频谱及功率谱密度。(7)对采样前后、插值前插值后信号进行比较。观察信号频谱的变化及周期延拓性。(8)讨论的自相关函数、功率谱密度与的自相关函数、功率谱密度之间的关系。6.3实验三随机信号通过线性和非线性系统后的特性分析6.3.1实验目的1、了解随机信号的均值、均方值、方差、自相关函数、互相关函数、概率密度、频谱及功率谱特性。2、研究随机信号通过线性系统和非线性系统后的均值、均方值、方差、自相关函数、互相关函数、概率密度、频谱及功率谱有何变化，分析随机信号通过线性系统和非线性系统后的特性。6.3.2实验内容1随机信号通过线性系统和非线性系统后的特性分析一、实验原理(1)随机信号的分析方法在信号系统中，可以把信号分成两大类：确定信号和随机信号。确定信号具有一定的变化规律，而随机信号无一定的变化规律，需要用统计特性进行分析。在这里引入了随机过程的概念。所谓随机过程，就是随机变量的集合，每个随机变量都是随机过程的一个取样序列。随机过程可分为平稳的和非平稳的、遍历的和非遍历的。如果随机信号的统计特性不随时间的推移而变化，则随机过程是平稳的。如果一个平稳的随机过程它的任意一个样本都具有相同的统计特性，则随机过程是遍历的。下面讨论的随机过程都认为是平稳的遍历的随机过程，因此，可以取随机过程的一个样本来描述随机过程的统计特性。随机过程的统计特性一般采用主要的几个平均统计特性函数描述，包括均值、均方值、方差、自相关函数、互相关函数、概率密度、频谱及功率谱密度等。①随机过程的均值均值表示集合平均值或数学期望值。基于随机过程的各态历经性，可用时间间隔T内的幅值平均值表示，即：均值表达了信号变化的中心趋势，或称之为直流分量。

②随机过程的均方值信号的均方值，或称为平均功率，其表达式为：均方值表示信号的强度，其正平方根，又称有效值，也是信号平均能量的一种表达。③随机信号的方差信号的方差定义为：称为均方差或标准差。可以证明，其中：描述了信号的波动量；描述了信号的静态量，方差反映了信号绕均值的波动程度。④随机过程的互相关函数信号的相关性是指客观事物变化量之间的相依关系。对于平稳随机过程和在两个不同时刻和时刻幅值的关联程度，可以用互相关函数表示。在离散情况下，信号和的互相关函数定义为：。⑤随机过程的频谱：信号频谱分析是采用傅立叶变换将时域信号变换为频域信号，从另一个角度来了解信号的特征。时域信号的付氏变换为：⑥随机过程的功率谱密度随机信号的功率谱密度是随机信号的各个样本在单位频带内的频谱分量消耗在一欧姆电阻上的平均功率的统计均值，是从频域描述随机信号的平均统计参量，表示的平均功率在频域上的分布。它只反映随机信号的振幅信息，而未反映相位信息。随机过程的功率谱密度为：(2)线性系统特性系统的数学模型满足叠加原理。对于任意常数和，输入信号和，有则称系统为线性系统，线性系统有下面的一些重要性质：叠加性、比例性、微分性、积分性、频率保持性等。(3)非线性系统特性在一般电子设备中，除了线性电路之外，通常还包括一些非线性电路，例如检波器、限幅器、鉴频器等。非线性电路具有下述特点：①叠加原理已不适用，当信号与噪声共同通过非线性电路时，不能像线性电路那样将它们分开研究。②会发生频谱变换，其输出产生了输入信号中没有的新频率分量，例如输入信号的各次谐波。二、实验任务与要求(1)通过实验要求掌握线性系统、非线性系统基本特性，比较未通过系统的随机信号与通过系统后的随机信号的特性。自选MATLAB或C/C++仿真软件之一。实验系统框图如图6-3图6-3-1线性系统、非线性系统测试(2)输入信号、噪声的测试与分析输入信号，其中：、、为1KHz、2KHz、3KHz，如图6-3-2所示。（a）输入信号时域波形图（b）输入信号频域波形图图6-3-2输入信号的时域、频域波形图噪声为高斯白噪声，其时域、频域波形图如图6-3-3所示。（a）高斯白噪声信号时域波形（b）高斯白噪声信号频域波形图图6-3-3高斯白噪声的时域、频域波形要求测试的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱并用波形图表示。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱的物理意义。(3)线性系统设计及测试研究随机信号经过线性系统后的线性变换问题，需要设计一个线性系统。线性系统设计成一个滤波器。要求信号经滤波器后只剩下2KHz频率的信号。滤波器的结构及参数和所涉及到的采样频率取值根据滤波器在这里所起的作用、输入信号本身的特点来确定，因此这里不作规定。滤波器设计好之后，要求测试它的幅频特性，以验证是否符合要求。经滤波器后随机信号的时域、频域特性如图6-3-4所示。（a）输入信号通过带通滤波器波形图（b）输入信号通过带通滤波器频域波形图图6-3-4信号通过滤波器后的时域、频域(4)非线性系统设计及测试研究随机信号经非线性系统后的非线性变换问题，需要设计一个非线性系统。在这里非线性系统分别设计成一个限幅器、一个平方律器件。平方律器件硬件实现一般利用二极管的特性曲线如图6-3-5所示图6-3-5二极管的特性曲线时域特性：当时,当时,。频域特性：。软件实现也是利用这个平方特性。限幅器是一种波形变换或整形电路。当输入信号在一定范围内变化时，输出电压跟随输入电压相应变化，完成信号的传输。而当输入电压超过这一范围时，其超过的部分就被削去，输出电压保持不变。由于限幅器能将一定范围以外的输入波形削去，所以限幅器又称削波器。限幅器、平方律器件设计好之后，要求测试它的频率特性并画出频率特性曲线，观察是否符合要求。经平方律器件后随机信号的时域、频域特性如图6-3-6所示。（a）经平方律检波器后时域波形图（b）经平方律检波器后频域波形图图6-3-6经平方律检波器后的时频、频域(5)系统测试与分析按以上要求编写好程序，计算、、、、的均值、均方值、方差，自相关函数、频谱及功率谱，并对两路系统输入和输出信号进行比较。2随机信号经线性、非线性系统并行叠加后分析一、实验原理随机信号经过并行系统的原理框图如图6-3-7图6-3-7二、实验任务与要求(1)通过实验掌握线性系统的特点以及分析方法，重点在于系统测试与分析。自选MATLAB或C/C++仿真软件之一。实验系统框图如图6-3-7(2)输入信号，其中：、为1KHz、2KHz，如图6-3-8所示。（a）时域波形图（b）频域波形图图6-3-8输入信号的时域、频域波形图高斯白噪声的时域、频域波形图如图6-3-9所示（a）白噪声信号时域波形图（b）白噪声信号频域波形图图6-3-9高斯白噪声的时域、频域要求测试的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱密度并用波形图表示。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱的物理意义。(3)、是线性系统，是非线性系统。系统：设计成滤波器，要求输入信号x(t)经后输出为1Khz正弦信号。系统：设计成滤波器，要求输入信号x(t)经后输出为2Khz正弦信号。滤波器的结构及参数和所涉及到的采样频率取值根据滤波器在这里所起的作用、输入信号本身的特点来确定，这里不作规定。滤波器设计好之后，要求测试它的频率特性并画出频率特性曲线，以验证是否符合要求。系统：设计成平方律器件。的时域、频域波形图如图6-3-10所示。（a）时域波形图（b）的功率谱密度图图6-3-10y(t)的时域、频域(4)系统测试与分析计算测试图6-3-7中、、、、的均值、均方值、方差、频谱、功率谱密度，自相关函数，并绘出函数曲线。分析信号经三个系统并行叠加后频谱发生的变化；分析线性、非线性系统分别在这里的作用。3随机信号经线性、非线性系统串行后性能分析一、实验原理随机信号经过串行系统的原理框图如图6-3-11所示图6-3-11串行系统的原理二、实验任务与要求(1)通过实验掌握线性系统、非线性系统的特点以及分析方法。自选MATLAB或C/C++仿真软件之一。(2)输入信号：。其正弦波的频率为2000Hz，噪声为高斯分布的白噪声。图6-3-12为输入信号的时域、频域波形图；图6-3-13为的时域、频域波形图。（a）时域波形图（b）频域波形图图6-3-12的时域、频域波形图（a）白噪声信号时域波形图（b）白噪声信号频域波形图图6-3-13高斯白噪声的时域、频域波形要求测试的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱密度并绘图。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱的物理意义。(3)、是线性系统，是非线性系统。系统：设计成滤波器，要求输入信号x(t)经后输出为2Khz正弦信号。系统：设计成滤波器，输出信号经后输出为4Khz正弦信号。滤波器的结构及参数和所涉及到的采样频率取值根据滤波器在这里所起的作用、输入信号本身的特点来确定，因此这里不作规定。滤波器设计好之后，要求测试它的频率特性并画出频率特性曲线，以验证是否符合要求。系统：设计成平方律器件。(4)信号通过、、后，使信号频率得到改变，如图6-3-14所示，仔细观察其结果。（a）信号时域波形图（b）信号频域波形图图6-3-14y(t)的时域、频域波形图(5)计算图6-3-11中、、、的均值、均方值、方差、频谱、功率谱，自相关函数，并绘出函数曲线。分析这个系统的功能，分析线性、非线性系统分别在这里的作用。4随机信号经线性、非线性系统串并后性能分析一、实验原理随机信号进行串并行系统的原理框图如图6-3-15所示。图6-3-15串并系统的原理框图二、实验任务与要求(1)通过实验掌握线性、非线性系统的特点以及分析方法，重点在于系统测试与分析。自选MATLAB或C/C++仿真软件之一编写程序和仿真。(2)输入信号：x(t)=方波+n(t删除x(t)=方波+n(t)，以前重复)。其方波的基频为1000Hz，为高斯分布的白噪声。图6-3-16为的时域、频域波形图，图6-3-17为的时域、频域波形图。删除x(t)=方波+n(t)，以前重复（a）输入信号时域波形图（b）输入信号频域波形图图6-3-16x(t)的时域、频域波形（a）白噪声信号时域波形图（b）白噪声信号频域波形图图6-3-17高斯白噪声的时域、频域波形图要求测试的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱密度并绘图。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱的物理意义。(3)、、是线性系统，是非线性系统。系统：设计成滤波器，要求输入信号x(t)经后输出为1Khz正弦信号。系统：设计成滤波器，要求输入信号x(t)经后输出为1KHz、3KHz混合信号。滤波器的结构及参数和所涉及到的采样频率取值根据滤波器在这里所起的作用、输入信号本身的特点来确定，因此这里不作规定。滤波器设计好之后，要求测试它的频率特性并画出频率特性曲线，以验证是否符合要求。系统：设计成平方律器件。要求输出信号的频率为3个频率成分的叠加，即1KHz、2KHz、3KHz。图6-3-18为的时域、频域波形图。（a）信号时域波形图（b）信号频域波形图图6-3-18y(t)的时域、频域波形(5)计算图6-3-15中、、、、的均值、均方值、方差、频谱、自相关函数、功率谱密度，并用波形图表示。分析线性、非线性系统分别在这里的作用。6.4实验四弱信号的提取方法研究6.4.1实验目的了解随机信号分析理论如何在实践中应用，掌握随机弱信号的检测及分析的几种方法。6.4.2实验内容1弱信号的提取方法研究一、实验原理(1)弱随机信号的检测及提取方法因为噪声总会影响信号检测的结果，所以信号检测是信号处理的重要内容之一，低信噪比下的信号检测是目前检测领域的热点，而强噪声背景下弱信号的提取又是信号检测的难点，其目的就是消除噪声，将弱信号从背景噪声中提取出来，或者用一些新技术和新方法来提高检测系统输出信号的信噪比。噪声主要来自于检测系统本身的电子电路，系统外的空间或其它设备产生的高频电磁信号等，通常从两种不同的途径来解决：①降低系统本身的噪声，提高信噪比。②采用相关接收技术，可以保证在被测信号功率小于噪声功率的情况下，仍能检测出信号。在电子学系统中，采用低噪声放大技术，选取适当的滤波器限制系统带宽，以抑制内部噪声和外部干扰，大大改善系统的信噪比，当信号较弱时，也能得到的结果。但当信号非常弱，比噪声小几个数量级甚至完全被噪声深深淹没时，上述方法就不会有效。当已知噪声中的有用信号波形时，利用信号和噪声在时间特性上的差别，可以用匹配滤波的方法进行检测。但当弱信号是未知信号时，则无法利用匹配滤波的方法进行检测。经过分析，白噪声为一个具有零均值的平稳随机过程，所以选取任一时间点，在该点前一段时间内将信号按时间分成若小段后，然后在选取时间点处将前面所分的每小段信号累加，若为白噪声信号，则时间均值依然为零，但当噪声中存在有用信号时，则时间均值不为零，由此特性，就可对强噪声背景中是否存在弱信号进行判定。对弱信号检测与提取有很多方法，常采用以下方法进行检测，这些检测方法都可以在与信号处理的参考文献中查找到。①自相关检测方法传统的自相关检测技术是应用信号周期性和噪声随机性的特点，通过自相关运算达到去除噪声的检测方法。由于信号和噪声是相互独立的过程，根据自相关函数的定义，信号只与信号本身相关与噪声不相关，而噪声之间一般也是不相关的。假设信号为，噪声为，则输入信号其相关函数：对于具有各态历经性的过程，可以利用样本函数的时间相关函数来替代随机过程的自相关函数。②多重自相关法多重自相关法是在传统自相关检测法的基础上，对信号的自相关函数再多次作自相关。即令：式中，是和的叠加；是和的叠加。信号经过相关运算后增加了信噪比，但其改变程度是有限的，因而限制了检测弱信号的能力。多重相关法将当作，重复自相关函数检测方法步骤，自相关的次数越多，信噪比提高的越多，因此可检测出淹没于强噪声中的弱信号，如图6-4-1所示。图6-4-1多重自相关③双谱估计理论及算法双谱变换是对信号的三阶累积量进行二维傅立叶变换，假定为零均值，三阶实平稳随机序列，其三阶相关函数为：则其双谱就定义为：对于经典的双谱估计方法，可分直接法和间接法两种。④时域方法主要是叠加平均技术，它对时域信号进行多次叠加取平均值以得到信号，其算法误差较大。采用对混合信号反复取样，累加平均的方法，使噪声信号自相削弱，从而再现有用信号。根据采样定理，如果信号在处连续，则在处的取样可表示为。如果在的各处均连接，每隔时间对取样一次，这个取样用表示，则式中为单位强度的周期性冲激函数序列，里删除的傅立叶变换为：删除若的频谱函数用表示，则由频域的卷积定理得：，应用冲激函数的取样性质，考虑到得；，式中是的周期.只要，即（为信号最高频率），亦即。因此只要满足条件，即可获得信号的全部信息，当然这里面也含有噪声。由此可见，系统工作频带越窄，叠加次数越多，等效噪声带宽越小，则系统的输出信噪比越高。但经过足够次数的采样、累加平均后，信噪比会大大提高。当增加时，系统工作频带变窄，这样就抑制了噪声，提高了信噪比。(2)小波变换算法针对于实际应用中完全被噪声淹没情况下的弱信号提取的问题，依据白噪声信号的小波变换系数相对于有用信号的小波变换系数小的特点，利用小波变换对信号进行消噪来提取弱信号，小波变换能够有效的消除噪声，将有用弱信号从受噪声污染的信号中提取出来。目前，众多新的方法正在研究中，有兴趣的初学者可以通过查找相关的参考文献加以关注。二、实验任务与要求(1)通过实验掌握一种弱信号的检测提取方法，重点在于系统测试与分析。(2)输入信号：白噪声加上弱周期信号组成（信噪比）。由信号发生器产生的正弦波信号，周期信号的频率自己确定，从计算机的音频输入端输入。要注意的是输入信号与计算机的声卡采样频率之间的关系。或者用MATLAB或C/C++仿真软件之一来实现正弦波信号和高斯白噪声。输入数据和高斯白噪声的时域、频域波形图如图6-4-2和图6-4-3所示。（a）读取的信号时域波形图（b）混合信号的频域波形图图6-4-2输入数据的时域、频域波形图（a）白噪声信号时域波形图（b）白噪声信号频域波形图图6-4-3高斯白噪声的时域、频域波形要求测试白噪声的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱密度并绘图。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱密度的物理意义。(3)弱信号提取过程框图如图6-4-4所示。图6-4-4弱信号提取过程框图(4)滤波器设计：设计一个滤波器。要求保留有用信号，滤波器的结构及参数和所涉及到的采样频率取值根据滤波器在这里所起的作用、输入信号本身的特点所以确定，因此这里不作规定。滤波器设计好之后，要求测试它的频率特性并画出频率特性曲线，以验证是否符合要求。(5)信号提取：选择一种算法，将弱信号从噪声中提取出来，提取信号的时域、频域波形图如图6-4-5所示。（a）信号的二重自相关法恢复信号（b）提取信号的频域波形图图6-4-5提取信号的时域、频域波形图(6)首先计算点的均值、均方值、方差、频谱及功率谱密度，确定输入信号中包含着有用信号；然后计算、的均值、均方值、方差、频谱及功率谱密度，并画出曲线；确定信号的周期。(7)将输入信号改为方波、三角波时选择那些方法更接近还原原始信号，结果如何？2维纳滤波器去除噪声的性能分析一、实验原理通常，真实信号在其输出的过程中会受到外界噪声的干扰，因此所得到的输出信号就会发生信号失真。用滤波器就可以处理噪声，从而得到一个较接近原来真实信号的输出信号。如果信号和噪声的频带不重叠，用经典滤波器就可以很好的实现滤波；但是如果信号和噪声的频带相互重叠，则要采用现代滤波器。其中以维纳滤波的功能更为突出，它是利用平稳随机过程的相关特性和频谱特性对混有噪声的信号进行滤波的方法。假设维纳滤波器的输入信号是，噪声是，信号加噪声是，维纳滤波器的冲激响应是，维纳滤波器原理如图6-4-6所示。图6-4-6维纳滤波器信号加噪声通过线性滤波器得到的目标信号估值：，估计误差:。估计误差为可正可负的随机变量,用它的均方值描述误差的大小显然更为合理。均方误差为：维纳滤波器的设计就是要确定，而可以从求解维纳-霍普方程得到。维纳-霍普方程为：式中，表示噪声加信号的自相关函数，表示信号和的互相关函数。求解维纳-霍普方程得到维纳滤波最优冲激响应：式中是求逆矩阵，是求转置矩阵。二、实验任务与要求(1)通过实验掌握维纳滤波器的设计方法并与传统滤波器进行比较。自选MATLAB或C/C++仿真软件之一编程并仿真。(2)信号经过维纳滤波器与经典滤波器的系统方框图如图6-4-7所示。图6-4-7维纳滤波器与经典滤波器的系统方框图产生输入信号：，其中：、、为1KHz、2KHz、3KHz，幅值为1v，为高斯白噪声。的时域、频域波形图如图6-4-8所示，高斯白噪声的时域、频域波形图如图6-4-9所示。（a）输入信号时域波形图（b）输入信号频域波形图图6-4-8输入信号的时域、频域波形图（a）白噪声信号时域波形图（b）白噪声信号频域波形图图6-4-9高斯白噪声的时域、频域波形要求测试白噪声的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱密度并绘图。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱密度的物理意义。(3)滤波器。要求信号经滤波器后保留有用信号。滤波器的结构及参数和所涉及到的采样频率取值根据滤波器在这里所起的作用、输入信号本身的特点所确定，因此这里不作规定。滤波器设计好之后，要求测试它的频率特性并画出频率特性曲线，以验证是否符合要求。(4)维纳滤波程序设计：要求：设计维纳滤波时，输入信号与滤波器输入相同。①改变噪声强度（信噪比分别0dB,-5dB,-10dB)，观察维纳滤波器输出数据。②改变维纳滤波器长度N（N=10、50、100、500），观察维纳滤波器输出波形的变化，观察维纳滤波器滤除噪声的效果。(5)计算、、的均值、均方值、方差、自相关函数、频谱及功率谱密度。(6)比较滤波器、维纳滤波器对输入信号的影响。3利用正交系统检测窄带信号方法一、实验原理当正态窄带噪声叠加正弦信号时，其合成过程为：令：则：式中：为合成过程的包络，为合成过程的相位，、均含有信息分量。窄带信号检测系统的框图如图6-4-10所示。图6-4-10窄带信号检测系统二、实验任务与要求(1)通过实验掌握如何将弱的周期信号检测出来的方法，重点在于系统测试与分析。自选MATLAB或C/C++仿真软件之一编程并仿真。(2)输入信号：为：其中，包络频率为1KHz,幅值为1v，在[]上均匀分布，为正态分布高斯噪声。图6-4-11为的时域、频域波形图。（a）的时域波形图（b）的频域波形图图6-4-11的时域、频域波形图要求测试白噪声的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱，并绘图。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱的物理意义。(3)滤波器设计：要求信号经滤波器后保留有用信号。滤波器的结构及参数和所涉及到的采样频率取值根据滤波器在这里所起的作用、输入信号本身的特点所确定，因此这里不作规定。滤波器设计好之后，要求测试它的频率特性并画出频率特性曲线，以验证是否符合要求。(4)输入信号通过希尔伯特转换器。希尔伯特转换(Hilberttransform)：是将信号与做卷积。因此，可以将的希尔伯特转换看成是将通过一个冲激响应为的线性滤波器，希尔伯特转换相当于一个正交滤波器。(5)用波形图表示和，其表达式为，。包络信号的时域、频域波形图如图6-4-12所示。（a）时域波形图（b）频域波形图图6-4-12时域、频域波形图(6)计算图6-4-10中的、、、、、、、、各点的信号均值、均方值、方差、自相关函数、频谱及功率谱密度，并用波形图表示。4利用系统的积累特性检测窄带信号方法一、实验原理用包络检测法来检测淹没在噪声中的周期信号，可以采用信号积累的方法来检测周期信号(信号周期已知)。信号积累的检测系统框图如图6-4-13所示。图6-4-13信号积累的检测系统二、实验任务与要求(1)通过实验掌握如何将弱的周期信号检测出来的方法。自选MATLAB或C/C++仿真软件之一编程并仿真。(2)输入信号：为：其中，包络频率为1KHz,幅值为1v，在上均匀分布，为正态高斯噪声。图6-4-14为的时域、频域波形图。（a）时域波形图（b）频域波形图图6-4-14的时域、频域波形图测试白噪声的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱，并绘图。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱的物理意义。(3)窄带滤波器：设计成一个低通滤波器。要求信号经低通滤波器后保留有用信号。低通滤波器的结构及参数和所涉及到的采样频率取值根据滤波器在这里所起的作用、输入信号本身的特点自己确定，因此这里不作规定。低通滤波器设计好之后，要求测试它的频率特性并画出频率特性曲线，以验证是否符合要求。(4)平方律检波器系统设计成一个平方律检波器。平方律检波器的特性曲线如图6-4-15所示。图6-4-15平方律检波器的特性曲线设计小信号平方律检波器，频率范围<100KHz，输入信号幅度小于1v。(5)归一化参数：，为输入信号的方差。(6)独立取样：取10个长度为N的序列，序列N长度必须是信号周期的整数倍数。这里需要求出一个信号周期有多少点，这样才能知道取得的长度是否为整周期。如果不为整周期则误差就会加大，在这里通过仿真实验进行验证，如N=1024的情况。图6-4-16为的时域、频域波形图，从图6-4-16中可以看出，包络信号已恢复。（a）的时域波形图（b）的频域波形图图6-4-16y(t)的时域、频域波形图(7)计算、低通滤波器输出、平方律检波输出、输出信号的均值、方差、频谱、功率谱密度，自相关函数，并绘出函数曲线。6.5实验五系统辨识的基本方法在实际工程应用中，常出现已知系统的输入和输出，期望得到该系统特性的场合，而系统辨识所研究的问题恰好是根据系统的输入信号和输出信号来确定系统的传输特性。6.5.1实验目的通过查阅相关的文献资料，了解分析电子系统的测量和识别方法，掌握随机信号的特性，包括均值、均方值、方差、自相关函数、互相关函数、概率密度、频谱及功率谱密度等在分析电子系统特性方面的应用，利用系统的输入和输出对系统进行分析，以得到系统的特性。6.5.2实验内容1实验原理系统辨识是在对某系统不知其结构的情况下，如果要确定该系统的冲激响应，可以通过将白噪声作为系统的输入，测量系统的输出为，于是可求得和的互相关函数，由于白噪声的自相关函数，而所以，。将通过一个低通滤波器，获得线性系统单位冲激响应，其原理框图如图6-5-1所示。图6-5-1利用互相关测量线性系统单位冲激响应掌握白噪声的特性，以及探讨这种测试方法的意义，重点在于系统测试与分析。输入信号为高斯白噪声，其时域、频域波形图如图6-5-2所示。（a）白噪声信号时域波形图（b）白噪声信号频域波形图图6-5-2高斯白噪声的时域、频域波形图要求测试白噪声的均值、均方值、方差，自相关函数、概率密度、频谱及功率谱密度并绘图。分析实验结果，掌握均值、均方值、方差，自相关函数、频谱及功率谱的物理意义。被测系统：①被测系统是一个低通滤波器。低通滤波器的通带为0KHz-1KHz、通带衰减小于1dB、阻带衰减大于35dB。②被测系统是一个带通滤波器。带通滤波器的通带为1KHz-2KHz、通带衰减小于1dB、阻带衰减大于35dB。③被测系统是一个微分器。如果被测系统是低通滤波器，输出是它的冲激响应为函数。实验结果应该为图6-5-3所示：同理，如果被测系统是带通滤波器，输出是它的冲激响应，如图6-5-4