模拟退火算法在矩阵实特征值求解中的应用.pdfVIP

模拟退火算法在矩阵实特征值求解中的应用 梁芳，祝庆 武汉理工大学理学院，武汉（430070 ） E-mail ：liangfang4022@ 摘 要：给出了矩阵特征值的定义及确定特征值范围的圆盘定理。基于模拟退火算法提出了 一种新算法，新算法用于解决实特征值的求解问题。具有通用性，算法的稳定性也大大提高。 数值实验表明，新算法的求解精度有了明显的提高。 关键词：特征值问题；圆盘定理；模拟退火算法 中图分类号：TP301.6 文献标志码：A 1. 引言 物理学、力学和工程技术中的很多问题在数学上都归结为求矩阵的特征值问题。许多 学者已经对此进行了长期的研究和尝试工作，常用的方法主要有幂法、子空间迭代法、 Lanczos 方法等。幂法只能求得端特征值，且算法的稳定性依赖于特征值的分布情况，当模 最大特征值和模次大特征值的模接近时收敛得非常慢；子空间迭代法是幂法的直接推广，收 敛速度有所提高，但其它缺陷依然得不到弥补；Lanczos 方法是目前被认为是求解高阶矩阵 特征值的最有效方法之一，但它对随机初始向量的选取以及舍入误差敏感，正交性受损,算 法稳定性不太理想。本文试用简单、快速、稳定性强的模拟退火算法来求解特征值问题。 2. 矩阵特征值的定义及相关结果 设矩阵A (a ) ∈M (R) ，如果存在λ∈R ，X ∈R n ，且X ≠0 满足AX λX ，则 ij n λ称为方阵 A 的特征值。X 为对应λ的特征向量。 定义 1 设矩阵A (a ) ∈M (R) ，称R a 为矩阵A 的第 行半径，以a 为圆心，R ij n i ∑ ij i ii i j ≠i j 1 为半径的圆称为矩阵的第i 行圆盘s ，即s {λ λ−a ≤R } 。 i i ii i 定义 2 设矩阵A (a ) ∈M (R) ，称R ∑a 为矩阵 A 的第i 列半径，以a 为圆心，R ij n i ij ii i j ≠i i 1 为半径的圆称为矩阵的第 列圆盘 ，即s {λ λ−a ≤R } 。 i s i i ii i 定理 1 若圆盘的第 行圆盘与第 列圆盘至少有一个是孤立圆盘，则有一个特征值属于它们 i i 的交集中，即半径取行半径与列半径的较小的一个。 k n