EWT经验小波matlab应用实现信号处理

完整的EWT工具箱 一维test测试m文件 具有二维工具箱 本测试主要应用一维信号处理

load（‘x.mat‘）;

params.SamplingRate =-1;%put -1 if you don‘t know the sampling rate
%不知道采样频率的情况下取-1
params.log=0;%Perform the detection on the log spectrum instead the spectrum
%是否利用对数谱

params.globtrend = ‘none‘;%Some preprocessing are available in order to remove a global trend of the spectrum
%‘none‘‘plaw‘‘poly‘‘morpho‘tophat‘
%为了消除频谱的全局趋势，可以进行一些预处理
params.degree=0; %degree for the polynomial interpolation  ‘poly‘预处理方法的degree

params.reg = ‘none‘;%Moreover regularized version of the spectrum can be obtained by the following methods
%‘none‘‘gaussian‘‘average‘‘closing‘
%频谱的正则化方法
params.lengthFilter = 0;%width of the above filters
params.sigmaFilter = 0;%Gaussian filter standard deviation

params.detect = ‘locmax‘;%This function segments f into a certain amount of supports by using different technics
%‘locmax‘‘locmaxmin‘‘locmaxminf‘‘adaptive‘ ‘adaptivereg‘‘scalespace‘ 
%把信号的频谱分成特定数量边界检测方法
params.N = 5;%maximum number of supports （needed for thelocmax and locmaxmin methods） 
%仅在‘locmax‘‘locmaxmin‘‘locmaxminf‘方法下使用
params.completion = 0; %choose if you want to force to have params.N modes in case the algorithm found less ones （0 or 1）
%边界检测数量等于0或1时，是否重置数量为N
params.typeDetect=‘otsu‘; %（for scalespace only）:‘otsu‘‘halfnormal‘‘empiricallaw‘‘mean‘‘kmeans‘
params.InitBounds = [0 0 0];%vector of initial bounds （in index domain） needed for the （adaptive and adaptivereg methods）
%（索引域）adaptive和adaptivereg方法的初始边界向量

% Choose the results you want to display 选择显示结果
Bound=1;   % Display the detected boundaries on the spectrum 显示频谱检测边界
Comp=1;    % Display the EWT components 显示信号的EWT分解部分
Rec=1;     % Display the reconstructed signal 显示原始信号和重建信号
TFplane=1;% Display the time-frequency plane （by using the Hilbert transform）. 
%You can decrease the frequency resolution by changing the subresf variable below.
subresf=1;

%We perform the empirical transform and its inverse compute the EWT 
%（and get the corresponding filter bank and list of boundaries）
[ewtmfbboundaries]=EWT1D（xparams）;

%Show the results

if Bound==1 %Show the boundaries on the spectrum
    div=1;
    if （strcmp（params.detect‘adaptive‘）||strcmp（params.detect‘adaptivereg‘））
        Show_EWT_Boundaries（abs（fft（x））boundariesdivparams.SamplingRateInitBounds）;
    else
        Show_EWT_Boundaries（abs（fft（x））boundariesdivparams.SamplingRate）;
    end
end

if Comp==1 %Show the EWT components and the reconstructed signal
    if Rec==1
        %compute the reconstruction
        rec=iEWT1D（ewtmfb）;
        Show_EWT（ewtxrec）;
    else
        Show_EWT（ewt）;
    end    
end

if TFplane==1 %Show the time-frequency plane by using the Hilbert transform
    TF=EWT_TF_Plan（ewtboundariesparams.SamplingRatex[][]subresf[]）;
end

%边际谱的绘制marginal specturm
fs=1;
mar_spec = mean（TF 2）;
ff = （1: size（TF 1）） / size（TF 1） * fs/2;
figure（3）
plot（ff mar_spec）
xlabel（‘频率/Hz‘）;
ylabel（‘信号幅值‘）;
title（‘边际谱‘）;

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

     文件       2825  2016-12-01 01:53  EWTtests\1D\EWT1D.m

     文件        162  2016-12-01 01:53  EWTtests\1D\EWT_beta.m

     文件       1730  2018-03-29 18:01  EWTtests\1D\EWT_InstantaneousComponents.m

     文件       1474  2016-12-01 01:53  EWTtests\1D\EWT_Meyer_FilterBank.m

     文件        920  2016-12-01 01:53  EWTtests\1D\EWT_Meyer_Scaling.m

     文件       1325  2016-12-01 01:53  EWTtests\1D\EWT_Meyer_Wavelet.m

     文件         93  2016-12-01 01:53  EWTtests\1D\EWT_Single_filter.m

     文件        841  2016-12-01 01:53  EWTtests\1D\iEWT1D.m

     文件       1628  2018-03-29 17:09  EWTtests\1D\IFcleaning.m

     文件        749  2016-12-01 01:53  EWTtests\1D\Modes_EWT1D.m

     文件       3365  2018-04-11 21:45  EWTtests\1D_Utilities\EWT_TF_Plan.m

     文件       1158  2018-04-02 16:53  EWTtests\1D_Utilities\Show_EWT.m

     文件       2432  2016-12-01 01:53  EWTtests\1D_Utilities\Show_EWT_Boundaries.m

     文件       1629  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\Adaptive_Bounds_Adapt.m

     文件       2009  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\EpsNeighLocalMaxMin.m

     文件       1223  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\LocalMaxima\LocalMax.m

     文件       2113  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\LocalMaxima\LocalMaxMin.m

     文件       2653  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\LocalMaxima\LocalMaxMin2.m

     文件       1020  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\ScaleSpace\EmpiricalLaw.m

     文件        881  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\ScaleSpace\GSS_BoundariesDetect.m

     文件       1080  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\ScaleSpace\HalfNormalLaw.m

     文件       1215  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\ScaleSpace\kmeansDetect.m

     文件       3406  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\ScaleSpace\LengthScaleCurve.m

     文件        407  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\ScaleSpace\maxcheckplateau.m

     文件       1280  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\ScaleSpace\MeaningfulScaleSpace.m

     文件        867  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\ScaleSpace\MeanTh.m

     文件       1270  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\ScaleSpace\OtsuMethod.m

     文件       1286  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\ScaleSpace\PlanGaussianScaleSpace.m

     文件       4037  2016-12-01 01:53  EWTtests\Boundaries\BoundariesDetectMethod\ScaleSpace\RemoveMerge.m

     文件        835  2016-12-01 01:53  EWTtests\Boundaries\EWT_Boundaries_Completion.m

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