广义逆在矩阵理论分析中有着重要的作用。文中讨论了中心对称矩阵A的广义中心对称{1, 4}逆的一种迭代算法，首先将广义{1, 4}逆转化为单变量线性矩阵方程组，然后建立求线性矩阵方阵组中心对称{1, 4}逆的修正共轭梯度算法(MCG算法)，证明了MCG算法的收敛性。数值算例表明，该算法具有很高的计算效率。

广义{1, 4}逆，中心对称解，修正共轭梯度算法，Moore-Penrose

Shijun Chen

Department of Basic Teaching and Research, Yango University, Fuzhou Fujian

Received: Mar. 19^th, 2021; accepted: Apr. 6^th, 2021; published: Apr. 22^nd, 2021

Generalized inverses play an important role in the analysis of matrix theory. In this paper, an iterative algorithm for generalized centrosymmetric {1, 4} inverse of centrosymmetric matrix is discussed. Firstly, the generalized {1, 4} inverse is transformed into a system of univariate linear matrix equations. Then, a modified conjugate gradient algorithm (MCG algorithm) is established for solving centrosymmetric {1, 4} inverses of linear matrices. The convergence of MCG algorithm is proved. Numerical examples show that the algorithm has high computational efficiency.

Keywords:Generalized {1, 4} Inverse, Centrosymmetric Solution, Modified Conjugate Gradient Algorithm, Moore-Penrose

Copyright © 2021 by author(s) and beplay安卓登录

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

E. H. Moore在1920年首次提出了广义逆矩阵的概念，R. Penrose明确给出了Moore广义逆矩阵定义，广义逆矩阵的研究进入了全新的时期。伴随计算科学的发展，广义逆矩阵在数理统计、系统理论、优化计算和控制论领域应用逐渐为人们所认识，因而大大推动了对广义逆的理论与应用的研究，使得广义逆理论得到迅速发展 [ 1 ] - [ 9 ]，刘红伟在文 [ 1 ] 中提出一种新算法求广义逆，该算法在某种条件下运算速度优于奇异值分解求广义逆算法；郭志荣在文 [ 2 ] 中从纯代数角度研究了代数广义逆的可加性与广义逆表示问题；梁少辉在文 [ 3 ] 中讨论了Quantale矩阵M-P广义逆的若干性质，给出了Quantale矩阵M-P广义逆的具体形式；柯圆圆在文 [ 4 ] 中给出环上(b, c)可逆元的一些等价刻画，并考虑环上2 × 2矩阵的(B, C)逆的存在性和表达式；欧阳光在文 [ 5 ] 中讨论了广义逆矩阵及其性质且给出加号逆矩阵的一种切实可用的计算方法。MCG算法不同于通常的共轭梯度算法，它不要求系数矩阵正定、可逆或者列满秩，因此总是可行的。本文讨论实中心对称矩阵Moore-Penrose广义中心对称{1, 4}逆问题，将Moore-Penrose广义中心对称{1, 4}逆问题转换为求单变量线性矩阵方程组的中心对称解，然后建立MCG算法求线性矩阵方程组的中心对称解。

定义1 设n阶单位矩阵 I = ( e 1 , e 2 , ⋯ , e n ) ，称n矩阵 S n = ( e n , e n − 1 , ⋯ , e 1 ) 为次单位矩阵。

若矩阵 X ∈ R n × n 满足 S n X S n = X ，称 X 为中心对称矩阵，用 C S R n × n 表示中心对称矩阵集合。

定义2 设 A ∈ R m × n ，若 X ∈ R n × m 满足

A X A = A , ( X A ) T = X A (1)

则称 X 为 A 的一个{1, 4}逆，记为 A ( 1 , 4 ) 。特别地，若 A ∈ C S R n × n ，且 S n X S n = X ，则称 X 为中心对称矩阵 A 的一个中心对称{1, 4}逆。根据广义逆矩阵的性质，矩阵的{1, 4}逆不是唯一的。

令

A 1 = B 1 = F 1 = A , C 1 = O m × m , D 1 = O n × n ,

A 2 = I n , B 2 = A , C 2 = − A T , D 2 = I m , F 2 = O m × m ,

则(1)式可以改为

