denoising

小波变换、EMD以及EMD结合小波变换几种常用的时频去噪方法的对比分析

load 97.mat;
x=X097_DE_time;
Fs=1024;
t=[1:1024]/Fs;
s=x（1:1024）;
N=1024;
n=0:N-1;
figure;
plot（ts）;
y=0.2*（randn（11024）+j*randn（11024））;
figure;
plot（ty）;
z=s+y‘;
figure;
plot（tz）;
m1=snr（sy）;
window=boxcar（length（s））;
si=periodogram（swindowNFs）;
window=boxcar（length（z））;
zi=periodogram（zwindowNFs）;
figure;
plot（si）;
figure;
plot（zi）;
r0=corrcoef（sz）;
%用db1小波对原始信号进行3层分解并提取系数
[cl]=wavedec（z2‘db3‘）;
%用默认阈值对信号进行消噪处理并图示结果
%用ddencmp函数获得信号的默认阈值
[thrsorhkeepapp]=ddencmp（‘den‘‘wv‘z）;
xd=wdencmp（‘gbl‘cl‘db3‘2thrsorhkeepapp）;
figure;
plot（txd）;
m2=snr（s（xd-s））;
window=boxcar（length（xd））;
xdi=periodogram（xdwindowNFs）;
figure;
plot（xdi）;
r1=corrcoef（sxd）;
%EMD去噪
Z=z;
[imfortnbits] = emd（Z）;
emd_visu（Ztimf）;
figure;
C1=im

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

     文件       2871  2009-12-07 19:27  新建文件夹\quzao.m

     文件       2676  2009-11-20 17:18  新建文件夹\QUZAO1.m

     文件       1183  2009-10-15 18:17  新建文件夹\rl.m

     目录          0  2010-11-09 22:16  新建文件夹

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

                 6730                    4