一类二次矩阵方程的条件数和后向误差

主要讨论一类二次矩阵方程X^2-EX-F=0的条件数和后向误差,其中E是一个对角矩阵,F是一个M矩阵.这类二次矩阵方程来源于Markov链的噪声Wiener-Hopf问题.实际问题中人们感兴趣的是它的M矩阵的解.应用Rice创立的基于Frobenius范数下的条件数理论,导出此类二次矩阵方程的M矩阵解的条件数的显式表达式.同时,也给出近似解的后向误差的定义以及一个可计算的表达式.最后,通过数值例子验证理论结果是有效的.     

Condition numbers and perturbation analysis for the Tikhonov regularization of discrete ill‐posed problems

One of the most successful methods for solving the least‐squares problem min_x∥Ax?b∥₂ with a highly ill‐conditioned or rank deficient coefficient matrix A is the method of Tikhonov regularization. In this paper, we derive the normwise, mixed and componentwise condition numbers and componentwise perturbation bounds for the Tikhonov regularization. Our results are sharper than the known results. Some numerical examples are given to illustrate our results. Copyright © 2010 John Wiley & Sons, Ltd.     

New perturbation bounds and condition numbers for the hyperbolic QR factorization

Using the modified matrix-vector equation approach, the technique of Lyapunov majorant function and the Banach fixed point theorem, we obtain some new rigorous perturbation bounds for R factor of the hyperbolic QR factorization under normwise perturbation. These bounds are always tighter than the one given in the literature. Moreover, the optimal first-order perturbation bounds and the normwise condition numbers for the hyperbolic QR factorization are also presented.     

矩阵方程ATXA=D的条件数与向后扰动分析

讨论矩阵方程ATXA=D,该方程源于振动反问题和结构模型修正.本文利用Moore-Penrose广义逆的性质,给出该方程解的条件数的上、下界估计.同时,利用Schauder不动点理论给出该方程的向后扰动界,这些结果可用于该矩阵方程的数值计算.     

Mixed,componentwise condition numbers and small sample statistical condition estimation of Sylvester equations

We present a componentwise perturbation analysis for the continuous‐time Sylvester equations. Componentwise, mixed condition numbers and new perturbation bounds are derived for the matrix equations. The small sample statistical method can also be applied for the condition estimation. These condition numbers and perturbation bounds are tested on numerical examples and compared with the normwise condition number. The numerical examples illustrate that the mixed condition number gives sharper bounds than the normwise one. Copyright © 2011 John Wiley & Sons, Ltd.     

Perturbation theory for generalized and constrained linear least squares

The perturbation analysis of weighted and constrained rank‐deficient linear least squares is difficult without the use of the augmented system of equations. In this paper a general form of the augmented system is used to get simple perturbation identities and perturbation bounds for the general linear least squares problem both for the full‐rank and rank‐deficient problem. Perturbation identities for the rank‐deficient weighted and constrained case are found as a special case. Interesting perturbation bounds and condition numbers are derived that may be useful when considering the stability of a solution of the rank‐deficient general least squares problem. Copyright © 2000 John Wiley & Sons, Ltd.     

Systems of matrix linear differential equations of first order

Sufficient conditions for the solvability in quadratures of systems of matrix linear ordinary differential equations of first order with one-sided multiplication by variable matrix coefficients are given in this paper. These conditions are stated in terms of the theory of Lie algebras. We consider matrix equations of higher orders that are equivalent to such systems. An illustrative example is considered. Translated fromMatematicheskie Zametki, Vol. 66, No. 1, pp. 63–75, July, 1999.     

An exponential matrix method for solving systems of linear differential equations

This paper presents an exponential matrix method for the solutions of systems of high‐order linear differential equations with variable coefficients. The problem is considered with the mixed conditions. On the basis of the method, the matrix forms of exponential functions and their derivatives are constructed, and then by substituting the collocation points into the matrix forms, the fundamental matrix equation is formed. This matrix equation corresponds to a system of linear algebraic equations. By solving this system, the unknown coefficients are determined and thus the approximate solutions are obtained. Also, an error estimation based on the residual functions is presented for the method. The approximate solutions are improved by using this error estimation. To demonstrate the efficiency of the method, some numerical examples are given and the comparisons are made with the results of other methods. Copyright © 2012 John Wiley & Sons, Ltd.     

New upper bound formulas with parameters for Ramsey numbers

In this paper, we obtain some new results R(5,12)?848, R(5,14)?1461, etc., and we obtain new upper bound formulas for Ramsey numbers with parameters.     

Delayed perturbation of Mittag‐Leffler functions and their applications to fractional linear delay differential equations

In this paper, we propose a delayed perturbation of Mittag‐Leffler type matrix function, which is an extension of the classical Mittag‐Leffler type matrix function and delayed Mittag‐Leffler type matrix function. With the help of the delayed perturbation of Mittag‐Leffler type matrix function, we give an explicit formula of solutions to linear nonhomogeneous fractional delay differential equations.     

The homogeneous projective transformation of general quadratic matrix equations

