matlab中平稳性检验,基于Matlab的信号平稳性检验系统
信号的平稳性检验在随机信号处理中起着十分基础的作用。由于平稳信号和非平稳信号的性质差别显著,因此在处理信号之前先行判断它的平稳性就显得尤为重要。虽然信号平稳性的定义十分明确,但是实际判断过程却是复杂的,例如观察尺度对信号平稳性判断就有很大的影响。
这一领域的研究已经取得了一定的成果。一些人提出了受限和带参数的非平稳性判定方法,而另一些人则将他们的平稳性判定建立在对原始数据的一些假设上。而对于更一般信号的平稳性检验的研究还没有取得太多成果。文献 [2,3]中又提到了这一问题,并且提出了一种新的检验平稳性框架。这一框架混合了时频透视法和有名的替代数据法。它的基本思想是引入“可控噪声”,即替代数据。并且由于替代数据的一些特性,它可以作为平稳性的评判标准。本文参考了文献[6]中的平稳性检验方法,设计了一个信号平稳性检验系统,并在 Matlab的GUI开发环境下实现了图形用户界面的设计。实践表明,本系统不但提供了友好的用户界面,并且可以方便地完成信号的平稳性检验。
1 平稳性检验原理
1.1 平稳性定义及其检验的重要性
假设有一个高斯过程{xl(t))(- ∞
为任意确定t时刻的全体平均。同时:
被称为自相关函数(ACF)。
对于一个弱平稳过程,它的μx(t) 和r(t1,t2)都是时不变的或者说与时间无关的。因此有:
式中:τ=t1-t2被称为时延。因此,对于平稳高斯过程{xl(t)},它的自相关函数或者它的功率谱密度函数(PSD)为:
足以确定它的性质。
另一方面,如果{xl(t))是非平稳的,它的μx(t)和r(t1,t2)就是时变的或者说和时间相关的。这样它的PSD就应该放在时频域分析。
由此可见,平稳性检验是任何信号处理前必不可少的一步,它决定了后续处理可以使用何种方法。
1.2 替代数据
替代数据的概念最初是由Theiler和其合作作者提出的,这种技术是用来产生一种所谓的“替代数据”,这种替代数据是平稳的,同时保持了原数据的一些相关的统计特性。
Theiler在文献[4]中提出了一种具体的产生替代数据的方法。由这种方法产生的替代数据是平稳的,同时保持了原数据的二阶统计特性。具体地说,替代数据保持了原数据功率谱的幅度值不变。
根据Wiener- Khintchin理论,信号的功率谱等于其傅里叶变换的幅值平方。因此保持信号的功率谱幅度值不变,就是保持其傅里叶变换的幅度值不变。因此,假设原数据为x(t),它的傅里叶变换为
产生。其中,φf是在[-π,π]上均匀分布的随机相位。这样就保证了s(t)和x(t)有相同的傅里叶变换幅值。在下面的例子中也可以看到,这样产生的s(t)也是平稳的。
1.3时频分布
时频分布主要用于分析非平稳随机信号的功率谱。由于非平稳随机信号的功率谱是时变的,因此在原来功率谱的基础上再引入时间轴,成为时频分布(TFD)。TFD可以显示出信号的功率谱随时间的变化情况。
具体来说,根据文献[9]中的定义,信号x(t)的时频分布Sx,K(t,f)可以表示为:
1.4 平稳性检验
平稳性可以体现在频谱随时间的波动上。具体来说,对于平稳信号,其频谱不随时间变化;而对于非平稳信号,其频谱会随时间改变。因此,可以通过比较不同时间点上频谱的相似程度来判断信号的平稳性。
按照文献[4]中的检验方法,定义不同时间点上的频谱与频谱平均值的距离cn(x)为:
2用户界面生成
2.1 Matlab中用户界面的生成
Matlab 为用户设计图形界面提供了一个高效、方便的集成环境。在Matlab中,基本的图形对象主要包括坐标轴、控件、下拉菜单和内容菜单。用户可以通过这些对象设计出界面友好,功能强大,操作简单的图形用户界面。图形用户界面的生成主要分为以下几个步骤:
(1)规划所设计的图形用户界面,主要包括:确定需要哪些窗口,每个窗口怎样布局,窗口中的各个对象各有什么功能,对象之间如何配合工作,以及相应的异常处理;
(2)在 Matlab提示行下输入GUIDE,载人用户界面开发环境;
(3)利用Layout Editor,完成用户面板以及界面的制作,并对相应的按钮及控件属性进行适当的设置;
(4)在Programme Editot中编辑各个对象的回调函数,实现各个对象韵具体功能;
(5)利用Mfile编译器生成客户端,完成随机数据仿真系统的设计。
2.2 用户界面介绍
本文所实现的用户界面主要包括两个窗口,分别是主窗口和数据生成窗口。由于Matlab对保存绘图区域有限制,因此设计时没有在窗口中设置固定的绘图区域。窗口只相当于一个命令菜单,所有的绘图将会以独立窗口的形式根据用户需求动态产生。这样便于用户对比和保存图片。下面对主要窗口分别加以介绍。
2.2.1 主窗口介绍
主窗口如图1所示。主窗口主要用于绘制原数据和替代数据的各种波形以及显示平稳性检验结果。其中,“Create/Open Original Data”按钮用来打开数据生成窗口。
“View/Change Parameters”按钮用来查看或改变当前仿真参数,它在原始数据存在的情况下才有效。主要的仿真参数有:
“Time Scale of TFD”和“Frequency Scale of TFD”用于确定绘制TFD图片时的时间/频率轴采样周期,由于计算和显示时频分布图比较耗费时间,将采样周期设大,可以提高速度,但是相应的时频分布图的分辨率会下降。
“Max Level of Hermite Function”用于确定求TFD时所使用Hermite函数的最高阶数。最高阶数越高,则分辨率越高,但是相应的计算时间会加长。
“The Number of Surrogates”用于确定平稳性检验时所用的参考替代数据个数。个数越多,则检验结果越精确,但是会极大地延长计算时间。
“Create Surrogate”按钮用于产生替代数据,其在原始数据存在的情况下才有效。由于替代数据具有随机性,因此用户可以多次产生不同的替代数据,观察它们的性质。
右上方的下拉菜单用于选择需要绘图或者保存数据的对象,主要包括原始数据的时域、频域和时频域图,替代数据的时域、频域和时频域图,以及替代数据的平稳度分布。它在原始数据存在的情况下才有效。
“Show Selected Plot”按钮用于在新窗口中绘制下拉菜单所选图线,它在原始数据存在的情况下才有效。
“Save Selected Data”按钮用于保存下拉菜单所选图线对应的数据,它在原始数据存在的情况下才有效。
2.2.2 数据生成窗口
数据生成窗口如图2所示。数据生成窗口主要用来产生实验用数据或者打开已经存在的实验数据。
“Creat Original Data”按钮用来产生测试用数据。按下此按钮后会提示输入产生数据用的参数。由于数据是通过公式:
产生的调频信号,因此需要确定参数P1和P2,另外还要确定t的区间和采样周期。数据成功产生后会在新建窗口中显示该数据时域波形。如果当前存在数据波形,将会覆盖它。
“Open Original Data”按钮用来打开已经存在的数据文件。选择好文件后会提示输人参数。主要包括读人数据的时间起点、时间采样周期和数据长度。数据成功读入后会在新建窗口中显示该数据时域波形。如果当前存在数据波形,将会覆盖它。
“Confirm”按钮用于确认新建窗口显示的数据就是用户想要的数据,并返回主窗口。它在创建或打开的数据存在的情况下才有效。
3数据仿真和分析
将实验数据取为调频信号 x(t)=sin(sin(t/8)πt)。t起始为0,采样周期为0.1 s,数据长度为400个点,其时域波形和频域波形如图3所示。
由式(6)产生的替代数据s(t)的时域波形和频域波形如图4所示。
由图3和图4不难看出,替代数据与原数据的傅里叶变换幅值相同,但替代数据傅里叶变换的相位是随机的。
图5显示了由式(7)计算得到的原数据和替代数据的时频分布图。由图5中可见,原数据的时频分布图有明显的结构性。它表明是非平稳的,而替代数据的时频分布图的结构性较原数据有明显减弱,表明替代数据的平稳性增加。
由式(11)计算得到的0的概率密度函数f(0) 如图6所示,其中一共计算了1 000次替代数据。
由图6可见,替代数据的平稳度主要分布在0.02附近。数据的平稳度落在0~0.04 之间可以认为是平稳的,而在此之外可以认为是非平稳的。
由式(11)计算得到的1=0.046。位于上述区间之外,因此判为非平稳。这一结果也与图5所示的结果相吻合。
4结 语
利用替代数据法和时频透视法,并采用Matlab 的GUI开发环境,设计了一个信号平稳性检验系统。该系统的用户界面友好。利用该系统可以观察信号及其替代数据的频域和时频域波形,检验信号的平稳性。
matlab中平稳性检验,基于Matlab的信号平稳性检验系统相关推荐
- matlab功率谱的峰值,基于Matlab功率谱密度估计方法.doc
基于Matlab功率谱密度估计方法 基于Matlab功率谱密度估计方法要在实际情况下, 许多平稳信号无法导出数学表达式, 要准确获取这些信号的功率谱密度存在一定的困难.根据维纳辛钦 (Wiener K ...
- Matlab之DNN:基于Matlab利用神经网络模型(epochs=10000000)预测勒布朗詹姆斯的2018年总决赛(骑士VS勇士)第一场得分、篮板、助攻
Matlab之DNN:基于Matlab利用神经网络模型(epochs=10000000)预测勒布朗詹姆斯的2018年总决赛(骑士VS勇士)第一场得分.篮板.助攻 目录 输出结果 预测勒布朗詹姆斯的20 ...
- 基于matlab的漏电保护器仿真,基于MATLAB仿真的高压断路器振动信号故障诊断
内容介绍 原文档由会员 第二波打卡 发布 基于MATLAB仿真的高压断路器振动信号故障诊断 16300字 44页 原创作品,已通过查重系统 摘 要 断路器在分合闸操作中所产生的振动信号是一系列无规则的 ...
- matlab模拟频谱仪,毕业设计基于MATLAB的声音信号频谱分析仪源代码V3.1(手机版)...
<[毕业设计]基于MATLAB的声音信号频谱分析仪源代码.doc>由会员分享,可免费在线阅读全文,更多与<毕业设计基于MATLAB的声音信号频谱分析仪源代码(V3.1)>相关文 ...
- 基于matlab的数字调制,基于MATLAB的多功能数字调制系统信号源仿真
内容简介: 毕业设计 基于MATLAB的多功能数字调制系统信号源仿真,共60页,27225字 摘要 数字通信与模拟通信系统相比有着灵活性.高效性和保密性等突出特点.数字通信系统仿真对于系统分析起着越来 ...
- wcdma matlab程序,一种基于MATLAB的WCDMA信号仿真方法.doc
目 录 1 引言 ----------------------- 6 2 WCDMA系统 ------------------- 7 2.1 WCDMA通信系统的模型与结构 ----------- 7 ...
- matlab 双轴图,基于MATLAB图像识别的双轴寻光系统的制作方法
本发明涉及新能源与图像处理与识别技术领域,具体涉及一种基于MATLAB图像处理的双轴寻光系统. 背景技术: 随着现在新能源技术的发展,如何高效地开发和利用好太阳能成为人们在新能源领域探索的焦点之一.光 ...
- 基于matlab fdma传输系统设计,基于MATLAB的LTE系统仿真研究
摘 要: 根据LTE系统的原理和模块构成,分析并建立了一个基于MATLAB的系统级仿真平台.针对LTE系统中被广泛认可的轮询.比例公平和最大载干比三种经典调度算法,利用仿真平台产生的数据对三者的性能 ...
- matlab自适应逆控制,基于matlab仿真模块的自适应有源噪声逆控制研究
ATE 农业技术与装备 X2 c(q-1) p(q-1) u2 C( x) p(q-1) LMS M(q-1) p(q-1) w( q-1) w( q-1) ! ! e2 d2 + - - + p(q ...
- 用matlab编程简单电子琴,基于MATLAB的数字信号发生器及简易电子琴设计论文.doc...
基于MATLAB的数字信号发生器及简易电子琴设计论文 基于MATLAB的数字信号发生器及简易电子琴设计 摘 要 数字信号发生器是一种基于软硬件结合实现的函数波形产生仪器.在工程实践中需要检测和分析的各 ...
最新文章
- ios app上架App Store需要多少费用?
- 两个计算机系统安装,如何在一台电脑上同时重装两个系统|戴尔电脑怎么安装两个系统...
- 关于Vmware下NAT模式物理机无法ping通虚拟机但是可以上网的解决方法
- SAP Spartacus ComponentData的提前subscription
- 深入理解Java ClassLoader及在 JavaAgent 中的应用
- 外呼机器人起名_智能外呼机器人,目前都有哪些公司做产品?
- Spring 注解AOP 入门
- mssql访问 oracle
- 基于智能网卡的vhost-user优化接收方法-学习笔记
- FFmpeg源代码简单分析:av_write_trailer()
- java中sql 语句格式转换_Java字符串转Sql语句In格式工具方法
- 利用IDLE对 dem进行批量拼接处理
- wps 甘特图_强烈推荐5款在线甘特图工具,又简单又清爽又好用
- 2022南京商业贷款提前还款
- javascript---对象和函数的引用、浅拷贝、深拷贝、递归
- 【Linux】Linux关闭防火墙、关机重启和查看系统运行级别
- Android界面编程之简单的图片浏览器
- 【时间函数】gettimeofday
- 极验滑块轨迹生成 3、4通用
- 桥牌记分规则 双人赛计分表