滚动轴承故障诊断

打开实时脚本

这个例子展示了如何基于加速度信号进行滚动轴承的故障诊断，特别是在存在来自其他机器部件的强屏蔽信号的情况下。该实例将演示如何应用包络谱分析和谱峰度来诊断轴承故障，并能够扩展到大数据应用。

问题概述

滚动元件轴承中的局部故障可能发生在外圈、内圈、保持架或滚动元件中。当滚动元件撞击外圈或内圈上的局部故障，或滚动元件上的故障撞击外圈或内圈[1]时，轴承和响应传感器之间会产生高频共振。下面的图片显示了一个滚动的元素撞击内部的一个局部断层。问题是如何检测和识别各种类型的故障。

机械故障预防技术(MFPT)挑战数据

MFPT挑战数据[4]包含从不同故障条件下的机器收集的23个数据集。前20组数据采集自轴承试验台，其中3组处于良好状态，3组处于恒载下的外圈故障，7组处于各种负载下的外圈故障，7组处于各种负载下的内圈故障。剩下的3个数据集来自真实的机器:一个油泵轴承，一个中速轴承和一个行星轴承。故障位置未知。在本例中，仅使用从已知条件的试验台收集的数据。

每个数据集包含一个加速度信号“gs”，采样率“sr”，轴转速“rate”，负载重量“load”，以及四个代表不同故障位置的关键频率:球通频率外场(BPFO)，球通频率内场(BPFI)，基本列车频率(FTF)和球自旋频率(BSF)。这是临界频率[1]的公式。

外场球通频率(BPFO)

$B P F O ＝ \frac{n f_{r}}{2} （ 1 - \frac{d}{D} 因为 ϕ ）$

球通频率，内圈(BPFI)

$B P F 我＝ \frac{n f_{r}}{2} （ 1 + \frac{d}{D} 因为 ϕ ）$

列车基本频率(FTF)，也称为笼速

$F T F ＝ \frac{f_{r}}{2} （ 1 - \frac{d}{D} 因为 ϕ ）$

球(辊)自旋频率

$B 年代 F ＝ \frac{D}{2 d} ［ 1 - {（ \frac{d}{D} 因为 ϕ ）}^{2} ］$

如图所示， $d$ 是球的直径， $D$ 为节径。的变量 $f_{r}$ 为轴转速， $n$ 是滚动元素的个数， $ϕ$ 为轴承接触角[1]。

轴承诊断的包络谱分析

在MFPT数据集中，轴转速是恒定的，因此不需要执行顺序跟踪作为预处理步骤来去除轴转速变化的影响。

当滚动单元撞击外圈或内圈局部故障时，或滚动单元上的故障撞击外圈或内圈时，冲击会调制相应的临界频率，如BPFO、BPFI、FTF、BSF。因此，由幅度解调产生的包络信号传递了更多的诊断信息，这是无法从原始信号的频谱分析中获得的。以MFPT数据集中的内部种族故障信号为例。

dataInner = load(fullfile(matlabroot)，“工具箱”，“predmaint”，.．.“predmaintdemos”，“bearingFaultDiagnosis”，.．.“train_data”，“InnerRaceFault_vload_1.mat”));

将原始的内部赛车故障数据在时域上可视化。

xInner = dataInner.bearing.gs;f罪人= dataInner.bearing.sr;tInner = (0:length(xInner)-1)/ f罪人;图(tInner, xInner) xlabel(“时间,(s)”) ylabel (“加速度(g)”)标题(“原始信号:内部种族问题”) xlim([0 0.1])

在频域中可视化原始数据。

figure [pInner, fpInner] = pspectrum(xInner, f罪人);pInner = 10*log10(pInner);plot(fpInner, pInner) xlabel(的频率(赫兹)) ylabel (“功率谱(dB)”)标题(“原始信号:内部种族问题”)传说(的功率谱）

现在放大低频范围内原始信号的功率谱，仔细观察BPFI处的频率响应及其前几个谐波。

figure plot(fpInner, pInner) ncomb = 10;helperPlotCombs(ncomb, dataInner.BPFI)的频率(赫兹)) ylabel (“功率谱(dB)”)标题(“原始信号:内部种族问题”)传说(的功率谱，“BPFI谐波”) xlim([0 1000])

在BPFI及其谐波中看不到清晰的模式。对原始信号的频率分析不能提供有用的诊断信息。

