黄金比例分割法确定对称逐次超松弛迭代法最佳松弛因子.pdfVIP

黄金比例分割法确定对称逐次超松弛迭代法的最佳 松弛因子 张德宣，杜成斌，孙立国 河海大学工程力学系，南京 (210098) E-mail：zdx.1234@163.com 摘 要：本文提出了将黄金分割法确定松弛因子与对称逐次超松弛法的改进迭代格式相结合 的迭代算法。算法应用黄金比例分割法确定最佳松弛因子，成功的将其与运行速度和效率很 高的对称逐次超松弛法的改进迭代格式相结合，并给出了迭代收敛性证明，编写了相应的程 序，对一实际结构进行的算例计算表明，与大型商业软件的计算结果相比较，本文所提算法 具有精度高，收敛快的优点。 关键词：黄金分割；对称逐次超松弛法；松弛因子 1．引言 线性代数方程组Ax b ( 其中A 为非奇异系数矩阵,b 为列向量) 的迭代解法通常有 Jacobi迭代法,Gauss Seidel迭代法和对称逐次超松弛迭代（SSOR ）法三种，此外还有USAOR[1] 迭代法，USSOR[2]迭代法等。 但是当结构的自由度达到较大的数目（如10000个自由度以上），Jacobi法和Gauss Seidel 法存在着运算速率很慢这一明显缺点，而且这些迭代方法的结果的精度不能保证。因此本文 针对运算速率很高的对称逐次超松弛法的改进迭代格式提出了用黄金分割法确定松弛因子 的思路，编写了黄金分割计算松弛因子与相应的迭代法相结合的程序，实际算例表明，本文 方法可提高求解大型线性方程组的求解效率和精度。 2 ．SSOR-PCG法改进的迭代格式 SSOR法在每步迭代中，需计算系数矩阵与方向向量的乘积，下面引用了该法的改进迭 代格式，避免了该乘积的计算，比原迭代格式节省计算量8 ％—50 ％。 SSOR法的改进迭代格式[3]为： 置初值 0 0 0 −1 0 0 0 0 −T 0 δ , g Kδ −R ,y W g , z −Vy , d W z ,k 0 ， (1) 0 R ： ( k k ) α y , Vy 式中( ， 为内积表示，下同。) 如果α ≤ε ，则停止，否则 -1- ⎧τk (y k , Vyk ) /(d k ,2z k −Vdk ), ⎪ k +1 k k δ δ τ d , ⎪ + k ⎪ k +1 k k −1 k k y y +τk (d +W (z −Vd )), ⎪ ⎪⎪βk (y k +1 , Vyk +1 ) /(y k , Vyk ), ⎨ (2) z k +1 −Vyk +1 +β z k , ⎪ k ⎪ k +1 −T k +1 d W z ⎪