故障诊断Matlab常用技巧

1 巴特沃斯低通滤波+傅立叶变换

2 离散傅立叶变换

3 小波变换

4 短时傅立叶等等

具体做法见下处代码

clc
close all
%%for example
fs = 12800;
t = (0 :52121)/fs;
data = sin(t).*cos(t);
figure
FFT(data,fs)
figure
BFFT(data,fs)
figure
XFFT(data,fs)
figure
xiaobo(data(30000:40000),fs)wlen = 512;
hop = 2;figure
shortFouier(wlen,hop,data(30000:40000),fs);
%巴特沃斯低通滤波+傅里叶变换function BFFT(data,fs)datay = data(30000:40000);  %取中间稳定数据lengthNormal = length(datay);datay = datay - mean(datay);wc = 300/(fs/2);        %低道滤波截止频率为300Hz[a,b] = butter(4,wc);datay = filter(a,b,datay);rifft = fft(datay);fz = (0:lengthNormal-1)*fs/lengthNormal;P2 = 2*abs(rifft/lengthNormal);plot(fz(1:lengthNormal/2),P2(1:lengthNormal/2),'r');ax = gca;ax.XLim = [0 300];xlabel('频率Hz');ylabel('幅值');title('频谱图（频域）');
end%小波去噪
function data = loseNoise(data)lev = 15;data = wden(data,'heursure','s','one',lev,'db4');
end%傅里叶
function FFT(data,fs)datay = data(30000:40000);  %取中间稳定数据lengthD = length(datay);rifft = fft(datay);fz = (0:lengthD - 1)*fs/lengthD;P2 = 2*abs(rifft)/lengthD;plot(fz(1:lengthD/2),P2(1:lengthD/2),'m');ax = gca;ax.XLim = [0 300];xlabel('频率(Hz)');ylabel('幅值');title('频谱图');
end%小波去噪+傅里叶
function XFFT(data,fs)datay = loseNoise(data);FFT(datay,fs);
end
%短时傅里叶function shortFouier(wlen,hop,data,fs)lengthData = length(data);Normal1 = wextend(1,'sym',data,round(lengthData/2));% 镜像延拓Normal2 = wkeep1(Normal1,lengthData + 1*wlen);% 短时 Fouierh = hamming(wlen);% 海明窗的窗长f = 0:1:200;[tfr2,f,t2] = spectrogram(Normal2,h,wlen-hop,f,fs);tfr2 = tfr2 * 2/wlen*2;imagesc(t2-wlen/fs/2,f,abs(tfr2));%[S,F,T,P] = spectrogram(Normal2,h,wlen-hop,f,frequency);%surf(T,F,10*log10(P),'edgecolor','none');%view(0,1000);ylim([-10 200])xlabel('Time, s');ylabel('Frequency, Hz');colorbar
end%小波去噪+小波变换
function xiaobo(data,fs)data = loseNoise(data);t = (0:length(data)-1)/fs;plot(t,data);wavename = 'cmor3-3';totalscal = 1024;Fc = centfrq(wavename);c = 2*Fc*totalscal;scals = c./(1:totalscal);f = 8*scal2frq(scals,wavename,1/fs);coefs = cwt(data,scals,wavename);figureimagesc(t,f,abs(coefs));set(gca,'YDir','normal');colorbar;xlabel('时间 t/s');ylabel('频率 f/Hz');title('小波时频图');
end

故障诊断Matlab常用技巧相关推荐

Matlab常用技巧
1.常用标点的功能逗号:用作要显示计算结果的指令与其后面的指令之间的分隔符:用作输入量与输出量之间的分隔符:用作数组元素分隔符. 分号:不显示计算结果的分隔符:数组行间的分隔符. 冒号:生成一维数值 ...
MATLAB常用技巧之复数判断，矩阵中空值判断，空矩阵判断
判断矩阵或向量中的空值(NAN) returnValue = isnan(matrix); 其中,matrix表示需要判断矩阵(向量),returnValue 是返回值,矩阵matrix中是空值(na ...
matlab快速入门案例及常用技巧 | 《matlab数学建模方法与实践（第三版）》学习笔记
目录快速入门案例: 解决流程: 具体实现: 一.获取数据二.数据探索和建模三.分享结果常用技巧一.常用标点功能二.常用操作指令三.指令编辑操作键四.matlab数据类型五.开发模式 ...
MATLAB小技巧（9）图片合成视频与视频分帧
MATLAB小技巧(9)图片合成视频与视频分帧前言一. MATLAB仿真-序列合成二. MATLAB仿真-视频分帧三. 小结前言 MATLAB进行图像处理相关的学习是非常友好的,可以从零开始 ...
MATLAB小技巧（25）竞争神经网络与SOM神经网络
MATLAB小技巧(25)竞争神经网络与SOM神经网络前言一. MATLAB仿真二. 仿真结果三. 小结前言 MATLAB进行图像处理相关的学习是非常友好的,可以从零开始,对基础的图像处理都 ...
MATLAB小技巧（32）FCM聚类
MATLAB小技巧(32)FCM聚类前言一. MATLAB仿真二. 仿真结果三. 小结前言 MATLAB进行图像处理相关的学习是非常友好的,可以从零开始,对基础的图像处理都已经有了封装好的许 ...
MATLAB小技巧（20）矩阵分析--主成分回归
MATLAB小技巧(20)矩阵分析--主成分回归前言一. MATLAB仿真二. 仿真结果三. 小结前言 MATLAB进行图像处理相关的学习是非常友好的,可以从零开始,对基础的图像处理都已经有 ...
Matlab常用操作入门及电力电子系统仿真
引言刚刚比较系统地学了一遍Matlab(主要是基础知识和电力电子方面的使用),专门撰写此文,写一下自己对于matlab及Power electronic方面的应用的总结和心得. Matlab入门因 ...
MATLAB小技巧（1）绘制图像能量图
MATLAB小技巧(1)绘制图像能量图前言一. MATLAB仿真二. 仿真结果三. 小结前言 MATLAB进行图像处理相关的学习是非常友好的,可以从零开始,对基础的图像处理都已经有了封装好的 ...

故障诊断Matlab常用技巧

故障诊断Matlab常用技巧

1 巴特沃斯低通滤波+傅立叶变换

2 离散傅立叶变换

3 小波变换

4 短时傅立叶等等

故障诊断Matlab常用技巧相关推荐

最新文章

热门文章