复平面分解T-S波和声波的FORTRAN代码

复平面分解可以用来将一个信号分解成不同频率的分量。对于地震信号，根据其传播方式的不同，通常可以将其分解为P波（纵波）、S波（横波）和T波（面波）等不同的分量。下面是一份将地震信号分解成T波和声波信号的FORTRAN代码：

! T-S波分解
subroutine ts_decomposition(x, t, n, t_wave, s_wave, dt)
implicit none
real, intent(in) :: x(n)       ! 原始地震信号
real, intent(in) :: t(n)       ! 时间轴
integer, intent(in) :: n       ! 数据点数
real, intent(out) :: t_wave(n) ! T波信号
real, intent(out) :: s_wave(n) ! 声波信号
real, intent(in) :: dt         ! 采样时间间隔
real :: omega, k, c, a, b, theta, r
integer :: i
do i = 1, n
    omega = 2.0 * pi / (n*dt)  ! 角频率
    k = omega / 2.0            ! 波数
    c = 2.0 * sqrt(2.0) / (3.0*dt)  ! T波速度
    a = x(i)
    b = (x(i+1) - x(i-1)) / (2.0*dt)
    theta = atan(b, a)
    r = sqrt(a**2 + b**2)
    if (theta >= -pi/2.0 .and. theta <= pi/2.0) then
        t_wave(i) = r*exp(cmplx(0.0, theta))
        s_wave(i) = 0
    else
        t_wave(i) = 0
        s_wave(i) = r*exp(cmplx(0.0, theta))
    end if
end do
end subroutine ts_decomposition

其中，`x` 表示原始地震信号，`t` 表示时间轴，`n` 表示数据点数，`t_wave` 表示T波信号，`s_wave` 表示声波信号，`dt` 表示采样时间间隔。该代码中将地震信号按照时间点逐一处理，计算每个时间点对应的T波和声波信号的复数值，再根据复数值的幅角将其分类。其中 $\omega$ 表示角频率，$k$ 表示波数，$c$ 表示T波速度，$a$ 和 $b$ 分别表示原始信号在当前时间点和相邻两个时间点之间的差值。

需要注意的是，这份代码只是一种基本的实现思路，具体实现还需要考虑数据的特点和算法的效率等问题。