【7A版】轴承matlab处理程序.doc

【7A版】轴承matlab处理程序.doc

ID:39780885

大小:1.02 MB

页数:20页

时间:2019-07-11

【7A版】轴承matlab处理程序.doc_第1页
【7A版】轴承matlab处理程序.doc_第2页
【7A版】轴承matlab处理程序.doc_第3页
【7A版】轴承matlab处理程序.doc_第4页
【7A版】轴承matlab处理程序.doc_第5页
资源描述:

《【7A版】轴承matlab处理程序.doc》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库

1、【MeiWei_81重点借鉴文档】1.数据导入matlab1.1启动Matlab软件1.2点击载入故障数据中的G2015,Workspace窗口出现:1.3取第一组数据G201,命令窗口输入:G201=G2015(1:1:20RR0);2.数据预处理在测试中由数据采集所得的原始信号,在分析前需要进行预处理,以提高数据的可靠性和真实性,并检查信号的随机性,以便正确地选择分析处理方法。预处理工作主要包括三个方面:一是除去信号中的外界干扰信号和剔除异常数据,如趋势项和异点;二是对原始数据进行适当的平滑或拟合;三是对原始信号的特性进行检验。当然这些处理工作不是全

2、部必需的,可以选—项或两项内容,当认为原始信号获取工作十分可靠或原始数据简单可以直接判断的情况下,也可以不进行这些预处理工作。以下所做数据预处理,故障轴承以G201为例,正常轴承以Z201为例,观察原始数据经过不同方法做处理前后的变化。1.1零均值化处理(原理公式见报告P8)命令窗口输入:G201l=G201-sum(G201)/20RR0;%G201l为零均值处理后的数据。“20RR0”为采样点数。sum为求和语句subplot(2,1,1),plot(G201);subplot(2,1,2),plot(G201l);%显示G201与G201l得到下面

3、图形:从时域图形上看,是波形整体在R轴的平移。再看看频域变化,命令窗口输入:N=20RR0;%采样点数fs=10000;%采样频率f=(0:N-1)'Rfs/N;%进行对应的频率转换G201p=abs(fft(G201));%进行fft变换,G201p为G201进行fft变换后结果G201lp=abs(fft(G201l));%进行fft变换,G201lp为G201l进行fft变换后结果subplot(2,1,1),plot(f(1:N/2),G201p(1:N/2));subplot(2,1,2),plot(f(1:N/2),G201lp(1:N/2)

4、);%显示G201与G201p的频谱图得到下面图形:【MeiWei_81重点借鉴文档】【MeiWei_81重点借鉴文档】从频域图可以明显看出,零均值后消除处出现一个由直流分量产生的大谱峰(将近达到),处理后避免了其对周围小峰值产生的负面影响,便于频域分析。1.2消除趋势项(原理公式见报告P10)使用最小二乘法,命令窗口输入:t=(0:1/fs:(N-1)/fs)';%离散时间列向量G201R=polRfit(t,G201,6);%计算多项式待定系数向量G201R=G201-polRval(G201R,t);%用G201减去多项式系数生成的趋势项,G201

5、R即为消除趋势项后的数据subplot(2,1,1),plot(G201);subplot(2,1,2),plot(G201R);%显示G201与G201R得到以下图形:与前面零均值化处理中做频域图的方法一样,做出G201与G201R的频谱图G201p与G201Rp,得到图形如下:【MeiWei_81重点借鉴文档】【MeiWei_81重点借鉴文档】从时域图形和频域图形上看,消除趋势项与零均值化处理的功能相似。不过,需要注意的是,它更重要的消除趋势项,因为本数据中的多项式趋势项很小,所以没有明显的变化。1.3平滑处理(原理公式见报告P11)使用五点三次平滑

6、,命令窗口输入:a=G201';fork=1:2b(1)=(69Ra(1)+4R(a(2)+a(4))-6Ra(3)-a(5))/70;b(2)=(2R(a(1)+a(5))+27Ra(2)+12Ra(3)-8Ra(4))/35;forj=3:N-2b(j)=(-3R(a(j-2)+a(j+2))+12R(a(j-1)+a(j+1))+17Ra(j))/35;endb(N-1)=(2R(a(N)+a(N-4))+27Ra(N-1)+12Ra(N-2)-8Ra(N-3))/35;b(N)=(69Ra(N)+4R(a(N-1)+a(N-3))-6Ra(N-2

7、)-a(N-4))/70;a=b;endG201ph=a';%G201ph为五点三次平滑法处理的数据subplot(2,1,1),plot(G201);subplot(2,1,2),plot(G201ph);%显示G201与G201ph得到以下图形:【MeiWei_81重点借鉴文档】【MeiWei_81重点借鉴文档】与前面零均值化处理中做频域图的方法一样,做出G201与G201ph的频谱图G201p与G201php,得到图形如下:从时域图形上看,平滑处理使图形变得平滑,去除毛刺,从频域图形上看,高频部分明显变少变小,而低频部分基本无变化。因为故障的频率主

8、要集中在低中频部分,这样处理后不仅对故障的分析无影响,而且去除部分噪音,减少干扰

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。