从时域数据可以看出，原始信号的幅值被调制到一定的频率，调制的主频在1/0.009 Hz左右 $\approx$ 111赫兹。已知滚动元件撞击内圈局部故障即BPFI的频率为118.875 Hz。这表明轴承可能有一个内圈故障。

图subplot(2,1,1) plot(tInner, xInner) xlim([0.04 0.06]) title(“原始信号:内部种族问题”) ylabel (“加速度(g)”)注释(“doublearrow”， [0.37 0.71]， [0.8 0.8]) text(0.047, 20， ['0.009秒\约1/BPFI, BPFI = 'num2str (dataInner.BPFI)))

为了提取调制振幅，计算原始信号的包络线，并将其可视化在底部子图上。

subplot(2,1,2) [pEnvInner, fEnvInner, xEnvInner, tEnvInner] = envspectrum(xInner, f罪人);plot(tEnvInner, xEnvInner) xlim([0.04 0.06]) xlabel(“时间(s)”) ylabel (“加速度(g)”)标题(“包络信号”）

现在计算包络信号的功率谱，看看在BPFI处的频率响应及其谐波。

figure plot(fEnvInner, pEnvInner) xlim([0 1000]) ncomb = 10;helperPlotCombs(ncomb, dataInner.BPFI)的频率(赫兹)) ylabel (“峰值振幅”)标题(“包络谱:内部种族错误”)传说(“包络谱”，“BPFI谐波”）

结果表明，大部分能量集中在BPFI及其谐波上。这表明轴承的内圈故障，这与数据的故障类型相匹配。

包络谱分析在其他故障类型中的应用

现在对正常数据和外部赛车故障数据重复同样的包络谱分析。

datannormal = load(fullfile(matlabroot)，“工具箱”，“predmaint”，.．.“predmaintdemos”，“bearingFaultDiagnosis”，.．.“train_data”，“baseline_1.mat”));xNormal = datnormal .bearing.gs;fsNormal = datnormal .bearing.sr;tNormal = (0:length(xNormal)-1)/fsNormal;[pEnvNormal, fEnvNormal] = envspectrum(xNormal, fsNormal);图(fEnvNormal, pEnvNormal) ncomb = 10;helperPlotCombs (ncomb [dataNormal。BPFOdataNormal.BPFI]) xlim([0 1000]) xlabel(的频率(赫兹)) ylabel (“峰值振幅”)标题(“包络谱:正常”)传说(“包络谱”，“BPFO谐波”，“BPFI谐波”）

正如预期的那样，正常轴承信号的包络谱在BPFO或BPFI处不显示任何显著的峰值。

dataOuter = load(fullfile(matlabroot)，“工具箱”，“predmaint”，.．.“predmaintdemos”，“bearingFaultDiagnosis”，.．.“train_data”，“OuterRaceFault_2.mat”));xOuter = dataOuter.bearing.gs;fsOuter = dataOuter.bearing.sr;tOuter = (0:length(xOuter)-1)/fsOuter;[pEnvOuter, fEnvOuter, xEnvOuter, tEnvOuter] = envspectrum(xOuter, fsOuter);figure plot(fEnvOuter, pEnvOuter) ncomb = 10;helperPlotCombs(ncomb, dataOuter.BPFO) xlim([0 1000]) xlabel(的频率(赫兹)) ylabel (“峰值振幅”)标题(“包络谱:外部种族故障”)传说(“包络谱”，“BPFO谐波”）

对于外部故障信号，在BPFO谐波处也没有明显的峰值。包络谱分析是否不能区分外圈故障轴承和健康轴承?让我们退一步看一下在不同条件下的时域信号。

首先让我们在时域中再次可视化信号并计算它们的峰度。峰度是随机变量的第四个标准化矩。它的特征是信号的冲动性或随机变量的尾部的重量。

图subplot(3,1,1) kurtInner = kurtosis(xInner);plot(tInner, xInner)“加速度(g)”)标题(“内部种族断层，峰度=”num2str(kurtInner)]) xlim([0 0.1]) subplot(3,1,2) kurtNormal =峰度(xNormal);plot(tNormal, xNormal)“加速度(g)”)标题(“正常，峰度=”num2str(kurtNormal)]) xlim([0 0.1]) subplot(3,1,3) kurtOuter = kurtosis(xOuter);plot(tOuter, xOuter) xlabel(“时间(s)”) ylabel (“加速度(g)”)标题(“外部种族断层，峰度=”num2str(kurtOuter))) xlim([0 0.1])

