血红细胞时滞模型的谱及其解的结构

研究一个带有时滞的血红细胞模型的解展开问题.对模型在平衡点处线性化,并利用泛函分析方法,将线性化模型写成抽象发展方程.借助半群理论证明了方程的适定性.对系统算子细致的谱分析,得到了本征值的渐近表达式.通过对算子的Riesz谱投影范数的渐近估计,证明系统的本征向量不能构成状态空间的基,但我们仍给出了方程的解在平衡点附近按照本征向量的的渐近展开.     

带时滞机器人模型解的展开

研究一个带有时滞机器人模型解的性质, 其中机器人的动力学行为由一组含有时滞的微分方程描述.通过引入Hilbert状态空间将其写成一个发展方程, 利用半群理论得到抽象发展方程的适定性. 通过对系统算子的谱作细致分析,得到谱的渐近表达式, 并证明系统的本征函数并不构成空间的基.但我们证明了对时间大于5倍时滞时, 方程的解按照本征向量的展开式.     

Riesz basis property of root vectors of non‐self‐adjoint operators generated by aircraft wing model in subsonic airflow

This paper is the third in a series of several works devoted to the asymptotic and spectral analysis of a model of an aircraft wing in a subsonic air flow. This model has been developed in the Flight Systems Research Center of UCLA and is presented in the works by Balakrishnan. The model is governed by a system of two coupled integro‐differential equations and a two‐parameter family of boundary conditions modeling the action of the self‐straining actuators. The differential parts of the above equations form a coupled linear hyperbolic system; the integral parts are of the convolution type. The system of equations of motion is equivalent to a single operator evolution–convolution equation in the energy space. The Laplace transform of the solution of this equation can be represented in terms of the so‐called generalized resolvent operator, which is an operator‐valued function of the spectral parameter. This generalized resolvent operator is a finite‐meromorphic function on the complex plane having the branch cut along the negative real semi‐axis. Its poles are precisely the aeroelastic modes and the residues at these poles are the projectors on the generalized eigenspaces. In the first two papers (see [33, 34]) and in the present one, our main object of interest is the dynamics generator of the differential parts of the system. This generator is a non‐self‐adjoint operator in the energy space with a purely discrete spectrum. In the first paper, we have shown that the spectrum consists of two branches, and have derived their precise spectral asymptotics with respect to the eigenvalue number. In the second paper, we have derived the asymptotical approximations for the mode shapes. Based on the asymptotical results of the first two papers, in the present paper, we (a) prove that the set of the generalized eigenvectors of the aforementioned differential operator is complete in the energy space; (b) construct the set of vectors which is biorthogonal to the set of the generalized eigenvectors in the case when there might be not only eigenvectors but associate vectors as well; and (c) prove that the set of the generalized eigenvectors forms a Riesz basis in the energy space. To prove the main result of the paper, we made use of the Nagy–Foias functional model for non‐self‐adjoint operators. The results of all three papers will be important for the reconstruction of the solution of the original initial‐boundary‐value problem from its Laplace transform in the forthcoming papers. Copyright © 2000 John Wiley & Sons, Ltd.     

Multidimensional parametrization and numerical solution of systems of nonlinear equations

The numerical solution of a system of nonlinear algebraic or transcendental equations is examined within the framework of the parameter continuation method. An earlier result of the author according to which the best parameters should be sought in the tangent space of the solution set of this system is now refined to show that the directions of the eigenvectors of a certain linear self-adjoint operator should be used for finding these parameters. These directions correspond to the extremal values of the quadratic form associated with the above operator. The parametric approximation of curves and surfaces is considered.     

System of nonsmooth variational inequalities with applications

In this paper, we consider a system of vector variational inequalities and a system of nonsmooth variational inequalities defined by means of Clarke directional derivative. We also consider the Nash equilibrium problem with vector pay-offs and its scalarized form. We present some relations among these systems and problems. The existence results for a solution of system of nonsmooth variational inequalities are given. As a consequence, we derive an existence result for a solution of Nash equilibrium problem with vector pay-offs.     

Quadratic model updating with no spill-over and incomplete measured data: Existence and computation of solution

