Hz工频信陷波器设计
50Hz陷波器设计(matlab)

源代码:%陷波器的设计%陷波器的传输函数为% B(1/z) (z-exp(j*2*pi*f0))*(z-exp(-j*2*pi*f0))%H(z) = -------- = --------------------------------------------% A(1/z) (z-a*exp(j*2*pi*f0))*(z-a*exp(-j*2*pi*f0))%其中f0为陷波器要滤除信号的频率,a为与陷波器深度相关的参数,a越大,深度越深。
%%已知信号中50Hz工频干扰,信号为x=sin(2*pi*50*n*Ts)+sin(2*pi*125*n*Ts);%要求通过陷波器滤除50Hz干扰信号%参数设置:采样率Ts=0.001s,采样长度:512点clf;clear;%设置初值f0=50;Ts=0.001;fs=1/Ts;NLen=512;n=0:NLen-1;%陷波器的设计apha=-2*cos(2*pi*f0*Ts);beta=0.96;b=[1 apha 1];a=[1 apha*beta beta^2];figure(1);freqz(b,a,NLen,fs);%陷波器特性显示x=sin(2*pi*50*n*Ts)+sin(2*pi*125*n*Ts);%原信号y=dlsim(b,a,x);%陷波器滤波处理%对信号进行频域变换。
xfft=fft(x,NLen);xfft=xfft.*conj(xfft)/NLen;y1=fft(y,NLen);y2=y1.*conj(y1)/NLen;figure(2);%滤除前后的信号对比。
subplot(2,2,1);plot(n,x);grid;xlabel('Time (s)');ylabel('Amplitude');title('Input signal');subplot(2,2,3);plot(n,y);grid;xlabel('Time (s)');ylabel('Amplitude');title('Filter output');subplot(2,2,2);plot(n*fs/NLen,xfft);axis([0 fs/2 min(xfft) max(xfft)]);grid;xlabel('Frequency (Hz)');ylabel('Magnitude (dB)');title('Input signal');subplot(2,2,4);plot(n*fs/NLen,y2);axis([0 fs/2 min(y2) max(y2)]);grid;xlabel('Frequency (Hz)');ylabel('Magnitude (dB)');title('Filter output');100200300400500600-2-112Time (s)A m p l i t u d eInput signal100200300400500600-2-112Time (s)A m p l i t u d eFilter output010020406080100120M a g n i t u d e (d B )010020406080100120M a g n i t u d e (d B )。
陷波器设计

