不变子空间上特征值的扰动

1引言设矩阵A∈C~(n×n),B∈C~(m×m),Q∈C~(n×m)为列满秩矩阵,令R=AQ-QB．当R的范数很小的时候,我们分析矩阵B的特征值对A的特征值的逼近性．当A,B都是Hermite阵时,上述问题已经被Kahan解决．近年来,对可对角化矩阵的情形,取得了一些新的成果．[4][5][6]中给出了几个范数不等式,并应用于矩阵特征值     

关于Kahan定理的一个注记

设A∈C~(n×n),B∈C~(k×k)均为Hermite矩阵,它们的特征值分别为{λ_j}_(j=1)~n和{μ_j}_(j=1)~k(k≤n);Q∈~(n×k)为列满秩矩阵.令 (1) 则存在A的k个特征值λ_(j_2),λ_(j_2),…,λ_(j_k),使得 (2) 其中σ_k为Q的最小奇异值,||·||_2表示矩阵的谱范数.这是著名的Kahan定理·1996年曹志浩等在[2]中将(2)加强为 (3) 这是Kahan的猜想.在本文中,我们讨论将Kahan定理中“B为k阶Hermite矩阵”改为B为k阶(任意)方阵后,特征值的扰动估计,有以下结果. 定理 设A∈C~(n×n)为Hermite矩阵,其特征值为{λ_j}_(j=1)~n,B∈C~(k×k)的特征值为{μ_j}_(j=1)~k,而Q∈C~(n×k)为列满秩矩阵.则存在A的k个特征值λ_(j_1),λ_(j_2),…,λ_(j_k),使得     

非线性规划中一个梯度投影法的收敛性质

我们考虑非线性规划问题(P)■f(x),其中R={x|Ax=a,Bx≤b},A是p×n矩阵,其秩为p,B是q×n矩阵,x∈E~n,a∈E~p,b∈E~q,f(x)∈C~1.我们以R~*表示(P)的最优解集合,并假定R非空.最近,M.S.Bazaraa与J.J.Goode     

求解Sylvester方程的广义非对称PMHSS算法

<正>1引言考虑如下Sylvester方程:AX+XB=F(1)这里A∈C~(m×m),B∈C~(n×n),F∈C~(m×n)是复数矩阵.令A=W+iT,B=U+iV,Q,T∈R~(m×m),U,V∈R~(n×n)都是实对称矩阵,且W,U是不定的,T,V是正定的.我们假定-TW≤T,-VU≤V.对于任意矩阵W和T,WT(W≤T)意味着T-W是     

广义特征值的扰动界限(Ⅱ)

本节将利用广义特征多项式的概念来研究广义特征值扰动界的上界估计.设A,B∈C~(n×n),首先定义一列算子:     

正则化HSS预处理鞍点矩阵的特征值估计

最近,Bai和Benzi针对鞍点问题提出了一类正则化HSS（Regularized Hermitian and skew-Hermitian splitting,RHSS）预处理子（BIT Numer.Math.,57（2017）287-311）.为了进一步分析RHSS预处理子的效果,本文重点研究了RHSS预处理鞍点矩阵特征值的估计,分析了复特征值实部和模的上下界、实特征值的上下界,还给出了特征值均为实数的充分条件.当正则化矩阵取为零矩阵时,RHSS预处理子退化为HSS预处理子,分析表明本文给出的复特征值实部的界比已有的结果更精确.数值算例验证了本文给出的理论结果.     

等式约束加权线性最小二乘问题的解法

1 引言 在实际应用中常会提出解等式约束加权线性最小二乘问题 min||b-Ax||_M,(1.1) x∈C~n s.t.Bx=d, 其中B∈C~(p×n),A∈C~(q×n),d∈C~p,b∈C~q,M∈C~(q×q)为Hermite正定阵. 对于问题(1.1),目前已有多种解法,见文[1—3).本文将利用广义逆矩阵的知识,给出(1.1)的通解及迭代解法.本文中关于矩阵广义逆与投影算子(矩阵)的记号基本上与文[4]的相同.例如,A~+表示A的MP逆,P_L表示到子空间L上的正交投影算子,λ_(max)(MAY)表示矩阵M~(1/2)AY的最大特征值.我们还要用到广义BD逆的概念: 设A∈C~(n×n),L为C~n的子空间,则称A_(L)~(+)=P_L(AP_L+P_L⊥)~+为A关于L的广义BD逆.     

