graphics.py 文件源码

def __lowpass_filter(self,ite_data,fs,order,highcut,zerophase,btype,ftype,rps=None):
        fe = fs/2.0
        high = highcut/fe
        if high>1:
            high=1
        if ftype == "cheby1":
            rp = rps
            z,p,k = signal.iirfilter(order,high,btype=btype,ftype=ftype,output="zpk",rp=rp)
        elif ftype == "cheby2":
            rs = rps
            z,p,k = signal.iirfilter(order,high,btype=btype,ftype=ftype,output="zpk",rs=rs)
        elif ftype == "ellip":
            rp = rps[0]
            rs = rps[1]
            z,p,k = signal.iirfilter(order,high,btype=btype,ftype=ftype,output="zpk",rp=rp,rs=rs)
        else:
            z,p,k = signal.iirfilter(order,high,btype=btype,ftype=ftype,output="zpk")
        sos = signal.zpk2sos(z,p,k)
        ite_data = signal.sosfilt(sos,ite_data)
        if zerophase:
            ite_data = signal.sosfilt(sos,ite_data[::-1])[::-1]
        return ite_data