关于Toeplitz矩阵问题的数值解法.pdf

摘 要 Toeplitz 矩阵在科学与工程的众多领域中均有着广泛的应用, 比如数字图 像与信号的处理, 微分方程数值解以及排队网络等等. 本文考虑两个与Toeplitz 矩阵有关的问题, 对称 Toeplitz 系统上的弱非线性方程组和空间分数阶非线性 Schrödinger 方程. 具体如下: 对于含对称Toeplitz 系统的弱非线性方程组, 通过分离线性项与非线性项, 利 用逆无关的预处理共轭梯度法（AIPCG）, 建立了以Picard 迭代法作为外迭代的 Picard-AIPCG 迭代方法. 其优势在于无需精确计算和存储雅各比矩阵, 只需解常 系数矩阵的线性子系统. 因此在实际应用中, 大大缩减了计算量和存储量. 而且利 用AIPCG 迭代作为内迭代的收敛速度非常快且不依赖于参数. 理论分析证明在 给定条件下该方法是全局收敛的. 数值实验表明在特定情形下, Picard-AIPCG 方 法是可行且有效的. 对于空间分数阶非线性Schrödinger 方程, 使用隐式守恒差分离散后, 原问题 转化成线性方程组, 其系数矩阵为非负对角阵、对称正定的Toeplitz 矩阵和复的 单位阵之和. 基于交替迭代方法, 本文提出了不对称的分裂(LS）迭代方法和改进 的平均移位分裂(MMS) 迭代方法, 并证明了其收敛性. 此两种迭代方法均需分别 求解以对角矩阵和Toeplitz 矩阵为系数矩阵的子系统, 从而可分别直接求解和快 速求解上述两个子系统. 此外, MMS 迭代方法不需要额外计算最优参数. 最后数 值例子验证了这两种迭代方法的可行性与高效性. 关键词: Toeplitz 矩阵, 弱非线性方程组, 分数阶, 非线性Schrödinger 方程, 数值 解法 i ABSTRACT Toeplitz matrix has been widely used in many fields of science and engineering, such as digital image and signal processing, numerical solutions of differential equa- tionsandqueuingnetwork,etc. Inthispaper,wewillconsidertwoproblemsrelatedto Toeplitz matrix, one is the weakly nonlinear equations on symmetric Toeplitz system, the otheristhespatialfractionalSchrödingerequations,detailsareasfollows: For solving symmetric Toeplitz systems of weakly nonlinear equations, based on the fact that the linear term is the strongly dominant to the nonlinear term, by using theapproximateinverse-freepreconditionedconjugategradient(AIPCG)iterationtech- niqueandPicarditerativemethodasouteriterativemethod,thePicard-AIPCGiteration method has been established. The advantage of this method is that the storage and the accurate computation of Jacobian matrix aren’t necessary and the only need is to solve theconstantToeplitzsub-systems,whichmaybeagreatsavingincomputerstorageand workload when facing the practical problems. Moreover, the inner iteration–AIPCG method converges very