与M-矩阵相关的一类二次矩阵方程的新迭代法（英文）

本文研究与M-矩阵相关的一类二次矩阵方程的数值解法.这类方程源于马尔可夫链的带噪Wiener-Hopf问题,其解中具有实际意义的是M-矩阵解.通过简单的变换,将该二次矩阵方程转化为M-矩阵代数Riccati方程.提出一种新的迭代方法,并对其进行收敛性分析.数值实验表明,新的迭代方法是可行的,且在一定条件下比现有的一些方法更为有效.     

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

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

摄动离散矩阵Lyapunov方程向后误差分析

研究摄动离散矩阵Lyapunov方程解的向后误差,利用矩阵Kronecker积的性质以及矩阵范数的性质,给出方程近似解的向后误差界,最后通过数值例子说明解的向后稳定性.     

矩阵方程X-A*XqA=I(0＜q＜1)Hermite正定解的扰动分析

首先证明了非线性矩阵方程X-A~*X~qA=I(0相似文献   

一类四元数矩阵三次方程的可解性(英文)

确立了某类分块矩阵[M（₁1） M₁₂ XM₂₁ Y M₂₃Z M₃₂ M₃₃]的最大秩公式,其中,X,Y和Z是三个受限于四元数线性矩阵方程A₁X=C₁,XB₁=C₂,A₂Y=D₁,YB₂=D₂,A₃Z=E₁,ZB₃=E₂的变量矩阵.作为该公式的一项应用,我们推导出上述矩阵方程解集等同于某类四元数三次矩阵方程组A₁X=C₁,XB₁=C₂,A₂Y=D₁,YB₂=D₂,A₃Z=E₁,ZB₃=E₂,XYZ=J解集的条件.     

二次矩阵方程X~2-bX-C=O的正定解

研究二次矩阵方程X2-bX-C=O(b>0,C为n×n阶正定阵)的正定解,证明了解的存在唯一性并且给出了求解方法.     

代数Riccati方程可稳解的条件数

１．引言 关于代数Ｒｉｃｃａｔｉ方程（ＡＲＥ）的研究是大量的．从数值角度看,有关数值方法,扰动理论的研究已比较深入．而关于条件数理论的研究则还不多［３］,［６］． Ｒｙｅｒｓ［１］研究了时连续代数 Ｒｉｃｃａｔｉ方程可稳解的条件数; Ｋｅｎｎｅｙ和 Ｈｅｗｅｒ［３］讨论了时连续代数 Ｒｉｃｃａｔｉ方程（以下简称 ＣＴＡＲＥ）可稳解的敏度分析,给出了一阶扰动界,引进了条件数; Ｓｕｎ［６］从最佳向后扰动理论角度研究了时离散代数Ｒｉｃｃａｔｉ方程（以下简称ＤＴＡＲＥ）可稳解的条件数;徐树方［８］针对 ＣＴＡ…     

M-矩阵Sylvester方程的一类Smith-like迭代法（英文）

本文研究了M-矩阵Sylvester方程的数值解法,这类矩阵方程广泛出现在科学计算和工程应用的许多领域.利用M-矩阵的性质和Smith方法的思想,提出了一类Smith-like迭代法以求解M-矩阵Sylvester方程,并给出了新方法的收敛性分析.数值实验表明,新方法是可行的,而且在一定条件下也是较为有效的.     

一类线性方程组的结构向后误差分析

考虑如下结构线性方程组(A B C 0)(x y)=(a b),其中A∈R~(m×m),B∈R~(m×n),C∈R~(n×m).本文给出该类结构方程组的结构向后扰动误差的显式表达式.数值例子表明求解该类问题稳定的算法得到的解不必是强稳定的.     

带多重右边的不定最小二乘问题的条件数

本文研究了带多重右边的不定最小二乘问题的条件数,给出了范数型、混合型及分量型条件数的表达式,同时,也给出了相应的结构条件数的表达式.所考虑的结构矩阵包含Toeplitz 矩阵、Hankel矩阵、对称矩阵、三对角矩阵等线性结构矩阵与Vandermonde矩阵、Cauchy矩阵等非线性结构矩阵.数值例子显示结构条件数总是紧于非结构条件数.     

Condition numbers and backward perturbation bound for linear matrix equations

This paper deals with the normwise perturbation theory for linear (Hermitian) matrix equations. The definition of condition number for the linear (Hermitian) matrix equations is presented. The lower and upper bounds for the condition number are derived. The estimation for the optimal backward perturbation bound for the Hermitian matrix equations is obtained. Copyright © 2010 John Wiley & Sons, Ltd.     

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).     

Backward error and perturbation bounds for high order Sylvester tensor equation

In this article, we investigate the backward error and perturbation bounds for the high order Sylvester tensor equation (STE). The bounds of the backward error and three types of upper bounds for the perturbed STE with or without dropping the second order terms are presented. The classic perturbation results for the Sylvester equation are extended to the high order case.     

On the condition number of linear least squares problems in a weighted Frobenius norm

LetA be anm × n, m n full rank real matrix andb a real vector of sizem. We give in this paper an explicit formula for the condition number of the linear least squares problem (LLSP) defined by min Ax–b₂,x ⁿ . Let and be two positive real numbers, we choose the weighted Frobenius norm [A, b] _F on the data and the usual Euclidean norm on the solution. A straightforward generalization of the backward error of [9] to this norm is also provided. This allows us to carry out a first order estimate of the forward error for the LLSP with this norm. This enables us to perform a complete backward error analysis in the chosen norms.Finally, some numerical results are presented in the last section on matrices from the collection of [5]. Three algorithms have been tested: the QR factorization, the Normal Equations (NE), the Semi-Normal Equations (SNE).     

Numerical analysis of a quadratic matrix equation

The quadratic matrix equation AX²+ BX + C = 0in n x nmatricesarises in applications and is of intrinsic interest as oneof the simplest nonlinear matrix equations. We give a completecharacterization of solutions in terms of the generalized Schurdecomposition and describe and compare various numerical solutiontechniques. In particular, we give a thorough treatment offunctional iteration methods based on Bernoulli’s method.Other methods considered include Newton’s method with exact line searches, symbolic solution and continued fractions.We show that functional iteration applied to the quadraticmatrix equation can provide an efficient way to solve the associated quadratic eigenvalue problem (²A + B + C)x = 0.     

Perturbation analysis of the matrix equation

Consider the nonlinear matrix equation X-A^*X^-pA=Q with 0<p1. This paper shows that there exists a unique positive definite solution to the equation. A perturbation bound and the backward error of an approximate solution to this solution is evaluated. We also obtain explicit expressions of the condition number for the unique positive definite solution. The theoretical results are illustrated by numerical examples.     

A modified Tikhonov regularization method for an axisymmetric backward heat equation

This paper deals with the inverse time problem for an axisymmetric heat equation. The problem is ill-posed. A modified Tikhonov regularization method is applied to formulate regularized solution which is stably convergent to the exact one. estimate between the approximate solution and exact technical inequality and improving a priori smoothness Meanwhile, a logarithmic-HSlder type error solution is obtained by introducing a rather assumption.     

Inaccurate linear equation system with a restricted-rank error matrix

It is well-known that the solution set of an interval linear equation system is a union of convex polyhedra the number of which increases, in general, exponentially with the problem size. As a consequence, the problem of finding the interval hull of the solution set is NP-hard as J. Rohn and V. Kreinovich proved in [13]. The purpose of this paper is to show that the solution set analysis can be simplified substantially provided the rank of the error matrix is restricted even if the assumption of interval character of data errors is replaced by a more general one. Especially, in the case of a rank-one error matrix we have to look into at most two convex subsets. Besides, a dual approach to describing the solution set is discussed. The original version of this approach was suggested in [7].