matlab_小波分解_希尔伯特解调_快速傅里叶变换

clc clear all SampleFrequency=20000;%采样频率 N=20480;%数据点数 data=load('C:\Program Files\MATLAB71\work\1st_test\2003.11.25.23.13.21');%读数据 fs=SampleFrequency; x=data(:,5); x=x-mean(x); %画时域波形 t=(1:N)/fs; figure('name','时域波形'); plot(t,x);%时域波形 xlabel('时间'); ylabel('幅值'); % [Fx,Fy] = ZOOMFFT(x,fs,25,40); % figure; % plot(Fx,Fy); % xlabel('频率'); % ylabel('幅值'); % 线性相关分析1：直接调用相关计算的函数 % y=xcorr(x); % y=y(1:N); % y=flipud(y); % for k=1:N % y(k)=y(k)/(N-k+1); % end % t=(1:N)/fs; % figure; % plot(t, y); % xlabel('time/s'); % ylabel('amp/mv'); % title('相关函数：调用函数计算'); % y=fft(y,N)/N; % y=2*abs(y(1:N/2)); % fx=(1:N/2)/N*fs; % figure('name','自功率谱'); % plot(fx,y); % xlabel('频率'); % ylabel('幅值'); %画频域波形 y1=fft(x,N)/N; y=2*abs(y1(1:N/2)); fx=(1:N/2)/N*fs; figure('name','频域波形'); plot(fx,y); xlabel('频率'); ylabel('幅值'); %zoom后 h=abs(hilbert(x)); H=h-mean(h); y1=fft(H,N)/N; y=2*abs(y1(1:N/2)); fx=(1:N/2)/N*fs; figure('name','hilbert解调频域波形'); plot(fx,y); xlabel('频率'); ylabel('幅值'); % [Fx,Fy] = ZOOMFFT(H,fs,150,350); % figure('name','hilbert解后zoomFFT频域波形'); % plot(Fx,Fy); % xlabel('频率'); % ylabel('幅值'); % % %hilbert解调后画频域波形 % h=abs(hilbert(x)); % H=h-mean(h); % y1=fft(H,N)/N; % y=2*abs(y1(1:N/2)); % fx=(1:N/2)/N*fs; % figure('name','hilbert解调频域波形'); % plot(fx,y); % xlabel('频率'); % ylabel('幅值'); % % %% 3层小波包分解 % 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); % myfft(dr31); % % dr32 % [data32]=general2ndre(d310,d32); % [data32]=general2ndre(data32,d220); % [dr32]=general2ndre(data32,d10); % myfft(dr32); % % dr33 % [data33]=general2ndre(d33,d340); % [data33]=general2ndre(d210,data33); % [dr33]=general2ndre(data33,d10); % myfft(dr33); % % dr34 % [data34]=general2ndre(d330,d34); % [data34]=general2ndre(d210,data34); % [dr34]=general2ndre(data34,d10); % myfft(dr34); % % dr35 % [data35]=general2ndre(d35,d360); % [data35]=general2ndre(data35,d240); % [dr35]=general2ndre(c10,data35); % myfft(dr35); % % dr36 % [data36]=general2ndre(d350,d36); % [data36]=general2ndre(data36,d240); % [dr36]=general2ndre(c10,data36); % myfft(dr36); % % dr37 % [data37]=general2ndre(d37,d380); % [data37]=general2ndre(d230,data37); % [dr37]=general2ndre(c10,data37); % myfft(dr37); % % dr38 % [data38]=general2ndre(d370,d38); % [data38]=general2ndre(d230,data38); % [dr38]=general2ndre(c10,data38); % myfft(dr38); % figure() % subplot(4,2,1) % plot(dr31) % ylabel('d31') % subplot(4,2,3) % plot(dr32) % ylabel('d32') % subplot(4,2,5) % plot(dr33) % ylabel('d33') % subplot(4,2,7) % plot(dr34) % ylabel('d34') % subplot(4,2,2) % plot(dr35) % ylabel('d35') % subplot(4,2,4) % plot(dr36) % ylabel('d36') % subplot(4,2,6) % plot(dr37) % ylabel('d37') % subplot(4,2,8) % plot(dr38) % ylabel('d38') % % % -----------求解信息熵-------------- % % c=[dr31;dr32;dr33;dr34;dr35;dr36;dr37;dr38]; % % n=10; % % m=size(c,2)/n; % % for i=1:size(c,1) % % for j=1:n % % xx=c(i,(m*(j-1)+1):(m*j)); % % Ex(i,j)=sum(xx.^2); % % end % % Exn(i)=sum(Ex(i,:)); % % for j=1:n % % p(j)=Ex(i,j)/Exn(i); % % end % % E(i)=-sum(p.*log(p)); % % end % % E=E'; % % Et(:,k)=E; % % end % % tsc % % fsc % % Et % % fsca=fsc/1000 % %