dfttest.py 文件源码

def parse_DFT_output(FmList, threshold=0.001):
    outputs = []
    for (i, Fm) in enumerate(FmList):
        if abs(Fm) > threshold:
            frequency = i
            amplitude = abs(Fm) * 2.0
            phase_angle = int(((cmath.phase(Fm) + pi2 + pi2 / 4.0) % pi2) / pi2 * 360 + 0.5)            
            outputs.append((amplitude, frequency, phase_angle))
    return outputs