advanced-harmonic-wavelet-packet 自己编的改进谐波小波包实现程序

自己编的改进谐波小波包实现程序，实现了信号经过谐波小波包变换后幅值不变，适合于机械故障特征提取

clear all;
close all;%得第一层信号时用的滤波器不插零，得第二层信号时用的滤波器插1个零，得第三层信号时用的滤波器插3个零。
N=1024;n=1:N;fs=1024;t=n*（1/fs）;M=2;J=5000;fh=fs/2;M1=2^M;detf=fs/N;
%x=sin（2*pi*15*t）+sin（2*pi*20*t）+sin（2*pi*75*t）+sin（2*pi*90*t）+sin（2*pi*140*t）+sin（2*pi*180*t）+sin（2*pi*200*t）+sin（2*pi*270*t）+sin（2*pi*311*t）+sin（2*pi*325*t）;
   % +sin（2*pi*375*t）;%fs=800，M=1不能重构200hz信号.fs=800，M=2不能重构50hz信号
%shuju212=load（‘F:\Bearing Data\Bearing Data\121.mat‘）;
%x=randn（1N）;
%shuju212=load（‘F:\Bearing Data\Bearing Data\212.mat‘）;
%x=shuju212.X212_DE_time（1:N）‘;%x=awgn（x1）;
%x=2*cos（2*pi*160*t）.*[zeros（1350） ones（140） zeros（1N-390）]+2*cos（2*pi*80*t）+5*cos（2*pi*10*t）;
%x=x+0*randn（1N）;
%x=2*cos（2*pi*160*t）.*[zeros（1350） ones（140） zeros（1N-390）]+2*cos（2*pi*80*t）+5*cos（2*pi*10*t）;
x=（1+0.3*cos（10*pi*t））.*sin（2*pi*

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5612  2012-09-17 13:09  advanced harmonic wavelet packet.m

----------- ---------  ---------- -----  ----

                 5612                    1