轴承过盈量程序matlab,轴承matlab处置程序.doc
轴承matlab处置程序
1. 数据导入matlab
1.1启动Matlab软件
1.2点击载入故障数据中的G2015,Workspace窗口出现:
1.3取第一组数据G201,命令窗口输入:
G201=G2015(1:1:20000);
2. 数据预处理
在测试中由数据采集所得的原始信号,在分析前需要进行预处理,以提高数据的可靠性和真实性,并检查信号的随机性,以便正确地选择分析处理方法。预处理工作主要包括三个方面:一是除去信号中的外界干扰信号和剔除异常数据,如趋势项和异点;二是对原始数据进行适当的平滑或拟合;三是对原始信号的特性进行检验。当然这些处理工作不是全部必需的,可以选—项或两项内容,当认为原始信号获取工作十分可靠或原始数据简单可以直接判断的情况下,也可以不进行这些预处理工作。以下所做数据预处理,故障轴承以G201为例,正常轴承以Z201为例,观察原始数据经过不同方法做处理前后的变化。
1.1零均值化处理(原理公式见报告P8)
命令窗口输入:
G201l=G201-sum(G201)/20000;%G201l为零均值处理后的数据。 “20000”为采样点数。sum为求和语句
subplot(2,1,1),plot(G201);subplot(2,1,2),plot(G201l);%显示G201与G201l
得到下面图形:
从时域图形上看,是波形整体在Y轴的平移。再看看频域变化,命令窗口输入:
N=20000; %采样点数
fs=10000; %采样频率
f=(0:N-1)'*fs/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)); %显示G201与G201p的频谱图
得到下面图形:
从频域图可以明显看出,零均值后消除处出现一个由直流分量产生的大谱峰(将近达到),处理后避免了其对周围小峰值产生的负面影响,便于频域分析。
1.2 消除趋势项(原理公式见报告P10)
使用最小二乘法,命令窗口输入:
t=(0:1/fs:(N-1)/fs)'; %离散时间列向量
G201x=polyfit(t,G201,6); %计算多项式待定系数向量
G201x=G201-polyval(G201x,t); %用G201减去多项式系数生成的趋势项,G201x即为消除趋势项后的数据
subplot(2,1,1),plot(G201);subplot(2,1,2),plot(G201x);%显示G201与G201x
得到以下图形:
与前面零均值化处理中做频域图的方法一样,做出G201与G201x的频谱图G201p与G201xp,得到图形如下:
从时域图形和频域图形上看,消除趋势项与零均值化处理的功能相似。不过,需要注意的是,它更重要的消除趋势项,因为本数据中的多项式趋势项很小,所以没有明显的变化。
1.3 平滑处理(原理公式见报告P11)
使用五点三次平滑,命令窗口输入:
a=G201';
for k=1:2
b(1)=(69*a(1)+4*(a(2)+a(4))-6*a(3)-a(5))/70;
b(2)=(2*(a(1)+a(5))+27*a(2)+12*a(3)-8*a(4))/35;
for j=3:N-2
b(j)=(-3*(a(j-2)+a(j+2))+12*(a(j-1)+a(j+1))+17*a(j))/35;
end
b(N-1)=(2*(a(N)+a(N-4))+27*a(N-1)+12*a(N-2)-8*a(N-3))/35;
b(N)=(69*a(N)+4*(a(N-1)+a(N-3))-6*a(N-2)-a(N-4))/70;
a=b;
end
G201ph=a'; %G201ph为五点三次平滑法处理的数据
subplot(2,1,1),plot(G201);subplot(2,1,2),plot(G201ph);%显示G201与G201ph
得到以下图形:
与前面零均值化处理中做频域图的方法一样,做出G201与G201ph的频谱图G201p与G201php,得到图形如下:
从时域图形上看,平滑处理使图形变得平滑,去除毛刺,从频域图形上看,高频部分明显变少变小,而低频部分基本无变化。因为故障的频率主要集中在低中频部分,这样处理后不仅对故障的分析无影响,而且去除部分噪音,减少干扰。
1.4 滤
轴承过盈量程序matlab,轴承matlab处置程序.doc相关推荐
- matlab控制算法C语言,PID算法Matlab仿真程序和C程序
<PID算法Matlab仿真程序和C程序>由会员分享,可在线阅读,更多相关<PID算法Matlab仿真程序和C程序(6页珍藏版)>请在人人文库网上搜索. 1.增量式PID控制算 ...
- matlab水印剪切攻击程序,可以运行的水印matlab程序(嵌入,提取,攻击测试等).doc
可以运行的水印matlab程序(嵌入,提取,攻击测试等).doc IMREAD YUANSHITUXIANGBMP SUBPLOT2,2,1IMSHOWTITLE 原始图像 WATERIMREAD W ...
- matlab绘制频散曲线,Matlab绘制频散曲线程序代码.docx
Matlab绘制频散曲线程序代码.docx 下载提示(请认真阅读)1.请仔细阅读文档,确保文档完整性,对于不预览.不比对内容而直接下载带来的问题本站不予受理. 2.下载的文档,不会出现我们的网址水印. ...
- c 传图片数据给matlab,c++ - 如何通过UDP将数据从C ++应用程序发送到Matlab并进行绘制 - 堆栈内存溢出...
我想使用Winsock通过UDP将数据从我的C ++应用程序发送到Matlab,并实时绘制数字. 例如:我得到了从1到10的数字,我想从C ++发送1,在Matlab中接收它,将其放在图形上(图,条等 ...
- matlab有意思程序,matlab有意思的小程序
10个C++趣味小程序,很有意思的.VIP专享文档 VIP专享文档是百度文库认... 现在很多人使用微信的时间已经非常长了,他们注册的微信号往上可能已经是5年前的事情了,正是由于不少使用者在这个过程当 ...
- matlab的exec程序,C++调用Matlab画图的一段程序
劳动节闲来无事,写了一天程序,just for fun. 看,这是C++调用Matlab画图的一段程序.暂时不想多解释了,有兴趣的话,看看下面的代码吧. 以下几段代码由上到下,越来越旧.最上面的是最新 ...
- dmc matlab程序,matlab编的DMC程序.doc
matlab编的DMC程序 clear all; % close all; %系统模型建立 num=[0.8]; den=[225 1]; [a,b,c,d]=tf2ss(num,den); % st ...
- matlab 万年历,matlab时钟万年历 matlab时钟万年历(附word说明文档):
程序主要 联合开发网 - pudn.com...
matlab时钟万年历 所属分类:图形图像处理 开发工具:matlab 文件大小:54KB 下载次数:6 上传日期:2018-08-24 21:54:23 上 传 者:小小葡萄123 说明: mat ...
- 精馏塔matlab,MATLAB图解精馏塔理论塔板数程序代码
<MATLAB图解精馏塔理论塔板数程序代码>由会员分享,可在线阅读,更多相关<MATLAB图解精馏塔理论塔板数程序代码(6页珍藏版)>请在人人文库网上搜索. 1.MATLAB图 ...
- 车牌识别与计算机编程,基于MATLAB的车牌识别程序详解.ppt
基于MATLAB的车牌识别程序详解 自定义一个字符函数,用来从车牌区域中提取出7个字符,其中利用切割函数来进行切割. 程序:function [word,result]=getword(d) word ...
最新文章
- python 命令-python常见命令
- 《深入理解 Java 内存模型》读书笔记(上)(干货,万字长文)
- 在阿里,我们这样帮助用户实现业务云原生化迁云
- 如何在CDH5上部署Dolphin Scheduler 1.3.1
- 产品经理必知的三个数据分析框架
- python矩阵乘法分治_分治法实现矩阵乘法
- 牛津美女硕士放弃百万年薪,用废弃物做轻奢包,马斯克妈妈也来捧场
- mysql查询的是问号_MySQL客户端查询中文显示为问号(linux) | 学步园
- 跑步进入全站 HTTPS ,这些经验值得你看看
- android 随身无线网卡,让小锐WiFi支持USB无线网卡/随身WiFi(附各种“随身wifi”芯片型号)...
- 小米便签源码分析——UI包
- Android Studio 报错:Incompatible types
- 2022年第十二届APMCM亚太地区数学建模竞赛1月加赛E题翻译以及思路
- APP遇到大量的真实手机号刷注册用户该如何应对?
- js如何截取某个字符串前面所有的字符串
- 详解“因果效应估计”
- 小程序学习:电影周周看
- 零代码—0代码—无代码平台是什么?怎么选?
- 了解一下chirp,扫频信号的生成
- 服务器怎么当虚拟主机,服务器怎么做虚拟主机
热门文章
- 论文阅读:2021 | CE-FPN: Enhancing Channel Information for Object Detection
- 面试题大全2-技能提升-项目部分-运维
- linux simg2img,simg2img工具
- 医疗时鲜资讯:移动医疗 or 互联网医疗 or 远程医疗?
- apt-get安装软件:依赖冲突问题及解决
- python 使用 reportlab 生成 pdf
- 计算机怎么删除表格,怎么快速删除电脑word文档中不想要的表格
- 时频分析matlab实例,Matlab时频分析工具箱及仿真案例
- 切比雪夫多项式(Chebyshev Polynomials)
- 从认知盈余说起,也谈分享精神