求高手帮我详细解释一下这个matlab程序，是关于地震波数值模拟的，事后另有重谢。
通过有限差分来解波动方程模拟波场，画的两个图一个是所用的子波，一个是最后时刻的波场快照。
常规地震
万有引力 的测量大致可分为地球物理学方法测量、空间测量、实验室内测量等三大类. 地球物理学方法测量G 是利用大的自然物体(如形状规则的山体、矿井和湖泊等)作为吸引质量。
该方法的主要优点是作为吸引质量的自然物体很大, 引力效应明显. 但由于吸引质量的尺度、密度及其分布等都不能精确测量, 所以实验的精度比较低. 随着航天技术的发展, 人们期望在太空开展测G 实验。
空间测量方法可以避免地面实验室中遇到的两大难题:一个是地面实验环境中的附加背景引力场作用, 另一个是地面振动噪声的干扰, 就目前的情况来看,空间测量G 的方法面临着很多新的技术难题, 仍在探索之中.实验室内测量万有引力常数G 的常用工具是精密扭秤和天平。
与地球物理学方法相比, 精密扭秤的最大优点是将待测的检验质量与吸引质量之间的万有引力相互作用置于与地球重力场方向正交的水平面内, 这样就在实验设计上极大地减少了重力及其波动的影响。
天平可以绕刀口在垂直面内上下倾斜以探垂直方向的引力作用。
常用的测量方法有: 直接倾斜法、补偿法、共振法、周期法和自由落体法等
