EWT 经验小波分解

经验小波分解，通过ewt变换分解重构信号，实现对信号的降噪

clear all
clc
fs=5;
Ts=1/fs;
h=1/fs;   %时间步长
t=0:Ts:4095*Ts;
D=0.31;
a=1;
b=1;
s=0.3*sin（2*pi*0.01*t）;
x1=s+sqrt（2*D）*randn（size（t））;


y=fft（x14096）; 
pyy=y.*conj（y）/4096;
ff=fs*（0:2048）/4096;


figure（1）
subplot（221）
plot（tx1‘k‘）
title（‘Input signal sn‘）
xlabel（‘Time‘）
ylabel（‘Amplitude‘）
subplot（222）
plot（ffpyy（1:2049）‘k‘）
xlabel（‘Frequency ‘）
ylabel（‘Amplitude‘）
title（‘Spectrum of input signal sn‘）
axis（[-0.1 0.1 0 1000]）

x=sr（abhx1）;

y=fft（x4096）; 
py=y.*conj（y）/4096;
ff=fs*（0:2048）/4096;


subplot（223）
plot（tx‘k‘）
title（‘Output signal x‘）
xlabel（‘Time‘）
ylabel（‘Amplitude‘）
subplot（224）
plot（ffpy（1:2049）‘k‘）
xlabel（‘Frequency‘）
ylabel（‘Amplitude‘）
title（‘Spectrum of output signal x‘）
axis（[-0.1 0.1 0 1000]）

% [t

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

     文件       2184  2016-04-28 14:31  EWT\1D\disp_hhs2.asv

     文件       2220  2016-04-28 14:31  EWT\1D\disp_hhs2.m

     文件       2246  2016-04-28 14:31  EWT\1D\EWT1D.m

     文件        162  2016-04-28 14:31  EWT\1D\EWT_beta.m

     文件       1474  2016-04-28 14:31  EWT\1D\EWT_Meyer_FilterBank.m

     文件        920  2016-04-28 14:31  EWT\1D\EWT_Meyer_Scaling.m

     文件       1325  2016-04-28 14:31  EWT\1D\EWT_Meyer_Wavelet.m

     文件       1318  2016-04-28 14:31  EWT\1D\Hilbert_EWT.m

     文件        854  2016-04-28 14:31  EWT\1D\iEWT1D.m

     文件        708  2016-04-28 14:31  EWT\1D\Modes_EWT1D.m

     文件       1346  2016-05-10 17:04  EWT\1D\Show_EWT.m

     文件        977  2016-04-28 14:31  EWT\1D\Show_EWT_Boundaries.m

     文件       1442  2016-04-28 14:31  EWT\2D\Curvelet\AnglesLocalMax.m

     文件       2182  2016-04-28 14:31  EWT\2D\Curvelet\AnglesLocalMaxMin.m

     文件       9804  2016-04-28 14:31  EWT\2D\Curvelet\Angular_sector.m

     文件       1270  2016-04-28 14:31  EWT\2D\Curvelet\CreateAngleGrid.m

     文件       4892  2016-04-28 14:31  EWT\2D\Curvelet\EWT2D_Curvelet.m

     文件       8118  2016-04-28 14:31  EWT\2D\Curvelet\EWT2D_Curvelet_FilterBank.m

     文件        971  2016-04-28 14:31  EWT\2D\Curvelet\EWT2D_Curvelet_Scaling.m

     文件       3819  2016-04-28 14:31  EWT\2D\Curvelet\EWT_Angles_Detect.m

     文件        856  2016-04-28 14:31  EWT\2D\Curvelet\iEWT2D_Curvelet.m

     文件        736  2016-04-28 14:31  EWT\2D\Curvelet\ShowCurveletFilters.m

     文件       4532  2016-04-28 14:31  EWT\2D\Curvelet\Show_Curvelets_boundaries.m

     文件        713  2016-04-28 14:31  EWT\2D\Curvelet\Show_EWT2D_Curvelet.m

     文件        338  2016-04-28 14:31  EWT\2D\drawArcEllipse.m

     文件        179  2016-04-28 14:31  EWT\2D\drawCircle.m

     文件        182  2016-04-28 14:31  EWT\2D\drawEllipse.m

     文件       2323  2016-04-28 14:31  EWT\2D\Littlewood-Paley\EWT2D_LittlewoodPaley.m

     文件       1529  2016-04-28 14:31  EWT\2D\Littlewood-Paley\EWT2D_Meyer_FilterBank.m

     文件       1011  2016-04-28 14:31  EWT\2D\Littlewood-Paley\EWT2D_Meyer_Scaling.m

............此处省略79个文件信息