%?fs=12000;
%?load?X156_DE_time;
%?xdata=X156_DE_time;
%?xdata=xdata（1:10240）;
%?N=10240;
%?
%?
%?xdata=（xdata-mean（xdata））/std（xdata1）;%mean平均值?std標準偏差???對原始信號進行歸一處理
%?
%?%時域波形
%?figure（1）;
%?subplot（211）;
%?plot（1:Nxdata）;
%?xlabel（‘時間?t/n‘）;
%?ylabel（‘電壓?V/v‘）;

%采樣頻率
fs=12000;
%軸承外環故障信號
%?fid=fopen（‘b.mat‘‘r‘）;%故障
N=10240;
%?xdata=fread（fidN‘int16‘）;%fread：以二進制形式讀取文件中數據，int16：16字節數，整數型
%?fclose（fid）;
load?X156_DE_time;
xdata=X156_DE_time;
xdata=xdata（1:10240）;
xdata=（xdata-mean（xdata））/std（xdata1）;%mean平均值?std標準偏差???對原始信號進行歸一處理

%時域波形
figure（1）;
subplot（211）;
plot（1:Nxdata）;
xlabel（‘時間?t/n‘）;
ylabel（‘電壓?V/v‘）;
%原始信號的功率譜
nfft=10240;
Y1=abs（fft（xdatanfft））;
subplot（212）;
plot（（0:nfft/2-1）/nfft*fsY1（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;

%db10小波進行4層分解
%一維小波分解
[cl]?=?wavedec（xdata4‘db8‘）;

%重構第1～4層細節信號
d4?=?wrcoef（‘d‘cl‘db10‘4）;
d3?=?wrcoef（‘d‘cl‘db10‘3）;
d2?=?wrcoef（‘d‘cl‘db10‘2）;
d1?=?wrcoef（‘d‘cl‘db10‘1）;

%顯示細節信號
figure（2）
title（‘細節信號‘）;
subplot（411）;
plot（d1‘LineWidth‘2）;
ylabel（‘d1‘）;
subplot（412）;
plot（d2‘LineWidth‘2）;
ylabel（‘d2‘）;
subplot（413）;
plot（d3‘LineWidth‘2）;
ylabel（‘d3‘）;
subplot（414）;
plot（d4‘LineWidth‘2）;
ylabel（‘d4‘）;
xlabel（‘時間?t/s‘）;

%細節信號功率譜
figure（3）
title（‘細節信號功率譜‘）;
%第一層細節信號的功率譜
s1=abs（fft（d1nfft））;
subplot（411）;
plot（（0:nfft/2-1）/nfft*fss1（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;
%第二層細節信號的功率譜
s2=abs（fft（d2nfft））;
subplot（412）;
plot（（0:nfft/2-1）/nfft*fss2（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;
%第三層細節信號的功率譜
s3=abs（fft（d3nfft））;
subplot（413）;
plot（（0:nfft/2-1）/nfft*fss3（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;
%第四層細節信號的功率譜
s4=abs（fft（d4nfft））;
subplot（414）;
plot（（0:nfft/2-1）/nfft*fss4（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;ylabel（‘功率譜?P/W‘）;


%第1層細節信號的包絡譜
figure（4）;
title（‘細節信號包絡圖‘）;
y=hilbert（d1）;
?y=y-mean（y）;
ydata=abs（y）;?
%ydata=abs（y）;?%求絕對值
%y=y-mean（y）;
nfft=10240;
p=abs（fft（ydatanfft））;
subplot（411）;
plot（（0:nfft/2-1）/nfft*fsp（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;

%第2層細節信號的包絡譜
y=hilbert（d2）;
?ydata=abs（y）;?%求絕對值
%?y=y-mean（y）;
nfft=10240;
p=abs（fft（ydatanfft））;
%figure（5）;
subplot（412）;
plot（（0:nfft/2-1）/nfft*fsp（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;

%第3層細節信號的包絡譜
y=hilbert（d3）;
?y=y-mean（y）;
?ydata=abs（y）;?%求絕對值
%?y=y-mean（y）;
nfft=10240;
p=abs（fft（ydatanfft））;
%figure（6）;
subplot（413）;
plot（（0:nfft/2-1）/nfft*fsp（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;

%第4層細節信號的包絡譜
y=hilbert（d4）;
?ydata=abs（y）;?%求絕對值
%?y=y-mean（y）;
nfft=10240;
p=abs（fft（ydatanfft））;
%figure（7）;
subplot（414）;
plot（（0:nfft/2-1）*0.1/nfft*fsp（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????3153??2011-12-25?22:17??zhouchengxinhaochuli.m