Adaptive Fourier decomposition‐based Dirac type time‐frequency distribution

The Dirac‐type time‐frequency distribution (TFD), regarded as ideal TFD, has long been desired. It, until the present time, cannot be implemented, due to the fact that there has been no appropriate representation of signals leading to such TFD. Instead, people have been developing other types of TFD, including the Wigner and the windowed Fourier transform types. This paper promotes a practical passage leading to a Dirac‐type TFD. Based on the proposed function decomposition method, viz., adaptive Fourier decomposition, we establish a rigorous and practical Dirac‐type TFD theory. We do follow the route of analytic signal representation of signals founded and developed by Garbo, Ville, Cohen, Boashash, Picinbono, and others. The difference, however, is that our treatment is theoretically throughout and rigorous. To well illustrate the new theory and the related TFD, we include several examples and experiments of which some are in comparison with the most commonly used TFDs. Copyright © 2016 John Wiley & Sons, Ltd.     

Composite generalized Laguerre spectral method for nonlinear Fokker–Planck equations on the whole line

In this paper, we propose a composite Laguerre spectral method for the nonlinear Fokker–Planck equations modelling the relaxation of fermion and boson gases. A composite Laguerre spectral scheme is constructed. Its convergence is proved. Numerical results show the efficiency of this approach and coincide well with theoretical analysis. Some results on the Laguerre approximation and techniques used in this paper are also applicable to other nonlinear problems on the whole line. Copyright © 2016 John Wiley & Sons, Ltd.     

A fast algorithm for soil dynamics calculations by wavelet decomposition

Summary The article presents a fast numerical algorithm for calculating the response of a halfspace under any surface loads. Under certain conditions there exists an analytical solution to the problem in the Fourier domain. To get the desired response, a numerical inverse Fourier transform of this analytic solution has to be made. By using a wavelet decomposition, the proposed algorithm can reduce the calculation time significantly, thus allowing the computation of complex problems. As an example, the response of the beam-halfspace coupled system under moving load is presented. Received 6 March 1997; accepted for publication 20 May 1997     

A SINGULAR VALUE DECOMPOSITION BASED TRUNCATION ALGORITHM IN SOLVING THE STRUCTURAL DAMAGE EQUATIONS

The structural damage identification through modal data often leads to solving a set of linear equations. Special numerical treatment is sometimes required for an accurate and stable solution owing to the ill conditioning of the equations. Based on the singular value decomposition (SVD) of the coefficient matrix, an error based truncation algorithm is proposed in this paper. By rejection of selected small singular values, the influence of noise can be reduced. A simply-supported beam is used as a simulation example to compare the results to other methods.Illustrative numerical examples demonstrate the good efficiency and stability of the algorithm in the nondestructive identification of structural damage through modal data.     

基于区域分解算法的地下水耦合模型及其应用

根据岩体中不同区域裂隙发育规模的差异分为连续区域和离散区域,在不同区域中使用不同的地下水运动数学模型,应用区域分解算法来解决这类问题。其中连续区域采用了等效连续介质模型,离散区域采用了随机裂隙网络模型,通过区域公共边界上水位和流量连续的条件将两模型耦合求解。将基于区域分解算法的耦合模型应用于锦屏水电站坝址区三维渗流场的模拟中,通过钻孔观测水位和计算水位的对比发现,该方法是有效的,能够应用于实际工程。     

THE THEORETICAL COST OF SEQUENTIAL AND PARALLEL ALGORITHMS FOR SOLVING LINEAR SYSTEMS OF EQUATIONS

ＴＨＥ　ＴＨＥＯＲＥＴＩＣＡＬ　ＣＯＳＴ　ＯＦ　ＳＥＱＵＥＮＴＩＡＬ　ＡＮＤ　ＰＡＲＡＬＬＥＬ　ＡＬＧＯＲＩＴＨＭＳ　ＦＯＲ　ＳＯＬＶＩＮＧ　ＬＩＮＥＡＲ　ＳＹＳＴＥＭＳ　ＯＦ　ＥＱＵＡＴＩＯＮＳＳａｌｍａｎＨ．Ａｂｂａｓ（ＲｅｃｅｉｖｅｄＭａｙ２４...     

Interface crack between two shear deformable elastic layers

