Time-Frequency Analysis of Nonlinear Systems: The Skeleton Linear Model and the Skeleton Curves

The joint time-frequency analysis method is adopted to study the nonlinear behavior varying with the instantaneous response for a class of S.D.O.F nonlinear system. A time-frequency masking operator, together with the conception of effective time-frequency region of the asymptotic signal are defined here. Based on these mathematical foundations, a so-called skeleton linear model (SLM) is constructed which has similar nonlinear characteristics with the nonlinear system. Two skeleton curves are deduced which can indicate the stiffness and damping in the nonlinear system. The relationship between the SLM and the nonlinear system, both parameters and solutions, is clarified. Based on this work a new identification technique of nonlinear systems using the nonstationary vibration data will be proposed through time-frequency filtering technique and wavelet transform in the following paper.

Identification of nonlinear dynamic systems: Time-frequency filtering and skeleton curves

The nonlinear behavior varying with the instantaneous response was analyzed through the joint time-frequency analysis method for a class of S. D. O. F nonlinear system. A masking operator an definite regions is defined and two theorems are presented. Based on these, the nonlinear system is modeled with a special time-varying linear one, called the generalized skeleton linear system (GSLS). The frequency skeleton curve and the damping skeleton curve are defined to describe the main feature of the non-linearity as well. Moreover, an identification method is proposed through the skeleton curves and the time-frequency filtering technique.

Identification of Nonlinear Systems Through Time-Frequency Filtering Technique

In the previous paper, a class of nonlinear system is mapped to a so-called skeleton linear model (SLM) based on the joint time-frequency analysis method. Behavior of the nonlinear system may be indicated quantitatively by the variance of the coefficients of SLM versus its response. Using this model we propose an identification method for nonlinear systems based on nonstationary vibration data in this paper. The key technique in the identification procedure is a time-frequency filtering method by which solution of the SLM is extracted from the response data of the corresponding nonlinear system. Two time-frequency filtering methods are discussed here. One is based on the quadratic time-frequency distribution and its inverse transform, the other is based on the quadratic time-frequency distribution and the wavelet transform. Both numerical examples and an experimental application are given to illustrate the validity of the technique.

Electron quantum path tuning and isolated attosecond pulse emission driven by a waveform-controlled multi-cycle laser field

We investigate high-order harmonic emission and isolated attosecond pulse (IAP) generation in atoms driven by a two-colour multi-cycle laser field consisting of an 800 nm pulse and an infrared laser pulse at an arbitrary wavelength. With moderate laser intensity, an IAP of similar to 220 as can be generated in helium atoms by using two-colour laser pulses of 35 fs/800 nm and 46 fs/1150 nm. The discussion based on the three-step semiclassical model, and time-frequency analysis shows a clear picture of the high-order harmonic generation in the waveform-controlled laser field which is of benefit to the generation of XUV IAP and attosecond electron pulses. When the propagation effect is included, the duration of the IAP can be shorter than 200 as, when the driving laser pulses are focused 1 mm before the gas medium with a length between 1.5 mm and 2 mm.

Theoretical study on the mechanism of the reaction of CH4+MgO

The reactions of (1) CH4 + MgO --> MgOH. + CH3. and (2) CH4 + MgO --> Mg + CH3OH have been studied on the singlet spin state potential energy surface at the MP2/6-311+G(2d,2p) level. These two reaction channels, both involving intermediates and transition states, have been rationalized by the structures of the species involved, natural bond orbital (NBO), and vibrational frequency analysis. We have considered two initial interacting models between CH4 and MgO: a collinear C-H approach to the O end of the MgO forming the MgOCH4 complex with C-3nu symmetry and three hydrogen atoms of the methane point to the Mg end of the MgO forming the OMgCH4 complex with C-1 symmetry. The calculations predict that reactions 1 and 2 are exothermic by 39.8 and 86.5 kJ mol(-1), respectively. Also, the former reaction proceeds more easily than the latter, and the complex HOMgCH3 is energetically preferred in the reaction of MgO + CH4.

Using otolith weight to predict the age of Pennahia macrocephalus in the mouth of the Beibu Gulf

