经验模式分解的数学理论研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
信号分析是对信号基本性质的研究和表征。从传统观点来看,信号的不同表示通常是将信号在函数空间的完备正交基集上展开。但是,由于基函数的形式固定,往往难以适应现实信号的复杂性,因此需要一种更加灵活的信号分解方法。于是,经验模式分解(EMD)方法应运而生,它不是采用某种正交基对信号进行分解,而是通过特殊的迭代规则将信号分解为一系列本征模态函数(IMF)之和,这种方法摆脱了传统信号分解方法中基函数的束缚,从而能够更加灵活的发掘信号特征。但是,这种灵活性是以迭代规则的复杂性为代价的,因此经验模式分解数学理论研究的重点在于对迭代规则性质的分析。
     经验模式分解中的迭代过程称为筛过程,其基本步骤如下:首先找到信号的极值点,然后构造信号的上包络和下包络,最后从原有信号中减去上下包络的均值。重复这个迭代过程,直到信号的上下包络均值为零,就得到了信号中包含的本征模态函数。从前面的步骤可以看出,筛过程包括两个关键步骤:(1)信号极值点的定位;(2)上下包络的构造。本文正是从这两个关键点出发,提出了一种精确的信号极值点定位方法和一种解析性质良好的包络构造方法。
     信号极值点定位算法主要建立在插值理论的基础上,分别提出了基于代数插值和三角插值的极值点定位算法。其中,代数插值采用了三次样条插值和分段Hermite插值,分别介绍了它们的插值公式的构造方法和插值误差估计的表达式。对于三角插值,不仅给出了插值公式,并首次给出了三角插值误差估计公式及其上界的推导方法。
     包络构造算法首先定义了一组信号包络必须满足的规则,然后采用二次多项式和三次多项式来构造满足规则的包络,最后将这种方法推广到n次多项式。同时,证明了这种包络的一个重要性质:采用这种包络进行迭代时,信号包络的波动会逐渐变小,最后收敛到恒包络信号。因此,采用这种包络的经验模式分解算法,是一种恒包络信号分解方法。在分析包络性质的过程中,发现这种包络通过若干关键节点,从而提出了一种基于分段Hermite插值的等效包络构造方法。
     本文针对目前经验模式分解中存在的理论空白开展了研究,解决了低采样率情况下信号中极值点定位的问题,同时设计了一种产生恒包络本征模态函数的包络构造方法。仿真结果表明,在对于恒包络组成的信号进行分解时,该方法优于现有的经验模式分解算法。因此,本文的研究工作具有较高的理论价值和一定的实际意义。
Signal analysis is used to research and represent the basic property of the signal. In traditional view, the different representations of the signal are made by the signal's extention onto different orthogonal basises. However, the basis function is fixed in an application so that it can not accommodate the special need of the practical signal. Empirical Mode Decomposition (EMD) is just to solve the problem. The signal decomposites into Intrinsic Mode Functions (IMFs) not by some orthogonal basis but by some iterative rule. The new method throws off the chains of the basis and finds the feature of the signal adaptively. However, the flexibility is at the cost of the complexity of the iterative rule. Therefore, the mathematical theory research of EMD focuses on the analysis to the property of iterative rule.
     The iterative process in EMD is called the sifting process, which has the following steps: finding the extrema of the signal, constructing the upper and lower envelopes of the signal by the extrema and subtracting the mean of the two envelopes from original signal. The steps are repeated until the mean envelope is close to zero and the residue signal is called IMF. From the steps, the two key parts are the location of the extream and the construction of the envelope. This dissertation focuses on these two points and proposes a method to locate the extrema and a theory to design the envelopes.
     The location algorithm is based on the interpolation theory including algebraic interpolation and trigonometric interpolation. Algebraic interpolation mainly composes of cubic interpolation and piecewise Hermite interpolation, and their interpolation formula and errors are introduced. For the trigonometric interpolation, the interpolation error is discussed and its the upper bound is also given.
     The envelope design defines the property of the undetermined envelope, and then constructs the envelope by quadratic polynomial, cubic polynomial and n-order polynomial. The important property of these envelopes is that the fluctuation of the envelope will be decreased in the iterative process and the iterative process will converge to a constant-envelope signal. Therefore, the EMD using the new envelope is a method that decomposites the signal into the constant-envelope signals. In the analysis of the envelopes, we find the envelopes pass through some special points and propose an equivalent construction method by piece-wise Hermite interpolation.
     The dissertation intends tofill the gaps in the theoretical research of Empirical Mode Decomposition. We solve the extrema location problem under low sampling rate and propose an envelope construction method that leads to a constant-envelope signal decomposition. Simulations show that the new algorithm is better than the traditional algorithm in the constant-envelope signal decomposition. Therefore, the works of this dissertation have strong theoretical and practical significance.
