FIR数字带通滤波器(Python实现)

数字信号综合设计题目简述:

（1）基带信号的产生及调制

def py_init(jidai_1,sin1_hz,zxb1_hz):    #（常数，载波频率，基带信号频率）
	y_sin1=10*np.sin(2*sin1_hz*np.pi*x)  #sinw1产生
	a=0.5*np.sin(2*zxb1_hz*np.pi*x)	#基带方波产生
	y_zxb1=[]
	for i in a:
		if i>0:
			y_zxb1.append(0.5+jidai_1)
		else:
			y_zxb1.append(-0.5+jidai_1)
	x_1=np.multiply(y_sin1,y_zxb1)		#相乘
	return x_1

（2）带通滤波器

def fda(x_1,Fstop1,Fstop2): #（输入的信号，截止频率下限，截止频率上限）
	b, a = signal.butter(8, [2.0*Fstop1/fs,2.0*Fstop2/fs], 'bandpass')
	filtedData = signal.filtfilt(b,a,x_1)
	return filtedData

（3）低通滤波器

def fda2(x_1,Fc):  #（输入的信号，限制频率）
	b, a = signal.butter(8,2.0*Fc/1000, 'lowpass')  
	filtedData = signal.filtfilt(b, a,x_1) 
	return filtedData

（4）设计系统各参数

采样频率为1000hz

fs=1000
ysc=0.2
x=np.arange(0,ysc,ysc/(ysc*fs))

基带信号1的频率为10hz，峰峰值为1V，sinw1的频率为100hz，峰峰值为20V

py_init(1.1,100,10)

基带信号2的频率为10hz，峰峰值为1V，sinw2的频率为300hz，峰峰值为20V

py_init(1.1,300,10)

带通滤波器1的截止频率下限为50hz，截止频率上限位250hz

fda(y_add,50,250)

带通滤波器2的截止频率下限为200hz，截止频率上限位490hz

fda(y_add,200,490)

低通滤波器1的截止频率为100hz

fda2(y_mul_1,100)

低通滤波器2的截止频率为300hz

fda2(y_mul_1,300)

（5）结果

基带信号与100HZ载波调制之后的频谱图

基带信号与300HZ载波调制之后的频谱图 

混合调制波通过带通滤波器1后的频谱图

 

混合调制波通过带通滤波器2后的频谱图

通过低通滤波器1后的频谱图

通过低通滤波器2后的频谱图

 300Hz频谱分析(由上到下:未通过通带滤波器,通过通带滤波器后与sinw2相乘,通过低通滤波器)

 

完整python代码:https://download.csdn.net/download/weixin_39739342/10841741