Small-amplitude limit cycles of polynomial Liénard systems

In this paper,we study the number of limit cycles appeared in Hopf bifurcations of a Linard system with multiple parameters.As an application to some polynomial Li’enard systems of the form x=y,y=gm(x)-fn(x)y,we obtain a new lower bound of maximal number of limit cycles which appear in Hopf bifurcation for arbitrary degrees m and n.     

A class of hybrid formulae for the numerical integration of stiff systems

A class of hybrid formulae suitable for the numerical integration of stiff systems of first order ordinary differential equations is presented. These formulae are defined in a predictor-corrector mode and contain an off-step pointx _n+v, which is distinct from the points defined by the step used. The parameterv is examined in an attempt to obtain better stability properties and higher order formulae of one step only. We obtain order five and eitherA-stability or stiff-stability according to the values ofv used. Some numerical results are presented.     

Elliptic equation with a singular potential in a domain with a conic point

This paper deals with the behavior of the nonnegative solutions of the problem $$- \Delta u = V(x)u, \left. u \right|\partial \Omega = \varphi (x)$$ in a conical domain Ω ? ? ⁿ , n ≥ 3, where 0 ≤ V (x) ∈ L₁(Ω), 0 ≤ ?(x) ∈ L₁(?Ω) and ?(x) is continuous on the boundary ?Ω. It is proved that there exists a constant C _*(n) = (n ? 2)²/4 such that if V ₀(x) = (c + λ ₁)|x|^?2, then, for 0 ≤ c ≤ C _*(n) and V(x) ≤ V ₀(x) in the domain Ω, this problem has a nonnegative solution for any nonnegative boundary function ?(x) ∈ L ₁(?Ω); for c > C _*(n) and V(x) ≥ V ₀(x) in Ω, this problem has no nonnegative solutions if ?(x) > 0.     

Normal Vectors on Manifolds of Critical Points for Parametric Robustness of Equilibrium Solutions of ODE Systems

Summary. {Equilibrium solutions of systems of parameterized ordinary differential equations \dot x = f(x, α) , x ∈ R ⁿ , α∈ R ^m can be characterized by their parametric distance to manifolds of critical solutions at which the behavior of the system changes qualitatively. Critical points of interest are bifurcation points and points at which state variable constraints or output constraints are violated. We use normal vectors on manifolds of critical points to measure the distance between these manifolds and equilibrium solutions as suggested in I. Dobson [J. Nonlinear Sci., 3:307-327, 1993], where systems of equations to calculate normal vectors on codimension-1 bifurcations were presented. We present a scheme to derive systems of equations to calculate normal vectors on manifolds of critical points which (i) generalizes to bifurcations of arbitrary codimension, (ii) can be applied to state variable constraints and output constraints, (iii) implies that the normal vector defining system of equations is of size c ₁ n+ c ₂ m+ c ₃ , c _i ∈ R , i.e., no bilinear terms nm or higher-order terms occur, (iv) reduces the number of equations for normal vectors on Hopf bifurcation manifolds compared to previous work, and (v) simplifies the proof of regularity of the normal vector system. As an application of this scheme, we present systems of equations for normal vectors to manifolds of output/state variable constraints, to manifolds of saddle-node, Hopf, cusp, and isola bifurcations, and we give illustrative examples of their use in engineering applications.} Received September 27, 2000; accepted December 10, 2001 Online publication March 11, 2002 Communicated by Y. G. Kevrekidis Communicated by Y. G. Kevrekidis rid="     

Mixed quadratic-cubic autocatalytic reaction–diffusion equations: Semi-analytical solutions

Semi-analytical solutions for autocatalytic reactions with mixed quadratic and cubic terms are considered. The kinetic model is combined with diffusion and considered in a one-dimensional reactor. The spatial structure of the reactant and autocatalyst concentrations are approximated by trial functions and averaging is used to obtain a lower-order ordinary differential equation model, as an approximation to the governing partial differential equations. This allows semi-analytical results to be obtained for the reaction–diffusion cell, using theoretical methods developed for ordinary differential equations. Singularity theory is used to investigate the static multiplicity of the system and obtain a parameter map, in which the different types of steady-state bifurcation diagrams occur. Hopf bifurcations are also found by a local stability analysis of the semi-analytical model. The transitions in the number and types of bifurcation diagrams and the changes to the parameter regions, in which Hopf bifurcations occur, as the relative importance of the cubic and quadratic terms vary, is explored in great detail. A key outcome of the study is that the static and dynamic stability of the mixed system exhibits more complexity than either the cubic or quadratic autocatalytic systems alone. In addition it is found that varying the diffusivity ratio, of the reactant and autocatalyst, causes dramatic changes to the dynamic stability. The semi-analytical results are show to be highly accurate, in comparison to numerical solutions of the governing partial differential equations.     

