Toeplitz 方程组的拟对称化快速算法.pdfVIP

第 31卷 第 5期 周 口师范学院学报 2014年 9月 Vo1．31No．5 JournalofZhoukouNormalUniversity Sep．2014 Toeplitz方程组的拟对称化快速算法 王 健 ，朱 变 ，徐 仲 (1．周 口师范学院计算机科学与技术学院，河南 周 口466001； 2．西北工业大学 应用数学系，陕西 西安 710072) 摘 要：利用 由 阶 Toeplitz矩 阵构造 的 2n阶拟对称矩 阵，给 出了求解 Toeplitz方程组的拟对称化快速算 法．该算法与 已有算法相 比，或者减少 了计算量 ，或者提高了计算精度． 关键词 ：Toeplitz矩 阵；线性方程组 ；位移结构；快速算法 中图分类号 ：O241．6 文献标志码 ：A 文章编号 ：1671—9476(2014)05—0008一O4 Toeplitz矩阵及 Toeplitz方程组的求解在线性预测 、数字信号处理、谱分析 、误差控制码、控制论、系 统识别及系统理论等领域内起着重要的作用口 ]．在求微分方程数值解和三次样条插值 中常需求解 以三 对角 Toeplitz矩阵为系数矩阵的线性方程组．1985年，Kumar将 Toeplitz矩阵镶嵌成一个大的循环矩阵， 进而得到求解 Toeplitz方程组 的快速算法 ]．1986年，Ben—Artzi给出求非奇异 Toeplitz矩阵逆 的显示表 达式，间接得到了求解 Toeplitz方程组的快速算法 ]．之后 1986到 1995年，Gohberg，Hsue等人相继给 出了求 Toeplitz方程组的快速算法 ]．笔者通过构造特殊分块矩阵，间接得到了求解线性方程组 Tx=b 的解的快速算法 ，该算法所需计算量为 O(n)． 1 算法推导 考虑求解 Toeplitz线性方程组 Tx=b，其中T：：=(￡)， 是 n阶Toeplitz矩阵，X一 (z ”， ) ， b一 (6，… ，b)T．容易验证 Toeplitz矩阵 T满足 T一 ===pel +elq ， (1) 其中P一 (￡。，t一 一，t- )，口： (O，t 一， )，el是第 i个分量为 1其余分量为0的 维列 向量 ， 一 (P ，el，…，：，0)是 阶移位矩阵． 构造 2 阶方阵 M ： fLTT]， (2) IT 0J 如果求得线性方程组 My— f (3) 的解向量Y一 ( ，Y ，…，Y ) ，这里f一 (0 ，b)，则Y的前n个分量构成的列向量即是Toeplitz方程组 一 b的解 向量．利用式(1)易知式(2)的矩阵M 满足 M 一 (]Mf ]一骞 ， c4 其中g — g =hT—hT一 (P{” ，O)，g —h4~一 (O ，P )，g ===hT一 (0 ，e’)，g —hT一 (口，0)． 引理Ⅲ 设 阶矩阵A一 (a) 的顺序主子阵A (愚一 1，2，…，)均可逆．给定线性方程组Ax—b， 其 中X一 (z1，…， )，b一 (6 ”，b)．设 b 一 (6 一，b) ，又设线性方程组A女 =：=b 和A女U 一 P 收稿 日期 ：2014一O2—28 基金项 目：国家 自然科学基金资助项 目(No；周 口师范学院青年科研基金 (No．ZKNUC0204) 作者简介：王健 (1984一)，男，河南郸城人，助教，硕士，主要从事数值计算、数值代数研究． 第 31卷 第 5期 王 健 ，等：Toeplitz方程组的拟对称化快速算法