An improved method based on the first-order shear deformable plate theory is developed to calculate the energy release rate and stress intensity factor for a crack at the interface of a bi-layer structure. By modeling the uncracked region of the structure as two separate Reissner-Mindlin plates bonded perfectly along the interface, this method is able not only to take into account the shear deformation in the cracked region, but also to capture the shear deformation in the uncracked region of the structure. A closed form solution of energy release rate and mode decomposition at the interface crack is obtained for a general loading condition, and it indicates that the energy release rate and stress intensity factor are determined by two independent loading parameters. Compared to the approach based on the classical plate theory, the proposed method provides a more accurate prediction of energy release rate as well as mode decomposition. The computational procedures introduced are relatively straightforward, and the closed form solution can be used to predict crack growth along the layered structures.     

Nonlinear Coherence in Multivariate Research: Invariants and the Reconstruction of Attractors

First, some linear techniques in multivariate time-series analysis in EEG research are reviewed to highlight the problem of estimating the dimensionality of the state space (embedding dimension), the reconstruction of an attractor, and the evaluation of invariant properties of the attractor. The traditional linear techniques included the usual spectral and cospectral measures of power, phase, and coherence to which stepwise discriminant analysis was applied for canonical representation of the attractor. Then, some traditional nonlinear techniques of attractor reconstruction and dimensional analysis which use the time-lagged univariate approach of Ruelle and Takens (Takens, 1981) are reviewed. Next, updates and multivariate generalizations that use singular-value decomposition (Broomhead & King, 1986) are reviewed. Finally, Stewart's (1995, 1996) multivariate generalization of the method of false nearest neighbors (Abarbanel, Brown, Sidorowich, & Tsimring, 1993; Kennel, Brown, & Abarbanel, 1992) is reviewed. These are particularly relevant for evaluating multivariate coherence in research on the complex cooperative dynamical systems found in neuroscience, psychology, and social science when time series of sufficient length are investigated.     

Finite Precision Representation of the Conley Decomposition

We present a theoretical basis for a novel way of studying and representing the long-time behavior of finite-dimensional maps. It is based on a finite representation of -pseudo orbits of a map by the sample paths of a suitable Markov chain based on a finite partition of phase space. The use of stationary states of the chain and the corresponding partition elements in approximating the attractors of maps and differential equations was demonstrated in Refs. 7 and 3 and proved for a class of stable attracting sets in Ref. 8. Here we explore the relationship between the communication classes of the Markov chain and basic sets of the Conley decomposition of a dynamical system. We give sufficient conditions for the existence of a chain transitive set and show that basic sets are isolated from each other by neighborhoods associated with closed communication classes of the chain. A partition element approximation of an isolating block is introduced that is easy to express in terms of sample paths. Finally, we show that when the map supports an SBR measure there is a unique closed communication class and the Markov chain restricted to those states is irreducible.     

A spectral element method using the modal basis and its application in solving second‐order nonlinear partial differential equations

We present a high‐order spectral element method (SEM) using modal (or hierarchical) basis for modeling of some nonlinear second‐order partial differential equations in two‐dimensional spatial space. The discretization is based on the conforming spectral element technique in space and the semi‐implicit or the explicit finite difference formula in time. Unlike the nodal SEM, which is based on the Lagrange polynomials associated with the Gauss–Lobatto–Legendre or Chebyshev quadrature nodes, the Lobatto polynomials are used in this paper as modal basis. Using modal bases due to their orthogonal properties enables us to exactly obtain the elemental matrices provided that the element‐wise mapping has the constant Jacobian. The difficulty of implementation of modal approximations for nonlinear problems is treated in this paper by expanding the nonlinear terms in the weak form of differential equations in terms of the Lobatto polynomials on each element using the fast Fourier transform (FFT). Utilization of the Fourier interpolation on equidistant points in the FFT algorithm and the enough polynomial order of approximation of the nonlinear terms can lead to minimize the aliasing error. Also, this approach leads to finding numerical solution of a nonlinear differential equation through solving a system of linear algebraic equations. Numerical results for some famous nonlinear equations illustrate efficiency, stability and convergence properties of the approximation scheme, which is exponential in space and up to third‐order in time. Copyright © 2014 John Wiley & Sons, Ltd.