引文
[1]L.Cohen.时-频分析:理论与应用.白居宪译.西安:西安交通大学出版社,1998
    [2]张贤达,保铮.非平稳信号分析与处理.北京:国防工业出版社,1998.
    [3]王宏禹.非平稳随机信号分析与处理.北京:国防工业出版社,1999.
    [4]Patrick J.Loughlin,Leon Cohen.The uncertainty principle:Global,Local,or,Both?.IEEE Transactions Signal Processing.2004,52(5):1218-1227.
    [5]Masoud Karimi-Ghartemani,Alirera K.Ziarani.A nonlinear time-frequency analysis method.IEEE Transactions on Signal Processing.2004,52(6):1585-1595.
    [6]Antonl Homs-Corbera,Jose Antonio,et al.Time-frequency detection and analysis of wheezes during forced exhalation.IEEE Transactions on Biomedical Engineering.2004,51(1):182-186.
    [7]Cedric Fevotte,Christian Doncarli.Two contributions to blind source separation using time-frequency distributions.IEEE Signal Processing Letters.2004,11(3):386-389.
    [8]Y.T.Chars.Wavelet basics.Boston:Kluwer.1995.
    [9]N.E.Huang,Z.Shen,S.R.Long,et al.The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis.Proc.Roy.Soc.,London.A.1998,454:903-995.
    [10]N.E.Huang,S.R.Long and Z.Shen.The mechanism for frequency downshift in nonlinear wave evolution.Adv.Appl.Mech.,1996,32:59-111.
    [11]X.Zhu,Z.Shen,S.D.Eckermann et al.Gravity wave characteristics in the middle atmosphere derived from the empirical mode decomposition method.J.Geophys.Res.,1997,102(D14):16545-16561.
    [12]W.Huang,Z.Shen,N.E.Huang and Y.C.Fung.Engineering analysis of biological variables:an example of blood pressure over 1 day.Proc.Natl.Acad.Sci.USA,1998,95(9):4816-4821.
    [13]N.E.Huang.Computer implemented empirical mode decomposition method,apparatus and article of manufacture,U.S.Patent 5983162,1999.
    [14]N.E.Huang,Z.Shen and S.R.Long.A new view of nonlinear water waves:The Hilbert Spectrum.Annual Review of Fluid Mechanics,1999,31(1):417-457.
    [15]M.L.C.Wu,S.Schubert and N.E.Huang.The Development of the South Asian Summer Monsoon and the Intraseasonal Oscillation.Journal of Climate,1999,12(7):2054-2075.
    [16]W. Huang, Z. Shen, N. E. Huang et al. Nonlinear indicial response of complex nonsta-tionary oscillations as pulmonary hypertension responding to step hypoxia. Proc. Natl. Acad. Sci. USA, 1999, 96(5): 1834-1839.
    
    [17]N.E. Huang, H.H. Shih and Z Shen. The Ages of Large Amplitude Coastal Seiches on the Caribbean Coast of Puerto Rico. Journal of Physical Oceanography, 2000, 30(8): 2001-2012.
    
    [18]N.E. Huang, C.C. Chem and K. Huang. A New Spectral Representation of Earthquake Data: Hilbert Spectral Analysis of Station TCU129, Chi-Chi, Taiwan, 21 September 1999. Bulletin of the Seismological Society of America, 2001,91(5): 1310-1338.
    
    [19]N. E. Huang, M. L. C. Wu, S. R. Long et al. A confidence limit for the empirical mode decomposition and Hilbert spectral analysis. Proc. Roy. Soc, London. A., 2003, 459:2317-2345.
    
    [20] P. Gloersen and N.E. Huang. Comparison of interannual intrinsic modes in hemispheric sea ice covers and other geophysical parameters. IEEE Transactions on Geoscience and Remote Sensing, 2003,41(5): 1062-1074
    
    [21]N. E. Huang, M. L. Wu, W. Qu et al. Applications of Hilbert-Huang transform to non-stationary financial time series analysis. Applied Stochastic Models in Business and Industry, 2003, 19(3): 245-268
    
    [22] Z. Wu and N. E. Huang. A study of the characteristics of white noise using the empirical mode decomposition method. Proc. R. Soc. London A, 2004,460: 1597-1611.
    
    [23]Z. Wu and N.E. Huang. Ensemble empirical mode decomposition: A noise assisted data analysis method. Center for Ocean-Land-Atmosphere Studies Technical Report, 2005, 193:1-51.
    
    [24] N. E. Huang. Computing Instantaneous Frequency by Normalizing Hilbert Transform, U.S. Patent 6901353,2005.
    
    [25]N. E. Huang and S. Shen. Hilbert-Huang Transform and Its Applications, 1st ed. New York: World Scientific, 2005.
    
    [26] L. F. Bliven, N. E. Huang and S. R. Long. Experimental study of the influence of wind on Benjamin-Feir sideband instability. Journal of Fluid Mechanics Digital Archive, 2006, 162: 237-260
    
    [27]N. E. Huang and S. R. Long. An experimental study of the surface elevation probability distribution and statistics of wind-generated waves. Journal of Fluid Mechanics Digital Archive, 2006,101(1): 179-200
    
    [28]N. E. Huang, S. R. Long, C. C. Tung et al. A unified two-parameter wave spectral model for a general sea state. Journal of Fluid Mechanics Digital Archive, 2006, 112: 203-224
    [29]N.E.Huang.A plea for adaptive data analysis.Proceedings of SPIE,2007,6576:65760P1-6.
    [30]Z.Wu,N.E.Huang,S.R.Long et al.On the trend,detrending,and variability of nonlinear and nonstationary time series.Proceeding of the National Academy of Sciences of the USA,2007,104(38):14889-14894
    [31]N.E.Huang and Z.Wu.A review on Hilbert-Huang transform:Method and its applications to geophysical studies.Reviews of Geophysics,2008,46:1-23.
    [32]N.E.Huang and L.W.Salvino.System and method of analyzing vibrations and identifying failure signatures in the vibrations.U.S.Patent 7346461,2008.
    [33]N.E.Huang and P.It.Analyzing nonstationary financial time series via hilbert-huang transform(HHT).U.S.Patent 7464006,2008.
    [34]G.Rilling,P.Flandrin,and P.Goncalves.On empirical mode decomposition and its algorithms.The IEEE Workshop Nonlinear Signal Image Process(NSIP),San Francisco,CA,2003.
    [35][Online].Available:http://perso.ens-lyon.fr/patrick.flandrin/emd.html
    [36]Q.Gai,X.Ma,H.Zhang et al.Processing time-varying signals by a new method.CIE International Conference on Radar,2001:1011-1014
    [37]钟佑明,秦树人,汤宝平.一种振动信号新变换法的研究.振动工程学报.2002,15(2):233-238
    [38]于德介,程军圣,杨宇.机械故障诊断的HILBERT-HUANG变换方法.北京:科学出版社,2006
    [39]Z.Xu,B.Huang and S.Xu.Exact location of extrema for empirical mode decomposition.Electronics Lett.2008,44(8):551-552.
    [40]J.Wang,Y.Peng,X.Peng.Similarity searching based boundary effect processing method for empirical mode decomposition.Electronics Letters,Jan.2007,43(1):58-59
    [41]Q.H.Chen,N.E.Huang,et al.A B-spline approach for empirical mode decompositions.Advances in Computational Mathematics,2006(24):171-195.
    [42]Y.Kopsinis,S.McLaughlin.Investigation and Performance Enhancement of the Empirical Mode Decomposition Method Based on a Heuristic Search Optimization Approach.IEEE Trans.Signal Process.,2008,56(1):1-13.
    [43]R.C.Sharpley and V.Vatchev.Analysis of the intrinsic mode functions.Constr.Approx.,2006,24:17-47.
    [44]S.Kizhner,K.Blank,T.Flatley et.al.On certain theoretical developments underlying the Hilbert-Huang transform.IEEE Aerospace Conference,2006:1-14.
    [45]A.O.Boudraa and J.C.Cexus.EMD-Based Signal Filtering.IEEE Transactions on Instrumentation and Measurement,2007,56(6):2196-2202.
    [46]J.C.Nunes,Y.Bouaoune,et al.Texture analysis based on the bidimensional empirical mode decomposition with gray-level co-occurrence models.IEEE,Machine Vision and Application,2003(2):633-635.
    [47]N.Bi,Q.Sun,D.Huang et al.Robust Image Watermarking Based on Multiband Wavelets and Empirical Mode Decomposition.IEEE Transactions on Image Processing,2007,16(8):1956-1966
    [48]K.I.Molla,K.Hirose,N.Minematsu.Localization Based Separation of Mixed Audio Signals with Binary Masking of Hilbert Spectrum,Acoustics,Speech and Signal Processing,2006.ICASSP 2006:V85-V88
    [49]K.I.Molla,K.Hirose,N.Minematsu.On the empirical mode decomposition in time frequency representation of audio signals considering disjoint orthogonality.Nonlinear Signal and Image Processing,2005.NSIP 2005:44-45
    [50]R.Balocchi,D.Menicucci,E.Santarcangelo,et al.Deriving the respiratory sinus arrhythmia from the heartbeat time series using empirical mode decomposition.Chaos,Solitons and Fractals,2004,20(1):171-177.
    [51]J.C.Echeverria,J.A.Crowe,M.S.Woolfson,et al.Application of empirical mode decomposition to heart rate variability analysis.Medical and Biological Engineering and Computing,2001(39):471-479.
    [52]Z.K.Peng,P.W.Tse,F.L.Chu.An improved Hilbert Huang transform and its application in vibration signal analysis.Journal of Sound and Vibration,2005(286):187-200.
    [53]高强,杜小山,范虹,孟庆丰.Hilbert-Huang Based Approach for Structural Damage Detection.振动工程学报,2007(1):15-18.
    [54]赵知劲,尹霆等.一种基于HHT的数字调制信号识别方法.微波学报,2006(3):66-70.
    [55]毛炜,金荣洪等.基于HHT变换的时频分析法及其在2FSK系统解调中的应用.电子与信息学报,2006(12):2318-2322.
    [56]D.Cummings,et al.Travelling Waves in Dengue Hemorrhagic Fever Incidence in Thailand,Nature,2004,427:344-347.
    [57]Hualou Liang,Qiuhua Lin and J.D.Z.Chen,Application of the Empirical Mode Decomposition to the Analysis of Esophageal Manometric Data in Gastroesophageal Reflux Disease,IEEE Transactions on Biomedical Engineering,2005,52(10),1692-1701.
    [58]G.Rilling and P.Flandrin,on the Influence of Sampling on the Empirical Mode Decomposition.Acoustics,Speech and Signal Processing,2006.ICASSP 2006,3:444-447.
    [59]Y.Kopsinis,S.MeLaughlin.Investigation of the Empirical Mode Decomposition Based on Genetic Algorithm Optimization Schemes.Acoustics,Speech and Signal Processing,2006.ICASSP 2007,3:1397-1400.
    [60]G.Rilling and P.Flandrin.One or two frequencies? The empirical mode decomposition answers.IEEE Transactions on Signal Processing,2007,56(1):85-95.
    [61]许宝杰,张建民,徐小力等.抑制EMD端点效应方法的研究.北京理工大学学报.2006,26(3):196-200
    [62]胡维平,莫家玲,龚英姬等.经验模态分解中多种边界处理方法的比较研究.电子与信息学报,2007,29(6):1394-1398
    [63]罗奇峰,石春香.Hilbert-Huang变换理论及其计算中的问题.同济大学学报,2003,31(6):637-640
    [64]黄大吉,赵进平,苏纪兰等.希尔伯特-黄变换的端点延拓.海洋学报,2003,25(1):1-11.
    [65]邓拥军,王伟,钱成春等.EMD方法及Hilbert变换中边界问题的处理.科学通报,2001,46(3):257-263.
    [66]张郁山,梁建文,胡垏贤.应用自回归模型处理EMD方法中的边界问题.自然科学进展,2003,13(10):1054-1059.
    [67]刘慧婷,张旻,程家兴.基于多项式拟合算法的EMD端点问题的处理.计算机工程与应用.2004,40(16):84-86.
    [68]龙思胜,张铁宝,龙峰.希尔伯特-黄变换中拟合过冲和端点飞翼的原因及解决办法.地震学报,2005,27(5):561-567
    [69]杜爱明,王彬,杨润海.Hilbert-Huang变换中的一种端点处理方法.地震研究,2007,30(1):54-58
    [70]李庆扬,王能超,易大义.数值分析(第4版).北京:清华大学出版社,2001.45-46
    [71]李岳生,黄友谦.数值逼近.北京:人民教育出版社,1979.47-48
    [72]李岳生.样条与插值.上海:上海科学技术出版社,1983.120-121
    [73]李岳生,齐东旭.样条函数方法.北京:科学出版社,1979.54-55
    [74]J.H.Ahlberg,E.N.Nilson,J.L.Walsh.样条理论及其应用.赵根榕,程正兴译.北京:科学出版社,1981.10-12
    [75]M.H.Schultz.样条分析.赵根榕译.上海:上海科学技术出版社,1979.71-71
    [76]徐利治,王仁宏,周蕴时.函数逼近的理论与方法.上海:上海科学技术出版社,1983.54-57
    [77]G.G.Lorentz.函数逼近论.谢庭藩,施咸亮译.上海:上海科学技术出版社,1981.48-48
    [78]A.A.G.Requicha.The zeros of entire functions:theory and engineering applications.Proceedings of the IEEE,1980,68(3):308-328
    [79]T.Schanze.Sine interpolation of discrete periodic signals.IEEE Trans.Signal Process.,1995,43(6):1502-1503.
    [80]F.Candocia and J.C.Principe.Comments on "sinc interpolation of discrete periodic signals".IEEE Trans.Signal Process.,1998,46(7):2044-2047.
    [81]S.R.Dooley and A.K.Nandi.Notes on the interpolation of discrete periodic signals using sine function related approaches.IEEE Trans.Signal Process.,2000,48(4):1201-1203.
    [82]G.Helmberg.A limit function for equidistant Fourier interpolation.J.Approx.Theory,1995,81:389-396.
    [83]G.Helmberg and P.Wagner.Manipulating Gibbs' Phenomenon for Fourier Interpolation.J.Approx.Theory,1997,89:308-320.
    [84]K.Yao and J.B.Thomas.On truncation error bounds for sampling representations of band-limited signals.IEEE Trans.Aerospace and Electronic Systems,1966,AES-2(6):640-647.
    [85]杨波尔斯基.双曲函数.邢富冲译.北京:中央民族学院出版社,1987
    [86]R.J.Duffin and A.C.Schaeffer.Some properties of functions of exponential type.Bull.~Am.~Math.~Soc.,1938,44:236-240.
    [87]盖广洪.经验模态分解的一种改进算法.西安交通大学学报.2004,38(11):1199-1202
    [88]Y.Kopsinis and S.McLaughlin.Improved EMD using doubly-iterative sifting and high order spline interpolation.EURASIP Journal on Advances in Signal Processing,2008,Article ID 128293.
    [89]杨世锡,胡劲松,吴昭同等.基于高次样条插值的经验模态分解方法研究.浙江大学学报.2004,38(3):267-270
    [90]钟佑明,金涛,秦树人.希尔伯特_黄变换中的一种新包络线算法.数据采集与处理.2005,20(1):13-17
    [91]刘霖雯,刘超,江成顺.EMD新算法及其应用.系统仿真学报.2007,19(2):446-447
    [92]S.D.Hawley,L.E.Atlas and H.J.Chizeck.Some properties of an empirical mode type signal decomposition algorithm.ICASSP 2008:3625-3628
    [93]张勇,陈滨.黄变换的穿越筛分法.电子科技大学学报.2007,36(1):8-10
    [94]钟佑明,秦树人,汤宝平.Hilbert-Huang变换中的理论研究.振动与冲击.2002,21(4):13-17
    [95]钟佑明,秦树人.HHT的理论依据探讨——Hilbert变换的局部乘积定理.振动与冲击.2006,25(2):12-19
    [96]钟佑明,秦树人.希尔伯特-黄变换的统一理论依据研究.振动与冲击.2006,25(3):40-43

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700