The relation between otolith weight (OW) and the age of marine fish is studied. A total of 222 individuals of bighead white croaker, Pennahia macrocephalus were sampled seasonally in the mouth of the Beibu Gulf, the South China Sea, in 2007. Since there are no significant differences in sagittal OW between otolith in pairs (Pa parts per thousand yen0.05), the undamaged left sagittal otolith is used for age determination. The highest correlations among standard length, OW and fish ages are confirmed by linear, exponential and multinomial regression. Results show that sagittal OW overlaps only occasionally among age groups, and to individuals with similar standard length, the older and slower-growing fish has a heavier otolith because of the continued otolith material deposition. There are differences in sagittal OW among different age groups and significant positive linear relationship with age (P < 0.05). The age readings can be verified by plotting the sagittal OW versus the standard length for age groups, and the individuals with similar standard length but in different ages can be separated by sagittal OW frequency analysis. Mostly, the predicted ages using the regression between sagittal OW and ages are closed to the observed ages by counting annulus on scale. It indicates that the sagittal OW analysis is a useful technique for validating the accuracy of age determination by annuli counts, especially for individuals of similar size. Furthermore, the technique is applied for Pennahia macrocephalus with discussion in this paper.

THE COMPOUND EYE AS AN INDICATOR OF AGE AND SHRINKAGE IN ANTARCTIC KRILL

Laboratory studies have shown that Antarctic krill (Euphausia superba) shrink if maintained in conditions of low food availability. Recent studies have also demonstrated that E. superba individuals may be shrinking in the field during winter. If krill shrink during the winter, conclusions reached by length-frequency analysis may be unreliable because smaller animals may not necessarily be younger animals. In this study, the correlation between the body-length and the crystalline cone number of the compound eye was examined. Samples collected in the late summer show an apparent linear relationship between crystalline cone number and body-length. From a laboratory population, it appears that when krill shrink the crystalline cone number remains relatively unchanged. If crystalline cone number is little affected by shrinking, then the crystalline cone number may be a more reliable indicator of age than body-length alone. The ratio of crystalline cone number to body-length offers a method for detecting the effect of shrinking in natural populations of krill. On the basis of the crystalline cone number count, it appears from a field collection in early spring that E. superba do shrink during winter.

基于匹配追踪方法的地震信号时频分析研究

Seismic signal is a typical non-stationary signal, whose frequency is continuously changing with time and is determined by the bandwidth of seismic source and the absorption characteristic of the media underground. The most interesting target of seismic signal’s processing and explaining is to know about the local frequency’s abrupt changing with the time, since this kind of abrupt changing is indicating the changing of the physical attributes of the media underground. As to the seismic signal’s instantaneous attributes taken from time-frequency domain, the key target is to search a effective, non-negative and fast algorithm time-frequency distribution, and transform the seismic signal into this time-frequency domain to get its instantaneous power spectrum density, and then use the process of weighted adding and average etc. to get the instantaneous attributes of seismic signal. Time-frequency analysis as a powerful tool to deal with time variant non-stationary signal is becoming a hot researching spot of modern signal processing, and also is an important method to make seismic signal’s attributes analysis. This kind of method provides joint distribution message about time domain and frequency domain, and it clearly plots the correlation of signal’s frequency changing with the time. The spectrum decomposition technique makes seismic signal’s resolving rate reach its theoretical level, and by the method of all frequency scanning and imaging the three dimensional seismic data in frequency domain, it improves and promotes the resolving abilities of seismic signal vs. geological abnormal objects. Matching pursuits method is an important way to realize signal’s self-adaptive decomposition. Its main thought is that any signal can be expressed by a series of time-frequency atoms’ linear composition. By decomposition the signal within an over completed library, the time-frequency atoms which stand for the signal itself are selected neatly and self-adaptively according to the signal’s characteristics. This method has excellent sparse decomposition characteristics, and is widely used in signal de-noising, signal coding and pattern recognizing processing and is also adaptive to seismic signal’s decomposition and attributes analysis. This paper takes matching pursuits method as the key research object. As introducing the principle and implementation techniques of matching pursuits method systematically, it researches deeply the pivotal problems of atom type’s selection, the atom dictionary’s discrete, and the most matching atom’s searching algorithm, and at the same time, applying this matching pursuits method into seismic signal’s processing by picking-up correlative instantaneous messages from time-frequency analysis and spectrum decomposition to the seismic signal. Based on the research of the theory and its correlative model examination of the adaptively signal decomposition with matching pursuit method, this paper proposes a fast optimal matching time-frequency atom’s searching algorithm aimed at seismic signal’s decomposition by frequency-dominated pursuit method and this makes the MP method pertinence to seismic signal’s processing. Upon the research of optimal Gabor atom’s fast searching and matching algorithm, this paper proposes global optimal searching method using Simulated Annealing Algorithm, Genetic Algorithm and composed Simulated Annealing and Genetic Algorithm, so as to provide another way to implement fast matching pursuit method. At the same time, aimed at the characteristics of seismic signal, this paper proposes a fast matching atom’s searching algorithm by means of designating the max energy points of complex seismic signal, searching for the most optimal atom in the neighbor area of these points according to its instantaneous frequency and instantaneous phase, and this promotes the calculating efficiency of seismic signal’s matching pursuit algorithm. According to these methods proposed above, this paper implements them by programmed calculation, compares them with some open algorithm and proves this paper’s conclusions. It also testifies the active results of various methods by the processing of actual signals. The problems need to be solved further and the aftertime researching targets are as follows: continuously seeking for more efficient fast matching pursuit algorithm and expanding its application range, and also study the actual usage of matching pursuit method.

