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

本文链接：https://blog.csdn.net/weixin_39739342/article/details/84948308

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

A. Mishka and Contest_x_flyx的博客-程序员秘密

A. Mishka and Contesttime limit per test1 secondmemory limit per test256 megabytesinputstandard inputoutputstandard outputMishka started participating in a programming contest. Ther...

前端开发和设计必备的Chrome插件_前端设计协作插件_ejinxian的博客-程序员秘密

Web设计、开发、调试工具集 Google Chrome是最好用的几个浏览器之一，自从2011年11月份赶超Firefox之后，已成为当今互联网的主流浏览器。今天，我来分享下自己收集的一系列Chrome插件，希望能够提高大家的工作效率。Devtools Terminal—嵌在浏览器中的终端。开发调试利器！LiveReload—集成LiveReload

Nginx Proxy buffer 工作原理_weixin_34306593的博客-程序员秘密

为什么80%的码农都做不了架构师？>>> ...

递归算法时间复杂度分析_递归的时间复杂度_AcceptedLin的博客-程序员秘密

递归算法时间复杂度分析时间复杂度：一般情况下，算法中基本操作重复的次数就是问题规模n的某个函数f（n），进而分析f（n）随n的变化情况并确定T（n）的数量级。这里用‘o’来表示数量级，给出算法时间复杂度。T（n）=o（f（n））；它表示随问题规模n的增大，算法的执行时间增长率和f（n）增长率成正比，这称作算法的渐进时间复杂度...

复习JavaWeb入门与基础知识&Java面试题HashMap在JDK1.8的改进&生活【记录一个咸鱼大学生三个月的奋进生活】030_Aspiriln的博客-程序员秘密

记录一个咸鱼大学生三个月的奋进生活030JavaWeb基础知识Tomcat相关知识JavaWeb基础知识在正式开始之前我们先来了解一下软件的分类？系统软件： Windows、Unix、Linux、Mac 等驱动程序：主板驱动、显卡驱动、声卡驱动、网卡驱动等工具软件：办公软件、杀毒软件、聊天软件、下载软件等商业应用软件：教务管理系统、售票系统、银行转账系统、电子商务网站等而我们现在要做的就是应用软件程序员，开发商业应用软件商业应用软件又分为 C/S（客户端/服务器）和 B/S（浏览

程序员难道不加班就没有前途吗?_Java后端技术的博客-程序员秘密

一、有没有不加班的程序员？有，我就是。我目前工作与一家500强外企（欧美），职称为高级软件开发专家。我在这家公司（500强外企）已经服务了10年，所有加班次数合起来不超过...