[16]刘崇春，裘正定，杜锡钰. 小波变换理论及其在信号处理中的应用[J]. 北方交通大学
学报, 1997,30(5):67～86.
[17]Gabor.D. Theory of communication.[J].IEEE,1946,93(3):429～457. [18]Wexler J. Raz S. Discrete Gabor expansion. Signal Processing[J]. 1990,21(3): 20 7～220.
[19]询撇利，何昭水，高鹰.信号处理的自适应理论.北京:科学出版社，2田6:130页 [20]吴湘淇.信号、系统与信号处理.北京:电子工业出版礼1996:78页
[21]T.C.M.Classen，W.F.G.Mecklenbrauker.The Wigner Distribution Journal of Rese arch，1980，35:276一303P
[22]T.C.M.Classen，W.F.G.Mecklenbrauker.The Wigner Distribution Journal of Rese arch，1980，35:375一390P
[23]杨福生.小波变换的工程分析与应用.北京:科学出版社，1999:78页
[24]马孝江,王凤利，蔡悦等.局域波时频分布在转子系统早期故障诊断中的应用研究.中国 电机工程，24(3):161一164页
[25]廖庆斌，李舜酩.一种旋转机械振动信号特征提取的新方法.中国机械工程.2006:167一 168页
[26] 何永勇,褚福磊,王庆禹.小波分析在应用中的两个问题研究.振动工程学报.2002.15 (2):228-232
[27] 赵建林, 代月明, 张冬梅, 朱习军.基于小波包的机械故障信号消噪处理2005.9 : 141-142
[28]杨立华，戴道清，等译.信号处理的小波导引.北京:机械工业出版社，9:159一188页 [29]黄晶，阂沛文.小波分析在管道缺陷超声检测中的应用.传感技术学报[J]，2003，(3):263-266
[30]邹晓红.小波分析在超声回波测量中的应用.传感器技术[J]，2005，24(7):78-81 [31]杜文正，黄先祥，谢建等.小波奇异性检测在超声信号处理中的应用.NDT无损检测[J]，2005， 26(6):319-321
23
附件： (1)rigrsure
y=Sheet1;%载入波形数据 subplot(2,1,1);plot(y);%原始信号
title('原始信号');xlabel('样本');ylabel('信号');
[C,L]=wavedec(y,3,'db6');%用小波函数db6对信号进行3层分解 thr=thselect(y,'rigrsure');%用rigrsure选择阈值 keepapp=1;%低频系数不进行阈值量化
xd=wdencmp('gbl',C,L,'db6',3,thr,'s',keepapp);%软阈值去噪 subplot(2,1,2);plot(xd);%去噪后的信号
title('去噪后的信号');xlabel('样本');ylabel('信号'); (2) heursure
y=Sheet1;%载入波形数据 subplot(2,1,1);plot(y);%原始信号
title('原始信号');xlabel('样本');ylabel('信号');
[C,L]=wavedec(y,3,'db6');%用小波函数db6对信号进行3层分解
24
thr=thselect(y,'heursure');%用heursure选择阈值 keepapp=1;%低频系数不进行阈值量化
xd=wdencmp('gbl',C,L,'db6',3,thr,'s',keepapp);%软阈值去噪 subplot(2,1,2);plot(xd);%去噪后的信号
title('去噪后的信号');xlabel('样本');ylabel('信号'); (3) sqtwolog
y=Sheet1;%载入波形数据 subplot(2,1,1);plot(y);%原始信号
title('原始信号');xlabel('样本');ylabel('信号');
[C,L]=wavedec(y,3,'db6');%用小波函数db6对信号进行3层分解 thr=thselect(y,'sqtwolog');%用sqtwolog选择阈值 keepapp=1;%低频系数不进行阈值量化
xd=wdencmp('gbl',C,L,'db6',3,thr,'s',keepapp);%软阈值去噪 subplot(2,1,2);plot(xd);%去噪后的信号
title('去噪后的信号');xlabel('样本');ylabel('信号'); (4) minimaxi
y=Sheet1;%载入波形数据 subplot(2,1,1);plot(y);%原始信号
title('原始信号');xlabel('样本');ylabel('信号');
[C,L]=wavedec(y,3,'db6');%用小波函数db6对信号进行3层分解 thr=thselect(y,'minimaxi');%用minimaxi选择阈值 keepapp=1;%低频系数不进行阈值量化
xd=wdencmp('gbl',C,L,'db6',3,thr,'s',keepapp);%软阈值去噪 subplot(2,1,2);plot(xd);%去噪后的信号
title('去噪后的信号');xlabel('样本');ylabel('信号');
25