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

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

%時(shí)域波形
figure（1）;
subplot（211）;
plot（1:Nxdata）;
xlabel（‘時(shí)間?t/n‘）;
ylabel（‘電壓?V/v‘）;
%原始信號(hào)的功率譜
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小波進(jìn)行4層分解
%一維小波分解
[cl]?=?wavedec（xdata4‘db8‘）;

%重構(gòu)第1～4層細(xì)節(jié)信號(hào)
d4?=?wrcoef（‘d‘cl‘db10‘4）;
d3?=?wrcoef（‘d‘cl‘db10‘3）;
d2?=?wrcoef（‘d‘cl‘db10‘2）;
d1?=?wrcoef（‘d‘cl‘db10‘1）;

%顯示細(xì)節(jié)信號(hào)
figure（2）
title（‘細(xì)節(jié)信號(hào)‘）;
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（‘時(shí)間?t/s‘）;

%細(xì)節(jié)信號(hào)功率譜
figure（3）
title（‘細(xì)節(jié)信號(hào)功率譜‘）;
%第一層細(xì)節(jié)信號(hào)的功率譜
s1=abs（fft（d1nfft））;
subplot（411）;
plot（（0:nfft/2-1）/nfft*fss1（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;
%第二層細(xì)節(jié)信號(hào)的功率譜
s2=abs（fft（d2nfft））;
subplot（412）;
plot（（0:nfft/2-1）/nfft*fss2（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;
%第三層細(xì)節(jié)信號(hào)的功率譜
s3=abs（fft（d3nfft））;
subplot（413）;
plot（（0:nfft/2-1）/nfft*fss3（1:nfft/2））;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘功率譜?P/W‘）;
%第四層細(xì)節(jié)信號(hào)的功率譜
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層細(xì)節(jié)信號(hào)的包絡(luò)譜
figure（4）;
title（‘細(xì)節(jié)信號(hào)包絡(luò)圖‘）;
y=hilbert（d1）;
?y=y-mean（y）;
ydata=abs（y）;?
%ydata=abs（y）;?%求絕對(duì)值
%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層細(xì)節(jié)信號(hào)的包絡(luò)譜
y=hilbert（d2）;
?ydata=abs（y）;?%求絕對(duì)值
%?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層細(xì)節(jié)信號(hào)的包絡(luò)譜
y=hilbert（d3）;
?y=y-mean（y）;
?ydata=abs（y）;?%求絕對(duì)值
%?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層細(xì)節(jié)信號(hào)的包絡(luò)譜
y=hilbert（d4）;
?ydata=abs（y）;?%求絕對(duì)值
%?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‘）;

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