摘要
首先将一维人工地震的数学模型化为偏微分方程组,离散此方程组的空间偏导数得到一个Hamilton系统,再用梯形公式,它对应ex的对角Pade逼近,来求该系统的数值解即得到这个问题的辛差分算法.证明了该算法的稳定性,给出了正演计算实例,及计算解和理论解的比较.
A Hamilton system is obtained by converting the mathematical model of artificial seism in one dimension to partial differential equations and discreting space partial derivative.The symplectic difference scheme is formed by using trapezoid formula for finding the numerical solution of this system, the formula corresponding to the diagonal Pade approximation of exp (x).The stability of this argorithm is proved. Some numerical.experiments of the positive problem are given and the numerical solutions are compared with the theoretical ones.AMS Classification 65M06
引文
[1]黄绪德.石油勘探对计算数学的新需求.数值计算与计算机应用,1993,14(1):71—81.[2]BubeKP,BurrigeR.Theone-dimensioninverseProblemofreflectionseismology.SIAMReview,1983,25(4):497-559.[3]FengKang.Ondifferenceschemesandsymplecticgeometry.InEd.FengKang,Procedingofthe1984BeijinSymposiumonDifferentealGeometryandDifferentialequation-ComputationofpartialequationsBeigin:SciencePress,1985,42—58.[4]汪家.分析力学.北京:高等教育出版社,1983,106—110.[5]FengKang,WuHuamo,QinMengZhao.SymplecticdiffereceschemesforlinearHamiltoncanolicalsystems.JournalofCompuationMathematics,1990,8(4):371—380.[6]LancasterP,TismenerskyM.TheTheoryofMatrices,ScondEdition,NewYork:1985,377—380.