请问在emd分解过程提到的间歇检测准则怎么实现?
网上下载的emd程序好像对这个没有介绍,谁做过这方面的工作吗? 本帖最后由 wdhd 于 2016-9-10 14:34 编辑原帖由 tangaoming 于 2007-4-28 20:31 发表
网上下载的emd程序好像对这个没有介绍,谁做过这方面的工作吗?
请仔细看最新版的emd.m中帮助信息 eight,你说的如下吗?我看了看,它只介绍了mask,但好像不是间歇检测准则,你说的那部分介绍了
Description
%
%
% IMF = EMD(X) where X is a real vector computes the Empirical Mode
% Decomposition of X, resulting in a matrix IMF containing 1 IMF per row, the
% last one being the residue. The default stopping criterion is the one proposed
% in :
%
% at each point, mean_amplitude < THRESHOLD2*envelope_amplitude
% &
% mean of boolean array {(mean_amplitude)/(envelope_amplitude) > THRESHOLD} < TOLERANCE
% &
% |#zeros-#extrema|<=1
%
% where mean_amplitude = abs(envelope_max+envelope_min)/2
% and envelope_amplitude = abs(envelope_max-envelope_min)/2
%
% IMF = EMD(X) where X is a complex vector computes Bivariate Empirical Mode
% Decomposition of X, resulting in a matrix IMF containing 1 IMF per row, the
% last one being the residue. The default stopping criterion is similar to the
% one proposed in :
%
% at each point, mean_amplitude < THRESHOLD2*envelope_amplitude
% &
% mean of boolean array {(mean_amplitude)/(envelope_amplitude) > THRESHOLD} < TOLERANCE
%
% where mean_amplitude and envelope_amplitude have definitions similar to the
% real case
%
% IMF = EMD(X,...,'Option_name',Option_value,...) sets options Option_name to
% the specified Option_value (see Options)
%
% IMF = EMD(X,OPTS) is equivalent to the above syntax provided OPTS is a struct
% object with field names corresponding to option names and field values being the
% associated values
%
% = EMD(...) returns an index of orthogonality
% ________
% _|IMF(i,:).*IMF(j,:)|
% ORT = \ _____________________
% /
% ? || X ||?
% i~=j
%
% and the number of iterations to extract each mode in NB_ITERATIONS
%
%
% Options
%
%
%stopping criterion options:
%
% STOP: vector of stopping parameters
% if the input vector's length is less than 3, only the first parameters are
% set, the remaining ones taking default values.
% default:
%
% FIX (int): disable the default stopping criterion and do exactly <FIX>
% number of sifting iterations for each mode
%
% FIX_H (int): disable the default stopping criterion and do <FIX_H> sifting
% iterations with |#zeros-#extrema|<=1 to stop
%
%bivariate/complex EMD options:
%
% COMPLEX_VERSION: selects the algorithm used for complex EMD ()
% COMPLEX_VERSION = 1: "algorithm 1"
% COMPLEX_VERSION = 2: "algorithm 2" (default)
%
% NDIRS: number of directions in which envelopes are computed (default 4)
% rem: the actual number of directions (according to ) is 2*NDIRS
%
%other options:
%
% T: sampling times (line vector) (default: 1:length(x))
%
% MAXITERATIONS: maximum number of sifting iterations for the computation of each
% mode (default: 2000)
%
% MAXMODES: maximum number of imfs extracted (default: Inf)
%
% DISPLAY: if equals to 1 shows sifting steps with pause
% if equals to 2 shows sifting steps without pause (movie style)
% rem: display is disabled when the input is complex
%
% INTERP: interpolation scheme: 'linear', 'cubic', 'pchip' or 'spline' (default)
% see interp1 documentation for details
%
% MASK: masking signal used to improve the decomposition according to
%
%
% Examples
%
%
%X = rand(1,512);
%
%IMF = emd(X);
%
%IMF = emd(X,'STOP',,'MAXITERATIONS',100);
%
%T=linspace(0,20,1e3);
%X = 2*exp(i*T)+exp(3*i*T)+.5*T;
%IMF = emd(X,'T',T);
%
%OPTIONS.DISLPAY = 1;
%OPTIONS.FIX = 10;
%OPTIONS.MAXMODES = 3;
% = emd(X,OPTIONS); 本帖最后由 wdhd 于 2016-9-10 14:35 编辑
原帖由 tangaoming 于 2007-4-28 21:35 发表
eight,你说的如下吗?我看了看,它只介绍了mask,但好像不是间歇检测准则,你说的那部分介绍了
Description
%
%
% IMF = EMD(X) where X is a real vector computes the Empirical Mode
% Decompositi ...
mask就是用来处理间歇检测信号(即模式混叠信号)的。另外,本版上有不少这个问题的讨论,请搜索一下
可不可以不用mask信号?用截止频率实现?就是间歇检测准则?
本帖最后由 wdhd 于 2016-9-10 14:35 编辑原帖由 eight 于 2007-4-28 22:34 发表
mask就是用来处理间歇检测信号(即模式混叠信号)的。另外,本版上有不少这个问题的讨论,请搜索一下 [/quo
我不想用mask,只是想辨识某一频段的特征。能否用截止频率(间歇检测准则实现)?
截止频率实现
具体的处理方法可能和所采用的程序有关。是采用将截止频率折算为对应的周期内的峰值点数的方法来实现间歇频率的,点数不是整数时,按照就近的原则取整。
页:
[1]