轴承故障特征计算MATLAB代码

---

前言

写的一些计算轴承故障特征频率的MATLAB代码，方便计算。

---

一、轴承特征频率计算公式

轴承特征频率计算公式

轴承示意

二、MATLAB代码

代码如下（示例）：

clc;
clear;
close all;
%% XJTU-SY轴承数据集为例
fs = 25600; %采样频率
N = 40;     %转频 Hz
n = 8;      %滚动体数量
d = 7.92;   %滚动体直径
D = 34.55;  %轴承中径   
alfe = 0;   %滚动体接触角


fprintf('\n转频：%5.4f Hz ',N); %转频

FTF = 1/2*N*(1-d*cos(alfe)/D);  
fprintf('\nFTF保持架故障：%5.4f Hz',FTF); 
fprintf('\nFTF保持架故障周期点数：%5.4f ',fs/FTF);

BSF = 1/2*N*D/d*(1-(d*cos(alfe)/D)^2);
fprintf('\nBSF滚动体旋转故障：%5.4f Hz',BSF);
fprintf('\nBSF滚动体旋转故障周期点数：%5.4f',fs/BSF);

BPFO = 1/2*N*n*(1-d*cos(alfe)/D);
fprintf('\nBPFO外圈故障：%5.4f Hz',BPFO);
fprintf('\nBPFO外圈故障周期点数：%5.4f ',fs/BPFO);

BPFI = 1/2*N*n*(1+d*cos(alfe)/D);
fprintf('\nBPFI内圈故障：%5.4f Hz',BPFI);
fprintf('\nBPFI内圈故障周期点数：%5.4f ',fs/BPFI);

---

提示

若有错误，请批评指正！