CN101188410A - 用于智能变送器的一阶滞后数字滤波算法 - Google Patents
用于智能变送器的一阶滞后数字滤波算法 Download PDFInfo
- Publication number
- CN101188410A CN101188410A CNA2007101708756A CN200710170875A CN101188410A CN 101188410 A CN101188410 A CN 101188410A CN A2007101708756 A CNA2007101708756 A CN A2007101708756A CN 200710170875 A CN200710170875 A CN 200710170875A CN 101188410 A CN101188410 A CN 101188410A
- Authority
- CN
- China
- Prior art keywords
- formula
- digital filtering
- filtering algorithm
- value
- unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Measuring Fluid Pressure (AREA)
Abstract
一种用于智能变送器的一阶滞后数字滤波算法,涉及信号处理技术领域;所要解决的是智能变送器的输出稳定、采样值无丢失、平滑性好的技术问题;该数字滤波算法的步骤如下:1)与前一次计算值相减,按照(式1)),计算出变化量ΔYn;2)按照(式2))的函数关系f(ΔY)对滤波时间Tf进行计算;3)然后按照(式3)),计算滤波系数a;4)按照(式4)),计算最后的采样滤波值n;本发明的一阶滞后数字滤波算法能使智能变送器在全量程范围内输出稳定、采样值无丢失、平滑性好及输出具有快速响应的特性。
Description
技术领域
本发明涉及信号处理技术,特别是涉及一种适用于带微处理器的工业测量用智能变送器的,利用汇编或高级语言编写的一阶滞后数字滤波算法的技术。
背景技术
在工业过程控制***中,大量的现场物理信号(压力、温度、流量等)需要被测量,这主要是通过变送器来完成的。变送器包括传感器和电子线路两个部分,传感器实现物理量的测量并转换为电信号,电子线路将传感器输出的电信号进行处理,转变为标准电信号(如4~20mA、1~5V等)并输出。目前广泛应用的智能变送器的电子线路部分是以微处理器为数据处理核心,还包括了信号处理、A/D转换、存储器、D/A转换、数字通信等单元,如图1所示。
传感器输出的模拟信号需要进行滤波,智能变送器一般使用的是一阶滞后数字滤波算法,如式①所示:
式中: Yn——未经滤波的第n次采样值;
a——滤波系数,在区间[0,1]内取值。
式①中的滤波系数a,表示滤波时间常数与采样周期的比值,即为式②:
式中:Tf——滤波时间常数;
T——采样周期;
当采样周期T远小于滤波时间常数Tf的时候,该算法等价于一般的一阶滞后RC模拟滤波器,对周期性干扰具有良好的抑制作用,适用于波动频率较高参数的滤波。其不足之处是带来了相位滞后,灵敏度低。滞后的程度取决于滤波时间常数与采样周期的比值(即a):a越大,滞后越多。实际应用的时候,常常根据经验设定滤波时间常数Tf,即阻尼调整。
随着对变送器要求的不断提高,希望变送器在大量程比、小信号情况下仍保持稳定输出,需要把滤波时间常数Tf调大,即增大a;但又希望保持输出的快速响应,即滞后较小(a减小),因此普通的一阶滞后滤波算法无法满足需要。一些智能变送器采用小信号切除技术,使得在零位或部分输出点保持稳定,但这种技术不能保证全量程较高的稳定性,同时会使输出平滑性不好。
发明内容
针对上述现有技术中存在的缺陷,本发明所要解决的技术问题是提供一种对模拟信号具有一阶滞后滤波器的良好抑制作用,在全量程范围内输出稳定、采样值无丢失、平滑性好及输出具有快速响应特性的用于智能变送器的一阶滞后数字滤波算法。
为了解决上述技术问题,本发明所提供的一种用于智能变送器的一阶滞后数字滤波算法,如式①所示:
式中:Yn——未经滤波的第n次采样值;
a——滤波系数,在区间[0,1]内取值。
式①中的滤波系数a,表示滤波时间常数与采样周期的比值,即为式②:
式中:Tf——滤波时间常数;
T——采样周期;
Tf=f(ΔY)
利用式③通过变化量ΔYn计算得到的滤波时间常数Tf,再代入到式②中计算滤波系数a;最后代入通用的一阶滞后数字滤波算法,即式①,计算得到滤波后的输出Yn。
进一步的,所述智能变送器的结构为:变送器的被测信号输入端经传感器单元、信号调理单元、A/D转换单元连接微处理器单元的信号输入端,微处理器单元的信号输出端经D/A转换单元连接变送器的信号输出端,微处理器单元分别连接通信单元和存储单元,传感器单元连接信号激励单元的输出端,电源模块连接并供电给变送器的所述其它电路单元;信号激励单元产生激励信号,并作用于传感器单元;传感器单元对被测信号进行测量,经过信号调理单元,进行采样和A/D转换单元,并将数字采样值送给微处理器单元;微处理器单元根据采样值进行滤波、校准等数值计算,即一阶滞后数字滤波算法的计算,并计算出输出信号的大小,利用D/A转换单元进行实际的信号输出;通信单元实现微处理器与控制站的双向数字传输,存储单元保存智能变送器的组态参数。
进一步的,对于金属电容式智能差压变送器,算法中所述滤波时间Tf与ΔY之间的设定的函数关系采用式④所描述的分段离散函数来实现:
式中:TH——较大的滤波时间常数;
TL——较小的滤波时间常数;
TNorm——根据经验预设的滤波时间;
Δ1——需要增强滤波的变化量上限值;
Δ2——需要减弱滤波的变化量下限值。
为了解决上述技术问题,本发明所提供的一种智能变送器的计算软件流程,其特征在于,程序的步骤如下:
1)首先进行初始化,然后进入循环主程序;
2)循环主程序中首先进行AD采样;
3)与前一次计算值相减,按照 ,计算出变化量ΔYn;
4)按照Tf=f(ΔY),的函数关系f(ΔY)对滤波时间Tf进行计算;
5)然后按照 ,计算滤波系数a;
7)再进行其它(包括非线性修正、电流计算等)计算和最终的输出循环至步骤2)。
本发明的用于智能变送器的一阶滞后数字滤波算法,具有以下优点:
1)算法是对原有一阶滞后数字滤波算法的改进,较其他复杂的数字滤波算法简单,容易在微处理器内部实现;
2)使用该算法设计的智能变送器,具有快速响应、全量程范围内输出稳定、采样值无丢失、平滑性好等优点。
附图说明
图1是智能变送器的结构示意图;
图2是本发明算法中分段函数的示意图;
图3是采用本发明算法后的智能变送器软件流程图;
图4是本发明实施例中分段函数的示意图。
具体实施方式
以下结合附图说明对本发明的实施例作进一步详细描述,但本实施例并不用于限制本发明,凡是采用本发明的相似结构及其相似变化,均应列入本发明的保护范围。
图1是智能变送器的结构示意图。变送器的被测信号输入端经传感器单元103、信号调理单元104、A/D转换单元105连接微处理器单元108的信号输入端,微处理器单元108的信号输出端经D/A转换单元106连接变送器的信号输出端,微处理器单元108分别连接通信单元107和存储单元109,传感器单元103连接信号激励单元102的输出端,电源模块101首先从回路中获得电能,连接并供电给变送器的所述其它电路单元;信号激励单元102产生激励信号,并作用于传感器单元103;传感器单元103对被测信号进行测量,经过信号调理单元104,进行采样和A/D转换单元105,并将数字采样值送给微处理器单元108;微处理器单元108根据采样值进行滤波、校准等数值计算,并计算出输出信号的大小,利用D/A转换单元106进行实际的信号输出;通信单元107实现微处理器与控制站的双向数字传输,存储单元109保存智能变送器的组态参数。
受被测信号波动和外部干扰的影响,传感器输出的模拟信号不够稳定,存在周期性波动,需要进行滤波,智能变送器一般使用的是一阶滞后数字滤波算法,见式①。当T<<Tf时,该算法便等价于一般的模拟滤波器。一阶滞后滤波算法对周期性干扰具有良好的抑制作用,适用于波动频率较高参数的滤波。其不足之处是带来了相位滞后,灵敏度低。滞后的程度取决于a值的大小。同时,它不能滤出频率高于采样频率二分之一的干扰信号。在智能变送器***中,通过设定滤波时间Tf,即阻尼时间,来调整输出的稳定性。
本发明实施例所提供的一种用于智能变送器的一阶滞后数字滤波算法:首先根据当前测量值Yn与上一次计算值的变化量绝对值ΔYn,按照某种函数关系对滤波时间常数Tf进行修正,即Tf是一个关于ΔY的函数,如式③所示:
Tf=f(ΔY)
利用式③通过变化量ΔYn计算得到的滤波时间常数Tf,再代入到式②中计算滤波系数a;最后代入通用的一阶滞后数字滤波算法,即式①,计算得到滤波后的输出Yn。
通常采用的一阶滞后数字滤波算法的滤波时间常数Tf是固定的,其大小根据输入信号和噪声的强弱按照经验进行选择。改进后的一阶滞后滤波算法,其滤波时间常数Tf是关于采样值变化量ΔY的函数f(ΔY),该函数可以是任意形式的一元函数表达式,包括线性函数、非线性函数以及离散函数等。该函数一般情况下具有如下特征:在采样值变化不大的情况下(ΔY较小),表示输入信号处于某一平衡状态,此时滤波时间常数Tf取值较大,使输出较稳定;当采样值变化较明显的时候(ΔY较大),表示输入信号阶跃变化,滤波时间常数Tf取值较小,使输出随输入迅速变化。总而言之,Tf与ΔY存在反比例关系。
本发明的用于智能变送器的一阶滞后数字滤波算法,首先基于采样值变化量ΔY,利用函数f(ΔY)使滤波时间Tf随ΔY动态调整。该函数具有输出与输入成反比例的特征,可以表示为图2。在采样值变化不大的情况下(ΔY较小),表示输入信号处于某一平衡状态,此时滤波时间常数Tf取值较大,使输出较稳定;当采样值变化较明显的时候(ΔY较大),表示输入信号阶跃变化,滤波时间常数Tf取值较小,使输出随输入迅速变化。
图3是采用本发明算法后的智能变送器软件流程图。程序开始以后:
1)首先进行初始化,然后进入循环主程序;
2)循环主程序中首先进行AD采样;
3)与前一次计算值相减,计算出变化量ΔYn,
4)按照式③的函数关系f(ΔY)对滤波时间Tf进行计算,Tf=f(ΔY);
5)然后按照式②计算滤波系数a,
7)再进行其它(包括非线性修正、电流计算等)计算和最终的输出循环至步骤2)。
实施例:
下面结合实施例对本实用新型作进一步的说明。
实施例:对于金属电容式智能差压变送器,算法中滤波时间Tf与ΔY之间的函数关系采用式④所描述的分段离散函数来实现。
式中:TH——较大的滤波时间常数;
TL——较小的滤波时间常数;
TNorm——根据经验预设的滤波时间;
Δ1——需要增强滤波的变化量上限值;
Δ2——需要减弱滤波的变化量下限值;
Claims (3)
1.一种用于智能变送器的一阶滞后数字滤波算法,如式①所示:
式中:Yn——未经滤波的第n次采样值;
α——滤波系数,在区间[0,1]内取值。
式①中的滤波系数α,表示滤波时间常数与采样周期的比值,即为式②:
式中:Tf——滤波时间常数;
T——采样周期;
Tf=f(ΔY)
利用式③通过变化量ΔYn计算得到的滤波时问常数Tf,再代入到式②中计算滤波系数α:最后代入通用的一阶滞后数字滤波算法,即式①,计算得到滤波后的输出Yn。
2.根据权利要求1所述的用于智能变送器的一阶滞后数字滤波算法,其特征在于,所述智能变送器的结构为:变送器的被测信号输入端经传感器单元、信号调理单元、A/D转换单元连接微处理器单元的信号输入端,微处理器单元的信号输出端经D/A转换单元连接变送器的信号输出端,微处理器单元分别连接通信单元和存储单元,传感器单元连接信号激励单元的输出端,电源模块连接并供电给变送器的所述其它电路单元。
3.根据权利要求1所述的用于智能变送器的一阶滞后数字滤波算法,其特征在于,对于金属电容式智能差压变送器,算法中所述滤波时间Tf与ΔY之间的设定的函数关系采用式④所描述的分段离散函数来实现:
式中:TH——较大的滤波时间常数;
TL——较小的滤波时间常数;
TNorm——根据经验预设的滤波时间;
Δ1——需要增强滤波的变化量上限值;
Δ2——需要减弱滤波的变化量下限值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2007101708756A CN101188410A (zh) | 2007-11-27 | 2007-11-27 | 用于智能变送器的一阶滞后数字滤波算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2007101708756A CN101188410A (zh) | 2007-11-27 | 2007-11-27 | 用于智能变送器的一阶滞后数字滤波算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101188410A true CN101188410A (zh) | 2008-05-28 |
Family
ID=39480635
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2007101708756A Pending CN101188410A (zh) | 2007-11-27 | 2007-11-27 | 用于智能变送器的一阶滞后数字滤波算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101188410A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102536345A (zh) * | 2012-01-13 | 2012-07-04 | 杭州哲达科技股份有限公司 | 一种提高bprt***运行能效的方法 |
CN104617914A (zh) * | 2015-02-11 | 2015-05-13 | 珠海格力电器股份有限公司 | 一种电器设备的信号滤波方法及*** |
CN104660215A (zh) * | 2015-02-26 | 2015-05-27 | 沈阳东北电力调节技术有限公司 | 一种低资源消耗的数字滤波方法 |
CN104883158A (zh) * | 2015-05-18 | 2015-09-02 | 北汽福田汽车股份有限公司 | 电动汽车加速踏板信号的滤波方法及滤波装置 |
CN105571690A (zh) * | 2014-10-09 | 2016-05-11 | 梅特勒-托利多(常州)精密仪器有限公司 | 数字称重传感器及传感器网络 |
CN107560645A (zh) * | 2017-08-29 | 2018-01-09 | 北京航空航天大学 | 一种光纤布拉格光栅传感器波长解调寻峰方法 |
CN110209628A (zh) * | 2019-06-05 | 2019-09-06 | 杭州华塑加达网络科技有限公司 | 一阶滞后滤波生成装置及方法 |
CN110342405A (zh) * | 2019-07-25 | 2019-10-18 | 上海振华重工(集团)股份有限公司 | 一种集装箱跨运车起升机构精确定位的控制方法 |
CN111216566A (zh) * | 2020-02-21 | 2020-06-02 | 威马智慧出行科技(上海)有限公司 | 一种车辆电机扭矩的控制方法及装置 |
-
2007
- 2007-11-27 CN CNA2007101708756A patent/CN101188410A/zh active Pending
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102536345B (zh) * | 2012-01-13 | 2014-11-12 | 杭州哲达科技股份有限公司 | 一种提高bprt***运行能效的方法 |
CN102536345A (zh) * | 2012-01-13 | 2012-07-04 | 杭州哲达科技股份有限公司 | 一种提高bprt***运行能效的方法 |
CN105571690B (zh) * | 2014-10-09 | 2018-09-07 | 梅特勒-托利多(常州)精密仪器有限公司 | 数字称重传感器及传感器网络 |
CN105571690A (zh) * | 2014-10-09 | 2016-05-11 | 梅特勒-托利多(常州)精密仪器有限公司 | 数字称重传感器及传感器网络 |
CN104617914B (zh) * | 2015-02-11 | 2018-09-07 | 珠海格力电器股份有限公司 | 一种电器设备的信号滤波方法及*** |
CN104617914A (zh) * | 2015-02-11 | 2015-05-13 | 珠海格力电器股份有限公司 | 一种电器设备的信号滤波方法及*** |
CN104660215A (zh) * | 2015-02-26 | 2015-05-27 | 沈阳东北电力调节技术有限公司 | 一种低资源消耗的数字滤波方法 |
CN104660215B (zh) * | 2015-02-26 | 2017-08-25 | 沈阳东北电力调节技术有限公司 | 一种低资源消耗的数字滤波方法 |
CN104883158B (zh) * | 2015-05-18 | 2018-03-27 | 北汽福田汽车股份有限公司 | 电动汽车加速踏板信号的滤波方法及滤波装置 |
CN104883158A (zh) * | 2015-05-18 | 2015-09-02 | 北汽福田汽车股份有限公司 | 电动汽车加速踏板信号的滤波方法及滤波装置 |
CN107560645A (zh) * | 2017-08-29 | 2018-01-09 | 北京航空航天大学 | 一种光纤布拉格光栅传感器波长解调寻峰方法 |
CN110209628A (zh) * | 2019-06-05 | 2019-09-06 | 杭州华塑加达网络科技有限公司 | 一阶滞后滤波生成装置及方法 |
CN110342405A (zh) * | 2019-07-25 | 2019-10-18 | 上海振华重工(集团)股份有限公司 | 一种集装箱跨运车起升机构精确定位的控制方法 |
CN110342405B (zh) * | 2019-07-25 | 2020-10-02 | 上海振华重工(集团)股份有限公司 | 一种集装箱跨运车起升机构精确定位的控制方法 |
CN111216566A (zh) * | 2020-02-21 | 2020-06-02 | 威马智慧出行科技(上海)有限公司 | 一种车辆电机扭矩的控制方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101188410A (zh) | 用于智能变送器的一阶滞后数字滤波算法 | |
CN103256044B (zh) | 一种随钻声波信号处理装置 | |
CN111769819B (zh) | 一种采样频率自适应可调的数据采集方法及*** | |
CN108429286A (zh) | 一种基于自抗扰控制的并网电流调节器 | |
CN101141162A (zh) | 一种射频模块性能的软硬件补偿方法及改进的射频模块 | |
CN105958470A (zh) | 一种基于双线性凸优化理论电力***双线性抗差估计方法 | |
CN110289821A (zh) | 适应工况环境的射频信号输出功率控制电路、方法及装置 | |
CN104895611A (zh) | 一种测量矿井粉尘浓度的智能传感器 | |
CN201498355U (zh) | 电子式电流互感器的模拟信号处理装置 | |
CN204302377U (zh) | 一种电流测量装置及智能断路器 | |
CN107863940A (zh) | 一种带温补功能的大动态范围alc保护电路 | |
CN103487778B (zh) | 一种模拟量校正方法及基于该方法的单cpu低压保护装置 | |
CN101813959B (zh) | 一种功率因数校正电路的控制装置及控制方法 | |
CN102769436A (zh) | 适用于便携式仪器设备的低成本微功耗微型化的信号调理模块 | |
CN102801151B (zh) | 电感量连续可调型消弧线圈的优化调谐***及方法 | |
CN106291096A (zh) | 一种用于非线性负荷电能表的多级重采样方法 | |
CN102394657A (zh) | 实时电能质量数据压缩算法 | |
CN106556724A (zh) | 电表及其时钟校准方法 | |
CN116953347A (zh) | 一种提高基波计量精度的方法及装置 | |
CN203301361U (zh) | 光伏并网逆变器的智能电流及电压隔离变送器 | |
CN109507624A (zh) | 一种电容模拟器及模拟电容的方法 | |
CN101902213B (zh) | 一种信号衰减控制方法及衰减式平衡控制电路 | |
CN2755839Y (zh) | 基于裂芯式或磁阀式可控电抗器的部分可调谐滤波装置 | |
CN201178414Y (zh) | 一种利用软件补偿方法改进的射频模块 | |
CN218496322U (zh) | 一种运算放大器为主的动态压力传感器变送器电路 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20080528 |