Okada：由有限矩形源引起的表面变形：计算 Okada 1985 年针对断层错位引起的位移、倾斜和应变的解。-matlab开发

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