CN104502675A - 电力信号的基波幅值测量方法和系统 - Google Patents

电力信号的基波幅值测量方法和系统 Download PDF

Info

Publication number
CN104502675A
CN104502675A CN201410848173.9A CN201410848173A CN104502675A CN 104502675 A CN104502675 A CN 104502675A CN 201410848173 A CN201410848173 A CN 201410848173A CN 104502675 A CN104502675 A CN 104502675A
Authority
CN
China
Prior art keywords
sequence
vector
filtering
real number
wave
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
Application number
CN201410848173.9A
Other languages
English (en)
Other versions
CN104502675B (zh
Inventor
李军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Electric Power Research Institute of Guangdong Power Grid Co Ltd filed Critical Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority to CN201410848173.9A priority Critical patent/CN104502675B/zh
Publication of CN104502675A publication Critical patent/CN104502675A/zh
Application granted granted Critical
Publication of CN104502675B publication Critical patent/CN104502675B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种电力信号的基波幅值测量方法和系统，所述方法包括：通过对虚数向量序列和实数向量序列数字滤波，生成虚数向量滤波序列和实数向量滤波序列，分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，生成实数向量滤波前段序列、实数向量滤波后段序列、虚数向量滤波前段序列和虚数向量滤波后段序列；依次对滤波前段序列和虚数向量滤波前段序列、以及实数向量滤波后段序列和虚数向量滤波后段序列进行积分、求相位、相位差；将所述相位差和所述参考频率转换为所述电力信号的基波频率。以所述基波频率为参考频率，进行基波幅值测量。实施本发明，可抑制虚数向量序列和实数向量序列中的混频干扰成分生成高精度的基波幅值。

Description

电力信号的基波幅值测量方法和系统
【技术领域】
本发明涉及电力技术领域，特别是涉及电力信号的基波幅值测量方法和系统。
【背景技术】
电力系统的频率测量、谐波测量、功率测量等在本质上均为正弦参数的测量。傅里叶变换等是实现正弦参数测量的基本方法，在电力系统有广泛的应用。但随着正弦测量技术的发展，傅里叶变换存在的问题也越显突出，难以进一步满足正弦参数高精度计算的要求。
电力系统功率计算首先是电压电流幅值和相位的计算，而电压电流幅值和相位的计算又首先是频率的计算，可认为频率测量是正弦参数计算的基础。在电力系统频率测量方面，有形式各样的频率测量或计算方法，如零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法、基于相位差的频率算法等。
但是，电网运行额定工频为50Hz，属于较低的频率，以上所述的频率测量方法对低频信号的频率测量精度不高，且抗噪声干扰性差，易导致基波幅值的测量精度低、抗噪声干扰性差。
【发明内容】
基于此，有必要针对以上所述的频率测量方法对低频信号的频率测量精度不高，且抗噪声干扰性差，易导致基波幅值的测量精度低、抗噪声干扰性差的问题，提供一种电力信号的基波幅值测量方法和系统。
一种电力信号的基波幅值测量方法，包括以下步骤：
根据预设时间长度和预设采样频率，对电力信号进行采样获得采样数据序列；
对所述采样数据序列的基波频率进行初测，获得初步基波频率，并以初步基波频率为参考频率；
将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列；
对所述实数向量序列进行数字滤波，获得实数向量滤波序列；
将所述参考频率的正弦函数与所述采样数据序列相乘，生成虚数向量序列；
对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列；
分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，生成实数向量滤波前段序列、实数向量滤波后段序列、虚数向量滤波前段序列和虚数向量滤波后段序列；
对所述实数向量滤波前段序列和所述虚数向量滤波前段序列分别进行积分运算，生成前段序列实数积分值和前段序列虚数积分值；
根据预设的相位转换规则，将所述前段序列实数积分值与所述前段序列虚数积分值转换为第一相位；
对所述实数向量滤波后段序列和所述虚数向量滤波后段序列分别进行积分运算，生成后段序列实数积分值和后段序列虚数积分值；
根据所述预设的相位转换规则，将所述后段序列实数积分值和所述后段序列虚数积分值转换为第二相位；
将所述第二相位减去所述第一相位，生成相位差；
根据预设的频率转换规则，将所述相位差和所述参考频率转换为所述电力信号的基波频率；
以所述基波频率为参考频率，并获取与所述基波频率对应的数字滤波参数；
将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列；
以所述数字滤波参数对所述实数向量序列进行数字滤波，生成实数向量滤波序列；
对所述实数向量滤波序列进行积分运算，生成实数向量积分值；
将所述参考频率的正弦函数与所述采样数据序列相乘，获得虚数向量序列；
以所述数字滤波参数对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列；
对所述虚数向量滤波序列进行积分运算，生成虚数向量积分值；
根据预设的幅值转换规则，将所述实数向量积分值和所述虚数向量积分值转换为基波幅值。
一种电力信号的基波幅值测量系统，包括：
采样模块，用于根据预设时间长度和预设采样频率，对电力信号进行采样获得采样数据序列；
初测模块，用于对所述采样数据序列的基波频率进行初测，获得初步基波频率，并以初步基波频率为参考频率；
第一实数向量序列模块，用于将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列；
第一实数向量滤波模块，用于对所述实数向量序列进行数字滤波，获得实数向量滤波序列；
第一虚数向量序列模块，用于将所述参考频率的正弦函数与所述采样数据序列相乘，生成虚数向量序列；
第一虚数向量滤波模块，用于对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列；
序列等分模块，用于分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，生成实数向量滤波前段序列、实数向量滤波后段序列、虚数向量滤波前段序列和虚数向量滤波后段序列；
前段序列积分模块，用于对所述实数向量滤波前段序列和所述虚数向量滤波前段序列分别进行积分运算，生成前段序列实数积分值和前段序列虚数积分值；
第一相位模块，用于根据预设的相位转换规则，将所述前段序列实数积分值与所述前段序列虚数积分值转换为第一相位；
后段序列积分模块，用于对所述实数向量滤波后段序列和所述虚数向量滤波后段序列分别进行积分运算，生成后段序列实数积分值和后段序列虚数积分值；
第二相位模块，用于根据所述预设的相位转换规则，将所述后段序列实数积分值和所述后段序列虚数积分值转换为第二相位；
相位差模块，用于将所述第二相位减去所述第一相位，生成相位差；
基波频率模块，用于根据预设的频率转换规则，将所述相位差和所述参考频率转换为所述电力信号的基波频率；
参考频率重置模块，用于以所述基波频率为参考频率，并获取与所述基波频率对应的数字滤波参数；
第二实数向量序列模块，用于将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列；
第二实数向量滤波模块，用于以所述数字滤波参数对所述实数向量序列进行数字滤波，生成实数向量滤波序列；
实数向量积分模块，用于对所述实数向量滤波序列进行积分运算，生成实数向量积分值；
第二虚数向量序列模块，用于将所述参考频率的正弦函数与所述采样数据序列相乘，获得虚数向量序列；
第二虚数向量滤波模块，用于以所述数字滤波参数对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列；
虚数向量积分模块，用于对所述虚数向量滤波序列进行积分运算，生成虚数向量积分值；
基波幅值模块，用于根据预设的幅值转换规则，将所述实数向量积分值和所述虚数向量积分值转换为基波幅值。
上述电力信号的基波幅值测量方法和系统，通过对虚数向量序列和实数向量序列数字滤波，生成虚数向量滤波序列和实数向量滤波序列，数字滤波可抑制虚数向量序列和实数向量序列中的混频干扰成分，得到高精度的虚数向量滤波序列和实数向量滤波序列，分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，生成实数向量滤波前段序列、实数向量滤波后段序列、虚数向量滤波前段序列和虚数向量滤波后段序列；依次对滤波前段序列和虚数向量滤波前段序列、以及实数向量滤波后段序列和虚数向量滤波后段序列进行积分、求相位、相位差；将所述相位差和所述参考频率转换为所述电力信号的基波频率。将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，目的是根据前后两段序列的相位差，生成高精度的基波频率。以所述基波频率为参考频率，进行基波幅值测量，可进一步生成高精度的基波幅值。
【附图说明】
图1是本发明电力信号的基波幅值测量方法的流程示意图；
图2是本发明电力信号的基波幅值测量系统的结构示意图。
图3是本发明电力信号的基波幅值测量系统中基波幅值相对误差随基波频率变化示意图。
【具体实施方式】
为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步地详细描述。
本发明中的步骤虽然用标号进行了排列，但并不用于限定步骤的先后次序，除非明确说明了步骤的次序或者某步骤的执行需要其他步骤作为基础，否则步骤的相对次序是可以调整的。
请参阅图1，图1是本发明电力信号的基波幅值测量方法的流程示意图。
本实施方式的所述电力信号的基波幅值测量方法，可包括以下步骤：
步骤S101，根据预设时间长度和预设采样频率，对电力信号进行采样获得采样数据序列。
步骤S102，对所述采样数据序列的基波频率进行初测，获得初步基波频率，并以初步基波频率为参考频率。
步骤S103，将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列。
步骤S104，对所述实数向量序列进行数字滤波，获得实数向量滤波序列。
步骤S105，将所述参考频率的正弦函数与所述采样数据序列相乘，生成虚数向量序列。
步骤S106，对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列。
步骤S107，分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，生成实数向量滤波前段序列、实数向量滤波后段序列、虚数向量滤波前段序列和虚数向量滤波后段序列。
步骤S108，对所述实数向量滤波前段序列和所述虚数向量滤波前段序列分别进行积分运算，生成前段序列实数积分值和前段序列虚数积分值。
步骤S109，根据预设的相位转换规则，将所述前段序列实数积分值与所述前段序列虚数积分值转换为第一相位。
步骤S110，对所述实数向量滤波后段序列和所述虚数向量滤波后段序列分别进行积分运算，生成后段序列实数积分值和后段序列虚数积分值。
步骤S111，根据所述预设的相位转换规则，将所述后段序列实数积分值和所述后段序列虚数积分值转换为第二相位。
步骤S112，将所述第二相位减去所述第一相位，生成相位差。
步骤S113，根据预设的频率转换规则，将所述相位差和所述参考频率转换为所述电力信号的基波频率。
步骤S114，以所述基波频率为参考频率，并获取与所述基波频率对应的数字滤波参数。
步骤S115，将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列。
步骤S116，以所述数字滤波参数对所述实数向量序列进行数字滤波，生成实数向量滤波序列。
步骤S117，对所述实数向量滤波序列进行积分运算，生成实数向量积分值。
步骤S118，将所述参考频率的正弦函数与所述采样数据序列相乘，获得虚数向量序列。
步骤S119，以所述数字滤波参数对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列。
步骤S120，对所述虚数向量滤波序列进行积分运算，生成虚数向量积分值。
步骤S121，根据预设的幅值转换规则，将所述实数向量积分值和所述虚数向量积分值转换为基波幅值。
本实施方式，通过对虚数向量序列和实数向量序列数字滤波，生成虚数向量滤波序列和实数向量滤波序列，数字滤波可抑制虚数向量序列和实数向量序列中的混频干扰成分，得到高精度的虚数向量滤波序列和实数向量滤波序列，分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，生成实数向量滤波前段序列、实数向量滤波后段序列、虚数向量滤波前段序列和虚数向量滤波后段序列；依次对滤波前段序列和虚数向量滤波前段序列、以及实数向量滤波后段序列和虚数向量滤波后段序列进行积分、求相位、相位差；将所述相位差和所述参考频率转换为所述电力信号的基波频率。将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，目的是根据前后两段序列的相位差，生成高精度的基波频率。以所述基波频率为参考频率，进行基波幅值测量，可进一步生成高精度的基波幅值。
其中，对于步骤S101，所述电力信号包括正弦基波信号、正弦1/3谐波成分、正弦1/2谐波成分、正弦2次谐波成分、正弦3次谐波成分、正弦4次谐波成分和正弦5次谐波成分。优选地可通过电网领域的采样设备对所述电力信号进行采样，获得采样数据序列。
优选地，可根据在额定频率50Hz，采样频率远大于电力系统额定频率的原则设置预设数的采样频率。
进一步地，为了保证一定的频率测量实时性，可取信号时间长度等于0.25s。
更进一步地，电力系统额定频率50Hz，为了提高性能，采样频率应远大于50Hz，优选地，设置采样频率等于fn＝5000Hz，采样间隔表达为式(1)：
T n = 1 f n - - - ( 1 ) ;
式中，Tn为采样间隔，单位s；fn为所述预设采样频率，单位Hz。
在一个实施例中，采样数据序列为式(2)：
式(2)中，ω为信号基波频率，单位rad/s；Tn为采样间隔，单位s；为初相位，单位rad。
采样数据序列的基波频率与参考频率的频率差为式(3)：
Ω＝ω-ωs   (3)；
式，Ω为频率差，单位rad/s；ωs为参考频率，单位rad/s。
对于步骤S102，可通过零交法对所述采样数据序列进行频率初测，获取所述初步频率。还可通过本领域技术人员惯用的其他频率测量方法对所述采样数据序列进行频率初测。
对于步骤S103，优选地，可通过乘法器将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列。所述乘法器是一种混频器。
在一个实施例中，在不考虑混频干扰时，所述实数向量序列，或实数混频序列如式(4)所示：
式中，XR(n)为实数向量序列。
对于步骤S104，可通过数字滤波器对所述实数向量序列进行多级数字滤波，生成实数向量滤波序列，滤除混频干扰成分。
在一个实施例中，数字滤波对所述实数向量序列中NT个连续离散值相加，然后取其算术平均值作为滤波输出。
在NT取值为二分之一参考频率的单位周期序列长度时，可以对1/2分次谐波和所有次谐波影响进行抑制，而在NT取值为三分之二参考频率的单位周期序列长度时，可以对1/3分次谐波影响进行抑制。因此，数字滤波由两种滤波参数的滤波器所构成，为了提高混频干扰的抑制性能，每种滤波参数的滤波器均由参数完全相同的三级数字滤波组成，其中用NT1表达滤波参数1，用NT2表达滤波参数2。数字滤波输出序列长度相对输入信号序列长度N减小了3NT1+3NT2。具体实数向量滤波为式(15)和式(16)。
在另一个实施例中，所述(有效的)实数向量滤波序列为式(5)：
式中，XRL(n)为实数向量滤波序列，K(Ω)为数字滤波在频差Ω的增益，单位无量纲，其中K(0)＝1；β(Ω)为数字滤波在频差Ω的移相，单位rad，其中β(0)＝0。
对于步骤S105，可通过乘法器将所述参考频率的正弦函数与所述采样数据序列相乘，获得虚数向量序列。所述乘法器是一种混频器。
在一个实施例中，在不考虑混频干扰时，所述虚数向量序列，或虚数混频序列如式(6)所示：
式中，XI(n)为虚数向量序列。
对于步骤S106，可通过数字滤波器对所述虚数向量序列进行多级数字滤波，生成虚数向量滤波序列，滤除混频干扰成分。
在一个实施例中，数字滤波对所述虚数向量序列中NT个连续离散值相加，然后取其算术平均值作为滤波输出。
在NT取值为二分之一参考频率的单位周期序列长度时，可以对1/2分次谐波和所有次谐波影响进行抑制，而在NT取值为三分之二参考频率的单位周期序列长度时，可以对1/3分次谐波影响进行抑制。因此，数字滤波由两种滤波参数的滤波器所构成，为了提高混频干扰的抑制性能，每种滤波参数的滤波器均由参数完全相同的三级数字滤波组成，其中用NT3表达滤波参数3，用NT4表达滤波参数4。数字滤波输出序列长度相对输入信号序列长度N减小了3NT3+3NT4。具体虚数向量滤波为式(18)和式(19)。
在另一个实施例中，所述(有效的)虚数向量滤波序列为式(7)：
式中，XIL(n)为虚数向量滤波序列，K(Ω)为数字滤波在频差Ω的增益，单位无量纲，其中K(0)＝1；β(Ω)为数字滤波在频差Ω的移相，单位rad，其中β(0)＝0。
对于步骤S107，在分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列之前，可先测量所述实数向量滤波序列或所述虚数向量滤波序列的序列长度，然后进行序列等分。
在一个实施例中，所述实数向量滤波前段序列(或所述虚数向量前段序列)与所述实数向量滤波后段序列(或所述虚数向量后段序列长度)的序列长度分别为式(8)：
M = N - N T 1 - N T 2 2 = N - N T 3 - N T 4 2 - - - ( 8 ) ;
对于步骤S108，可通过积分器对所述实数向量滤波前段序列和所述虚数向量滤波前段序列分别进行积分运算，生成前段序列实数积分值和前段序列虚数积分值。
在一个实施例中，通过以下公式(9)对所述实数向量滤波前段序列和所述虚数向量滤波前段序列分别进行积分运算，生成前段序列实数积分值和前段序列虚数积分值：
式中，RA(ωs)为前段序列实数积分值，IA(ωs)为前段序列虚数积分值。
对于步骤S109，优选地，所述预设的相位转换规则如公式(10)所示：
其中，PHA(ωs)为所述第一相位，单位rad。
在一个实施例中，根据预设的相位转换规则，将所述前段序列实数积分值与所述前段序列虚数积分值转换为第一相位的步骤包括以下步骤：
将所述前段序列虚数积分值除以所述前段序列实数积分值，生成第一比值。
获取所述第一比值的反余切函数值的相反数，生成所述第一相位。
对于步骤S110，可通过积分器对所述实数向量滤波后段序列和所述虚数向量滤波后段序列分别进行积分运算，生成后段序列实数积分值和后段序列虚数积分值。
在一个实施例中，通过以下公式(11)对所述实数向量滤波后段序列和所述虚数向量滤波后段序列分别进行积分运算，生成后段序列实数积分值和后段序列虚数积分值：
式中，RB(ωs)为后段序列实数积分值，IB(ωs)为后段序列虚数积分值。
对于步骤S111，优选地，所述预设的相位转换规则如公式(12)所示：
其中，PHB(ωs)为所述第二相位，单位rad。
在一个实施例中，根据所述预设的相位转换规则，将所述后段序列实数积分值和所述后段序列虚数积分值转换为第二相位的步骤包括以下步骤：
将所述后段序列虚数积分值除以所述后段序列实数积分值，生成第二比值。
获取所述第二比值的反余切函数值的相反数，生成所述第二相位。
对于步骤S112，优选地，所述相位差如式(13)所示：
ΔPH(ωs)＝PHB(ωs)-PHA(ωs)＝ΩTnM   (13)；
其中，ΔPH(ωs)为所述相位差。
对于步骤S113，优选地，预设的频率转换规则如式(14)所示：
ω = ΔPH ( ω s ) T n M + ω s - - - ( 14 ) ;
其中，ω为基波频率，单位rad。
在一个实施例中，根据预设的频率转换规则，将所述相位差和所述参考频率转换为所述电力信号的基波频率的步骤包括以下步骤：
检测所述虚数向量滤波序列或所述实数向量滤波序列的序列长度。
获取所述序列长度与采样间隔的乘积。
获取所述相位差与所述乘积的比值的二倍，生成第一频率，数量上等于所述采样数据序列的基波频率与所述参考频率的频率差。
将所述第一频率与所述参考频率相加，生成所述电力信号的基波频率。
对于步骤S114，以所述基波频率为参考频率，根据所述参考频率重设数字滤波的参数，可进一步减小测量误差。
对于步骤S115，优选地，可通过乘法器将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列。所述乘法器是一种混频器。
对于步骤S116，可通过数字滤波器对所述实数向量序列进行多级数字滤波，生成实数向量滤波序列，滤除混频干扰成分。
在一个实施例中，数字滤波对所述实数向量序列中NT个连续离散值相加，然后取其算术平均值作为滤波输出。
在NT取值为二分之一参考频率的单位周期序列长度时，可以对1/2分次谐波和所有次谐波影响进行抑制，而在NT取值为三分之二参考频率的单位周期序列长度时，可以对1/3分次谐波影响进行抑制。因此，数字滤波由两种滤波参数的滤波器所构成，为了提高混频干扰的抑制性能，每种滤波参数的滤波器均由参数完全相同的三级数字滤波组成，其中用NT1表达滤波参数1，用NT2表达滤波参数2。数字滤波输出序列长度相对输入信号序列长度N减小了3NT1+3NT2。
滤波参数NT取值为二分之一参考频率的单位周期序列长度时，第一数字滤波为式(15)：
X 1 ( n ) = 1 N T 1 Σ n N T 1 - 1 { 1 N T 1 Σ n N T 1 - 1 [ 1 N T 1 Σ n N T 1 - 1 X R ( n ) ] } n = 0,1,2,3 , . . . . . . , N - 1 - - - ( 15 ) ;
式中，X1(n)为第一数字滤波输出序列，XR(n)为第一数字滤波输入序列(所述实数向量序列)，NT1为第一滤波参数,即连续离散值相加数量。
NT取值为三分之二参考频率的单位周期序列长度时，第二数字滤波为式(16)：
X 2 ( n ) = 1 N T 2 Σ n N T 2 - 1 { 1 N T 2 Σ n N T 2 - 1 [ 1 N T 2 Σ n N T 2 - 1 X 1 ( n ) ] } n = 0,1,2,3 , . . . . , N - 3 N T 1 - 1 X RL ( n ) = X 2 ( n ) n = 0,1,2,3 , . . . , N - 3 N T 1 - 3 N T 2 - 1 - - - ( 16 ) ;
式中，X2(n)为第二数字滤波输出序列，XRL(n)为实数向量滤波序列，X1(n)为所述第一数字滤波输出序列，NT2为第二滤波参数，即连续离散值相加数量。
在另一个实施例中，对所述实数向量序列进行数字滤波，生成实数向量滤波序列的步骤包括以下步骤：
通过第一滤波器对所述实数向量序列进行三级数字滤波，生成第一滤波数据序列，其中，所述第一滤波器的每级数字滤波的滤波参数均为所述参考频率单位周期序列长度的2倍。
通过第二滤波器对所述第一滤波数据序列进行三级数字滤波，生成所述实数向量滤波序列，其中，所述第二滤波器的每级数字滤波的滤波参数为所述参考频率单位周期序列长度的1.5倍。
在其他实施例中，还可通过第一滤波器对实数向量序列进行三级以上数字滤波，通过第二滤波器对所述第一滤波数据序列进行三级以上数字滤波。
对于步骤S117，可通过积分器对所述实数向量滤波序列进行积分运算，生成实数向量积分值。
在混频干扰得到完全抑制前提下，实数向量积分为式(17)
式中，Re为实数向量积分值，A为输入信号的基波幅值。
对于步骤S118，可通过乘法器将所述参考频率的正弦函数与所述采样数据序列相乘，获得虚数向量序列。所述乘法器是一种混频器。
对于步骤S119，可通过数字滤波器对所述虚数向量序列进行多级数字滤波，生成虚数向量滤波序列，滤除混频干扰成分。
在一个实施例中，数字滤波对所述虚数向量序列中NT个连续离散值相加，然后取其算术平均值作为滤波输出。
在NT取值为二分之一参考频率的单位周期序列长度时，可以对1/2分次谐波和所有次谐波影响进行抑制，而在NT取值为三分之二参考频率的单位周期序列长度时，可以对1/3分次谐波波影响进行抑制。因此，数字滤波由两种滤波参数的滤波器所构成，为了提高混频干扰的抑制性能，每种滤波参数的滤波器均由参数完全相同的三级数字滤波组成，其中用NT3表达滤波参数3，用NT4表达滤波参数4。数字滤波输出序列长度相对输入信号序列长度N减小了3NT3+3NT4。
滤波参数NT取值为二分之一参考频率的单位周期序列长度时，第三数字滤波为式(18)：
X 3 ( n ) = 1 N T 3 Σ n N T 3 - 1 { 1 N T 3 Σ n N T 3 - 1 [ 1 N T 3 Σ n N T 3 - 1 X I ( n ) ] } n = 0,1,2,3 , . . . . , N - 1 - - - ( 18 ) ;
式中，X3(n)为第三数字滤波输出序列，XI(n)为第三数字滤波输入序列(所述虚数向量序列)，NT3为第三滤波参数，即连续离散值相加数量。
NT取值为三分之二参考频率的单位周期序列长度时，第四数字滤波为式(19)：
X 4 ( n ) = 1 N T 4 Σ n N T 4 - 1 { 1 N T 4 Σ n N T 4 - 1 [ 1 N T 4 Σ n N T 4 - 1 X 3 ( n ) ] } n = 0,1,2,3 , . . . . , N - 3 N T 3 - 1 X IL ( n ) = X 4 ( n ) n = 0,1,2,3 , . . . , N - 3 N T 3 - 3 N T 4 - 1 - - - ( 19 ) ;
式中，X4(n)为第四数字滤波输出序列，XIL(n)为虚数向量滤波序列，X3(n)为所述第三数字滤波输出序列，NT4为第四滤波参数，即连续离散值相加数量。
在另一个实施例中，对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列的步骤包括以下步骤：
通过第三滤波器对所述虚数向量序列进行三级数字滤波，生成第二滤波数据序列，其中，所述第三滤波器的每级数字滤波的滤波参数均为所述参考频率单位周期序列长度的2倍。
通过第四滤波器对所述第二滤波数据序列进行三级数字滤波，生成所述虚数向量滤波序列，其中，所述第四滤波器的每级数字滤波的滤波参数为所述参考频率单位周期序列长度的1.5倍。
在其他实施例中，还可通过第三滤波器对虚数向量序列进行三级以上数字滤波，通过第四滤波器对所述第二滤波数据序列进行三级以上数字滤波。
对于步骤S120，可通过积分器对所述虚数向量滤波序列进行积分运算，生成虚数向量积分值。
在混频干扰得到完全抑制前提下，虚数向量积分为式(20)：
式中，Im为虚数向量积分值，A为输入信号的基波幅值。
对于步骤S121，优选地，所述预设的幅值转换规则可为将所述实数向量积分值和所述虚数向量积分值转换为基波幅值的转换公式。
基波幅值的转换为式(21)：
X ( ω ) = R e 2 + I m 2 = A - - - ( 21 ) ;
式中，X(ω)为基波幅值，量程和单位由输入信号的量程和单位决定；A为输入信号的基波幅值。
在一个实施例中，根据预设的幅值转换规则，将所述实数向量积分值和所述虚数向量积分值转换为基波幅值的步骤包括以下步骤：
获取所述实数向量积分值与所述虚数向量积分值的平方和。
对获取的平方和开根，生成基波幅值。
请参阅图2，图2是本发明电力信号的基波幅值测量系统的结构示意图。
本实施方式的所述电力信号的基波幅值测量系统，可包括采样模块1010、初测模块1020、第一实数向量序列模块1030、第一实数向量滤波模块1040、第一虚数向量序列模块1050、第一虚数向量滤波模块1060、序列等分模块1070、前段序列积分模块1080、第一相位模块1090、后段序列积分模块1100、第二相位模块1110、相位差模块1120、基波频率模块1130、参考频率重置模块1140、第二实数向量序列模块1150、第二实数向量滤波模块1160、第二实数向量积分模块1170、第二虚数向量序列模块1180、第二虚数向量滤波模块1190、第二虚数向量积分模块1200和基波幅值模块1210，其中：
采样模块1010，用于根据预设时间长度和预设采样频率，对电力信号进行采样获得采样数据序列。
初测模块1020，用于对所述采样数据序列的基波频率进行初测，获得初步基波频率，并以初步基波频率为参考频率。
第一实数向量序列模块1030，用于将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列。
第一实数向量滤波模块1040，用于对所述实数向量序列进行数字滤波，获得实数向量滤波序列。
第一虚数向量序列模块1050，用于将所述参考频率的正弦函数与所述采样数据序列相乘，生成虚数向量序列。
第一虚数向量滤波模块1060，用于对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列。
序列等分模块1070，用于分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，生成实数向量滤波前段序列、实数向量滤波后段序列、虚数向量滤波前段序列和虚数向量滤波后段序列。
前段序列积分模块1080，用于对所述实数向量滤波前段序列和所述虚数向量滤波前段序列分别进行积分运算，生成前段序列实数积分值和前段序列虚数积分值。
第一相位模块1090，用于根据预设的相位转换规则，将所述前段序列实数积分值与所述前段序列虚数积分值转换为第一相位。
后段序列积分模块1100，用于对所述实数向量滤波后段序列和所述虚数向量滤波后段序列分别进行积分运算，生成后段序列实数积分值和后段序列虚数积分值。
第二相位模块1110，用于根据所述预设的相位转换规则，将所述后段序列实数积分值和所述后段序列虚数积分值转换为第二相位。
相位差模块1120，用于将所述第二相位减去所述第一相位，生成相位差。
基波频率模块1130，用于根据预设的频率转换规则，将所述相位差和所述参考频率转换为所述电力信号的基波频率。
参考频率重置模块1140，用于以所述基波频率为参考频率，并获取与所述基波频率对应的数字滤波参数。
第二实数向量序列模块1150，用于将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列。
第二实数向量滤波模块1160，用于以所述数字滤波参数对所述实数向量序列进行数字滤波，生成实数向量滤波序列。
实数向量积分模块1170，用于对所述实数向量滤波序列进行积分运算，生成实数向量积分值。
第二虚数向量序列模块1180，用于将所述参考频率的正弦函数与所述采样数据序列相乘，获得虚数向量序列。
第二虚数向量滤波模块1190，用于以所述数字滤波参数对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列。
虚数向量积分模块1200，用于对所述虚数向量滤波序列进行积分运算，生成虚数向量积分值。
基波幅值模块1210，用于根据预设的幅值转换规则，将所述实数向量积分值和所述虚数向量积分值转换为基波幅值。
本实施方式，通过对虚数向量序列和实数向量序列数字滤波，生成虚数向量滤波序列和实数向量滤波序列，数字滤波可抑制虚数向量序列和实数向量序列中的混频干扰成分，得到高精度的虚数向量滤波序列和实数向量滤波序列，分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，生成实数向量滤波前段序列、实数向量滤波后段序列、虚数向量滤波前段序列和虚数向量滤波后段序列；依次对滤波前段序列和虚数向量滤波前段序列、以及实数向量滤波后段序列和虚数向量滤波后段序列进行积分、求相位、相位差；将所述相位差和所述参考频率转换为所述电力信号的基波频率。将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，目的是根据前后两段序列的相位差，生成高精度的基波频率。以所述基波频率为参考频率，进行基波幅值测量，可进一步生成高精度的基波幅值。
其中，对于采样模块1010，所述电力信号包括正弦基波信号、正弦1/3谐波成分、正弦1/2谐波成分、正弦2次谐波成分、正弦3次谐波成分、正弦4次谐波成分和正弦5次谐波成分。优选地可通过电网领域的采样设备对所述电力信号进行采样，获得采样数据序列。
优选地，可根据在额定频率50Hz，采样频率远大于电力系统额定频率的原则设置预设数的采样频率。
进一步地，为了保证一定的频率测量实时性，可取信号时间长度等于0.25s。
更进一步地，电力系统额定频率50Hz，为了提高性能，采样频率应远大于50Hz，优选地，设置采样频率等于fn＝5000Hz，采样间隔表达为式(1)：
T n = 1 f n - - - ( 1 ) ;
式中，Tn为采样间隔，单位s；fn为所述预设采样频率，单位Hz。
在一个实施例中，采样数据序列为式(2)：
式(2)中，ω为信号基波频率，单位rad/s；Tn为采样间隔，单位s；为初相位，单位rad。
采样数据序列的基波频率与参考频率的频差为式(3)：
Ω＝ω-ωs   (3)；
式，Ω为频率差，单位rad/s；ωs为参考频率，单位rad/s。
对于初测模块1020，可通过零交法对所述采样数据序列进行频率初测，获取所述初步频率。还可通过本领域技术人员惯用的其他频率测量方法对所述采样数据序列进行频率初测。
对于实数向量序列模块1030，优选地，可通过乘法器将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列。所述乘法器为一种混频器。
在一个实施例中，在不考虑混频干扰时，所述实数向量序列，或实数混频序列如式(4)所示：
式中，XR(n)为实数向量序列。
对于实数向量滤波模块1040，可通过数字滤波器对所述实数向量序列进行多级数字滤波，生成实数向量滤波序列，滤除混频干扰成分。
在一个实施例中，数字滤波对所述实数向量序列中NT个连续离散值相加，然后取其算术平均值作为滤波输出。
在NT取值为二分之一参考频率的单位周期序列长度时，可以对1/2分次谐波和所有次谐波影响进行抑制，而在NT取值为三分之二参考频率的单位周期序列长度时，可以对1/3分次谐波影响进行抑制。因此，数字滤波由两种滤波参数的滤波器所构成，为了提高混频干扰的抑制性能，每种滤波参数的滤波器均由参数完全相同的三级数字滤波组成，其中用NT1表达滤波参数1，用NT2表达滤波参数2。数字滤波输出序列长度相对输入信号序列长度N减小了3NT1+3NT2。具体实数向量滤波为式(15)和式(16)。
在另一个实施例中，所述(有效的)实数向量滤波序列为式(5)：
式中，XRL(n)为实数向量滤波序列，K(Ω)为数字滤波在频差Ω的增益，单位无量纲，其中K(0)＝1；β(Ω)为数字滤波在频差Ω的移相，单位rad，其中β(0)＝0。
对于虚数向量序列模块1050，可通过乘法器将所述参考频率的正弦函数与所述采样数据序列相乘，获得虚数向量序列。
在一个实施例中，在不考虑混频干扰时，所述虚数向量序列，或虚数混频序列如式(6)所示：
式中，XI(n)为虚数向量序列。
对于虚数向量滤波模块1060，可通过数字滤波器对所述虚数向量序列进行多级数字滤波，生成虚数向量滤波序列，滤除混频干扰成分。
在一个实施例中，数字滤波对所述虚数向量序列中NT个连续离散值相加，然后取其算术平均值作为滤波输出。
在NT取值为二分之一参考频率的单位周期序列长度时，可以对1/2分次谐波和所有次谐波影响进行抑制，而在NT取值为三分之二参考频率的单位周期序列长度时，可以对1/3分次谐波影响进行抑制。因此，数字滤波由两种滤波参数的滤波器所构成，为了提高混频干扰的抑制性能，每种滤波参数的滤波器均由参数完全相同的三级数字滤波组成，其中用NT3表达滤波参数3，用NT4表达滤波参数4。数字滤波输出序列长度相对输入信号序列长度N减小了3NT3+3NT4。具体虚数向量滤波为式(18)和式(19)。
在另一个实施例中，所述(有效的)虚数向量滤波序列为式(7)：
式中，XIL(n)为虚数向量滤波序列，K(Ω)为数字滤波在频差Ω的增益，单位无量纲，其中K(0)＝1；β(Ω)为数字滤波在频差Ω的移相，单位rad，其中β(0)＝0。
对于序列等分模块1070，在分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列之前，可先测量所述实数向量滤波序列或所述虚数向量滤波序列的序列长度，然后进行序列等分。
在一个实施例中，所述实数向量滤波前段序列(或所述虚数向量前段序列)与所述实数向量滤波后段序列(或所述虚数向量后段序列长度)的序列长度分别为式(8)：
M = N - N T 1 - N T 2 2 = N - N T 3 - N T 4 2 - - - ( 8 ) ;
对于前段序列积分模块1080，可通过积分器对所述实数向量滤波前段序列和所述虚数向量滤波前段序列分别进行积分运算，生成前段序列实数积分值和前段序列虚数积分值。
在一个实施例中，前段序列积分模块1080可通过以下公式(9)对所述实数向量滤波前段序列和所述虚数向量滤波前段序列分别进行积分运算，生成前段序列实数积分值和前段序列虚数积分值：
式中，RA(ωs)为前段序列实数积分值，IA(ωs)为前段序列虚数积分值。
对于第一相位模块1090，优选地，所述预设的相位转换规则如公式(10)所示：
其中，PHA(ωs)为所述第一相位，单位rad。
在一个实施例中，第一相位模块1090还可用于：
将所述前段序列虚数积分值除以所述前段序列实数积分值，生成第一比值。
获取所述第一比值的反余切函数值的相反数，生成所述第一相位。
对于后段序列积分模块1100，可通过积分器对所述实数向量滤波后段序列和所述虚数向量滤波后段序列分别进行积分运算，生成后段序列实数积分值和后段序列虚数积分值。
在一个实施例中，后段序列积分模块1100可通过以下公式(11)对所述实数向量滤波后段序列和所述虚数向量滤波后段序列分别进行积分运算，生成后段序列实数积分值和后段序列虚数积分值：
式中，RB(ωs)为后段序列实数积分值，IB(ωs)为后段序列虚数积分值。
对于第二相位模块1110，优选地，所述预设的相位转换规则如公式(12)所示：
其中，PHB(ωs)为所述第二相位，单位rad。
在一个实施例中，第二相位模块1110还可用于：
将所述后段序列虚数积分值除以所述后段序列实数积分值，生成第二比值。
获取所述第二比值的反余切函数值的相反数，生成所述第二相位。
对于相位差模块1120，优选地，所述相位差如式(13)所示：
ΔPH(ωs)＝PHB(ωs)-PHA(ωs)＝ΩTnM(13)；
其中，ΔPH(ωs)为所述相位差。
对于基波频率模块1130，优选地，预设的频率转换规则如式(14)所示：
ω = ΔPH ( ω s ) T n M + ω s - - - ( 14 ) ;
其中，ω为基波频率，单位rad。
在一个实施例中，基波频率模块还包括序列长度检测模块、乘积模块、第一频率模块和相加模块，其中：
所述序列长度检测模块用于检测所述虚数向量滤波序列或所述实数向量滤波序列的序列长度。
所述乘积模块用于获取所述序列长度与采样间隔的乘积。
所述第一频率模块用于获取所述相位差与所述乘积的比值的二倍，生成第一频率，数量上等于所述采样数据序列的基波频率与所述参考频率的频率差。
所述相加模块用于将所述第一频率与所述参考频率相加，生成所述电力信号的基波频率。
对于参考频率重置模块1140，以所述基波频率为参考频率，根据所述参考频率重设数字滤波的参数，可进一步减小测量误差。
对于第二实数向量序列模块1150，优选地，可通过乘法器将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列。所述乘法器为一种混频器。
对于第二实数向量滤波模块1160，可通过数字滤波器对所述实数向量序列进行多级数字滤波，生成实数向量滤波序列，滤除混频干扰成分。
在一个实施例中，数字滤波对所述实数向量序列中NT个连续离散值相加，然后取其算术平均值作为滤波输出。
在NT取值为二分之一参考频率的单位周期序列长度时，可以对1/2分次谐波和所有次谐波影响进行抑制，而在NT取值为三分之二参考频率的单位周期序列长度时，可以对1/3分次谐波影响进行抑制。因此，数字滤波由两种滤波参数的滤波器所构成，为了提高混频干扰的抑制性能，每种滤波参数的滤波器均由参数完全相同的三级数字滤波组成，其中用NT1表达滤波参数1，用NT2表达滤波参数2。数字滤波输出序列长度相对输入信号序列长度N减小了3NT1+3NT2。
滤波参数NT取值为二分之一参考频率的单位周期序列长度时，第一数字滤波为式(15)：
X 1 ( n ) = 1 N T 1 Σ n N T 1 - 1 { 1 N T 1 Σ n N T 1 - 1 [ 1 N T 1 Σ n N T 1 - 1 X R ( n ) ] } n = 0,1,2,3 , . . . . . . , N - 1 - - - ( 15 ) ;
式中，X1(n)为第一数字滤波输出序列，XR(n)为第一数字滤波输入序列(所述实数向量序列)，NT1为第一滤波参数,即连续离散值相加数量。
NT取值为三分之二参考频率的单位周期序列长度时，第二数字滤波为式(16)：
X 2 ( n ) = 1 N T 2 Σ n N T 2 - 1 { 1 N T 2 Σ n N T 2 - 1 [ 1 N T 2 Σ n N T 2 - 1 X 1 ( n ) ] } n = 0,1,2,3 , . . . . , N - 3 N T 1 - 1 X RL ( n ) = X 2 ( n ) n = 0,1,2,3 , . . . , N - 3 N T 1 - 3 N T 2 - 1 - - - ( 16 ) ;
式中，X2(n)为第二数字滤波输出序列，XRL(n)为实数向量滤波序列，X1(n)为所述第一数字滤波输出序列，NT2为第二滤波参数，即连续离散值相加数量。
在另一个实施例中，第二实数向量滤波模块1160可用于：
通过第一滤波器对所述实数向量序列进行三级数字滤波，生成第一滤波数据序列，其中，所述第一滤波器的每级数字滤波的滤波参数均为所述参考频率单位周期序列长度的2倍。
通过第二滤波器对所述第一滤波数据序列进行三级数字滤波，生成所述实数向量滤波序列，其中，所述第二滤波器的每级数字滤波的滤波参数为所述参考频率单位周期序列长度的1.5倍。
在其他实施例中，还可通过第一滤波器对实数向量序列进行三级以上数字滤波，通过第二滤波器对所述第一滤波数据序列进行三级以上数字滤波。
对于第二实数向量积分模块1170，可通过数字滤波器对所述实数向量序列进行多级数字滤波，生成实数向量滤波序列，可通过积分器对所述实数向量滤波序列进行积分运算，生成实数向量积分值。
在混频干扰得到完全抑制前提下，实数向量积分为式(17)
式中，Re为实数向量积分值，A为输入信号的基波幅值。
对于第二虚数向量序列模块1180，可通过乘法器将所述参考频率的正弦函数与所述采样数据序列相乘，获得虚数向量序列。
对于第二虚数向量滤波模块1190，可通过数字滤波器对所述虚数向量序列进行多级数字滤波，生成虚数向量滤波序列，滤除混频干扰成分。
在一个实施例中，数字滤波对所述虚数向量序列中NT个连续离散值相加，然后取其算术平均值作为滤波输出。
在NT取值为二分之一参考频率的单位周期序列长度时，可以对1/2分次谐波和所有次谐波影响进行抑制，而在NT取值为三分之二参考频率的单位周期序列长度时，可以对1/3分次谐波影响进行抑制。因此，数字滤波由两种滤波参数的滤波器所构成，为了提高混频干扰的抑制性能，每种滤波参数的滤波器均由参数完全相同的三级数字滤波组成，其中用NT3表达滤波参数3，用NT4表达滤波参数4。数字滤波输出序列长度相对输入信号序列长度N减小了3NT3+3NT4。
滤波参数NT取值为二分之一参考频率的单位周期序列长度时，第三数字滤波为式(18)：
X 3 ( n ) = 1 N T 3 Σ n N T 3 - 1 { 1 N T 3 Σ n N T 3 - 1 [ 1 N T 3 Σ n N T 3 - 1 X I ( n ) ] } n = 0,1,2,3 , . . . . , N - 1 - - - ( 18 ) ;
式中，X3(n)为第三数字滤波输出序列，XI(n)为第三数字滤波输入序列(所述虚数向量序列)，NT3为第三滤波参数,即连续离散值相加数量。
NT取值为三分之一参考频率单位周期序列长度时，第四数字滤波为式(19)：
X 4 ( n ) = 1 N T 4 Σ n N T 4 - 1 { 1 N T 4 Σ n N T 4 - 1 [ 1 N T 4 Σ n N T 4 - 1 X 3 ( n ) ] } n = 0,1,2,3 , . . . . , N - 3 N T 3 - 1 X IL ( n ) = X 4 ( n ) n = 0,1,2,3 , . . . , N - 3 N T 3 - 3 N T 4 - 1 - - - ( 19 ) ;
式中，X4(n)为第四数字滤波输出序列，XIL(n)为虚数向量滤波序列，X3(n)为所述第三数字滤波输出序列，NT4为第四滤波参数，即连续离散值相加数量。
在另一个实施例中，第二虚数向量滤波模块1190可用于：
通过第三滤波器对所述虚数向量序列进行三级数字滤波，生成第二滤波数据序列，其中，所述第三滤波器的每级数字滤波的滤波参数均为所述参考频率单位周期序列长度的2倍。
通过第四滤波器对所述第二滤波数据序列进行三级数字滤波，生成所述虚数向量滤波序列，其中，所述第四滤波器的每级数字滤波的滤波参数为所述参考频率单位周期序列长度的1.5倍。
在其他实施例中，还可通过第三滤波器对虚数向量序列进行三级以上数字滤波，通过第四滤波器对所述第二滤波数据序列进行三级以上数字滤波。
对于第二虚数向量积分模块1200，可通过积分器对所述虚数向量滤波序列进行积分运算，生成虚数向量积分值。
在混频干扰得到完全抑制前提下，虚数向量积分为式(20)：
式中，Im为虚数向量积分值，A为输入信号的基波幅值。
对于基波幅值模块1210，优选地，所述预设的幅值转换规则可为将所述实数向量积分值和所述虚数向量积分值转换为基波幅值的转换公式。
基波幅值的转换为式(21)：
X ( ω ) = R e 2 + I m 2 = A - - - ( 21 ) ;
式中，X(ω)为基波幅值，量程和单位由输入信号的量程和单位决定；A为输入信号的基波幅值。
在一个实施例中，基波幅值模块1210可用于：
获取所述实数向量积分值与所述虚数向量积分值的平方和。
对获取的平方和开根，生成基波幅值。
请参阅图3，图3是本发明电力信号的基波幅值测量系统中基波幅值相对误差随基波频率变化示意图。
对50Hz工频条件下的电力信号的基波幅值测量方法进行仿真实验，仿真实验条件为：信号基波频率变化范围45Hz-55Hz，信号采样频率10KHz，信号离散数据量化位数24bit，频率初测单元所测初步频率相对误差±0.25％。
在基波频率45Hz-55Hz变化、信号窗口时间0.25s、参考频率误差0.25％，得到基波幅值相对误差随基波频率变化的实验结果如图3所示。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种电力信号的基波幅值测量方法，其特征在于，包括以下步骤：
根据预设时间长度和预设采样频率，对电力信号进行采样获得采样数据序列；
对所述采样数据序列的基波频率进行初测，获得初步基波频率，并以初步基波频率为参考频率；
将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列；
对所述实数向量序列进行数字滤波，获得实数向量滤波序列；
将所述参考频率的正弦函数与所述采样数据序列相乘，生成虚数向量序列；
对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列；
分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，生成实数向量滤波前段序列、实数向量滤波后段序列、虚数向量滤波前段序列和虚数向量滤波后段序列；
对所述实数向量滤波前段序列和所述虚数向量滤波前段序列分别进行积分运算，生成前段序列实数积分值和前段序列虚数积分值；
根据预设的相位转换规则，将所述前段序列实数积分值与所述前段序列虚数积分值转换为第一相位；
对所述实数向量滤波后段序列和所述虚数向量滤波后段序列分别进行积分运算，生成后段序列实数积分值和后段序列虚数积分值；
根据所述预设的相位转换规则，将所述后段序列实数积分值和所述后段序列虚数积分值转换为第二相位；
将所述第二相位减去所述第一相位，生成相位差；
根据预设的频率转换规则，将所述相位差和所述参考频率转换为所述电力信号的基波频率；
以所述基波频率为参考频率，并获取与所述基波频率对应的数字滤波参数；
将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列；
以所述数字滤波参数对所述实数向量序列进行数字滤波，生成实数向量滤波序列；
对所述实数向量滤波序列进行积分运算，生成实数向量积分值；
将所述参考频率的正弦函数与所述采样数据序列相乘，获得虚数向量序列；
以所述数字滤波参数对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列；
对所述虚数向量滤波序列进行积分运算，生成虚数向量积分值；
根据预设的幅值转换规则，将所述实数向量积分值和所述虚数向量积分值转换为基波幅值。
2.根据权利要求1所述的电力信号的基波幅值测量方法，其特征在于，所述电力信号包括正弦基波信号、正弦1/3谐波成分、正弦1/2谐波成分、正弦2次谐波成分、正弦3次谐波成分、正弦4次谐波成分和正弦5次谐波成分。
3.根据权利要求1所述的电力信号的基波幅值测量方法，其特征在于，根据预设的相位转换规则，将所述前段序列实数积分值与所述前段序列虚数积分值转换为第一相位的步骤包括以下步骤：
将所述前段序列虚数积分值除以所述前段序列实数积分值，生成第一比值；
获取所述第一比值的反余切函数值的相反数，生成所述第一相位。
4.根据权利要求1所述的电力信号的基波幅值测量方法，其特征在于，根据所述预设的相位转换规则，将所述后段序列实数积分值和所述后段序列虚数积分值转换为第二相位的步骤包括以下步骤：
将所述后段序列虚数积分值除以所述后段序列实数积分值，生成第二比值；
获取所述第二比值的反余切函数值的相反数，生成所述第二相位。
5.根据权利要求1所述的电力信号的基波幅值测量方法，其特征在于，根据预设的幅值转换规则，将所述实数向量积分值和所述虚数向量积分值转换为基波幅值的步骤包括以下步骤：
获取所述实数向量积分值与所述虚数向量积分值的平方和；
对获取的平方和开根，生成基波幅值。
6.根据权利要求1至5中任意一项所述的电力信号的基波幅值测量方法，其特征在于，根据预设的频率转换规则，将所述相位差和所述参考频率转换为所述电力信号的基波频率的步骤包括以下步骤：
检测所述虚数向量滤波序列或所述实数向量滤波序列的序列长度；
获取所述序列长度与采样间隔的乘积；
获取所述相位差与所述乘积的比值的二倍，生成第一频率，数量上等于所述采样数据序列的基波频率与所述参考频率的频率差；
将所述第一频率与所述参考频率相加，生成所述电力信号的基波频率。
7.一种电力信号的基波幅值测量系统，其特征在于，包括：
采样模块，用于根据预设时间长度和预设采样频率，对电力信号进行采样获得采样数据序列；
初测模块，用于对所述采样数据序列的基波频率进行初测，获得初步基波频率，并以初步基波频率为参考频率；
第一实数向量序列模块，用于将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列；
第一实数向量滤波模块，用于对所述实数向量序列进行数字滤波，获得实数向量滤波序列；
第一虚数向量序列模块，用于将所述参考频率的正弦函数与所述采样数据序列相乘，生成虚数向量序列；
第一虚数向量滤波模块，用于对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列；
序列等分模块，用于分别将所述实数向量滤波序列和所述虚数向量滤波序列等分为两段序列，生成实数向量滤波前段序列、实数向量滤波后段序列、虚数向量滤波前段序列和虚数向量滤波后段序列；
前段序列积分模块，用于对所述实数向量滤波前段序列和所述虚数向量滤波前段序列分别进行积分运算，生成前段序列实数积分值和前段序列虚数积分值；
第一相位模块，用于根据预设的相位转换规则，将所述前段序列实数积分值与所述前段序列虚数积分值转换为第一相位；
后段序列积分模块，用于对所述实数向量滤波后段序列和所述虚数向量滤波后段序列分别进行积分运算，生成后段序列实数积分值和后段序列虚数积分值；
第二相位模块，用于根据所述预设的相位转换规则，将所述后段序列实数积分值和所述后段序列虚数积分值转换为第二相位；
相位差模块，用于将所述第二相位减去所述第一相位，生成相位差；
基波频率模块，用于根据预设的频率转换规则，将所述相位差和所述参考频率转换为所述电力信号的基波频率；
参考频率重置模块，用于以所述基波频率为参考频率，并获取与所述基波频率对应的数字滤波参数；
第二实数向量序列模块，用于将所述参考频率的余弦函数与所述采样数据序列相乘，生成实数向量序列；
第二实数向量滤波模块，用于以所述数字滤波参数对所述实数向量序列进行数字滤波，生成实数向量滤波序列；
实数向量积分模块，用于对所述实数向量滤波序列进行积分运算，生成实数向量积分值；
第二虚数向量序列模块，用于将所述参考频率的正弦函数与所述采样数据序列相乘，获得虚数向量序列；
第二虚数向量滤波模块，用于以所述数字滤波参数对所述虚数向量序列进行数字滤波，生成虚数向量滤波序列；
虚数向量积分模块，用于对所述虚数向量滤波序列进行积分运算，生成虚数向量积分值；
基波幅值模块，用于根据预设的幅值转换规则，将所述实数向量积分值和所述虚数向量积分值转换为基波幅值。
8.根据权利要求7所述的电力信号的基波幅值测量系统，其特征在于，所述第一相位模块还用于将所述前段序列虚数积分值除以所述前段序列实数积分值，生成第一比值；获取所述第一比值的反余切函数值的相反数，生成所述第一相位。
9.根据权利要求7所述的电力信号的基波幅值测量系统，其特征在于，所述第二相位模块还用于将所述后段序列虚数积分值除以所述后段序列实数积分值，生成第二比值；获取所述第二比值的反余切函数值的相反数，生成所述第二相位。
10.根据权利要求7至9中任意一项所述的电力信号的基波幅值测量系统，其特征在于，所述基波频率模块还包括序列长度检测模块、乘积模块、第一频率模块和相加模块，其中：
所述序列长度检测模块用于检测所述虚数向量滤波序列或所述实数向量滤波序列的序列长度；
所述乘积模块用于获取所述序列长度与采样间隔的乘积；
所述第一频率模块用于获取所述相位差与所述乘积的比值的二倍，生成第一频率，数量上等于所述采样数据序列的基波频率与所述参考频率的频率差；
所述相加模块用于将所述第一频率与所述参考频率相加，生成所述电力信号的基波频率。
CN201410848173.9A 2014-12-29 2014-12-29 电力信号的基波幅值测量方法和系统 Active CN104502675B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410848173.9A CN104502675B (zh) 2014-12-29 2014-12-29 电力信号的基波幅值测量方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410848173.9A CN104502675B (zh) 2014-12-29 2014-12-29 电力信号的基波幅值测量方法和系统

Publications (2)

Publication Number Publication Date
CN104502675A true CN104502675A (zh) 2015-04-08
CN104502675B CN104502675B (zh) 2017-05-24

Family

ID=52944086

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410848173.9A Active CN104502675B (zh) 2014-12-29 2014-12-29 电力信号的基波幅值测量方法和系统

Country Status (1)

Country Link
CN (1) CN104502675B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104991104A (zh) * 2015-07-09 2015-10-21 广东电网有限责任公司电力科学研究院 电力信号的幅值检测方法和系统
CN105092970A (zh) * 2015-09-18 2015-11-25 广东电网有限责任公司电力科学研究院 获取电力信号序列正弦函数零初相位基准点的方法和系统
CN105116202A (zh) * 2015-07-14 2015-12-02 云南电网有限责任公司玉溪供电局 一种快速获取电信号基波幅值和多次谐波幅值的方法
CN105203840A (zh) * 2015-09-18 2015-12-30 广东电网有限责任公司电力科学研究院 对电力信号序列进行零初相位正弦函数调制的方法和系统
CN105372471A (zh) * 2015-07-09 2016-03-02 深圳市科润宝实业有限公司 正弦信号的幅值检测方法和系统
CN113514686A (zh) * 2021-07-13 2021-10-19 北京英博电气股份有限公司 电压基波幅值的检测方法、装置、设备和存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1493002A (zh) * 2002-03-25 2004-04-28 三菱电机株式会社 电子式电能计和功率关联量运算电路
CN102435844A (zh) * 2011-11-01 2012-05-02 南京磐能电力科技股份有限公司 一种频率无关的正弦信号相量计算方法
CN103185837A (zh) * 2013-03-25 2013-07-03 华中科技大学 一种电力系统频率测量的方法
CN103353550A (zh) * 2013-04-24 2013-10-16 武汉大学 一种测量电力系统信号频率及谐波参数的方法
CN103412171A (zh) * 2013-07-12 2013-11-27 江苏大学 一种基于极端学习机的电网谐波电压信号检测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1493002A (zh) * 2002-03-25 2004-04-28 三菱电机株式会社 电子式电能计和功率关联量运算电路
CN102435844A (zh) * 2011-11-01 2012-05-02 南京磐能电力科技股份有限公司 一种频率无关的正弦信号相量计算方法
CN103185837A (zh) * 2013-03-25 2013-07-03 华中科技大学 一种电力系统频率测量的方法
CN103353550A (zh) * 2013-04-24 2013-10-16 武汉大学 一种测量电力系统信号频率及谐波参数的方法
CN103412171A (zh) * 2013-07-12 2013-11-27 江苏大学 一种基于极端学习机的电网谐波电压信号检测方法

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104991104A (zh) * 2015-07-09 2015-10-21 广东电网有限责任公司电力科学研究院 电力信号的幅值检测方法和系统
CN105372471A (zh) * 2015-07-09 2016-03-02 深圳市科润宝实业有限公司 正弦信号的幅值检测方法和系统
CN104991104B (zh) * 2015-07-09 2017-09-08 广东电网有限责任公司电力科学研究院 电力信号的幅值检测方法和系统
CN105372471B (zh) * 2015-07-09 2018-05-25 深圳市科润宝实业有限公司 正弦信号的幅值检测方法和系统
CN105116202A (zh) * 2015-07-14 2015-12-02 云南电网有限责任公司玉溪供电局 一种快速获取电信号基波幅值和多次谐波幅值的方法
CN105092970A (zh) * 2015-09-18 2015-11-25 广东电网有限责任公司电力科学研究院 获取电力信号序列正弦函数零初相位基准点的方法和系统
CN105203840A (zh) * 2015-09-18 2015-12-30 广东电网有限责任公司电力科学研究院 对电力信号序列进行零初相位正弦函数调制的方法和系统
CN105203840B (zh) * 2015-09-18 2017-11-14 广东电网有限责任公司电力科学研究院 对电力信号序列进行零初相位正弦函数调制的方法和系统
CN105092970B (zh) * 2015-09-18 2018-01-30 广东电网有限责任公司电力科学研究院 获取电力信号序列正弦函数零初相位基准点的方法和系统
CN113514686A (zh) * 2021-07-13 2021-10-19 北京英博电气股份有限公司 电压基波幅值的检测方法、装置、设备和存储介质
CN113514686B (zh) * 2021-07-13 2024-02-06 北京英博电气股份有限公司 电压基波幅值的检测方法、装置、设备和存储介质

Also Published As

Publication number Publication date
CN104502675B (zh) 2017-05-24

Similar Documents

Publication Publication Date Title
CN104502706A (zh) 电力信号的谐波幅值测量方法和系统
CN104535836A (zh) 电力信号的基波频率测量方法和系统
CN104502700A (zh) 电力信号的正弦参数测量方法和系统
CN104635044A (zh) 基于幅值调制的电力信号频率检测方法和系统
CN104635045A (zh) 基于相位调制的电力信号频率检测方法和系统
CN104459321A (zh) 电力信号的基波相位测量方法和系统
CN104459320A (zh) 电力信号的谐波相位测量方法和系统
CN104502675A (zh) 电力信号的基波幅值测量方法和系统
CN105067880A (zh) 对电力信号进行正交调制的方法和系统
CN104502698B (zh) 电力信号的频率测量方法和系统
Salcic et al. An improved Taylor method for frequency measurement in power systems
CN105372471A (zh) 正弦信号的幅值检测方法和系统
CN104991104B (zh) 电力信号的幅值检测方法和系统
CN105182077B (zh) 根据余弦函数调制的电力信号全相位差检测方法和系统
CN105445547A (zh) 正弦信号的频率检测方法和系统
CN105445551A (zh) 正弦信号的截止相位检测方法和系统
CN104977467A (zh) 电力信号的初相位检测方法和系统
CN105445552A (zh) 正弦信号的初相位检测方法和系统
CN105092967A (zh) 电力信号的频率检测方法和系统
CN105372490B (zh) 基于零初相位基准正交倍频序列的频率测量的方法和系统
CN105548692A (zh) 根据零初相位基准正交序列进行频率测量的方法和系统
CN105301355B (zh) 根据正弦函数调制的电力信号频率检测方法和系统
CN105044462A (zh) 电力信号的截止相位检测方法和系统
CN105425032B (zh) 从电力信号中获取正弦函数倍频序列的方法和系统
CN105158564A (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
GR01 Patent grant
GR01 Patent grant