^** Email: vassilios.tsachouridis{at}ieee.org^*** Email: basil.kouvaritakis{at}eng.ox.ac.uk Algebraic quadratic equations are a special case of a singlegeneralized algebraic quadratic matrix equation (GQME). Hence,the importance of that equation in science and engineering isevident. This paper focus on the study of solutions of thatGQME and a unified framework for the characterization and identificationof solutions at infinity and of finite solutions of generalquadratic algebraic matrix equations is presented. The analysisis based on the concept of homogeneous projective transformationfor general polynomial systems (Morgan, 1986). In addition,a numerical error analysis for the computed solutions is providedfor the assessment of numerical accuracy, stability and conditioningof the computed solutions. The proposed framework is independentof any numerical method and therefore it can be used along withvarious possible numerical methods for the GQME solution, especiallymatrix flow-based algorithms (Chu, 1994) (e.g. continuation/homotopy,Morgan, 1989).     

Small perturbation solutions for nonlocal elliptic equations

We present a small perturbation result for nonlocal elliptic equations, which says that for a class of nonlocal operators, the solutions are in C^σ+α for any α∈(0,1) as long as the solutions are small. This is a nonlocal generalization of a celebrated result of Savin in the case of second order equations.     

On Frobenius normwise condition numbers for Moore–Penrose inverse and linear least‐squares problems

Condition numbers play an important role in numerical analysis. Classical condition numbers are normwise: they measure the size of both input perturbations and output errors using norms. In this paper, we give explicit, computable expressions depending on the data, for the normwise condition numbers for the computation of the Moore–Penrose inverse as well as for the solutions of linear least‐squares problems with full‐column rank. Copyright © 2007 John Wiley & Sons, Ltd.     

伴有边界摄动的Volterra型积分微分方程组的奇摄动

讨论了伴有边界摄动的二阶非线性Volterra型积分微分方程组的奇摄动.在适当的条件下,利用对角化技巧证明了解的存在性,构造出解的渐近展开式并给出余项的一致有效估计.     

Structured Backward Error and Condition Number for Linear Systems of the Type A* Ax = b

We present a formulation for the structured condition number and for the structured backward error for the linear system A* Ax = b, when the rectangular matrix A is subjected to normwise perturbations. Perturbations on the data A and the solution x are measured in the Frobenius norm. Numerical experiments are provided that show the relevance of this condition number in the prediction of the computing error when solving such systems.     

Comparison between variational iteration method and homotopy perturbation method for linear and nonlinear partial differential equations with the nonhomogeneous initial conditions

In this study, linear and nonlinear partial differential equations with the nonhomogeneous initial conditions are considered. We used Variational iteration method (VIM) and Homotopy perturbation method (HPM) for solving these equations. Both methods are used to obtain analytic solutions for different types of differential equations. Four examples are presented to show the application of the present techniques. In these schemes, the solution takes the form of a convergent series with easily computable components. The present methods perform extremely well in terms of efficiency and simplicity. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

一类KKT系统的结构敏度分析

本文讨论一类KKT系统的敏度分析,这类KKT系统产生于用有限元方法离散Stokes方程,有结构特性．首先给出了最佳向后扰动界,接下来定义了偏条件数并导出了表达式．最后给出了新的扰动界．     

THE STRUCTURE SENSITIVITY ANALYSIS FOR TRIANGULAR SYSTEMS

Triangular systems play a fundamental role in matrix computations. It has become commonplace that triangular systems are solved to be more accurate even if they are ill-conditioned. In this paper, we define structured condition number and give structured (forward) perturbation bound. In addition, we derive the representation of optimal structured backward perturbation bound.     

正倒向随机微分方程组的数值解法

1990年,Pardoux和Peng(彭实戈)解决了非线性倒向随机微分方程(backward stochastic differential equation,BSDE)解的存在唯一性问题,从而建立了正倒向随机微分方程组(forward backward stochastic differential equations,FBSDEs)的理论基础;之后,正倒向随机微分方程组得到了广泛研究,并被应用于众多研究领域中,如随机最优控制、偏微分方程、金融数学、风险度量、非线性期望等.近年来,正倒向随机微分方程组的数值求解研究获得了越来越多的关注,本文旨在基于正倒向随机微分方程组的特性,介绍正倒向随机微分方程组的主要数值求解方法.我们将重点介绍讨论求解FBSDEs的积分离散法和微分近似法,包括一步法和多步法,以及相应的数值分析和理论分析结果.微分近似法能构造出求解全耦合FBSDEs的高效高精度并行数值方法,并且该方法采用最简单的Euler方法求解正向随机微分方程,极大地简化了问题求解的复杂度.文章最后,我们尝试提出关于FBSDEs数值求解研究面临的一些亟待解决和具有挑战性的问题.     

A class of backward doubly stochastic differential equations with discontinuous coefficients

In this work the existence of solutions of one-dimensional backward doubly stochastic differential equations （BDSDEs） with coefficients left-Lipschitz in y （may be discontinuous） and Lipschitz in z is studied. Also, the associated comparison theorem is obtained.