HHT在电磁信号处理中的应用

In modem signal Processing，non－linear，non－Gaussian and non－stable signals are usually the analyzed and Processed objects，especially non－stable signals. The convention always to analyze and Process non－stable signals are： short time Fourier transform，Wigner－Ville distribution，wavelet Transform and so on. But the above three algorithms are all based on Fourier Transform，so they all have the shortcoming of Fourier Analysis and cannot get rid of the localization of it. Hilbert－Huang Transform is a new non－stable signal processing technology，proposed by N. E. Huang in 1998. It is composed of Empirical Mode Decomposition (referred to as EMD) and Hilbert Spectral Analysis (referred to as HSA). After EMD Processing，any non－stable signal will be decomposed to a series of data sequences with different scales. Each sequence is called an Intrinsic Mode Function (referred to as IMF). And then the energy distribution plots of the original non－stable signal can be found by summing all the Hilbert spectrums of each IMF. In essence，this algorithm makes the non－stable signals become stable and decomposes the fluctuations and tendencies of different scales by degrees and at last describes the frequency components with instantaneous frequency and energy instead of the total frequency and energy in Fourier Spectral Analysis. In this case，the shortcoming of using many fake harmonic waves to describe non－linear and non－stable signals in Fourier Transform can be avoided. This Paper researches in the following parts： Firstly，This paper introduce the history and development of HHT，subsequently the characters and main issues of HHT. This paper briefly introduced the basic realization principles and algorithms of Hilbert－Huang transformation and confirms its validity by simulations. Secondly, This paper discuss on some shortcoming of HHT. By using FFT interpolation, we solve the problem of IMF instability and instantaneous frequency undulate which are caused by the insufficiency of sampling rate. As to the bound effect caused by the limitation of envelop algorithm of HHT, we use the wave characteristic matching method, and have good result. Thirdly, This paper do some deeply research on the application of HHT in electromagnetism signals processing. Based on the analysis of actual data examples, we discussed its application in electromagnetism signals processing and noise suppression. Using empirical mode decomposition method and multi-scale filter characteristics can effectively analyze the noise distribution of electromagnetism signal and suppress interference processing and information interpretability. It has been founded that selecting electromagnetism signal sessions using Hilbert time-frequency energy spectrum is helpful to improve signal quality and enhance the quality of data.

各向异性与非线性介质中弹性波场传播特征研究

The topic of this study is about the propagation features of elastic waves in the anisotropic and nonlinear media by numerical methods with high accuracy and stability. The main achievements of this paper are as followings: Firstly, basing on the third order elastic energy formula, principle of energy conservation and circumvolved matrix method, we firstly reported the equations of non-linear elastic waves with two dimensions and three components in VTI media. Secondly, several conclusions about some numerical methods have been obtained in this paper. Namely, the minimum suitable sample stepth in space is about 1/8-1/12 of the main wavelength in order to distinctly reduce the numerical dispersion resulted from the numerical mehtod, at the same time, the higher order conventional finite difference (CFD) schemes will give little contribution to avoid the numerical solutions error accumulating with time. To get the similar accuracy with the fourth order center finite difference method, the half truncation length of SFFT should be no less than 7. The FDFCT method can present with the numerical solutions without obvious dispersion when the paprameters of FCT is suitable (we think they should be in the scope from 0.0001 to 0.07). Fortunately, the NADM method not only can reported us with the higher order accuracy solutions (higher than that of the fourth order finite difference method and lower than that of the sixth order finite difference method), but also can distinctly reduce the numerical dispersion. Thirdly, basing on the numerial and theoretical analysis, we reported such nonlinear response accumulating with time as waveform aberration, harmonic generation and resonant peak shift shown by the propagation of one- and two-dimensional non-linear elasticwaves in this paper. And then, we drew the conclusion that these nonlinear responses are controlled by the product between nonlinear strength (SN) and the amplitude of the source. At last, the modified FDFCT numerical method presented by this paper is used to model the two-dimensional non-linear elastic waves propagating in VTI media. Subsequently, the wavelet analysis and polarization are adopted to investigate and understand the numerical results. And then, we found the following principles (attention: the nonlinear strength presented by this paper is weak, the thickness of the -nonlinear media is thin (200m), the initial energy of the source is weak and the anisotropy of the media is weak too): The non-linear response shown by the elastic waves in VTI media is anisotropic too; The instantaneous main frequency sections of seismic records resulted from the media with a non-linear layer have about 1/4 to 1/2 changes of the initial main frequency of source with that resulted from the media without non-linear layer; The responses shown by the elasic waves about the anisotropy and nonlinearity have obvious mutual reformation, namely, the non-linear response will be stronger in some directions because of the anisotropy and the anisotropic strength shown by the elastic waves will be stronger when the media is nonlinear.

