adjoint.py 文件源码

python
阅读 21 收藏 0 点赞 0 评论 0

项目:SeisFlows_tunnel 作者: DmBorisov 项目源码 文件源码
def Envelope3(syn, obs, nt, dt, eps=0.):
    # envelope lag
    # (Yuan et al 2015, eqs B-2, B-5)
    esyn = abs(_analytic(syn))
    eobs = abs(_analytic(obs))

    erat = _np.zeros(nt)
    erat[1:-1] = (esyn[2:] - esyn[0:-2])/(2.*dt)
    erat[1:-1] /= esyn[1:-1]
    erat *= misfit.Envelope3(syn, obs, nt, dt)

    wadj = -erat*syn + _hilbert(erat*_hilbert(esyn))
    return wadj
评论列表
文章目录


问题


面经


文章

微信
公众号

扫码关注公众号