clear?all;
close?all;

%----------------信號說明—————————————————%
%???normal880-1????正常齒輪箱??實測轉速880
%???normal880-2????加載0.2A下正常齒輪箱?數據記錄?實測轉速800
%???normal880-3????加載0.1A下正常齒輪箱據記錄?實測轉速820
%???normal880-4????加載0.05A下正常齒輪箱據記錄?實測轉速852?
%???傳感器:
%???1------------TACH1??轉速，光電，?每周1個脈沖
%???2------------ch1??輸入軸?X方向位移????
%???3------------ch2???Y方向位移
%???4------------CH3???加速度?輸入軸電機側軸承Y??
%???5------------CH4??加速度?輸出軸電機側軸承Y
%???6------------CH5????加速度?輸入軸負載側軸承Y
%???7------------CH6?加速度?輸出軸負載?側軸承Y
%???8------------CH7??加速度?輸出軸負載?側軸承X
%???9------------CH8??磁電式速度??輸出軸負載?側軸承X?
%-------------------------------------------------------%
fs=2000*2.56;????????%?采樣率
load?data.mat
t=0:1/fs:（length（signal（1:））-1）/fs;?????%以時間為橫坐標

%%?去噪法二：小波分解消噪
%使用一維信號的自動消噪函數wden對信號進行消噪
%將信號使用小波函數‘sym5‘分解到第5層?
%使用minimaxi閾值選擇對系數進行處理，消除噪聲信號?
n=length（signal（1:））;
for?i=1:1:4
????signalnew（i:）=wden（signal（i:）‘minimaxi‘‘s‘‘mln‘5‘sym5‘）;???%小波分解去噪
end

%?時域波形圖---去噪前，去噪后
figure（1）
subplot（221）;
plot（tsignal（1:）‘-b‘）;hold?on;plot（tsignalnew（1:）‘-r‘）;
xlim（[011]）;ylabel（‘加速度‘）;xlabel（‘時間（s）‘）;title（‘1.正常齒輪箱（880r/min）‘）;
subplot（222）;
plot（tsignal（2:）‘-b‘）;hold?on;plot（tsignalnew（2:）‘-r‘）;
xlim（[011]）;ylabel（‘加速度‘）;xlabel（‘時間（s）‘）;title（‘2.加載0.2A下正常齒輪箱（800r/min）‘）;
subplot（223）;
plot（tsignal（3:）‘-b‘）;hold?on;plot（tsignalnew（3:）‘-r‘）;
xlim（[011]）;xlabel（‘加速度‘）;xlabel（‘時間（s）‘）;title（‘3.加載0.1A下正常齒輪箱（820r/min）‘）;
subplot（224）;
plot（tsignal（4:）‘-b‘）;hold?on;plot（tsignalnew（4:）‘-r‘）;
xlim（[011]）;ylabel（‘加速度‘）;xlabel（‘時間（s）‘）;title（‘4.加載0.05A下正常齒輪箱（852r/min）‘）;
legend（‘去噪前‘‘去噪后‘）;

figure;
plot（abs（fft（signal（1:）））1024）;

%時域特征計算----結果顯示在命令框
for?i=1:1:4
????feature（i1）=var（signalnew（i:））;????????%特征1?方差
????feature（i2）=skewness（signalnew（i:））;???%特征2?偏度
????feature（i3）=kurtosis（signalnew（i:））;???%特征3?峰度
????feature（i4）=std（signalnew（i:））;????????%特征4?均方根
????disp（[‘信號‘num2str（i）‘的方差、偏度、峰度、均方根分別為‘]）;???%輸出
????disp（feature（i:））;
end

%小波變換
percent=zeros（433）;
for?i=1:1:4
????percent（i1）=i;
????c=cwt（signalnew（i:）1:32‘morl‘）;???%morlet小波??32維
????a=zeros（321）;
????for?j=1:1:32
????????for?jj=1:1:n
????????????a（j1）=a（j1）+（c（jjj））.^2;
????????end?
????end????
????%計算每個尺度對應能量占總能量的百分比
????sum1=0;
????for?j=1:1:32
????????sum1=sum1+a（j）;
????end
????b=zeros（321）;
????for?j=1:1:32
????????b（j1）=a（j1）./sum1;
????end
????b=b‘;
????percent（i2:1:33）=b（1:）;
????figure（2）;
????color=[0?0?0;
???????????0?0?1;
???????????0?1?0;
???????????1?0?0];
????for?i=1:1:4
????????plot（percent（i2:33）‘color‘color（i:））;hold?on;grid?on;
????end
????xlabel（‘尺度‘）;ylabel（‘各尺度能量百分比‘）;title（‘齒輪振動信號能量-尺度圖‘）;
end
legend（‘信號1‘‘信號2‘‘信號3‘‘信號4‘）

%功率譜計算?Welch法
nfft=128;
window=hamming（100）;
noverlap=20;
figure（3）;
for?i=1:1:4
????[Pxxf]=pwelch（signalnew（i:）windownoverlapnfftfs）;
????Pxx_log=10*log（Pxx）;
????plot（fPxx_log‘color‘color（i:））;hold?on;grid?on;
e

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????486691??2017-08-29?16:17??2017082902\data.mat
?????文件????????3599??2017-08-29?16:48??2017082902\process.m
?????文件????????8239??2017-08-30?14:35??2017082902\process_gui.fig
?????文件???????14743??2017-08-30?14:36??2017082902\process_gui.m
?????目錄???????????0??2017-08-30?14:34??2017082902\