Quadratic finite element model updating problem (QFEMUP), to be studied in this paper, is concerned with updating a symmetric nonsingular quadratic pencil in such a way that, a small set of measured eigenvalues and eigenvectors is reproduced by the updated model. If in addition, the updated model preserves the large number of unupdated eigenpairs of the original model, the model is said to be updated with no spill-over. QFEMUP is, in general, a difficult and computationally challenging problem due to the practical constraint that only a very small number of eigenvalues and eigenvectors of the associated quadratic eigenvalue problem are available from computation or measurement. Additionally, for practical effectiveness, engineering concerns such as nonorthogonality and incompleteness of the measured eigenvectors must be considered. Most of the existing methods, including those used in industrial settings, deal with updating a linear model only, ignoring damping. Only in the last few years a small number of papers been published on the quadratic model updating; several of the above issues have been dealt with both from theoretical and computational point of views. However, mathematical criterion for existence of solution has not been fully developed. In this paper, we first (i) prove a set of necessary and sufficient conditions for the existence of a solution of the no spill-over QFEMUP, then (ii) present a parametric representation of the solution, assuming a solution exists and finally, (iii) propose an algorithm for QFEMUP with no spill-over and incomplete measured eigenvectors. Interestingly, it is shown that the parametric representation can be constructed with the knowledge of only the few eigenvalues and eigenvectors that are to be updated and the corresponding measured eigenvalues and eigenvectors—complete knowledge of eigenvalues and eigenvectors of the original pencil is not needed, which makes the solution readily applicable to real-life structures.     

二阶耗散动力系统的降维对解长期行为的误差估计

基于非线性动力学理论,对一类高维二阶耗散自治动力系统的降维及其对解的长期行为的影响进行了理论分析．该分析将方程的解投影到控制方程的线性算子的特征向量所张成的完备空间中,并在相空间中引入一距离的概念,方便地解决了缩减后系统与原始系统解之间的误差或距离的描述．基于此距离定义,首先,分析了由于高阶模态的截取对解的长期行为的影响,并推导出了相应的误差估计,该估计表明由于降维对系统长期行为的影响不仅与系统的高阶子空间中的固有频率和阻尼比乘积的最小值有关,并且与高阶子空间中的某一最大固有频率有关．然后,将一般的模态截取视为对原系统的解的一个扰动,对一些文献中由于降维程度的不同而造成解的拓扑性质发生变化的现象进行了定性的解释．     

Gap Functions and Lyapunov Functions

Equilibrium problems play a central role in the study of complex and competitive systems. Many variational formulations of these problems have been presented in these years. So, variational inequalities are very useful tools for the study of equilibrium solutions and their stability. More recently a dynamical model of equilibrium problems based on projection operators was proposed. It is designated as globally projected dynamical system (GPDS). The equilibrium points of this system are the solutions to the associated variational inequality (VI) problem. A very popular approach for finding solution of these VI and for studying its stability consists in introducing the so-called "gap-functions", while stability analysis of an equilibrium point of dynamical systems can be made by means of Lyapunov functions. In this paper we show strict relationships between gap functions and Lyapunov functions.     

Extension of the solution of nonlinear equations in the neighborhood of a bifurcation point

Using the method of continuous extension with respect to a parameter we develop a method of constructing the load trajectory of a structure having both limit points and bifurcation points. The method is applicable for the systems of nonlinear algebraic equations that describe the family of extremals that minimize the value of the total potential strain energy of the structure, and makes it possible to find all the branches of the load trajectory emanating from a bifurcation point and extend the solution along any of them. The method is based on the fact that the eigenvectors of the augmented Jacobian of the system of equations in the extended space of variables that correspond to zero eigenvalues on the main branch of the load trajectory are bifurcation vectors and form the active subspace of solutions of the equations of the extension. Meanwhile the other eigenvectors form the passive subspace that contains the extension vector with respect to the main branch of the load. As a result the entire process of computing the extension vector of the solution at any point of the load trajectory reduces to determining the eigenvectors of the augmented Jacobian of the original system of nonlinear algebraic equations, identifying them according as they belong to the active or passive subspace, and forming the extension vector of the solution using them and analytic relations Translated fromMatematychni Metody ta Fizyko-Mekhanichni Polya, Vol. 41, No. 1, 1998, pp. 35–46.     

Application of finite difference method of lines on the heat equation

In this article, we apply the method of lines (MOL) for solving the heat equation. The use of MOL yields a system of first–order differential equations with initial value. The solution of this system could be obtained in the form of exponential matrix function. Two approaches could be applied on this problem. The first approach is approximation of the exponential matrix by Taylor expansion, Padé and limit approximations. Using this approach leads to create various explicit and implicit finite difference methods with different stability region and order of accuracy up to six for space and superlinear convergence for time variables. Also, the second approach is a direct method which computes the exponential matrix by applying its eigenvalues and eigenvectors analytically. The direct approach has been applied on one, two and three‐dimensional heat equations with Dirichlet, Neumann, Robin and periodic boundary conditions.     

Analytic solution to the problem of the temperature jump in a gas with rotational degrees of freedom

An exact solution is obtained for the first time for the problem of the temperature jump in a gas with allowance for internal (rotational) degrees of freedom. The treatment is based on a model collision integral proposed by the authors. The problem reduces to the solution of a boundary-value problem for a linear vector transport equation with matrix kernel. Separation of the variable leads to a characteristic equation for which eigenvectors are found in a space of generalized functions and the eigenvalue spectrum is investigated. An expansion of the solution to the problem with respect to eigenvectors of the continuous and discrete spectra is established. On the basis of the conditions of solvability of the vector Riemann-Hilbert boundary-value problem which arises in the process of the proof, an exact (in closed form) expression is obtained for the temperature jump.Moscow Pedagogical University. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 95, No. 3, pp. 530–540, June, 1993.     

用奇异值分解方法计算具有重特征值矩阵的特征矢量

若当(Jordan)形是矩阵在相似条件下的一个标准形,在代数理论及其工程应用中都具有十分重要的意义．针对具有重特征值的矩阵,提出了一种运用奇异值分解方法计算它的特征矢量及若当形的算法．大量数值例子的计算结果表明,该算法在求解具有重特征值的矩阵的特征矢量及若当形上效果良好,优于商用软件MATLAB和MATHEMATICA．     

Limit periodic travelling wave solution of a model for biological invasions

In this paper, we consider a class of biological invasion model with density-dependent migrations and Allee effect, which is reduced to one ordinary differential form via the travelling wave solution ansatz. For the corresponding planar system, we firstly obtain the first several weak focal values of its one equilibrium by computing the singular point quantities, then determine the existence of one stable limit cycle from its Hopf bifurcation. Thus a special periodic travelling wave solution which is isolate as a limit is obtained, and it corresponds to the particular real patterns of spread during biological invasions, which is an interesting discovery.     

Stochastic processes and special functions: On the probabilistic origin of some positive kernels associated with classical orthogonal polynomials

We describe and investigate a class of Markovian models based on a form of “dynamic occupancy problem” originating in statistical mechanics. The most fundamental of these gives rise to a transition-probability matrix over (N + 1) discrete states, which proves to have the Hahn polynomials as eigenvectors. The structure of this matrix, which is a convolution of two negative hypergeometric distributions, leads to a factorization into finite-difference sumoperators having forms analogous to the Erdelyi-Kober operators for the continuous variable. These make possible the exact solution of the corresponding eigenvalue problem and hence the spectral representation of the transition matrix. By taking suitable limits, further families of Markov processes can be generated having other classical polynomials as eigenvectors; these, like the polynomials, inherit their properties from the original Hahn system. The Meixner, Jacobi and Laguerre systems arise in this way, having their origin in variants of the basic model. In the last of these cases, the spectral resolution of the continuous transition kernel proves to be identical with Erdelyi's (1938) bilinear formula, which is thus both generalized and given a physical interpretation. Various symmetry and “duality” properties are explored and a number of interesting formulas are obtained as by-products. The use of statistical models to generate kernels, which are thereby guaranteed to be both positive and positive definite, appears to be mathematically fruitful, while the models themselves seem likely to have application to a variety of topics in applied probability.     

A method for the approximate solution of quasi-static problems for hardening bodies

A method for the approximate solution of quasi-static problems for hardening elastoplastic bodies is proposed. The constitutive relation of the model is taken in the form of a variational inequality. An approximate solution of the initial problem is constructed in time steps and, by means of the finite element method, is reduced to the solution of a system of two variational inequalities in corresponding finite-dimensional space. It is shown that the solution of this system is equivalent to finding the saddle point of the corresponding quadratic functional. To find the saddle point, Udzawa's algorithm is used, by means of which the process of finding the velocity vector and stress tensor reduces to the successive calculation of these quantities: the velocity vector is determined from the variational inequality corresponding to the equilibrium equation, and the stress tensor is determined from the variational inequality corresponding to the constitutive relation. The latter inequality is reduced to a certain non-linear equation containing the operation of projection onto a closed convex set corresponding to the elastic strains of the medium. In turn, the solution of the non-linear equation is constructed using the method of successive approximations. To illustrate the use of the proposed method, the one-dimensional problem of the quasi-static deformation of a cylindrical tube under a load applied to its internal surface is considered.     

Backward error of polynomial eigenvalue problems solved by linearization

