A direct algebraic method applied to obtain complex solutions of some nonlinear partial differential equations

By using some exact solutions of an auxiliary ordinary differential equation, a direct algebraic method is described to construct the exact complex solutions for nonlinear partial differential equations. The method is implemented for the NLS equation, a new Hamiltonian amplitude equation, the coupled Schrodinger–KdV equations and the Hirota–Maccari equations. New exact complex solutions are obtained.     

The dynamic complexity of a three-species Beddington-type food chain with impulsive control strategy

In this paper, by using theories and methods of ecology and ordinary differential equation, the dynamics complexity of a prey–predator system with Beddington-type functional response and impulsive control strategy is established. Conditions for the system to be extinct are given by using the Floquet theory of impulsive equation and small amplitude perturbation skills. Furthermore, by using the method of numerical simulation with the international software Maple, the influence of the impulsive perturbations on the inherent oscillation is investigated, which shows rich dynamics, such as quasi-periodic oscillation, narrow periodic window, wide periodic window, chaotic bands, period doubling bifurcation, symmetry-breaking pitchfork bifurcation, period-halving bifurcation and crises, etc. The numerical results indicate that computer simulation is a useful method for studying the complex dynamic systems.     

A Multiresolution Method for Numerical Reduction and Homogenization of Nonlinear ODEs

The multiresolution analysis (MRA) strategy for the reduction of a nonlinear differential equation is a procedure for constructing an equation directly for the coarse scale component of the solution. The MRA homogenization process is a method for building a simpler equation whose solution has the same coarse behavior as the solution to a more complex equation. We present two multiresolution reduction methods for nonlinear differential equations: a numerical procedure and an analytic method. We also discuss one possible appproach to the homogenization method.     

Ray Method Expansions for Surface and Internal Waves in Inhomogeneous Oceans of Variable Depth

A ray method formalism is developed for the analysis of surface and internal waves in an inhomogeneous ocean of variable depth. In this method, we deduce from the governing system of equations a system of first order ordinary differential equations, for the group lines (rays of the ray method) and the propagation of phase and amplitude on them. The dispersion relation for these waves arises as an eigen-condition on an eigen-value problem involving an ordinary differential equation in the depth variable. The deduced equation for amplitude propagation has the interpretation of a statement of conservation of action.     

Solitary waves of the EW and RLW equations

Eight finite difference methods are employed to study the solitary waves of the equal-width (EW) and regularized long–wave (RLW) equations. The methods include second-order accurate (in space) implicit and linearly implicit techniques, a three-point, fourth-order accurate, compact operator algorithm, an exponential method based on the local integration of linear, second-order ordinary differential equations, and first- and second-order accurate temporal discretizations. It is shown that the compact operator method with a Crank–Nicolson discretization is more accurate than the other seven techniques as assessed for the three invariants of the EW and RLW equations and the L₂-norm errors when the exact solution is available. It is also shown that the use of Gaussian initial conditions may result in the formation of either positive or negative secondary solitary waves for the EW equation and the formation of positive solitary waves with or without oscillating tails for the RLW equation depending on the amplitude and width of the Gaussian initial conditions. In either case, it is shown that the creation of the secondary wave may be preceded by a steepening and an narrowing of the initial condition. The creation of a secondary wave is reported to also occur in the dissipative RLW equation, whereas the effects of dissipation in the EW equation are characterized by a decrease in amplitude, an increase of the width and a curving of the trajectory of the solitary wave. The collision and divergence of solitary waves of the EW and RLW equations are also considered in terms of the wave amplitude and the invariants of these equations.     

Complex Lie symmetries for scalar second-order ordinary differential equations