半张量积下矩阵方程组AX=B,XC=D的最小二乘解

本文研究了半张量积下矩阵方程组AX=B,XC=D在不同情况下的最小二乘解X*∈R~(p×q),其中矩阵A∈R~(m×n),B∈R~(h×k),C∈R~(a×b),D∈R~(l×d)给定.根据半张量积的定义将其转变为普通乘积下的矩阵方程组,再结合矩阵奇异值分解及矩阵微分给出该方程组在不同情况下最小二乘解的解析表达式,并用数值算例加以验证.     

幂等矩阵下Yang-Baxter-like方程解的新的显式表达式

正1引言设C~(m×n)表示m×n阶复矩阵的集合,I_n表示n阶单位矩阵.对于矩阵A∈C~(m×n),A~*表示它的共轭转置矩阵.设矩阵A∈C~(n×n),如果A~2=A,则称矩阵A为幂等矩阵;如果A~2=A=A~*,则称矩阵A为正交投影矩阵.设A∈C~(n×n)本文主要研究下面的二次矩阵方程AXA=XAX,(1.1)称之为Yang-Baxter-like方程,因为其与统计物理中分别由Yang[1]和Baxter[2]独立得到的经典Yang-Baxter方程相似.     

界约束下算子方程最小二乘问题的条件梯度法

研究如下界约束下算子方程最小二乘问题:min x∈Ω‖L(X:A_1,…,At;B_1,…,B_t)-T‖~2,其中‖.‖为Frobenius范数,L(X:A_1…A_t;B_1,…,B_t)为关于X的线性矩阵算子(或齐次线性变换),Ai∈R~(p×m),B_j∈R~(n×q)i,j=1,…,n为算子L的系数矩阵,丁为右端矩阵,ΩR~(m×n)为界约束凸集合.提出了求解问题的条件梯度迭代算法及其简要收敛性分析,并给出条件梯度算法的几类加速形式.随机数据和图像恢复模型数据的实验结果表明说明算法是可行高效的.     

Numerical enclosure for each eigenvalue in generalized eigenvalue problem

An algorithm for enclosing all eigenvalues in generalized eigenvalue problem Ax=λBx is proposed. This algorithm is applicable even if A∈C^n×n is not Hermitian and/or B∈C^n×n is not Hermitian positive definite, and supplies nerror bounds while the algorithm previously developed by the author supplies a single error bound. It is proved that the error bounds obtained by the proposed algorithm are equal or smaller than that by the previous algorithm. Computational cost for the proposed algorithm is similar to that for the previous algorithm. Numerical results show the property of the proposed algorithm.     

Congruence of Hermitian matrices by Hermitian matrices

Two Hermitian matrices A,B∈M_n(C) are said to be Hermitian-congruent if there exists a nonsingular Hermitian matrix C∈M_n(C) such that B=CAC. In this paper, we give necessary and sufficient conditions for two nonsingular simultaneously unitarily diagonalizable Hermitian matrices A and B to be Hermitian-congruent. Moreover, when A and B are Hermitian-congruent, we describe the possible inertias of the Hermitian matrices C that carry the congruence. We also give necessary and sufficient conditions for any 2-by-2 nonsingular Hermitian matrices to be Hermitian-congruent. In both of the studied cases, we show that if A and B are real and Hermitian-congruent, then they are congruent by a real symmetric matrix. Finally we note that if A and B are 2-by-2 nonsingular real symmetric matrices having the same sign pattern, then there is always a real symmetric matrix C satisfying B=CAC. Moreover, if both matrices are positive, then C can be picked with arbitrary inertia.     

Bergman-Hartogs型域的全纯自同构群

我们考虑一类以有界对称域D为底的Bergman-Hartogs型域Ω={(wm(1),...,w(r),z)∈C1×···×Cmr×D:∥w(1)∥2p1+···+∥w(r)∥2prKD(z,z)-q},其中KD(z,z)是D上的Bergman核函数,r 1且为正整数,参数p1,...,pr1和q0为实数.我们给出它的全纯自同构群,并且证明当r=1时此自同构群为最大全纯自同构群;当r1时,若Ω的全纯自同构变换F将(0,z)∈{0}×D映到(0,z*)∈{0}×D,则F在我们给出的全纯自同构群中.     

Perturbation Analysis for the Eigenvalue Problem of a Formal Product of Matrices

We study the perturbation theory for the eigenvalue problem of a formal matrix product A ₁ ^{s 1} ··· A _p ^{s p}, where all A _k are square and s _k {–1, 1}. We generalize the classical perturbation results for matrices and matrix pencils to perturbation results for generalized deflating subspaces and eigenvalues of such formal matrix products. As an application we then extend the structured perturbation theory for the eigenvalue problem of Hamiltonian matrices to Hamiltonian/skew-Hamiltonian pencils.     

上三角算子矩阵的谱

设X,y是Banach空间,对A∈B(X),B∈B(y),C∈B(Y,X),以M_C记X⊕Y上的算子(ACOB).本文给出了算子M_C的20种谱的结构表示,18种谱的填洞性质以及关于这些问题的有趣例子.     

Nearest southeast submatrix that makes multiple a prescribed eigenvalue. Part 1

Given four complex matrices A,B,C and D, where A∈C^n×n and D∈C^m×m, and given a complex number z₀: What is the (spectral norm) distance from D to the set of matrices X∈C^m×m such that z₀ is a multiple eigenvalue of the matrix

     

Neville elimination for rank-structured matrices

In this paper it is shown that Neville elimination is suited to exploit the rank structure of an order-r quasiseparable matrix A∈C^n×n by providing a condensed decomposition of A as product of unit bidiagonal matrices, all together specified by O(nr) parameters, at the cost of O(nr³) flops. An application of this result for eigenvalue computation of totally positive rank-structured matrices is also presented.     

The determinant bundle on the moduli space of stable triples over a curve

We construct a holomorphic Hermitian line bundle over the moduli space of stable triples of the form (E₁, E₂,?), where E₁ and E₂ are holomorphic vector bundles over a fixed compact Riemann surfaceX, and?: E₂ → E₁ is a holomorphic vector bundle homomorphism. The curvature of the Chern connection of this holomorphic Hermitian line bundle is computed. The curvature is shown to coincide with a constant scalar multiple of the natural Kähler form on the moduli space. The construction is based on a result of Quillen on the determinant line bundle over the space of Dolbeault operators on a fixed C^∞ Hermitian vector bundle over a compact Riemann surface.     

About a construction and some analysis of time inhomogeneous diffusions on monotonely moving domains

We construct and analyze in a very general way time inhomogeneous (possibly also degenerate or reflected) diffusions in monotonely moving domains E⊂R×R^d, i.e. if E_t?{x∈R^d|(t,x)∈E}, t∈R, then either E_s⊂E_t, ∀s?t, or E_s⊃E_t, ∀s?t, s,t∈R. Our major tool is a further developed L²(E,m)-analysis with well chosen reference measure m. Among few examples of completely different kinds, such as e.g. singular diffusions with reflection on moving Lipschitz domains in R^d, non-conservative and exponential time scale diffusions, degenerate time inhomogeneous diffusions, we present an application to what we name skew Bessel process on γ. Here γ is either a monotonic function or a continuous Sobolev function. These diffusions form a natural generalization of the classical Bessel processes and skew Brownian motions, where the local time refers to the constant function γ≡0.     

“Localized” self-adjointness of Schrödinger type operators on Riemannian manifolds

We prove self-adjointness of the Schrödinger type operator , where ∇ is a Hermitian connection on a Hermitian vector bundle E over a complete Riemannian manifold M with positive smooth measure dμ which is fixed independently of the metric, and V∈L_loc¹(EndE) is a Hermitian bundle endomorphism. Self-adjointness of H_V is deduced from the self-adjointness of the corresponding “localized” operator. This is an extension of a result by Cycon. The proof uses the scheme of Cycon, but requires a refined integration by parts technique as well as the use of a family of cut-off functions which are constructed by a non-trivial smoothing procedure due to Karcher.