基于AVR单片机的数字滤波器滤除工频干扰的快速算法

发布者:RadiantRiver最新更新时间:2013-12-09 来源: eefocus关键字:AVR单片机  数字滤波器  工频干扰 手机看文章 扫描二维码
随时随地手机看文章

工频干扰广泛存在各种工业现场中,其产生的途径主要包括输电馈线、照明设备、发动机以及各种电子仪器设备等。一般可以通过滤波电路消除工频干扰,但这必将增加硬件结构的复杂程度。实际上,还可以采用数字信号处理的相关算法,通过软件滤波器滤除工频干扰。软件滤波算法的采用,无疑会在简化电路结构的同时,使系统的硬件资源得到更加充分的利用,并达到降低产品成本的要求。

AVR 单片机是Atmel 公司生产的8位精简指令集(RISC)单片机。与同类单片机相比,在运算速度、外设资源、灵活性等方面性能均衡,性价比较高。 AVR 单片机适合C语言开发,Mega 系列AVR 单片机还有一个内部硬件乘法器单元。这些特点都为软件滤波器的实现提供了极大的便利。

1 滤波算法

常用的单片机滤波算法包括中值滤波、均值滤波、滑动平均以及复合滤波算法等。工频干扰的频率范围在50Hz附近,可以采用一个截止频率远低于50Hz的低通滤波器来滤除工频干扰。

假设输入信号为x(t),输出信号为y(t),则一个RC低通滤波器表达式为:,其中t为时间常数。

连续时间信号经过采样后成为离散时间信号,低通滤波器的表达式也变为:

 

为采样时间间隔,k为归一化时间。

由(1)式可得:

 

所以系统的传递函数为:

假设采样频率Fs=500Hz,a分别取0.8, 0.85, 0.9, 0.95, 代入(4)式,利用matlab 画出频率响应曲线,如图1。其中50Hz频率对应的幅度衰减见表1。

[page]

2 定点小数表达方式

为了精确构造数字滤波器,经常要用到浮点数据和系数。在进行浮点数乘法运算时,针对AVR 单片机设计的C编译器例如AVR-GCC,需要加入额外的数学库函数进行编译。而这会使编译后程序的代码量增加、处理时间加长、处理器的开销也随之增加。为了更大限度地降低系统开销,提高程序效率,采用定点小数表示形式进行乘法运算是最佳选择。

AVR 单片机是一种8位精简指令集(RISC)单片机。其中megaAVR 系列内部都带有一个硬件乘法器,计算一次8位乘8位的定点乘法只需2个时钟周期。因此采用8位定点采样数据乘以8位系数的定点乘法方式完成滤波算法是最高效的。

低通滤波处理涉及的运算形式为一个纯小数系数和一个已知数据相乘再相加。因此将系数采用定点小数的表示形式,对于提高算法速度是至关重要的,可以定义一种8 位定点小数表示形式--Q8 数,其各位权系数如下:

Q8 数的表示范围从0到1-2-8=0.99609375,每两个数之间的间隔是2-8(0.00390625),其所能表达的纯小数共有28=256个。例如11011000就表示2-1+2-2+2-4+2-5=0.84375,而11011001 就是表示2-1+2-2+2-4+2-5+2-8=0.84765625 ,因此0.84375 和0.84765625 之间的纯小数只能用这两个数中的一个近似表示了。这对于乘法计算的精度有一定的影响,但是由于滤波公式(3)中的系数a和(1-a)都是常数,在整体性能稳定的情况下,系数微小的不确定性对滤波器整体性能并没有太大的影响。

3 分配系数法原理

从(3)式可知,滤波算法可以用迭代计算实现,为保证每个新的输出值都可以作为下次计算的输入值,必须使输出值和输入值的位宽度一致。 AVR单片机内部硬件乘法器的输出结果为16 位,两次乘法运算的结果还要进行加法运算,其结果很有可能超过16 位宽度。如果要进行迭代计算,就要将乘加运算的结果转化成8 位表示方式。一种解决方法是用查表法实现乘法计算,这样运算结果就直接表达成8 位定点数形式,不用进行表示方式的转化,但是这种方案要占用额外的硬件存储空间构造一张查找表。

