适用于含跳变边沿信号的平滑滤波方法.pdf

《适用于含跳变边沿信号的平滑滤波方法.pdf》由会员分享，可在线阅读，更多相关《适用于含跳变边沿信号的平滑滤波方法.pdf（12页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 202010885672.0 (22)申请日 2020.08.28 (71)申请人 电子科技大学 地址 611731 四川省成都市高新区 （西区） 西源大道2006号 (72)发明人 马春光陈元鄢然周逸 吕洪光罗勇 (74)专利代理机构 电子科技大学专利中心 51203 代理人 邓黎 (51)Int.Cl. H03H 17/02(2006.01) (54)发明名称 一种适用于含跳变边沿信号的平滑滤波方 法 (57)摘要 本发明公开了一种适用于含跳变边沿信号 的平滑滤波方法， 属于。

2、数字滤波系统中对数据波 形平滑处理领域。 本发明通过求取信号一阶差分 的方法对跳变边沿进行判断， 能够实时自动判断 含有跳变边沿的信号的跳变边沿， 从而在进行平 滑处理时能够自动保留跳变边沿， 相比传统的直 接平滑滤波能够有效防止该类信号跳变边沿信 息丢失， 因此相比传统的平滑滤波算法能够适用 的范围更广。 权利要求书1页 说明书5页 附图5页 CN 112104339 A 2020.12.18 CN 112104339 A 1.一种适用于含跳变边沿信号的平滑滤波方法， 包括以下步骤： S1.采集一段包含完整跳变边沿的待平滑处理数据； S2.将待平滑处理数据输入嵌入式微处理器中， 对于新进入队。

3、列的数据， 在进行取平均 值运算前先抓取新进入队列的数据以及此数据的前一个数据， 若新进入队列的数据为第一 个数据， 则该数据的前一个数据取0； S3.计算两个数据的差值的绝对值， 该绝对值即为待平滑处理数据的一阶后向差分数 据的绝对值， 将该绝对值与判断阈值进行比较， 并输出判断信号； S4.根据判断信号决定是否进行平滑滤波处理， 如果判断信号有效则表示当前数据为 跳变边沿状态， 不进行平均运算； 如果判断信号无效， 则判断当前数据为刚脱离跳变边沿状 态， 或者为非刚脱离状态， 对于前一种状态需要进行重新求平均运算， 对于后一种状态则直 接进行求平均运算。 2.如权利要求1所述的一种适用于含。

4、跳变边沿信号的平滑滤波方法， 其特征在于， 步骤 S3中所述判断阈值通过以下步骤得到： (1)求出步骤S1中的待平滑处理数据的一阶后向差分数据的绝对值； (2)再求得待平滑处理数据的二阶后向差分数据的绝对值， 即对待平滑处理数据一阶 后向差分数据再求一次一阶后向差分， 并求取绝对值； (3)求取二阶后向差分数据的绝对值的最大值， 并找出该最大值对应的时刻； 然后求出 该时刻对应的一阶后向差分数据的绝对值， 选取该绝对值为判断阈值。 3.如权利要求2所述的一种适用于含跳变边沿信号的平滑滤波方法， 其特征在于， 步骤 S4中所述判断信号无效时， 通过以下步骤判断是否为刚脱离跳变边沿状态： 在每一轮。

5、取平均运算时监测步骤S3中所述判断信号是否产生， 如果在上一轮取平均运 算时监测到判断信号产生， 而本轮取平均运算时未监测到判断信号产生， 则表示当前数据 为刚脱离跳变边沿状态， 从而进行重新取平均运算； 如果上一轮取平均运算时未监测到判 断信号产生， 而本轮取平均运算时也未监测到判断信号产生， 则表示当前数据为非刚脱离 状态， 从而进行直接取平均运算。 4.如权利要求2或3所述的一种适用于含跳变边沿信号的平滑滤波方法， 其特征在于， 步骤S4中， 对待滤波数据取平均运算的平滑滤波公式为： 其中， N为平滑窗宽， xn-i、 xn分别为第n-i个、 第n个待平滑处理数据， i1 , 2.n， 。