Estimation of states and parameters in continuous non-linear systems with discrete observations

We consider a class of continuous non-linear systems defined by the ordinary differential equation x = f(x, t) + g(x, t)u, where u is an unknown input representing noise or disturbances. The object is to estimate states and parameters in these systems by means of a fixed number of discrete observations y_i = h(x(t_i), t_i) + v_i, 1 ? i ? m, where the v_i represents unknown errors in the measurements y_i. No statistical assumptions are made concerning the nature of the unknown input u or the unknown measurement errors v_i. A weighted least squares criterion is defined as a measure of the optimal estimate. A result concerning the existence of solutions of the differential equation which minimize the criterion is presented. The necessary conditions for an optimal estimate, a set of Euler-Lagrange equations and multi-point discontinuous non-linear boundary conditions, are given. The multi-point problem is converted to an equivalent continuous two-point boundary value problem of larger dimension in the case in which the observations are assumed to be linear functions of the state. A pair of equivalent quasilinearization algorithms is defined for the two-point system and the multi-point system. Quadratic convergence for these algorithms is proved.     

Elementary approximations for zeros of Bessel functions

Let j_vk, y_vk and c_vk denote the kth positive zeros of the Bessel functions J_v(x), Y_v(x) and of the general cylinder function C_v(x) = cos αJ_v(x)?sin αY_v(x), 0 ? α < π, respectively. In this paper we extend to c_vk, k = 2, 3,..., some linear inequalities presently known only for j_vk. In the case of the zeros y_vk we are able to extend these inequalities also to k = 1. Finally in the case of the first positive zero j_v1 we compare the linear enequalities given in [9] with some other known inequalities.     

Two new approaches to computing Hopf bifurcation problems

We consider the computation of Hopf bifurcation for ordinary differential equations. Two new extended systems are given for the calculation of Hopf bifurcation problems: the first is composed of differential-algebraic equations with index 1, the other consists of differential equations by using a symmetry inherited from the autonomous system of ordinary differential equations. Both methods are especially suitable for calculating bifurcating periodic solutions since they transform the Hopf bifurcation problem into regular nonlinear boundary value problems which are very easy to implement. The bifurcation solutions become isolated solutions of the extended system so that our methods work both in the subcritical and supercritical case. The extended systems are based on an additional parameter ε; practical experience shows that one gets convergence for ε sufficiently large so that a substantial part of the bifurcating branch can be computed. The two methods are illustrated by numerical examples and compared with other procedures.     

Algorithm for the computation of Bessel function integrals

A fortran subroutine is given for the computation of integrals of the form ∫^c₀f(x)J_v(αx)dx, where v = 0, 1,…,10.     

Spatially nonhomogeneous equilibrium in a reaction-diffusion system with distributed delay

We consider a reaction-diffusion system with general time-delayed growth rate and kernel functions. The existence and stability of the positive spatially nonhomogeneous steady-state solution are obtained. Moreover, taking minimal time delay τ as the bifurcation parameter, Hopf bifurcation near the steady-state solution is proved to occur at a critical value τ=τ₀. Especially, the Hopf bifurcation is forward and the bifurcated periodic solutions are stable on the center manifold. The general results are applied to competitive and cooperative systems with weak or strong kernel function respectively.     

Elementary divisors of induced transformations on symmetry classes of tensors

Let V_χ(G) denote the symmetry class of tensors over the vector space V associated with the permutation group G and irreducible character χ. Write v₁*v₂*...*v_m for the decomposable symmetrized product of the indicated vectors (m=degG). If T is a linear operator on V, let K(T) denote the associated operator on V_χ(G), i.e., K(T)v₁*v₂*...*v_m=Tv₁*Tv₂*...*Tv_m. Denote by D(T) the derivation operator D(T)v₁*v₂*...*v_m=Tv₁*v₂...*v_m+v₁*Tv₂*v₃* ...*v_m+...+v₁*v₂*...*v_m–1*Tv_m. The article concerns the elementary divisors of K(T) and D(T).     