{ A 1 X B 1 + C 1 X T D 1 = F 1 A 2 X B 2 + C 2 X T D 2 = F 2 (2)

问题的提出：任意给出中心对称矩阵 A ，求方程组(2)的中心对称解 X 和 X T ，则 A 的中心对称广义{1, 4}逆就是矩阵 X 。

本文通过建立求线性矩阵方程组(2)中心对称解的MCG算法，给出了MCG算法的性质证明该算法是收敛的，给出了算法收敛性定理。最后给出了算例，验证了MCG算法求中心对称矩阵 A 的中心对称{1, 4}逆是有效的。

结合共轭梯度算法原理以及中心对称矩阵的特点，建立求矩阵方程组(2)的MCG算法如下：

第1步：任意给定初始矩阵 X 1 ∈ C S R n × n ，置 k : = 1 ，计算

R k ( i ) = F i − A i X k B i − C i X k T D i   ( i = 1 ,   2 ) , R k = ( R k ( 1 ) R k ( 2 ) ) ,

R ˜ k = ∑ i = 1 2 ( A i T R k ( i ) B i T + D i ( R k ( i ) ) T C i ) , P k = 1 2 ( R ˜ k + S R ˜ k S )

第2步：若 R k ( 1 ) 与 R k ( 2 ) 均为零矩阵，停止计算；否则，计算

α k = ‖ R k ‖ 2 ‖ P k ‖ 2 , X k + 1 = X k + α k P k

第3步：计算

R k + 1 ( i ) = F i − A i X k + 1 B i − C i X k + 1 T D i   ( i = 1 ,   2 ) ,

R ˜ k + 1 = ∑ i = 1 2 ( A i T R k + 1 ( i ) B i T + D i ( R k + 1 ( i ) ) T C i ) ,

β k + 1 = ‖ R k + 1 ‖ 2 ‖ R k ‖ 2 , P k + 1 = 1 2 ( R ˜ k + 1 + S R ˜ k + 1 S ) + β k + 1 P k

易见，MCG算法中的矩阵 X k , P k ∈ C S R n × n ，下面给出MCG算法的基本性质，证明MCG算法在有限步计算之后停止。

性质1 MCG算法中的矩阵 R i ， P i ， R ˜ i 满足

t r ( R i + 1 T R j ) = t r ( R i T R j ) − β j − 1 t r ( P i T R ˜ j )

证明 由MCG算法可得

t r ( R i + 1 T R j ) = t r ( ( R i + 1 ( 1 ) R i + 1 ( 2 ) ) T R j ) = t r ( ( F 1 − A 1 X i + 1 B 1 − C 1 X i + 1 T D 1 F 2 − A 2 X i + 1 B 2 − C 2 X i + 1 T D 2 ) T R j ) = t r ( ( R i ( 1 ) − α i ( A 1 P i B 1 + C 1 P i T D 1 ) R i ( 2 ) − α i ( A 2 P i B 2 + C 2 P i T D 2 ) ) T R j ) = t r ( R i T R j ) − α i t r ( ( A 1 P i B 1 + C 1 P i T D 1 A 2 P i B 2 + C 2 P i T D 2 ) T R j )

= t r ( R i T R j ) − α i t r ( ( A 1 P i B 1 + C 1 P i T D 1 A 2 P i B 2 + C 2 P i T D 2 ) T ( R j ( 1 ) R j ( 2 ) ) ) = t r ( R i T R j ) − α i t r ( ∑ k = 1 2 ( B k T P i T A k + D k T P i C k T ) R j ( k ) ) = t r ( R i T R j ) − α i t r ( ∑ k = 1 2 ( P i T A k T R j ( k ) B k + P i C k T ( R j ( k ) ) T D k T ) ) = t r ( R i T R j ) − α i t r ( P i T R ˜ j )

性质2 当 k ≥ 2 时，对于MCG算法中的矩阵 R i ， P i 和 R ˜ i ，有

t r ( R i T R j ) = 0 , t r ( P i T P j ) = 0 , ( i ≠ j ,   i , j = 1 , 2 , 3 , ⋯ , k ) (3)

证明 采用数学归纳法，对于 k = 2 ，由性质1可知，

t r ( R 2 T R 1 ) = t r ( R 1 T R 1 ) − ‖ R 1 ‖ 2 ‖ P 1 ‖ 2 t r ( P 1 T R ˜ 1 ) = ‖ R 1 ‖ 2 − ‖ R 1 ‖ 2 ‖ P 1 ‖ 2 t r ( P 1 T R ˜ 1 + S R ˜ 1 S 2 ) = ‖ R 1 ‖ 2 − ‖ R 1 ‖ 2 ‖ P 1 ‖ 2 t r ( P 1 T P 1 ) = 0

t r ( P 2 T P 1 ) = t r ( ( R ˜ 2 + S R ˜ 2 S 2 + ‖ R 2 ‖ 2 ‖ R 1 ‖ 2 P 1 ) T P 1 ) = t r ( R ˜ 2 + S R ˜ 2 S 2 P 1 ) + ‖ R 2 ‖ 2 ‖ R 1 ‖ 2 t r ( P 1 T P 1 ) = t r ( R ˜ 2 P 1 T ) + ‖ R 2 ‖ 2 ‖ R 1 ‖ 2 ‖ P 1 ‖ 2 = ‖ P 1 ‖ 2 ‖ R 1 ‖ 2 ( t r ( R 1 T R 2 ) − t r ( R 2 T R 2 ) ) + ‖ R 2 ‖ 2 ‖ R 1 ‖ 2 ‖ P 1 ‖ 2 = − ‖ P 1 ‖ 2 ‖ R 1 ‖ 2 t r ( R 2 T R 2 ) + ‖ R 2 ‖ 2 ‖ R 1 ‖ 2 ‖ P 1 ‖ 2 = 0

假设 k = s   ( s ≥ 2 ) 时(3)式成立，则当 k = s + 1   ( s ≥ 2 ) 时，由性质1可得

t r ( R s + 1 T R s ) = t r ( R s T R s ) − ‖ R s ‖ 2 ‖ P s ‖ 2 t r ( P s T R ˜ s ) = t r ( R s T R s ) − ‖ R s ‖ 2 ‖ P s ‖ 2 t r ( P s T R ˜ s + S R ˜ s S 2 ) = ‖ R s ‖ 2 − ‖ R s ‖ 2 ‖ P s ‖ 2 t r ( P s T ( P s − ‖ R s ‖ 2 ‖ R s − 1 ‖ 2 P s − 1 ) ) = ‖ R s ‖ 2 − ‖ R s ‖ 2 ‖ P s ‖ 2 ( t r ( P s T P s ) − ‖ R s ‖ 2 ‖ R s − 1 ‖ 2 t r ( P s T P s − 1 ) ) = ‖ R s ‖ 2 − ‖ R s ‖ 2 ‖ P s ‖ 2 t r ( P s T P s ) = 0

t r ( P s + 1 T P s ) = t r ( ( R ˜ s + 1 + S R ˜ s + 1 S 2 + ‖ R s + 1 ‖ 2 ‖ R s ‖ 2 P s ) T P s ) = t r ( ( R ˜ s + 1 + S R ˜ s + 1 S ) T 2 P s ) + ‖ R s + 1 ‖ 2 ‖ R s ‖ 2 t r ( P s T P s ) = t r ( R ˜ s + 1 P s T ) + ‖ R s + 1 ‖ 2 ‖ R s ‖ 2 t r ( P s T P s ) = 0

当 j = 1 时，有

t r ( R s + 1 T R 1 ) = t r ( R s T R 1 ) − ‖ R s ‖ 2 ‖ P s ‖ 2 t r ( P s T R ˜ 1 ) = 0 − ‖ R s ‖ 2 ‖ P s ‖ 2 t r ( P s T R ˜ 1 + S R ˜ 1 S 2 ) = − ‖ R s ‖ 2 ‖ P s ‖ 2 t r ( P s T P 1 ) = 0

当 j = 2 , 3 , 4 , ⋯ , s − 1 时，有

t r ( R s + 1 T R j ) = t r ( R s T R j ) − ‖ R s ‖ 2 ‖ P s ‖ 2 t r ( P s T R ˜ j ) = 0 − ‖ R s ‖ 2 ‖ P s ‖ 2 t r ( P s T R ˜ j + S R ˜ j S 2 ) = − ‖ R s ‖ 2 ‖ P s ‖ 2 t r ( P s T ( P j − ‖ R j ‖ 2 ‖ R j − 1 ‖ 2 P j − 1 ) ) = − ‖ R s ‖ 2 ‖ P s ‖ 2 ( t r ( P s T P j ) − ‖ R j ‖ 2 ‖ R j − 1 ‖ 2 t r ( P s T P j − 1 ) ) = 0

t r ( P s + 1 T P j ) = t r ( ( R ˜ s + 1 + S R ˜ s + 1 S 2 + ‖ R s + 1 ‖ 2 ‖ R s ‖ 2 P s ) T P j ) = t r ( R ˜ s + 1 P j T ) + ‖ R s + 1 ‖ 2 ‖ R s ‖ 2 t r ( P s T P j ) = t r ( R ˜ s + 1 P j T ) = t r ( P j T R ˜ s + 1 ) = ‖ P j ‖ 2 ‖ R j ‖ 2 ( t r ( R j T R s + 1 ) − t r ( R j + 1 T R s + 1 ) ) = 0

根据矩阵迹的性质可得

t r ( R j T R s + 1 ) = t r ( R s + 1 T R j ) = 0 , t r ( P j T P s + 1 ) = t r ( P s + 1 T P j ) = 0

所以当 k = s + 1 时(3)式也成立。由归纳法原理可得 1 ≤ j < i ≤ k 时(3)式成立。

性质3 设 X 是问题的任意一组解，则由算法得到的矩阵 X ( k ) ， P k 及 R k 满足

t r ( P k T ( X − X ( k ) ) ) = ‖ R k ‖ 2   ( k = 1 , 2 , 3 , ⋯ ) .

定理1 若中心对称矩阵 A 存在中心对称{1, 4}逆，则对任意初始矩阵 X 1 ∈ C S R n × n ，MCG算法可在有限步计算后得到 A 中心对称{1, 4}逆。若方阵 A 无中心对称{1, 4}逆，则在MCG算法中存在正整数k，使得 R k ≠ O 且 Q k = O 。

下面给出两个算例，在例1中，给出两个方阵说明文中建立的MCG算法是可行的。在例2中给出不可逆方阵 A 在不同阶数下计算结果，说明MCG算法具有很高的效率，且都能求出矩阵方程组(2)广义中心对称{1, 4}逆解。

文中计算时间(秒)、矩阵阶数、MCG算法迭代次数、MCG算法中断次数、实际误差的范数依次用time、n、 k 1 、 k 12 和error表示。为避免迭代次数过分增加，设定MCG算法迭代次数上限为2999次，MCG算法终止准则为 ε = 10 − 12 ，初始矩阵 X 1 均为零矩阵。

例1 用MCG算法求中心对称矩阵 A 的一个广义中心对称{1, 4}逆，在本例中，中心对称矩阵 A 分别取可逆方阵和不可逆方阵。结果如下

1) 若 A = ( 6 5 13 3 − 4 3 13 5 6 ) ， A ( 1 , 4 ) = A − 1 = ( − 0.0526 0.0472 0.0903 0.0283 − 0.1792 0.0283 0.0903 0.0472 − 0.0526 ) 。