可以从逆向进行思考:由(3)式可知,每个新的输出值y(k)都与上一次的输出值y(k-1)和新的输入值x(k)有关。 y(k-1)和x(k)都是8 位的,因此最大值为0xFF。为了使a×y(k-1)+(1-a)×x(k)不超过0xFFFF,两个系数a和(1-a)的和不能超过0xFFFF/0xFF=0x101 。实际上,a+(1-a)等于“1” ,因此这里的0x101 就可以看作“1” 。如果取a=0.9,那么对应地将0x101 平均分成10份,取其中的9 份,即0x101×0.9 近似等于0xE7 ,相应地0.1 就等于0x101-0xE7=0x1A。这里的0xE7 可以近似被认为是0.9 的一种定点Q8 数表示形式,而0.1 的定点Q8 数表示形式就是0x1A。由于滤波器系数a和(1-a)采用了Q8 数的表示形式,这种将16 位乘加运算结果转化为8 位定点数表示形式的工作就变得简单了,只需通过移位运算,取y(k)的高8 位即可,对应的C语言代码为:

y(k)=(char)(y(k)>>8)

在C语言编程处理中,并不需要建立一个数组来储存y(k)的值,而只需定义两个unsigned char型的变量分别储存y(k-1)和x(k) 。当乘加计算a×y(k-1)+(1-a)×x(k),将结果转化成8 位定点数形式,再将其赋值给y(k-1)所对应的变量即可。因此采用迭代方式进行乘加运算后,整个运算过程只需要两个变量和两个常数参加即可。

通过这种处理,y(k)就可以作为计算下一次输出值y(k+1)的一个已知量,并继续与Q8 数形式的滤波器系数相乘,得到新的输出值。这种处理方式简化了乘加运算的完成过程,节省了系统硬件资源,并降低了处理器开销。

4 采样时间的控制

采用单片机进行数字信号处理,一种有效而准确的数据采集方式就是通过计数器中断服务程序(ISR)控制AD对输入信号进行精确采样。但是(图2)中断服务程序(ISR)的开销影响了AD采样时间间隔的精确度,同时如果中断服务程序(ISR)的开销过大,必然导致AD的最高采样频率的降低。因此,要想获得精确的采样频率,就必须在尽量减少中断服务程序开销的前提下,适当调整计数器中断的时间间隔。这可以通过调整OCRO的预置数来完成。

5 算法流程图

滤波算法是通过中断服务程序(ISR)来完成的,整个应用程序的主函数main()主要负责初始化计数器中断,并处理其它应用。整个程序的流程图如图3 所示。

[page]

本算法的C语言代码(附录A)经过AVR-GCC编译器的编译后,“.text”段只有310个字节,大大节省了单片机的flash 空间。

6 基于VMLAB的滤波系统仿真实现

VMLAB的全称为:Visual Micro lab。它针对AVR系列单片机和ST62 系列单片机设计,是一个单片机的虚拟原型框架,可以提供给用户一个真正意义上的虚拟微控制器(MCU)设计实验室。它具有强大的多窗口、多文件的编辑器,微控制器的集成开发环境,拥有一系列的集成开发工具,图形界面的调试器,混合模式的模拟-数字电路仿真器,代码质量检查器等。基于MCU,它可以仿真出包括模拟元器件在内的更多外围设备,并具有交互式器件模拟仿真功能。

假设有用信号是2V大小的直流信号,工频干扰是峰峰值为1V,频率为50Hz的正弦波,建立单片机AD的输入信号表示形式如下:

VMLAB 通过工程文件来管理和控制各种仿真信息、硬件连接以及显示I/O 电压波形等。根据本算法的特点,采用Atmega16作为目标单片机,时钟选为8MHz ,建立工程文件。恰当设置OCRO寄存器,使计数器比较匹配中断的时间间隔约为2ms ,这样AD的采样频率近似认为Fs等于500Hz 。经过仿真,对比结果如表3 。

