CN111181526B - 一种用于信号处理的滤波方法 - Google Patents

一种用于信号处理的滤波方法 Download PDF

Info

Publication number
CN111181526B
CN111181526B CN202010005848.9A CN202010005848A CN111181526B CN 111181526 B CN111181526 B CN 111181526B CN 202010005848 A CN202010005848 A CN 202010005848A CN 111181526 B CN111181526 B CN 111181526B
Authority
CN
China
Prior art keywords
data set
filtering method
data
signal processing
filtering
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.)
Active
Application number
CN202010005848.9A
Other languages
English (en)
Other versions
CN111181526A (zh
Inventor
王鹏鹏
聂一雄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN202010005848.9A priority Critical patent/CN111181526B/zh
Publication of CN111181526A publication Critical patent/CN111181526A/zh
Application granted granted Critical
Publication of CN111181526B publication Critical patent/CN111181526B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0201Wave digital filters
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0211Frequency selective networks using specific transformation algorithms, e.g. WALSH functions, Fermat transforms, Mersenne transforms, polynomial transforms, Hilbert transforms

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Algebra (AREA)
  • Analogue/Digital Conversion (AREA)
  • Emergency Protection Circuit Devices (AREA)

Abstract

本申请提供了一种用于信号处理的滤波方法,通过将模拟信号转换为数字信号后得到原始数据组,根据原始数据组形成差值数组,将差值数组按计算规则进行多次迭代计算,对数据进行滤波处理。本申请的滤波方法能去除各种干扰造成的波动,使滤波后的数据得到较好的平滑度,并且数据不失真,仍在前后时刻的变化范围内,最大限度接近实际值,而且本申请的方法运算量小,适合大部分实时性要求较高的场合。

Description

一种用于信号处理的滤波方法
技术领域
本发明属于信号处理技术领域,具体涉及一种用于信号处理的滤波方法。
背景技术
随着电子技术与数字技术的飞速发展,信号的传递速度与CPU的处理速度越来越快,因此对数据采集和处理的要求也越来越高。模拟量是指变量在一定范围内连续变化的量,分为电流型和电压型,经常用于单片机或PLC中。由于在单片机主程序结束一个扫描周期时才更新对模拟量的采集,因此会出现采集数据异常突变的现象。而且电压型和电流型模拟量数据都容易受到外部电路的干扰,因此需要对采集的模拟量进行一定的滤波处理。现有技术中对采集到的模拟量数据进行处理使用的软件滤波方法主要分为两种,第一种为除极限值法,主要思想为通过比较计算把最大值和最小值舍去,把剩余值作为有效采集值进行处理;第二种为平均值法,主要思想为对一定时间内所有采集到的数据进行求平均值,以达到滤波效果。其它方法都是这两种基本方法的衍生方法。
在以上方法中,除极限值法能够有效克服偶然因素引起的波动,但是无法抑制周期性的干扰,而且平滑度差。平均值法不适用于检测速度慢或者实时性要求高的场合。在自然界中,有些变量实际上的确是在突变的,而且没有规律,一旦除去极限值或求平均值,可能会使处理后的数据远离模拟量的真实范围,出现人为造成的失真。因此提供一种能去除各种干扰造成的波动而且数据不会失真的用于信号处理的滤波方法尤为重要。
发明内容
基于此,本发明旨在提供一种用于信号处理的滤波方法,通过对原始数组的差值数组按计算规则进行多次迭代计算,实现对数据进行滤波处理。用于解决现有技术存在的问题,通过本发明的滤波方法去除各种干扰造成的波动,得到较好的平滑度,而且数据不失真,仍在前后时刻的变化范围内,最大限度接近实际值,而且本发明的方法运算量小,适合大部分实时性要求较高的场合。
本发明提供一种用于信号处理的滤波方法,包括:
步骤1、将采集的模拟量信号转换为数字量信号;
步骤2、获得原始数据数组A=[ai],其中,ai为原始数据组中的第i个元素,i=0,1,2,......,k,k为原始数据组中元素的个数;
步骤3、根据原始数组计算差值数组
Figure BDA0002355258240000021
j为迭代次数,
Figure BDA0002355258240000022
为第j次迭代时差值数组中的第i个元素,设目标迭代次数为d,其中,
Figure BDA0002355258240000023
步骤4、进行迭代计算,当j为奇数时,
Figure BDA0002355258240000024
当j为偶数时,
Figure BDA0002355258240000025
当迭代计算次数达到d次时,退出迭代计算;
步骤5、形成处理后数据组C=[ci],i=0,1,2,......,k,ci为处理后数据组中的第i个元素,
Figure BDA0002355258240000026
步骤6、将处理后数据组的数字量转换为模拟量并输出显示波形图;
进行迭代计算的目的是通过迭代计算弱化数据量之间的差值,对波动进行缓冲;奇数次迭代与偶数次迭代使用不同的计算方法的目的是防止处理后数组超前或滞后,使数据不失真,若统一将当前数值与前一项数值求平均值会造成处理后数组超前,若统一将当前数值与后一项数值求平均值会造成处理后数组滞后,而奇数次迭代将当前数值与后一项数值求平均值,偶数次迭代将当前数值与前一项数值求平均值,可以解决此问题。
优选地,所述目标迭代次数为偶数;目的是当奇数次迭代使数组向前偏移一次后,通过偶数次迭代使数组向后偏移一位,保持数组不超前也不滞后。
从以上技术方案可以看出,本发明具有以下优点:
本发明提供了一种用于信号处理的滤波方法,通过将采集得到的模拟信号转换为数字信号后得到原始数据组,根据原始数据组形成差值数组,将差值数组按计算规则进行多次迭代计算,对数据进行滤波处理。本发明的滤波方法能去除各种干扰造成的波动,使滤波后的数据得到较好的平滑度,并且数据不失真,仍在前后时刻的变化范围内,最大限度接近实际值,而且本发明的方法运算量小,适合大部分实时性要求较高的场合。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1本发明一种用于信号处理的滤波方法的实施流程图;
图2本发明一个实施例中一种用于信号处理的滤波方法滤波前信号原始波形图;
图3本发明一个实施例中一种用于信号处理的滤波方法滤波后信号效果波形图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本实施例提供一种用于信号处理的滤波方法,应用于智能断路器中对电流数据采集后进行滤波处理,实施流程如下:
将智能断路器采集的电流信号模拟量转换成数字量,并保存为原始数据组A=[ai],i=0,1,2,......,2000,原始数据组元素包括从0至2000共2001个数据值;
设置迭代次数d=10。并形成差值数组
Figure BDA0002355258240000031
其中,
Figure BDA0002355258240000041
开始第1次迭代:得到差值数组
Figure BDA0002355258240000042
开始第2次迭代:得到差值数组
Figure BDA0002355258240000043
开始第3次迭代:得到差值数组
Figure BDA0002355258240000044
开始第4次迭代:得到差值数组
Figure BDA0002355258240000045
开始第5次迭代:得到差值数组
Figure BDA0002355258240000046
开始第6次迭代:得到差值数组
Figure BDA0002355258240000047
开始第7次迭代:得到差值数组
Figure BDA0002355258240000048
开始第8次迭代:得到差值数组
Figure BDA0002355258240000049
开始第9次迭代:得到差值数组
Figure BDA00023552582400000410
开始第10次迭代:得到差值数组
Figure BDA00023552582400000411
形成处理后数据组C=[ci],i=0,1,2,......,2000,
Figure BDA0002355258240000051
滤波处理结束。
将处理后数据组中的数字量数据转换成模拟量信号并输出显示波形图。
图2为采集的模拟量信号在滤波前的信号原始波形图,图3为通过本实施例的方法进行滤波后输出的模拟量信号效果波形图。可以看出本发明提供的一种用于信号处理的滤波方法的滤波效果显著,处理后的数据波纹明显减少。本方法能够有效克服偶然因素引起的波动,平滑度较好,而且本方法的运算量小,可以适合大部分实时性要求较强的场合。
本发明也可以用于对智能断路器中采集电压的模拟量信号进行滤波处理,还可用于对电子凸轮示教采集的模拟量数据进行滤波处理。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。

Claims (2)

1.一种用于信号处理的滤波方法,其特征在于,包括:
步骤1、将采集的模拟量信号转换为数字量信号;
步骤2、获得原始数据数组A=[ai],其中,ai为原始数据组中的第i个元素,i=0,1,2,......,k,k为原始数据组中元素的个数;
步骤3、根据原始数据组计算差值数组
Figure FDA0004076327190000011
j为迭代次数,
Figure FDA0004076327190000012
为第j次迭代时差值数组中的第i个元素,设目标迭代次数为d,其中,
Figure FDA0004076327190000013
步骤4、进行迭代计算,当j为奇数时,
Figure FDA0004076327190000014
当j为偶数时,
Figure FDA0004076327190000015
当迭代计算次数达到d次时,退出迭代计算;
步骤5、形成处理后数据组C=[ci],i=0,1,2,......,k,ci为处理后数据组中的第i个元素,
Figure FDA0004076327190000016
步骤6、将处理后数据组的数字量转换为模拟量并输出显示波形图。
2.根据权利要求1所述的一种用于信号处理的滤波方法,其特征在于,所述目标迭代次数d为偶数。
CN202010005848.9A 2020-01-03 2020-01-03 一种用于信号处理的滤波方法 Active CN111181526B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010005848.9A CN111181526B (zh) 2020-01-03 2020-01-03 一种用于信号处理的滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010005848.9A CN111181526B (zh) 2020-01-03 2020-01-03 一种用于信号处理的滤波方法

Publications (2)

Publication Number Publication Date
CN111181526A CN111181526A (zh) 2020-05-19
CN111181526B true CN111181526B (zh) 2023-03-17

Family

ID=70649408

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010005848.9A Active CN111181526B (zh) 2020-01-03 2020-01-03 一种用于信号处理的滤波方法

Country Status (1)

Country Link
CN (1) CN111181526B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114577115B (zh) * 2022-03-29 2023-11-10 广州德亚机械制造有限公司 一种基于plc的高速运动检测边沿位置方法及装置
CN115173555A (zh) * 2022-07-05 2022-10-11 北京虹林泰电气工程有限公司 一种基于低压智能断路器的线路监控保护方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104572076A (zh) * 2014-12-12 2015-04-29 烟台智慧云谷云计算有限公司 一种模拟量采集的软滤波算法及滤波系统
WO2015192696A1 (zh) * 2014-06-20 2015-12-23 深圳大学 分布式网络中基于gossip算法的多目标DOA估计系统及估计方法
CN106603076A (zh) * 2016-11-14 2017-04-26 江苏兆伏爱索新能源有限公司 一种ad采样信号的处理方法及ad采样信号装置
CN110515063A (zh) * 2019-08-13 2019-11-29 青岛海洋科学与技术国家实验室发展中心 基于迭代平稳离散小波变换的水声信号处理方法和装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6221257B2 (ja) * 2013-02-26 2017-11-01 沖電気工業株式会社 信号処理装置、方法及びプログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015192696A1 (zh) * 2014-06-20 2015-12-23 深圳大学 分布式网络中基于gossip算法的多目标DOA估计系统及估计方法
CN104572076A (zh) * 2014-12-12 2015-04-29 烟台智慧云谷云计算有限公司 一种模拟量采集的软滤波算法及滤波系统
CN106603076A (zh) * 2016-11-14 2017-04-26 江苏兆伏爱索新能源有限公司 一种ad采样信号的处理方法及ad采样信号装置
CN110515063A (zh) * 2019-08-13 2019-11-29 青岛海洋科学与技术国家实验室发展中心 基于迭代平稳离散小波变换的水声信号处理方法和装置

Also Published As

Publication number Publication date
CN111181526A (zh) 2020-05-19

Similar Documents

Publication Publication Date Title
CN111181526B (zh) 一种用于信号处理的滤波方法
Namin et al. Efficient hardware implementation of the hyperbolic tangent sigmoid function
CN110688088B (zh) 一种面向神经网络的通用非线性激活函数计算装置和方法
CN109188069B (zh) 一种用于负载开关事件检测的脉冲噪声滤除方法
CN107203135B (zh) 一种基于动态反馈-前馈的自适应控制系统
CN112737711A (zh) 一种自适应噪声基底估计方法及其宽带载波检测方法
CN107565930B (zh) Ad采样的滤波方法及装置
CN112104339B (zh) 一种适用于含跳变边沿信号的平滑滤波方法
CN105281708B (zh) 一种基于分段并行处理的高速fir滤波实现方法
CN108984783A (zh) 一种使用大数据进行递减处理的方法
JPS62274336A (ja) 平方根の概算装置
WO2017012410A1 (zh) 一种信号时延调整方法、装置及计算机存储介质
CN112766076A (zh) 一种电力负荷超短期预测方法、系统、设备及存储介质
CN108805011B (zh) 一种数字滤波方法及系统
CN116384455A (zh) 一种非均匀分段线性化激活函数硬件实现方法
WO2019127357A1 (zh) 一种正交编码器信号滤波方法、装置及系统
CN111062116B (zh) 一种特性未知的模拟量按键信息精确识别方法
CN109100441B (zh) 一种液相色谱曲线去除脉冲干扰的方法
CN109900922B (zh) 转速确定方法、装置、电子设备及可读存储介质
CN107800405B (zh) 基于fpga设计的径向基函数神经网络自适应增强器电路
CN112803428A (zh) 一种受端主网架动态无功电源配置节点选取方法及终端
CN113360852B (zh) 一种道路载荷谱高效清洗方法、系统、装置及存储介质
CN111723677A (zh) 一种基于自适应阈值的小波去噪方法
CN108416311B (zh) 一种基于可编程门阵列和坐标旋转处理的方位角获取方法
CN110620564A (zh) 一种节约fpga资源的抛物线插值滤波器

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant