波动方程辛算法的迭代求解  

ITERATIVE METHODS FOR SYMPLECTIC ALGORITHM OF WAVE EQUATION

作　　者：蒋长锦[1] 

机构地区：[1]中国科学技术大学数学系,安徽合肥230026

出　　处：《计算物理》2002年第1期13-16,共4页Chinese Journal of Computational Physics

摘　　要：对 2x2 利用中心差商算子 ,对expt作对角Pad啨逼近 ,由波动偏微分方程可得到两类具有O(Δx2 +Δt2l)和O(Δx4+Δt2l)精度的辛格式 .对由此类辛格式产生的线性方程组构造了两种迭代解法 ,并对l=1,2 ,3,4给出了它们的收敛条件 .By using the central difference quotient operator for  2x 2 and the diagonal Padé approximation of exp t , two kinds of symplectic schemes which have accuracy O( Δ x 2+ Δ t 2l ) and O( Δ x 4+ Δ t 2l ) , respectively, can be attained for wave partial differential equation. Two iterative methods are described for the linear systems formed from the above schemes. Their conditions of convergence are given for l =1,2,3,4. The numerical experiments demonstrate that the symplectic algorithm have efficiency and both methods are convergent. [

关 键 词：HAMILTON系统 辛差分格式 迭代解法 收敛条件 波动方程 

分 类 号：O413.1[理学—理论物理] O241.6[理学—物理]