matlab_信号DB值计算，提取整体趋势

close all clear all clc cd('C:\Program Files\MATLAB71\work\2nd_test');%%%%%%%%%%%读数据 Allname=struct2cell(dir); %%%%%%%%%%%结构体 [m,n]=size(Allname); for i=1:n-2 Imdata=importdata(Allname{1,i+2}); DATA=Imdata;%%%%%%%%%%%%%%%%%%%%%%%%将第i个文件数据赋给x 1表示第一个通道 N=20480; DATA=DATA(1:N ,:);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%取x中的前10240个点进行分析 fs=20000;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%采样频率 x=DATA(:,1)-mean(DATA(:,1)); s=x; [c1,d1]=general2nddecm(s); [d21,d22]=general2nddecm(c1); [d23,d24]=general2nddecm(d1); [d31,d32]=general2nddecm(d21); [d33,d34]=general2nddecm(d22); [d35,d36]=general2nddecm(d23); [d37,d38]=general2nddecm(d24); %% 三层小波包信号重构 d310(1:length(d31))=0; d320(1:length(d32))=0; d330(1:length(d33))=0; d340(1:length(d34))=0; d350(1:length(d35))=0; d360(1:length(d36))=0; d370(1:length(d37))=0; d380(1:length(d38))=0; d210(1:length(d21))=0; d220(1:length(d22))=0; d230(1:length(d23))=0; d240(1:length(d24))=0; c10(1:length(c1))=0; d10(1:length(d1))=0; % dr31 [data31]=general2ndre(d31,d320); [data31]=general2ndre(data31,d220); [dr31]=general2ndre(data31,d10); [maxB,maxFx]=getdbzhi(dr31); MAXB(1)=maxB; MAXFX(1)=maxFx; % dr32 [data32]=general2ndre(d310,d32); [data32]=general2ndre(data32,d220); [dr32]=general2ndre(data32,d10); [maxB,maxFx]=getdbzhi(dr32); MAXB(2)=maxB; MAXFX(2)=maxFx; % dr33 [data33]=general2ndre(d33,d340); [data33]=general2ndre(d210,data33); [dr33]=general2ndre(data33,d10); [maxB,maxFx]=getdbzhi(dr33); MAXB(3)=maxB; MAXFX(3)=maxFx; % dr34 [data34]=general2ndre(d330,d34); [data34]=general2ndre(d210,data34); [dr34]=general2ndre(data34,d10); [maxB,maxFx]=getdbzhi(dr34); MAXB(4)=maxB; MAXFX(4)=maxFx; % dr35 [data35]=general2ndre(d35,d360); [data35]=general2ndre(data35,d240); [dr35]=general2ndre(c10,data35); [maxB,maxFx]=getdbzhi(dr35); MAXB(5)=maxB; MAXFX(5)=maxFx; % dr36 [data36]=general2ndre(d350,d36); [data36]=general2ndre(data36,d240); [dr36]=general2ndre(c10,data36); [maxB,maxFx]=getdbzhi(dr36); MAXB(6)=maxB; MAXFX(6)=maxFx; % dr37 [data37]=general2ndre(d37,d380); [data37]=general2ndre(d230,data37); [dr37]=general2ndre(c10,data37); [maxB,maxFx]=getdbzhi(dr37); MAXB(7)=maxB; MAXFX(7)=maxFx; % dr38 [data38]=general2ndre(d370,d38); [data38]=general2ndre(d230,data38); [dr38]=general2ndre(c10,data38); [maxB,maxFx]=getdbzhi(dr38); MAXB(8)=maxB; MAXFX(8)=maxFx; maxB=max(MAXB); maxFx=mean(MAXFX); BB(i)=maxB; CC(i)=maxFx; BB=BB'; CC=CC'; end