CN102798748A - 一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法 - Google Patents
一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法 Download PDFInfo
- Publication number
- CN102798748A CN102798748A CN2012102431807A CN201210243180A CN102798748A CN 102798748 A CN102798748 A CN 102798748A CN 2012102431807 A CN2012102431807 A CN 2012102431807A CN 201210243180 A CN201210243180 A CN 201210243180A CN 102798748 A CN102798748 A CN 102798748A
- Authority
- CN
- China
- Prior art keywords
- phase
- voltage
- sampling
- phase voltage
- nth harmonic
- 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
- Measuring Phase Differences (AREA)
- Measurement Of Current Or Voltage (AREA)
Abstract
本发明公开了一种基于迭代傅里叶变换计算交流电压幅值和相位的方法,属于电力供电技术领域。包括:确定采样频率,初始化参数;确定欲提取的交流电压频率;计算正弦函数和余弦函数;采样交流电压值;判断采样序号,计算交流电压次谐波余弦幅值和余弦幅值;获得交流电压次谐波电压幅值和相位。本发明利用每次采样数据进行一次加减法计算即可获得准确的基波电压幅值和相位,较传统的傅里叶变换减少大量的运算过程,能够满足电力系统中信号检测的快速性和准确性要求。
Description
技术领域
本发明属于电力供电技术领域,涉及电力供电系统中动态无功补偿装置设计,具体为一种基于迭代傅里叶变换计算交流电压信号幅值相位的方法。
背景技术
钢铁冶金行业中诸如轧机、电弧炉等冲击负载给电网电压造成闪变,因此需要利用动态无功补偿技术进行实时控制以稳定母线电压。在应用动态无功补偿技术时,需要检测母线电压的基波幅值和相位,而采用传统的基于傅里叶变换获得基波幅值相位的方法需要在每一采样点进行大量的计算才能获得精确的结果。
发明内容
针对现有技术中存在的问题,本发明提出一种基于迭代傅里叶变换计算交流电压幅值和相位的方法,本发明利用每次采样数据进行一次加减法计算即可获得准确的基波电压幅值和相位,较传统的傅里叶变换减少大量的运算过程,能够满足电力系统中信号检测的快速性和准确性要求。
本发明提出一种基于迭代傅里叶变换计算交流电压幅值和相位的方法,其特征在于:包括以下几个步骤:
(1)确定采样频率fs,初始化采样序号k为1,初始化A、B或C相电压n次谐波余弦幅值an为0,初始化A、B或C相电压n次谐波正弦幅值bn为0;
(3)计算正弦函数sin(wnkTs)和余弦函数cos(wnkTs),将计算结果分别保存为两个数组;若采样序号k满足k>N时,从k=N+1开始,令k=1重新开始计数,使sin(wnkTs)和cos(wnkTs)一直反复利用N点采样值循环计算;其中Ts为采样周期,Ts=1/fs;
(4)采样交流电压A、B或C相电压值u(k),若采样序号k满足N≥k>0,则进入步骤(5),若采样序号k满足k>N,则进入步骤(6);
(5)根据下式得到A、B或C相电压n次谐波余弦幅值an和A、B或C相电压n次谐波正弦幅值bn,然后进入步骤(7);
其中u(t)表示任意时刻t时A、B或C相电压,wn=nw为A、B或C相电压n次谐波的角速度,w为A、B或C相基波电压信号的角速度,w=2πf,f为基波频率,t0表示大于等于零的任意时刻,Ts为采样周期,Ts=1/fs,k为采样序号;
(6)根据下式得到A、B或C相电压n次谐波余弦幅值an和A、B或C相电压n次谐波正弦幅值bn,然后进入步骤(7);
得到A、B或C相电压n次谐波电压初始相位;然后跳转至步骤(4),采样下一时刻的A、B或C相电压信号。
本发明具有的优点在于:
1、本发明提出一种基于迭代傅里叶变换计算交流电压幅值和相位的方法,在每次采样间隔只需要利用新的采样数据进行一次加减法计算即可获得准确的基波电压幅值和相位。
2、本发明提出一种基于迭代傅里叶变换计算交流电压幅值和相位的方法,较传统的基于过零法的相位检测方法具有更快的锁相速度。
3、本发明提出一种基于迭代傅里叶变换计算交流电压幅值和相位的方法,能够在一个周期内获得准确的幅值和相位变化情况。
附图说明
图1为本发明提出的基于迭代傅里叶变换计算交流电压幅值和相位的方法流程图;
图2为实施例1中交流电压信号中A相电压的实测采样曲线;
图3为实施例1中计算的正弦函数sin(w1kTs)和余弦函数cos(w1kTs)cos(w1kTs)曲线;
图4为实施例1中A相电压幅值变化曲线;
图5为实施例1中A相电压相位变化曲线。
具体实施方式
下面结合附图对本发明作进一步说明,以使本领域的技术人员可以更好的理解本发明并能予以实施。
本发明提出的一种基于迭代傅里叶变换计算交流电压幅值和相位的方法,该方法通过采样交流电压A、B、C三相电压信号,计算A、B、C三相的电压信号幅值和相位,交流电压每一相电压信号幅值和相位的计算方法完全相同,因此本发明给出的一种基于迭代傅里叶变换计算交流电压幅值和相位的方法,为计算交流电压信号中任意一相电压幅值和相位的方法。
电力系统中由于存在非线性元件,交流电压中往往含有大量谐波,因此检测电力系统电压基波分量时需要采用傅里叶变换才能获得准确的交流电压幅值和相位,设被检测交流电压的A、B或C相电压u的表达式为:
其中An、nw、分别为A、B或C相n次谐波的电压幅值、角频率、初始相位,w=2πf,w为A、B或C相基波电压信号的角速度,n为大于1的正整数,f为A、B或C相基波电压信号频率。在离散系统中,若采样频率为fs,则采样周期Ts=1/fs,每一个基波周期的采样点数为N=fs/f,按照离散傅里叶变换的基本原理,在任意时刻t,A、B或C相电压u(t)表示为:
其中:
a0为A、B或C相电压直压分量,an为A、B或C相电压n次谐波余弦幅值,bn为A、B或C相电压n次谐波正弦幅值,t0表示大于等于零的任意时刻。根据公式(3)、公式(4)和公式(5),在采样时刻t0+kTs,且N≥k>0时,an和bn可以分别表示为:
其中k为采样序号。
比较公式(4)、(5)与公式(6)、(7)可以发现,公式(6)、(7)中等式右边前半部分为利用时刻t0+kTs之前的采样数据进行计算的结果,后半部分为当前时刻采样数据计算结果,因此将A、B或C相电压n次谐波余弦幅值an和A、B或C相电压n次谐波正弦幅值bn的计算均转化为迭代形式的单步计算。
当k>N时,在采样时刻t0+kTs,根据迭代原理,A、B或C相电压n次谐波余弦幅值an和A、B或C相电压n次谐波正弦幅值bn可以分别表示为:
an(k)、bn(k)分别为第k次采样获得的A、B或C相电压的n次谐波余弦幅值、正弦幅值,根据(6)~(9),时刻t0为任意值,对于实际采样数据无法确定时刻t0,不失一般意义下可设t0=0,因此,cos(nwkTs)和sin(nwkTs)可转化为一个信号周期中的采样点。以基波电压信号为例,正弦函数sin(wkTs)和余弦函数cos(wkTs)分为N个采样点,由于sin(wkTs)和cos(wkTs)是以2π/w为周期的函数,因此在k>N,从k=N+1开始,可令k=1重新开始计数,从而使得sin(wkTs)和cos(wkTs)一直反复利用N点采样值循环计算。
根据公式(2),n次谐波A、B或C相电压un(t)为:
综上,本发明提出一种基于迭代傅里叶变换计算交流电压幅值和相位的方法,如图1所示,具体为:
(1)确定采样频率fs,初始化参数,令k=1,an=0,bn=0;
(2)确定需要提取的A、B或C相电压n次谐波的信号频率fn,A、B或C相电压n次谐波的角速度为wn=2πfn,计算提取信号一个周期的采样点数
(3)计算正弦函数sin(wnkTs)和余弦函数cos(wnkTs),其中k≤N,将计算结果分别保存为两个数组;当k>N时,从k=N+1开始,令k=1重新开始计数,从而使得sin(wnkTs)和cos(wnkTs)一直反复利用N点采样值循环计算。
(4)采样A、B或C相电压值u(k),判断k是否满足N≥k>0,若满足则进入步骤(5),否则进入步骤(6);
(5)根据公式(6)、(7)计算an、bn,然后进入步骤(7);
(6)根据公式(8)、(9)计算an、bn,然后进入步骤(7);
(7)根据公式(11)、(12)计算A、B或C相电压n次谐波的幅值和初始相位,然后跳转至步骤(4),采样下一时刻的A、B或C相电压信号。
实施例:
某钢厂6.5kV母线上接有一负荷系统,通过实际测量获得如图2所示的A相电压曲线,此时采样频率fs=10000Hz,采样周期Ts=0.0001s。若需要从检测电压数据中获得基波电压(即n=1)幅值和相位,则本实施例中需要提取的基波电压信号的频率为f1=50Hz,角频率为w1=2πf1=314.15rad/s,一个周期采样点数
当N≥k>0时,计算正弦函数sin(w1kTs)和余弦函数cos(w1kTs)值,具体数值曲线如图3所示。
按照本发明方法中步骤(4)~(7)进行判断计算,可得最终稳定状态时A相基波电压正弦幅值a1=284020,余弦幅值b1=-874120,则A相电压幅值初始相位交流电压中B、C两相电压的幅值和相位的计算方法与A相电压的幅值和相位的计算方法完全相同,只需将采样数据分别换成B、C相的电压数据即可。本实施例中A相电压的幅值变化曲线如图4所示,A相电压相位的变化曲线如图5所示。从图4中可以看出本发明提出的电压幅值计算方法只需要一个周期采样数据即能收敛于实际电压幅值,从图5中可以看出本发明提出的电压相位计算方法也只需要一个周期采样数据即能收敛于实际电压相位值。因此本发明提出一种基于迭代傅里叶变换计算交流电压幅值和相位的方法较传统的基于过零法的相位检测方法具有更快的锁相速度。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。
Claims (1)
1.一种基于迭代傅里叶变换计算交流电压幅值和相位的方法,其特征在于:包括以下几个步骤:
(1)确定采样频率fs,初始化采样序号k为1,初始化A、B或C相电压n次谐波余弦幅值an为0,初始化A、B或C相电压n次谐波正弦幅值bn为0;
(3)计算正弦函数sin(wnkTs)和余弦函数cos(wnkTs),将计算结果分别保存为两个数组;若采样序号k满足k>N时,从k=N+1开始,令k=1重新开始计数,使sin(wnkTs)和cos(wnkTs)一直反复利用N点采样值循环计算;其中Ts为采样周期,Ts=1/fs;
(4)采样交流电压A、B或C相电压值u(k),若采样序号k满足N≥k>0,则进入步骤(5),若采样序号k满足k>N,则进入步骤(6);
(5)根据下式得到A、B或C相电压n次谐波余弦幅值an和A、B或C相电压n次谐波正弦幅值bn,然后进入步骤(7);
其中u(t)表示任意时刻t时A、B或C相电压,wn=nw为A、B或C相电压n次谐波的角速度,w为A、B或C相基波电压信号的角速度,w=2πf,f为A、B或C相基波电压信号频率,t0表示大于等于零的任意时刻,Ts为采样周期,Ts=1/fs,k为采样序号;
(6)根据下式得到A、B或C相电压n次谐波余弦幅值an和A、B或C相电压n次谐波正弦幅值bn,然后进入步骤(7);
得到A、B或C相电压n次谐波电压初始相位;然后跳转至步骤(4),采样下一时刻的A、B或C相电压信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210243180.7A CN102798748B (zh) | 2012-07-13 | 2012-07-13 | 一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210243180.7A CN102798748B (zh) | 2012-07-13 | 2012-07-13 | 一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102798748A true CN102798748A (zh) | 2012-11-28 |
CN102798748B CN102798748B (zh) | 2014-08-13 |
Family
ID=47197915
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210243180.7A Active CN102798748B (zh) | 2012-07-13 | 2012-07-13 | 一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102798748B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103344815A (zh) * | 2013-06-08 | 2013-10-09 | 中国农业大学 | 一种宽范围变化的电气参数测量方法及系统 |
CN105067882A (zh) * | 2015-07-24 | 2015-11-18 | 许继电气股份有限公司 | 一种电量幅值的测定方法 |
CN105372471A (zh) * | 2015-07-09 | 2016-03-02 | 深圳市科润宝实业有限公司 | 正弦信号的幅值检测方法和系统 |
CN107167658A (zh) * | 2017-05-11 | 2017-09-15 | 南京国电南自电网自动化有限公司 | 一种高精度抗干扰的电力系统基波频率及相位测量方法 |
CN107328987A (zh) * | 2017-06-27 | 2017-11-07 | 国网福建省电力有限公司 | 一种光电式电流传感器波形失真的处理方法及装置 |
CN108414816A (zh) * | 2018-01-31 | 2018-08-17 | 许继集团有限公司 | 一种交流单相电压的解算方法及装置 |
CN109164293A (zh) * | 2018-07-10 | 2019-01-08 | 深圳市华宝新能源股份有限公司 | 一种市电掉电检测方法、电子设备、存储介质和储能系统 |
CN111134637A (zh) * | 2019-12-06 | 2020-05-12 | 深圳绿米联创科技有限公司 | 信号处理方法、装置、电子设备及存储介质 |
CN112255457A (zh) * | 2020-09-22 | 2021-01-22 | 天津电气科学研究院有限公司 | 适用于自动准同期装置的相角差测量方法 |
CN112268501A (zh) * | 2020-09-03 | 2021-01-26 | 西安铁路信号有限责任公司 | 一种可应用于物体线性位移或转角位置的检测方法 |
CN112710890A (zh) * | 2020-12-15 | 2021-04-27 | 特变电工西安电气科技有限公司 | 一种单相正弦交流相量实时计算方法及装置 |
CN113180689A (zh) * | 2021-04-08 | 2021-07-30 | 上海数创医疗科技有限公司 | 一种基于高斯函数拟合的p波波峰及起始位置提取方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008120282A1 (ja) * | 2007-02-27 | 2008-10-09 | Mitsubishi Electric Corporation | 交流電気量測定装置 |
CN101408568A (zh) * | 2007-10-12 | 2009-04-15 | 深圳科士达科技股份有限公司 | 一种测量交流电的电信号相位的方法及装置 |
CN102495282A (zh) * | 2011-12-27 | 2012-06-13 | 国电南瑞科技股份有限公司 | 一种交流电压频率的测量方法 |
CN102508029A (zh) * | 2011-09-26 | 2012-06-20 | 北京东标电气股份有限公司 | 一种用于电网相位角跟踪的方法 |
-
2012
- 2012-07-13 CN CN201210243180.7A patent/CN102798748B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008120282A1 (ja) * | 2007-02-27 | 2008-10-09 | Mitsubishi Electric Corporation | 交流電気量測定装置 |
CN101408568A (zh) * | 2007-10-12 | 2009-04-15 | 深圳科士达科技股份有限公司 | 一种测量交流电的电信号相位的方法及装置 |
CN102508029A (zh) * | 2011-09-26 | 2012-06-20 | 北京东标电气股份有限公司 | 一种用于电网相位角跟踪的方法 |
CN102495282A (zh) * | 2011-12-27 | 2012-06-13 | 国电南瑞科技股份有限公司 | 一种交流电压频率的测量方法 |
Non-Patent Citations (1)
Title |
---|
张同尊: "《一种基于离散傅里叶变换的频率测量算法》", 《电力系统自动化》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103344815B (zh) * | 2013-06-08 | 2015-08-19 | 中国农业大学 | 一种宽范围变化的电气参数测量方法及系统 |
CN103344815A (zh) * | 2013-06-08 | 2013-10-09 | 中国农业大学 | 一种宽范围变化的电气参数测量方法及系统 |
CN105372471B (zh) * | 2015-07-09 | 2018-05-25 | 深圳市科润宝实业有限公司 | 正弦信号的幅值检测方法和系统 |
CN105372471A (zh) * | 2015-07-09 | 2016-03-02 | 深圳市科润宝实业有限公司 | 正弦信号的幅值检测方法和系统 |
CN105067882A (zh) * | 2015-07-24 | 2015-11-18 | 许继电气股份有限公司 | 一种电量幅值的测定方法 |
CN105067882B (zh) * | 2015-07-24 | 2017-11-07 | 许继电气股份有限公司 | 一种电量幅值的测定方法 |
CN107167658B (zh) * | 2017-05-11 | 2019-06-28 | 南京国电南自电网自动化有限公司 | 一种高精度抗干扰的电力系统基波频率及相位测量方法 |
CN107167658A (zh) * | 2017-05-11 | 2017-09-15 | 南京国电南自电网自动化有限公司 | 一种高精度抗干扰的电力系统基波频率及相位测量方法 |
CN107328987A (zh) * | 2017-06-27 | 2017-11-07 | 国网福建省电力有限公司 | 一种光电式电流传感器波形失真的处理方法及装置 |
CN108414816A (zh) * | 2018-01-31 | 2018-08-17 | 许继集团有限公司 | 一种交流单相电压的解算方法及装置 |
CN108414816B (zh) * | 2018-01-31 | 2020-08-11 | 许继集团有限公司 | 一种交流单相电压的解算方法及装置 |
CN109164293A (zh) * | 2018-07-10 | 2019-01-08 | 深圳市华宝新能源股份有限公司 | 一种市电掉电检测方法、电子设备、存储介质和储能系统 |
CN111134637A (zh) * | 2019-12-06 | 2020-05-12 | 深圳绿米联创科技有限公司 | 信号处理方法、装置、电子设备及存储介质 |
CN112268501B (zh) * | 2020-09-03 | 2022-03-04 | 西安铁路信号有限责任公司 | 一种可应用于物体线性位移或转角位置的检测方法 |
CN112268501A (zh) * | 2020-09-03 | 2021-01-26 | 西安铁路信号有限责任公司 | 一种可应用于物体线性位移或转角位置的检测方法 |
CN112255457A (zh) * | 2020-09-22 | 2021-01-22 | 天津电气科学研究院有限公司 | 适用于自动准同期装置的相角差测量方法 |
CN112255457B (zh) * | 2020-09-22 | 2022-06-07 | 天津电气科学研究院有限公司 | 适用于自动准同期装置的相角差测量方法 |
CN112710890A (zh) * | 2020-12-15 | 2021-04-27 | 特变电工西安电气科技有限公司 | 一种单相正弦交流相量实时计算方法及装置 |
CN112710890B (zh) * | 2020-12-15 | 2024-02-27 | 特变电工西安电气科技有限公司 | 一种单相正弦交流相量实时计算方法及装置 |
CN113180689A (zh) * | 2021-04-08 | 2021-07-30 | 上海数创医疗科技有限公司 | 一种基于高斯函数拟合的p波波峰及起始位置提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102798748B (zh) | 2014-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102798748B (zh) | 一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法 | |
CN102818921B (zh) | 一种基于迭代傅里叶变换计算交流电流信号幅值和相位的方法 | |
CN101479612B (zh) | 频率测量装置 | |
CN101957398B (zh) | 一种基于机电与电磁暂态混合仿真技术检测计算电网一次时间常数的方法 | |
CN101806832B (zh) | 一种低频率信号的频率测量方法 | |
US20140062457A1 (en) | Three-phase ac phase sequence detecting method and apparatus | |
RU2014120171A (ru) | Способ и устройство для ввода электрического тока в электрическую сеть | |
CN102435844A (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN102735938A (zh) | 一种电网电压基波正序相角的快速检测方法 | |
CN103558436A (zh) | 基于单相锁相环算法的检测电网电压幅值、频率和相角的方法 | |
CN108808669A (zh) | 高压直流输电系统换流器的动态相量建模方法 | |
CN104535890A (zh) | 一种电网电压跌落的快速检测方法 | |
CN105158592A (zh) | 电网不对称条件下三相电网阻抗实时检测方法 | |
KR101952063B1 (ko) | 누설 전류 검출 장치 | |
CN104360137A (zh) | 一种适用于动态电压恢复器的电压暂降检测方法 | |
CN102590618A (zh) | 一种电网基波正序电压相位的检测方法 | |
CN102508022A (zh) | 采用最优乘子牛顿算法检测电网频率的方法 | |
CN104101781B (zh) | 变电站母线电压相位角瞬时值测量方法 | |
CN102809688B (zh) | 一种基于迭代傅里叶变换计算三相电压实时相位的方法 | |
CN104868909B (zh) | 一种基于电压正交谐振器qr无静差的锁频锁相环及其测量方法 | |
CN107703358B (zh) | 一种基于改进二阶广义积分器的锁相算法 | |
CN107797017B (zh) | 一种电力变压器带电检测损耗特征参数的方法 | |
CN103293379B (zh) | 基于有效值的apf谐波检测方法及其直流侧电压的控制方法 | |
CN105137186A (zh) | 一种微机自动准同期装置的同期电压相位差测量法 | |
CN104090143A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |