块循环矩阵方程组的新算法

1　基本概念形如 A=a1 a2 … a Na N a1 … a N- 1?彙?廰2 a3 … a1的矩阵称为由 a1 ,a2 ,… ,a N 生成的循环矩阵 .力学和工程中的轴对称结构的计算产生上述循环矩阵 [2 - 3] .以循环矩阵A为系数矩阵的方程组 ,称为循环矩阵方程组 .已有的求解循环矩阵方程组的办法主要是各种迭代法 ,如递推法及 SOR,SSOR,SAOR超松弛迭代法[2 - 6] 等 .定义 1　形如A =A1 A2 … ANAN A1… AN- 1?彙?廇2 A3… A1　 (Ai,i =1 ,2 ,… ,N为 m阶矩阵 )的矩阵称为由 A1 ,A2 ,… ,AN 生成的块循环矩阵 .定义 2　系数矩阵 A为块循环矩阵的方程组AX …     

方程组x′=Ax的新解法与计算机实现

在介绍B.VAN ROOT SELAAR的解方程组x′=Ax的一种新方法的基础上,对矩阵F(0)求法作了补充,对照以往通常的解法,分析了它的优越性.文章用完全开放性的Maple语言程序在计算机上实现了这种方法的应用,并通过生动的例子说明了同样是借助计算机强大的计算功能,新的解法在速度上要提高上百倍,更有实用价值.     

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

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

