几类结构矩阵的快速算法及其应用-计算数学专业论文.docxVIP

摘 摘 要 众所周知，在工程计算和实际应用中有许多问题最终都归结为矩阵计算 问题，而且不同的应用会导出一些具有特殊结构的矩阵计算．最常见的一些结 构矩阵有Toeplitz矩阵k-Jl，Hankel矩阵k+Jl，Vandermonde矩阵陋一-11， Cauchy矩阵【士】等等．处理与这些结构矩阵有关的矩阵计算问题(例如， 求解线性方程组，计算特征值等)，若矩阵的阶数较小时，通常的经典算法是 可行的(例如LU分解算法．QR算法等)．然而，在许多实际应用当中，矩 阵的阶效n很大∽一108一109)或某个线性方程组需要多次计算直到得到一 个满意的结果(例如用迭代法时)，此时这些经典的算法由于代价太大而失去 了实际意义． 因此，针对这些结构矩阵的特点而设计一些能利用它们的结构的，数值稳 定的快速算法，具有非常重要的意义．正因为结构矩阵在实际应用中所具有的 重要意义，国内外众多的学者将目光投入到这—领域．结构矩阵的快速算法中 最著名的莫过于快速傅里叶变换(即FF丁)，有许多快速算法均是由快速傅里 叶变换导出的．因此，著名数学家Charles Van Loan曾这样评价快速傅里叶 变换算法t。从计算的角度看，快速傅里叶变换是本世纪最杰出的成就之一， 毫不夸张地说，快速傅里叶变换改变了科学与工程计算的面貌，如果没有它， 生活将会是另一种景象”． 本论文主要研究了Toeplitz矩阵．Hankel矩阵、Pascal矩阵以及合流 Cauchy—Vandermonde矩阵的一些性质及相关的快速算法，同时还给出这 些快速算法的数值实验和在一些问题中的应用．理论和数值实验显示，这些快 速算法是行之有效的． 第一章，我们简单介绍了研究结构矩阵快速算法的现实意义、研究概况以 及常用的研究方法，同时也给出了与本论文有关的几类结构矩阵的定义及其 基本性质． 在第二章和第三章，我们主要是利用Toeplitz矩阵和Hankel矩阵的特 殊结构，导出相应的递推关系式，然后再利用快速傅里叶变换(FFr)，给出 了计算Toeplitz矩阵的正弦变换和Hankel矩阵的余弦变换的快速算法(算 法计算复杂度为O(ntogn))．该算法不仅快而且存贮有效，因为在执行该快速 算法的过程中，不需要存贮任何矩阵．同时在第二章中，我们还给出了该快速 算法在利用Jacobi旋转变换计算Toeplitz矩阵的特征值中的应用．数值实验 显示。对变换后的矩阵利用Jacobi旋转变换方法求其特征值所需要的扫描次 显示。对变换后的矩阵利用Jacobi旋转变换方法求其特征值所需要的扫描次 数大约为对原矩阵利用Jacobi旋转变换方法求其特征值所需要的扫描次数的 三分之一． 第四章，我们研究了Pascal矩阵的一些性质，并且给出了计算系数矩阵 为Pascal矩阵的线性方程组的解的一个快速算法，该算法的计算复杂度为 O(ntogn)，比文[221所给的快速算法运算量要少([221中所给的快速算法的计 算复杂度为O(no))．此外，我们还将我们所设计的快速算法用到了常系数非 同态常微分方程的求解中． 第五章，我们讨论了一类在有理插值问题中有着广泛应用的矩阵，即合流 Cauchy—Vandermonde矩阵．利用块Neville消去，给出了合流Cauchy— Vandermonde矩阵及其逆矩阵的块双对角分解，推广了文[521的结果． 第六章，我们讨论了K．R．Drieasel在文【19】中所提出的—个问题，即 Problem A：任给n个实数，若它们是一个n阶实Hankel矩阵的特征 值，则这n个实数必须满足何条件? 对Problem A，我们给出了几个必要条件，并且得到了Problem A的必要 条件的一个通用表达式。推广了文【26l的部分结果． 关键词：结构矩阵，快速算法，Toeplitz矩阵，Hankel矩阵，Pascal矩 阵，合流Cauchy—Vandermonde矩阵． II ABSTRACTMany ABSTRACT Many important problems in engineering and applied sciences ultimately can be reduced to matrix computation problems．Moreover，varioll8 appli- cations often introduce a special structure into the corresponding matrices． Among classical examples are Toeplitz matrices【ai—J】，Hankel matrices k钾】， Vandermonde matrices【a,S。1】，Cauchy matrices【六】