浅谈数字滤波抗干扰技术有哪些应用？

　　1 引言
　　仪器仪表设备在现场测试过程中，由于生产变量的测试数据对生产过程具有重要的意义，因此对各种物理量测试数据精度要求是比较高的。
　　在前向测试通道上采用的抗干扰措施中，滤波方法是抑制干扰的一种有效途径。在工业现场中，可利用硬件滤波器电路或软件滤波器算法提高测试数据的准确性。硬件滤波措施是使用较多的一种方法，技术比较成熟，但同时也增加了设备，提高了成本，而且电子设备的增加有可能带来新的干扰源。而采用软件滤波算法不需增加硬件设备，可靠性高，功能多样，使用灵活，具有许多硬件滤波措施所不具备的优点，当然它需要占一定的运行时间。
　　2 常用的几种软件滤波方法
　　（1）中值滤波法：即每次取N个AD值，去除其中的最大值和最小值而取剩余的N-2个A/D转换值的平均值。
　　（2）程序判断滤波法：即根据经验确定出两次采样的最大偏差ΔY，若先后两次采样的信号相减数值大于ΔY，表明输入为干扰信号，应去除；用上次采样值与本次采样值比较，若小于或等于Δ Y，表明没有受到干扰，此时本次采样值有效，这样可以滤去随机干扰和传感器不稳定而引起的误差。
　　（3）利用格拉布斯（Grubbs）准则进行处理：根据误差理论，要有效地剔除偶然误差，一般要测量10次以上，兼顾到精度和响应速度，取15次为一个单位。在取得的15个数据中，有些可能含有较大的误差，需要对它们分检，剔除可疑值，提高自适应速度。对可疑值的剔除有多种准则，如莱以达准则、肖维勒（Chauvenet）准则、格拉布斯（Grubbs）准则等。以Grubbs准则为例，它认为若某测量值 xi对应的残差Vi满足下式
　　|Vi|=| xi-|》g（n，a）× σ（X）
　　时应将该数据舍去。式中，为n次采集到的AD 值的平均值，=（∑xi）/n ；σ（X）为测量数据组的标准差，由贝塞尔函数可得： σ（X）=［（∑Vi2 ）/（n-1）］1/2；g（n， a）是取决于测量次数n和显著性水平a （相当于犯“弃真” 错误的概率系数），a通常取0.01或0.05。通过查表可得：当 n=15时，a=0.05， g（n，a）=2.41。
　　把15次采集到的AD值存入一个数组中然后求平均值，计算残差，求标准差σ（X）。将残差绝对值与2.41倍的标准差σ（X）比较。剔除可疑值以后，再求平均值，求出新的平均值以后，应再重复以上过程，验证是否还有可疑值存在。据我们对测量装置大量的实际测试结果看，这样做没有什么必要，因为一般只有第一遍即可达到要求。
　　然而这种方法也有它的不足， 利用Grubbs准则需要处理大量的数据，而在一般的工业现场测试设备中，仪表结构大多采用嵌入式结构，如AVR单片机。这些MCU程序空间和数据空间有限，若处理大量数据，难以满足资源要求。而且，由于Grubbs准则要求MCU进行大量数据处理，使得系统降低了信号采集速率，影响实时性。
　　3 AD7705内置的数字滤波器
　　上述的几种软件滤波方法对抑制个别的异常数据方面具有一定作用，然而对于由工频干扰引起的平稳随机干扰信号不太理想。随着近几年来集成电路的快速发展，一些数字滤波器也集成到数据处理芯片当中，而且数字滤波器的设计参数通过软件可编程实现，具有极大的灵活性和实用性。
　　3.1 芯片介绍及应用
　　目前市场上由美国AD公司新近推出的AD7705 就是一例。它是内置了数字滤波器的带信号调理电路的16位A/D转换器件，可应用于低频测量的2通道的模拟前端。该器件可以接受直接来自传感器的低电压输入信号，然后产生串行的数字输出，利用∑-Δ转换技术实现16位无丢失代码性能。选定的输入信号被送到一个基于模拟调制器的增益可编程专用前端，片内数字滤波器处理调制器的输出信号，通过片内控制寄存器可调节滤波器的截止点和输出更新率，从而对数字滤波器的第一个陷波进行编程。
　　AD7705有以下几个特点：
　　（1）差分输入通道ADC，具有16位无丢失代码和0.003%非线形误差；
　　（2）MCU的接口为串行的三线接口模式；
　　（3）可编程增益1～128，信号极性以及更新速率配置有工作寄存器，通过对寄存器的操作可以选择等；
　　（4） 器件还包括自标定和系统校准功能，以消除器件本身或系统的增益和编程误差；
　　（5） 内带数字信号处理电路，由软件可编程进行设置。
　　引脚排列如图1。
　　
　　主要管脚说明如下：
　　MCLK IN：为转换器提供主时钟信号，能以晶体/谐振器或外部时钟的形式提供。
　　MCLK OUT：当主时钟为晶体/谐振器时，晶体谐振器被接在MCLK IN和MCLK OUT之间；如果在MCLK IN引脚处连接一个外部时钟，MCLK OUT将提供一个反向时钟信号。
　　CS ：片选，低电平有效的逻辑输入。
　　A1N2（+）［A1N1］：对于AD7705，差分模拟输入通道2的正输入端；对于AD7706，模拟输入通道1的输入端。
　　A1N1（+）［A1N2］：差分模拟输入通道1的正输入端。
　　A1N1（-）［COMMON］：差分模拟输入通道1的负输入端。
　　A1N2（-）［A1N3］：差分输入模拟通道2的负输入端。
　　DRDY ：逻辑输出。这个输出端上的逻辑低电平表示可从AD7705的数字寄存器获取新的输出字。
　　DOUT：串行数据输出端。
　　DIN：串行数据输入端。
　　SCLK：串行时钟，施密特逻辑输入。
　　从AD转换器件发展趋势来看，目前的16位A/ D转换器除具有较高的转换精度外，还具有高度的集成化、可编程等特点 ［2，3］。AD7705符合这些特点，其简化的应用电路如图2。
　　
　　图2给出了一个完整的前向测试通道，图中 AD7705可直接接收传感器输入信号，而不用外接信号调理电路，模拟量输出增益可实现软件可编程调节。输出的AD值通过SPI三线口模式发送给MCU。图中，MCU采用的是P87LPC764，它是 PHILIPS公司新推出的增强型8位单片机，内核与MCS51系列兼容，内置电源检测、看门狗WDT、 IIC总线等资源。
　　AD7705的数据接口采用了SPI接口标准。SPI 接口利用4根线可完成主从之间的数据通信。接口非常方便。这4根口线分别为时钟线（SCLK）、数据输入线（SDI）、数据输出线（SDO）和片选线（CS）。DIN线用来向片内寄存器传输数据，而DOUT线用来访问寄存器里的数据。SCLK是串行时钟输入，所有的数据传输都和SCLK信号有关。DRDY线作为状态信号，以提示数据什么时候已准备好从寄存器读数据。输出寄存器中有新的数据时，DRDY变为低电平，在数据寄存器更新前，若DRDY变为高电平，则提示这个时候不读数据，以免在寄存器更新的过程中读数据。CS作为片选信号，在本电路中由于只用一片AD7705故CS端接地。
　　3.2 低通滤波及陷波处理
　　AD7705包含一个片内低通滤波器，用它处理器件的∑-Δ调制器的输出信号。所以，该器件不仅提供模数转换功能，而且还具备一定的滤波能力。AD7705内嵌了数字低通滤波器，数字滤波发生在模-数转换之后，它能消除模-数转换过程中产生的噪音。此外，数字滤波器容易实现可编程性。依靠数字滤波器设计，用户可以编程截断频率和输出更新频率。
　　AD7705的数字滤波器是一个低通（sin x/x）3 滤波器（也称为sinc3），其在Z 域的转换功能可描述为，
　　
　　在3次域可描述为，
　　
　　图3是截止频率为15.72Hz的滤波器频率响应，这一频率与滤波器的第一个陷波频率（50Hz）对应，这曲线显示从直流到300Hz。在数字滤波器采样频率的两边这个频率曲线各重复一次。这种滤波器响应与一个平均滤波器的相似。因此，图中的输出速率为50Hz，滤波器的第一个陷波频率为50Hz。这种（sin x/x）3 滤波器的陷波位置多次在第一个陷波的倍数处重复。在这些陷波处滤波器提供大于100dB 的衰减。
　　
　　AD7705内部带有数字滤波器，具有抑制干扰功能。对于50Hz的工频干扰在第一陷波位置已衰减180dB，这样有效抑制工频干扰。另外，我们在每次采集AD值时，每采集10个AD值作为一组数据，去除其中的最大值，最小值而取剩余的8个A/D转换值的平均值，这样可有效去除数据随机误差。