6、Njp为跳变边沿开始的时刻。 权利要求书 1/1 页 2 CN 112104339 A 2 一种适用于含跳变边沿信号的平滑滤波方法 技术领域 0001 本发明属于数字滤波系统中对数据波形平滑处理领域， 特别涉及一种适用于含跳 变边沿信号的平滑滤波方法。 背景技术 0002 在数字滤波系统中常会出现经一次低通滤波后所得波形仍含有小幅度高频噪声 从而导致波形整体不平整， 会产生小幅度抖动的现象， 这种干扰现象对于需要尽可能保持 平稳的波形来说是无法容忍的， 而由于这些经过一次低通滤波仍然无法滤除的小幅度高频 噪声其频率处在更低的范围， 因此除非进一步减少低通滤波器的通带频率， 否则无法进行 滤除，。

7、 但对于数字滤波器尤其是最常用的FIR(Finite Impulse Response， 有限脉冲响应) 数字滤波器来说， 在保持信号采样频率不变的情况下， 如果要进一步减小其通带频率则意 味着要增加滤波器阶数， 而数字滤波器阶数决定了实现数字滤波器所需要的硬件资源(乘 法器、 加法器等)的多少。 因此为了尽可能节省硬件资源， 并且能够滤除信号本身所含有的 小幅高频噪声， 平滑滤波算法就是最常用最有效的一种手段。 0003 传统的平滑滤波算法的做法是设定合适长度的平滑窗宽N， 并直接将采集到的N个 待平滑处理数据进行平均， 得到的平均值就是这一次的平滑滤波结果， 比如由西安理工大 学李强等人提。

8、出的一种基于双窗宽滑动滤波的ip-iq谐波检测算法中其滑动滤波算法部分 就是采用此算法李强， 王博强， 王子豪.一种基于双窗宽滑动滤波的ip-iq谐波检测算法， 电 力电容器与无功补偿， 40(3):23-28,2019.， 但是这种算法有一个缺点， 即对于一些含有跳 变边沿的信号， 如方波、 三角波、 锯齿波或者脉冲信号等， 使用这种传统平滑滤波方法将会 导致在平滑处理的过程中将这些信号具有的跳变边沿进行平滑处理， 从而致使其跳变边沿 严重失真， 而这种情况对于需要特意保留其跳变边沿的信号(如采样时钟信号、 触发信号、 异步复位信号等)来说是致命的， 因此传统的平滑滤波算法对于此类信号是不适。

9、用的。 发明内容 0004 鉴于上述现有技术中平滑滤波算法无法应对具有跳变边沿信号的不足， 本发明的 目的在于改良平滑滤波算法， 使其能够实时分辨待处理信号是否处于跳变边沿状态， 从而 自动选择是否进行平滑处理。 本发明提出的平滑滤波算法能够适用于含有跳变边沿信号的 平滑处理。 0005 为达到上述目的， 本发明采用的技术方案如下： 0006 一种适用于含跳变边沿信号的平滑滤波方法， 包括以下步骤： 0007 S1.采集一段包含完整跳变边沿(包含上升沿和下降沿)的待平滑处理数据； 0008 S2.将待平滑处理数据输入嵌入式微处理器中， 对于新进入队列的数据， 在进行取 平均值运算前先抓取新进入。

10、队列的数据以及此数据的前一个数据， 若新进入队列的数据为 第一个数据， 则该数据的前一个数据取0； 0009 S3.计算两个数据差值的绝对值， 该绝对值即为待平滑处理数据的一阶后向差分 说明书 1/5 页 3 CN 112104339 A 3 数据的绝对值， 将该绝对值与判断阈值进行比较， 并输出判断信号； 0010 S4.根据判断信号决定是否进行平滑滤波处理， 如果判断信号有效则表示当前数 据为跳变边沿状态， 不进行平均运算； 如果判断信号无效， 则判断当前数据为刚脱离跳变边 沿状态， 或者为非刚脱离状态， 对于前一种状态需要进行重新求平均运算， 对于后一种状态 则直接进行求平均运算。 00。