A scalar complex ordinary differential equation can be considered as two coupled real partial differential equations, along with the constraint of the Cauchy–Riemann equations, which constitute a system of four equations for two unknown real functions of two real variables. It is shown that the resulting system possesses those real Lie symmetries that are obtained by splitting each complex Lie symmetry of the given complex ordinary differential equation. Further, if we restrict the complex function to be of a single real variable, then the complex ordinary differential equation yields a coupled system of two ordinary differential equations and their invariance can be obtained in a non-trivial way from the invariance of the restricted complex differential equation. Also, the use of a complex Lie symmetry reduces the order of the complex ordinary differential equation (restricted complex ordinary differential equation) by one, which in turn yields a reduction in the order by one of the system of partial differential equations (system of ordinary differential equations). In this paper, for simplicity, we investigate the case of scalar second-order ordinary differential equations. As a consequence, we obtain an extension of the Lie table for second-order equations with two symmetries.     

Solitary waves of the EW and RLW equations

Eight finite difference methods are employed to study the solitary waves of the equal-width (EW) and regularized long–wave (RLW) equations. The methods include second-order accurate (in space) implicit and linearly implicit techniques, a three-point, fourth-order accurate, compact operator algorithm, an exponential method based on the local integration of linear, second-order ordinary differential equations, and first- and second-order accurate temporal discretizations. It is shown that the compact operator method with a Crank–Nicolson discretization is more accurate than the other seven techniques as assessed for the three invariants of the EW and RLW equations and the L₂-norm errors when the exact solution is available. It is also shown that the use of Gaussian initial conditions may result in the formation of either positive or negative secondary solitary waves for the EW equation and the formation of positive solitary waves with or without oscillating tails for the RLW equation depending on the amplitude and width of the Gaussian initial conditions. In either case, it is shown that the creation of the secondary wave may be preceded by a steepening and an narrowing of the initial condition. The creation of a secondary wave is reported to also occur in the dissipative RLW equation, whereas the effects of dissipation in the EW equation are characterized by a decrease in amplitude, an increase of the width and a curving of the trajectory of the solitary wave. The collision and divergence of solitary waves of the EW and RLW equations are also considered in terms of the wave amplitude and the invariants of these equations.     

Meromorphic solutions of an auxiliary ordinary differential equation using complex method

In this paper, we employ the complex method to obtain first all meromorphic solutions of an auxiliary ordinary differential equation and then find all meromorphic exact solutions of the classical Korteweg–de Vries equation, Boussinesq equation, ( 3 + 1)‐dimensional Jimbo–Miwa equation, and Benjamin–Bona–Mahony equation. Our results show that the method is more simple than other methods. Copyright © 2013 John Wiley & Sons, Ltd.     

Terminal-Dependent Statistical Inferences for FBSDE

Backward stochastic differential equation (BSDE) has been well studied and widely applied in mathematical finance. The main difference from the original stochastic differential equation (OSDE) is that the BSDE is designed to depend on a terminal condition, which plays key roles in certain financial and ecological circumstances. However, to the best of our knowledge, the terminal-dependent statistical inference for such model has not been explored in the existing literature. This article proposes two terminal-dependent estimation methods via terminal control variable the integral form models of forward-backward stochastic differential equation (FBSDE). We take these measures because the resulting models contain terminal condition as model variable, and therefore, the corresponding estimators inherit the terminal-dependent characteristic. In this article, the FBSDE is first rewritten as regression versions and then two semi-parametric estimation approaches are proposed. Because of the control variable and integral form, our regression versions are more complex than the classical ones, and the inference methods are somewhat different from which designed for the OSDE. Even so, the statistical properties of the terminal-dependent methods are similar to the classical ones. Simulations are conducted to demonstrate finite sample behaviors.     

具有Holling Ⅲ功能反应和阶段结构的Gompertz生态系统的持久性

以生态学与微分方程的理论和方法为基础,建立了一类具有HollingⅢ功能反应和阶段结构的生态Gompertz模型.利用频闪映射,获得了捕食者灭绝周期解,分析了此周期解的全局吸引性.在对食饵进行脉冲收获和捕食者具有成长期时滞条件下,运用脉冲微分方程比较定理和小振幅扰动技巧,获得了系统一致持续生存的条件.     

Green's function-based integral approaches to linear transient boundary-value problems and their stability characteristics (I)

