fortran写一个傅里叶变换的程序

下面是一个简单的Fortran程序，用于计算离散傅里叶变换（DFT）：

PROGRAM dft
    IMPLICIT NONE
    
    INTEGER :: N, k, n
    REAL :: x(100), re(100), im(100)
    
    ! 读取信号长度 N 和信号 x(n)
    WRITE(*, *) '请输入信号长度 N：'
    READ(*, *) N
    WRITE(*, *) '请输入信号 x(n)：'
    DO n = 1, N
        READ(*, *) x(n)
    END DO
    
    ! 计算 DFT
    DO k = 1, N
        DO n = 1, N
            re(k) = re(k) + x(n) * COS(2 * PI() * (k-1) * (n-1) / N)
            im(k) = im(k) - x(n) * SIN(2 * PI() * (k-1) * (n-1) / N)
        END DO
    END DO
    
    ! 输出结果
    WRITE(*, *) 'DFT结果：'
    DO k = 1, N
        WRITE(*, *) 're(', k, ') = ', re(k), ', im(', k, ') = ', im(k)
    END DO

END PROGRAM dft

这个程序使用了双重循环计算DFT，其中re(k)和im(k)分别存储了第k个频率分量的实部和虚部。注意，这个程序没有使用FFT算法，因此它的时间复杂度是O(N^2)，不适用于大规模的信号处理。如果需要处理更大的信号，可以使用现成的FFT库。