本章先分析说明伪谱法求解波动方程的算法原理，再给出弹性波场伪谱法数值模拟的改进算法；然后，通过直接引入P波波场变量和S波波场变量，给出可实现弹性波场的P波和S波分解的波动方程，并用伪谱法实现弹性波场的P波和S波分解的数值模拟；最后用伪谱法对井间地震的复杂波场作数值模拟，分析井间地震波场中各种波的传播规律。

CSDN海神之光上传的全部代码均可运行，亲测可用，尽我所能，为你服务； 1、代码压缩包内容 主函数：Forward_modeling_PS.m； 调用函数：其他m文件；无需运行 运行结果效果图； 2、代码运行版本 Matlab 2019b；若运行有误，根据提示修改；若不会，可私信博主； 3、运行操作步骤 步骤一：将所有文件放到Matlab的当前文件夹中； 步骤二：双击打开Forward_modeling_PS.m文件； 步骤三：点击运行，等程序运行完得到结果； 4、物理应用 仿真：导航、地震、电磁、电路、电能、机械、工业控制、水位控制、直流电机、平面电磁波、管道瞬变流 光学：光栅、杨氏双缝、单缝、多缝、圆孔、矩孔衍射、夫琅禾费、干涉、拉盖尔高斯、光束、光波、涡旋 定位问题：chan、taylor、RSSI、music、卡尔曼滤波UWB 气动学：弹道、气体扩散、龙格库弹道 运动学：倒立摆、泊车 天体学：卫星轨道、姿态

1、程序功能： 1.1调整地震波格式 程序主要针对太平洋地震中心的地震波文件格式，该类格式一般为5列，读取顺序为从左至右，从上至下，如下文操作方法中所示。将下载文件导入程序后，可生成一列的地震波文件，导出保存为txt、excel。 1.2生成反应谱、与设计谱进行对比 导入地震波加速度时程后，会显示各类反应谱。因建筑抗震设计规范的设计反应谱采用绝对加速度反应谱标定，所以绝对加速度反应谱与设计反应谱进行对比。 1.3生成人工波

Hach地震波原始数据

North地震波原始数据

天津波-南北向.txt

从改进BISQ 模型双相介质所对应的一阶速度―应力运动方程出发，构建2×2N 阶交错网格有限差分模拟算法，为了尽可能地减小或消除数值模拟中由人工边界引起的虚假反射，建立完全匹配层(PML)吸收边界的2×2N 阶交错网格有限差分算法。详细地讨论了PML 吸收边界条件的构建及其有限差分算法的实现。通过MATLAB编程进行波场模拟，将加入 PML 吸收边界、常规指数衰减吸收边界及未加吸收边界的 3 种数值模拟结果进行对比，论证PML 吸收边界能十分有效地吸收边界反射。

谱元法属广义有限元法，可用于地震波场模拟与震源反演，但几乎不能在单台普通电脑上用于大尺度地震波场模拟。Microwulf是一种低成本、高效率的主从式便携型计算集群系统，能有效地利用个人电脑组建并行计算环境。基于三维弹性波动方程研究了谱元法基本算法，并利用谱元法和自搭建Microwulf系统对2008年汶川5．12地震波场进行模拟。模拟结果表明：(1)搭建的Microwulf系统能成功进行大尺度地震波场的谱元法数值模拟；(2)计算区域地形能引起地震波场畸变；(3)震中附近区域峰值地震动分布模式大体与发震断裂

利用代码合成地震波，更快更好的去接触地震波及了解合成过程

raytrace_射线追踪法_地震成像_地震波_射线_射线层析