振动信号EMD分析matlab算法-abbr_96cb0047277bd4af6526156ae39ba1bc.rar

振动信号EMD分析matlab算法-abbr_96cb0047277bd4af6526156ae39ba1bc.rar
振动信号EMD分析matlab算法，介绍了EMD分析处理的基本原理，可在此基础上改进算法应用到具体问题上
程序如下:
clc
clear
%z=[4319.0,4358.5,4392.2,4418.9,4438.5,4453.1,4465.1,4473.0,4476.5,4477.0,4477.5,4478.8,4477.7,4472.1,4462.2,4449.9,4435.9,4419.0,4399.5,4380.5,4364.3,4351.5,4341.5,4334.7,4332.7,4335.2,4338.7,4339.5,4336.7,4332.5,4327.9,4322.4,4316.7,4312.3,4310.9,4313.0,4317.4,4324.3,4335.3,4351.0,4369.4,4386.5,4398.6,4405.0,4408.6,4411.6,4415.5,4421.3,4430.1,4441.2,4450.6,4456.0,4458.0,4460.1,4464.2,4467.7,4468.2,4465.7,4461.9,4458.3,4454.2,4448.6,4442.0,4434.9,4427.0,4417.2,4405.9,4396.0,4388.8,4383.3,4377.6,4370.3,4363.2,4357.3,4350.8,4341.8,4330.9,4321.2,4314.9,4311.8,4311.2,4313.7,4320.7,4330.7,4340.5,4347.3,4351.9,4354.9,4355.6,4352.8,4346.9,4341.4,4338.9,4338.9,4340.1,4341.7,4345.0,4350.3,4355.0,4357.3,4356.6,4355.0,4353.5,4351.0,4348.0,4345.9,4347.3,4352.1,4357.5,4361.9,4365.6,4370.0,4374.3,4377.3,4378.8,4380.3,4383.3,4386.1,4387.3,4386.9,4386.8,4388.6,4390.8,4392.2,4391.6,4389.7,4387.0,4382.0,4375.8,4370.6,4368.2,4369.2,4371.0,4371.0,4369.4,4367.2,4365.8,4364.9,4365.2,4368.3,4374.3,4380.6,4383.5,4381.1,4375.9,4371.0,4366.8,4364.1,4363.6,4366.5,4371.9,4376.6,4379.0,4380.4,4382.7,4385.9,4387.8,4387.2,4386.0,4385.8,4387.4,4389.7,4392.4,4396.9,4403.7,4410.8,4416.1,4419.2,4421.6,4424.7,4427.2,4428.5,4429.3,4431.3,4435.1,4437.9,4438.4,4437.0,4435.2,4433.2,4429.2,4422.9,4416.7,4412.6,4410.5,4408.5,4405.4,4402.1,4399.1,4394.5,4385.7,4372.2,4354.5,4332.6,4303.7];
%z=[5210.1,5291.7,5362.0,5424.2,5481.8,5534.9,5582.4,5624.1,5662.4,5698.9,5732.0,5759.4,5779.6,5793.0,5800.1,5799.7,5792.1,5780.1,5766.4,5751.9,5734.5,5712.3,5685.5,5653.8,5614.0,5562.3,5496.4,5418.3,5331.7,5238.9];
%z=[5585.5,5650.5,5699.2,5737.4,5765.4,5782.2,5787.8,5785.3,5778.3,5767.8,5755.5,5743.7,5735.3,5731.1,5728.3,5723.9,5716.2,5703.8,5684.7,5655.8,5616.0,5566.8,5509.7,5444.9,5372.3,5294.4,5214.8,5136.8,5061.1,4985.6,4907.9,4828.8,4749.6,4671.7,4597.6,4530.5,4473.1,4426.7,4388.7

%EMD  computes Empirical Mode Decomposition
%
%
%   Syntax
%
%
% IMF = EMD（X）
% IMF = EMD（X...‘Option_name‘Option_value...）
% IMF = EMD（XOPTS）
% [IMFORTNB_ITERATIONS] = EMD（...）
%
%
%   Description
%
%
% IMF = EMD（X） where X is a real vector computes the Empirical Mode
% Decomposition [1] 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 [2]:
%
%   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 [3] 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 [2]:
%
%   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（XOPTS） 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 
%
% [IMFORTNB_ITERATIONS] = 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 [THRESHOLDTHRESHOLD2TOLERANCE]
% if the input vector‘s length is less than 3 only the first parameters are
% set the remaining ones taking default values.
% default: [0.050.50.05]
%
% FIX （int）: disable the default stopping criterion and do exactly  
% number of sifting iterations for each mode
%
% FIX_H （int）: disable the default stopping criterion and do  sifting 
% iterations with |#zeros-#extrema|<=1 to stop [4]
%
%  bivariate/complex EMD options:
%
% COMPLEX_VERSION: selects the algorithm used for complex EMD （[3]）
% 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 [3]） 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

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

     文件      22273  2009-04-21 10:47  EMD 可运行\emd.m

     文件       3593  2009-04-21 10:51  EMD 可运行\fang02.m

     目录          0  2009-04-21 10:47  EMD 可运行

     文件        183  2009-02-02 09:47  Matlab中文论坛--助努力的人完成毕业设计.url

     文件       3543  2009-02-01 16:19  使用帮助：新手必看.htm

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

                29592                    5