方程组x''''=Ax的新解法与计算机实现

在介绍B.VAN ROOT SELAAR的解方程组x'=Ax的一种新方法的基础上,对矩阵F(0)求法作了补充,对照以往通常的解法,分析了它的优越性.文章用完全开放性的Maple语言程序在计算机上实现了这种方法的应用,并通过生动的例子说明了同样是借助计算机强大的计算功能,新的解法在速度上要提高上百倍,更有实用价值.     

欧氏环上一次不定方程组的矩阵解法

利用欧几里德算法从理论上对多元一次不定方程组在欧氏环上的可逆线性变换下的解进行深入的研究 ,并提出用矩阵的初等变换求解欧氏环上多元一次不定方程组的算法 ,即矩阵解法 .     

接触问题中的快速凝缩消元法

用分离解法求解弹性接触问题时，在增量加载和迭代过程中，由于接触区某些节点的状态发生改变而导致方程组的系数矩阵某些行和列元素随之变化。根据此特点，本推导了一种新的自适应迭代算法-快速凝缩消元法，并给出具体的迭代步骤，避免了系数矩阵变化时必须重新形成矩阵的重复计算。     

相容线性方程组的Huang方法及其推广

本文从计算广义逆矩阵的Greville算法出发,给出并推广了Huang所提出的算法。这种新的算法可以用来求亚定、超定方程组的极小最小二乘解。计算实践表明,此法对解坏条件方程组具有良好的效果。     

离散椭圆型方程的矩阵方程表示及其直接解法

在[1]中讨论了解离散椭圆型方程的直接法,[2]、[3]中在椭圆型方程的类型和区域上得到了进一步的推广.它们的优点是计算起来快速且仅需要极小的存贮量.本文把求解的线性方程组写成矩阵方程形式,进一步发展矩阵分解法.首先,为了简单,如[1]那样考虑方程组     

非对称稀疏线性方程组的快速外存解法及其在无网格法计算中的应用

针对局部Petrov-Galerkin无网格法(MLPG)等无网格方法的计算所产生的大型非对称稀疏线性方程组,介绍了一种新的直接解法．与一般非对称求解过程不同,该解法从现有的对称正定解法中演变出来,其分解过程在矩阵的上、下三角阵中对称进行．新的矩阵分解算法可以通过修改对称矩阵分解算法的代码来实现,这提供了从对称解法到非对称解法的快捷转换．还针对MLGP法以及有限元法所产生的方程组开发了多块外存算法(multi-blocked out-of-core strategy)来扩大求解规模．测试结果证明该方法大幅度提高了大型非对称稀疏线性方程组的求解速度．     

一类线性隐式A稳定的单步法

1.引言 对于Stiff方程组初值问题的数值解法,Dahlquist在[1]中引进了 A稳定的概念,并且证明了显式的线性多步法(包括显式的Runge-Kutta方法)不可能是A稳定的.现在已经有许许多多隐式A稳定或Stiff稳定的方法,但绝大多数在数值解的过程中必须解由于隐式方法所产生的非线性方程组,而非线性方程组的求解过程往往又要采用Newton-Raphson迭代方法,因此需要计算方程y’=f(x,y)的右函数f(x,y)的Jacobi矩阵以及与此有关的逆矩阵.本文的主要思想是:既然在数值解过程中要计算f(x,y)的Jacobi矩阵,那么不妨在数值公式中明显的出现f(x,y)的一阶偏导数.我们将A稳定公式     

Circulant Wavelet Preconditioners for Solving Elliptic Differential Equations and Boundary Integral Equations

In this paper is discussed solving an elliptic equation and a boundary integral equation of the second kind by representation of compactly supported wavelets. By using wavelet bases and the Galerkin method for these equations, we obtain a stiff sparse matrix that can be ill-conditioned. Therefore, we have to introduce an operator which maps every sparse matrix to a circulant sparse matrix. This class of circulant matrices is a class of preconditioners in a Banach space. Based on having some properties in the spectral theory for this class of matrices, we conclude that the circulant matrices are a good class of preconditioners for solving these equations. We called them circulant wavelet preconditioners (CWP). Therefore, a class of algorithms is introduced for rapid numerical application.     

Parallel algorithms for solving tridiagonal and near-circulant systems

Many problems in mathematics and applied science lead to the solution of linear systems having circulant coefficient matrices. This paper presents a new stable method for the exact solution of non-symmetric tridiagonal circulant linear systems of equations. The method presented in this paper is quite competitive with Gaussian elimination both in terms of arithmetic operations and storage requirements. It is also competitive with the modified double sweep method. This method can be applied to solve the near-circulant tridiagonal system. In addition, the method is modified to allow for parallel processing.     

A note on inversion of Toeplitz matrices

It is shown that the invertibility of a Toeplitz matrix can be determined through the solvability of two standard equations. The inverse matrix can be denoted as a sum of products of circulant matrices and upper triangular Toeplitz matrices. The stability of the inversion formula for a Toeplitz matrix is also considered.     

在两循环矩阵类中求解线性方程组反问题的快速算法

本文利用多项式最大公因式 ,给出了线性方程组的反问题在 r-循环矩阵类和对称 r-循环矩阵类中有唯一解的充要条件 ,进而得到线性方程组在 r循环矩阵类和对称 r-循环矩阵类中的反问题求唯一解的算法 .最后给出了应用该算法的数值例子 .     

The structured distance to normality of Toeplitz matrices with application to preconditioning

A formula for the distance of a Toeplitz matrix to the subspace of {e^i?}‐circulant matrices is presented, and applications of {e^i?}‐circulant matrices to preconditioning of linear systems of equations with a Toeplitz matrix are discussed. Copyright © 2010 John Wiley & Sons, Ltd.     

A Note on Parallel Preconditioning for the All-at-Once Solution of Riesz Fractional Diffusion Equations

The $p$-step backward difference formula (BDF) for solving systems of ODEs can be formulated as all-at-once linear systems that are solved by parallel-in-time preconditioned Krylov subspace solvers (see McDonald et al. [36] and Lin and Ng [32]). However, when the BDF$p$ (2 ≤ $p$ ≤ 6) method is used to solve time-dependent PDEs, the generalization of these studies is not straightforward as $p$-step BDF is not selfstarting for $p$ ≥ 2. In this note, we focus on the 2-step BDF which is often superior to the trapezoidal rule for solving the Riesz fractional diffusion equations, and show that it results into an all-at-once discretized system that is a low-rank perturbation of a block triangular Toeplitz system. We first give an estimation of the condition number of the all-at-once systems and then, capitalizing on previous work, we propose two block circulant (BC) preconditioners. Both the invertibility of these two BC preconditioners and the eigenvalue distributions of preconditioned matrices are discussed in details. An efficient implementation of these BC preconditioners is also presented, including the fast computation of dense structured Jacobi matrices. Finally, numerical experiments involving both the one- and two-dimensional Riesz fractional diffusion equations are reported to support our theoretical findings.     

包含切比雪夫多项式的循环矩阵行列式的计算

行首加r尾r右循环矩阵和行尾加r首r左循环矩阵是两种特殊类型的矩阵,这篇论文中就是利用多项式因式分解的逆变换这一重要的技巧以及这类循环矩阵漂亮的结构和切比雪夫多项式的特殊的结构,分别讨论了第一类、第二类切比雪夫多项式的关于行首加r尾r右循环矩阵和行尾加r首r左循环矩阵的行列式,从而给出了行首加r尾r右循环矩阵和行尾加r首r左循环矩阵的行列式显式表达式.这些显式表达式与切比雪夫多项式以及参数r有关.这一问题的应用背景主要在循环编码,图像处理等信息理论方面.     

The computation of the square roots of circulant matrices

In the first part of this paper, we investigate the reduced forms of circulant matrices and quasi-skew circulant matrices. By using their properties we present two efficient algorithms to compute the square roots of circulant matrices and quasi-skew circulant matrices, respectively. Those methods are faster than the traditional algorithm which is based on the Schur decomposition. In the second part, we further consider circulant H-matrices with positive diagonal entries and develop two algorithms for computing their principal square roots. Those two algorithms have the common advantage that is they only need matrix-matrix multiplications in their iterative sequences, an operation which can be done very efficiently on modern high performance computers.     

Double circulant matrices

Double circulant matrices are introduced and studied. By a matrix-theoretic method, the rank r of a double circulant matrix is computed, and it is shown that any consecutive r rows of the double circulant matrix are linearly independent. As a generalization, multiple circulant matrices are also introduced. Two questions on square double circulant matrices are posed.     

具有位数3和4的双随机循环矩阵中的素元

本文研究了双随机循环矩阵中素元的分类问题．由于任一n阶双随机循环矩阵都可以唯一地表示为移位的n-1次一元多项式,从而可把双随机循环矩阵中素元的分类问题简化为解双随机循环矩阵上的一个方程．应用此原理,本文完全解决了判别具有位数3的n阶双随机循环矩阵是否为素元的问题,并给出了n阶双随机循环矩阵中一类具有位数4的素元．