关于Toeplitz矩阵计算_11_15.pdf

第二章 Toeplitz 矩阵的逆矩阵的研究 2.1 引言 Toeplitz 矩阵在数值分析、优化理论、概率统计、自动控制、数字信号处理、 系统辨识、最小二乘估计、自回归滤波器设计等众多科学领域中都有广泛的应用， 是应用最广泛的特殊矩阵之一。由Toeplitz 矩阵的定义可以看出，Toeplitz 矩阵T 是 特殊的次对称矩阵，若T 是可逆的，由次对称矩阵的性质可知T −1 仍然是次对称矩 阵，但是它一般不再是 Toeplitz 矩阵，因此求Toeplitz 矩阵的逆矩阵有一定的困难， 所以，Toeplitz 矩阵的逆矩阵受到不少研究者的关注。 早在 1964 年和 1969 年 Trench 和 Zohar 就给出了一种 Trench-Zohar 算法[1， 2]求 Toeplitz 矩阵的逆矩阵，这种算法的前提条件是 Toeplitz 矩阵T 的各阶顺序主 子阵T (k 1,2, , n) 均非奇异，Zohar 将其推广到强飞奇异（各阶顺序主子阵均非 k 奇异）Toeplitz 矩阵，Trench-Zohar 算法还适合于求分块 Toeplitz 矩阵T 的逆矩阵。 在 1973 年，Akaike 做了进一步改进，得到了分块 Toeplitz 矩阵T 的逆矩阵的Akaike 算法[3]。这些内容将在本章第二节中阐述。 Toeplitz 矩阵的逆矩阵一般不再是 Toeplitz 矩阵，但是，其逆矩阵可以表示成 一些三角 Toeplitz 矩阵的乘积。1972 年，Gohberg 和 Semencul 在俄罗斯提出来了 这个结论，并给出了 Gohberg-Semencul 公式[4]。该Gohberg-Semencul 定理表明， 当 Toeplitz 矩阵T 及n −1阶顺序主子阵Tn−1 可逆时，Toeplitz 矩阵T 的逆矩阵T −1 由 它的第一列及第n 列的2n 个元素确定。当年，Gohberg 和 Krupnik 也提出了一个类 似的结论[5]：只要Toeplitz 矩阵的逆矩阵T −1 的第一列的最后一个分量不为零，逆 矩阵T −1 由它的第一列及第二列的元素确定。在 2004 年，Gohberg、Kaashoek 和 Schagen 又给出了 Gohberg-Semencul 定理的另外两种证明[9]，一种证明用的是 Schur 补和 LU 分解来证明的；另外一种证明利用状态空间技术。因为 Gohberg-Semencul 定理和 Gohberg-Krupnik 定理分别要求逆矩阵T −1 第一列的第一 个分量、最后一个分量不为零。如果此分量恰好为零，就无法将逆矩阵T −1 表示成 一些三角 Toeplitz 矩阵的乘积之和。因此，在当时，有些 Toeplitz 矩阵的逆矩阵根 本没办法由它的两列确定。之后，研究者就沿着 Gohberg、Semencul 和 Krupnik 这 个方向继续往下做，试图解决这个问题。1984 年，Heinig-Rost 提出了 Toeplitz 矩 5 阵的逆矩阵T −1 可以由T −1 的第一列n 元素和 T 的第一行n −1元素确定，并用这 2n −1个元素形成的三角 Toeplitz 矩阵的乘积之和表示逆矩阵T −1 [6]。并且，引入 T-Bezout 矩阵的概念，利用 T-Bezout 矩阵的性质也可以给出 Heinig-Rost 定理完整 的证明，还可以得出 Heinig-Ros 公式。而且，得出非常著名的结论：任意一个可 逆的n 阶 Toeplitz 矩阵的逆矩阵T −1 都是 T-Bezout 矩阵。 在 1986 年，Ben-Artzi 和 Shalom 得到一个结论[7]：只要适当地选取 Toeplitz 矩阵之逆矩阵T −1 的三列，就可以构造出任意一个Toeplitz 矩阵之逆矩阵T −1 ，并给 出了这三个列向量和逆矩阵T −1 的关系，即逆矩阵T −1 可以表示成由这三个列向量 形成的三角 Toeplitz 矩阵的乘积之和，他们还得到了一些非常经典的推论。既然一 个n 阶 Toeplitz 矩阵是由2n −1个元素确定的，Ben-Artzi 和 Shalom 的定理里的三 个列向量（3n 个元素）构造逆矩阵T −1 能否减到2n −1个元素来构造逆矩阵T −1 ， 而且给出T −1 具体的表达式呢？这时，还存在这样一些问题：最少需要多少个T −1 的 元素才能构造出