It is commonplace in many application domains to utilize polynomial eigenvalue problems to model the behaviour of physical systems. Many techniques exist to compute solutions of these polynomial eigenvalue problems. One of the most frequently used techniques is linearization, in which the polynomial eigenvalue problem is turned into an equivalent linear eigenvalue problem with the same eigenvalues, and with easily recoverable eigenvectors. The eigenvalues and eigenvectors of the linearization are usually computed using a backward stable solver such as the QZ algorithm. Such backward stable algorithms ensure that the computed eigenvalues and eigenvectors of the linearization are exactly those of a nearby linear pencil, where the perturbations are bounded in terms of the machine precision and the norms of the matrices defining the linearization. Although we have solved a nearby linear eigenvalue problem, we are not certain that our computed solution is in fact the exact solution of a nearby polynomial eigenvalue problem. Here, we perform a backward error analysis for the solution of a specific linearization for polynomials expressed in the monomial basis. We use a suitable one-sided factorization of the linearization that allows us to map generic perturbations of the linearization onto structured perturbations of the polynomial coefficients. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Exponential and Super Stability of a Wave Network

In this paper, we investigate the spectral distribution and stability of a star-shaped wave network with N edges, of which the feedback gain constants fail to satisfy the assumptions for Riesz basis generation. By a detailed spectral analysis, we present the explicit expressions of the spectra, which consist of simple eigenvalues located on a vertical line in the complex left half-plane. In addition we show that the eigenvectors are not complete in the state space. Further, we decompose the state space into the spectral-subspace and another invariant subspace of infinite dimension, which form a topological direct sum. We prove that, in the spectral-subspace, the solution can be expanded according to the eigenvectors, and hence the solution is exponentially stable; in the other subspace, the associated semigroup is super-stable, i.e., the solution is identical to zero after finite time. In particular, we give the explicit decay rate and the maximum existence time of the nonzero part of the solution.     

Antidiffusive Lagrangian‐remap schemes for models of polydisperse sedimentation

One‐dimensional models of gravity‐driven sedimentation of polydisperse suspensions with particles that belong to N size classes give rise to systems of N strongly coupled, nonlinear first‐order conservation laws for the local solids volume fractions. As the eigenvalues and eigenvectors of the flux Jacobian have no closed algebraic form, characteristic‐wise numerical schemes for these models become involved. Alternative simple schemes for this model directly utilize the velocity functions and are based on splitting the system of conservation laws into two different first‐order quasi‐linear systems, which are solved successively for each time iteration, namely, the Lagrangian and remap steps (so‐called Lagrangian‐remap [LR] schemes). This approach was advanced in (Bürger, Chalons, and Villada, SIAM J Sci Comput 35 (2013), B1341–B1368) for a multiclass Lighthill–Whitham‐Richards traffic model with nonnegative velocities. By incorporating recent antidiffusive techniques for transport equations a new version of these Lagrangian‐antidiffusive remap (L‐AR) schemes for the polydisperse sedimentation model is constructed. These L‐AR schemes are supported by a partial analysis for N = 1. They are total variation diminishing under a suitable CFL condition and therefore converge to a weak solution. Numerical examples illustrate that these schemes, including a more accurate version based on MUSCL extrapolation, are competitive in accuracy and efficiency with several existing schemes. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1109–1136, 2016     

Incremental stochastic model for the temporal distribution of peak traffic demand

This paper presents the incremental stochastic user equilibrium (ISUE) model for predicting how travellers select their departure times from an origin in a single origin-destination pair system if they have desired times of arrival at the destination. The temporal distribution of the peak traffic is the result of commuters' selections of departure times. The stochastic user equilibrium (SUE) model is one of the techniques for estimating this distribution, but is computationally cumbersome to apply. In addition, existence and uniqueness of the solution to the SUE formulation and its approximations have not been proven. The ISUE is another approximation to the SUE, and it is based on an approach for which existence and uniqueness of the solution have been established.     

An extension of Krugman’s core–periphery model to the case of a continuous domain: Existence and uniqueness of solutions of a system of nonlinear integral equations in spatial economics

We consider a model that is an extension of Krugman’s core–periphery model to the case of a bounded closed domain included in a Euclidean space. We can describe the relation of the density of workers, the density of nominal wages, and the density of real wages by the system of nonlinear integral equations of the model. If we obtain a solution of the system under the condition that the density of workers is given, then the solution is called a short-run equilibrium. In this paper we prove that this model has a short-run equilibrium, and we obtain a sufficient condition for its uniqueness. Moreover we obtain upper and lower estimates for short-run equilibria, and we construct a useful iteration scheme to numerically obtain short-run equilibria.