可以使用的EWT经验小波变换

如果要运行所有功能，则需要在计算机上正确安装以下Matlab工具箱： - Flandrin的EMD工具箱（在一维变换中需要执行希尔伯特变换并且可视化时频平面） 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱（2D变换除了基于张量变换之外） 可从http://www.cs.technion.ac.il/~elad/software/获得 如果要运行所有功能，则需要在计算机上正确安装以下Matlab工具箱： - Flandrin的EMD工具箱（在一维变换中需要执行希尔伯特变换并且可视化时频平面） 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱（2D变换除了基于张量变换之外） 可从http://www.cs.technion.ac.il/~elad/software/获得 如果要运行所有功能，则需要在计算机上正确安装以下Matlab工具箱： - Flandrin的EMD工具箱（在一维变换中需要执行希尔伯特变换并且可视化时频平面） 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱（2D变换除了基于张量变换之外） 可从http://www.cs.technion.ac.il/~elad/software/获得 如果要运行所有功能，则需要在计算机上正确安装以下Matlab工具箱： - Flandrin的EMD工具箱（在一维变换中需要执行希尔伯特变换并且可视化时频平面） 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱（2D变换除了基于张量变换之外） 可从http://www.cs.technion.ac.il/~elad/software/获得 如果要运行所有功能，则需要在计算机上正确安装以下Matlab工具箱： - Flandrin的EMD工具箱（在一维变换中需要执行希尔伯特变换并且可视化时频平面） 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱（2D变换除了基于张量变换之外） 可从http://www.cs.technion.ac.il/~elad/software/获得 这个工具箱组织如下： EWT ?| ?| -1D：1D EWT功能 ?| -2D：2D EWT功能 ?| | - 小波：经验曲线变换 ?| | -Littlewood-Paley：经验的Littlewood-Paley小波变换 ?| | -Ridgelet：经验Ridgelet变换 ?| | - 张量：经验张量小波变换 ? - 边界：用于执行傅里叶支持的函数 ?| | - LocalMaxima：根据当地最大值，中途或当地最小值执行检测的功能 ?| | -MorphoMath：执行形态学操作符对谱进行预处理的功能 ?| | -PowerLaw：通过去除其幂律近似来预处理谱 ?| | -ScaleSpace：基于尺度空间方法执行检测的函数 ?| - 文档：工具箱文档 ?| -Tests ?| | -1D：对几个1D信号执行基本测试的功能 ?| | -2D：用于在不同图像上执行几个2D变换的基本测试的功能 ?| -utilities ?| | -1D：在1D情况下绘制结果的有用函数（时频平面，分量，边界） ?| | -2D：用于在2D情况下绘制结果的有用函数（不同类型的组件，2D边界，...） 如果要运行所有功能，则需要在计算机上正确安装以下Matlab工具箱： - Flandrin的EMD工具箱（在一维变换中需要执行希尔伯特变换并且可视化时频平面） 可从http://perso.ens-lyon.fr/patrick.flandrin/emd.html获得 - Elad的伪极谱FFT工具箱（2D变换除了基于张量变换之外） 可从http://www.cs.technion.ac.il/~elad/software/获得 已经包含了

function [ewtmfbboundaries]=EWT1D（fparams）

% =========================================================================
% function [ewtmfbboundaries]=EWT1D（fparams）
% 
% Perform the Empirical Wavelet Transform of f over Nscale scales. See 
% also the documentation of EWT_Boundaries_Detect for more details about
% the available methods and their parameters.
%
% Inputs:
%   -f: the input signal
%   -params: structure containing the following parameters:
%       -params.log: 0 or 1 to indicate if we want to work with
%                    the log spectrum
%       -params.preproc: ‘none‘‘plaw‘‘poly‘‘morpho‘tophat‘
%       -params.method: ‘locmax‘‘locmaxmin‘‘locmaxminf‘‘adaptive‘
%                       ‘adaptivereg‘‘scalespace‘
%       -params.reg: ‘none‘‘gaussian‘‘average‘‘closing‘
%       -params.lengthFilter: width of the above filters
%       -params.sigmaFilter: standard deviation of the above Gaussian
%                            filter
%       -params.N: maximum number of supports
%       -params.degree: degree of the polynomial （needed for the
%                       polynomial approximation preprocessing）
%       -params.completion: 0 or 1 to indicate if we try to complete
%                           or not the number of modes if the detection
%                           find a lower number of mode than params.N
%       -params.InitBounds: vector of initial bounds （in index domain）
%                           needed for the adaptive and adaptivereg methods
%				-params.typeDetect: （for scalespace method only） ‘otsu‘
%                           ‘halfnormal‘‘empiricallaw‘‘mean‘‘kmeans‘
%
% Outputs:
%   -ewt: cell containing first the low frequency component and
%         then the successives frequency subbands
%   -mfb: cell containing the filter bank （in the Fourier domain）
%   -boundaries: vector containing the set of boundaries corresponding
%                to the Fourier line segmentation （normalized between
%                0 and Pi）
%
% Author: Jerome Gilles
% Institution: UCLA - Department of Mathematics
% Year: 2013
% Version: 2.0
% =========================================================================

