okada地震模型，地球物理/地质等数学模型，matlab实现代码

Okada [1985] 模型计算了弹性半空间中由于剪切和拉伸断层引起的表面变形的解析解。 该模型广泛用于模拟由构造断层（地震）或火山岩脉（岩浆侵入）等局部扰动产生的地面变形。 给定矩形断层几何形状（长度、宽度、深度、走向、倾角）和 3 分量位错幅度（前倾、滑动和张开），它计算自由表面的位移、倾斜和应变。 建议的 Matlab 脚本是冈田方程的文字转录，除了它在地理参考（东、北、上）中转置，其中断层由相对于北的走向角定义，位错参数由下式给出: rake, slip and opening (而不是 U1, U2, U3)，遵循 Aki & Richards [1980] 的定义。 所有坐标和深度都是相对于断层质心的。 拉梅常数 λ 和 μ 由泊松比 ν 代替（各向同性介质的默认值为 0.25），因为方程与其他弹性参数无关。 这些方程还针对 (x,y) 坐标和除倾角之外的所有输入参数进行了

半空间中位错型变形引起的静态同震位移的计算 描述：这个精确解通常用于构建由于水下地震引起的同震位移。 然后使用这些位移以某种方式构建海啸传播代码的初始条件。 有关更多详细信息，请参阅Okada原始文件： * 冈田，Y. (1985)。 由于半空间中的剪切和拉伸断层导致的表面变形。 公牛。 地震。 社会。 上午，75，1135-1154。 或者我免费提供的预印本： * D. Dutykh, F. Dias，移动底部产生的水波。 书中：“海啸和非线性波”，Kundu, A.（编辑），Springer Verlag 2007 http://hal.archives-ouvertes.fr/hal-00115875/ 在作者的其他出版物中可以找到更复杂的海啸生成问题方法： http://www.denys-dutykh.com/

OKada弹性位错模型的matlab源代码