matlab 时域同步平均,时域平均法解析.doc
时域平均法
图2所示是截取不同的段数N进行时域平均的效果。由图可见,虽然原来信号(N=1)的信噪比很低,SNR=0.5,但经过多段平均后,信噪比大大提高。当N=256段时,可以得到几乎接近理想的正弦信号。而原始信号中的正弦分量,几乎完全被其它信号和随机噪声所淹没。
对于旋转机械或往复机械运行中所产生的带有周期重复性的机械信号x(t),如果x(t)由周期信号f(t)和白噪声n(t)组成,即
(1)
以f(t)的周期T去截取信号x(t),共截得N段,然后将各段对应点相加,由于白噪声的不相关性,可得到
(2)
再对x(ti)平均,便得到输出信号y(ti)
(3)
此时输出的白噪声是原来输入信号x(t)
中的白噪声的,因此信噪比提高
倍。下面是对此结论的分析。
采用时域平均法提取周期分量f(t),
也可看作是梳状滤波器进行滤波的过
程。以为时间间隔对),n=0,1,2,…
…。设N为叠加、平均的周期段数目,
M为一个周期中的采样数目,为采样间隔。x)为滤波器的输入,则
滤波器的输出y(n)为
,),x(n)的z变换为Y(z),有
(5)
滤波器的传递函数H(z)为
(6)
令,且有周期==,传递函数
(7)
传递函数H(z)的相位为
= (8)
当时,在的值可由洛必达法则求得: (9)
可见时域平均过程等价于具有中心频率,也就是的梳状滤波器。图是滤波器关于频率的幅值响应曲线的形状。由
曲线形状。,的值称为滤波器主瓣峰值,此时= (10)
通过滤波器后的噪音能量的改变,即ENB为
(11)
由于时域平均时,截取信号段的周期为T,即频率为f0,所以,的范围为, (x=)
(12)
通过时域平均,从方差意义上讲,使平均后的噪音缩小了N倍,相当于输出的噪音是输入噪音的1/。因此梳状滤波器抑制了白噪音,提高了信噪比。
没有经过时域平均的振动信号如图4所示,经过时域平均的振动信号的时域频谱如图5所示。
经过对比可以发现,没有经过时域平均的振动信号的相位不够明显,里面包含很多的杂质,对这种信号进行分析处理后所得结果不够准确。经过时域平均时后信号中的与基频无关的信号会被过滤掉。图中所显示的谱图相位明显清晰,符合实验要求。
把两种时域波形同时进行FFT分析处理,得到结果如图6所示。
图中上曲线表示没有进行时域平均的信号经过快速傅里叶变换后的结果,下曲线表示振动信号经过时域平均后做FFT 处理得到的结果。我们可以观察到是否经过时域平均对信号的基频没有任何影响,而在其它区域,上曲线的值都要比下曲线高出很多,也就是说,如果有故障发生,其振动特点更能在时域同步平均后的曲线中被发现,而如果没有进行时域同步平均处理,故障就很容易被较高频率的噪声淹没掉。所以通过这个实例说明经过时域平均处理能够很好地去除高频的干扰和噪声对整体的影响。
3.误差分析
离散序列x(n)时域平均时的截取长度M,由其中人们感兴趣的周期分量的周期T和采样间隔Δt确定,确切地说当M取T /Δt时就近取整值。感兴趣分量的周期T在实际中有办法事先精确确定,一般情况下总有M≠T/Δt ,即总有截取误差ΔT= T- T0 (T0=ΜΔt )存在。ΔT的存在,对用时域平均方法提取信号中感兴趣的周期分量及倍频分量将有影响。
下面首先讨论时域平均方法抑制非感兴趣周期分量的原理。在此基础上,探讨截取误差ΔT对时域平均结果的影响及合理平均段数的确定。
式( 1)为离散序列的时域平均公式,离散信号与连续信号时域平均的本质一样,为了方便起见,这里用连续信号来研究。设连续信号x(t)中含有非感兴趣的周期分量xi(t) =sin(2πfit+ Qi ) ,其周期记为Ti ,角频率记为ki ,现把信号x(t)按其中的感兴趣周期分量的周期T (角频率记为k)连续截取N 段进行时域平均处理,则分量xi(t)平均后结果为
(13)
令,,则 (14)
其中
matlab 时域同步平均,时域平均法解析.doc相关推荐
- 时域同步平均(TSA)降噪原理
项目说明 旋转机械设备出现故障损伤时,振动信号中通常包含有规律性的故障特征,其提取对于设备故障诊断具有重要意义.但由于噪声干扰.传递路径衰减.多振源耦合等因素,往往需要对振动信号进行降噪处理,如滤波器 ...
- Python 信号处理——时域同步平均技术(TSA)
文章目录 前言 一.TSA实现方式 二.使用步骤 1.主要代码 2.示例Demo 总结 前言 时域同步平均技术是一种信号降噪的手段,前提假设是在某个测量节点上(暂且这么命名)主频信号的幅值随时间是不变 ...
- 时域同步平均TSA-学习笔记
谈谈自己实际工作中对于旋转机械故障诊断的感受,大家一起交流见解.信号处理技术可以做的很深,有些方法很简单,但很实用,工作中大部分还是基于信号的前处理.降噪.频域技术分析来进行故障诊断(特征频域或者时域 ...
- MATLAB观察系统响应,基于MATLAB的系统的时域响应分析精选.doc
基于MATLAB的系统的时域响应分析精选 摘 要 系统的时域响应分析是衡量控制系统的重要指标.所谓控制系统的时域响应分析就是在时间域内求解系统的微分方程,然后根据绘制出来的曲线分析系统的性能和各主要参 ...
- matlab 线性系统响应,MATLAB实验二线性系统时域响应分析.doc
MATLAB实验二线性系统时域响应分析 武汉工程大学实验报告 专业 班号 组别 01 教师 姓名 同组者 ( 个人 ) 实验名称 实验二 线性系统时域响应分析 实验日期 2011-11-24 第 2 ...
- 《数字语音处理》- 1. 基于MATLAB的语音信号时域特征分析
声明 本文仅在CSDN发布,其他均为盗版.请支持正版! 正版链接: https://blog.csdn.net/meenr/article/details/117629753 基于MATLAB的语音信 ...
- matlab施加阶跃相应,(基于MATLAB的系统的时域响应分析)(最终版)最新版
<基于MATLAB的系统的时域响应分析.doc>由会员分享,可免费在线阅读全文,更多与<(基于MATLAB的系统的时域响应分析)(最终版)>相关文档资源请在帮帮文库(www.w ...
- 用matlab验证时域卷积,用MATLAB实现信号的时域卷积
用MATLAB实现信号的时域卷积 用 MATLAB 实现信号的时域卷积摘要:利用抽样方法绘制连续信号和离散信号.编制函数计算连续时间信号,离散时间信号的卷积.关键字:MATLAB:离散卷积:连续卷积实 ...
- 频域串联滞后校正matlab,自动控制原理实验七基于MATLAB控制系统频域法串联校正设计.doc...
自动控制原理实验七基于MATLAB控制系统频域法串联校正设计.doc 实验七 基于MATLAB控制系统频域法串联校正设计 一.实验目的 (1)对给定系统设计满足频域或时域指标的串联校正装置: (2)掌 ...
最新文章
- SQL SERVER全面优化-------写出好语句是习惯
- 解决Uncaught SyntaxError: Unexpected token var报错问题
- java 注解 属性 类型_跟光磊学Java开发-Java注解
- SpringMVC框架、Spring boot框架、SSM區別
- 在SQL SERVER里面用命令查包含某字段的表
- fiddler限速/弱网模拟
- My PaintBrush Pro for mac(专业的绘图画板)
- struts2拦截器定义
- Python生信练习
- 软件工程阶段性总结(一)——概述
- 目标检测YOLO实战应用案例100讲-基于深度学习的显著性目标检测研究与应用(论文篇)
- ATM维护人员教大家正确使用银行卡和取款机
- python工程师怎么考-【一个合格的Python工程师需要达到怎样的编程水平】
- 群辉nas虚拟linux,UNRAID教程:3分钟 用unraid自带的虚拟机 安装 黑群晖NAS DSM系统 很强大!...
- Centos在NAT模式下的设置
- composer windows 下面内存溢出 Allowed memory size of 1610612736 bytes exhausted
- Android Jetpack架构组件之Navigation
- 仿抖音--音乐裁剪控件,android实现IOS版本效果
- 两个复数相加减c语言编程,用c语言 如何编写两个复数的运算啊 都含有虚部 谢谢啊...
- 【MySql】MySql事务隔离级别