必威体育Betway必威体育官网
当前位置:首页 > IT技术

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

时间:2019-07-16 21:41:03来源:IT技术作者:seo实验室小编阅读:73次「手机版」
 

带通滤波器

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

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 

相关阅读

支付宝数字证书mac安装控件方法和操作流程?

关于支付宝数字证书mac可能有些人会感觉很陌生,今天小编就专题讲解关于电脑安全设置差异还有系统的不同,会出现不同的安装方法,还有

事件分析的正确姿势:探源头、落定义、看数字、用理性(以

下面本文将从探源头、用第一性原理思考(落定义)、看数字、理性推演这四个角度,来聊聊事件分析的正确姿势。这两天朋友圈被腐国星巴克

动态规划挖金矿问题——python实现

最近在复习动态规划问题,在处理挖金矿问题的时候发现网上以python实现的代码很少,于是自己整理一份。 问题描述:漫画图解 公式和讲解

matlab数字图像处理(一) 基本操作

1.1  MATLAB图像处理基本操作 本文中对于大多数的操作,是对数字图像处理领域中最为著名的“lena”图片进行操作的。原图如下(Figur

python实现K_mean算法

# encoding: utf-8 ''' #!/usr/bin/env python @author: yudian @contact: [email protected] @file: k_means.py @time: 2018/12

分享到:

栏目导航

推荐阅读

热门阅读