11、11 进一步地， 步骤S3中所述判断阈值通过以下步骤得到： 0012 (1)求出步骤S1中的待平滑处理数据的一阶后向差分数据的绝对值， 即将原数据 的每一个数据减去前一个数据再取绝对值， 若为第一个数据， 则该数据的前一个数据取0； 0013 (2)再求得待平滑处理数据的二阶后向差分数据的绝对值， 即对待平滑处理数据 一阶后向差分数据再求一次一阶后向差分， 并求取绝对值； 0014 (3)求取二阶后向差分数据的绝对值的最大值， 并找出该最大值对应的时刻； 然后 求出该时刻对应的一阶后向差分数据的绝对值， 选取该绝对值为判断阈值。 0015 进一步地， 步骤S4中所述判断信号无效时， 通过以下步。

12、骤判断是否为刚脱离跳变 边沿状态： 0016 在每一轮取平均运算时监测S3中所述判断信号是否产生， 如果在上一轮取平均运 算时监测到判断信号产生， 而本轮取平均运算时未监测到判断信号产生， 则表示当前数据 为刚脱离跳变边沿状态， 从而进行重新取平均运算； 如果上一轮取平均运算时未监测到判 断信号产生， 而本轮取平均运算时也未监测到判断信号产生， 则表示当前数据为非刚脱离 状态， 从而进行直接取平均运算。 0017 本发明提出的一种适用于含跳变边沿信号的平滑滤波方法， 能够实时自动判断含 有跳变边沿信号的跳变边沿， 从而在进行平滑处理时能够自动保留跳变边沿， 相比传统的 直接平滑滤波能够有效防止。

13、该类信号跳变边沿信息的丢失， 因此本专利提出的平滑滤波算 法适用范围更广。 附图说明 0018 图1是本发明实施例的流程图。 0019 图2是本发明实施例中待平滑处理信号波形图。 0020 图3是本发明实施例中基于一阶差分的阈值选取示意图。 0021 图4是本发明实施例中FPGA中MAC FIR滤波器结构图。 0022 图5是本发明实施例中改进后的平滑滤波器模块框图。 0023 图6是本发明实施例中改进的平滑滤波器工作时序图； 其中图6(a)是信号跳变标 志jump_edge产生时序图， 图6(b)是平滑滤波运算时序图。 0024 图7是本发明实施例中平滑滤波结果波形图； 其中图7(a)是传统。

14、平滑滤波处理结 果图， 图7(b)是改进后平滑滤波处理结果图。 具体实施方式 0025 本发明适用于跳变边沿信号的平滑滤波方法可在不同的可编程器件上设计实现， 如单片机、 PLC、 FPGA、 ARM等。 本实施例以较为主流的FPGA进行硬件实现， 下面结结合附图和 具体实施例对本发明作进一步的详细描述。 说明书 2/5 页 4 CN 112104339 A 4 0026 本实施例所基于的硬件实施平台为FPGA(Field Programmable Gate Array， 可编 程逻辑门阵列)， 由于其具有集成度高、 编译更加灵活， 并行运算处理更加高效等特点， 其从 硬件电路设计层面可以简化。

15、电路复杂程度， 去除电路中冗余逻辑控制电路元件， 从而能够 降低电路设计的开发困难及缩短开发周期， 并且FPGA内部自带功能强大的数字信号处理IP 核以及存储器IP核， 在设计数字滤波器、 实现数字滤波算法上有非常大的优势， 因此使用 FPGA作为数字滤波器实现的硬件平台， 可以不需改变硬件环境就可重新配置参数， 且具有 廉价、 灵活、 方便等优势。 0027 如附图1所示， 本发明实施例提供了一种适用于含跳变边沿信号的平滑滤波方法， 其中对于待平滑处理信号的分析以及判断阈值计算选取部分在MATLAB软件中进行， 算法的 具体硬件实现在FPGA中进行。 0028 首先采集一段包含完整跳变边沿(。

16、包含上升沿和下降沿)的待平滑处理数据， 将其 导入能够进行复杂、 精确计算的MATLAB中， 在MATLAB中作出该段信号波形图。 为能够更直观 说明本实施例， 选取的待平滑处理数据波形如附图2所示， 可以看到该段信号具有上升沿与 下降沿， 在平滑处理中需要保留防止其失真， 同时需要说明的是该段信号已经经过一次FIR 低通滤波， 其中仍然包含小幅高频噪声导致信号波形整体起伏不稳定。 0029 接下来针对传统直接平滑滤波算法进行改进说明， 传统的直接平滑滤波算法的思 路是设定合适长度的平滑窗宽N， 并将依次采集到的N个待平滑处理数据看为一个长度为N 的队列， 该队列中数据随着新数据的输入不断进行。