小波域地震资料的处理方法（WSSP）研究

The modeling formula based on seismic wavelet can well simulate zero - phase wavelet and hybrid-phase wavelet, and approximate maximal - phase and minimal - phase wavelet in a certain sense. The modeling wavelet can be used as wavelet function after suitable modification item added to meet some conditions. On the basis of the modified Morlet wavelet, the derivative wavelet function has been derived. As a basic wavelet, it can be sued for high resolution frequency - division processing and instantaneous feature extraction, in acoordance with the signal expanding characters in time and scale domains by each wavelet structured. Finally, an application example proves the effectiveness and reasonability of the method. Based on the analysis of SVD (Singular Value Decomposition) filter, by taking wavelet as basic wavelet and combining SVD filter and wavelet transform, a new de - noising method, which is Based on multi - dimension and multi-space de - noising method, is proposed. The implementation of this method is discussed the detail. Theoretical analysis and modeling show that the method has strong capacity of de - noising and keeping attributes of effective wave. It is a good tool for de - noising when the S/N ratio is poor. To give prominence to high frequency information of reflection event of important layer and to take account of other frequency information under processing seismic data, it is difficult for deconvolution filter to realize this goal. A filter from Fourier Transform has some problems for realizing the goal. In this paper, a new method is put forward, that is a method of processing seismic data in frequency division from wavelet transform and reconstruction. In ordinary seismic processing methods for resolution improvement, deconvolution operator has poor part characteristics, thus influencing the operator frequency. In wavelet transform, wavelet function has very good part characteristics. Frequency - division data processing in wavelet transform also brings quite good high resolution data, but it needs more time than deconvolution method does. On the basis of frequency - division processing method in wavelet domain, a new technique is put forward, which involves 1) designing filter operators equivalent to deconvolution operator in time and frequency domains in wavelet transform, 2) obtaining derivative wavelet function that is suitable to high - resolution seismic data processing, and 3) processing high resolution seismic data by deconvolution method in time domain. In the method of producing some instantaneous characteristic signals by using Hilbert transform, Hilbert transform is very sensitive to high - frequency random noise. As a result, even though there exist weak high - frequency noises in seismic signals, the obtained instantaneous characteristics of seismic signals may be still submerged by the noises. One method for having instantaneous characteristics of seismic signals in wavelet domain is put forward, which obtains directly the instantaneous characteristics of seismic signals by taking the characteristics of both the real part (real signals, namely seismic signals) and the imaginary part (the Hilbert transfom of real signals) of wavelet transform. The method has the functions of frequency division and noise removal. What is more, the weak wave whose frequency is lower than that of high - frequency random noise is retained in the obtained instantaneous characteristics of seismic signals, and the weak wave may be seen in instantaneous characteristic sections (such as instantaneous frequency, instantaneous phase and instantaneous amplitude). Impedance inversion is one of tools in the description of oil reservoir. one of methods in impedance inversion is Generalized Linear Inversion. This method has higher precision of inversion. But, this method is sensitive to noise of seismic data, so that error results are got. The description of oil reservoir in researching important geological layer, in order to give prominence to geological characteristics of the important layer, not only high frequency impedance to research thin sand layer, but other frequency impedance are needed. It is difficult for some impedance inversion method to realize the goal. Wavelet transform is very good in denoising and processing in frequency division. Therefore, in the paper, a method of impedance inversion is put forward based on wavelet transform, that is impedance inversion in frequency division from wavelet transform and reconstruction. in this paper, based on wavelet transform, methods of time - frequency analysis is given. Fanally, methods above are in application on real oil field - Sansan oil field.

