%正常工作状况下
%用db5小波对信号进行6层分解
wpt=wpdec(iz,6,'db5','shannon');
%提取频段系数
cc1=wprcoef(wpt,[6 0]);
cc2=wprcoef(wpt,[6 1]);
cc3=wprcoef(wpt,[6 2]);
cc4=wprcoef(wpt,[6 3]);
cc5=wprcoef(wpt,[6 4]);
cc6=wprcoef(wpt,[6 5]);
cc7=wprcoef(wpt,[6 6]);
cc8=wprcoef(wpt,[6 7]);
cc9=wprcoef(wpt,[6 8]);
cc10=wprcoef(wpt,[6 9]);
cc11=wprcoef(wpt,[6 10]);
cc12=wprcoef(wpt,[6 11]);
cc13=wprcoef(wpt,[6 12]);
cc14=wprcoef(wpt,[6 13]);
cc15=wprcoef(wpt,[6 14]);
cc16=wprcoef(wpt,[6 15]);
cc17=wprcoef(wpt,[6 16]);
cc18=wprcoef(wpt,[6 17]);
cc19=wprcoef(wpt,[6 18]);
cc20=wprcoef(wpt,[6 19]);
cc21=wprcoef(wpt,[6 20]);
cc22=wprcoef(wpt,[6 21]);
cc23=wprcoef(wpt,[6 22]);
cc24=wprcoef(wpt,[6 23]);
cc25=wprcoef(wpt,[6 24]);
cc26=wprcoef(wpt,[6 25]);
cc27=wprcoef(wpt,[6 26]);
cc28=wprcoef(wpt,[6 27]);
cc29=wprcoef(wpt,[6 28]);
cc30=wprcoef(wpt,[6 29]);
cc31=wprcoef(wpt,[6 30]);
cc32=wprcoef(wpt,[6 31]);
cc33=wprcoef(wpt,[6 32]);
cc34=wprcoef(wpt,[6 33]);
cc35=wprcoef(wpt,[6 34]);
cc36=wprcoef(wpt,[6 35]);
cc37=wprcoef(wpt,[6 36]);
cc38=wprcoef(wpt,[6 37]);
cc39=wprcoef(wpt,[6 38]);
cc40=wprcoef(wpt,[6 39]);
cc41=wprcoef(wpt,[6 40]);
cc42=wprcoef(wpt,[6 41]);
cc43=wprcoef(wpt,[6 42]);
cc44=wprcoef(wpt,[6 43]);
cc45=wprcoef(wpt,[6 44]);
cc46=wprcoef(wpt,[6 45]);
cc47=wprcoef(wpt,[6 46]);
cc48=wprcoef(wpt,[6 47]);
cc49=wprcoef(wpt,[6 48]);
cc50=wprcoef(wpt,[6 49]);
cc51=wprcoef(wpt,[6 50]);
cc52=wprcoef(wpt,[6 51]);
cc53=wprcoef(wpt,[6 52]);
cc54=wprcoef(wpt,[6 53]);
cc55=wprcoef(wpt,[6 54]);
cc56=wprcoef(wpt,[6 55]);
cc57=wprcoef(wpt,[6 56]);
cc58=wprcoef(wpt,[6 57]);
cc59=wprcoef(wpt,[6 58]);
cc60=wprcoef(wpt,[6 59]);
cc61=wprcoef(wpt,[6 60]);
cc62=wprcoef(wpt,[6 61]);
cc63=wprcoef(wpt,[6 62]);
cc64=wprcoef(wpt,[6 63]);
c1=wpcoef(wpt,[6,0]);
c2=wpcoef(wpt,[6,1]);
c3=wpcoef(wpt,[6,2]);
c4=wpcoef(wpt,[6,3]);
c5=wpcoef(wpt,[6,4]);
c6=wpcoef(wpt,[6,5]);
c7=wpcoef(wpt,[6,6]);
c8=wpcoef(wpt,[6,7]);
c9=wpcoef(wpt,[6,8]);
c10=wpcoef(wpt,[6,9]);
c11=wpcoef(wpt,[6,10]);
c12=wpcoef(wpt,[6,11]);
c13=wpcoef(wpt,[6,12]);
c14=wpcoef(wpt,[6,13]);
c15=wpcoef(wpt,[6,14]);
c16=wpcoef(wpt,[6,15]);
c17=wpcoef(wpt,[6,16]);
c18=wpcoef(wpt,[6,17]);
c19=wpcoef(wpt,[6,18]);
c20=wpcoef(wpt,[6,19]);
c21=wpcoef(wpt,[6,20]);
c22=wpcoef(wpt,[6,21]);
c23=wpcoef(wpt,[6,22]);
c24=wpcoef(wpt,[6,23]);
c25=wpcoef(wpt,[6,24]);
c26=wpcoef(wpt,[6,25]);
c27=wpcoef(wpt,[6,26]);
c28=wpcoef(wpt,[6,27]);
c29=wpcoef(wpt,[6,28]);
c30=wpcoef(wpt,[6,29]);
c31=wpcoef(wpt,[6,30]);
c32=wpcoef(wpt,[6,31]);
c33=wpcoef(wpt,[6,32]);
c34=wpcoef(wpt,[6,33]);
c35=wpcoef(wpt,[6,34]);
c36=wpcoef(wpt,[6,35]);
c37=wpcoef(wpt,[6,36]);
c38=wpcoef(wpt,[6,37]);
c39=wpcoef(wpt,[6,38]);
c40=wpcoef(wpt,[6,39]);
c41=wpcoef(wpt,[6,40]);
c42=wpcoef(wpt,[6,41]);
c43=wpcoef(wpt,[6,42]);
c44=wpcoef(wpt,[6,43]);
c45=wpcoef(wpt,[6,44]);
c46=wpcoef(wpt,[6,45]);
c47=wpcoef(wpt,[6,46]);
c48=wpcoef(wpt,[6,47]);
c49=wpcoef(wpt,[6,48]);
c50=wpcoef(wpt,[6,49]);
c51=wpcoef(wpt,[6,50]);
c52=wpcoef(wpt,[6,51]);
c53=wpcoef(wpt,[6,52]);
c54=wpcoef(wpt,[6,53]);
c55=wpcoef(wpt,[6,54]);
c56=wpcoef(wpt,[6,55]);
c57=wpcoef(wpt,[6,56]);
c58=wpcoef(wpt,[6,57]);
c59=wpcoef(wpt,[6,58]);
c60=wpcoef(wpt,[6,59]);
c61=wpcoef(wpt,[6,60]);
c62=wpcoef(wpt,[6,61]);
c63=wpcoef(wpt,[6,62]);
c64=wpcoef(wpt,[6,63]);
xs=[c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11 c12 c13 c14 c15 c16 c17 c18 c19 c20 c21 c22 c23 c24 c25 c26 c27 c28 c29 c30 c31 c32 c33 c34 c35 c36 c37 c38 c39 c40 c41 c42 c43 c44 c45 c46 c47 c48 c49 c50 c51 c52 c53 c54 c55 c56 c57 c58 c59 c60 c61 c62 c63 c64];
Ed1=sum(c1.^2);
Ed2=sum(c2.^2);
Ed3=sum(c3.^2);
Ed4=sum(c4.^2);
Ed5=sum(c5.^2);
Ed6=sum(c6.^2);
Ed7=sum(c7.^2);
Ed8=sum(c8.^2);
Ed9=sum(c9.^2);
Ed10=sum(c10.^2);
Ed11=sum(c11.^2);
Ed12=sum(c12.^2);
Ed13=sum(c13.^2);
Ed14=sum(c14.^2);
Ed15=sum(c15.^2);
Ed16=sum(c16.^2);
Ed17=sum(c17.^2);
Ed18=sum(c18.^2);
Ed19=sum(c19.^2);
Ed20=sum(c20.^2);
Ed21=sum(c21.^2);
Ed22=sum(c22.^2);
Ed23=sum(c23.^2);
Ed24=sum(c24.^2);
Ed25=sum(c25.^2);
Ed26=sum(c26.^2);
Ed27=sum(c27.^2);
Ed28=sum(c28.^2);
Ed29=sum(c29.^2);
Ed30=sum(c30.^2);
Ed31=sum(c31.^2);
Ed32=sum(c32.^2);
Ed33=sum(c33.^2);
Ed34=sum(c34.^2);
Ed35=sum(c35.^2);
Ed36=sum(c36.^2);
Ed37=sum(c37.^2);
Ed38=sum(c38.^2);
Ed39=sum(c39.^2);
Ed40=sum(c40.^2);
Ed41=sum(c41.^2);
Ed42=sum(c42.^2);
Ed43=sum(c43.^2);
Ed44=sum(c44.^2);
Ed45=sum(c45.^2);
Ed46=sum(c46.^2);
Ed47=sum(c47.^2);
Ed48=sum(c48.^2);
Ed49=sum(c49.^2);
Ed50=sum(c50.^2);
Ed51=sum(c51.^2);
Ed52=sum(c52.^2);
Ed53=sum(c53.^2);
Ed54=sum(c54.^2);
Ed55=sum(c55.^2);
Ed56=sum(c56.^2);
Ed57=sum(c57.^2);
Ed58=sum(c58.^2);
Ed59=sum(c59.^2);
Ed60=sum(c60.^2);
Ed61=sum(c61.^2);
Ed62=sum(c62.^2);
Ed63=sum(c63.^2);
Ed64=sum(c64.^2);
Ed=[Ed1 Ed2 Ed3 Ed4 Ed5 Ed6 Ed7 Ed8 Ed9 Ed10 Ed11 Ed12 Ed13 Ed14 Ed15 Ed16 Ed17 Ed18 Ed19 Ed20 Ed21 Ed22 Ed23 Ed24 Ed25 Ed26 Ed27 Ed28 Ed29 Ed30 Ed31 Ed32 Ed33 Ed34 Ed35 Ed36 Ed37 Ed38 Ed39 Ed40 Ed41 Ed42 Ed43 Ed44 Ed45 Ed46 Ed47 Ed48 Ed49 Ed50 Ed51 Ed52 Ed53 Ed54 Ed55 Ed56 Ed57 Ed58 Ed59 Ed60 Ed61 Ed62 Ed63 Ed64];
Eall=Ed1+ Ed2 +Ed3+ Ed4 +Ed5 +Ed6 +Ed7 +Ed8 +Ed9 +Ed10 +Ed11 +Ed12 +Ed13 +Ed14 +Ed15 +Ed16 +Ed17 +Ed18 +Ed19 +Ed20 +Ed21 +Ed22 +Ed23 +Ed24 +Ed25 +Ed26 +Ed27 +Ed28 +Ed29 +Ed30 +Ed31 + Ed32 +Ed33 +Ed34 +Ed35 +Ed36 +Ed37 +Ed38 +Ed39 +Ed40 +Ed41 +Ed42 +Ed43 +Ed44 +Ed45 +Ed46 +Ed47 +Ed48 +Ed49 +Ed50 +Ed51 +Ed52 +Ed53 +Ed54 +Ed55 +Ed56 +Ed57 +Ed58 +Ed59 +Ed60 +Ed61 + Ed62 +Ed63 +Ed64;
power1=[Ed1 Ed2 Ed3 Ed4 Ed5 Ed6 Ed7 Ed8 Ed9 Ed10 Ed11 Ed12 Ed13 Ed14 Ed15 Ed16 Ed17 Ed18 Ed19 Ed20 Ed21 Ed22 Ed23 Ed24 Ed25 Ed26 Ed27 Ed28 Ed29 Ed30 Ed31 Ed32 Ed33 Ed34 Ed35 Ed36 Ed37 Ed38 Ed39 Ed40 Ed41 Ed42 Ed43 Ed44 Ed45 Ed46 Ed47 Ed48 Ed49 Ed50 Ed51 Ed52 Ed53 Ed54 Ed55 Ed56 Ed57 Ed58 Ed59 Ed60 Ed61 Ed62 Ed63 Ed64];
PE1=Ed1/Eall;
PE2=Ed2/Eall;
PE3=Ed3/Eall;
PE4=Ed4/Eall;
PE5=Ed5/Eall;
PE6=Ed6/Eall;
PE7=Ed7/Eall;
PE8=Ed8/Eall;
PE9=Ed9/Eall;
PE10=Ed10/Eall;
PE11=Ed11/Eall;
PE12=Ed12/Eall;
PE13=Ed13/Eall;
PE14=Ed14/Eall;
PE15=Ed15/Eall;
PE16=Ed16/Eall;
PE17=Ed17/Eall;
PE18=Ed18/Eall;
PE19=Ed19/Eall;
PE20=Ed20/Eall;
PE21=Ed21/Eall;
PE22=Ed22/Eall;
PE23=Ed23/Eall;
PE24=Ed24/Eall;
PE25=Ed25/Eall;
PE26=Ed26/Eall;
PE27=Ed27/Eall;
PE28=Ed28/Eall;
PE29=Ed29/Eall;
PE30=Ed30/Eall;
PE31=Ed31/Eall;
PE32=Ed32/Eall;
PE33=Ed33/Eall;
PE34=Ed34/Eall;
PE35=Ed35/Eall;
PE36=Ed36/Eall;
PE37=Ed37/Eall;
PE38=Ed38/Eall;
PE39=Ed39/Eall;
PE40=Ed40/Eall;
PE41=Ed41/Eall;
PE42=Ed42/Eall;
PE43=Ed43/Eall;
PE44=Ed44/Eall;
PE45=Ed45/Eall;
PE46=Ed46/Eall;
PE47=Ed47/Eall;
PE48=Ed48/Eall;
PE49=Ed49/Eall;
PE50=Ed50/Eall;
PE51=Ed51/Eall;
PE52=Ed52/Eall;
PE53=Ed53/Eall;
PE54=Ed54/Eall;
PE55=Ed55/Eall;
PE56=Ed56/Eall;
PE57=Ed57/Eall;
PE58=Ed58/Eall;
PE59=Ed59/Eall;
PE60=Ed60/Eall;
PE61=Ed61/Eall;
PE62=Ed62/Eall;
PE63=Ed63/Eall;
PE64=Ed64/Eall;
PE=[PE1 PE2 PE3 PE4 PE5 PE6 PE7 PE8 PE9 PE10 PE11 PE12 PE13 PE14 PE15 PE16 PE17 PE18 PE19 PE20 PE21 PE22 PE23 PE24 PE25 PE26 PE27 PE28 PE29 PE30 PE31 PE32 PE33 PE34 PE35 PE36 PE37 PE38 PE39 PE40 PE41 PE42 PE43 PE44 PE45 PE46 PE47 PE48 PE49 PE50 PE51 PE52 PE53 PE54 PE55 PE56 PE57 PE58 PE59 PE60 PE61 PE62 PE63 PE64];
T=PE;
disp('能量分布比')
%故障电弧状况下
%用db5小波对信号进行6层分解
wpt=wpdec(ic,6,'db5','shannon');
%提取频段系数
cc1=wprcoef(wpt,[6 0]);
cc2=wprcoef(wpt,[6 1]);
cc3=wprcoef(wpt,[6 2]);
cc4=wprcoef(wpt,[6 3]);
cc5=wprcoef(wpt,[6 4]);
cc6=wprcoef(wpt,[6 5]);
cc7=wprcoef(wpt,[6 6]);
cc8=wprcoef(wpt,[6 7]);
cc9=wprcoef(wpt,[6 8]);
cc10=wprcoef(wpt,[6 9]);
cc11=wprcoef(wpt,[6 10]);
cc12=wprcoef(wpt,[6 11]);
cc13=wprcoef(wpt,[6 12]);
cc14=wprcoef(wpt,[6 13]);
cc15=wprcoef(wpt,[6 14]);
cc16=wprcoef(wpt,[6 15]);
cc17=wprcoef(wpt,[6 16]);
cc18=wprcoef(wpt,[6 17]);
cc19=wprcoef(wpt,[6 18]);
cc20=wprcoef(wpt,[6 19]);
cc21=wprcoef(wpt,[6 20]);
评论1