结果表明，内部故障信号具有较大的冲动性，使得包络谱分析能够有效地捕捉到BPFI处的故障特征。对于外部赛车故障信号，BPFO处的幅值调制较为明显，但被强噪声所掩盖。正常信号不显示任何幅度调制。在包络谱分析之前，提取在BPFO(或提高信噪比)处振幅调制的脉冲信号是一个关键的预处理步骤。下一节将介绍峰度图和谱峰度，提取峰度最高的信号，并对滤波后的信号进行包络谱分析。

波段选择的峰度图和谱峰度

峰度图和谱峰度计算在频带内的局部峰度。它们是定位具有最高峰度(或最高信噪比)[2]的频段的强大工具。在确定峰度最高的频带后，可以对原始信号进行带通滤波器，以获得更脉冲的信号进行包络谱分析。

等级= 9;图形峭度图(xOuter, fsOuter, level)

峰度图表明，以2.67 kHz为中心，带宽为0.763 kHz的频段峰度最高，为2.71。

现在使用由峰度图所建议的最佳窗口长度来计算谱峰度。

图wc = 128;pkurtosis(xOuter, fsOuter, wc)

要在光谱图上可视化频带，计算光谱图并将光谱峰度放置在一侧。用另一种方式解释谱峰度，高谱峰度值表示相应频率下功率的高方差，这使得谱峰度成为定位信号[3]非平稳分量的有用工具。

help perspectrogramandspectralkurtosis (xOuter, fsOuter, level)

采用所建议的中心频率和带宽对信号进行带通滤波，可以增强信号的峰度，提取出外部故障的调制幅值。

[~， ~， ~， fc， ~， BW] = kurtogram(xOuter, fsOuter, level);BPF = designfilt(“bandpassfir”，“FilterOrder”, 200,“CutoffFrequency1”fc-BW / 2,.．.“CutoffFrequency2”fc + BW / 2,“SampleRate”, fsOuter);xOuterBpf = filter(bpf, xOuter);[pEnvOuterBpf, fEnvOuterBpf, xEnvOuterBpf, tEnvBpfOuter] = envspectrum(xOuter, fsOuter，.．.“FilterOrder”, 200,“乐队”， [fc-BW/2 fc+BW/2]);figure subplot(2,1,1) plot(tOuter, xOuter, tEnvOuter, xEnvOuter)“加速度(g)”)标题(“原始信号:外部种族故障，峰度=”， num2str(kurtOuter)]) xlim([0 0.1]) legend(“信号”，“信封”kurtOuterBpf =峰度(xOuterBpf);plot(tOuter, xOuterBpf, tEnvBpfOuter, xEnvOuterBpf)“加速度(g)”xlabel([0 0.1])“时间(s)”)标题(“带通滤波信号:外部种族故障，峰度=”， num2str(kurtOuterBpf)])“信号”，“信封”）

可以看出，带通滤波后峰度值增大。现在想象一下频域的包络信号。