陡坡带砂砾岩体油气成藏模式研究——以东营凹陷西北部陡坡带为例

Based on the study of sequence stratigraphy, modern sedimentary, basin analysis, and petroleum system in abrupt slop of depression, this paper builds sedimentary system and model, sandy bodies distribution, and pool-forming mechanism of subtle trap. There are some conclusions and views as follows. By a lot of well logging and seismic analysis, the author founded up the sequence stratigraphic of the abrupt slope, systematically illustrated the abrupt slope constructive framework, and pointed out that there was a special characteristics which was that south-north could be divided to several fault block and east-west could be carved up groove and the bridge in studying area. Based all these, the author divided the studying area to 3 fault block zone in which because of the groove became the basement rock channel down which ancient rivers breathed into the lake, the alluvial fan or fan delta were formed. In the paper, the author illustrated the depositional system and depositional model of abrupt slope zone, and distinguished 16 kinds of lithofacies and 3 kinds of depositional systems which were the alluvial fan and fan-delta system, lake system and the turbidite fan or turbidity current deposition. It is first time to expound completely the genetic pattern and distributing rule of the abrupt slope sandy-conglomeratic fan bodies. The abrupt slope sandy-conglomeratic fan bodies distribute around the heaves showing itself circularity shape. In studying area, the sandy-conglomeratic fan bodies mainly distribute up the southern slope of Binxian heave and Chenjiazhuang heave. There mainly are these sandy-conglomeratic fan body colony which distributes at a wide rage including the alluvial fan, sub-water fluvial and the turbidite fan or the other turbidity current deposition in the I fault block of the Wangzhuang area. In the II fault block there are fan-delta front and sub-water fluvial. And in the Binnan area, there mainly are those the alluvial fan (down the basement rock channel) and the sandy-conglomeratic fan body which formed as narrowband sub-water fluvial (the position of bridge of a nose) in the I fault block, the fan-delta front sandy-conglomeratic fan body in the H fault block and the fan-delta front and the turbidity current deposition sandy-conglomeratic fan body in the m fault block. Based on the reservoir outstanding characteristics of complex classic composition and the low texture maturity, the author comparted the reservoir micro-structure of the Sha-III and Sha-IV member to 4 types including the viscous crude cementation type, the pad cementation type, the calcite pore-funds type and the complex filling type, and hereby synthetically evaluated 4 types sandy- conglomeratic fan body reservoir. In the west-north abrupt slope zone of Dongying Depression, the crude oil source is belonging to the Sha-III and Sha-IV member, the deep oil of Lijin oilfield respectively come from the Sha-III and Sha-IV member, which belongs to the autogeny and original deposition type; and the more crude oil producing by Sha-IV member was migrated to the Wangzhuan area and Zhengjia area. The crude oil of Binnan oil-field and Shanjiasi oil-field belongs to mixed genetic. It is the first time to illustrate systematically the genetic of the viscous crude that largely being in the studying area, which are that the dissipation of the light component after pool-forming, the biological gradation action and the bath-oxidation action, these oil accumulation belonging to the secondary viscous crude accumulation. It is also the first time to compart the studying area to 5 pool-forming dynamical system that have the characteristic including the common pressure and abnormal pressure system, the self-fountain and other-fountain system and the closing and half-closing system etc. The 5 dynamical systems reciprocally interconnected via the disappearance or merger of the Ethology and the fluid pressure compartment zone, the fault and the unconformity surface, hereby formed duplicated pattern oil-gas collecting zone. Three oil-gas pool-forming pattern were founded, which included the self-fountain side-direction migrated collecting pattern, the self-fountain side-direction ladder-shape pool-forming pattern and the other-fountain pressure releasing zone migrated collecting pattern. A series of systemic sandy-conglomeratic fan bodies oil-gas predicting theory and method was founded, based on the groove-fan corresponding relation to confirm the favorable aim area, according as the characteristic of seismic-facies to identify qualitatively the sandy-conglomeratic fan bodies or its scale, used the temporal and frequency analysis technique to score the interior structure of the sandy- conglomeratic fan bodies, applied for coherent-data system analysis technology to describe the boundary of the sandy-conglomeratic fan bodies, and utilized the well logging restriction inversion technique to trace quantificational and forecast the sandy-conglomeratic fan bodies. Applied this technique, totally 15 beneficial sandy-conglomeratic fan bodies were predicted, in studying area the exploration was preferably guided, and the larger economic benefit and social benefit was acquired.