从表3可以看出:随着a的增大,算法收敛的时间变长,同时50Hz对应的衰减幅度增加,衰减的幅度值和理论推导基本一致。另外,当a=0.95时,DA输出的均值变小。这主要是进行循环迭代运算时,需要将16位的变量转化为8位表示形式所导致的。在有用信号失真较小的情况下,为使滤波器达到降低工频干扰的最佳效果,必须恰当选择a 值。经过以上的仿真试验可以发现,当a=0.9 时,衰减幅度、DA输出均值和算法收敛时间表现比较均衡,可以作为一般情况下的选择值。

将VMLAB 中虚拟示波器的显示数据导出到一个*.cvs文件中,用matlab 读出这些数据,并画出不同a值对应的输出响应,如图4 。从图4 可以清晰看出不同a值下算法的性能变化的大致走向。

将AD的采样间隔设置为4ms ,对应的采样频率Fs就变为250Hz ,其它条件不变。通过VMLAB 进行仿真,对比结果如表4 、图5。

对比Fs = 500Hz 的情况,随着采样频率Fs降低,50Hz频率的幅度衰减值会逐渐增加。这主要是因为随着采样频率降低,低通滤波器的截至频率fc也随之降低,相应的滤波器在50Hz处的衰减也就越来越低。根据奈奎斯特低通采样定理,当采样频率小于100Hz时,由于信号频谱混叠,滤波器对50Hz信号的滤波效果将会变差。如果只是对缓变信号进行采样,采样频率比100Hz稍大即可。但是随着采样频率的降低,滤波算法的收敛时间也会增加。因此必须在算法的滤波性能和收敛时间上进行折衷考虑。

本文提出的分配系数法设计数字滤波器,算法速度快、代码效率高、滤波效果理想,是一种实用的数字滤波器设计方法,体现了将算法嵌入到具体硬件的思想。另一方面,将定点小数的表示形式进行适当扩展,这个算法还可以用于10位或16位AD转换精度的应用场合。

关键字:AVR单片机  数字滤波器  工频干扰 引用地址:基于AVR单片机的数字滤波器滤除工频干扰的快速算法

上一篇:可编程数字电位器与AVR单片机的通信说明
下一篇:AVR 内部EEPROM读写式例

推荐阅读最新更新时间:2024-03-16 13:31

