waveletana.m

内燃机转子仿真

%wavelet 
clear all;
clc;
close all;
format long;
load 'bneiranji';
load 'bneiranji1';
load 'bneiranji2';
load 'bneiranji3';
load 'bneiranji4';
ip=7;
sig=xyz(ip,:);
sig1=xyz1(ip,:);
sig2=xyz2(ip,:);
sig3=xyz3(ip,:);
sig4=xyz4(ip,:);
[l_sig,h_sig]=wavedec(sig,5,'db3');
[l_sig1,h_sig1]=wavedec(sig1,5,'db3');
[l_sig2,h_sig2]=wavedec(sig2,5,'db3');
[l_sig3,h_sig3]=wavedec(sig3,5,'db3');
[l_sig4,h_sig4]=wavedec(sig4,5,'db3');
typewave=1;
if typewave==1

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

dsig1_5=wrcoef('d',l_sig,h_sig,'db3',5);
dsig1_4=wrcoef('d',l_sig,h_sig,'db3',4);
dsig1_3=wrcoef('d',l_sig,h_sig,'db3',3);
dsig1_2=wrcoef('d',l_sig,h_sig,'db3',2);
dsig1_1=wrcoef('d',l_sig,h_sig,'db3',1);

dsig2_5=wrcoef('d',l_sig,h_sig,'db3',5);
dsig2_4=wrcoef('d',l_sig,h_sig,'db3',4);
dsig2_3=wrcoef('d',l_sig,h_sig,'db3',3);
dsig2_2=wrcoef('d',l_sig,h_sig,'db3',2);
dsig2_1=wrcoef('d',l_sig,h_sig,'db3',1);

dsig3_5=wrcoef('d',l_sig,h_sig,'db3',5);
dsig3_4=wrcoef('d',l_sig,h_sig,'db3',4);
dsig3_3=wrcoef('d',l_sig,h_sig,'db3',3);
dsig3_2=wrcoef('d',l_sig,h_sig,'db3',2);
dsig3_1=wrcoef('d',l_sig,h_sig,'db3',1);

dsig4_5=wrcoef('d',l_sig,h_sig,'db3',5);
dsig4_4=wrcoef('d',l_sig,h_sig,'db3',4);
dsig4_3=wrcoef('d',l_sig,h_sig,'db3',3);
dsig4_2=wrcoef('d',l_sig,h_sig,'db3',2);
dsig4_1=wrcoef('d',l_sig,h_sig,'db3',1);
else
 dsig_5=wrcoef('a',l_sig,h_sig,'db3',5);
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_5=wrcoef('a',l_sig,h_sig,'db3',5);
dsig1_4=wrcoef('a',l_sig,h_sig,'db3',4);
dsig1_3=wrcoef('a',l_sig,h_sig,'db3',3);
dsig1_2=wrcoef('a',l_sig,h_sig,'db3',2);
dsig1_1=wrcoef('a',l_sig,h_sig,'db3',1);

dsig2_5=wrcoef('a',l_sig,h_sig,'db3',5);
dsig2_4=wrcoef('a',l_sig,h_sig,'db3',4);
dsig2_3=wrcoef('a',l_sig,h_sig,'db3',3);
dsig2_2=wrcoef('a',l_sig,h_sig,'db3',2);
dsig2_1=wrcoef('a',l_sig,h_sig,'db3',1);

dsig3_5=wrcoef('a',l_sig,h_sig,'db3',5);
dsig3_4=wrcoef('a',l_sig,h_sig,'db3',4);
dsig3_3=wrcoef('a',l_sig,h_sig,'db3',3);
dsig3_2=wrcoef('a',l_sig,h_sig,'db3',2);
dsig3_1=wrcoef('a',l_sig,h_sig,'db3',1);

dsig4_5=wrcoef('a',l_sig,h_sig,'db3',5);
dsig4_4=wrcoef('a',l_sig,h_sig,'db3',4);
dsig4_3=wrcoef('a',l_sig,h_sig,'db3',3);
dsig4_2=wrcoef('a',l_sig,h_sig,'db3',2);
dsig4_1=wrcoef('a',l_sig,h_sig,'db3',1);
end

tu=4;
if tu==0
figure(1);
plot(t,dsig_5);
figure(2);
plot(t,dsig_4);
figure(3);
plot(t,dsig_3);
figure(4);
plot(t,dsig_2);
figure(5);
plot(t,dsig_1);
figure(6);
plot(t,sig);
elseif tu==1
 figure(1);
plot(t,dsig1_5);
figure(2);
plot(t,dsig1_4);
figure(3);
plot(t,dsig1_3);
figure(4);
plot(t,dsig1_2);
figure(5);
plot(t,dsig1_1);
figure(6);
plot(t,sig1);   
elseif tu==2
 figure(1);
plot(t,dsig2_5);
figure(2);
plot(t,dsig2_4);
figure(3);
plot(t,dsig2_3);
figure(4);
plot(t,dsig2_2);
figure(5);
plot(t,dsig2_1);
figure(6);
plot(t,sig2);   
elseif tu==3
 figure(1);
plot(t,dsig3_5);
figure(2);
plot(t,dsig3_4);
figure(3);
plot(t,dsig3_3);
figure(4);
plot(t,dsig3_2);
figure(5);
plot(t,dsig3_1);
figure(6);
plot(t,sig3);   
elseif tu==4
 figure(1);
plot(t,dsig4_5);
figure(2);
plot(t,dsig4_4);
figure(3);
plot(t,dsig4_3);
figure(4);
plot(t,dsig4_2);
figure(5);
plot(t,dsig4_1);
figure(6);
plot(t,sig4);   
end