A FAST SINE TRANSFORM ALGORITHM FOR TOEPLITZ MATRICES AND ITS APPLICATIONS

In this paper, a fast algorithm for the discrete sine transform(DST) of a Toeplitz matrix of order N is derived. Only O(N log N) O(M) time is needed for the computation of M elements. The auxiliary storage requirement is O(N). An application of the new fast algorithm is also discussed.     

FAST PARALLELIZABLE METHODS FOR COMPUTING INVARIANT SUBSPACES OF HERMITIAN MATRICES

We propose a quadratically convergent algorithm for computing the invariant subspaces of an Hermitian matrix. Each iteration of the algorithm consists of one matrix-matrix multiplication and one QR decomposition. We present an accurate convergence analysis of the algorithm without using the big O notation. We also propose a general framework based on implicit rational transformations which allows us to make connections with several existing algorithms and to derive classes of extensions to our basic algorithm with faster convergence rates. Several numerical examples are given which compare some aspects of the existing algorithms and the new algorithms.     

ITERATIVE ILU PRECONDITIONERS FOR LINEAR SYSTEMS AND EIGENPROBLEMS

Iterative ILU factorizations are constructed,analyzed and applied as preconditioners to solve both linear systems and eigenproblems.The computational kernels of these novel Iterative ILU factorizations are sparse matrix-matrix multiplications,which are easy and efficient to implement on both serial and parallel computer architectures and can take full advantage of existing matrix-matrix multiplication codes.We also introduce level-based and threshold-based algorithms in order to enhance the accuracy of the proposed Iterative ILU factorizations.The results of several numerical experiments illustrate the efficiency of the proposed preconditioners to solve both linear systems and eigenvalue problems.     

BLOCK-TRIANGULAR PRECONDITIONERS FOR SYSTEMS ARISING FROM EDGE-PRESERVING IMAGE RESTORATION

Signal and image restoration problems are often solved by minimizing a cost function consisting of an l2 data-fidelity term and a regularization term. We consider a class of convex and edge-preserving regularization functions. In specific, half-quadratic regularization as a fixed-point iteration method is usually employed to solve this problem. The main aim of this paper is to solve the above-described signal and image restoration problems with the half-quadratic regularization technique by making use of the Newton method. At each iteration of the Newton method, the Newton equation is a structured system of linear equations of a symmetric positive definite coefficient matrix, and may be efficiently solved by the preconditioned conjugate gradient method accelerated with the modified block SSOR preconditioner. Our experimental results show that the modified block-SSOR preconditioned conjugate gradient method is feasible and effective for further improving the numerical performance of the half-quadratic regularization approach.     

分块K—循环Toeplitz矩阵求逆的快速付氏变换法

1算法描述及推导 Toeplitz矩阵及Toeplitz系统的求解在谱分析、线性预测、误差控制码、自回归滤波器设计等领域内起着重要的作用~[1-3]，而分块Toeplitz矩阵在计算机的时序分析、自回归时序模型滤波中也经常出现~[4]。对一般Toeplitz矩阵求逆，其算术复杂性为O(n~2)~[5]-[6]，其中n为Toepleitz矩阵的阶，而K-循环Toeplitz矩阵的求逆，其算术复杂性可降为O(nlog_2n)，本文提供了mn附分块K-循环Toeplitz矩阵求逆的一种快速付氏变换算法，其算术复杂性为O（mnlog_2mn）.     

AN INVERSE EIGENVALUE PROBLEM FOR JACOBI MATRICES

Let T1,n be an n x n unreduced symmetric tridiagonal matrix with eigenvaluesand is an (n - 1) x (n - 1) submatrix by deleting the kth row and kth column, k = 1, 2,be the eigenvalues of T1,k andbe the eigenvalues of Tk+1,nA new inverse eigenvalues problem has put forward as follows: How do we construct anunreduced symmetric tridiagonal matrix T1,n, if we only know the spectral data: theeigenvalues of T1,n, the eigenvalues of Ti,k-1 and the eigenvalues of Tk+1,n?Namely if we only know the data: A1, A2, An,how do we find the matrix T1,n? A necessary and sufficient condition and an algorithm ofsolving such problem, are given in this paper.     