图(fEnvOuterBpf, pEnvOuterBpf);Ncomb = 10;helperPlotCombs(ncomb, dataOuter.BPFO) xlim([0 1000]) xlabel(的频率(赫兹)) ylabel (“峰值振幅”)标题(带通滤波信号的包络谱:外部种族故障)传说(“包络谱”，“BPFO谐波”）

结果表明，通过将原始信号带通滤波到由峰态图和谱峰度指示的频带，包络谱分析能够揭示出故障在BPFO处的特征及其谐波。

批处理

现在，让我们使用文件集成数据存储将该算法应用于一批训练数据。

工具箱中提供了数据集的有限部分。将数据集复制到当前文件夹，并启用写权限:

拷贝文件(.．.fullfile (matlabroot“工具箱”，“predmaint”，“predmaintdemos”，.．.“bearingFaultDiagnosis”)，.．.“RollingElementBearingFaultDiagnosis-Data-master”) fileattrib (fullfile (“RollingElementBearingFaultDiagnosis-Data-master”，“train_data”，‘* .mat‘)，' + w ') fileattrib (fullfile (“RollingElementBearingFaultDiagnosis-Data-master”，“test_data”，‘* .mat‘)，' + w '）

要查看完整的数据集，请转到此链接https://github.com/mathworks/RollingElementBearingFaultDiagnosis-Data下载整个存储库为zip文件，并将其保存到与实时脚本相同的目录中。使用以下命令解压文件:

如果存在(“RollingElementBearingFaultDiagnosis-Data-master.zip”，“文件”)解压缩(“RollingElementBearingFaultDiagnosis-Data-master.zip”）结束

本例中的结果是从完整的数据集生成的。完整数据集包含一个包含14个mat文件的训练数据集(2个正常文件，4个内部比赛故障，7个外部比赛故障)和一个包含6个mat文件的测试数据集(1个正常文件，2个内部比赛故障，3个外部比赛故障)。

通过将函数句柄分配给ReadFcn而且WriteToMemberFcn，文件集成数据存储将能够导航到文件以检索所需格式的数据。例如，MFPT数据有一个结构轴承它储存振动信号gs，采样率老，等等。而不是返回承载结构本身readMFPTBearing函数，使文件集成数据存储返回振动信号gs在内部轴承数据结构。

fileLocation = fullfile(“。”，“RollingElementBearingFaultDiagnosis-Data-master”，“train_data”）;fileExtension =“.mat”；集合训练= fileEnsembleDatastore(fileLocation, fileExtension);ensembleTrain。ReadFcn＝@readMFPTBearing; ensembleTrain.DataVariables = [“gs”，“老”，“速度”，“负载”，“BPFO”，“BPFI”，““保障未来粮食供给””，“巴斯夫”];ensembleTrain。条件变量= [“标签”，“文件名”];ensembleTrain。WriteToMemberFcn = @writeMFPTBearing;ensembleTrain。年代electedVariables = [“gs”，“老”，“速度”，“负载”，“BPFO”，“BPFI”，““保障未来粮食供给””，“巴斯夫”，“标签”，“文件名”］

ReadFcn: @readMFPTBearing WriteToMemberFcn: @writeMFPTBearing DataVariables: [8×1 string] IndependentVariables: [0×0 string] ConditionVariables: [2×1 string] SelectedVariables: [10×1 string] ReadSize: 1 NumMembers: 14 LastMemberRead: [0×0 string] Files: [14×1 string]

integrbletraintable = tall(集合火车)

使用“本地”配置文件启动并行池(parpool)…连接6个工人。ensembleTrainTable = M×10高表gs sr速率加载BPFO BPFI“保障未来粮食供给”BSF标签文件名  _________________ _____ ____ ____ ______ ______ ______ _____ __________________ ________________________ [ 25 146484×1双)48828 0 81.125 118.88 14.838 63.91“内座圈的错”“InnerRaceFault_vload_1”(146484×1双)48828 25 50 81.125 118.88 14.838 63.91“内座圈的错”“InnerRaceFault_vload_2”(146484×1双)48828 25 100 81.125 118.88 14.838 63.91“内座圈的错”“InnerRaceFault_vload_3”[146484×1 double] 48828 25 150 81.125 118.88 14.838 63.91 "InnerRaceFault " "InnerRaceFault_vload_4" [146484×1 double] 48828 25 200 81.125 118.88 14.838 63.91 "InnerRaceFault " "InnerRaceFault_vload_5" [585936×1 double] 97656 25 270 81.125 118.88 14.838 63.91 "OuterRaceFault " "OuterRaceFault_1" [585936×1 double] 97656 25 270 81.125 118.88 14.838 63.91 "OuterRaceFault " "OuterRaceFault_2" [146484×1 double] 48828 25 25 81.125 118.88 14.838 63.91 "OuterRaceFault " "OuterRaceFault_vload_1":: : : : : : : : : : : : : : : : : : :

从最后一节的分析中，注意到带通滤波包络谱振幅在BPFO和BPFI是轴承故障诊断的两个条件指标。因此，下一步是从所有的训练数据中提取这两个条件指标。为了使算法更健壮，可以设置一个较窄的带宽(bandwidth = $10 Δ f$ ,在那里 $Δ f$ 为功率谱的频率分辨率)在BPFO和BPFI附近，然后求出该窄带内的最大振幅。算法包含在bearingFeatureExtraction下面列出的函数。注意，在本例的其余部分中，围绕BPFI和BPFO的包络谱振幅称为“BPFIAmplitude”和“BPFOAmplitude”。

要并行处理数据，请使用以下代码% ppool = gcp;% n = numpartitions(ensemble bletrain, ppool);% parfor ct = 1:n% subEnsembleTrain =分区(ensembleTrain, n, ct);%重置(subEnsembleTrain);% while hasdata(subEnsembleTrain)% bearingFeatureExtraction (subEnsembleTrain);%结束%结束ensembleTrain。DataVariables = [ensembleTrain.DataVariables;“BPFIAmplitude”；“BPFOAmplitude”];重置(ensembleTrain)而hasdata (ensembleTrain) bearingFeatureExtraction (ensembleTrain)结束

在将新的条件指示器添加到文件集成数据存储中时，指定SelectedVariables从文件集成数据存储中读取相关数据，并创建包含提取的条件指示器的特征表。

ensembleTrain。年代electedVariables = [“BPFIAmplitude”，“BPFOAmplitude”，“标签”];featureTableTrain = tall(集合火车);featureTableTrain = gather(featureTableTrain);

使用并行池'local'计算高表达式:-通过1:0%完成计算0%完成

3秒内完成评估

featureTableTrain

featureTableTrain =14×3表BPFIAmplitude BPFOAmplitude Label _____________ _____________ __________________ 0.33918 0.082296“内部种族故障”0.31488 0.026599“内部种族故障”0.52356 0.036609“内部种族故障”0.52899 0.028381“内部种族故障”0.13515 0.012337“内部种族故障”0.004024 0.03574“外部种族故障”0.0044918 0.1835“外部种族故障”0.0074993 0.30166“外部种族故障”0.013662 0.12468“外部种族故障”0.0070963 0.28241“外部种族故障”0.011244 0.17975“外部种族故障”0.0036798 0.0050208“正常”0.00359 0.0069449“正常”

可视化已经创建的特性表。

图gscatter (featureTableTrain。BPFIAmplitude featureTableTrain。BPFOAmplitude featureTableTrain。标签、[]“牛+”)包含(“BPFI振幅”) ylabel (“BPFO振幅”）

BPFI振幅和BPFO振幅的相对值可能是不同故障类型的有效指标。这里创建了一个新特征，它是两个现有特征的对数比，并在一个按不同故障类型分组的直方图中可视化。

featureTableTrain。我OLogRatio = log(featureTableTrain.BPFIAmplitude./featureTableTrain.BPFOAmplitude); figure hold在直方图(featureTableTrain.IOLogRatio (featureTableTrain.Label = =“内部种族错误”)，“BinWidth”， 0.5)直方图(featureTableTrain.IOLogRatio(featureTableTrain.Label==“外部种族断层”)，“BinWidth”， 0.5)直方图(featureTableTrain.IOLogRatio(featureTableTrain.Label==“正常”)，“BinWidth”情节,0.5)(-1.5[-1.5南0.5 - 0.5],[0 3南0 3],“k——”)举行从ylabel (“数”)包含(“日志(BPFIAmplitude / BPFOAmplitude)”)传说(“内部种族错误”，“外部种族错误”，“正常”，分类边界的）

直方图显示了三种不同承载条件之间的清晰分离。BPFI和BPFO振幅之间的对数比是对轴承故障进行分类的有效特征。为了简化示例，导出了一个非常简单的分类器:if $日志（ \frac{BPFIAmp litude}{BPFOAmp litude} ） \leq - 1 ． 5$ ，轴承有外圈故障;如果 $- 1 ． 5 < 日志（ \frac{BPFIAmp litude}{BPFOAmp litude} ） \leq 0 ． 5$ ，轴承正常;如果 $日志（ \frac{BPFIAmp litude}{BPFOAmp litude} ） > 0 ． 5$ ，轴承有内圈故障。

使用测试数据集进行验证

现在，让我们将工作流应用到一个测试数据集，并验证上一节中获得的分类器。这里的测试数据包含1个正常数据集，2个内部赛车故障数据集，3个外部赛车故障数据集。

fileLocation = fullfile(“。”，“RollingElementBearingFaultDiagnosis-Data-master”，“test_data”）;fileExtension =“.mat”；= fileEnsembleDatastore(fileLocation, fileExtension);ensembleTest。ReadFcn＝@readMFPTBearing; ensembleTest.DataVariables = [“gs”，“老”，“速度”，“负载”，“BPFO”，“BPFI”，““保障未来粮食供给””，“巴斯夫”];ensembleTest。条件变量= [“标签”，“文件名”];ensembleTest。WriteToMemberFcn = @writeMFPTBearing;ensembleTest。年代electedVariables = [“gs”，“老”，“速度”，“负载”，“BPFO”，“BPFI”，““保障未来粮食供给””，“巴斯夫”，“标签”，“文件名”］

ReadFcn: @readMFPTBearing WriteToMemberFcn: @writeMFPTBearing DataVariables: [8×1 string] IndependentVariables: [0×0 string] ConditionVariables: [2×1 string] SelectedVariables: [10×1 string] ReadSize: 1 NumMembers: 6 LastMemberRead: [0×0 string] Files: [6×1 string]

ensembleTest。DataVariables = [ensembleTest.DataVariables;“BPFIAmplitude”；“BPFOAmplitude”];重置(ensembleTest)而hasdata (ensembleTest) bearingFeatureExtraction (ensembleTest)结束

ensembleTest。年代electedVariables = [“BPFIAmplitude”，“BPFOAmplitude”，“标签”];featureTableTest = tall(ensemble bletest);featureTableTest = gather(featureTableTest);

使用并行池“local”计算高表达式:-通过1 / 1:在1秒内完成计算

featureTableTest。我OLogRatio = log(featureTableTest.BPFIAmplitude./featureTableTest.BPFOAmplitude); figure hold在直方图(featureTableTrain.IOLogRatio (featureTableTrain.Label = =“内部种族错误”)，“BinWidth”， 0.5)直方图(featureTableTrain.IOLogRatio(featureTableTrain.Label==“外部种族断层”)，“BinWidth”， 0.5)直方图(featureTableTrain.IOLogRatio(featureTableTrain.Label==“正常”)，“BinWidth”， 0.5)直方图(featureTableTest.IOLogRatio(featureTableTest.Label==“内部种族错误”)，“BinWidth”， 0.1)直方图(featureTableTest.IOLogRatio(featureTableTest.Label==“外部种族断层”)，“BinWidth”， 0.1)直方图(featureTableTest.IOLogRatio(featureTableTest.Label==“正常”)，“BinWidth”情节,0.1)(-1.5[-1.5南0.5 - 0.5],[0 3南0 3],“k——”)举行从ylabel (“数”)包含(“日志(BPFIAmplitude / BPFOAmplitude)”)传说(“内部种族问题-火车”，“列车外轨故障”，“正常-列车”，.．.“内部种族错误-测试”，“外部种族故障-测试”，“正常-测试”，.．.分类边界的）

测试数据集的BPFI和BPFO振幅的对数比与训练数据集的对数比分布一致。最后一节得到的朴素分类器在测试数据集上达到了完美的精度。

应该注意的是，单一的特征通常不足以得到一般化的分类器。更复杂的分类器可以通过将数据分成多个片段(创建更多的数据点)，提取多个与诊断相关的特征，根据其重要级别选择特征的子集，并使用统计和机器学习工具箱中的分类学习者应用程序训练各种分类器来获得。有关此工作流程的更多细节，请参考示例“使用Simulink生成故障数据”。万博1manbetx

总结

这个例子展示了如何使用峰度图、谱峰度和包络谱来识别滚动轴承中不同类型的故障。将该算法应用于磁盘上的一批数据集，结果表明，带通滤波包络谱在BPFI和BPFO处的振幅是轴承诊断的两个重要条件指标。

参考文献

兰德尔，罗伯特·B，杰罗姆·安东尼。滚动轴承诊断教程。机械系统与信号处理“，．2011年第2期第25卷第485-520页。

[2]安东尼，Jérôme。"瞬态故障检测用峭度图的快速计算"机械系统与信号处理“，．2007年第1期，第21卷，108-124页。

[3]安东尼，Jérôme。光谱峰度:描述非平稳信号的有用工具。机械系统与信号处理“，．2006年第2期第20卷第282-307页。

[4]贝赫弗，埃里克。用于测试诊断和预测算法的基于条件的维护故障数据库。2013.https://www.mfpt.org/fault-data-sets/

辅助函数

函数bearingFeatureExtraction(整体)从轴承数据中提取条件指标数据=读取(集成);X = data.gs{1};Fs = data.sr;%临界频率BPFO = data.BPFO;BPFI = data.BPFI;等级= 9;[~， ~， ~， fc， ~， BW] =峰值图(x, fs, level);带通滤波包络谱[pEnvpBpf, fEnvBpf] = envspectrum(x, fs，)“FilterOrder”, 200,“乐队”， [max([fc-BW/2 0]) min([fc+BW/2 0.999*fs/2])]);deltaf = fEnvBpf(2) - fEnvBpf(1);BPFIAmplitude = max(pEnvpBpf((fEnvBpf > (bpf -5*deltaf)) & (fEnvBpf < (BPFI+5*deltaf))));bpfo振幅= max(pEnvpBpf((fEnvBpf > (bpf -5*deltaf)) & (fEnvBpf < (bpf +5*deltaf))));集合，表(BPFIAmplitude, BPFOAmplitude, BPFIAmplitude)“VariableNames”, {“BPFIAmplitude”，“BPFOAmplitude”}));结束

另请参阅

fileEnsembleDatastore