2) 若 A = ( 5 3 5 7 8 7 5 3 5 ) ，则 A ( 1 , 4 ) = X 1 = ( 0.1053 − 0.0789 0.1053 − 0.1842 0.2632 − 0.1842 0.1053 − 0.0789 0.1053 ) 。

从(1)中可以看出，当 A 可逆时，由MCG算法求得的广义中心对称{1, 4}逆就是 A 的逆矩阵；从(2)中可以看出，当 A 不可逆时，则MCG算法能在2步迭代计算求得 A 广义中心对称{1, 4}逆。

例2 取矩阵 A 如下，按MCG算法计算步骤求得矩阵方程组(2)的广义中心对称{1, 4}逆，结果如表1 (Matlab软件2016版-PIV3.0GHZ微机)：

A ˜ = ( a i j ) ，其中 a i j = ( i − j ) 2 ,   ( i , j = 1 , 2 , ⋯ , n ) ， A = A ˜ + S A ˜ S

表1. 方程组(2)的广义中心对称{1, 4}逆计算结果

从表1结果可以看出，文中建立的MCG算法求广义中心对称{1, 4}逆具有很高的计算效率。

广义逆求解算法多种多样，如对矩阵 A 满秩分解、奇异值分解等，也可以建立求广义逆递归计算公式，各种算法在满足某些条件下均能求出矩阵的广义逆。本文建立的求中心对称矩阵{1, 4}逆的MCG算法无意与其他算法做比较，只是给出一种求中心对称矩阵中心对称{1, 4}逆的迭代算法，该算法使用范围广，对矩阵 A 限制条件少，修改算法中某些矩阵，也可以建立方阵其他特殊解的MCG算法。

2019年福建省教育厅中青年教师教育科研项目：(JAT190410)。

陈世军. 中心对称矩阵的广义中心对称{1, 4}逆的迭代算法Iterative Algorithm for Generalized Centrosymmetric {1, 4} Inverse of Centrosymmetric Matrix[J]. 应用数学进展, 2021, 10(04): 1032-1038. https://doi.org/10.12677/AAM.2021.104112