python实现 快速傅里叶变换(fft)、滤波器设计(butter/chebv)

频谱分析

fft的意义在于频谱分析，基于python实现的fft、filter代码如下。其中滤波器分别用了巴特沃斯和契比雪夫二个设计，参数见源码的解释。比较复杂的是契比雪夫的第二个参数的意义。官方原文的解释为The maximum ripple allowed below unity gain in the passband. Specified in decibels, as a positive number.

---

from scipy import signal
import numpy as np
import math
import matplotlib.pyplot as pl

def fft(yt, sampling_rate, fft_size=None):  
    if fft_size is None:  
        fft_size = len(yt)  
    yt = yt[:fft_si