基于时间抽取的基2—FFT算法和音频文件的读取

**功能实现：**读取音频文件的数据，取其中一部分数据，生成幅度谱，相位谱。

import numpy as np
import pylab as plt
import librosa
from math import *
import wave
import time

time_start=time.time()


def fft(xn):
    '''
    计算xn的傅里叶变换,要求xn的长度为2的整数次幂
    '''
    N = len(xn)
    M = log2(N)
    assert isinstance(M, float)
    reverseOrder(xn)
    for L in range(1, int(M) + 1):#第L级蝶形运算
        B = 2 ** (L - 1)#每级蝶形运算有B个旋转因子
        for J in range(B):
            W = e ** (-1j *