Global stability of traveling wavefronts for nonlocal reaction-diffusion equations with time delay

This paper is concerned with the stability of traveling wavefronts for a population dynamics model with time delay. Combining the weighted energy method and the comparison principle, the global exponential stability of noncritical traveling wavefronts(waves with speeds c c_*, where c = c~* is the minimal speed) is established, when the initial perturbations around the wavefront decays to zero exponentially in space as x →-∞, but it can be allowed arbitrary large in other locations, which improves the results in [9, 18, 21].     

On the behavior of the oscillatory solutions of first or second order delay differential equations

Some new results are given concerning the behavior of the oscillatory solutions of first or second order delay differential equations. These results establish that all oscillatory solutions x of a first or second order delay differential equation satisfy x(t)=O(v(t)) as t→∞, where v is a nonoscillatory solution of a corresponding first or second order linear delay differential equation. Some applications of the results obtained are also presented.     

Existence and continuity of an implicit function in the neighborhood of an abnormal point

Solutions to the equation F(x, ??) = 0 with unknown x and the parameter ?? in the neighborhood of the solution (x _*, ??_*) under the additional constraint x ?? U, where U is a closed convex set, are studied. The sufficient conditions for existence of an implicit function without prior assumption of the normalcy of point x _* are given. The obtained result is used to investigate the local solvability of controlled systems with mixed constraints.     

Hopf bifurcation in symmetric configuration of predator-prey-mutualist systems

In this paper we apply the equivariant degree method to study Hopf bifurcations in a system of differential equations describing a symmetric predator-prey-mutualist model with diffusive migration between interacting communities. A topological classification (according to symmetry types), of symmetric Hopf bifurcation in configurations of populations with D₈, D₁₂, A₄ and S₄ symmetries, is presented with estimation on minimal number of bifurcating branches of periodic solutions.     

The ring of an outer von Neumann frame in modular lattices

We prove the following theorem. Let (a ₁, . . . , a _m , c ₁₂, . . . , c _1m ) be a spanning von Neumann m-frame of a modular lattice L, and let (u ₁, . . . , u _n , v ₁₂, . . . , v _1n ) be a spanning von Neumann n-frame of the interval [0, a ₁]. Assume that either m ≥ 4, or L is Arguesian and m ≥ 3. Let R* denote the coordinate ring of (a ₁, . . . , a _m , c ₁₂, . . . , c _1m ). If n ≥ 2, then there is a ring S* such that R* is isomorphic to the ring of all n × n matrices over S*. If n ≥ 4 or L is Arguesian and n ≥ 3, then we can choose S* as the coordinate ring of (u ₁, . . . , u _n , v ₁₂, . . . , v _1n ).     

Complex hopf bifurcation of integral surfaces

The existence of the Hopf bifurcation of a complex ordinary differential equation system in the complex domain is studied in this paper by using the complex qualitative theory. In the complex domain, we conclude that the Hopf bifurcation appears for both directions of the parameter. The formulae of the Hopf bifurcation are also given in this paper.     

Hopf bifurcation and topological horseshoe of a novel finance chaotic system

This paper deals with the existence of both Hopf bifurcation and topological horseshoe for a novel finance chaotic system. First, through rigorous mathematical analysis, we show that a Hopf bifurcation occurs at systems’ three equilibriums S_0,1,2 and Hopf bifurcation at equilibrium S₀ is non-degenerate and supercritical. Second, the computer-assisted verifications for horseshoe chaos in the system are given. Simulation results are presented to support the analysis.     

Error-free computer solution of certain systems of linear equations

In this article we propose a procedure which generates the exact solution for the system Ax = b, where A is an integral nonsingular matrix and b is an integral vector, by improving the initial floating-point approximation to the solution. This procedure, based on an easily programmed method proposed by Aberth [1], first computes the approximate floating-point solution x^* by using an available linear equation solving algorithm. Then it extracts the exact solution x from x^* if the error in the approximation x^* is sufficiently small. An a posteriori upper bound for the error of x^* is derived when Gaussian Elimination with partial pivoting is used. Also, a computable upper bound for |det(A)|, which is an alternative to using Hadamard's inequality, is obtained as a byproduct of the Gaussian Elimination process.