元音和复合音识别中噪音分离机制的实验研究

Harmonicity is one of the important features of a vowel. It makes great contribution to pitch and quality of vowel. However, contribution of a mistuned harmonic will decrease as it is mistuned increasingly. A mistuned harmonic will be segregated as noise from complex by auditory system, which was called harmonic sieve (Duifhuis, 1982). According to Darwin (1986) and Moore et al (1985), the critical value of one mistuned harmonic would be segregated from vowel or complex is 3% to 8%－－Harmonic Mistuned Effect (HME). Further questions need to be answered. For example, how will the harmonic sieve separate noise or whether the critical value change when more than two harmonics are mistuned? And what affect the HME? Three experiments were conducted to these questions. Experiment one was dealt with the number of mistuned harmonics as a factor affecting the HME. The position effect of HME was concerned in experiment two. The last experiment considered the relationship between HME and phase of the mistuned harmonic. The results indicated that (1) the HME was much greater when more than two harmonics were mistuned than only one harmonic was mistuned; (2) harmonic position played an important role in HME, the higher the harmonic was, the less HME was found for the complex, and the closer to formant the harmonic stood, the more significant HME existed; and (3) phase did not affect the HME significantly, however, its indirect contribution still existed, which related to the starting amplitude of a mistuned harmonic.

Multistep Methods for Integrating the Solar System

High order multistep methods, run at constant stepsize, are very effective for integrating the Newtonian solar system for extended periods of time. I have studied the stability and error growth of these methods when applied to harmonic oscillators and two-body systems like the Sun-Jupiter pair. I have also tried to design better multistep integrators than the traditional Stormer and Cowell methods, and I have found a few interesting ones.

Franck-Condon simulation of the photoelectron spectrum of SO2- including Duschinsky effects

A theoretical method to calculate multidimensional Franck-Condon factors including Duschinsky effects is described and used to simulate the photoelectron spectrum of the anion SO. Geometry optimizations and harmonic vibrational frequency calculations have been performed on the XA(1) state of SO2 and (XB1)-B-2 state of SO2. Franck-Condon analyses and spectral simulation were carried out on the first photoelectron band of SO2. The theoretical spectra obtained by employing CCSD(T)/6-31 I+G(2d,p) values are in excellent agreement with the experiment. In addition, the equilibrium geometric parameters, r(c)(OS) = 0.1508 +/- 0.0005 nm and theta(e)(O-S-0) = 113.5 +/- 0.5 degrees, of the (XB1)-B-2 state of SO2, are derived by employing an iterative Franck-Condon analysis procedure in the spectral simulation. (c) 2005 Elsevier B.V. All rights reserved.

Calculation of the multimode Franck-Condon factors based on the coherent state method

The equivalence of two ways for the calculation of overlap integrals, i.e. the Sharp Rosenstock generating function method and the Doktorov coherent state method, has been proved. On the basis of the generating function of the overlap integrals, a new closed form expression for the Franck - Condon integrals for overlap multidimensional harmonic oscillators has been exactly derived. In addition, some useful analytical expressions for the calculations of the multimode Franck - Condon factors have been given.