【正文】 。數(shù)據(jù)\測試數(shù)據(jù)\good\后軸承故障\1號機\39。 %文件夾2路徑 39。數(shù)據(jù)\測試數(shù)據(jù)\good\前軸承故障\1號機\39。 %文件夾3路徑 39。數(shù)據(jù)\測試數(shù)據(jù)\good\異物\1號\39。 %文件夾4路徑 % 39。數(shù)據(jù)\測試數(shù)據(jù)\good\異物和軸承故障39。 %文件夾5路徑 }。 Title = {39。1號好電機39。 39。1號電機后軸承故障39。 39。1號電機前軸承故障39。 39。1號電機異物39。 % 39。1號電機異物和軸承39。 }。FileNum = [50 25 25 25]。 %文件夾中文件數(shù)目Fs = 50000。 %采樣頻率X = cell(4,1)。 %存儲波形數(shù)據(jù)NL = cell(4,1)。 %文件波形點數(shù) mag = cell(4,1)。 %fft變換數(shù)據(jù) C = zeros(4,50)。 haxes = zeros(1,4)。%% 讀取數(shù)據(jù)for i = 1 : 4 for j = 1 : FileNum(i) FilePath = [Folder{i} num2str(j)]。 %文件路徑 fid = fopen(FilePath, 39。r39。)。 %打開文件 X{i}{j} = fread(fid, 39。float3239。)。 %將數(shù)據(jù)存入cell數(shù)組 fclose(fid)。 %關(guān)閉文件 endend%% 繪制波形%{figure(39。name39。,39。波形39。)。 hold on 。 for j = 1 : FileNum(1) N = length(X{1}{j})。 t = 0 : (N1)。 t = t/Fs。 plot(t, X{1}{j} )。 %繪制圖形 % pause end %} %% FFT,繪制頻譜圖 for i = 1 : 4 forj=1 : FileNum(i) % fft變換 NL{i}{j}=length(X{i}{j})。 % N=NL{i}{j}。 % n=0:N1。 % f=n*Fs/N。 % Y=fft(X{i}{j},N)。 % mag{i}{j}=abs(Y)。 % end end %% 繪制頻譜圖 for i = 2 : 4 FigureName = [Title{1} 39。和39。 Title{i}]。 % figure(39。name39。, FigureName)。 % axes1=subplot(211)。 % box on。 % title(Title{1})。 % xlabel(39。頻率(Hz)39。)。 % ylabel(39。振幅39。)。 % hold on % MAG = zeros(40000,1)。 for j = 1: FileNum(1) % %N = NL{1}{j}。 MAG = MAG+mag{1}{j}(1:40000)。 % %plot( f(1200:1400), mag{1}{j}(1200:1400),39。color39。, rand(1,3) )。 %plot( f(1:N/2), mag{1}{j}(1:N/2),39。color39。, rand(1,3) )。 end plot( f(1:40000), MAG/FileNum(1),39。color39。, rand(1,3) )。 axes2 = subplot(212)。 box on。 title(Title{i})。 xlabel(39。頻率(Hz)39。)。 ylabel(39。振幅39。)。 hold on MAG = zeros(40000,1)。 for j = 1 : FileNum(i) %N = NL{i}{j}。 MAG = MAG+mag{i}{j}(1:40000)。 % plot( f(1200:1400), mag{i}{j}(1200:1400),39。color39。, rand(1,3) )。 % plot( f(1:N/2), mag{i}{j}(1:N/2),39。color39。, rand(1,3) )。 end plot( f(1:40000), MAG/FileNum(i),39。color39。, rand(1,3) )。 linkaxes([axes2 axes1],39。x39。)。 end %% 特征段 最大值提取 for i = 1 : 4 for j = 1 : FileNum(i) C(i,j) = max( mag{i}{j}(1200:1400) )。 end end figure(39。name39。,39。箱型圖39。)。 boxplot(C(1:4, 1:25)39。,39。label39。,{39。1號好電機39。,39。后軸承故障39。,39。前軸承故障39。,39。異物39。})。 box on。for i = 1 : 4 M(i) = mean(C(i,1:25))。 S(i) = std(C(i,1:25))。 %標準差 V(i) = var(C(i,1:25))。 %方差 SK(i) = skewness(C(i,1:25))。 %偏度 K(i) = kurtosis(C(i,1:25))。 %峰度 endfigure(39。name39。,39。平均值39。)。box on。hbm = bar(M,39。b39。)。title(39。幾種電機狀況特征段最大值的平均值39。)。set(gca,39。xticklabel39。,{39。1號好電機39。,39。后軸承故障39。,39。前軸承故障39。,39。異物39。})。figure(39。name39。,39。標準差39。)。box on。hbs = bar(S,39。m39。)。title(39。幾種電機狀況特征段最大值的標準差39。)。set(gca,39。xticklabel39。,{39。1號好電機39。,39。后軸承故障39。,39。前軸承故障39。,39。異物39。})。figure(39。name39。,39。偏度39。)。box on。hbsk = bar(SK,39。g39。)。title(39。幾種電機狀況特征段最大值的偏度39。)。set(gca,39。xticklabel39。,{39。1號好電機39。,39。后軸承故障39。,39。前軸承故障39。,39。異物39。})。figure(39。name39。,39。峰度39。)。box on。hbk = bar(K,39。k39。)。title(39。幾種電機狀況特征段最大值的峰度39。)。set(gca,39。xticklabel39。,{39。1號好電機39。,39。后軸承故障39。,39。前軸承故障39。,39。異物39。})