加权Toeplitz最小二乘问题的预处理算法.pdfVIP

第一章 引言 §1．1 研究背景和意义 在本论文巾，我们丰要考虑对加权Toeplitz矩阵最小二乘I、口J题的预处理迭代 算法．考虑下面的最d-乘问题 IIAz_6眩 (1．1) ra工in 其中矩阵A和右端项b分别为 A=阱6=阱 (1．2) X 这里D是mm阶对角矩阵，对角线上的元素都是正数，，是佗×礼阶单位矩阵． K是…个m×礼阶Toeplitz矩阵且是列满秩的，，是。。个给定的右端项，p0是 一个正则化参数【25】．这类问题通常出现在图像还原【15】和非线性图像恢复…II． 在这些应川中，问题的规模会很大，很容易上百万．但很多时候我们必须在有限 的时间里解决这些问题，冈此必须寻求有效的快速算法米求解这类I’口J题． 对于一般的Toeplitz线性方程组 Kx=b． (1．3) (1．3)的快速Cholesky分解法，也称Schur方法．以上两种直接算法的时问复杂 矩阵，这些算法均不适用．因为加权Toeplitz矩阵本身的性质和空间的特性，使得 A的置换秩会很大『l71．所以在求解这类加权Toeplitz最小二乘问题时，无法使 用快速直接算法，这时我们可以使用迭代算法，如Krylov子空间迭代算法． 我们知道，线性最d、乘问题可以化为等价的正规方程，然后用共轭梯度法 来求解f51．共轭梯度法在理论上是直接法．若方程阶数是n，则共轭梯度法的第 n步迭代解z。必定是精确解．但实际上，由于舍入误差的影响，Xn往往不能满足 原方程．另外，对于一般的迭代法(包括共轭梯度法)，当方程组的系数矩阵的条件 2 数较大或者谱分布较差时．迭代算法的收敛速度通常很慢．为了改进这一个问题， 人们提出了预处理技术，即在方程两边同时乘以一个非奇异矩阵的逆，然后再用 迭代法米求解．近年来．人们不断提出各利，新的预处理技术．预处理技术已经成 为解大型系数方程组的一个重要手段．但是由于预处理子与问题本身密切相关， 针对不同的l’日J题，预处理子也各不相同．冈此如何构造有效的预处理子是F1前数 值计算领域II I的热门研究课题之一． 子有着特殊的要求，即求解预处理系统时的运算量不能超过D(几logn)．F1前常川 的一类预处理子为循环矩阵．因为循环矩阵可以通过Fourier矩阵对角化【14】，所 以求解一个以循环矩阵为系数矩阵的线性方程组时，运算量为0(礼log礼)．对于 和T．Chan循环预处理子等．有学者将这些预处理子推广到了求解Toeplitz矩阵 最t]、--乘I、日J题，对于置换秩较小的情形，取得了较好的效果f8，9，,01．但对于置换 秩较大的加权Toeplitz矩阵最小二乘问题，几前研究结果并不多【7】． §1．2 本文的研究内容和论文框架 本论文首先将加权Toeplitz最小_乘f’日J题(1．1)转化为等价的鞍点问题．然 后研究基于对称与反对称分裂(Hs)的预处理子的构造及其性质，并通过数值算 例来说明这些预处理子的有效性． 本论文主要做了以下工作： (1)通过引入三个参数矩阵，构造出了基于Hs分裂的预处理子的一般形式． (2)通过对这三个参数矩阵的不同选取，构造出了三类有效的预处理子，并对这 些预处理子的性质进行了理论分析，给出了预处理后的系数矩阵的特征值分 布情况． (3)对提出的预处理子进行了数值算例，验证了新的预处理子的有效性． 本文主要结构如下：在第二章，我们对Toeplitz矩阵和循环矩阵作了介绍．在 第三章，我们给出了基于HS分裂的预处理子，并进行了理论分析．最后，在第四 章，我们给出了数值试验的结果． 第二章Toeplitz矩阵和循环矩阵 设An∈RnX-，如果‰可写成下面的形式 qji ％Ⅶ ～咄 知 ¨ ‰哪 ‰ = 鲫们；