AVR单片机全系列性能参数表
AVR单片机全系列性能参数表 包括:tiny11,tiny12,tiny13,tiny15L,tiny2313,tiny26,tiny26L,tiny28L,tiny28V,90S1200, 90S2313 Mega128,Mega128L,Mega16,Mega16L,Mega162,Mega162L,Mega162V, Mega169,Mega169L,Mega169V,Mega32,Mega32L 包括:Mega48,Mega88,Mega168,Mega64,Mega64L,Mega8,Mega8L,Mega8515, Mega8515L,Mega8535,Mega8535L
[单片机]
<font color='red'>AVR单片机</font>全系列性能参数表
AVR单片机电梯外呼板系统设计
随着现代社会的发展,科学技术的进步,出现了众多高层建筑和智能建筑。电梯,作为高层建筑内部一种重要的交通工具,其应用规模日益扩大。而作为电梯系统中必不可少的一部分,电梯召唤显示板(简称电梯外呼板)的应用也随之剧增。 电梯外呼板应用于每层楼的电梯门外,供乘客及电梯维保人员使用。电梯外呼板将乘客及维保人员的需求信息通过CAN总线传达给电梯主板,电梯主板接收信息并执行相应的操作。同时,电梯主板将电梯的实时运行信息通过CAN总线传递给电梯外呼板,通过外呼板LED显示出来,供乘客参阅。 AVR单片机具有高可靠性、功能强、高速度、低功耗和低价位的特点,本系统选用了高档ATmega列AVR单片机ATmegal6。它具有先进的RI-SC结构
[单片机]
<font color='red'>AVR单片机</font>电梯外呼板系统设计
AVR单片机用于太阳能电池控制器的设计
简介:在此设计的太阳能控制器性能稳定,具有过充过放保护和温度补偿。经过测试,系统显示出良好的控制效果,不仅提高了太阳电池的工作效率,同时也保护了所使用的蓄电池,在利用绿色能源方面,具有一定的社会效益和广泛的推广价值。 随着能源危机和环境污染的加深,太阳能的研究和利用受到广泛的关注。太阳能是人类取之不尽用之不竭的可再生能源,也是清洁能源,不产生任何的环境污染,在太阳能的有效利用中,太阳能充电是近些年发展最快,最具活力的研究领域,是其中最受瞩目的项目之一。太阳能电池发电是基于“光生伏打效应”原理,将太阳能转化为电能,利用充电效应将太阳辐射直接转化为电能。它具有永久性、清洁性和灵活性大的优点,是其他能源无法比拟的。 1 太阳能控
[单片机]
<font color='red'>AVR单片机</font>用于太阳能电池控制器的设计
基于AVR单片机的高精度频率调节器
1?引言? 在现代化工业生产与高精度测试中,我们需要相当精确的频率来帮助判断设备性能指标。而且我们希望能够微调该频率。采用压控振荡器得到的频率不够精确,微调频率步骤烦琐,耗时漫长,因此有些测试项目限制了压控振荡器的使用。 Mega系列单片机是Atmel公司于2002年起陆续推出的。这款AVR增强型单片机具有速度快,抗干扰能力强,价格低廉等诸多优点。为了加快AVR单片机的软件编制,Atmel以及第三方提供的开发工具多种多样,程序开发方便有效。该单片机内部FLASH结构功能灵活,加锁后很难解密,可以最大限度地保护知识产权。AVR单片机可以广泛应用于通信、野外测试、汽车电子、医疗器械等领域,并且适用于各种低电压、低功耗的场合。? 本文
[单片机]
AVR单片机程序设计架构
昨天在回家的火车上,带了一本《匠人手记》,看了几篇,受益匪浅。其中一篇讲到编程思路,也就是如何入手构建一个程序。我用C语言在原文的基础上稍微展开一下,以备日后查看。 简单的程序只要用简单的方式就可以了,我在此之前也只会这一种方式。也就是: void main() { init();//初始化 while (1) { module01();//模块1,或若干语句 } } 匠人在书中讲到事件驱动机制,例如: void main() { init(); while (1) { if (事件1条件) shijian1();//执行事件1 if (事件2条件)
[单片机]
AVR单片机晶振电路设计
Mega16已经内置RC振荡线路,可以产生1M、2M、4M、8M的振荡频率。不过,内置的毕竟是RC振荡,在一些要求较高的场合,比如要与RS232通信需要比较精确的波特率时,建议使用外部的晶振线路。 早期的90S系列,晶振两端均需要接22pF左右的电容。Mega系列实际使用时,这两只小电容不接也能正常工作。不过为了线路的规范化,我们仍建议接上。 重要说明:实际应用时,如果你不需要太高精度的频率,可以使用内部RC振荡。即这部分不需要任何的外围零件。
[单片机]
<font color='red'>AVR单片机</font>晶振电路设计
基于AVR单片机的新型防汽车追尾安全装置设计
1 引言 随着我国汽车保有量逐年增加,高速公路的发展使得汽车的行驶速度越来越快。车流量也越来越大,导致高速公路交通事故频频发生,其发生率甚至达到普通公路的4倍,且有上升趋势。在高速公路所发生的交通事故中尤以汽车追尾居多,约占30%~40%,而在这些事故中,80%以上的事故是由于司机反应不及时或者判断失误引起的。由于目前只有高档车型才配有防碰撞装置,而且以往设计中只单纯考虑在车前或车后安装防碰撞装置,不尽全面,所以设计一套低成本、通用性好的汽车防车前车后追尾的安全装置具有广阔的应用和市场前景。因此。这里提出一种基于AVR的新型防汽车追尾安全装置设计。 2 安全装置组成和硬件电路设计 安全装置硬件电路设计主要由单片
[单片机]
小广播
添点儿料...
无论热点新闻、行业分析、技术干货……
设计资源 培训 开发板 精华推荐

最新单片机文章
何立民专栏 单片机及嵌入式宝典

北京航空航天大学教授,20余年来致力于单片机与嵌入式系统推广工作。

换一换 更多 相关热搜器件
电子工程世界版权所有 京B2-20211791 京ICP备10001474号-1 电信业务审批[2006]字第258号函 京公网安备 11010802033920号 Copyright © 2005-2024 EEWORLD.com.cn, Inc. All rights reserved