17、右移操作， 即将队首数据更新为新输入 数据xn， 同时将队列中其他数据整体向右移动， 并将原先处于队尾的数据xn-N去除， 并 且每次新输入一次数据， 在移位操作的同时将该队列所有数据求平均值， 即作为该次平滑 处理后数据yn， 传统的直接平滑滤波算法表达式如式1所示： 0030 0031 其中xn-i表示第n-i个待平滑处理数据,i1,2.n。 0032 这种直接平滑滤波算法如果应用在含跳变边沿的信号中则会导致在平均计算的 过程中将信号的跳变边沿也进行平均处理， 从而丢失跳变边沿信息。 而本发明改进的平滑 滤波算法具有能够判断原信号跳变边沿的能力， 从而对跳变边沿进行保留处理。 0033 本。

18、发明通过求取信号一阶差分的方法对跳变边沿进行判断， 具体方案为： 求出该 段信号所对应的一阶后向差分数据的绝对值1st_differ， 即将原信号的每一个数据减去前 一位数据再取绝对值。 因为对于离散信号来说， 求其一阶差分信号即相当于求连续信号的 一阶导数， 一阶差分信号代表了原离散信号的前后数据变化率， 一阶差分信号越大， 则代表 其原信号前后数据变化越大， 则对于一个整体变化幅度不大但在特定时刻具有跳变边沿的 信号来说， 如果其一阶差分信号在某一时刻陡然增大则说明其跳变边沿发生在这一时刻； 然而对于一般的离散信号来说， 其跳变边沿通常会持续一段时间， 即其跳变边沿为一些离 散数据组成的，。

19、 这就意味在其跳变边沿处原信号的一阶差分信号为一个幅度很大的波峰， 因此如果要根据一阶差分信号决定判断阈值， 需要令该阈值能够区分出此波峰， 从而才能 将跳变边沿完整保留。 为达到此目的， 则需要再求原信号的二阶后向差分信号的绝对值 2nd_differ， 即对原信号的一阶后向差分数据再求一次一阶后向差分， 并求取绝对值。 二阶 差分信号代表了一阶差分信号前后数据的变化率， 因此二阶差分信号的最大值出现的时刻 说明书 3/5 页 5 CN 112104339 A 5 就代表了一阶差分信号的波峰出现时刻， 因此选取该时刻的一阶差分信号的值为判断阈值 thershold， 就可以合理地区分出原信号。

20、的整个跳变边沿， 判断阈值的表达式为式2、 式3、 式 4所示： 0034 xn|xn-xn-1| (式2) 0035 max(|2x|)|2xthershold_N| (式3) 0036 thershold|xthershold_N| (式4) 0037 其中xn表示第n个待滤波数据的一阶后向差分数据的绝对值； 2x表示待滤 波数据的二阶后向差分数据， max(.)表示求最大值函数， thershold_N为判断阈值 thershold所在位置。 0038 以此判断方法计算得到的判断阈值如附图3所示， 可以看到该选定的判断阈值较 为完整地区分了一阶差分信号的最大波峰， 对应原信号则为完整的跳。

