waveletana1.m

内燃机转子仿真

%wavelet 
clear all;
clc;
close all;
format long;
load 'bneiranji';
load 'bneiranji1';
load 'bneiranji2';
load 'bneiranji3';
load 'bneiranji4';
load 'bneiranji5';
ip=7;
sig=xyz(ip,:)*180/pi;
sig1=xyz1(ip,:)*180/pi;
sig2=xyz2(ip,:)*180/pi;
sig3=xyz3(ip,:)*180/pi;
sig4=xyz4(ip,:)*180/pi;
sig5=xyz5(ip,:)*180/pi;
[l_sig,h_sig]=wavedec(sig,4,'db1');
[l_sig1,h_sig1]=wavedec(sig1,4,'db1');
[l_sig2,h_sig2]=wavedec(sig2,4,'db1');
[l_sig3,h_sig3]=wavedec(sig3,4,'db1');
[l_sig4,h_sig4]=wavedec(sig4,4,'db1');
[l_sig5,h_sig5]=wavedec(sig5,4,'db1');

typewave=2;
if typewave==1


dsig_4=wrcoef('d',l_sig,h_sig,'db1',4);
dsig_3=wrcoef('d',l_sig,h_sig,'db1',3);
dsig_2=wrcoef('d',l_sig,h_sig,'db1',2);
dsig_1=wrcoef('d',l_sig,h_sig,'db1',1);


dsig1_4=wrcoef('d',l_sig1,h_sig1,'db1',4);
dsig1_3=wrcoef('d',l_sig1,h_sig1,'db1',3);
dsig1_2=wrcoef('d',l_sig1,h_sig1,'db1',2);
dsig1_1=wrcoef('d',l_sig1,h_sig1,'db1',1);


dsig2_4=wrcoef('d',l_sig2,h_sig2,'db1',4);
dsig2_3=wrcoef('d',l_sig2,h_sig2,'db1',3);
dsig2_2=wrcoef('d',l_sig2,h_sig2,'db1',2);
dsig2_1=wrcoef('d',l_sig2,h_sig2,'db1',1);


dsig3_4=wrcoef('d',l_sig3,h_sig3,'db1',4);
dsig3_3=wrcoef('d',l_sig3,h_sig3,'db1',3);
dsig3_2=wrcoef('d',l_sig3,h_sig3,'db1',2);
dsig3_1=wrcoef('d',l_sig3,h_sig3,'db1',1);


dsig4_4=wrcoef('d',l_sig4,h_sig4,'db1',4);
dsig4_3=wrcoef('d',l_sig4,h_sig4,'db1',3);
dsig4_2=wrcoef('d',l_sig4,h_sig4,'db1',2);
dsig4_1=wrcoef('d',l_sig4,h_sig4,'db1',1);


dsig5_4=wrcoef('d',l_sig5,h_sig5,'db1',4);
dsig5_3=wrcoef('d',l_sig5,h_sig5,'db1',3);
dsig5_2=wrcoef('d',l_sig5,h_sig5,'db1',2);
dsig5_1=wrcoef('d',l_sig5,h_sig5,'db1',1);
else

dsig_4=wrcoef('a',l_sig,h_sig,'db3',4);
dsig_3=wrcoef('a',l_sig,h_sig,'db3',3);
dsig_2=wrcoef('a',l_sig,h_sig,'db3',2);
dsig_1=wrcoef('a',l_sig,h_sig,'db3',1);


dsig1_4=wrcoef('a',l_sig1,h_sig1,'db3',4);
dsig1_3=wrcoef('a',l_sig1,h_sig1,'db3',3);
dsig1_2=wrcoef('a',l_sig1,h_sig1,'db3',2);
dsig1_1=wrcoef('a',l_sig1,h_sig1,'db3',1);


dsig2_4=wrcoef('a',l_sig2,h_sig2,'db3',4);
dsig2_3=wrcoef('a',l_sig2,h_sig2,'db3',3);
dsig2_2=wrcoef('a',l_sig2,h_sig2,'db3',2);
dsig2_1=wrcoef('a',l_sig2,h_sig2,'db3',1);


dsig3_4=wrcoef('a',l_sig3,h_sig3,'db3',4);
dsig3_3=wrcoef('a',l_sig3,h_sig3,'db3',3);
dsig3_2=wrcoef('a',l_sig3,h_sig3,'db3',2);
dsig3_1=wrcoef('a',l_sig3,h_sig3,'db3',1);


dsig4_4=wrcoef('a',l_sig4,h_sig4,'db3',4);
dsig4_3=wrcoef('a',l_sig4,h_sig4,'db3',3);
dsig4_2=wrcoef('a',l_sig4,h_sig4,'db3',2);
dsig4_1=wrcoef('a',l_sig4,h_sig4,'db3',1);

dsig5_4=wrcoef('a',l_sig5,h_sig5,'db3',4);
dsig5_3=wrcoef('a',l_sig5,h_sig5,'db3',3);
dsig5_2=wrcoef('a',l_sig5,h_sig5,'db3',2);
dsig5_1=wrcoef('a',l_sig5,h_sig5,'db3',1);
end

nfft=2048;

dsig_hil=hilbert(dsig_4);
adsig_hil=abs(dsig_hil);
adsig_hil=adsig_hil-mean(adsig_hil);
pxx=abs(fft(adsig_hil,nfft));

dsig1_hil=hilbert(dsig1_4);
adsig1_hil=abs(dsig1_hil);
adsig1_hil=adsig1_hil-mean(adsig1_hil);
pxx1=abs(fft(adsig1_hil,nfft));

dsig2_hil=hilbert(dsig2_4);
adsig2_hil=abs(dsig2_hil);
adsig2_hil=adsig2_hil-mean(adsig2_hil);
pxx2=abs(fft(adsig2_hil,nfft));

dsig3_hil=hilbert(dsig3_4);
adsig3_hil=abs(dsig3_hil);
adsig3_hil=adsig3_hil-mean(adsig3_hil);
pxx3=abs(fft(adsig3_hil,nfft));

dsig4_hil=hilbert(dsig4_4);
adsig4_hil=abs(dsig4_hil);
adsig4_hil=adsig4_hil-mean(adsig4_hil);
pxx4=abs(fft(adsig4_hil,nfft));

dsig5_hil=hilbert(dsig5_4);
adsig5_hil=abs(dsig5_hil);
adsig5_hil=adsig4_hil-mean(adsig5_hil);
pxx5=abs(fft(adsig5_hil,nfft));

figure(1);
subplot(6,1,1);
plot(t,dsig_4);
subplot(6,1,2);
plot(t,dsig1_4);
subplot(6,1,3);
plot(t,dsig2_4);
subplot(6,1,4);
plot(t,dsig3_4);
subplot(6,1,5);
plot(t,dsig4_4);
subplot(6,1,6);
plot(t,dsig5_4);

figure(2);
subplot(6,1,1);
plot(t,dsig_3);
subplot(6,1,2);
plot(t,dsig1_3);
subplot(6,1,3);
plot(t,dsig2_3);
subplot(6,1,4);
plot(t,dsig3_3);
subplot(6,1,5);
plot(t,dsig4_3);
subplot(6,1,6);
plot(t,dsig5_3);

figure(3);
subplot(6,1,1);
plot(t,dsig_2);
subplot(6,1,2);
plot(t,dsig1_2);
subplot(6,1,3);
plot(t,dsig2_2);
subplot(6,1,4);
plot(t,dsig3_2);
subplot(6,1,5);
plot(t,dsig4_2);
subplot(6,1,6);
plot(t,dsig5_2);

figure(4);
subplot(6,1,1);
plot(t,dsig_1);
subplot(6,1,2);
plot(t,dsig1_1);
subplot(6,1,3);
plot(t,dsig2_1);
subplot(6,1,4);
plot(t,dsig3_1);
subplot(6,1,5);
plot(t,dsig4_1);
subplot(6,1,6);
plot(t,dsig5_1);

figure(5)
subplot(6,1,1);
plot((0:nfft/2-1)/nfft/h,pxx(1:nfft/2));
ylim([0 10]);
subplot(6,1,2);
plot((0:nfft/2-1)/nfft/h,pxx1(1:nfft/2));
ylim([0 10]);
subplot(6,1,3);
plot((0:nfft/2-1)/nfft/h,pxx2(1:nfft/2));
ylim([0 10]);
subplot(6,1,4);
plot((0:nfft/2-1)/nfft/h,pxx3(1:nfft/2));
ylim([0 10]);
subplot(6,1,5);
plot((0:nfft/2-1)/nfft/h,pxx4(1:nfft/2));
ylim([0 10]);
subplot(6,1,6);
plot((0:nfft/2-1)/nfft/h,pxx5(1:nfft/2));
ylim([0 10]);