| Basemi I. Selim,杜磊,于波,朱宣儒.解一般矩阵方程的GPBiCG(m,l)法[J].数学研究及应用,2019,39(4):408~432 |
| 解一般矩阵方程的GPBiCG(m,l)法 |
| The GPBiCG($m,l$) Method for Solving General Matrix Equations |
| 投稿时间：2019-01-20 修订日期：2019-03-05 |
| DOI：10.3770/j.issn:2095-2651.2019.04.008 |
| 中文关键词: GPBiCG$(m,l)$法 Krylov子空间法 矩阵方程 Kronecker积 向量化算子 |
| 英文关键词:GPBiCG($m,l$) method Krylov Subspace method matrix equations Kronecker product vectorization operator |
| 基金项目:国家自然科学基金(Nos.11501079; 11571061), 埃及高等教育委员会. |
|
| 摘要点击次数: 1305 |
| 全文下载次数: 764 |
| 中文摘要: |
| 最近提出的用于解非对称线性方程组$Ax=b$的广义乘积型双共轭梯度法GPBiCG$(m,l)$是一种基于GPBiCG和BiCGSTAB法的混合型算法,该算法在许多数值试验中都有不错的收敛表现.借助Kronecker积和向量化算子,本文将GPBiCG$(m,l)$法做了推广并用之解一般矩阵方程Equation 1和一般离散时间周期矩阵方程组Equation 2, 其中包括出现在许多应用领域的Lyapunov, Stein 和Sylvester 等矩阵方程.通过数值试验与一些现有算法对比,检验了所提GPBiCG$(m,l)$法的准确性和有效性. |
| 英文摘要: |
| The generalized product bi-conjugate gradient (GPBiCG($m,l$)) method has been recently proposed as a hybrid variant of the GPBiCG and the BiCGSTAB methods to solve the linear system $Ax = b$ with non-symmetric coefficient matrix, and its attractive convergence behavior has been authenticated in many numerical experiments. By means of the Kronecker product and the vectorization operator, this paper aims to develop the GPBiCG($m,l$) method to solve the general matrix equation $$\sum^{p}_{i=1}{\sum^{s_{i}}_{j=1} A_{ij}X_{i}B_{ij}} = C,$$ and the general discrete-time periodic matrix equations $$\sum^{p}_{i=1}{\sum^{s_{i}}_{j=1} (A_{i,j,k}X_{i,k}B_{i,j,k}+C_{i,j, k}X_{i,k+1}D_{i,j,k})} = M_{k},~~k = 1, 2, \ldots,t,$$ which include the well-known Lyapunov, Stein, and Sylvester matrix equations that arise in a wide variety of applications in engineering, communications and scientific computations. The accuracy and efficiency of the extended GPBiCG($m,l$) method assessed against some existing iterative methods are illustrated by several numerical experiments. |
| 查看全文 查看/发表评论 下载PDF阅读器 |
|
|
|