Two Green's function-based formulations are applied to the governing differential equation which describes unsteady heat or mass transport in an isotropic homogeneous 1-D domain. In this first part of a two series of papers, the linear form of the differential equation is addressed. The first formulation, herein denoted the quasi-steady Green element (QSGE) formulation, uses the Laplace differential operator as auxiliary equation to obtain the singular integral representation of the governing equation, while the second, denoted the transient Green element (TGE), uses the transient heat equation as auxiliary equation. The mathematical simplicity of the Green's function of the first formulation enhances the ease of solution of the integral equations and the resultant discrete equations. From the point of computational convenience, therefore, the first formulation is preferred. The stability characteristics of the two formulations are evaluated by examining how they propagate various Fourier harmonics in speed and amplitude. We found that both formulations correctly reproduce the theoretical speed of the harmonics, but fail to propagate the amplitude of the small harmonics correctly for Courant value of about unity. The QSGE formulation with difference weighting values between 0.67 and 0.75, and the TGE formulation provide optimal performance in numerical stability.     

Properties of the solutions of the fourth-order Bessel-type differential equation

The structured Bessel-type functions of arbitrary even-order were introduced by Everitt and Markett in 1994; these functions satisfy linear ordinary differential equations of the same even-order. The differential equations have analytic coefficients and are defined on the whole complex plane with a regular singularity at the origin and an irregular singularity at the point of infinity. They are all natural extensions of the classical second-order Bessel differential equation. Further these differential equations have real-valued coefficients on the positive real half-line of the plane, and can be written in Lagrange symmetric (formally self-adjoint) form. In the fourth-order case, the Lagrange symmetric differential expression generates self-adjoint unbounded operators in certain Hilbert function spaces. These results are recorded in many of the papers here given as references. It is shown in the original paper of 1994 that in this fourth-order case one solution exists which can be represented in terms of the classical Bessel functions of order 0 and 1. The existence of this solution, further aided by computer programs in Maple, led to the existence of a linearly independent basis of solutions of the differential equation. In this paper a new proof of the existence of this solution base is given, on using the advanced theory of special functions in the complex plane. The methods lead to the development of analytical properties of these solutions, in particular the series expansions of all solutions at the regular singularity at the origin of the complex plane.     

扩充偏微分方程(组)守恒律和对称的辅助方程方法及微分形式吴方法的应用

首先,我们给出了引入伴随方程(组)扩充原方程(组)的策略使给定偏微分方程(组)的扩充方程组具有对应泛瓯即,成为Lagrange系统的方法,以此为基础提出了作为偏微分方程(组)传统守恒律和对称概念的一种推广-偏微分方程(组)扩充守恒律和扩充对称的概念;其次,以得到的Lagrange系统为基础给定了确定原方程(组)扩充守恒律和扩充对称的方法,从而达到扩充给定偏微分方程(组)的首恒律和对称的目的;第三,提出了适用于一般形式微分方程(组)的计算固有守恒律的方法;第四,实现以上算法过程中,我们先把计算(扩充)守恒律和对称问题均归结为求解超定线性齐次偏微分方程组(确定方程组)的问题.然后,对此关键问题我们提出了用微分形式吴方法处理的有效算法;最后,作为方法的应用我们计算确定了非线性电报方程组在内的五个发展方程(组)的新守恒律和对称,同时也说明了方法的有效性.     

Post-buckling and nonlinear free vibration analysis of geometrically imperfect functionally graded beams resting on nonlinear elastic foundation

In this paper, post-buckling and nonlinear vibration analysis of geometrically imperfect beams made of functionally graded materials (FGMs) resting on nonlinear elastic foundation subjected to axial force are studied. The material properties of FGMs are assumed to be graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents. The assumptions of a small strain and moderate deformation are used. Based on Euler–Bernoulli beam theory and von-Karman geometric nonlinearity, the integral partial differential equation of motion is derived. Then this partial differential equation (PDE) problem, which has quadratic and cubic nonlinearities, is simplified into an ordinary differential equation (ODE) problem by using the Galerkin method. Finally, the governing equation is solved analytically using the variational iteration method (VIM). Some new results for the nonlinear natural frequencies and buckling load of the imperfect functionally graded (FG) beams such as the effects of vibration amplitude, elastic coefficients of foundation, axial force, end supports and material inhomogeneity are presented for future references. Results show that the imperfection has a significant effect on the post-buckling and vibration response of FG beams.     

