Numerical methods for Hamilton Jacobi functional differential equations

Initial and initial boundary value problems for first order partial functional differential equations are considered. Explicit difference schemes of the Euler type and implicit difference methods are investigated. The following theoretical aspects of the methods are presented. Sufficient conditions for the convergence of approximate solutions are given and comparisons of the methods are presented. It is proved that assumptions on the regularity of given functions are the same for both the methods. It is shown that conditions on the mesh for explicit difference schemes are more restrictive than suitable assumptions for implicit methods. There are implicit difference schemes which are convergent and corresponding explicit difference methods are not convergent. Error estimates for both the methods are construted.     

On Time Discretizations for Spectral Methods

New methods are introduced for the time integration of the Fourier and Chebyshev methods of solution for dynamic differential equations. These methods are unconditionally stable, even though no matrix inversions are required. Time steps are chosen by accuracy requirements alone. For the Fourier method both leapfrog and Runge-Kutta methods are considered. For the Chebyshev method only Runge-Kutta schemes are tested. Numerical calculations are presented to verify the analytic results. Applications to the shallow water equations are presented.     

State models and asymptotic behaviour of 2D positive systems

Homogeneous 2D positive systems are 2D state-space models whosevariables are alwalys nonnegative and, consequently, are describedby a pair of nonnegative square matrices (A, B). In the paper,the properties of these pairs are discussed both in the generalcase and under particular assumptions like finite memory, separability,and property L. Various aspects of the positive asymptotic dynamic are considered;in particular, sufficient conditions are provided guaranteeingthat the local states are eventually strictly positive. Finally,some results on the convergence of the states towards a constantasymptotic distribution are presented.     

Depth notions for orthogonal regression

Global depth, tangent depth and simplicial depths for classical and orthogonal regression are compared in examples, and properties that are useful for calculations are derived. The robustness of the maximum simplicial depth estimates is shown in examples. Algorithms for the calculation of depths for orthogonal regression are proposed, and tests for multiple regression are transferred to orthogonal regression. These tests are distribution free in the case of bivariate observations. For a particular test problem, the powers of tests that are based on simplicial depth and tangent depth are compared by simulations.     

Estimating a survival function with incomplete cause-of-death data

We propose a random censorship model which permits uncertainty in the cause of death assessments for a subset of the subjects in a survival experiment. A nonparametric maximum likelihood approach and a “self-consistency” approach are considered. The solution sets corresponding to both approaches are found. They are infinite and identical. Only some of the solutions are consistent; i.e., the MLEs and self-consistent estimators are not consistent in general. Two estimates are thus proposed and their asymptotic properties are studied. It is shown that both estimates are strongly consistent and converge to Gaussian processes. The covariance structures of these Gaussian processes are derived.     

Some properties of singular integral operators over an open curve Dedicated to the 85th birthday of Prof.Lu Jinake

In this paper, the stability properties, the endpoint behavior and the invertible relations of Cauchy-type singular integral operators over an open curve are discussed. If the endpoints of the curve are not special, this type of operators are proved to be stable. At the endpoints, either the singularity or smoothness of the operators are exactly described. And the function sets or spaces on which the operators are invertible as well as the corresponding inverted operators are given. Meanwhile, some applications for the solution of Cauchy-type singular integral equations are illustrated.     

Some properties of singular integral operators over an open curve Dedicated to the 85th birthday of Prof. Lu Jianke

In this paper, the stability properties, the endpoint behavior and the invertible relations of Cauchy-type singular integral operators over an open curve are discussed. If the endpoints of the curve are not special, this type of operators are proved to be stable. At the endpoints, either the singularity or smoothness of the operators are exactly described. And the function sets or spaces on which the operators are invertible as well as the corresponding inverted operators are given. Meanwhile, some applications for the solution of Cauchy-type singular integral equations are illustrated.     

Three types of matrix stability

Three types of stability of real matrices are compared and necessary conditions are obtained in terms of the principal submatrices. For normal matrices and matrices whose off-diagonal elements are all positive, these conditions are sufficient, and the three types of stability are all equivalent. Necessary and sufficient conditions in terms of the elements of the matrix are proven for matrices of order 2 and 3.     

Periodic and solitary-wave solutions of the Degasperis–Procesi equation

Travelling-wave solutions of the Degasperis–Procesi equation are investigated. The solutions are characterized by two parameters. For propagation in the positive x-direction, hump-like, inverted loop-like and coshoidal periodic-wave solutions are found; hump-like, inverted loop-like and peakon solitary-wave solutions are obtained as well. For propagation in the negative x-direction, there are solutions which are just the mirror image in the x-axis of the aforementioned solutions. A transformed version of the Degasperis–Procesi equation, which is a generalization of the Vakhnenko equation, is also considered. For propagation in the positive x-direction, hump-like, loop-like, inverted loop-like, bell-like and coshoidal periodic-wave solutions are found; loop-like, inverted loop-like and kink-like solitary-wave solutions are obtained as well. For propagation in the negative x-direction, well-like and inverted coshoidal periodic-wave solutions are found; well-like and inverted peakon solitary-wave solutions are obtained as well. In an appropriate limit, the previously known solutions of the Vakhnenko equation are recovered.     

Some properties of complex matrix-variate generalized Dirichlet integrals

Dirichlet integrals and the associated Dirichlet statistical densities are widely used in various areas. Generalizations of Dirichlet integrals and Dirichlet models to matrix-variate cases, when the matrices are real symmetric positive definite or hermitian positive definite, are available [4]. Real scalar variables case of the Dirichlet models are generalized in various directions. One such generalization of the type-2 or inverted Dirichlet is looked into in this article. Matrix-variate analogue, when the matrices are hermitian positive definite, are worked out along with some properties which are mathematically and statistically interesting.     

