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
评论列表
文章目录