CN102435849A - 使用数字滤波器提高频率测量精度的方法 - Google Patents
使用数字滤波器提高频率测量精度的方法 Download PDFInfo
- Publication number
- CN102435849A CN102435849A CN2011103297340A CN201110329734A CN102435849A CN 102435849 A CN102435849 A CN 102435849A CN 2011103297340 A CN2011103297340 A CN 2011103297340A CN 201110329734 A CN201110329734 A CN 201110329734A CN 102435849 A CN102435849 A CN 102435849A
- Authority
- CN
- China
- Prior art keywords
- signal
- filter
- frequency
- digital filter
- sampling
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000005070 sampling Methods 0.000 claims abstract description 34
- 238000001914 filtration Methods 0.000 claims abstract description 15
- 230000004044 response Effects 0.000 claims abstract description 6
- 238000013461 design Methods 0.000 claims description 4
- 230000035939 shock Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
Images
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Measurement Of Current Or Voltage (AREA)
Abstract
一种使用数字滤波器提高频率测量精度的方法，外部信号依次通过抗混叠高频滤波器、高速信号采样模块传送至数字滤波器，该滤波器可以用差分方程表示：其中y(n)为滤波器的输出，x(n)为滤波器的输入，ai为滤波器系数，M为数字滤波器的阶数，滤波器系数即为滤波器的冲击响应系列h(n)，带有干扰和谐波的信号采样值x(n)，经过与数字滤波器h(n)卷积的结果就为我们希望的滤波之后的信号y(n)；经过数字滤波器滤波的信号y(n)与原始信号x(n)之间存在相位差，精确的相位差可以计算出来，但是不影响信号的周期，因此测量y(n)的频率也就是x(n)的真实频率，保证频率测量的高精度。
Description
技术领域
本发明涉及交流采样技术，尤其是使用数字滤波器提高频率测量精度的方法。
背景技术
随着交流采样技术的广泛应用，频率测量的精度显得越来越重要，现在一般在采样信号中查找相邻的两个过零点，用相邻过零点之间的时间计算工频信号的频率，对于频率跟踪算法，再用这个测量频率，修改采样信号的频率，以保证每一个周波的采样点一样；而对于使用准同步算法的设备，也要以这个频率为依据，对采样数据进行插值计算。可见频率的测量精度直接关系到交流采样设备的测量准确度。
在产品研发中，当信号中包含了多种谐波成分，并且谐波含量很高的时候，经常遇到频率测量不正确的问题，从而直接导致测量数据的错误，或者准确度达不到要求，导致这种结果的原因是因为在过零点附近密集分布多个过零点，干扰信号淹没了真实信号的过零点，从而产生过零点搜索错误。
发明内容
本发明所要解决的技术问题是提供一种使用数字滤波器提高频率测量精度的方法，即使信号中含有很多高频干扰或谐波成分，滤波后，仍然可以很准确的找到所需要的过零点，保证频率测量的高精度。
为解决上述技术问题，本发明的技术方案是：一种使用数字滤波器提高频率测量精度的方法，外部信号源依次通过抗混叠高频滤波器、高速信号采样模块传送至数字滤波器；所述数字滤波器为FIR滤波器，用差分方程表示为：
其中y(n)为数字滤波器的输出，x(n)为数字滤波器的输入，ai为数字滤波器的系数，M为数字滤波器的阶数；高速信号采样模块实现AD采样；具体的测量方法如下：
(1)外部信号经过电压互感器、电流互感器后，将一次高压、大电流信号转换为500V以下的交流电压信号，小于100A的交流电流信号，再经过一级精密测量电压互感器、测量电流互感器将信号变换为后面测量电路可以测量范围的小信号；
(2)外部信号经过处理后，还含有很多超出AD采样速率的频率成分，这些信号会严重干扰采样信号的采样，会使采样信号失真，因此需要抗混叠滤波器滤掉超出采样率2倍频率的信号；
(3)经过抗混叠滤波器滤波后的信号可以被AD采样，采样信号的频率与AD的采样率满足采样定理的要求，经过AD采样后的信号为数字信号，用符号x(n)表示，n取值为整数0、1…n-1，代表有n个采样点；
(4)根据我们所关心的信号频率，提出截止频率的频率值，对截止频率外的信号的衰减程度，算成理想滤波器的阶数和系数，再根据实际设备情况将理想滤波器的阶数和系数进行修正，就得到了实际使用的滤波器的阶数和系数；滤波器系数即为滤波器的冲击响应系列h(n)，带有干扰和谐波的信号采样值x(n)，经过与数字滤波器h(n)卷积的结果就为我们希望的滤波之后的信号y(n)；经过数字滤波器滤波的信号y(n)与原始信号x(n)之间存在相位差，精确的相位差可以计算出来，但是不影响信号的周期，因此测量y(n)的频率也就是x(n)的真实频率。
作为改进，所述步骤(4)中，根据50Hz左右的工频信号，设计数字滤波器，在-3dB增益点的截止频率为90Hz。
本发明与现有技术相比所带来的有益效果是：
设计一个滤波器，可以用差分方程表示：
其中y(n)为滤波器的输出，x(n)为滤波器的输入，ai为滤波器系数，M为数字滤波器的阶数，滤波器系数即为滤波器的冲击响应系列h(n)，带有干扰和谐波的信号采样值x(n)，经过与数字滤波器h(n)卷积的结果就为我们希望的滤波之后的信号y(n)；经过数字滤波器滤波的信号y(n)与原始信号x(n)之间存在相位差，精确的相位差可以计算出来，但是不影响信号的周期，因此测量y(n)的频率也就是x(n)的真实频率。
附图说明
图1为本发明硬件连接图。
具体实施方式
下面结合说明书附图对本发明作进一步说明。
如图1所示，一种使用数字滤波器提高频率测量精度的方法，外部信号源依次通过抗混叠高频滤波器、高速信号采样模块传送至数字滤波器；所述数字滤波器为FIR滤波器，用差分方程表示为：
其中y(n)为数字滤波器的输出，x(n)为数字滤波器的输入，ai为数字滤波器的系数，M为数字滤波器的阶数；高速信号采样模块实现AD采样；具体的测量方法如下：
(1)外部信号经过电压互感器、电流互感器后，将一次高压、大电流信号转换为500V以下的交流电压信号，小于100A的交流电流信号，再经过一级精密测量电压互感器、测量电流互感器将信号变换为后面测量电路可以测量范围的小信号；
(2)外部信号经过处理后，还含有很多超出AD采样速率的频率成分，这些信号会严重干扰采样信号的采样，会使采样信号失真，因此需要抗混叠滤波器滤掉超出采样率2倍频率的信号；
(3)经过抗混叠滤波器滤波后的信号可以被AD采样，采样信号的频率与AD的采样率满足采样定理的要求，经过AD采样后的信号为数字信号，用符号x(n)表示，n取值为整数0、1…n-1，代表有n个采样点；
(4)根据我们所关心的信号频率，提出截止频率的频率值，对截止频率外的信号的衰减程度，算成理想滤波器的阶数和系数，再根据实际设备情况将理想滤波器的阶数和系数进行修正，就得到了实际使用的滤波器的阶数和系数；滤波器系数即为滤波器的冲击响应系列h(n)，带有干扰和谐波的信号采样值x(n)，经过与数字滤波器h(n)卷积的结果就为我们希望的滤波之后的信号y(n)；经过数字滤波器滤波的信号y(n)与原始信号x(n)之间存在相位差，精确的相位差可以计算出来，但是不影响信号的周期，因此测量y(n)的频率也就是x(n)的真实频率。
用软件实现数字滤波器有很高的灵活性，可以设计硬件不可能实现的近似理想特性的各种滤波器，通过软件改变滤波器参数即可实现不同性能的滤波器。在实际应用中，我们一般只关心某一种频率，如：50Hz左右的工频信号，设计一个软件滤波器，在-3dB增益点的截止频率为90Hz，将采样值信号进行低通滤波后，滤除90Hz以上的频率信号，滤波后的信号就是我们最关注的工频信号，对滤除工频信号以外的信号查找相近的过零点，比使用滤波前的信号更准确，即使信号中含有很多高频干扰或谐波成分，滤波后，仍然可以很准确的找到所需要的过零点，保证频率测量的高精度。
用软件实现本实施例数字滤波器，用C语言描述如下：
假设滤波系数用变量a[m]表示，输入信号用x[n]表示，输出用y[n]表述，滤波器的阶数为m，则实现信号滤波的软件代码可以这么写：
Claims (2)
1.一种使用数字滤波器提高频率测量精度的方法，其特征在于：外部信号源依次通过抗混叠高频滤波器、高速信号采样模块传送至数字滤波器；所述数字滤波器为FIR滤波器，用差分方程表示为：
其中y(n)为数字滤波器的输出，x(n)为数字滤波器的输入，ai为数字滤波器的系数，M为数字滤波器的阶数；高速信号采样模块实现AD采样；具体的测量方法如下：
(1)外部信号经过电压互感器、电流互感器后，将一次高压、大电流信号转换为500V以下的交流电压信号，小于100A的交流电流信号，再经过一级精密测量电压互感器、测量电流互感器将信号变换为后面测量电路可以测量范围的小信号；
(2)外部信号经过处理后，还含有很多超出AD采样速率的频率成分，这些信号会严重干扰采样信号的采样，会使采样信号失真，因此需要抗混叠滤波器滤掉超出采样率2倍频率的信号；
(3)经过抗混叠滤波器滤波后的信号可以被AD采样，采样信号的频率与AD的采样率满足采样定理的要求，经过AD采样后的信号为数字信号，用符号x(n)表示，n取值为整数0、1…n-1，代表有n个采样点；
(4)根据我们所关心的信号频率，提出截止频率的频率值，对截止频率外的信号的衰减程度，算成理想滤波器的阶数和系数，再根据实际设备情况将理想滤波器的阶数和系数进行修正，就得到了实际使用的滤波器的阶数和系数；滤波器系数即为滤波器的冲击响应系列h(n)，带有干扰和谐波的信号采样值x(n)，经过与数字滤波器h(n)卷积的结果就为我们希望的滤波之后的信号y(n)；经过数字滤波器滤波的信号y(n)与原始信号x(n)之间存在相位差，精确的相位差可以计算出来，但是不影响信号的周期，因此测量y(n)的频率也就是x(n)的真实频率。
2.根据权利要求1所述的使用数字滤波器提高频率测量精度的方法，其特征在于：所述步骤(4)中，根据50Hz左右的工频信号，设计数字滤波器，在-3dB增益点的截止频率为90Hz。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2011103297340A CN102435849A (zh) | 2011-10-26 | 2011-10-26 | 使用数字滤波器提高频率测量精度的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2011103297340A CN102435849A (zh) | 2011-10-26 | 2011-10-26 | 使用数字滤波器提高频率测量精度的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN102435849A true CN102435849A (zh) | 2012-05-02 |
Family
ID=45984004
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2011103297340A Pending CN102435849A (zh) | 2011-10-26 | 2011-10-26 | 使用数字滤波器提高频率测量精度的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102435849A (zh) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104569532A (zh) * | 2015-01-29 | 2015-04-29 | 国家电网公司 | 一种电子式互感器采集合并一体化装置 |
| CN105973349A (zh) * | 2016-04-19 | 2016-09-28 | 温州科技职业学院 | 一种科氏流量计固有频率的在线预测方法 |
| CN106405230A (zh) * | 2016-08-23 | 2017-02-15 | 北京德威特继保自动化科技股份有限公司 | 频率测量方法和装置 |
| CN108026832A (zh) * | 2015-10-07 | 2018-05-11 | 株式会社电装 | 转速检测装置 |
| CN108535541A (zh) * | 2018-03-19 | 2018-09-14 | 湖北工业大学 | 一种带有数字工频滤波的相位比较装置及方法 |
| CN108594003A (zh) * | 2018-06-30 | 2018-09-28 | 宁波中车时代传感技术有限公司 | 一种电压传感器电路 |
| CN109889692A (zh) * | 2019-02-21 | 2019-06-14 | 无锡中感微电子股份有限公司 | 一种电视信号重采样处理方法及装置 |
| CN111049501A (zh) * | 2020-01-14 | 2020-04-21 | 西安清泰科新能源技术有限责任公司 | 一种电压信号采样滤波电路 |
| CN112162144A (zh) * | 2020-08-19 | 2021-01-01 | 国网江西省电力有限公司供电服务管理中心 | 一种基于卷积增益的电网频率快速测量方法及装置 |
| CN112557750A (zh) * | 2020-11-23 | 2021-03-26 | 国网湖南省电力有限公司 | 一种计量装置互感器二次回路信号频率测量方法及系统 |
| CN113466537A (zh) * | 2021-05-26 | 2021-10-01 | 国网江苏省电力有限公司淮安供电分公司 | 一种电力系统抗干扰电流电压采集方法 |
-
2011
- 2011-10-26 CN CN2011103297340A patent/CN102435849A/zh active Pending
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104569532A (zh) * | 2015-01-29 | 2015-04-29 | 国家电网公司 | 一种电子式互感器采集合并一体化装置 |
| CN108026832A (zh) * | 2015-10-07 | 2018-05-11 | 株式会社电装 | 转速检测装置 |
| CN105973349A (zh) * | 2016-04-19 | 2016-09-28 | 温州科技职业学院 | 一种科氏流量计固有频率的在线预测方法 |
| CN105973349B (zh) * | 2016-04-19 | 2019-04-09 | 温州科技职业学院 | 一种科氏流量计固有频率的在线预测方法 |
| CN106405230A (zh) * | 2016-08-23 | 2017-02-15 | 北京德威特继保自动化科技股份有限公司 | 频率测量方法和装置 |
| CN106405230B (zh) * | 2016-08-23 | 2019-04-12 | 北京德威特电气科技股份有限公司 | 频率测量方法和装置 |
| CN108535541A (zh) * | 2018-03-19 | 2018-09-14 | 湖北工业大学 | 一种带有数字工频滤波的相位比较装置及方法 |
| CN108594003A (zh) * | 2018-06-30 | 2018-09-28 | 宁波中车时代传感技术有限公司 | 一种电压传感器电路 |
| CN109889692A (zh) * | 2019-02-21 | 2019-06-14 | 无锡中感微电子股份有限公司 | 一种电视信号重采样处理方法及装置 |
| CN109889692B (zh) * | 2019-02-21 | 2021-05-04 | 无锡中感微电子股份有限公司 | 一种电视信号重采样处理方法及装置 |
| CN111049501A (zh) * | 2020-01-14 | 2020-04-21 | 西安清泰科新能源技术有限责任公司 | 一种电压信号采样滤波电路 |
| CN111049501B (zh) * | 2020-01-14 | 2023-04-25 | 西安清泰科新能源技术有限责任公司 | 一种电压信号采样滤波电路 |
| CN112162144A (zh) * | 2020-08-19 | 2021-01-01 | 国网江西省电力有限公司供电服务管理中心 | 一种基于卷积增益的电网频率快速测量方法及装置 |
| CN112162144B (zh) * | 2020-08-19 | 2023-07-21 | 国网江西省电力有限公司供电服务管理中心 | 一种基于卷积增益的电网频率快速测量方法及装置 |
| CN112557750A (zh) * | 2020-11-23 | 2021-03-26 | 国网湖南省电力有限公司 | 一种计量装置互感器二次回路信号频率测量方法及系统 |
| CN113466537A (zh) * | 2021-05-26 | 2021-10-01 | 国网江苏省电力有限公司淮安供电分公司 | 一种电力系统抗干扰电流电压采集方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102435849A (zh) | 使用数字滤波器提高频率测量精度的方法 | |
| CN101915874B (zh) | 一种基于傅立叶变换的谐波检测方法 | |
| Yao et al. | Measurement of power system harmonic based on adaptive Kaiser self‐convolution window | |
| CN103257271A (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
| CN111679236B (zh) | 一种直流暂态阶跃响应延时测试方法、系统及装置 | |
| CN101539596A (zh) | 电网频率监测的方法 | |
| CN203287435U (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置 | |
| Upadhyay | Class of Recursive Wideband Digital Differentiators and Integrators. | |
| CN105487034A (zh) | 一种0.05级电子式互感器校验方法及系统 | |
| CN103969508B (zh) | 一种实时高精密的电力谐波分析方法及装置 | |
| CN104502707A (zh) | 一种基于三次样条插值的电力系统同步相量测量方法 | |
| Li | Design and realization of FIR digital filters based on MATLAB | |
| CN103575987B (zh) | 基于dsp间谐波检测仪器及其检测方法 | |
| CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
| CN106680585A (zh) | 谐波/间谐波的检测方法 | |
| CN104111373A (zh) | 一种用于智能变电站数字化电能的计量方法 | |
| CN109239463B (zh) | 一种基于线性修正算法的介质损耗测量方法 | |
| CN103592514A (zh) | 一种新型谐波高精度检测方法 | |
| CN117169590B (zh) | 一种基于软件变采样率的电力谐波分析的方法和装置 | |
| CN110082642B (zh) | 基于全相位微分滤波的电网工况故障时刻检测方法及装置 | |
| CN111505375A (zh) | 一种频率测量方法及装置 | |
| CN104504272B (zh) | 基于理想重构滤波器冲激响应拟合的数字积分器实现方法 | |
| CN115825546A (zh) | 一种直流暂态电流宽频数字标准器实现方法 | |
| CN105811919A (zh) | 用于智能变电站数字化保护的低通滤波器设计方法 | |
| CN107104674A (zh) | 电参量计量过程中的噪声功率估计装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20120502 |