AN INVERSE EIGENVALUE PROBLEM FOR JACOBI MATRICES

In this paper, we discuss an inverse eigenvalue problem for constructing a 2n × 2n Jacobi matrix T such that its 2n eigenvalues are given distinct real values and its leading principal submatrix of order n is a given Jacobi matrix. A new sufficient and necessary condition for the solvability of the above problem is given in this paper. Furthermore, we present a new algorithm and give some numerical results.     

对称正定Toeplitz方程组的多级迭代求解

二级迭代法亦称内外迭代法. 多级迭代法由多个二级迭代嵌套而成.这些方法特别适合于并行计算,同时可以理解为古典迭代法的延伸或共轭梯度法的预处理子.本文讨论了对称正定Toeplitz线性方程组多级迭代法. 首先,基于Toeplitz矩阵的结构, 我们给出了多级块Jacobi分裂,然后证明了每一级分裂均为P-正则分裂, 并证明了当每一级内迭代次数均为偶数时,迭代法的收敛性. 最后通过数值实例验证了此方法的有效性.     

TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH

We present a sufficient and necessary condition for a so-called Cnk pattern to have positive semidefinite (PSD) completion. Since the graph of the Cnk pattern is composed by some simple cycles, our results extend those given in [1] for a simple cycle. We also derive some results for a partial Toeplitz PSD matrix specifying the Cnk pattern to have PSD completion and Toeplitz PSD completion.     

AN IMPROVEMENT ON THE QL ALGORITHM FOR SYMMETRIC TRIDIAGONAL MATRICES*

This paper establishes an improvement on the QL algorithm for a symmetric tridiagonal matrix T so that we can work out the eigenvalues of T faster. Meanwhile, the new algorithm don't worsen the stability and precision of the former algorithm.     

GENUINE-OPTIMAL CIRCULANT PRECONDITIONERS FOR WIENER-HOPF EQUATIONS

