attenuation.h

Author

Zhu Dengda (zhudengda@mail.iggcas.ac.cn)

Date

2024-07-24

Functions

cplx_t grt_attenuation_law(real_t Qinv, cplx_t omega)

品质因子Q 对 波速的影响, Futterman causal Q, 参考Aki&Richards, 2009, Chapter 5.5.1

\[ c(\omega) = c(2\pi)\times (1 + \frac{1}{\pi Q} \log(\frac{\omega}{2\pi}) + \frac{i}{2Q}) \]
其中虚数部分的正负号和书中不同,是因为书中使用的傅里叶变换的e指数符号和我们通常使用的相反。

参数:
  • Qinv[in] 1/Q

  • omega[in] 复数频率 \( \tilde{\omega} =\omega - i\zeta \)

返回:

atncoef 系数因子,作用在 \( k=\omega / c(\omega)\)的计算

void grt_py_attenuation_law(real_t Qinv, real_t omg[2], real_t atte[2])

attenuation_law函数在python中被调用的版本,长度2的数组分别表示复数的实部和虚部