CN112710890B - 一种单相正弦交流相量实时计算方法及装置 - Google Patents
一种单相正弦交流相量实时计算方法及装置 Download PDFInfo
- Publication number
- CN112710890B CN112710890B CN202011479222.8A CN202011479222A CN112710890B CN 112710890 B CN112710890 B CN 112710890B CN 202011479222 A CN202011479222 A CN 202011479222A CN 112710890 B CN112710890 B CN 112710890B
- Authority
- CN
- China
- Prior art keywords
- phase
- alternating current
- signal
- sinusoidal alternating
- real
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 56
- 238000005070 sampling Methods 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 15
- 230000010363 phase shift Effects 0.000 claims description 9
- 238000001914 filtration Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/04—Measuring peak values or amplitude or envelope of ac or of pulses
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R25/00—Arrangements for measuring phase angle between a voltage and a current or between voltages or currents
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Networks Using Active Elements (AREA)
Abstract
本发明公开了一种单相正弦交流相量实时计算方法及装置,装置包括移相滤波器和幅值、相位运算模块,其中低通滤波器用于单相正弦信号实时滤波和构建正交信号,幅值、相位运算模块根据输入正弦交流信号及其正交信号实时计算相量。本发明可以实时计算单相正弦交流相量,并且算法结构简单、易于实现。与传统的单相正弦交流相量计算方法相比,所需采样数据远少于1/4周期采样数,运算量小,计算速度快,可用于实时计算。
Description
技术领域
本发明属于电气或电子***测量控制技术领域,具体涉及一种单相正弦交流相量实时计算方法及装置。
背景技术
电气或电子***测量和控制设备中,通常需要检测正弦交流信号的幅值和相位,即获取正弦交流相量。三相正弦交流信号可通过基于dq变换的三相锁相环获取其幅值和相位;但单相正弦交流信号通常只能通过傅里叶变换或虚拟dq变换方式计算正弦交流相量。傅里叶变换通常可采用离散傅里叶变换(DFT)或快速傅里叶变换(FFT),根据输入信号整周期或半周期的采样数据计算,需要存储大量数据且计算复杂。虚拟dq变换的方式则需要对输入信号进行存储和延时来获取正交信号,一般至少要存储1/4周期的采样数据。然而实时***的存储资源和计算资源有限,当信号采样率较高时,这两类传统方法均需要耗费大量的存储空间和计算资源来满足单相正弦交流相量的计算需求,甚至无法进行实时计算。
发明内容
本发明所要解决的技术问题是针对现有技术所存在的上述不足,提供一种单相正弦交流相量实时计算方法,该计算方法所需采样数据远少于1/4周期采样数,运算量小,计算速度快,可用于实时计算。
为实现上述目的,采用如下技术方案:
一种单相正弦交流相量实时计算方法,包括:
构建输入正弦交流信号的正交信号;
根据输入正弦交流信号及其正交信号实时计算相量。
优选的,采用移相滤波器构建输入正弦交流信号的正交信号,其传递函数为:
将移相滤波器的传递函数离散化,得到实时计算方程为:
uf(k)=a0u(k)+a1u(k-1)+a2u(k-2)-b1uf(k-1)-b2uf(k-2)
式中:u(k)为输入正弦交流信号在当前时刻的采样值,uf(k)为输入正弦交流信号在当前时刻的正交信号计算值,a0、a1、a2、b1、b2为移相滤波器系数;
按照移相滤波器实时计算方程对输入正弦交流信号的采样值进行迭代计算,得到构建后的输入正弦交流信号的正交信号。
优选的,所述移相滤波器满足:在输入正弦交流信号的频率f0处幅频特性为0dB、相频特性为-90°。
优选的,所述移相滤波器采用二阶低通滤波器。
优选的,根据输入正弦交流信号及其正交信号实时计算相量,计算公式为:
式中:Uk为所述单相正弦交流相量,u(k)为输入正弦交流信号在当前时刻的采样值,uf(k)为输入正弦交流信号在当前时刻的正交信号计算值,Arg(-uf(k)+j·u(k))为复数(-uf(k)+j·u(k))的辐角。
本发明提供的另一个技术方案是:
一种单相正弦交流相量实时计算装置,包括:
移相滤波器,用于构建输入正弦交流信号的正交信号;
幅值相位运算模块,用于根据输入正弦交流信号及其正交信号实时计算相量。
优选的,所述移相滤波器满足:在输入正弦交流信号的频率f0处幅频特性为0dB、相频特性为-90°。
优选的,所述移相滤波器采用二阶低通滤波器。
优选的,所述移相滤波器的传递函数为:
将移相滤波器的传递函数离散化,得到实时计算方程为:
uf(k)=a0u(k)+a1u(k-1)+a2u(k-2)-b1uf(k-1)-b2uf(k-2)
式中:u(k)为输入正弦交流信号在当前时刻的采样值,uf(k)为输入正弦交流信号在当前时刻的正交信号计算值,a0、a1、a2、b1、b2为移相滤波器系数;
按照移相滤波器实时计算方程对输入正弦交流信号的采样值进行迭代计算,得到构建后的输入正弦交流信号的正交信号。
优选的,所述幅值相位运算模块根据输入正弦交流信号及其正交信号实时计算相量的计算公式为:
式中:Uk为所述单相正弦交流相量,u(k)为输入正弦交流信号在当前时刻的采样值,uf(k)为输入正弦交流信号在当前时刻的正交信号计算值,Arg(-uf(k)+j·u(k))为复数(-uf(k)+j·u(k))的辐角。
本发明的有益效果如下:
1、本发明实施例设计的移相滤波器可有效实现输入的单相正弦交流信号的正交信号,即相位相差90°的正弦交流信号,同时滤除了输入信号中的高频分量,优于传统方法中依靠延时实现正交信号。
2、本发明实施例与传统方法相比,在实时计算过程中只需要6个数据,占用的存储空间小。
3、本发明实施例与传统方法相比,本发明中移相滤波和幅值相位运算模块的公式简单,运算量小,占用计算资源少。
附图说明
构成本申请的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本发明实施例中单相正弦交流相量实时计算方法组成模块示意图;
图2为本发明实施例中移相滤波器的一种实施方案的频率特性图;
图3为本发明实施例中单相正弦交流相量实时计算方法的一种实施效果图。
具体实施方式
下面将参考附图并结合实施例来详细说明本发明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
以下详细说明均是示例性的说明,旨在对本发明提供进一步的详细说明。除非另有指明,本发明所采用的所有技术术语与本申请所属领域的一般技术人员的通常理解的含义相同。本发明所使用的术语仅是为了描述具体实施方式,而并非意图限制根据本发明的示例性实施方式。
(一)本发明实施例提供了一种单相正弦交流相量实时计算方法及装置,算法结构简单、易于实现。与传统的单相正弦交流相量计算方法相比,在采样率较高时所需采样数据远少于1/4周期采样数,运算量小,计算速度快,可用于实时计算。
如图1所示,该装置包括移相滤波器和幅值相位运算模块,该方法通过移相滤波器构建输入正弦交流信号的正交信号,通过幅值相位运算模块根据输入正弦交流信号及其正交信号实时计算相量。
移相滤波器满足在输入正弦交流信号的频率f0处幅频特性为0dB,相频特性为-90°,即:
式中:为所述移相滤波器的频率特性函数,A(ω)为所述移相滤波器的幅频特性函数,/>为所述移相滤波器的相频特性函数,f0为输入正弦交流信号的频率。
A1.根据所述的移相滤波器特征设计移相滤波器传递函数;典型的,所述移相滤波器可采用满足条件的二阶低通滤波器,其传递函数为:
A2.将移相滤波器传递函数离散化,得到实时计算方程;上述传递函数对应的实时计算方程为:
uf(k)=a0u(k)+a1u(k-1)+a2u(k-2)-b1uf(k-1)-b2uf(k-2)
式中:u(k)为输入正弦交流信号在当前时刻的采样值,uf(k)为输入正弦交流信号在当前时刻的正交信号计算值,a0、a1、a2、b1、b2为移相滤波器系数。
A3.存储输入正弦交流信号的采样值,按照移相滤波器实时计算方程进行迭代计算,存储计算结果;
A4.重复A3步骤。
所述幅值相位运算模块根据输入正弦交流信号及其正交信号实时计算相量。计算公式为:
式中:Uk为所述单相正弦交流相量,u(k)为输入正弦交流信号在当前时刻的采样值,uf(k)为输入正弦交流信号在当前时刻的正交信号计算值,Arg(-uf(k)+j·u(k))为复数(-uf(k)+j·u(k))的辐角。
(二)下面结合具体的实施例子对本发明实施例提供的技术方案进行解释和说明:
在本发明的一种实施方式中采样率设为20000采样点/秒,输入单相正弦交流信号的频率为f0=50Hz,满足如下方程:
所述的移相滤波器,按照以下步骤实施:
A1.根据所述的移相滤波器特征设计移相滤波器传递函数。所述移相滤波器的一种实施方案为满足条件的二阶低通滤波器,其传递函数为:
该滤波器满足所述的移相滤波器特征,频率特性图如图2所示,在输入正弦交流信号的频率f0=50Hz处幅频特性为0dB,相频特性为-90°,即:
A2.将移相滤波器传递函数离散化,得到实时计算方程。上述传递函数对应的实时计算方程为:
uf(k)=8.6272237×10-5u(k)+1.7254447×10-4u(k-1)+8.6272237×10-5u(k-2)+1.9777869uf(k-1)-0.9780310uf(k-2)
A3.存储输入正弦交流信号的采样值u(k)、u(k-1)、u(k-2),按照移相滤波器实时计算方程进行迭代计算,得到正交信号计算值uf(k),存储计算结果uf(k)用于下次计算;
A4.重复A3步骤。
u(k)和uf(k)传递给所述幅值相位运算模块,根据输入正弦交流信号及其正交信号实时计算相量。计算公式为:
如图3所示,自上而下依次为u(k)、uf(k)、Uk幅值、Uk相位的计算结果。由图3可知,本发明所述的移相滤波器有效实时输出了输入单相正弦交流信号的正交信号,所述幅值相位运算模块在有效实时计算出单相正弦交流相量的幅值和相位。整个计算过程只需要存储6个数据,且单次运算量小。传统计算方法在该实施实例设定信号下至少要存储100个采样点数据用于计算,且运算量大,本发明所述方法具有明显优势。
由技术常识可知,本发明可以通过其它的不脱离其精神实质或必要特征的实施方案来实现。因此,上述公开的实施方案,就各方面而言,都只是举例说明,并不是仅有的。所有在本发明范围内或在等同于本发明的范围内的改变均被本发明包含。
Claims (6)
1.一种单相正弦交流相量实时计算方法,其特征在于,包括:
构建输入正弦交流信号的正交信号;
根据输入正弦交流信号及其正交信号实时计算相量;
采用移相滤波器构建输入正弦交流信号的正交信号,其传递函数为:
将移相滤波器的传递函数离散化,得到实时计算方程为:
uf(k)=a0u(k)+a1u(k-1)+a2u(k-2)-b1uf(k-1)-b2uf(k-2)
式中:u(k)为输入正弦交流信号在当前时刻的采样值,uf(k)为输入正弦交流信号在当前时刻的正交信号计算值,a0、a1、a2、b1、b2为移相滤波器系数;
按照移相滤波器实时计算方程对输入正弦交流信号的采样值进行迭代计算,得到构建后的输入正弦交流信号的正交信号;
根据输入正弦交流信号及其正交信号实时计算相量,计算公式为:
式中:Uk为所述单相正弦交流相量,u(k)为输入正弦交流信号在当前时刻的采样值,uf(k)为输入正弦交流信号在当前时刻的正交信号计算值,Arg(-uf(k)+j·u(k))为复数(-uf(k)+j·u(k))的辐角。
2.根据权利要求1所述的单相正弦交流相量实时计算方法,其特征在于,所述移相滤波器满足:在输入正弦交流信号的频率f0处幅频特性为0dB、相频特性为-90°。
3.根据权利要求1所述的单相正弦交流相量实时计算方法,其特征在于,所述移相滤波器采用二阶低通滤波器。
4.一种单相正弦交流相量实时计算装置,其特征在于,包括:
移相滤波器,用于构建输入正弦交流信号的正交信号;
幅值相位运算模块,用于根据输入正弦交流信号及其正交信号实时计算相量;
所述移相滤波器的传递函数为:
将移相滤波器的传递函数离散化,得到实时计算方程为:
uf(k)=a0u(k)+a1u(k-1)+a2u(k-2)-b1uf(k-1)-b2uf(k-2)
式中:u(k)为输入正弦交流信号在当前时刻的采样值,uf(k)为输入正弦交流信号在当前时刻的正交信号计算值,a0、a1、a2、b1、b2为移相滤波器系数;
按照移相滤波器实时计算方程对输入正弦交流信号的采样值进行迭代计算,得到构建后的输入正弦交流信号的正交信号;
所述幅值相位运算模块根据输入正弦交流信号及其正交信号实时计算相量的计算公式为:
式中:Uk为所述单相正弦交流相量,u(k)为输入正弦交流信号在当前时刻的采样值,uf(k)为输入正弦交流信号在当前时刻的正交信号计算值,Arg(-uf(k)+j·u(k))为复数(-uf(k)+j·u(k))的辐角。
5.根据权利要求4所述的单相正弦交流相量实时计算装置,其特征在于,所述移相滤波器满足:在输入正弦交流信号的频率f0处幅频特性为0dB、相频特性为-90°。
6.根据权利要求4所述的单相正弦交流相量实时计算装置,其特征在于,所述移相滤波器采用二阶低通滤波器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011479222.8A CN112710890B (zh) | 2020-12-15 | 2020-12-15 | 一种单相正弦交流相量实时计算方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011479222.8A CN112710890B (zh) | 2020-12-15 | 2020-12-15 | 一种单相正弦交流相量实时计算方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112710890A CN112710890A (zh) | 2021-04-27 |
CN112710890B true CN112710890B (zh) | 2024-02-27 |
Family
ID=75543257
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011479222.8A Active CN112710890B (zh) | 2020-12-15 | 2020-12-15 | 一种单相正弦交流相量实时计算方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112710890B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004289576A (ja) * | 2003-03-24 | 2004-10-14 | Hitachi Kokusai Electric Inc | カーテシアン型送信機 |
CN1677844A (zh) * | 2004-04-02 | 2005-10-05 | 泰瑞达公司 | 高性能信号的发生 |
CN102798748A (zh) * | 2012-07-13 | 2012-11-28 | 中冶南方工程技术有限公司 | 一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法 |
JP2013195222A (ja) * | 2012-03-19 | 2013-09-30 | Daihen Corp | 交流信号測定装置 |
CN104181577A (zh) * | 2014-09-01 | 2014-12-03 | 中国科学技术大学 | 一种基于全数字化技术的束流位置和相位测量***及方法 |
CN106410810A (zh) * | 2016-11-24 | 2017-02-15 | 国网江苏省电力公司经济技术研究院 | 一种附加阻尼控制的upfc控制单元 |
JP2017215774A (ja) * | 2016-05-31 | 2017-12-07 | 一般財団法人電力中央研究所 | 過渡現象解析装置、方法、及びプログラム |
CN109214371A (zh) * | 2018-10-30 | 2019-01-15 | 南京航空航天大学 | 一种基于快速傅氏变换的正弦信号滤波*** |
CN109842143A (zh) * | 2019-02-19 | 2019-06-04 | 郑州电力高等专科学校 | 电压源换流器高压直流输电技术互连弱交流***控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6934654B2 (en) * | 2003-03-21 | 2005-08-23 | Schweitzer Engineering Laboratories, Inc. | System and method for exact compensation of fundamental phasors |
-
2020
- 2020-12-15 CN CN202011479222.8A patent/CN112710890B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004289576A (ja) * | 2003-03-24 | 2004-10-14 | Hitachi Kokusai Electric Inc | カーテシアン型送信機 |
CN1677844A (zh) * | 2004-04-02 | 2005-10-05 | 泰瑞达公司 | 高性能信号的发生 |
JP2013195222A (ja) * | 2012-03-19 | 2013-09-30 | Daihen Corp | 交流信号測定装置 |
CN102798748A (zh) * | 2012-07-13 | 2012-11-28 | 中冶南方工程技术有限公司 | 一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法 |
CN104181577A (zh) * | 2014-09-01 | 2014-12-03 | 中国科学技术大学 | 一种基于全数字化技术的束流位置和相位测量***及方法 |
JP2017215774A (ja) * | 2016-05-31 | 2017-12-07 | 一般財団法人電力中央研究所 | 過渡現象解析装置、方法、及びプログラム |
CN106410810A (zh) * | 2016-11-24 | 2017-02-15 | 国网江苏省电力公司经济技术研究院 | 一种附加阻尼控制的upfc控制单元 |
CN109214371A (zh) * | 2018-10-30 | 2019-01-15 | 南京航空航天大学 | 一种基于快速傅氏变换的正弦信号滤波*** |
CN109842143A (zh) * | 2019-02-19 | 2019-06-04 | 郑州电力高等专科学校 | 电压源换流器高压直流输电技术互连弱交流***控制方法 |
Non-Patent Citations (3)
Title |
---|
An adaptive quasi-notch filter for a biased sinusoidal signal estimation;Giuseppe Fedele 等;《2011 9th IEEE International Conference on Control and Automation (ICCA)》;第1060-1065页 * |
基于数学形态学滤波的电压暂降检测方法;康军胜 等;《电力科学与工程》;第32卷(第1期);第58-63页 * |
频率估计的多段差频正弦信号加权融合算法;肖玮 等;《中国科学技术大学学报》;第42卷(第02期);第124-132页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112710890A (zh) | 2021-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Abdollahi et al. | Frequency estimation: A least-squares new approach | |
CA2929382C (en) | Method of removing decaying dc component from power system fault signal | |
Sinha et al. | A pre-filter based PLL for three-phase grid connected applications | |
JP5855886B2 (ja) | 周波数検出装置 | |
CN108020721B (zh) | 一种基于IpDFT的非平衡电力***的频率估计方法 | |
Meral | Improved phase-locked loop for robust and fast tracking of three phases under unbalanced electric grid conditions | |
CN105529950A (zh) | 一种基于二阶广义积分器的单相并网逆变器控制方法 | |
CN107423261B (zh) | 非理想微电网条件下基于ovpr的正负序分量的分离方法 | |
Terriche et al. | Matrix pencil method‐based reference current generation for shunt active power filters | |
CN105048995B (zh) | 巴特沃斯数字滤波器及利用其实现频率自适应的滤波方法 | |
WO2014094317A1 (zh) | 基于模型预测控制的主动前端整流器滤波延迟补偿方法 | |
CN103534601A (zh) | 用于实时谐波频谱分析仪的装置和方法 | |
CN102095915B (zh) | 一种采用多同步参考坐标系变换的电压信号检测装置 | |
CN105203837A (zh) | 无功功率测量方法 | |
Duda et al. | FIR filters compliant with the IEEE standard for M class PMU | |
Vekić et al. | Multi–Resonant observer PLL with real-time estimation of grid unbalances | |
Zadeh et al. | Analysis of phasor measurement method in tracking the power frequency of distorted signals | |
CN112710890B (zh) | 一种单相正弦交流相量实时计算方法及装置 | |
CN110365018A (zh) | 一种自适应宽频带正负序分离方法 | |
Marčetić et al. | Unbalanced three‐phase distribution system frequency estimation using least mean squares method and positive voltage sequence | |
Destro et al. | Implementation aspects of adaptive window moving average filter applied to PLLs—Comparative study | |
Wu et al. | Effect of adding DC‐offset estimation integrators in there‐phase enhanced phase‐locked loop on dynamic performance and alternative scheme | |
Feng et al. | A novel DFT algorithm used in active power filter under frequency distortion | |
Radil et al. | Methods for estimation of voltage harmonic components | |
CN108957118A (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 |