上采样/上变频/半带滤波器设计仿真代码python

# -*- coding: utf-8 -*-
“““
Created on Sun Aug 26 23:42:40 2018

@author: Administrator
“““
import numpy as np
import matplotlib.pyplot as plt
import scipy.signal as signal

def halfBand（cutofffsripple1ripple2）:
    “““
    cutoff: 通带截止频率;
    fs: 采样频率;
    ripple1:通带最大波纹（+dB）;
    ripple2:阻带最小衰减（+dB）。
    “““
    if ripple1        ripple=ripple2
    else:
        ripple=ripple1#取最严格的限制
    #wc=2*np.pi*cutoff/fs
    f_a=0.5*fs-cutoff#阻带截止频率，半带滤波器特点，f_c+f_a=pi
    numtaps