fft

#!/usr/bin/env python

import sys

import typedbytes

import numpy as np

input = typedbytes.PairedInput(sys.stdin)

output = typedbytes.PairedOutput(sys.stdout)

clen = 1<<22

for (key, value) in input:

    v = np.fromstring(value, np.int16)

    if len(v) >= clen:

    fft_v = np.fft.rfft(v, clen) #np.fft库中提供了一个rfft函数，方便我们对实数信号进行FFT计算

    print >>sys.stderr, 'key:', key, 'value:',fft_v         #为了调试用的，看看输出的key和value是什么

    output.write((key, fft_v)) #输出的key是文件名，value是fft序列