带通滤波器
数字信号综合设计题目简述:
（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可能有些人会感觉很陌生，今天小编就专题讲解关于电脑安全设置差异还有系统的不同，会出现不同的安装方法，还有
下面本文将从探源头、用第一性原理思考（落定义）、看数字、理性推演这四个角度，来聊聊事件分析的正确姿势。这两天朋友圈被腐国星巴克
最近在复习动态规划问题，在处理挖金矿问题的时候发现网上以python实现的代码很少，于是自己整理一份。 问题描述：漫画图解 公式和讲解
1.1 MATLAB图像处理基本操作 本文中对于大多数的操作，是对数字图像处理领域中最为著名的“lena”图片进行操作的。原图如下（Figur
# encoding: utf-8 ''' #!/usr/bin/env python @author: yudian @contact: [email protected] @file: k_means.py @time: 2018/12