%% Boundary detection
% We compute the Fourier transform of f
ff=fft（f）;
% We extract the boundaries of Fourier segments
boundaries = EWT_Boundaries_Detect（abs（ff（1:round（length（ff）/2）））params）;
boundaries = boundaries*pi/round（length（ff）/2）;

%% Filtering
% We extend the signal by miroring to deal with the boundaries
l=round（length（f）/2）;
f=[f（l-1:-1:1）;f;f（end:-1:end-l+1）];
ff=fft（f）;

% We build the corresponding filter bank
mfb=EWT_Meyer_FilterBank（boundarieslength（ff））;

% We filter the signal to extract each subband
ewt=cell（length（mfb）1）;
for k=1:length（mfb）
    ewt{k}=real（ifft（conj（mfb{k}）.*ff））;
    ewt{k}=ewt{k}（l:end-l）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        3542  2017-12-13 20:09  说明.txt
     目录           0  2017-12-13 16:43  EWT20161130\
     目录           0  2017-12-13 16:43  EWT20161130\EWT\
     目录           0  2017-12-13 16:43  EWT20161130\EWT\1D\
     文件        2825  2016-12-01 01:53  EWT20161130\EWT\1D\EWT1D.m
     文件        1732  2016-12-01 01:53  EWT20161130\EWT\1D\EWT_InstantaneousComponents.m
     文件        1474  2016-12-01 01:53  EWT20161130\EWT\1D\EWT_Meyer_FilterBank.m
     文件         920  2016-12-01 01:53  EWT20161130\EWT\1D\EWT_Meyer_Scaling.m
     文件        1325  2016-12-01 01:53  EWT20161130\EWT\1D\EWT_Meyer_Wavelet.m
     文件          93  2016-12-01 01:53  EWT20161130\EWT\1D\EWT_Single_filter.m
     文件         162  2016-12-01 01:53  EWT20161130\EWT\1D\EWT_beta.m
     文件        1628  2016-12-01 01:53  EWT20161130\EWT\1D\IFcleaning.m
     文件         749  2016-12-01 01:53  EWT20161130\EWT\1D\Modes_EWT1D.m
     文件         841  2016-12-01 01:53  EWT20161130\EWT\1D\iEWT1D.m
     目录           0  2017-12-13 16:43  EWT20161130\EWT\2D\
     目录           0  2017-12-13 16:43  EWT20161130\EWT\2D\Curvelet\
     文件        1442  2016-12-01 01:53  EWT20161130\EWT\2D\Curvelet\AnglesLocalMax.m
     文件        2182  2016-12-01 01:53  EWT20161130\EWT\2D\Curvelet\AnglesLocalMaxMin.m
     文件        8620  2016-12-01 01:53  EWT20161130\EWT\2D\Curvelet\Angular_sector.m
     文件        1288  2016-12-01 01:53  EWT20161130\EWT\2D\Curvelet\CreateAngleGrid.m
     文件        8376  2016-12-01 01:53  EWT20161130\EWT\2D\Curvelet\EWT2D_Curvelet.m
     文件        7085  2016-12-01 01:53  EWT20161130\EWT\2D\Curvelet\EWT2D_Curvelet_FilterBank.m
     文件         971  2016-12-01 01:53  EWT20161130\EWT\2D\Curvelet\EWT2D_Curvelet_Scaling.m
     文件        2674  2016-12-01 01:53  EWT20161130\EWT\2D\Curvelet\EWT_Angles_Detect.m
     文件         856  2016-12-01 01:53  EWT20161130\EWT\2D\Curvelet\iEWT2D_Curvelet.m
     目录           0  2017-12-13 16:43  EWT20161130\EWT\2D\Littlewood-Paley\
     文件        3226  2016-12-01 01:53  EWT20161130\EWT\2D\Littlewood-Paley\EWT2D_LittlewoodPaley.m
     文件        1529  2016-12-01 01:53  EWT20161130\EWT\2D\Littlewood-Paley\EWT2D_Meyer_FilterBank.m
     文件        1011  2016-12-01 01:53  EWT20161130\EWT\2D\Littlewood-Paley\EWT2D_Meyer_Scaling.m
     文件        1410  2016-12-01 01:53  EWT20161130\EWT\2D\Littlewood-Paley\EWT2D_Meyer_Wavelet.m
     文件        1175  2016-12-01 01:53  EWT20161130\EWT\2D\Littlewood-Paley\EWT2D_UP_Meyer_Wavelet.m
............此处省略209个文件信息