xiaobobaochengxu.m

实现小波变换

%t=0.001:0.001:1; 
t=1:1000; 
s1=sin(2*pi*50*t*0.001)+sin(2*pi*120*t*0.001)+rand(1,length(t)); 

for t=1:500; 
s2(t)=sin(2*pi*50*t*0.001)+sin(2*pi*120*t*0.001)+rand(1,length(t)); 
end 

for t=501:1000; 
s2(t)=sin(2*pi*200*t*0.001)+sin(2*pi*120*t*0.001)+rand(1,length(t)); 
end 


subplot(9,2,1) 
plot(s1) 
title('原始信号') 
ylabel('S1') 

subplot(9,2,2) 
plot(s2) 
title('故障信号') 
ylabel('S2') 

wpt=wpdec(s1,3,'db1','shannon'); 

%plot(wpt); 
s130=wprcoef(wpt,[3,0]); 
s131=wprcoef(wpt,[3,1]); 
s132=wprcoef(wpt,[3,2]); 
s133=wprcoef(wpt,[3,3]); 
s134=wprcoef(wpt,[3,4]); 
s135=wprcoef(wpt,[3,5]); 
s136=wprcoef(wpt,[3,6]); 
s137=wprcoef(wpt,[3,7]); 
s10=norm(s130); 
s11=norm(s131); 
s12=norm(s132); 
s13=norm(s133); 
s14=norm(s134); 
s15=norm(s135); 
s16=norm(s136); 
s17=norm(s137); 

st10=std(s130); 
st11=std(s131); 
st12=std(s132); 
st13=std(s133); 
st14=std(s134); 
st15=std(s135); 
st16=std(s136); 
st17=std(s137); 

disp('正常信号的特征向量'); 
snorm1=[s10,s11,s12,s13,s14,s15,s16,s17] 
std1=[st10,st11,st12,st13,st14,st15,st16,st17] 

subplot(9,2,3);plot(s130); 
ylabel('S130'); 
subplot(9,2,5);plot(s131); 
ylabel('S131'); 
subplot(9,2,7);plot(s132); 
ylabel('S132'); 
subplot(9,2,9);plot(s133); 
ylabel('S133'); 
subplot(9,2,11);plot(s134); 
ylabel('S134'); 
subplot(9,2,13);plot(s135); 
ylabel('S135'); 
subplot(9,2,15);plot(s136); 
ylabel('S136'); 
subplot(9,2,17);plot(s137); 
ylabel('S137'); 


wpt=wpdec(s2,3,'db1','shannon'); 

%plot(wpt); 
s230=wprcoef(wpt,[3,0]); 
s231=wprcoef(wpt,[3,1]); 
s232=wprcoef(wpt,[3,2]); 
s233=wprcoef(wpt,[3,3]); 
s234=wprcoef(wpt,[3,4]); 
s235=wprcoef(wpt,[3,5]); 
s236=wprcoef(wpt,[3,6]); 
s237=wprcoef(wpt,[3,7]); 

s20=norm(s230); 
s21=norm(s231); 
s22=norm(s232); 
s23=norm(s233); 
s24=norm(s234); 
s25=norm(s235); 
s26=norm(s236); 
s27=norm(s237); 

st20=std(s230); 
st21=std(s231); 
st22=std(s232); 
st23=std(s233); 
st24=std(s234); 
st25=std(s235); 
st26=std(s236); 
st27=std(s237); 

disp('故障信号的特征向量'); 
snorm2=[s20,s21,s22,s23,s24,s25,s26,s27] 
std2=[st20,st21,st22,st23,st24,st25,st26,st27] 


subplot(9,2,4);plot(s230); 
ylabel('S230'); 
subplot(9,2,6);plot(s231); 
ylabel('S231'); 
subplot(9,2,8);plot(s232); 
ylabel('S232'); 
subplot(9,2,10);plot(s233); 
ylabel('S233'); 
subplot(9,2,12);plot(s234); 
ylabel('S234'); 
subplot(9,2,14);plot(s235); 
ylabel('S235'); 
subplot(9,2,16);plot(s236); 
ylabel('S236'); 
subplot(9,2,18);plot(s237); 
ylabel('S237'); 

%fft 

figure 
y1=fft(s1,1024); 
py1=y1.*conj(y1)/1024; 
y2=fft(s2,1024); 
py2=y2.*conj(y2)/1024; 

y130=fft(s130,1024); 
py130=y130.*conj(y130)/1024; 
y131=fft(s131,1024); 
py131=y131.*conj(y131)/1024; 
y132=fft(s132,1024); 
py132=y132.*conj(y132)/1024; 
y133=fft(s133,1024); 
py133=y133.*conj(y133)/1024; 
y134=fft(s134,1024); 
py134=y134.*conj(y134)/1024; 
y135=fft(s135,1024); 
py135=y135.*conj(y135)/1024; 
y136=fft(s136,1024); 
py136=y136.*conj(y136)/1024; 
y137=fft(s137,1024); 
py137=y137.*conj(y137)/1024; 

y230=fft(s230,1024); 
py230=y230.*conj(y230)/1024; 
y231=fft(s231,1024); 
py231=y231.*conj(y231)/1024; 
y232=fft(s232,1024); 
py232=y232.*conj(y232)/1024; 
y233=fft(s233,1024); 
py233=y233.*conj(y233)/1024; 
y234=fft(s234,1024); 
py234=y234.*conj(y234)/1024; 
y235=fft(s235,1024); 
py235=y235.*conj(y235)/1024; 
y236=fft(s236,1024); 
py236=y236.*conj(y236)/1024; 
y237=fft(s237,1024); 
py237=y237.*conj(y237)/1024; 

f=1000*(0:511)/1024; 
subplot(1,2,1); 
plot(f,py1(1:512)); 
ylabel('P1'); 
title('原始信号的功率谱') 
subplot(1,2,2); 
plot(f,py2(1:512)); 
ylabel('P2'); 
title('故障信号的功率谱') 

figure 

subplot(4,2,1); 
plot(f,py130(1:512)); 
ylabel('P130'); 
title('S130的功率谱') 
subplot(4,2,2); 
plot(f,py131(1:512)); 
ylabel('P131'); 
title('S131的功率谱') 
subplot(4,2,3); 
plot(f,py132(1:512)); 
ylabel('P132'); 
subplot(4,2,4); 
plot(f,py133(1:512)); 
ylabel('P133'); 
subplot(4,2,5); 
plot(f,py134(1:512)); 
ylabel('P134'); 
subplot(4,2,6); 
plot(f,py135(1:512)); 
ylabel('P135'); 
subplot(4,2,7); 
plot(f,py136(1:512)); 
ylabel('P136'); 
subplot(4,2,8); 
plot(f,py137(1:512)); 
ylabel('P137'); 
figure 

subplot(4,2,1); 
plot(f,py230(1:512)); 
ylabel('P230'); 
title('S230的功率谱') 
subplot(4,2,2); 
plot(f,py231(1:512)); 
ylabel('P231'); 
title('S231的功率谱') 
subplot(4,2,3); 
plot(f,py232(1:512)); 
ylabel('P232'); 
subplot(4,2,4); 
plot(f,py233(1:512)); 
ylabel('P233'); 
subplot(4,2,5); 
plot(f,py234(1:512)); 
ylabel('P234'); 
subplot(4,2,6); 
plot(f,py235(1:512)); 
ylabel('P235'); 
subplot(4,2,7); 
plot(f,py236(1:512)); 
ylabel('P236'); 
subplot(4,2,8); 
plot(f,py237(1:512)); 
ylabel('P237'); 
figure 
%plottree(wpt)