21、变边沿。 0039 通过判断阈值判断原信号的跳变边沿是否到来， 从而决定是否对数据进行平均处 理， 因此步骤S4本质上其实为分段进行平滑处理。 对于原待平滑处理信号， 可以根据其跳变 边沿到来的时刻Njp对其进行分段， 如果只将原信号分为未发生跳变的时刻， 即nNjp； 与发 生跳变的时刻， 即nNjp， 则会导致在对刚脱离跳变时刻的数据进行平滑处理时， 仍然将前 面包含跳变边沿的数据放入队列中进行平均从而导致失真， 因此对于刚脱离跳变时刻而仍 在一个平滑滤波窗宽N的时刻内， 即NjpnNjp+N， 其中的数据在进行平滑处理时必须清除队 列中之前包含跳变边沿数据的旧数据， 从Njp之后的数据开。

22、始重新进行平均运算， 同时由于 此时相当于该队列中非0有效数据从1开始不断增加到N， 因此此时平滑滤波器系数也不应 保持1/N不变， 而应该也从1不断减小到1/N， 因此得出改进后的平滑滤波算法的公式为： 0040 0041 FPGA负责实现平滑滤波系统硬件的实现， 由式1可以得知滑动平均滤波算法可以 看作一种特殊的FIR滤波算法， 即该FIR滤波器的系数全部为1/N， 因此可以在传统FIR滤波 器结构的基础上通过改进来实现本发明中的改进平滑滤波算法。 在FPGA中最常用的一种结 构为MAC结构， 该结构经典框图如附图4所示， 该结构采用两片RAM来分别控制采样数据X以 及滤波器系数H的读取，。

23、 并通过一个乘累加器对两路数据进行乘累加运算， 从而得到滤波后 信号数据Y。 本发明即在这种结构的基础上进行改进， 来实现上述式5所示改进后的平滑滤 波算法。 0042 本发明实施例中FPGA基于改进平滑滤波方法的模块框图如附图5所示， 为了在 FPGA中计算每次待平滑处理数据xn的一阶后向差分的绝对值需要在每次卷积的特定时 刻采集xn与xn-1的数据， 并额外使用一个浮点数减法器进行实时相减运算， 将得出数 据的绝对值x_sub_abs与设定好的判断阈值thershold进行比较， 如果大于则发出跳变标志 信号jump_edge， 此信号有效状态下不进行平滑处理， 如果此轮数据处理jump_。

24、edge信号刚 好无效， 则说明数据刚脱离跳变边沿， 则需要对存储队列数据的x_ram发出复位指令清空输 出端口的旧数据； 同时重新加载存放平均系数的H_average_rom中的数据， 其具体工作时序 图如附图6(a)、 6(b)所示； 其中图6(a)中， wrong表示错误数据， jump_edge不会根据此数据 说明书 4/5 页 6 CN 112104339 A 6 选择是否有效； 待平滑处理数据的一阶后向差分的绝对值x_subn的表达式为： 0043 x_subn|xn-xn-1| (式6) 0044 从附图5中的乘累加器模块输出的数据y就是经过该平滑滤波器平滑处理后的信 号数据。 。

25、0045 通过附图7(a)、 7(b)波形结果可以看出， 采用传统的平滑滤波算法处理该含跳变 边沿的信号后， 原信号所包含的跳变边沿信息基本完全丢失， 而采用本发明改进的平滑滤 波算法对该信号进行平滑处理后， 既能去除原信号中包含的小幅度高频噪声， 增强原信号 的稳定性， 同时也可以保留原信号的跳变边沿信息， 因此本发明设计的平滑滤波算法能够 适用于含跳变边沿信号的平滑处理。 0046 由于本发明所采用的硬件实施平台为FPGA， 因此本领域技术人员在采用本发明时 可以方便地根据实际情况调整该平滑滤波器的窗宽N而不需对本发明其他软硬件结构进行 改动， 同时本发明可以在FPGA中实现与前一级FIR滤波器的级联， 从而实现实时的FIR滤波 与平滑滤波处理。 说明书 5/5 页 7 CN 112104339 A 7 图1 说明书附图 1/5 页 8 CN 112104339 A 8 图2 图3 说明书附图 2/5 页 9 CN 112104339 A 9 图4 图5 图6(a) 说明书附图 3/5 页 10 CN 112104339 A 10 图6(b) 图7(a) 说明书附图 4/5 页 11 CN 112104339 A 11 图7(b) 说明书附图 5/5 页 12 CN 112104339 A 12 。