求解块下Hessenberg方程组的新算法.docVIP

求解块下Hessenberg方程组的新算法 　　A New Algorithm of Solving Block Lower Hessenberg Equations 　　 　　Wei Yanhong； Pan Lijing 　　（Weinan Normal University，Weinan 714000，China） 　　 　　摘要： 给出了求解块下Hessenberg方程组的新算法――线性插值法，该方法所需要的运算量较通常的Gauss消去法所需要的运算量小，而且有很高的计算精度． 　　Abstract: A new algorithm of solving block lower Hessenberg Equations is proposed, which is named linear interpolation method. Compared with direct method（Gauss expunction method）, the operational quantity of this method increases slower and the precision is higher. 　　关键词： 块下Hessenberg矩阵 线性代数方程组 线性插值法 　　Key words: blocked lower Hessenberg matrix；linear algebraic equations；linear interpolation method 　　中图分类号：G42文献标识码：A文章编号：1006-4311（2011）２7-0155-02 　　0引言 　　在工程领域，经常要求解以块状下Hessenberg矩阵 　　为系数矩阵的大型线性方程组Ax=f， A为可逆方阵，Ai，j是r阶方阵，Ai，i+1是r阶可逆方阵，xi和fi为r维列向量。采用迭代方法[1-3]求解这类问题时，需要讨论算法的收敛性；采用直接方法（Gauss消去法）求解时，运算量会随方程组阶数的增大而快速增加。本文基于块下Hessenberg方程组建立了一种线性插值法，该方法所需要的运算量较Gauss消去法所需要的乘除法运算量小，而且有很高的计算精度。 　　1算法推导 　　方程组Ax=f的m个子方程为 　　A■x■+A■x■=f■A■x■+A■x■+…+A■x■=f■（i=2，…，m-1）A■x■+A■x■+…+A■x■=f■（1） 　　令 　　则方程组■x=■等价于 　　A■x■+A■x■=f■A■x■+A■x■+…+A■x■=f■（i=2，3…，m-1）（2） 　　任意给定x1，由1式的第1个方程解出x2，由第i个方程解出xi+1（i=2，3…，m-1），可得方程组■x=■的一组解。 　　分别取 　　x■■=00■0，x■■=10■0，…，x■■=01■0，…，x■■=00■1 　　可求得方程组■x=■的r+1组解 　　x■=x■■x■■■x■■，x■=x■■x■■■x■■，x■=x■■x■■■x■■，…，x■=x■■x■■■x■■， 　　令y■=x■-x■，即y■■=x■■-x■■（i=1，…，m），则有■y■=0（k=1，2，…，r）任意选取参数t1，t2，…，tr，构造向量 　　x■=x■+t■y■+…t■y■=x■+y■，…，y■t■t■■t■（3） 　　易知■x■=■。要求x*满足1式的第m个方程，可得 　　A■x■■+（y■■，…，y■■）t■t■■t■+A■x■■+y■■，…，y■■t■t■■t■+…+A■x■■+y■■，…，y■■t■t■■t■=f■ 　　整理得 　　A■y■■，…，y■■+A■y■■，…，y■■+…+A■（y■■，…，y■■）t■t■t■=f■-A■x■■+Am，2x■■+…+A■x■■（4） 　　若4式有解t1，t2，…，tr，则由3式确定的x*是方程组Ax=f的唯一解，称这种方法为解块下Hessenberg方程组的线性插值方法。 　　2运算量分析 　　给定x1，求解第i个方程约需■r■+r■+ir■次运算，求得方程组■x=■的一组解约需（m-1）（■r■+r■）+■r■次运算，r+1组解约需（r+1）■r■（m-1）+■r■次运算，求解4式约需■+mr■+（m+1）r■次运算，由3式计算x*约需mr2次运算。总的运算次数大约为 　　■m■r■（r+1）+■mr■（4r■-2r+6）-■r■（r■-3r-6） 　　直接求解方程组Ax=f（Gauss消去法）约需■（mr）■+（mr）■次乘除法运算。 　　当r=3，4，6，9时，分别取m=100和m=1000，线性插值法与直接法所需的乘除法运算次数对比见表1和表2。