基于Hartley变换的三维真振幅偏移算子研究

在实际应用中，以快速Fourier变换为基础的偏移方法，将本来是实数的地震道转化为复数参加运算，导致了计算机内存的增加。本文把只有纯实数运算的Hartley变换引入到基于Fourier变换的偏移算法，再利用三维真振幅偏移单程波方程，结合Fourier变换与Hartley变换的内在关系，经过数学推理，具体导出了裂步Hartley变换真振幅偏移算子。与一般裂步Fourier法相比，裂步Hartley变换真振幅偏移算法既提高了计算效率又对球面扩散问题进行了振幅补偿。     

Homoclinic Orbit for the Cubic–Quintic Ginzburg–Landau Equation

The Ginzburg–Landau equation with small complex coefficients is discussed in this paper. A transformation is introduced to change the equation into a three order, ordinary differential equation and the existence of the homoclinic orbit for this system has been proved by analytical methods.     

Transverse vibration of viscoelastic sandwich beam with time-dependent axial tension and axially varying moving velocity

Nonlinearly parametric resonances of axially accelerating moving viscoelastic sandwich beams with time-dependent tension are investigated in this paper. Based on the Kelvin differential constitutive equation, the controlling equation of the transverse vibration of a beam with large deflection is established. The system has been subjected to a time varying velocity and a harmonic axial tension. Here the governing equation of motion contains linear parametric terms and two frequencies, one is the frequency of axially moving velocity and the other one is the frequency of varying tension. The method of multiple scales is applied directly to the governing equation to obtain the complex eigenfunctions and natural frequencies of the system. The elimination of secular terms leads to the steady-state response and amplitude of vibrations. The influence of various parameters such as initial tension on natural frequencies and the amplitude of axial fluctuation, the phase angle between the two frequencies on response curves has been investigated for two different resonance conditions. With the help of numerical results, it has been shown that by using suitable initial tension, the amplitude of axial fluctuation, the phase angle, the vibration of the sandwich beam can be significantly controlled.     

时滞Duffing型方程周期解的存在唯一性

该文利用拓扑度方法研究了一类时滞依赖状态的广义Duffing型泛函微分方程x'(t)$ 该文利用拓扑度方法研究了一类时滞依赖状态的广义Duffing型泛函微分方程x'(t)$ 该文利用拓扑度方法研究了一类时滞依赖状态的广义Duffing型泛函微分方程x'(t) g(x(t-τ(t,x(t))))=f(t)周期解的存在性,得到了方程周期解存在的充分条件和必要条件．研究了当滞量为常值时,方程周期解的存在唯一性．并且给出了所研究问题的一个应用实例．     

A random differential transform method: Theory and applications

In this article, a Differential Transform Method (DTM) based on the mean fourth calculus is developed to solve random differential equations. An analytical mean fourth convergent series solution is found for a nonlinear random Riccati differential equation by using the random DTM. Besides obtaining the series solution of the Riccati equation, we provide approximations of the main statistical functions of the stochastic solution process such as the mean and variance. These approximations are compared to those obtained by the Euler and Monte Carlo methods. It is shown that this method applied to the random Riccati differential equation is more efficient than the two above mentioned methods.     

如何不用复数讲解常系数线性微分方程

在高等数学课程中,复指数函数及其导数知识的严格讲解,通常要比微分方程知识的讲解晚很多.这使得微分方程的教学在逻辑上有些不足.用复值函数解的复系数线性组合推导出实值函数解,在教学实践中,学生经常感到迷惑.不以复数的任何知识作为前提,给出了常系数微分方程的一种自然的讲解方法.