- 博客(21)
- 资源 (8)
- 收藏
- 关注
原创 卷积神经网络西储大学轴承故障诊断(基于时频变换)
网上有不少方法,本人尝试了下. 一、选取西储大学轴承数据 二、利用短时傅立叶或小波变换为时频图,选取了cmor小波的方法,得到640张时频图如下 三、直接利用alexnet,做分类实验 得到结果 感觉总体正确率很高,可分性较强。 ...
2020-04-09 10:34:16
11491
129
转载 轴承故障诊断之迁移学习
本文来自本人公众号:** 轴承故障诊断与寿命预测 ** 迁移学习有人比做人工智能的最后一公里。 迁移学习的通俗解释为把一个已经解决一定问题的模型调整为解决近似问题的模型。 迁移学习的核心在于问题的相似性度量和调整,无论是数据迁移还是模型迁移等,核心在于问题之间的联系性,通过找到关联转移的方法,实现不同目的。比如轴承故障诊断过程中,可能工况具有差别,导致振动信号采集过程中的幅度具有不同,但本质的特征频率不会因为信号幅度的大小而改变,这就属于不同工况间的内在联系,通过抓住本质的过程,实现模型和方法解决不同的问题
2020-09-02 10:14:32
5790
8
原创 轴承故障诊断之时域指标
传统的轴承故障诊断方法有时域,频域,时频域等方法,从今天开始写一下这些传统方法。 时域指标作为最传统的信号分析方法,有独特的使用优势,同样采样数据驱动的机器学习方法时,也可以利用这些统计量作为样本进行判断学习。 ...
2020-05-10 12:22:32
5989
5
原创 研华PCI1716L的C#编程
新建windows窗体项目; 添加引用研华的库文件 将库添加using 整个代码如下: using System; using System.Collections.Generic; using System.ComponentModel; using System.Data; using System.Drawing; using System.Linq; using System.Text...
2019-10-24 09:27:03
3841
1
原创 matlab判断IIR传递函数稳定性
num = input('输入分子: '); den = input('输入分母: '); N = max(length(num),length(den)); x = 1; y0 = 0; S = 0;zi = zeros(1,N-1); for n = 1:1000 [y] = filter(num,den,x); if abs(y) < 0.000001, break, ...
2019-09-22 15:50:20
2643
原创 DFT与DTFT区别
k = 0:15; w = 0:511; x = cos(2*pi*k*3/16); X = fft(x); XE = fft(x,512); plot(k/16,abs(X),'o', w/512,abs(XE)) xlabel('\omega/\pi'); ylabel('Magnitude')
2019-09-22 09:45:07
596
原创 中值滤波频响
h1 = ones(1,5)/5; h2 = ones(1,14)/14; [H1,w] = freqz(h1, 1, 256); [H2,w] = freqz(h2, 1, 256); m1 = abs(H1); m2 = abs(H2); subplot(211) plot(w/pi,m1,'r-',w/pi,m2,'b--'); ylabel('Magnitude'); xlabel('\o...
2019-09-21 21:54:09
351
原创 中值滤波
R = 50; a = rand(1,R)-0.4; b = round(a); m = 0:R-1; s = 2*m.*(0.9.^m); x = s + b; y = medfilt1(x,3); subplot(2,1,1) stem(m,x);axis([0 50 -1 8]); xlabel('n');ylabel('Amplitude'); title('Impulse Nois...
2019-09-21 20:37:19
244
原创 MATLAB滑动平均滤波
示例 R = 50; d = rand(R,1)-0.5; m = 0:1:R-1; s = 2*m.*(0.9.^m); x = s + d'; M = input('输入滤波重复采样数 :'); subplot(211) plot(m,d,'r-',m,s,'b--',m,x,'k--') xlabel('Time index n'); ylabel('Amplitude') legend(...
2019-09-21 20:23:03
3052
原创 matlab互相关运算
求两个数组相关性 x = input('输入数组1 '); y = input('输入数组2 '); n1 = length(y)-1; n2 = length(x)-1; r = conv(x,fliplr(y)); k = (-n1):n2'; subplot(211) stem(k,r); xlabel('Lag index'); ylabel('Amplitude'); r1=xcorr(...
2019-09-21 17:55:18
8213
1
原创 matlab平均值消噪声
原始信号s=0.9(n(0.9^n) 加入随机噪声rand(n) 利用多次采样取平均值消除噪声,代码如下: clc clear all close all r=50; M=0:r-1; s=0.9*M.*(0.9.^M); k=rand(r,1)-0.5; s1=s+k'; subplot(311) stem(M,s),title('原始信号') subplot(312) stem(M,s1),t...
2019-09-21 10:53:26
1877
原创 matlab中短时傅里叶变换tfrstft函数用法
[TFR,T,F]=TFRSTFT(X,T,N,H,TRACE) X : 信号. T : 时间序列 (默认值 :1:length(X)). N : 频率点数 (默认值: length(X)). H : 选择的平滑窗函数 (默认值: Hamming(N/4)). TRACE : 如果不为零,跟踪计算过程 (默认为 : 0). TFR : 得到的时频谱值(复数) F : 归一化...
2019-09-18 20:25:18
16792
3
原创 matlab时频分析工具箱tftb-0.2详细安装教程
通俗点讲吧,自己装了一遍,截图讲下: 首先下载安装包 解压,将解压后文件夹拷到D:\Program Files\MATLAB\R2016b\toolbox,具体根据个人安装路径修改。 然后打开matlab,找到主页里设置路径点击 点击添加并包含子文件夹按钮 选择刚才拷到目录下的文件夹 保存关闭,点击设置路径上面预设按钮 点击常规里--------更新工具箱路径缓存,保存完关闭即可 ...
2019-09-18 19:20:18
8773
3
原创 matlab中wextend函数简单示例理解
学习短时傅里叶变换,发现对信号进行延拓时用到wextend函数 学习了一下: 初始函数COS函数一个周期,-pi到pi 进行sym延拓代码如下: clc clear all t=-pi:pi/240:pi t2=-2*pi:pi/240:2*pi+pi/120 z=cos(t) subplot(211) plot(t,z) z2=wextend(1,'sym',z,round(length(z)/...
2019-09-18 17:10:23
5856
原创 python模块time获取系统时间及暂停函数
import time for i in range(100): print(time.strftime('%Y-%m-%d %H:%M:%S',time.localtime(time.time()))) time.sleep(1) 可以按照格式得到系统时间
2019-09-16 10:27:12
1186
原创 接上一篇fft画图乱直线,频谱图横坐标频率设置问题
上一篇画出图: 明明函数只有三个频率为啥六条线,其实是横坐标设置问题及频谱搬移问题, 因为plot(k)没有x,只有一个变量,那x应该怎么取呢,freq=(0:N/2)*fs/N;这样,取采样点的一半作图,画的时候这样plot(freq,abs(k(1:N/2+1))) 再看看图,代码改为下面这样 close all; clear all; clc; %采样频率 fs=1024; %采样间...
2019-09-07 15:51:41
1536
原创 fft画图出现乱直线情况
信号处理小白肯定开始用matlab换图时,y=f(x),上来就k=fft(y)和plot(k)了, 但画出图来发现,,,,,,,,,, **What is are?**遇到的同学举个手。。。。。。。 这其实是横坐标t的取值以及abs问题,并不是FFT的问题,当t取值不当及abs没加时,确实会出现直线乱跳的图,那时间怎么取呢? %采样频率 fs=1024;/自己定 %采样间隔 Ts=1/fs; %...
2019-09-07 15:22:52
885
原创 PyTorch创建自己的图像分类数据集
制作个人图像分类数据集过程中,查阅资料发现几个问题,在解决的过程中总结! 个人认为这是很通俗易懂详细的过程了 制作个人分类用数据集具体步骤如下: 1、将个人收集的图片归到一个文件夹内如下图: 我个人做实验收集的文件夹data内图如下图 2、建立与图片对应的label文件(标签文件),我label中数据如下: 3、在文件夹内建立py代码文件,本实验代码如下: import dataset im...
2019-09-05 09:31:02
13187
11
原创 可视化学深度学习
深度学习是什么? 这不用说了,怎么开始学深度学习呢?当然会给你推荐深度学习的经典之作《“花书”》,当你满心欢悦拿到这本书的时候(大佬除外),你会发现HOW What is are? 然后就开始了入门到放弃。。。。。 因此为提高大家入坑深度学习的兴趣,我简单讲一下自己的学习过程(本人菜鸟) 首先我是学习下深度学习是啥,我总结最重要的是理解卷积神经网络!何为卷积神经网络呢,即为具有深度的神经网络,说起...
2019-09-04 16:17:05
710
1
原创 cifar10 数据集转换成图片中imsave()问题详解
`做cifar10 数据集转换成图片程序,发现许多现有程序运行过程中会有imsave()报错,目前许多解释未安装pillow库所致,本人经过验证发现,目前pip安装库后,仍会出现scipy.misc内没有imsave()的问题,研究发现根本问题在于scipy.misc已经将imsave(),imread弃用,现在代替使用为mageio库,有做cifar10 数据集转换成图片等过程的朋友可以尝试使用...
2019-09-04 15:49:02
647
1
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人
RSS订阅