轴承振动信号的归一化处理

这是我编写的关于轴承振动信号归一化处理的程序

% 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