1. IntroductionWienerHopf equations are integral equations defined on the haif line:where rr > 0, a(.) C L1(ro and g(.) E L2(at). Here R = (--oo,oo) and ty [0,oo). Inou-r discussions, we assume that a(.) is colljugate symmetric, i.e. a(--t) = a(t). WienerHop f equations arise in a variety of practical aPplicatiolls in mathematics and ellgineering, forinstance, in the linear prediction problems fOr stationary stochastic processes [8, pp.145--146],diffuSion problems and scattering problems […     

一类Toeplitz线性代数方程组的预处理GMRES方法

本文研究一类来源于分数阶特征值问题的Toeplitz线性代数方程组的求解.构造Strang循环矩阵作为预处理矩阵来求解该Toeplitz线性代数方程组,分析了预处理后系数矩阵的特征值性质.提出求解该线性代数方程组的预处理广义极小残量法（PGMRES）,并给出该算法的计算量.数值算例表明了该方法的有效性.     

基于均值的Toeplitz矩阵填充的子空间算法

本文提出一种基于均值的Toeplitz矩阵填充的子空间算法.通过在左奇异向量空间中对已知元素的最小二乘逼近,形成了新的可行矩阵;并利用对角线上的均值化使得迭代后的矩阵保持Toeplitz结构,从而减少了奇异向量空间的分解时间.理论上,证明了在一定条件下该算法收敛于一个低秩的Toeplitz矩阵.通过不同已知率的矩阵填充数值实验展示了Toeplitz矩阵填充的新算法比阈值增广Lagrange乘子算法在时间上和精度上更有效.     

THE SINE TRANSFORM OPERATOR IN THE BANACH SPACE OF SYMMETRIC MATRICES AND ITS APPLICATION IN IMAGE RESTORATION

In this paper, we study an operator s which maps every n-by-n symmetric matrix A, to a matrix s(A_n) that minimizes || B_n-A_n || F over the set of all matrices B_n, that can be diagonalized by the sine transform. The matrix s(A_n), called the optimal sine transform preconditioner, is defined for any n-by-n symmetric matrices A_n. The cost of constructing s(A_n) is the same as that of optimal circulant preconditioner c(A_n) which is defined in [8], The s(A_n) has been proved in [6] to be a good preconditioner in solving symmetric Toeplitz systems with the preconditioned conjugate gradient (PCG) method. In this paper, we discuss the algebraic and geometric properties of the operator s, and compute its operator norms in Banach spaces of symmetric matrices. Some numerical tests and an application in image restoration are also given.     

几乎各向同性的高维空间分数阶扩散方程的分块快速正则Hermite分裂预处理方法

一类空间分数阶扩散方程经过有限差分离散后所得到的离散线性方程组的系数矩阵是两个对角矩阵与Toeplitz型矩阵的乘积之和.在本文中,对于几乎各向同性的二维或三维空间分数阶扩散方程的离散线性方程组,采用预处理Krylov子空间迭代方法,我们利用其系数矩阵的特殊结构和具体性质构造了一类分块快速正则Hermite分裂预处理子.通过理论分析,我们证明了所对应的预处理矩阵的特征值大部分都聚集于1的附近.数值实验也表明,这类分块快速正则Hermite分裂预处理子可以明显地加快广义极小残量(GMRES)方法和稳定化的双共轭梯度(BiCGSTAB)方法等Krylov子空间迭代方法的收敛速度.     

SINE TRANSFORM MATRIX FOR SOLVING TOEPLITZ MATRIX PROBLEMS

1. IntroductionStrang[1] first studied the use of circulallt matrices C for solving systems of linear eqllationsTi x = b witha symmetric positive definite Toeplitz matrix.Numerous authors such as T.Chan[2],R.Chan,etc.[3],[4],[5], Tyrtyshnikov[6], Huckle[7] and T.Ku and C.Kuo[8] proposed differentfamilies of circulallt / skew- circulant precondit ioners.Appling the preconditioned conjugate gradient algorithm(PCGA) to solve the systems Ti x -b, we must find a preconditioner P such that P…     

INERTIA SETS OF SYMMETRIC SIGN PATTERN MATRICES

1　IntroductionIn qualitative and combinatorial matrix theory,we study properties ofa matrix basedon combinatorial information,such as the signs of entries in the matrix.A matrix whoseentries are from the set{ + ,-,0 } is called a sign pattern matrix ( or sign pattern,or pat-tern) .We denote the setof all n× n sign pattern matrices by Qn.For a real matrix B,sgn( B) is the sign pattern matrix obtained by replacing each positive( respectively,negative,zero) entry of B by+ ( respectively,-,0 )…     

求解带Toeplitz矩阵的线性互补问题的一类预处理模系矩阵分裂迭代法

针对系数矩阵为对称正定Toeplitz矩阵的线性互补问题,本文提出了一类预处理模系矩阵分裂迭代方法.先通过变量替换将线性互补问题转化为一类非线性方程组,然后选取Strang或T.Chan循环矩阵作为预优矩阵,利用共轭梯度法进行求解.我们分析了该方法的收敛性.数值实验表明,该方法是高效可行的.     

长方矩阵加权群逆的迭代法

1引 言与引理 最近,文[1]定义了长方矩阵的一种加权群逆:设A∈Cm×n,W∈Cn×m.称满足下列矩阵方程组的矩阵X∈Cm×n为A的加W权群逆:(W1)AWXWA=A, (W2)XWAWX=X, (W3)AWX=XWA通常记A的加W权群逆为A#W.若A#W存在,则它是唯一的.