中学生日常烦恼和缓冲因素的研究

This study is a research on the stress of the secondary school students. The whole work included the construction of the Daily Hassle Scale of the Secondary School Students and the research on the mediate factors of the stressors. Six kinds of the daily hassle were found by the factor analysis: problem with academic performance, problem with school life, problem with teachers, problem with family, problem with heterosexual relation, problem with self-esteem. The reliability and validity of the scale is high. By the regression analysis of the SCL-90 score, we found that the factors which influence the physical and mental health of the secondary school students are the personal character(neuroticism, locus of control), daily hassle (problem with academic performance, problem with heterosexual relation) and age.

大学生应激应付研究

In this study, we made research on coping with stress in undergraduate students. The whole work included two parts: (a) We developed the Coping Strategies Questionnaire and had it tested. Five kinds of coping strategies were found by factor analysis: Problem Solving, Passive Acceptance, Support Seeking, Emotion Regulation, and Unadaptation. The reliability and validity coefficients of the questionnaire was high. (b) We examined the relations among personality characteristics (Locus of control, Extraversion, Neuroticism), emotion state, cognitive appraisal (psychological pressure, controllability, resource of social support), and the specific coping processes in several specific stressful events. The results indicated that coping have some degree of consistency across different situations which suggested that there is a general coping tendency, and that this general coping tendency, in addition to the above considered variables, have different degree of effects on subjects' choice of specific coping strategies. Finally, theoretical model of the relationship was tested through path analysis.

Analysis of strip electric saturation model of crack problem in piezoelectric materials

This paper presents a fully anisotropic analysis of strip electric saturation model proposed by Gao et al. (1997) (Gao, H.J., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids, 45, 491-510) for piezoelectric materials. The relationship between the size of the strip saturation zone ahead of a crack tip and the applied electric displacement field is established. It is revealed that the critical fracture stresses for a crack perpendicular to the poling axis is linearly decreased with the increase of the positive applied electric field and increases linearly with the increase of the negative applied electric field. For a crack parallel to the poring axis, the failure stress is not effected by the parallel applied electric field. In order to analyse the existed experimental results, the stress fields ahead of the tip of an elliptic notch in an infinite piezoelectric solid are calculated. The critical maximum stress criterion is adopted for determining the fracture stresses under different remote electric displacement fields. The present analysis indicates that the crack initiation and propagation from the tip of a sharp elliptic notch could be aided or impeded by an electric displacement field depending on the field direction. The fracture stress predicted by the present analysis is consistent with the experimental data given by Park and Sun (1995) (Park, S., Sun, C.T., 1995. Fracture criteria for piezoelectric materials. J. Am. Ceram. Soc 78, 1475-1480).

Three-dimensional crack problem analysis using boundary element method with finite-part integrals

A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor solutions to several typical planar crack problems in a finite body are evaluated.