On the boundedness and unboundedness of external polyhedral estimates for reachable sets of linear differential systems

The boundedness and unboundedness properties of external polyhedral (paralle-lepiped-valued) estimates are investigated for reachable sets of linear differential systems with a stable matrix. Boundedness and unboundedness criteria on an infinite time interval are presented for two types of estimates (“touching” estimates, which were introduced earlier, and estimates with constant orientation matrix). Conditions for the system matrix and bounding sets are given under which there are bounded estimates among the estimates of the mentioned types, under which there are unbounded estimates, and under which all the estimates are bounded or all the estimates are unbounded. In terms of the exponents of the estimates, the possible rate of their growth is described. For two-dimensional systems, the classification and comparison of possible situations of the boundedness or unboundedness for estimates of both types are given and boundedness criteria for estimates with special (orthogonal and “quasi-orthogonal”) constant orientation matrices are found. Results of numerical modeling are presented.     

Relations Between Linear Equations and Painlevé’s Equations

Painlevé equations are studied on the basis of linear equations, which are generic for them. Different possible approaches are compared to each other. Formulas binding these approaches are derived. Symmetries demonstrated in the equations are also a subject of discussion.     

Parabolic flow over a flat plate with wave disturbance in the main stream

In this paper the effects of unsteady flow past a flat plate, which is at rest, are studied when the steady parabolic flow outside the boundary layer is in a wave disturbance. Solutions are developed for large and small times. The behaviour of the amplitudes and phase leads of skin friction components along chordwise and spanwise directions are studied in detail. Two cases, one when the main stream lines are straight and yawed, and the other when the main stream lines are parabolic and concave with respect to the points on the chordwise direction are studied separately and the results are interpreted graphically.     

关于Toeplitz+Hankel线性方程组的迭代解法

本文研究Toeplitz+Hankel线性方程组的预处理迭代解法.我们提出了几个新的预条件子,并分析了预处理矩阵的谱性质,当生成函数在Wiener类中时,预处理矩阵的特征值聚集在1附近.数值实验表明该预处理子比文[5]中的预处理子更有效.     

New perturbation bounds for denumerable Markov chains

This paper is devoted to perturbation analysis of denumerable Markov chains. Bounds are provided for the deviation between the stationary distribution of the perturbed and nominal chain, where the bounds are given by the weighted supremum norm. In addition, bounds for the perturbed stationary probabilities are established. Furthermore, bounds on the norm of the asymptotic decomposition of the perturbed stationary distribution are provided, where the bounds are expressed in terms of the norm of the ergodicity coefficient, or the norm of a special residual matrix. Refinements of our bounds for Doeblin Markov chains are considered as well. Our results are illustrated with a number of examples.     

Potential-based numerical solution of Dirichlet problems for the Helmholtz equation

Three-dimensional Dirichlet problems for the Helmholtz equation are considered in generalized formulations. By applying single-layer potentials, they are reduced to Fredholm boundary integral equations of the first kind. The equations are discretized using a special averaging method for integral operators with weak singularities in the kernels. As a result, the integral equations are approximated by systems of linear algebraic equations with easy-to-compute coefficients, which are solved numerically by applying the generalized minimal residual method. A modification of the method is proposed that yields solutions in the spectra of interior Dirichlet problems and integral operators when the integral equations are not equivalent to the original differential problems and are not well-posed. Numerical results are presented for assessing the capabilities of the approach.     

Some properties of singular integral operators over an open curve

In this paper, the stability properties, the endpoint behavior and the invertible relations of Cauchy-type singular integral operators over an open curve are discussed. If the endpoints of the curve are not special, this type of operators are proved to be stable. At the endpoints, either the singularity or smoothness of the operators are exactly described. And the function sets or spaces on which the operators are invertible as well as the corresponding inverted operators are given. Meanwhile, some applications for the solution of Cauchy-type singular integral equations are illustrated. This work was partially supported by the National Natural Science Foundation of China (Grant No. 10471048)     

Dynamics of a Stewart platform

The kinematics and dynamics of a Stewart platform supported by six pneumocylinders are considered. The differential equations of motion are written, and the forces providing the fulfillment of the formulated law of motion are calculated. The inertia and weight of pneumocylinders are introduced into consideration to refine the equations of motion. The obtained equations are used to study the motion of a loaded Stewart platform, providing the stability of this motion by means of feedback control. Some numerical examples are given.     

Multiple Attractors in Juvenile-adult Single Species Models

The role of age-structure and the Allee effect in generating multiple attractors in juvenile-adult single species single patch discrete-time models without dispersal are studied. In the presence of the Allee effect juvenile-adult single patch models support multiple attractors. However, in the absence of the Allee effect single attractors are supported when the dynamics are compensatory while multiple attractors are supported under overcompensatory dynamics. When the governing dynamics are compensatory, the boundaries of the basins of attraction have simple structure while complicated fractal basin boundaries are supported under overcompensatory dynamics.     

A classification of vector harmonizable processes

Hilbert space valued second order stochastic processes over the real line are considered. Various hmonizabilities and V-boundednesses are introduced and their interrelations are obtained as well as heir integral representations. Examples are given to distinguish most of the harmonizabilities. Stationary dilations of harmonizable processes are also discussed. Finally, for some harmonizabilities, some convergence theorems for sequences of processes are obtained