feature_extract_1.m

利用MATLAB实现一维信号时间序列的

load L05_G19;
middle=reshape(L05_G19,112,2001);
s1=0;
s2=s1;
s3=s1;
s4=s1;

t1=wpdec(s1,3,'db1','shannon');
s130=wprcoef(t1,[3,0]);
s131=wprcoef(t1,[3,1]);
s132=wprcoef(t1,[3,2]);
s133=wprcoef(t1,[3,3]);
s134=wprcoef(t1,[3,4]);
s135=wprcoef(t1,[3,5]);
s136=wprcoef(t1,[3,6]);
s137=wprcoef(t1,[3,7]);
hold on
figure(1)
subplot(4,2,1);plot(s130);
Ylabel('s130');
subplot(4,2,2);plot(s131);
Ylabel('s131');
subplot(4,2,3);plot(s132);
Ylabel('s132');
subplot(4,2,4);plot(s133);
Ylabel('s133');
subplot(4,2,5);plot(s134);
Ylabel('s134');
subplot(4,2,6);plot(s135);
Ylabel('s135');
subplot(4,2,7);plot(s136);
Ylabel('s136');
subplot(4,2,8);plot(s137);
Ylabel('s137');
hold off
s10=norm(s130);
s11=norm(s131);
s12=norm(s132);
s13=norm(s133);
s14=norm(s134);
s15=norm(s135);
s16=norm(s136);
s17=norm(s137);
ss1=[s10,s11,s12,s13,s14,s15,s16,s17]
%%
t2=wpdec(s2,3,'db1','shannon');
s230=wprcoef(t2,[3,0]);
s231=wprcoef(t2,[3,1]);
s232=wprcoef(t2,[3,2]);
s233=wprcoef(t2,[3,3]);
s234=wprcoef(t2,[3,4]);
s235=wprcoef(t2,[3,5]);
s236=wprcoef(t2,[3,6]);
s237=wprcoef(t2,[3,7]);
hold on
figure(2)
subplot(4,2,1);plot(s230);
Ylabel('s230');
subplot(4,2,2);plot(s231);
Ylabel('s231');
subplot(4,2,3);plot(s232);
Ylabel('s232');
subplot(4,2,4);plot(s233);
Ylabel('s233');
subplot(4,2,5);plot(s234);
Ylabel('s234');
subplot(4,2,6);plot(s135);
Ylabel('s235');
subplot(4,2,7);plot(s236);
Ylabel('s236');
subplot(4,2,8);plot(s237);
Ylabel('s237');
hold off
s20=norm(s230);
s21=norm(s231);
s22=norm(s232);
s23=norm(s233);
s24=norm(s234);
s25=norm(s235);
s26=norm(s236);
s27=norm(s237);
ss2=[s20,s21,s22,s23,s24,s25,s26,s27]
%%
t3=wpdec(s3,3,'db1','shannon');
s330=wprcoef(t3,[3,0]);
s331=wprcoef(t3,[3,1]);
s332=wprcoef(t3,[3,2]);
s333=wprcoef(t3,[3,3]);
s334=wprcoef(t3,[3,4]);
s335=wprcoef(t3,[3,5]);
s336=wprcoef(t3,[3,6]);
s337=wprcoef(t3,[3,7]);
hold on 
figure(3)
subplot(4,2,1);plot(s330);
Ylabel('s330');
subplot(4,2,2);plot(s331);
Ylabel('s331');
subplot(4,2,3);plot(s332);
Ylabel('s332');
subplot(4,2,4);plot(s333);
Ylabel('s333');
subplot(4,2,5);plot(s334);
Ylabel('s334');
subplot(4,2,6);plot(s335);
Ylabel('s335');
subplot(4,2,7);plot(s336);
Ylabel('s336');
subplot(4,2,8);plot(s337);
Ylabel('s337');
hold off
s30=norm(s330);
s31=norm(s331);
s32=norm(s332);
s33=norm(s333);
s34=norm(s334);
s35=norm(s335);
s36=norm(s336);
s37=norm(s337);
ss3=[s30,s31,s32,s33,s34,s35,s36,s37]
%%
t4=wpdec(s4,3,'db1','shannon');
s430=wprcoef(t4,[3,0]);
s431=wprcoef(t4,[3,1]);
s432=wprcoef(t4,[3,2]);
s433=wprcoef(t4,[3,3]);
s434=wprcoef(t4,[3,4]);
s435=wprcoef(t4,[3,5]);
s436=wprcoef(t4,[3,6]);
s437=wprcoef(t4,[3,7]);
hold on
figure(4)
subplot(4,2,1);plot(s430);
Ylabel('s430');
subplot(4,2,2);plot(s431);
Ylabel('s431');
subplot(4,2,3);plot(s432);
Ylabel('s432');
subplot(4,2,4);plot(s433);
Ylabel('s433');
subplot(4,2,5);plot(s434);
Ylabel('s434');
subplot(4,2,6);plot(s435);
Ylabel('s435');
subplot(4,2,7);plot(s436);
Ylabel('s436');
subplot(4,2,8);plot(s437);
Ylabel('s437');
hold off
s40=norm(s430);
s41=norm(s431);
s42=norm(s432);
s43=norm(s433);
s44=norm(s434);
s45=norm(s435);
s46=norm(s436);
s47=norm(s437);
ss4=[s40,s41,s42,s43,s44,s45,s46,s47]
%%
figure(5)
hold on 
plot(s137,'r--')
plot(s337,'b-.')
hold off