FIR数字带通滤波器(Python实现)_python 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 
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_39739342/article/details/84948308

智能推荐

Java Web开发_异步处理以及前端中Vue框架的简单使用（Day3）_后端异步前端怎么处理-程序员宅基地

文章浏览阅读304次。此博客主要记录在学习黑马程序员2023版JavaWeb开发课程的一些笔记，方便总结以及复习。_后端异步前端怎么处理

python数据类型01_python[80., 20., 1000, 200]-程序员宅基地

文章浏览阅读1w次。文章目录数值类型整型(int)long（长整型）浮点数复数不同进制表示数值类型转换数据类型信息获取math 模块、cmath 模块python数学函数abs(x)ceil()cmp()exp()fabs()floor()log()log10()max()min()modf()pow()round()sqrt()python随机数函数choice()randrange()random()seed()..._python[80., 20., 1000, 200]

机器视觉halcon轮廓线处理关键算子-常州龙熙机器视觉培训班_halcon中的轮廓线 导数-程序员宅基地

文章浏览阅读876次，点赞23次，收藏21次。halcon 轮廓线处理 关键算子_halcon中的轮廓线 导数

自动驾驶人机交互HMI产品技术方案_自动驾驶hmi用什么开发-程序员宅基地

文章浏览阅读544次。HMI产品是L4车辆的人机交互程序，为高速运营、港口单车、测试路测等提供状态可视化、任务交互、自动驾驶行车控制、编队控制功能。_自动驾驶hmi用什么开发

Matlab画散点图并拟合（使用cftool函数拟合）_matlab散点图拟合函数-程序员宅基地

文章浏览阅读4w次，点赞13次，收藏120次。Matlab根据坐标点进行绘制散点图并拟合成图像可以使用cftool函数，下面以二维数据拟合进行举例：（1）首先输入数据点x=[0.20,2,4.01,5.99,8.08,9.98,11.96,14.00,15.99,18.00,19.98,21.98,23.99,25.97,28.01,30.00,32.04,33.99,35.98,37.99,39.99,42.00,43.99,45...._matlab散点图拟合函数

Java 命令行运行参数大全_命令行运行java参数-程序员宅基地

文章浏览阅读6.8k次。javac 用法：javac 其中，可能的选项包括： -g 生成所有调试信息 -g:none 不生成任何调试信息 -g:{lines,vars,source} _命令行运行java参数

随便推点

阿里云mysql空间不足_阿里云MySQL 实例空间使用率过高的原因和解决方法-程序员宅基地

文章浏览阅读419次。用户在使用 MySQL 实例时，会遇到空间使用告警甚至超过实例限额被锁定的情况。在 RDS 控制台的实例基本信息中，即会出现如下信息：本文将介绍造成空间使用率过高的常见原因及其相应的解决方法。对于MySQL 5.6版本的实例，升级实例规格和存储空间后即可解锁实例，关于如何升级实例配置，请参见变更配置。•常见原因造成 MySQL 实例空间使用率过高，主要有如下四种原因：Binlog 文件占用高。数据..._阿里云m2实例数超过限制99999

JQuery信息提示框插件 jquery.toast.js 的使用-程序员宅基地

文章浏览阅读1.1w次，点赞5次，收藏13次。1.下载https://github.com/kamranahmedse/jquery-toast-plugin在线预览地址2.导入在页面中引入jquery.toast.css文件，jquery和jquery.toast.js文件。<link type="text/css" rel="stylesheet" href="css/jquery.toast.css">..._jquery.toast.js

vue2+vue3——1~35-程序员宅基地

文章浏览阅读271次。vue2+vue3

电脑远程控制软件哪个好用?(4款远程控制软件推荐)_安企神控制软件-程序员宅基地

文章浏览阅读940次，点赞12次，收藏19次。本文介绍了四款远程控制电脑的软件，这四款远程控制电脑软件操作方法都很简单，大家可以根据自己的需要选择合适的软件即可。在另一台电脑的Chrome浏览器中登录同一个谷歌账号，打开谷歌远程桌面选择要控制的电脑，再输入PIN码即可远程控制电脑。是一款好用的电脑远程控制软件，用户可以通过网络远程连接到其他计算机，轻松实现远程监控、远程技术支持。在两台电脑上都登录QQ账号，主控端电脑打开要控制的好友聊天窗口，单击右上角的更多按钮。，在管理者的电脑上安装管理端，在员工的电脑上安装员工端，安装好后会自动进行连接和上线。_安企神控制软件

80 Gbps 的 USB4 2.0 要来了！_usb4+2.0-程序员宅基地

文章浏览阅读1w次，点赞10次，收藏7次。整理 | 郑丽媛出品 | CSDN（ID：CSDNnews）两年前，英特尔在公布新一代 Thunderbolt 4（以下简称雷电 4）接口标准时曾说：“不是所有 USB4 都能和雷电 4 平起平坐。”如今看来，这句话的顺序可能要颠倒一下了：本月初，USB 推广组官宣了 USB4 v2.0，其可通过 USB Type-C 提供高达 80 Gbps（相当于 10GB/s）的数据传输速率——不仅是 U..._usb4+2.0

【jdk8 jdk17 jdk21 在线中文文档】-程序员宅基地

文章浏览阅读123次。jdk8中文文档jdk17在线文档jdk21在线文档

推荐文章

热门文章

相关标签