def DFT(sig): #离散傅里叶变换 t=np.linspace(0, 1.0, len(sig)) #创建等间隔时间序列 f = np.arange(len(sig)//2+1, dtype=complex) for index in range(len(f)): #计算公式2.1.2,重复上述步骤④～⑥获得傅里叶变换结果 f[index]=complex(np.sum(np.cos(2*np.pi*index*t)*sig), -np....