CN102692540A - 一种频率跟踪交流采样方法 - Google Patents
一种频率跟踪交流采样方法 Download PDFInfo
- Publication number
- CN102692540A CN102692540A CN2012100739404A CN201210073940A CN102692540A CN 102692540 A CN102692540 A CN 102692540A CN 2012100739404 A CN2012100739404 A CN 2012100739404A CN 201210073940 A CN201210073940 A CN 201210073940A CN 102692540 A CN102692540 A CN 102692540A
- Authority
- CN
- China
- Prior art keywords
- value
- sampling
- zero
- frequency
- current 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
Images
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开一种频率跟踪交流采样方法,(1)过零点的检测:对于正弦型模拟量信号,其值总是围绕零值上下周期性变化,每一个周期将经历一次由负到正和由正到负的变化,两次由负到正的变化点的时间差及为信号的周期,程序中通过2次采样的值由负(或零)到正的改变记为过零点;(2)采样周期的调整(频率跟踪):如果每个周期采样N个点,我们把采样点序号记为0,1,2,3……N-2,N-1,交流采样无处不在,本发明不但可以省去过零检测硬件电路以及系统开销,还进一步提高采样精度。本发明代码简单可靠,占用资源极小。
Description
技术领域
本发明涉及一种频率跟踪交流采样方法。
背景技术
对正弦型模拟信号的分析通常通过交流采样,再对采样数据进行傅里叶变换以及均方根等算法,计算模拟信号的谐波幅值相位和有效值等参数。但这要求交流采样时取一个整周期的数据进行计算,否则产生误差。
正弦交流电在任何1个周期内采用积分方法可求得电压、电流和功率参数[2],且三者计算方法基本相同,所以本文以电压为例,其积分公式为
式中U为电压有效值;Um为电压幅值;ω为角频率;θ0为电压初相角。
如果所求电压包含整数次谐波,则只要在式(1)的积分变量中加入谐波分量。计算结果就将包含谐波值。在实际交流采样计算中用均方根法代替式(1)求电压有效值.均方根法公式为
式中un=Umsin(θ0+nTs);n为每周期采样点数;Ts为采样周期。
上式可变换为
当采样频率不变,信号频率发生变化时,使用式(2)计算会产生误差。若电压频率为f,采样频率为fs,每周期的采样点数为fs/f=N+ m/M,0<m/M<1,不是一个整数,则用均方根法计算电压时,式(2)应该改为
式中un为最后1个采样点电压值;△m/M为un所占采样点数的比例。
在实际计算时,小数△m/M很难参与计算,N+△m/M只能取为整数N或者N+1,若四舍五入取采样点时,则最大会相差0.5个点。将式(4)改写为
假设△m/M<0.5,采样点数取N,则
近似代替U2,相对误差为
于是
由式(7)可以看出,误差和每周期采样点数N、最后1个采样点uN及uN所占采样点数的比例△m/M有关。显然,当|uN|=O或|uN|=2 U,且△m/M=0.5时误差绝对值最大,最大相对误差为
最大误差绝对值和每周期采样点数的关系如图1所示(uN=0时,相对误差为正;uN=2 U时为负)。
当N=128时。最大误差仍达到0.20%:同时由于各种原因每次采样的初相角吼总会发生变化,所以|uN|将在0~2 U之间变化,计算结果的误差将在±0.20%内晃动。如果不使用四舍五入法而使用截断算法计算采样点数,误差将更大。
以往的系统中,其一是假定被测信号的频率不变,如电力系统中就假定交流电的频率是50Hz恒定的,而直接按其整数倍固定采样频率,由于实际频率多少会有些变化以及定时器本身有些误差,这只能用在对精度要求不高的场合;其二采用过零触发电路检测信号周波的开始和结束时间,再通过定时器计算出信号周期和采样间隔。这需要引入硬件电路以及定时器本身带来的误差。
发明内容
为了克服现有技术中存在的上述不足之处,本发明的目的在于提供一种频率跟踪交流采样方法,能够对信号进行频率跟踪的交流采样,简化硬件的设计,提高频率跟踪精度从而提高测量精度。
为了达到上述之目的,发明采用如下具体技术方案:一种频率跟踪交流采样方法,其特征在于以下步骤:
(1)过零点的检测:对于正弦型模拟量信号,其值总是围绕零值上下周期性变化,每一个周期将经历一次由负到正和由正到负的变化,两次由负到正的变化点的时间差及为信号的周期,程序中通过2次采样的值由负(或零)到正的改变记为过零点;
(2)采样周期的调整(频率跟踪):如果每个周期采样N个点,我们把采样点序号记为0,1,2,3……N-2,N-1。
与现有的技术相比,发明具有以下突出优点和效果:交流采样无处不在,本发明不但可以省去过零检测硬件电路以及系统开销,还进一步提高采样精度。本发明代码简单可靠,占用资源极小。
附图说明
图1为现有技术中的最大误差与每周期采样点数的关系图。
具体实施方式
实施例1
一种频率跟踪交流采样方法,其特征在于以下步骤:
(1)过零点的检测:对于正弦型模拟量信号,其值总是围绕零值上下周期性变化,每一个周期将经历一次由负到正和由正到负的变化,两次由负到正的变化点的时间差及为信号的周期,程序中通过2次采样的值由负(或零)到正的改变记为过零点;
(2)采样周期的调整(频率跟踪):如果每个周期采样N个点,我们把采样点序号记为0,1,2,3……N-2,N-1。当采样频率与信号频率同步的时候,每次过零点都会出现在相同的序号位置,如果本次过零点的序号大于上次过零点的序号,则表示采样频率太高,则降低采样频率(加大采样周期,按定时器分辨率进行),反之则表示采样频率太低,则增加采样频率,通过一定时间的调整,采样频率则保证了与信号频率同步。如果信号频率按50Hz,定时器分辨率1uS,每个周波采样64个点。当频率由50Hz改变到49.5Hz时,其调整时间是3个周期。如果定时器分辨率是0.1uS则调整时间是31周期,当然其精度更高。当然为了提高跟踪速度,也可以通过2次过零点的时间差直接计算调整量。不过由于电力系统频率变化非常缓慢,这个就不是必须的。由于采样时间误差可以累积,因此即使只有很小的不同步通过一定的时间积累都会反映出来而得到调整。
Claims (1)
1.一种频率跟踪交流采样方法,其特征在于以下步骤:
(1)过零点的检测:对于正弦型模拟量信号,其值总是围绕零值上下周期性变化,每一个周期将经历一次由负到正和由正到负的变化,两次由负到正的变化点的时间差及为信号的周期,程序中通过2次采样的值由负(或零)到正的改变记为过零点;
(2)采样周期的调整(频率跟踪):如果每个周期采样N个点,我们把采样点序号记为0,1,2,3……N-2,N-1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100739404A CN102692540A (zh) | 2012-03-20 | 2012-03-20 | 一种频率跟踪交流采样方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100739404A CN102692540A (zh) | 2012-03-20 | 2012-03-20 | 一种频率跟踪交流采样方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102692540A true CN102692540A (zh) | 2012-09-26 |
Family
ID=46858132
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012100739404A Pending CN102692540A (zh) | 2012-03-20 | 2012-03-20 | 一种频率跟踪交流采样方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102692540A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103383412A (zh) * | 2013-07-10 | 2013-11-06 | 珠海许继芝电网自动化有限公司 | 一种自适应软硬件频率跟踪采样的方法 |
CN112485511A (zh) * | 2020-12-11 | 2021-03-12 | 国网辽宁省电力有限公司阜新供电公司 | 基于跟踪补偿算法的交流采样装置及其采样补偿方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0915032A (ja) * | 1995-07-03 | 1997-01-17 | Tech Res & Dev Inst Of Japan Def Agency | 信号周波数追尾方法 |
US6298309B1 (en) * | 1999-03-05 | 2001-10-02 | Abb Power T&D Company Inc. | Current based frequency tracking method and apparatus |
CN101350518A (zh) * | 2008-09-05 | 2009-01-21 | 国电南瑞科技股份有限公司 | 数字化保护测控装置采样值处理方法 |
CN101655519A (zh) * | 2009-09-14 | 2010-02-24 | 国电南京自动化股份有限公司 | 数字化变电站测控装置交流采样数据处理方法 |
-
2012
- 2012-03-20 CN CN2012100739404A patent/CN102692540A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0915032A (ja) * | 1995-07-03 | 1997-01-17 | Tech Res & Dev Inst Of Japan Def Agency | 信号周波数追尾方法 |
US6298309B1 (en) * | 1999-03-05 | 2001-10-02 | Abb Power T&D Company Inc. | Current based frequency tracking method and apparatus |
CN101350518A (zh) * | 2008-09-05 | 2009-01-21 | 国电南瑞科技股份有限公司 | 数字化保护测控装置采样值处理方法 |
CN101655519A (zh) * | 2009-09-14 | 2010-02-24 | 国电南京自动化股份有限公司 | 数字化变电站测控装置交流采样数据处理方法 |
Non-Patent Citations (1)
Title |
---|
曾院辉等: "一种软件频率跟踪方法", 《电力系统自动化》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103383412A (zh) * | 2013-07-10 | 2013-11-06 | 珠海许继芝电网自动化有限公司 | 一种自适应软硬件频率跟踪采样的方法 |
CN103383412B (zh) * | 2013-07-10 | 2015-08-26 | 珠海许继芝电网自动化有限公司 | 一种自适应软硬件频率跟踪采样的方法 |
CN112485511A (zh) * | 2020-12-11 | 2021-03-12 | 国网辽宁省电力有限公司阜新供电公司 | 基于跟踪补偿算法的交流采样装置及其采样补偿方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Phadke et al. | Synchronized phasor and frequency measurement under transient conditions | |
JP4987068B2 (ja) | 交流電気量測定装置 | |
CN103543335B (zh) | 一种同步相量测量方法 | |
JP6416072B2 (ja) | 同期フェーザ測定装置およびパルス生成装置 | |
CN102135568B (zh) | 交流电信号的周期检测方法及相位差检测方法 | |
CN102288821B (zh) | 三相电路相位差的测量方法、测量装置 | |
CN102236048B (zh) | 一种电力系统相量频率测量方法 | |
CN103941088A (zh) | 一种基于三相信号的电力系统频率快速测量方法 | |
CN102035554A (zh) | 电能计量和分析系统、方法和模数转换电路 | |
CN109633266B (zh) | 一种频率测量方法、系统、装置及计算机可读存储介质 | |
CN102879639A (zh) | 一种电力系统中频率的实时测量方法 | |
CN102135567A (zh) | 一种电力系统交流采样实时频率跟踪及谐波测量方法 | |
CN108333426A (zh) | 基于傅氏算法的电力系统频率测量方法 | |
CN102243272A (zh) | 一种高精度采样数据同步的谐波分析方法 | |
CN103018555A (zh) | 一种高精度的电力参数软件同步采样方法 | |
Mokeev | Filter synthesis for PMU | |
CN103063913A (zh) | 用于傅里叶变换的频率跟踪方法 | |
JP5281126B2 (ja) | レゾルバによる角度検出方法 | |
CN103969508A (zh) | 一种实时高精密的电力谐波分析方法及装置 | |
CN102692540A (zh) | 一种频率跟踪交流采样方法 | |
CN102608415A (zh) | 基于加权双拟合的软件频率跟踪算法 | |
CN103543331A (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN105550516B (zh) | 基于Adams算法的频率自适应谐波间谐波分析方法 | |
CN103592513B (zh) | 电力信号谐波分析方法和装置 | |
CN103884910B (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20120926 |