CN105004913A - 一种电力系统交流采样频率实时跟踪方法 - Google Patents
一种电力系统交流采样频率实时跟踪方法 Download PDFInfo
- Publication number
- CN105004913A CN105004913A CN201510429939.4A CN201510429939A CN105004913A CN 105004913 A CN105004913 A CN 105004913A CN 201510429939 A CN201510429939 A CN 201510429939A CN 105004913 A CN105004913 A CN 105004913A
- Authority
- CN
- China
- Prior art keywords
- analog
- digital converter
- alternating current
- electric power
- 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
Abstract
本发明涉及一种电力系统交流采样频率实时跟踪方法,将被测信号同时输入到固定采样率的第1模数转换器和可变采样率的第2模数转换器进行模数转换;将离散采样序列经过加窗FFT变换,得到幅值谱序列;根据幅值谱序列中的幅值最大值yk及其对应的谱线号k,以及相邻谱线号k+1对应的幅值yk+1对被测交流信号的频率进行校正,得到校正频率f;根据第2模数转换器每周期所需采样点数与校正频率f对第2模数转换器的采样率进行调整,第2模数转换器对电力系统交流信号进行整周期等分同步采样分析;本发明提供的电力系统交流采样频率实时跟踪方法不需使用锁相环,整体可靠性高,跟踪效果好。
Description
技术领域
本发明涉及一种工频交流采样频率实时跟踪技术,属于电力系统交流采样测量技术领域。
背景技术
电力系统中电压、电流信号的采集分析处理目前多数都采用交流采样,即对信号的瞬时值按照一定的规律进行采样,然后进行均方根有效值计算或者通过离散傅里叶变换进行谐波分析,无论是均方根计算还是离散傅里叶变换都需要对信号进行整周期等分采样,才能保证计算的准确性,而电网的频率并不是固定的50Hz,因而需要对信号频率进行实时跟踪。
目前多数的手段是通过硬件方式,通过硬件锁相环电路产生同步于被测交流信号的采样脉冲。这种方式硬件复杂,而且当信号谐波成分较多及含量较大时,需要更复杂的滤波等电路才能达到理想的跟踪效果。
发明内容
本发明的目的在于克服现有技术的不足,并提供一种整体可靠性高,跟踪效果好的电力系统交流采样频率实时跟踪方法。
实现本发明目的所采用的技术方案为,一种电力系统交流采样频率实时跟踪方法,包括如下步骤:
(1)将被测交流信号同时输入到第1模数转换器和第2模数转换器进行模数转换,并输出离散采样序列,其中第1模数转换器为固定采样率模数转换器,其采样率为fs,第2模数转换器为可变采样率模数转换器,其采样率为fc;
(2)取N点长度的第1模数转换器输出的离散采样序列,经过加窗并进行FFT变换,得到幅值谱序列{y0,y1,……,yn},其中n=N-1;
(3)遍历幅值谱序列,获取幅值谱序列中的幅值最大值yk及其对应的谱线号k,以及相邻谱线号k+1对应的幅值yk+1;
(4)根据k、yk和yk+1对被测交流信号的频率进行校正,得到校正频率f,校正公式为f=[k-(v-2)/(v+1)]*fs/N,其中v为校正系数,计算公式为v=yk/yk+1;
(5)调整第2模数转换器的采样率fc,调整公式为fc=f*M,其中M为第2模数转换器每周期所需采样点数,第2模数转换器对电力系统交流信号进行整周期等分同步采样分析。
N和M均为2的整次方。
步骤(2)中加窗使用的窗函数为Hanning窗。
由上述技术方案可知,本发明提供的电力系统交流采样频率实时跟踪方法,通过两个模数转换器实现电力系统交流信号的频率跟踪整周期等分同步采样,其中第1模数转换器的采样率为固定采样率fs,用于对电力系统交流信号进行频率测量校正,可得到接近电力系统交流信号真实频率的校正频率f;第2模数转换器的采样率为可变采样率fc,根据每周期信号所需采样点数M,调整采样率fc=f*M,实现频率跟踪,第2模数转换器输出的采样序列即为电力系统交流信号的整周期等分采样序列。
采用Hanning窗函数对电力系统交流信号进行截断,可达到很高的频率校正精度。
与现有技术相比,本发明提供的方法通过对第1模数转换器采样系列进行频率校正,得到接近电力系统交流信号真实频率的校正频率,校正频率直接乘以第2模数转换器每周期的采样点数,即是第2模数转换器的采样率,满足了第2模数转换器整周期等分采样的要求,因此不需增加额外的硬件锁相电路。该方法利用数字信号处理技术,降低了现有硬件电路的复杂度,节省了硬件设计成本,能够高效可靠地解决电力数据监测中的频率跟踪问题;由于无需额外频率跟踪锁相硬件电路,硬件的简化避免了硬件频率跟踪电路的长期稳定性差、故障隐患多等问题,提高了可靠性,非常适用于以高性能MCU为核心的测量设备。
具体实施方式
下面结合具体实施例对本发明进行详细具体说明,本发明的内容不局限于以下实施例;为便于对本发明的技术特点进行说明,实施例设定一个硬件系统,系统中的硬件包括一个独立外部16bit高精度频率可调整模数转换器(即第2模数转换器,其采样率为fc)与一个片上自带12bit固定频率模数转换器(即第1模数转换器,其采样率为fs)的高性能MCU。
本发明提供的电力系统交流采样频率实时跟踪方法,包括如下步骤:
(1)将被测交流信号同时输入到第1模数转换器和第2模数转换器进行模数转换,其中MCU片上12bit模数转换器采用固定采样频率fs=50*256Hz;
(2)获得N=256个MCU片上12bit模数转换器采样得到的采样点,并加Hanning窗进行快速傅里叶变换,得到幅值谱序列{y0,y1,……,y255};
(3)对得到幅值谱序列进行遍历,获取幅值谱序列中的幅值最大值yk及其对应的谱线号k,以及相邻谱线号k+1对应的幅值yk+1;
(4)MCU根据公式v=yk/yk+1计算校正系数v,随后根据校正公式计算校正频率f,校正公式为f=[k-(v-2)/(v+1)]*fs/N;
(5)MCU根据调整公式fc=f*M调整第2模数转换器的采样率fc,M为第2模数转换器每周期所需采样点数,取M=512,第2模数转换器对电力系统交流信号进行整周期等分同步采样分析。
MCU内部程序对独立外部16bit高精度频率可调整模数转换器获得的频率同步离散采样信号进行有效值、谐波等需要信息的分析与计算,上述实施例获得的频率校正数据如下表:
被测交流信号真实频率(Hz) | 校正频率f(Hz) | 误差精度 |
49.5 | 49.4953 | -0.0005875 |
49.8 | 49.7951 | -0.0006125 |
50 | 49.9982 | -0.000225 |
50.5 | 50.4949 | -0.0006375 |
50.8 | 50.7981 | -0.0002375 |
51 | 50.9999 | -0.0001 |
51.5 | 51.5028 | 0.00035 |
52 | 52.0037 | 0.0004625 |
由上表可得,本发明提供的电力系统交流采样频率实时跟踪方法,其同步采样频率误差精度均在0.001Hz范围内,满足了电力系统交流信号高精度采样的需求。
对于工频电压、电流信号的测量,精度要求较高时,一般都采用外部模数转换器(ADC),通过FPGA、MCU等控制ADC的采样率,而目前MCU的性能越来越高,成本越来越低,很多MCU内部还集成有10位以上的ADC。采用本发明方法可利用MCU内部ADC对信号进行频率测量校正,进而根据真实频率调整外部ADC的采样率,实现频率跟踪,均方根运算和傅里叶分析都采用外部ADC的采样数据。利用本发明提供的方法,可在不增加额外硬件电路的基础上,仅通过软件计算处理解决电网交流信号频率的同步跟踪问题。
Claims (3)
1.一种电力系统交流采样频率实时跟踪方法,其特征在于,包括如下步骤:
(1)将被测交流信号同时输入到第1模数转换器和第2模数转换器进行模数转换,并输出离散采样序列,其中第1模数转换器为固定采样率模数转换器,其采样率为fs,第2模数转换器为可变采样率模数转换器,其采样率为fc;
(2)取N点长度的第1模数转换器输出的离散采样序列,经过加窗并进行FFT变换,得到幅值谱序列{y0,y1,……,yn},其中n=N-1;
(3)遍历幅值谱序列,获取幅值谱序列中的幅值最大值yk及其对应的谱线号k,以及相邻谱线号k+1对应的幅值yk+1;
(4)根据k、yk和yk+1对被测交流信号的频率进行校正,得到校正频率f,校正公式为f=[k-(v-2)/(v+1)]*fs/N,其中v为校正系数,计算公式为v=yk/yk+1;
(5)调整第2模数转换器的采样率fc,调整公式为fc=f*M,其中M为第2模数转换器每周期所需采样点数,第2模数转换器对电力系统交流信号进行整周期等分同步采样分析。
2.根据权利要求1所述的电力系统交流采样频率实时跟踪方法,其特征在于:N和M均为2的整次方。
3.根据权利要求1所述的电力系统交流采样频率实时跟踪方法,其特征在于:步骤(2)中加窗使用的窗函数为Hanning窗。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510429939.4A CN105004913A (zh) | 2015-07-21 | 2015-07-21 | 一种电力系统交流采样频率实时跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510429939.4A CN105004913A (zh) | 2015-07-21 | 2015-07-21 | 一种电力系统交流采样频率实时跟踪方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105004913A true CN105004913A (zh) | 2015-10-28 |
Family
ID=54377656
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510429939.4A Pending CN105004913A (zh) | 2015-07-21 | 2015-07-21 | 一种电力系统交流采样频率实时跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105004913A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110031675A (zh) * | 2019-04-19 | 2019-07-19 | 南京大学 | 一种数采系统实际采样频率的测量方法 |
CN110702972A (zh) * | 2019-10-08 | 2020-01-17 | 浙江大学 | 一种模拟信号自适应采样方法及装置 |
CN110716167A (zh) * | 2019-11-12 | 2020-01-21 | 中电科仪器仪表有限公司 | 任意波形发生器幅度扫频自适应分辨率校准方法及系统 |
CN113466552A (zh) * | 2021-07-14 | 2021-10-01 | 南京海兴电网技术有限公司 | 一种定间隔采样下频率跟踪方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2091462U (zh) * | 1991-05-05 | 1991-12-25 | 中国矿业大学 | 便携式谐波分析仪 |
CN102128982A (zh) * | 2010-12-23 | 2011-07-20 | 东莞市创锐电子技术有限公司 | 一种基于加窗插值fft基频跟踪技术的谐波分析方法 |
CN102135567A (zh) * | 2011-02-25 | 2011-07-27 | 深圳市业通达实业有限公司 | 一种电力系统交流采样实时频率跟踪及谐波测量方法 |
EP1791314A3 (en) * | 2005-11-29 | 2012-10-03 | Samsung Electronics Co., Ltd. | Apparatus and method for carrier frequency synchronization in an OFDM system |
CN103454495A (zh) * | 2013-09-13 | 2013-12-18 | 电子科技大学 | 自适应高精度快速频谱分析方法 |
CN103969508A (zh) * | 2014-05-07 | 2014-08-06 | 江苏天浩达科技有限公司 | 一种实时高精密的电力谐波分析方法及装置 |
-
2015
- 2015-07-21 CN CN201510429939.4A patent/CN105004913A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2091462U (zh) * | 1991-05-05 | 1991-12-25 | 中国矿业大学 | 便携式谐波分析仪 |
EP1791314A3 (en) * | 2005-11-29 | 2012-10-03 | Samsung Electronics Co., Ltd. | Apparatus and method for carrier frequency synchronization in an OFDM system |
CN102128982A (zh) * | 2010-12-23 | 2011-07-20 | 东莞市创锐电子技术有限公司 | 一种基于加窗插值fft基频跟踪技术的谐波分析方法 |
CN102135567A (zh) * | 2011-02-25 | 2011-07-27 | 深圳市业通达实业有限公司 | 一种电力系统交流采样实时频率跟踪及谐波测量方法 |
CN103454495A (zh) * | 2013-09-13 | 2013-12-18 | 电子科技大学 | 自适应高精度快速频谱分析方法 |
CN103969508A (zh) * | 2014-05-07 | 2014-08-06 | 江苏天浩达科技有限公司 | 一种实时高精密的电力谐波分析方法及装置 |
Non-Patent Citations (3)
Title |
---|
余佳兵 等: "窗谱校正方法的实用峰值搜寻算法研究", 《振动工程学报》 * |
刘凤新 等: "基于改进DFT变换的高精度实时电网频率跟踪算法", 《电测与仪表》 * |
谢明 等: "离散频谱分析的一种新校正方法", 《重庆大学学报(自然科学版)》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110031675A (zh) * | 2019-04-19 | 2019-07-19 | 南京大学 | 一种数采系统实际采样频率的测量方法 |
CN110702972A (zh) * | 2019-10-08 | 2020-01-17 | 浙江大学 | 一种模拟信号自适应采样方法及装置 |
CN110716167A (zh) * | 2019-11-12 | 2020-01-21 | 中电科仪器仪表有限公司 | 任意波形发生器幅度扫频自适应分辨率校准方法及系统 |
CN110716167B (zh) * | 2019-11-12 | 2021-12-07 | 中电科思仪科技股份有限公司 | 任意波形发生器幅度扫频自适应分辨率校准方法及系统 |
CN113466552A (zh) * | 2021-07-14 | 2021-10-01 | 南京海兴电网技术有限公司 | 一种定间隔采样下频率跟踪方法 |
CN113466552B (zh) * | 2021-07-14 | 2024-02-02 | 南京海兴电网技术有限公司 | 一种定间隔采样下频率跟踪方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101915874B (zh) | 一种基于傅立叶变换的谐波检测方法 | |
CN102288807B (zh) | 一种测量电网电压闪变的方法 | |
CN103869162B (zh) | 一种基于时域准同步的动态信号相量测量方法 | |
CN103257271B (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
CN103543335B (zh) | 一种同步相量测量方法 | |
CN105004913A (zh) | 一种电力系统交流采样频率实时跟踪方法 | |
CN102236048B (zh) | 一种电力系统相量频率测量方法 | |
CN103983849B (zh) | 一种实时高精度的电力谐波分析方法 | |
CN203909138U (zh) | 一种通用三相交流信号同步采样装置 | |
CN102495282A (zh) | 一种交流电压频率的测量方法 | |
CN203287435U (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置 | |
CN103575984A (zh) | 基于凯塞窗双谱线插值fft的谐波分析方法 | |
CN102135567A (zh) | 一种电力系统交流采样实时频率跟踪及谐波测量方法 | |
CN103969508B (zh) | 一种实时高精密的电力谐波分析方法及装置 | |
CN112986775B (zh) | 一种局部放电脉冲相位主动校准方法 | |
CN104181391A (zh) | 数字功率计谐波检测的方法 | |
CN110389312A (zh) | 一种适用于现场pmu测试的校准器相量测量方法 | |
CN109655665A (zh) | 基于布莱克曼窗的全相位傅里叶谐波分析方法 | |
CN104502707A (zh) | 一种基于三次样条插值的电力系统同步相量测量方法 | |
CN101510687A (zh) | 电网中利用窗函数实现多采样率信号的频率转换方法 | |
CN111625769A (zh) | 一种基于拉格朗日插值和三次指数平滑的pmu-scada数据对时与融合方法 | |
CN105137186A (zh) | 一种微机自动准同期装置的同期电压相位差测量法 | |
CN105137241A (zh) | 一种自适应电网频率的电能质量数据采集方法与装置 | |
CN114966194A (zh) | 一种基于adc的低频计量三相电能表 | |
CN108037351A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151028 |
|
RJ01 | Rejection of invention patent application after publication |