陷波器设计陷波器是无限冲击响应(IIR)数字滤波器,该滤波器可以用以下常系数线性差分方程表示:∑∑==---=M i Ni i i i n y b i n x a n y 01)()()( (1)式中: x(n)和y(n)分别为输人和输出信号序列;i a 和i b 为滤波器系数。
对式(1)两边进行z 变换,得到数字滤波器的传递函数为:∏∏∑∑===-=---==N i i M i i N i i iM i ii p z z z z b z a z H 1100)()()( (2)式中:i z 和i p 分别为传递函数的零点和极点。
由传递函数的零点和极点可以大致绘出频率响应图。
在零点处,频率响应出现极小值;在极点处,频率响应出现极大值。
因此可以根据所需频率响应配置零点和极点,然后反向设计带陷数字滤波器。
考虑一种特殊情况,若零点i z 在第1象限单位圆上,极点i p 在单位圆内靠近零点的径向上。
为了防止滤波器系数出现复数,必须在z 平面第4象限对称位置配置相应的共轭零点*i z 、共轭极点*i p 。
这样零点、极点配置的滤波器称为单一频率陷波器,在频率ωo 处出现凹陷。
而把极点设置在零的的径向上距圆点的距离为l-μ处,陷波器的传递函数为:))1()()1(())(()(2121z z z z z z z z z H μμ------= (3) 式(3)中μ越小,极点越靠近单位圆,则频率响应曲线凹陷越深,凹陷的宽度也越窄。
当需要消除窄带干扰而不能对其他频率有衰减时,陷波器是一种去除窄带干扰的理想数字滤波器。
当要对几个频率同时进行带陷滤波时,可以按(2)式把几个单独频率的带陷滤波器(3)式串接在一起。
一个例子:设有一个输入,它由50Hz 信号和100Hz 信号组成。
50Hz 是一个干扰信号,要设计一个50 Hz 的带陷滤波器,采样频率为400Hz 。
4/400/5021ππω=⨯=因此z 平面上的零极点可设置为4/14/1999.0ππj j ep e z ±±== 展开式为70637064)707.0707.0(999.0)4sin 4(cos 999.0999.0707.0707.022224sin 4cos 4/14/1j j j e p j j j e z j j ±=±=±=±=±=±=±±ππππππ== 它的传递函数为2121221111999.04126.11414.11999.04126.11414.1)7063.07063.0)(7063.07063.0()707.0707.0)(707.0707.0())(())(()(----**+-+-=+-+-=+---+---=----=z z z z z z z z j z j z j z j z p z p z z z z z z H因此分子系数是[1 1.414 1];分母系数是[1 1.4126 0.999]。
低通50Hz陷波器

2.1电路工作原理.........................................................4
2.2元器件选择...........................................................5
学生姓名:学号:专业班级:
一、设计条件
1.可选元件
(1)选题规定的“可选、限选元件”
(2)电阻、电容、电感、电位器等,按需使用
(3)自备元件
2.可用仪器
万用表,示波器,交流毫伏表,信号发生器,直流稳压电源
二、设计任务及要求
1.设计任务
根据技术要求和已知条件,完成选题电路的设计、装配与调试。
2.设计要求
(1)选题规定的“设计内容和要求”;
(2)选择电路方案,完成对确定方案电路的设计。包括:计算电路元件参数、选择元件、画出总体电路原理图;
(3)用软件仿真整体或部分核心实验电路,得出适当结果;
(4)装配、调试作品,按规定格式写出课程设计报告书。
三、时间安排
1.第12周前:布置设计任务,讲解设计要求、实施计划、设计报告等要求。
2.1电路工作原理
陷波器也称带阻滤波器(窄带阻滤波器),它能在保证其他频率的信号不损失的情况下,有效的抑制输入信号中某一频率信息。所以当电路中需要滤除存在的某一特定频率的干扰信号时,就经常用到陷波器。本实验所用为双T带阻滤波网络。
2.2元器件选择
元器件选择清单
序号
名称
型号
数量
1
双运放
陷波器设计

数字信号处理课程设计报告书题目陷波器设计课程设计任务书课题题目摘要随着数字技术的发展,数字滤波器在许多领域得到广泛的应用。
它是通信、语言、图像、自动控制、雷达、航空航天、生物医学信号处理等领域中的一种基本处理部件,具有稳定性好、精度高、灵活性大等突出优点。
在信号采集时,往往受到50Hz电源频率干扰,尤其是在供电系统不稳定、外界环境适应性差时严重影响要采集信号的正确判断。
本设计研究一种在MATLAB语言环境下分别用IIR和FIR滤波器设计方法设计实现一个数字陷波器,并将设计的滤波器应用到混合的正弦信号,通过仿真测试,用两种方法设计的滤波器可以很好的消除50Hz的工频干扰,并分析比较了各种方法所设计的陷波器性能。
在设计IIR数字陷波器过程中,是用椭圆数字陷波滤波器的设计方法,而FIR数字陷波器的设计主要用窗函数法、频率采样法及等波纹逼近法。
FIR滤波器可以得到严格的线性相位,但它的传递函数的极点固定在原点,只能通过改变零点位置来改变性能,为了达到高的选择性,必须用较高的阶数,对于同样的滤波器设计指标,FIR滤波器要求的阶数可能比IIR 滤波器高5~10倍。
IIR滤波器的设计可以利用模拟滤波器的许多现成的设计公式、数据和表格,计算的工作量较小。
关键词数字陷波器;50Hz工频干扰;IIR和FIR滤波器目录课程设计任务书 (I)摘要 (II)1设计概述 (1)1.1设计背景 (1)1.2设计目的 (1)1.3设计内容及要求 (1)2设计方案及实现 (2)2.1总体方案设计 (2)2.2设计原理 (2)2.2.1数字陷波器原理 (2)2.2.2IIR数字陷波器原理 (3)2.2.3FIR数字陷波器原理 (3)3设计结果分析 (8)3.1IIR数字陷波器设计 (8)3.2FIR数字陷波器设计 (10)3.2.1用窗函数法设计陷波器 (10)3.2.2频率采样法设计陷波器 (12)3.2.3基于切比雪夫等波纹逼近法 (13)4总结 (16)1 设计概述1.1 设计背景在我国采用的是50Hz 频率的交流电,所以在平时需要对信号进行采集处理和分析时,常会存在50Hz 的工频干扰,对我们的信号处理造成很大干扰,因此50Hz 陷波器在日常成产生活中被广泛应用,其技术已基本成熟。
基于UAF42通用滤波芯片的50Hz陷波器设计

第27卷 第8期2006年8月仪器仪表学报Chinese Journal of Scientific InstrumentVol127No18Aug12006基于UAF42通用滤波芯片的50H z陷波器设计3颜 良1 陈儒军1,2 刘 石1 陈一平11(中南大学信息物理工程学院 长沙 410083)2(中南大学信息科学与工程学院 长沙 410083)摘要 为了压制工业电源的干扰(如50Hz干扰),经常要在输入电路中采用陷波电路。
由于传统的双T陷波电路对元件的精度和对称性要求比较高,在实际中很难调试。
文中介绍的由通用滤波器芯片UA F42构成的50Hz陷波电路克服了这方面的缺点。
实验结果表明:基于UA F42设计的50Hz陷波器的频率响应达到了设计要求,具有调试方便、对元件精度要求不高的优点。
关键词 陷波器 UAF42 FIL TER42中图分类号 TP337 文献标识码 A 国家标准学科分类代码 510110Design of a50H z notch f ilter using the universal active f ilter chip UAF42Yan Liang1 Chen Rujun1,2 Liu Shi1 Chen Yiping11(School of I nf o2p hysics and Geomatics Engineering,Cent ral S out h Universit y,Changsha410083,China)2(College of I nf ormation Science and Engineering,Cent ral S outh Universit y,Changsha410083,China)Abstract To supp ress t he50Hz interference,a notch filter is needed in pre2stage circuit.The conventional twin2T network requires high2p recision component s and st rict symmet ry,but it is difficult to t une.The50Hz notch filter using t he universal active filter chip(UA F42)eliminates above disadvantages.The experiment re2 sult s show t hat t he50Hz notch filter using UA F42meet s t he design requirement s.The design has t he advanta2 ges of easy t uning and lower component tolerance requirement.K ey w ords notch filter UA F42 FIL TER421 引 言 物探仪器工作频率低,且一般在野外工作,易于受外界信号干扰,特别是在此频率段内受到强烈的干扰源,如50Hz或60Hz市电电网信号影响。
陷波器设计

陷波器设计由传递函数的零点和极点可以大致绘出频率响应图。
在零点处,频率响应出 现极小值;在极点处,频率响应出现极大值。
因此可以根据所需频率响应配置零 点和极点,然后反向设计带陷数字滤波器。
考虑一种特殊情况,若零点 Z |在第1 象限单位圆上,极点P i 在单位圆内靠近零点的径向上。
为了防止滤波器系数出 现复数,必须在z 平面第4象限对称位置配置相应的共轭零点Z j 、共轭极点p i < 这样零点、极点配置的滤波器称为单一频率陷波器, 在频率①o 处出现凹陷 而把极点设置在零的的径向上距圆点的距离为1-卩处,陷波器的传递函数为:(Z Z 1)(Z Z 2) (z (1 )zj(z (1 )Z 2)式⑶ 中卩越小,极点越靠近单位圆,贝擞率响应曲线凹陷越深,凹陷的宽 度也越窄。
当需要消除窄带干扰而不能对其他频率有衰减时,陷波器是一种去除 窄带干扰的理想数字滤波器。
当要对几个频率同时进行带陷滤波时,可以按 (2)式把几个单独频率的带陷 滤波器(3)式串接在一起。
一个例子:设有一个输入,它由50Hz 信号和100Hz 信号组成。
50Hz 是一个干扰 信号,要设计一个50 Hz 的带陷滤波器,采样频率为400Hz=12 50/400/4因此z 平面上的零极点可设置为/4P 1 0.999ej /4=0.999(COS : j sin :) 0.999(0.707 j0.707) 7064 j7063它的传递函数为陷波器是无限冲击响应 差分方程表示: My(n) aX n式中:x(n)和y(n)分别为输人和输出信号序列;a :和b 为滤波器系数。
对式(1)两边进行z 变换,得到数字滤波器的传递函数为:Ma i z 1H(z) + b i z 1I 0式中:Z |和P i 分别为传递函数的零点和极点。
i)(IIR)数字滤波器,该滤波器可以用以下常系数线性by( n I) (1)I 1M(z Z |)■N ----------⑵(z P i ) I 1H(z)P 1 展开式为乙0.999e/4e j /4= cos_ 42 2j sin 40.707 j0.707(z Z i)(z Z1)(Z P i)(z p 1) (z 0.707 j0.707)(z 0.707 j 0.707)(z 0.7063 j0.7063)(z 0.7063 j0.7063)2z 1.414z 1 1 1.414z2z 1.4126z 0.999 1 1.4126z 0.999z因此分子系数是[1 1.414 1] 差分方程有y(n) a(2)y( n 1) a(3)y(n y(n) x(n) b(2)x( n 1) ;分母系数是[1 1.4126 0.999]。
自编IIR型50hz陷波器,Matlab源代码,有详细注释

Matlab源代码:%% 巴特沃斯型50hz陷波器,可改变陷波带宽和阶数function [Num,Den] = ZB_50_filter(f0,B1,N)%输入参数% f0 陷波器中心频率% B1 单边带宽% N 滤波器阶数fs=1000;%采样率T=1/fs;rp=3;%通带衰减rs=N/2*10; %阻带衰减wp1=((f0-B1)/1000)*2*pi;% 下通带截止频率wp2=((f0+B1)/1000)*2*pi;% 上通带截止频率ws1=((f0-B1/4)/1000)*2*pi;%阻带下限频率ws2=((f0+B1/4)/1000)*2*pi;%阻带上限频率wc1=(2/T)*tan(wp1/2);wc2=(2/T)*tan(wp2/2);wr1=(2/T)*tan(ws1/2);wr2=(2/T)*tan(ws2/2);w0=sqrt(wc2*wc1);%阻带中心频率B=wc2-wc1;%带宽wp=1;ws=wp*(wr1*B)/(w0^2-wr1^2);%归一化阻带截止频率[N,wc]=buttord(wp,ws,rp,rs,'s');[Z,P,K]=buttap(N);[Bd,Ad]=zp2tf(Z,P,K);%将零极点转化成传输形式[B,A]=lp2bs(Bd,Ad,w0,B);%对低通滤波器进行频率转换,为带阻滤波器[b,a]=bilinear(B,A,fs);%利用双线性变换法,转换成数字滤波器Num=b;%滤波器分子Den=a;%滤波器分母end举例:设计中心频率为50hz,带宽为0.6hz的IIR的陷波器:f0=50;B1=0.3; %单边带宽0.3hzN=2;[Num,Den] = ZB_50_filter(f0,B1,N);%调用函数计算结果:Num=[0.998118588556613 -1.898537748170150.998118588556613];Den=[1 -1.89853774817015 0.996237177113225]调用matlab中的fdatool滤波器工具箱验证:与fdatool自带IIR的butterworth滤波器的幅频曲线一致,滤波器系数几乎一致:。
50Hz工频陷波电路图

50Hz工频陷波电路图图2前置放大器电路A VCC LMC6464主放大电路<无线心电信号采集系统研究>李罗,卢建刚计算机测量与控制英文刊名:COMPUTER MEASUREMENT & CONTROL 年,卷(期): 2006,14(12)被引用次数: 2次图2 系统原理图图 1 nRF905 模块的高频头用户接口电路管脚[21] 李朝青.无线发送/接收IC 芯片及其数据通信技术选编.北京:北京航空航天大学出版社.2003[26] 张毅刚,彭喜元.单片机原理与应用设计.哈尔滨工业大学出版社.2008,7 NRF905 应用原理图心电信号接收电路图通信电路图3-4 单片机最小系统板显示模块无线射频收发芯片nRF905内置有天线,同时内部集成有调制,解调、编码/解码等功能,故在通信过程中能自动生成前导码和CRC校验,而不需要"接人"网络就能享受通信服务。
本设计根据nRF905的特点设计的无线数据收发系统,经过多次实验证明,其发射端能正确地将数据传送出去;同时,经nRF905发射后,接收端也能正确接收并显示数据,有效通信距离大于200米。
在有障碍物体的混凝土结构的建筑内测试,其有效直线通信距离大于50 m。
此外,该系统采用了比较完善的软件、硬件设计以及抗干扰措施,这样,就可以保证系统工作的安全性和可靠性,并具有通用性,便于投入实际应用,而且稍作改动就可以应用到小区传呼、工业数据采集、生物信号采集,无线遥控等其它一些短距离无线通信领域,以实现无线数据的双向传输,具有较好的市场应用价值。
为了减少电路板对无线射频部分的通信干扰,本系统中所用的nRF905将构成一个最小系统作为一个功能模块,直接以功能模块的形式引出与单片机的接口。
这样既敬爱年少了电路板对无线通信的干扰,又方便了无线模块的更换。
nRF905芯片内部集成了电源管理、晶振、低噪放大器、频率合成器、功放等模块,可以自动完成处理字头和CRC(循环冗余码校验)的工作,片内硬件自动完成曼彻斯特编码/解码。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第1章摘要本文介绍一种基于运算放大器的工频信号陷波器的设计与制作,用以消除叠加在频率为1kHz以上的测试信号中所包含的50Hz工频信号。
叙述内容包括工频信号陷波器的工作原理与设计思路,介绍了陷波器的参数计算及其选择,通过multisim 仿真,记录和分析了该陷波器的工作特性与陷波性能,论证了该陷波器的可行性。
此次设计的陷波器优点是:陷波性能良好,带宽较小,品质因数Q可调,即滤波性能便于调整,电路线路简单,具有实际应用价值。
缺点是:对于元器件的参数要求高,需要仔细调节。
第2章设计原理概述及设计要求
陷波器的基本原理及作用
陷波器也称带阻滤波器(窄带阻滤波器),它能在保证其他频率的信号不损失的情况下,有效的抑制输入信号中某一频率信息。
所以当电路中需要滤除存在的某一特定频率的干扰信号时,就经常用到陷波器。
在我国采用的是50hz频率的交流电,所以在平时需要对信号进行采集处理和分析时,常会存在50hz的工频干扰,对我们的信号处理造成很大干扰,因此50Hz陷波器在日常成产生活中被广泛应用,其技术已基本成熟。
工频陷波器不仅在通信领域里被大量应用,还在自动控制、雷达、声纳、人造卫星、仪器仪表测量及计算机技术等领域有着广泛的应用。
设计要求
1:完成题目的理论设计模型;
2:完成电路的multisim仿真;
3:完成一份设计说明书(其中包括理论设计的相关参数以及仿真结果);
4:提交一份电路原理图
第3章基于运算放大器的工频信号陷波器设计
理论分析
f和抑制带宽BW之间的关系为:
陷波器的实现方法有很多,本次设计采用的是电路比较简单,易于实现的双T 型陷波器。
双T型带阻滤波器的主体包括三部分内容:选频部分、放大器部分、反馈部分。
此陷波器具有良好的选频特性和比较高的Q值。
图双T 型陷波器电路
图中,2A 用作放大器,其输出端作为整个电路的输出。
1A 接成电压跟随器的形式。
因为双T 网络只有在离中心频率较远时才能达到较好的衰减特性,因此滤波器的Q 值不高。
加入电压跟随器是为了提高Q 值,此电路中,Q 值可以提高到50以上,调节1R 、2R 两个电阻的阻值,来控制陷波器的滤波特性,包括带阻滤波的频带宽度和Q 值的高低。
在图2中,O C U U =, 1C Z sC
=
, 2212O O R U U R R =+,令212R K R R =+,1
n R = 对节点A 列KCL 方程,得:
()()()2i A O A A O U U sC U U sC n U KU -+-=- (1)
同样,对节点B 列KCL 方程,得:
()()()2i B O B B O U U n U U n sC U KU -+-=- (2)
同样,对节点C 列KCL 方程,得:
()()A O O B U U sC U U n -=- (3)
由式(1)、(2)、(3)可得到电路的传递函数为:
()222
22
2
s 44o i U s C n H U s C n nsC nsCK
+==++- (4) 令s j ω=,
01RC
ω=
得:
()()2
02
001141H j j K ωωωωωωω⎛⎫- ⎪⎝⎭=⎛⎫⎛⎫-+- ⎪ ⎪⎝⎭⎝
⎭ (5)
由带阻滤波器的标准形式:()2
020011H j j Q ωωωωωωω⎛⎫- ⎪
⎝⎭=
⎛⎫ ⎪
⎛⎫⎝⎭-+ ⎪⎝⎭
(6)
可得:()
1
41Q K =
- ,其中212R K R R =+
当0ωω=时,()1H j ω≈
当ω远大于0ω,或者ω远小于0ω时,增益接近1。
令()0.707H ω=,可以求得:
()021H f f K ⎤=-⎥⎦
(7)
()021L f f K ⎤=-⎥⎦
(8)
()
01
41H L f Q f f K =
=-- (9)
()041H L BW f f K f =-=- (10)
因此,当 1K ≈时,Q 值极高,BW 接近于零,所以我们可以改变K 的值来调节带宽,Q 值越大,带阻曲线越窄,限波效果约好,但是实际应用中Q 值不能无限大,如果Q 值过大,会引起电路的振荡,不稳定。
图 双T 陷波器的幅频特性
其中心频率0f 的计算公式为:
01
2f RC
π=
(11) 因此T 型结构中的电容和电阻用来确定中心频率的值,可 以通过改变这些电容和电阻来选择需要滤除的频率值。
参数的具体计算及选择
本次设计要求消除叠加在频率为1kHz 以上的测试信号中所包含的50Hz 工频信
号,所以中心频率050z f H =,1
314
RC =,令R= Ω,则C=100nf 。
带宽BW 越小越好,取51H f Hz ≥,49L f Hz ≤,可得()0412BW K f =-≤,25Q ≥ 即2
12
0.99R K R R =
>+,2199R R >,取Ω=5001R Ω=K R 1002
为了防止中心频率漂移,要使用镀银云母电容或碳酸盐电 容和金属膜电阻。
常见衰减量为40—50dB ,如果要得到60dB 的 衰减量,必须要求电阻的误差小于%,电容误差小于%。
电路组成
双T 陷波器Multisim 仿真电路。
图双T 陷波器Multisim 仿真电路
仿真图中的元器件参数:
第4章基于运算放大器的工频信号陷波器性能测试
multisim 50Hz正弦频率信号源仿真
图f=50Hz时实测波形
图
图f=49Hz时实测波形图
图f=51Hz时实测波形图
仿真数据记录及简要分析
表f=49Hz、50Hz、51Hz时数据记录
频率/Hz 输入幅值/V 输出幅值/V 衰减/dB 相位差
49 π
50 π
分析:当频率为49Hz 的信号通过陷波器时,由于该频率小于第一窄带阻带的
49.505L f Hz =,因此衰减较少,所以该信号能通过陷波器。
当频率为50Hz 的信号通过
陷波器时,由于该频率等于第一窄带阻带的中心频率,所以该信号不能通过陷波器。
当频率为51Hz 的信号通过陷波器时,由于该频率大于第一窄带阻带的50.500H f Hz =,因此衰减较少,所以该信号能通过陷波器。
注:本小节multisim 仿真图中黄色(颜色较浅)是输入曲线,红色(颜色较深)是输出曲线。