CN112688666A - 一种仪器仪表二级均值滤波算法 - Google Patents
一种仪器仪表二级均值滤波算法 Download PDFInfo
- Publication number
- CN112688666A CN112688666A CN202011485247.9A CN202011485247A CN112688666A CN 112688666 A CN112688666 A CN 112688666A CN 202011485247 A CN202011485247 A CN 202011485247A CN 112688666 A CN112688666 A CN 112688666A
- Authority
- CN
- China
- Prior art keywords
- add
- reg
- filtering
- result
- 1order
- 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.)
- Granted
Links
Images
Landscapes
- Complex Calculations (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明公开了一种仪器仪表二级均值滤波算法，采用二级滤波，滤波参数设置为a＝11,b＝4，滤波周期131.072ms；4us采样率，采用2n作为累加次数，不采用除法器。基于FPGA高速AD采集下对采集结果进行滤波，该算法也适用于其他控制器在高速AD采集后进行有的效滤波处理。能保证仪器仪表表头的电压电流回读值的真实性和稳定性，可靠性高、滤波频率能在极宽宽度范围内灵活选择、执行效率高。
Description
技术领域
本发明涉及一种仪器仪表，尤其涉及一种仪器仪表二级均值滤波算法。
背景技术
为了保证仪器仪表表头的电压电流回读值的真实性和稳定性，需要通过一定的处理技术，满足仪器仪表设备中对于回读电压电流等参数采集结果的滤波效果。
现有技术中：
1.有限(FIR)滤波，如图1所示：
该算法在高频率采样下当低通滤波频率越高时滤波阶数增长过大，导致占用FPGA资源过多无法综合和布局。且由于系数为整数时归一化较困难，滤波结果出现衰减；而且浮点运算占用资源较多，执行效率较低(乘法器执行效率较慢)。
2、无限(IIR)滤波：
该算法克服了有限(FIR)滤波的大部分缺点，但是对于FPGA处理器无法实现。
3、均值滤波：
其中除法操作占用时间与资源较多，而且被除数不能超过57位左右(有上限)，导致滤波频率范围不够宽。当采样频率高时，滤波频率不够低。
4、平滑均值滤波，如图2所示：
该算法需要计算效率非常高，流水方式产生均值结果，缺点是占用RAM资源过多，如果滤波阶数过大时，将导致资源占用率过高。
5、二级均值滤波，如图3所示：
该算法为两级均值算法，极大的扩展了滤波范围。加法器采用分时复用方式，两阶求和次数为2a和2b，求均值只需要截取a和b以上的结果即可，省掉了除法器的使用，降低了加法器的位数，采用分时复用的方式让两阶加法计算用一个加法器，进一步降低了资源使用率。缺点是：均值结果不是流水方式输出。
上述现有技术的滤波不能够满足仪器仪表设备中对于回读电压电流等参数采集结果的滤波效果，总结缺点有：
高采样率，低滤波频率占用资源过大；
滤波频率调节不够灵活，滤波范围不够宽。
发明内容
本发明的目的是提供一种仪器仪表二级均值滤波算法。
本发明的目的是通过以下技术方案实现的：
本发明的仪器仪表二级均值滤波算法，采用二级滤波，滤波参数设置为a＝11,b＝4，滤波周期131.072ms；
4us采样率，采用2n作为累加次数，不采用除法器。
由上述本发明提供的技术方案可以看出，本发明实施例提供的仪器仪表二级均值滤波算法，基于FPGA高速AD采集下对采集结果进行滤波，该算法也适用于其他控制器在高速AD采集后进行有的效滤波处理，可靠性高、滤波频率可在极宽宽度范围内灵活选择、执行效率高。
附图说明
图1为现有技术中有限(FIR)滤波算法示意图；
图2为现有技术中平滑均值滤波算法示意图；
图3为现有技术中二级均值滤波算法示意图；
图4为本发明实施例提供的仪器仪表二级均值滤波算法示意图。
具体实施方式
下面将对本发明实施例作进一步地详细描述。本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
本发明的仪器仪表二级均值滤波算法，其较佳的具体实施方式是：
采用二级滤波，滤波参数设置为a＝11,b＝4，滤波周期131.072ms；
4us采样率，采用2n作为累加次数，不采用除法器。
包括步骤：
1)FPGA控制ADC每4us采集电压电流值，作为输入数据；
2)每次采集数据进入滤波器，与1阶累加结果Reg_1order_Add_Sum相加，同时1 阶累加计数器Reg_1order_Add_Num加1；
3)判断1阶累加计数器Reg_1order_Add_Num的第a位是否变为1：
若为0，则回到第1)步；
若为1，则截取1阶累加结果a位之上的数据位Reg_1order_Add_Sum[a+15:a]，作为2阶均值滤波器的输入数据Reg_1order_Add_Avg，同时清除1阶累加计数器 Reg_1order_Add_Num为0与1阶累加结果Reg_1order_Add_Sum为0，同时触发该结果与2阶累加结果Reg_2order_Add_Sum求和运算，2阶累加计数器Reg_2order_Add_Num 加1，判断2阶累加计数器Reg_2order_Add_Num的b位是否变为1：
若为0，则回到第1)步；
若为1，则截取2阶累加结果b位以上的数据位作为均值滤波的一次结果 Reg_2order_Add_Avg，并清除累加器结果Reg_2order_Add_Sum与累加计数器 Reg_2order_Add_Num。
本发明的仪器仪表二级均值滤波算法，能保证仪器仪表表头的电压电流回读值的真实性和稳定性，可靠性高、滤波频率能在极宽宽度范围内灵活选择、执行效率高。该算法是基于FPGA高速AD采集下对采集结果进行滤波，该算法也适用于其他控制器在高速 AD采集后进行有的效滤波处理。
具体实施例，如图4所示：
1、FPGA控制ADC每4us采集电压电流值，作为输入数据。
2、每次采集数据进入滤波器，与1阶累加结果Reg_1order_Add_Sum相加，同时1 阶累加计数器Reg_1order_Add_Num加1。
3、判断1级累加计数器Reg_1order_Add_Num的第a位是否变为1。若为0，则回到第1步。若为1，则截取1级累加结果a位之上的数据位Reg_1order_Add_Sum[a+15:a]，作为2阶均值滤波器的输入数据Reg_1order_Add_Avg，同时清除1级累加计数器 Reg_1order_Add_Num为0与1级累加结果Reg_1order_Add_Sum为0。同时触发该结果与2级累加结果Reg_2order_Add_Sum求和运算，2级累加计数器Reg_2order_Add_Num 加1，判断2级累加计数器Reg_2order_Add_Num的b位是否变为1，若为0，则回到第1 步。若为1，则截取2级累加结果b为以上的数据位作为均值滤波的一次结果 Reg_2order_Add_Avg，并清除累加器结果Reg_2order_Add_Sum与累加计数器 Reg_2order_Add_Num。
本发明采用二级滤波；采用2n作为累加次数，删除了除法器。能真实反映电源回读的电压电流平均值，目前15KW直流电源DH17800系列使用此方案，使得电源回读值成功滤除掉高频波形，非常稳定的反应了电源的回读电压电流平均值。4us采样率，滤波参数设置为a＝11,b＝4(131.072ms滤波周期)。
具体实施中，该帧格式内部改变地址字节格式定义，用于在复杂的多控制器之间传输数据；数据长度定义成1/2/3/4字节长度，单帧提高传输量。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。
Claims (2)
1.一种仪器仪表二级均值滤波算法，其特征在于：
采用二级滤波，滤波参数设置为a＝11,b＝4，滤波周期131.072ms；
4us采样率，采用2n作为累加次数，不采用除法器。
2.根据权利要求1所述的仪器仪表二级均值滤波算法，其特征在于，包括步骤：
1)FPGA控制ADC每4us采集电压电流值，作为输入数据；
2)每次采集数据进入滤波器，与1阶累加结果Reg_1order_Add_Sum相加，同时1阶累加计数器Reg_1order_Add_Num加1；
3)判断1阶累加计数器Reg_1order_Add_Num的第a位是否变为1：
若为0，则回到第1)步；
若为1，则截取1阶累加结果a位之上的数据位Reg_1order_Add_Sum[a+15:a]，作为2阶均值滤波器的输入数据Reg_1order_Add_Avg，同时清除1阶累加计数器Reg_1order_Add_Num为0与1阶累加结果Reg_1order_Add_Sum为0，同时触发该结果与2阶累加结果Reg_2order_Add_Sum求和运算，2阶累加计数器Reg_2order_Add_Num加1，判断2阶累加计数器Reg_2order_Add_Num的b位是否变为1：
若为0，则回到第1)步；
若为1，则截取2阶累加结果b位以上的数据位作为均值滤波的一次结果Reg_2order_Add_Avg，并清除累加器结果Reg_2order_Add_Sum与累加计数器Reg_2order_Add_Num。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011485247.9A CN112688666B (zh) | 2020-12-16 | 2020-12-16 | 一种仪器仪表二级均值滤波算法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011485247.9A CN112688666B (zh) | 2020-12-16 | 2020-12-16 | 一种仪器仪表二级均值滤波算法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112688666A true CN112688666A (zh) | 2021-04-20 |
| CN112688666B CN112688666B (zh) | 2022-08-19 |
Family
ID=75448187
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011485247.9A Active CN112688666B (zh) | 2020-12-16 | 2020-12-16 | 一种仪器仪表二级均值滤波算法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112688666B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6175849B1 (en) * | 1998-02-10 | 2001-01-16 | Lucent Technologies, Inc. | System for digital filtering in a fixed number of clock cycles |
| CN101540595A (zh) * | 2009-04-14 | 2009-09-23 | 清华大学 | 一种多阶数的均值数字滤波器及滤波方法 |
| CN107016217A (zh) * | 2017-04-28 | 2017-08-04 | 茂硕电源科技股份有限公司 | 一种采样滤波电流值的计算方法及其系统 |
-
2020
- 2020-12-16 CN CN202011485247.9A patent/CN112688666B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6175849B1 (en) * | 1998-02-10 | 2001-01-16 | Lucent Technologies, Inc. | System for digital filtering in a fixed number of clock cycles |
| CN101540595A (zh) * | 2009-04-14 | 2009-09-23 | 清华大学 | 一种多阶数的均值数字滤波器及滤波方法 |
| CN107016217A (zh) * | 2017-04-28 | 2017-08-04 | 茂硕电源科技股份有限公司 | 一种采样滤波电流值的计算方法及其系统 |
Non-Patent Citations (2)
| Title |
|---|
| 周翠: "基于FPGA的高斯滤波器的快速实现", 《空军雷达学院学报》 * |
| 杨澜等: "基于FPGA的车辆振动信号提取方法", 《吉林大学学报(信息科学版)》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112688666B (zh) | 2022-08-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101257288B (zh) | 一种可配置参数的有限冲激响应数字滤波器 | |
| CN104020398B (zh) | 一种换流变压器局部放电波形特征提取方法 | |
| CN105913835B (zh) | 一种基于数学形态学的自适应滤波方法 | |
| CN102098004A (zh) | 一种变带宽数字下变频器及实现方法 | |
| CN104467739B (zh) | 一种带宽、中心频点可调的数字滤波器及其实现方法 | |
| CN102035503B (zh) | 滤波装置及级联积分梳状滤波器的滤波方法 | |
| CN108918931A (zh) | 一种负荷分解中功率信号自适应滤波方法 | |
| CN105281708B (zh) | 一种基于分段并行处理的高速fir滤波实现方法 | |
| CN109391246A (zh) | 一种可配置的数字抽取滤波器 | |
| CN102891662B (zh) | 一种通用的速率下变换、上变换装置及方法 | |
| CN112688666B (zh) | 一种仪器仪表二级均值滤波算法 | |
| CN101316102A (zh) | 用于数字基带接收机的滤波器及其设计方法 | |
| US20140005963A1 (en) | Power Quality Diagnosis for Power Conditioning | |
| CN201860303U (zh) | 数字滤波器电路 | |
| CA2393445A1 (en) | Method and apparatus for recursive data filtering | |
| CN106802293A (zh) | 波形峰值检测方法及装置 | |
| CN107102196A (zh) | 一种复杂电力信号的噪声滤除装置及方法 | |
| CN105811919A (zh) | 用于智能变电站数字化保护的低通滤波器设计方法 | |
| Jayaprakasan et al. | Design of CIC based decimation filter structure using FPGA for WiMAX applications | |
| CN111585542A (zh) | 一种功率信号脉冲噪声滤除的方法及滤波器 | |
| CN111896807B (zh) | 基波频率测量方法、测量终端及存储介质 | |
| CN107979355A (zh) | 一种fir滤波器及其滤波方法 | |
| CN104000578A (zh) | 用于心电信号qrs波实时检测的asic芯片 | |
| CN207518557U (zh) | 一种激光陀螺小抖动稳频无延迟带通滤波器 | |
| CN101473535A (zh) | 用于在数字域处理数字输入信号的方法及用于处理数字输入信号的数字滤波器电路 |
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 |