CN105372471B - 正弦信号的幅值检测方法和系统 - Google Patents
正弦信号的幅值检测方法和系统 Download PDFInfo
- Publication number
- CN105372471B CN105372471B CN201510400897.1A CN201510400897A CN105372471B CN 105372471 B CN105372471 B CN 105372471B CN 201510400897 A CN201510400897 A CN 201510400897A CN 105372471 B CN105372471 B CN 105372471B
- Authority
- CN
- China
- Prior art keywords
- frequency
- sequence
- phase
- preset
- amplitude
- 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
Landscapes
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
本发明涉及一种正弦信号的幅值检测方法和系统,对采样所得的信号序列进行截短处理,获得截短信号序列;以所测参考频率的余弦函数和参考频率的正弦函数分别与信号序列和截短信号序列相乘,生成两组实频向量序列和虚频向量序列;通过对两组虚频向量序列和实频向量序列数字陷波,生成两组虚数向量陷波序列和实数向量陷波序列,进而积分生成两组虚数向量积分值和实数向量积分值;将两组实数向量积分值和虚数向量积分值转换为序列幅值和两个相位,再将两个相位进行扩展,获得扩展相位;将两个扩展相位转换为频差;将所述序列幅值、所述频差转换为所述正弦信号的基波幅值。实施本发明,得到准确度较高的基波幅值。
Description
技术领域
本发明涉及电力技术领域,特别是涉及一种正弦信号的幅值检测方法和系统。
背景技术
电力系统的频率测量、相位测量、幅值测量等在本质上均为正弦信号参数的测量。傅里叶变换是实现正弦信号参数测量的基本方法,在电力系统中有广泛的应用。但随着正弦参数测量技术的发展,傅里叶变换存在的问题也越显突出,其难以进一步满足电力系统对正弦参数高准确度计算的要求。
在电力系统正弦信号参数测量方面,还有一些改进的参数测量方法,如零交法、基于滤波的测量法、基于小波变换法、基于神经网络的测量法、基于DFT变换的测量法等。但是,以上所述方法对低频信号(比如,电网运行额定工频在50Hz附近)的测量精度低,且抗谐波和噪声干扰性差。
发明内容
本发明所要解决的技术问题在于,提供一种一种正弦信号的幅值检测方法,其对电力系统的低频正弦信号的幅值测量精度高,且抗谐波和噪声干扰性好。
为解决上述技术问题,本发明采用如下技术方案:
一种正弦信号的幅值检测方法,包括以下步骤:
根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度;
根据所述初步采样序列长度,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列;
对所述初步采样序列进行频率初测,生成所述正弦信号的初步频率,以所述初步频率给定参考频率;
根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度;
根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度;
根据预设序列长度,从所述初步采样序列中,获得信号序列;
对所述信号序列进行截短处理,获得截短信号序列;
以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,生成第一实频向量序列和第一虚频向量序列;
以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,生成第二实频向量序列和第二虚频向量序列;
分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波,生成第一实频向量陷波序列和第一虚频向量陷波序列;
分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算,生成第一实频向量积分值和第一虚频向量积分值;
分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波,生成第二实频向量陷波序列和第二虚频向量陷波序列;
分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算,生成第二实频向量积分值和第二虚频向量积分值;
根据预设的相位和幅值转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位和序列幅值;
根据预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位;
根据预设的相位扩展规则,对所述第一相位进行扩展,获得第一扩展相位;
根据预设的相位扩展规则,对所述第二相位进行扩展,获得第二扩展相位;
根据预设的频差转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号频率与所述参考频率的频差;
根据预设的幅值转换规则,将所述序列幅值、所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的基波幅值。
一种正弦信号的幅值检测系统,包括:
初步序列长度模块,用于根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度;
初步采样模块,用于根据所述初步采样序列长度,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列;
频率初测模块,用于对所述初步采样序列进行频率初测,生成所述正弦信号的初步频率,以所述初步频率给定参考频率;
周期序列模块,用于根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度;
序列长度模块,根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度;
信号序列模块,用于根据所述预设序列长度,从所述初步采样序列中,获得信号序列;
截短序列模块,用于对所述信号序列进行截短处理,获得截短信号序列,;
第一混频模块,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,生成第一实频向量序列和第一虚频向量序列;
第二混频模块,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,生成第二实频向量序列和第二虚频向量序列;
第一陷波模块,用于分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波,生成第一实频向量陷波序列和第一虚频向量陷波序列;
第一积分模块,用于分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算,生成第一实频向量积分值和第一虚频向量积分值;
第二陷波模块,用于分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波,生成第二实频向量陷波序列和第二虚频向量陷波序列;
第二积分模块,用于分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算,生成第二实频向量积分值和第二虚频向量积分值;
第一相位模块,用于根据预设的相位和幅值转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位和序列幅值;
第二相位模块,用于根据所述预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位;
第一相位扩展模块,用于根据预设的相位扩展规则,将所述第一相位扩展为第一扩展相位;
第二相位扩展模块,用于根据预设的相位扩展规则,将所述第二相位扩展为第二扩展相位;
第一相位扩展模块,用于根据预设的相位扩展规则,将所述第一相位扩展为第一扩展相位;
第二相位扩展模块,用于根据预设的相位扩展规则,将所述第二相位扩展为第二扩展相位;
频差转换模块,用于根据预设的频差转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号频率与所述参考频率的频差;
幅值转换模块,用于根据预设的幅值转换规则,将所述序列幅值、所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的基波幅值。
本发明与现有技术相比的有益效果是:上述正弦信号的幅值检测方法和系统,对采样所得的信号序列进行截短处理,获得截短信号序列;以所测参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列和所述截短信号序列相乘,生成两组实频向量序列和虚频向量序列;通过对两组虚频向量序列和实频向量序列数字陷波,生成两组虚数向量陷波序列和实数向量陷波序列,进而积分生成两组虚数向量积分值和实数向量积分值;再根据预设的相位转换规则,将两组实数向量积分值和虚数向量积分值转换为序列幅值和两个相位;再将两个相位进行扩展,获得扩展相位;根据预设的频差转换规则,将两个扩展相位转换为所述正弦信号频率与所述参考频率的频差;根据预设的幅值转换规则,将所述序列幅值、所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的基波幅值。
由于数字陷波可以在某一个频率点迅速衰减输入信号,以达到阻碍此频率信号通过的效果,因此,当该数字陷波的陷波频率点设为对应混频干扰频率点时,该数字陷波对所述混频干扰频率具有深度的抑制作用。如此,采用本发明的正弦信号的幅值检测方法和设备,其对低频正弦信号的幅值测量精度高,且抗谐波和噪声干扰性好。
附图说明
图1是本发明正弦信号的幅值检测方法在一些实施方式中的流程示意图。
图2是本发明正弦信号的幅值检测系统在一些实施方式中的结构示意图。
图3是本发明正弦信号的幅值检测方法的信号序列与截短信号序列示意图。
图4是本发明正弦信号的幅值检测系统的基波幅值检测相对误差实验结果示意图。
图5是本发明正弦信号的幅值检测系统的三角窗函数1外形与三角窗算术平均陷波器1的频域特性示意图。
图6是本发明正弦信号的幅值检测系统的三角窗函数2外形与三角窗算术平均陷波器2的频域特性示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
本发明中的步骤虽然用标号进行了排列,但并不用于限定步骤的先后次序,除非明确说明了步骤的次序或者某步骤的执行需要其他步骤作为基础,否则步骤的相对次序是可以调整的。
在一些实施例中,如图1所示,正弦信号的幅值检测方法包括以下步骤:
S101根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度;
S102根据所述初步采样序列长度,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列;
S103对所述初步采样序列进行频率初测,生成所述正弦信号的初步频率,以所述初步频率给定参考频率;
S104根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度;
S105根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度;
S106根据预设序列长度,从所述初步采样序列中,获得信号序列;
S107对所述信号序列进行截短处理,获得截短信号序列;
S108以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,生成第一实频向量序列和第一虚频向量序列;
S109以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,生成第二实频向量序列和第二虚频向量序列;
S110分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波,生成第一实频向量陷波序列和第一虚频向量陷波序列;
S111分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算,生成第一实频向量积分值和第一虚频向量积分值;
S112分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波,生成第二实频向量陷波序列和第二虚频向量陷波序列;
S113分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算,生成第二实频向量积分值和第二虚频向量积分值;
S114根据预设的相位和幅值转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位和序列幅值;
S115根据预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位;
S116根据预设的相位扩展规则,对所述第一相位进行扩展,获得第一扩展相位;
S117根据预设的相位扩展规则,对所述第二相位进行扩展,获得第二扩展相位;
S118根据预设的频差转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号频率与所述参考频率的频差;
S119根据预设的幅值转换规则,将所述序列幅值、所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的基波幅值。
本实施方式,对采样所得的信号序列进行截短处理,获得截短信号序列,;以所测参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列和所述截短信号序列相乘,生成两组实频向量序列和虚频向量序列;通过对两组虚频向量序列和实频向量序列数字陷波,生成两组虚数向量陷波序列和实数向量陷波序列,进而积分生成两组虚数向量积分值和实数向量积分值;再根据预设的相位转换规则,将两组实数向量积分值和虚数向量积分值转换为序列幅值和两个相位;再将两个相位进行扩展,获得扩展相位;根据预设的频差转换规则,将两个扩展相位转换为所述正弦信号频率与所述参考频率的频差;根据预设的幅值转换规则,将所述序列幅值、所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的基波幅值。
由于数字陷波可以在某一个频率点迅速衰减输入信号,以达到阻碍此频率信号通过的效果,因此,当该数字陷波的陷波频率点设为对应混频干扰频率点时,该数字陷波对所述混频干扰频率具有深度的抑制作用。如此,采用本发明的正弦信号的幅值检测方法和设备,其对低频正弦信号的幅值测量精度高,且抗谐波和噪声干扰性好。
之后如不加说明,所述正弦信号的幅值均指正弦信号的基波幅值。
其中,对于步骤S101,电力系统频率范围在45Hz-55Hz,优选的,取正弦信号频率下限fmin为45Hz;根据实际需要设置所述预设信号周期数C2π,优选的,所述预设信号周期数C2π取整数。
在一些实施例中,取C2π为整数12。
初步采样序列长度计算,为式(1):
式中,Nstart为所述初步采样序列长度,单位无量纲;(int)表示取整;C2π为所述预设信号周期数,单位无量纲;fmin为正弦信号频率范围的下限,单位Hz;fn为所述预设采样频率,单位Hz。
对于步骤S102,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列。对单基波频率正弦信号,获得所述正弦信号的初步采样序列为式(2):
其中,Xstart(n)为初步采样序列;A为信号幅值,单位v;ω为信号频率,单位rad/s;Tn为采样间隔,单位s;fn为所述采样频率,单位Hz;n为序列离散数,单位无量纲;为信号初相位,单位rad;Nstart为初步采样序列长度,单位无量纲。
对于步骤S103,可通过零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法对所述初步采样序列进行频率初测,获取所述初步频率。
所述初步频率表达为式(3):
ωo (3);
其中,ωo为初步频率,单位rad/s;
优选地,所述参考频率表达为式(4):
ωs=ωo (4);
其中,ωs为参考频率,单位rad/s;ωo为初步频率,单位rad/s。
对于步骤S104,根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度:
在一些实施例中,所述正弦信号的单位周期序列长度计算,为式(5):
式中,N2π为所述单位周期序列长度,单位无量纲;(int)为取整数;fn为预设采样频率,单位Hz;fs为Hz单位的参考频率;ωs为rad/s单位的参考频率。
所述单位周期序列长度整数化存在1个采样间隔内的误差。
对于步骤S105,根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度:
在一些实施例中,所述预设序列长度为所述单位周期序列长度的12倍,所述预设序列长度计算,为式(6):
N=(int)(C2πN2π) (6);
其中,N为预设序列长度,单位无量纲;(int)为取整数;N2π为所述单位周期序列长度,单位无量纲。由于存在误差,所述预设序列长度所包含的信号周期整数是大约的。
对于步骤S106,根据预设序列长度,从所述初步采样序列中,获得信号序列:
在一些实施例中,所述预设序列长度的信号序列,为式(7):
其中,Xi(n)为信号序列;Xstart(n)为初步采样序列;A为信号幅值,单位v;ω为信号频率,单位rad/s;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;为信号初相位,单位rad;N信号序列长度,单位无量纲,信号序列长度等于所述预设序列长度;Nstart为初步采样序列长度,单位无量纲。
所述信号序列的图形表达,如图3所示。
对于步骤S107,对所述信号序列进行截短处理,获得截短信号序列:
在一些实施例中,所述截短信号序列长度为所述单位周期序列长度的11.25倍,所述截短信号序列长度计算为式(8):
NS=N-0.75N2π (8);
式中,NS为截短信号序列长度,单位无量纲;N为信号序列长度,单位无量纲;N2π为信号单位周期序列长度,单位无量纲。系数0.75代表所述截短值为所述单位周期序列长度的0.75倍;
优选地,所述截短信号序列表达为式(9):
式中,X2(n)为截短信号序列;NS为截短信号序列长度,单位无量纲。
所述截短信号序列的图形表达,如图3所示。
对于步骤S108,优选地,所述参考频率的余弦函数和所述参考频率的正弦函数可分别为以所述参考频率为频率、以Tn为间隔离散变量的正弦函数和余弦函数。
在一些实施例中,以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,得到所述第一实频向量序列和所述第一虚频向量序列为式(10):
其中,R1(n)为所述第一实频向量序列;I1(n)为所述第一虚频向量序列;Ω为信号频率与参考频率的频差,单位rad/s;和为有效分量;和为混频干扰频率成分。
对于步骤S109,以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,得到所述第二实频向量序列和所述第二虚频向量序列为式(11):
式中,R2(n)为所述第二实频向量序列;I2(n)为所述第二虚频向量序列;Ω为信号频率与参考频率的频差,单位rad/s;和为有效分量;和为混频干扰频率。
对于步骤S110,所述第一实频向量序列和所述第一虚频向量序列中包含混频干扰频率。当输入信号中包含直流成分、次谐波成分及分次谐波成分时,所述混频干扰频率将更加复杂,这些混频干扰频率严重影响计算准确度。虽然窗口函数和积分运算本身对混频干扰频率具有良好的衰减作用,但没有针对性,不能够对所述复杂的混频干扰频率产生深度的抑制作用,不能满足参数的高准确度计算需要。为了有针对性的抑制所述混频干扰频率的影响,可以对该混频干扰进行数字陷波器,由于数字陷波可以在某一个频率点迅速衰减输入信号,以达到阻碍此频率信号通过的效果,因此,当该数字陷波的陷波频率点设为对应混频干扰频率点时,该数字陷波对所述混频干扰频率具有完全的抑制作用。
优选地,数字陷波具体采用滑动三角窗算术平均陷波,即将若干个连续离散值乘以三角窗函数后再相加,然后取其算术平均值作为本次陷波值输出。滑动三角窗算术平均陷波器需要设置三角窗参数,所述三角窗参数具体指三角窗函数序列的长度NW。在三角窗参数NW值为信号周期序列长度的3倍,可以对1/3分次谐波产生的混频干扰频率进行抑制。在三角窗参数NW值为信号周期序列长度的4倍,可以对直流、1/2分次、1次、2次、3次、4次、5次谐波等产生的混频干扰频率进行抑制。
考虑到实际存在误差等因数,如所述参数存在1个采样间隔以内的误差,为了深度抑制混频干扰频率影响,还可再增加一级滑动矩形窗算术平均陷波器,即将若干个连续离散值直接相加,然后取其算术平均值作为本次陷波值输出。滑动矩形窗算术平均陷波器需要设置矩形窗参数,所述矩形窗参数具体指矩形窗函数序列的长度ND。矩形窗参数ND取值为信号周期序列长度的1.5倍,可以对1/3分次谐波产生的混频干扰频率进行抑制。而ND取值为信号周期序列长度的2倍,可以对直流、1/2分次、1次、2次、3次、4次、5次谐波等产生的混频干扰频率进行抑制。
优选地,数字陷波包括二级滑动三角窗算术平均陷波和二级滑动矩形窗算术平均陷波,该四级数字陷波式可为式(12):
ND1=(int)(1.5N2π+0.5)
ND2=(int)(2N2π+0.5)
对X(n)n=0,1,2,3,....,N-1
对W1(n)n=0,1,2,3,....,NW1-1
对W2(n)n=0,1,2,3,....,NW2-1
对XD(n)n=0,1,2,3,....,N-NW1-ND1-NW2-ND2-1
式中,X(n)为四级数字陷波输入序列,输入序列长度N;XD(n)为四级数字陷波输出序列,输出序列长度N-NW1-ND1-NW2-ND2;W1(n)为三角函数1,其中函数峰值为1,零频率增益为0.5;W2(n)为三角函数2,其中函数峰值为1,零频率增益为0.5;NW1为三角窗参数一,即三角窗函数1序列长度,单位无量纲,要求为奇数,以保证三角窗函数的外形为等腰三角形(如图5所示),(int)代表取整数;ND1为矩形窗参数一,即矩形窗函数1序列长度,单位无量纲,(int)代表取整数;NW2为三角窗参数二,即三角窗函数2序列长度,单位无量纲,要求为奇数,以保证三角窗函数的外形为等腰三角形(如图6所示),(int)代表取整数;ND2矩形窗参数二,即矩形窗函数2序列长度,单位无量纲,式中给出了计算公式,(int)代表取整数。
在一些实施例中,三角窗参数NW1取值为所述参考频率的单位周期序列长度的3倍,矩形窗参数ND1取值为所述参考频率的单位周期序列长度的1.5倍,三角窗参数NW2取值为所述参考频率的单位周期序列长度的4倍,矩形窗参数ND1取值为所述参考频率的单位周期序列长度的2倍。四级数字陷波需要使用10.5倍信号周期序列长度。
对上述实施例,在所述正弦信号的基波频率100π、单位rad/s,得到三角窗算术平均陷波算器1的频域特性,如图5所示。得到三角窗算术平均陷波算器2的频域特性,如图6所示。
优选地,在所述混频干扰频率成分得到完全抑制前提下,所述第一实频向量陷波序列和所述第一虚频向量陷波序列为(13):
对R1(n)I1(n) n=0,1,2,3,....,N-1
对RD1(n)ID1(n)n=0,1,2,3,....,N-NW1-ND1-NW2-ND2-1
式中,RD1(n)为所述第一实频向量陷波序列;ID1(n)为所述第一虚频向量陷波序列;K(Ω)为数字陷波在频差Ω的无量纲增益;α(Ω)为数字陷波在频差Ω的移相,单位rad。
对于步骤S111,优选地,可通过本领域技术人员惯用的积分器进行积分运算。
积分运算式为(14):
式中,R1为第一实频向量积分值;I1为第一虚频向量积分值。L1为第一积分长度,单位无量纲,L1为1.5倍所述单位周期序列长度。
对于步骤S112,同理和优选地,在所述混频干扰频率成分得到完全抑制前提下,所述第二实频向量陷波序列和所述第二虚频向量陷波序列为式(15):
对R2(n)I2(n)n=0,1,2,3,....,Ns-1
对RD2(n)ID2(n)n=Ns-NW1-ND1-NW2-ND2-1
其中,RD2(n)为所述第二实频向量陷波序列;ID2(n)为所述第二虚频向量陷波序列;K(Ω)为数字陷波在频差Ω的无量纲增益;α(Ω)为数字陷波在频差Ω的移相,单位rad。
对于步骤S113,优选地,积分运算式可为(16):
式中,R2为第二实频向量积分值;I2为第二虚频向量积分值。L2为第二积分长度,单位无量纲,L2为0.75倍所述单位周期序列长度。
对于步骤S114,优选地,预设的相位转换规则对应于虚频向量积分值和实频向量转换为相位的转换式,可通过以下公式(17)将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位,可通过以下公式(18)将所述第一虚频向量积分值与所述第一实频向量积分值转换为序列幅值;
式中,PH1为第一相位,单位rad。
式中,Amp为序列幅值,单位v。
在一些实施例中,根据预设的相位和幅值转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位和序列幅值的步骤包括以下步骤:
获取所述第一虚频向量积分值与所述第一实频向量积分值的比值。
获取所述比值的反正切函数值的相反数,生成所述第一相位。
获取所述第一虚频向量积分值与所述第一实频向量积分值的平方和。
对获取的平方和开根,生成所述序列幅值。
对于步骤S115,优选地,可通过以下公式(19)将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位:
式中,PH2为第二相位,单位rad。
在一些实施例中,根据所述预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位的步骤包括以下步骤:
获取所述第二虚频向量积分值与所述第二实频向量积分值的比值;
获取所述比值的反正切函数值的相反数,生成所述第二相位。
对于步骤S116,所述第一相位在0~±0.5πrad范围,但实际序列相位可能会超出±0.5πrad范围,因此必须根据相位扩展规则对所述第一相位进行扩展,扩展后的相位范围在0~±πrad范围,第一扩展相位,为式(20):
式中,Ph1为第一扩展相位,范围在0~±πard;&代表与逻辑。
在一些实施例中,根据预设的相位扩展规则,对所述第一相位进行扩展,获得第一扩展相位的步骤包括以下步骤:
如果所述第一实频向量积分值大于等于零的同时,所述第一虚频向量积分值的相反数大于等于零,则第一扩展相位等于所述第一相位;
如果所述第一实频向量积分值大于等于零的同时,所述第一虚频向量积分值的相反数小于零,则第一扩展相位等于所述第二相位;
如果所述第一实频向量积分值小于零的同时,所述第一虚频向量积分值的相反数大于等于零,则第一扩展相位等于所述第二相位加πrad;
如果所述第一实频向量积分值小于零的同时,所述第一虚频向量积分值的相反数小于零,则第一扩展相位等于所述第二相位减πrad;
对于步骤S117,所述第二相位在0~±0.5πrad范围,但实际序列相位可能会超出±0.5πrad范围,因此必须根据相位扩展规则对所述第二相位进行扩展,扩展后的相位范围在0~±πrad范围,第二扩展相位,为式(21):
式中,PH2为所述第二相位,范围在0~±0.5πard;Ph2为第二扩展相位,范围在0~±πard;&代表与逻辑。
在一些实施例中,根据预设的相位扩展规则,对所述第二相位进行扩展,获得第二扩展相位的步骤包括以下步骤:
如果所述第二实频向量积分值大于等于零的同时,所述第二虚频向量积分值的相反数大于等于零,则第二扩展相位等于所述第二相位;
如果所述第二实频向量积分值大于等于零的同时,所述第二虚频向量积分值的相反数小于零,则第二扩展相位等于所述第二相位;
如果所述第二实频向量积分值小于零的同时,所述第二虚频向量积分值的相反数大于等于零,则第二扩展相位等于所述第二相位加πrad;
如果所述第二实频向量积分值小于零的同时,所述第二虚频向量积分值的相反数小于零,则第二扩展相位等于所述第二相位减πrad;
对于步骤S118,优选地,所述预设的频差转换规则与将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号频率与所述参考频率的频差的转换式相应。
根据式(20)式(21),所述正弦信号频率与所述参考频率的频差计算式可为式(22):
其中,△ω为所述正弦信号频率与所述参考频率的频差,单位rad/s。
在一些实施例中,根据预设的频差转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号频率与所述参考频率的频差的步骤包括以下步骤:
获取对所述第一实频向量陷波序列进行积分运算的积分计算长度,生成第一积分长度。
获取对所述第二实频向量陷波序列进行积分运算的积分计算长度,生成第二积分长度。
获取所述第一扩展相位与所述第二扩展相位的差值,生成第一差值。
获取所述第一差值与所述预设采样频率的乘积的2倍,生成第一乘积。
获取所述第一积分长度与所述第二积分长度的差值,生成第二差值。
获取所述第一乘积与所述第二差值的比值,生成所述频差。
对于步骤S119,优选地,所述预设的幅值转换规则与将所述序列幅值和所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的幅值的转换式相应。所述正弦信号的幅值的计算式可为式(23):
其中
式中,AMP为所述正弦信号的基波幅值检测值,单位v;Amp为所述序列幅值,单位v;△ω为所述频差,单位rad/s;KL1(△ω)为所述第一积分在所述频差的幅值增益,单位无量纲;KW1(△ω)为所述滑动三角窗一算术平均陷波器在所述频差的幅值增益,单位无量纲;KD1(△ω)为所述滑动矩形窗一算术平均陷波器在在所述频差的幅值增益,单位无量纲;KW2(△ω)为所述滑动三角窗二算术平均陷波器在所述频差的幅值增益,单位无量纲;KD2(△ω)为所述滑动矩形窗二算术平均陷波器在所述频差的幅值增益,单位无量纲。
在一个实施例中,根据预设的幅值转换规则,将所述序列幅值、所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的基波幅值的步骤包括以下步骤:
获取所述序列幅值;
获取所述正弦信号频率与所述参考频率的频差;
获取所述采样间隔;
获取所述第一积分长度;
获取所述滑动三角窗一算术平均陷波器的三角窗参数一;
获取所述滑动矩形窗一算术平均陷波器的矩形窗参数一;
获取所述滑动三角窗二算术平均陷波器的三角窗参数二;
获取所述滑动矩形窗二算术平均陷波器的矩形窗参数二;
获取所述频差与所述采样间隔及所述第一积分长度的乘积的0.5倍,生成第一乘积;
获取所述频差与所述采样间隔及所述三角窗参数一的乘积的0.25倍,生成第二乘积;
获取所述频差与所述采样间隔及所述矩形窗参数一的乘积的0.5倍,生成第三乘积;
获取所述频差与所述采样间隔及所述三角窗参数二的乘积的0.25倍,生成第四乘积;
获取所述频差与所述采样间隔及所述矩形窗参数二的乘积的0.5倍,生成第五乘积;
获取所述第一乘积的正弦函数值与所述第一乘积的比值,生成所述第一积分运算在所述频差的幅值增益;
获取所述第二乘积的正弦函数值与所述第二乘积的比值的平方,生成所述滑动三角窗一算术平均陷波器在所述频差的幅值增益;
获取所述第三乘积的正弦函数值与所述第三乘积的比值,生成所述滑动矩形窗一算术平均陷波器在所述频差的幅值增益;
获取所述第四乘积的正弦函数值与所述第四乘积的比值的平方,生成所述滑动三角窗二算术平均陷波器在所述频差的幅值增益;
获取所述第五乘积的正弦函数值与所述第五乘积的比值,生成所述滑动矩形窗二算术平均陷波器在所述频差的幅值增益;
获取所述第一积分在所述频差的幅值增益、所述滑动三角窗一算术平均陷波器在所述频差的幅值增益、所述滑动矩形窗一算术平均陷波器在所述频差的幅值增益、所述滑动三角窗二算术平均陷波器在所述频差的幅值增益、所述滑动矩形窗二算术平均陷波器在所述频差的幅值增益的乘积,生成在所述频差的总幅值增益;
本发明还公开了正弦信号的幅值检测系统,在一些实施例中,如图2所示,该正弦信号的幅值检测系统,包括:初步序列长度模块1010、初步采样模块1020、频率初测模块1030、周期序列模块1040、序列长度模块1050、信号序列模块1060、截短序列模块1070、第一混频模块1080、第二混频模块1090、第一陷波模块1100、第一积分模块1110、第二陷波模块1120、第二积分模块1130、第一相位模块1140、第二相位模块1150、第一相位扩展模块1160、第二相位扩展模块1170、频差转换模块1180、幅值转换位模块1190。其中:
初步序列长度模块1010,用于根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度;
初步采样模块1020,用于根据所述初步采样序列长度,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列;
频率初测模块1030,用于对所述初步采样序列进行频率初测,生成所述正弦信号的初步频率,以所述初步频率给定参考频率;
周期序列模块1040,用于根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度;
序列长度模块1050,根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度;
信号序列模块1060,用于根据所述预设序列长度,从所述初步采样序列中,获得信号序列;
截短序列模块1070,用于对所述信号序列进行截短处理,获得截短信号序列;
第一混频模块1080,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,生成第一实频向量序列和第一虚频向量序列;
第二混频模块1090,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,生成第二实频向量序列和第二虚频向量序列;
第一陷波模块1100,用于分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波,生成第一实频向量陷波序列和第一虚频向量陷波序列;
第一积分模块1110,用于分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算,生成第一实频向量积分值和第一虚频向量积分值;
第二陷波模块1120,用于分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波,生成第二实频向量陷波序列和第二虚频向量陷波序列;
第二积分模块1130,用于分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算,生成第二实频向量积分值和第二虚频向量积分值;
第一相位模块1140,用于根据预设的相位和幅值转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位和序列幅值;
第二相位模块1150,用于根据所述预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位;
第一相位扩展模块1160,用于根据预设的相位扩展规则,将所述第一相位扩展为第一扩展相位;
第二相位扩展模块1170,用于根据预设的相位扩展规则,将所述第二相位扩展为第二扩展相位;
频差转换模块1180,用于根据预设的频差转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号频率与所述参考频率的频差;
幅值转换模块1190,用于根据预设的幅值转换规则,将所述序列幅值、所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的基波幅值。
本实施方式,对采样所得的信号序列进行截短处理,获得截短信号序列;以所测参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列和所述截短信号序列相乘,生成两组实频向量序列和虚频向量序列;通过对两组虚频向量序列和实频向量序列数字陷波,生成两组虚数向量陷波序列和实数向量陷波序列,进而积分生成两组虚数向量积分值和实数向量积分值;再根据预设的相位转换规则,将两组实数向量积分值和虚数向量积分值转换为序列幅值和两个相位;再将两个相位进行扩展,获得扩展相位;根据预设的频差转换规则,将两个扩展相位转换为所述正弦信号频率与所述参考频率的频差;根据预设的幅值转换规则,将所述序列幅值、所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的基波幅值。
由于数字陷波可以在某一个频率点迅速衰减输入信号,以达到阻碍此频率信号通过的效果,因此,当该数字陷波的陷波频率点设为对应混频干扰频率点时,该数字陷波对所述混频干扰频率具有深度的抑制作用。如此,采用本发明的正弦信号的幅值检测方法和设备,其对低频正弦信号的幅值测量精度高,且抗谐波和噪声干扰性好。
其中,初步序列长度模块1010,电力系统频率范围在45Hz-55Hz,优选的,取正弦信号频率下限fmin为45Hz;并根据实际需要设置所述预设信号周期数C2π,优选的,所述预设信号周期数C2π取整数。
在一些实施例中,取C2π为整数12。
初步采样序列长度计算,为上述式(1)。
对于初步采样模块1020,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列。对单基波频率正弦信号,获得所述正弦信号的初步采样序列为上述式(2)。
对于频率初测模块1030,可通过零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法对所述初步采样序列进行频率初测,获取所述初步频率。
所述初步频率表达为上述式(3)。
优选地,所述参考频率表达为上述式(4)。
对于周期序列模块1040,根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度:
在一些实施例中,所述正弦信号的单位周期序列长度计算,为上述式(5)。
所述单位周期序列长度整数化存在1个采样间隔内的误差。
对于序列长度模块1050,根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度:
在一些实施例中,所述预设序列长度为所述单位周期序列长度的12倍,所述预设序列长度计算,为上述式(6)。
对于信号序列模块1060,根据预设序列长度,从所述初步采样序列中,获得信号序列:
在一些实施例中,所述预设序列长度的信号序列,为上述式(7)。
所述信号序列的图形表达,如图3所示。
对于截短序列模块1070,对所述信号序列进行截短处理,获得截短信号序列:
在一些实施例中,所述截短信号序列长度为所述单位周期序列长度的11.25倍,所述截短信号序列长度计算为上述式(8)。
优选地,所述截短信号序列表达为上述式(9)。
所述截短信号序列的图形表达,如图3所示。
对于第一混频模块1080,优选地,所述参考频率的余弦函数和所述参考频率的正弦函数可分别为以所述参考频率为频率、以Tn为间隔离散变量的正弦函数和余弦函数。
在一些实施例中,以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,得到所述第一实频向量序列和所述第一虚频向量序列为上述式(10)。
对于第二混频模块1090,以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,得到所述第二实频向量序列和所述第二虚频向量序列为上述式(11)。
对于第一陷波模块1100,所述第一实频向量序列和所述第一虚频向量序列中包含混频干扰频率。当输入信号中还在直流成分、次谐波成分及分次谐波成分时,所述混频干扰频率将更加复杂,这些混频干扰频率严重影响计算准确度。虽然窗口函数和积分运算本身对混频干扰频率具有良好的衰减作用,但没有针对性,不能够对所述复杂的混频干扰频率产生深度的抑制作用,不能满足参数的高准确度计算需要。
为了有针对性的抑制所述混频干扰频率的影响,可以对该混频干扰进行数字陷波器,由于数字陷波可以在某一个频率点迅速衰减输入信号,以达到阻碍此频率信号通过的效果,因此,当该数字陷波的陷波频率点设为对应混频干扰频率点时,该数字陷波对所述混频干扰频率具有完全的抑制作用。
优选地,数字陷波器具体采用滑动三角窗算术平均陷波器,即将若干个连续离散值乘以三角窗函数后再相加,然后取其算术平均值作为本次陷波值输出。滑动三角窗算术平均陷波器需要设置三角窗参数,所述三角窗参数具体指三角窗函数序列的长度NW。在三角窗参数NW值为信号周期序列长度的3倍,可以对1/3分次谐波产生的混频干扰频率进行抑制。在三角窗参数NW值为信号周期序列长度的4倍,可以对直流、1/2分次、1次、2次、3次、4次、5次谐波等产生的混频干扰频率进行抑制。
考虑到实际存在误差等因数,如所述参数存在1个采样间隔以内的误差,为了深度抑制混频干扰频率影响,还可再增加一级滑动矩形窗算术平均陷波器,即将若干个连续离散值直接相加,然后取其算术平均值作为本次陷波值输出。滑动矩形窗算术平均陷波器需要设置矩形窗参数,所述矩形窗参数具体指矩形窗函数序列的长度ND。矩形窗参数ND取值为信号周期序列长度的1.5倍,可以对1/3分次谐波产生的混频干扰频率进行抑制。而ND取值为信号周期序列长度的2倍,可以对直流、1/2分次、1次、2次、3次、4次、5次谐波等产生的混频干扰频率进行抑制。
优选地,数字陷波器由二级滑动三角窗算术平均陷波器加二级滑动矩形窗算术平均陷波器共四级数字陷波器所构成该四级数字陷波器式可为上述式(12)。
在一些实施例中,三角窗参数NW1取值为所述参考频率的单位周期序列长度的3倍,矩形窗参数ND1取值为所述参考频率的单位周期序列长度的1.5倍,三角窗参数NW2取值为所述参考频率的单位周期序列长度的4倍,矩形窗参数ND1取值为所述参考频率的单位周期序列长度的2倍。四级数字陷波需要使用10.5倍信号周期序列长度。
对上述实施例,在所述正弦信号的基波频率100π、单位rad/s,得到三角窗算术平均陷波算器1的频域特性,如图5所示。得到三角窗算术平均陷波算器2的频域特性,如图6所示。
优选地,在所述混频干扰频率成分得到完全抑制前提下,所述第一实频向量陷波序列和所述第一虚频向量陷波序列为上述式(13)。
对于第一积分模块1110,优选地,可通过本领域技术人员惯用的积分器进行积分运算。积分运算式为上述式(14)。
对于第二陷波模块1120,同理和优选地,在所述混频干扰频率成分得到完全抑制前提下,所述第二实频向量陷波序列和所述第二虚频向量陷波序列为上述式(15)。
对于第二积分模块1130,优选地,积分运算式可为上述式(16)。
对于第一相位模块1140,优选地,预设的相位转换规则对应于虚频向量积分值和实频向量转换为相位的转换式,可通过上述公式(17)将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位,可通过上述公式(18)将所述第一虚频向量积分值与所述第一实频向量积分值转换为序列幅值。
在一些实施例中,第一相位模块1140可用于:
获取所述第一虚频向量积分值与所述第一实频向量积分值的比值。
获取所述比值的反正切函数值的相反数,生成所述第一相位。
获取所述第一虚频向量积分值与所述第一实频向量积分值的平方和。
对获取的平方和开根,生成所述序列幅值。
对于第二相位模块1150,优选地,可通过上述公式(19)将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位。
在一些实施例中,第二相位模块1150可用于:
获取所述第二虚频向量积分值与所述第二实频向量积分值的比值;
获取所述比值的反正切函数值的相反数,生成所述第二相位。
对于第一相位扩展模块1160,所述第一相位在0~±0.5πrad范围,但实际序列相位可能会超出±0.5πrad范围,因此必须根据相位扩展规则对所述第一相位进行扩展,扩展后的相位范围在0~±πrad范围,第一扩展相位,为上述式(20):
在一些实施例中,第一相位扩展模块1160可用于:
如果所述第一实频向量积分值大于等于零的同时,所述第一虚频向量积分值的相反数大于等于零,则第一扩展相位等于所述第一相位;
如果所述第一实频向量积分值大于等于零的同时,所述第一虚频向量积分值的相反数小于零,则第一扩展相位等于所述第二相位;
如果所述第一实频向量积分值小于零的同时,所述第一虚频向量积分值的相反数大于等于零,则第一扩展相位等于所述第二相位加πrad;
如果所述第一实频向量积分值小于零的同时,所述第一虚频向量积分值的相反数小于零,则第一扩展相位等于所述第二相位减πrad;
对于第二相位扩展模块1170,所述第二相位在0~±0.5πrad范围,但实际序列相位可能会超出±0.5πrad范围,因此必须根据相位扩展规则对所述第二相位进行扩展,扩展后的相位范围在0~±πrad范围,第二扩展相位,为上述式(20)。
在一些实施例中,第二相位扩展模块1170可用于:
如果所述第二实频向量积分值大于等于零的同时,所述第二虚频向量积分值的相反数大于等于零,则第二扩展相位等于所述第二相位;
如果所述第二实频向量积分值大于等于零的同时,所述第二虚频向量积分值的相反数小于零,则第二扩展相位等于所述第二相位;
如果所述第二实频向量积分值小于零的同时,所述第二虚频向量积分值的相反数大于等于零,则第二扩展相位等于所述第二相位加πrad;
如果所述第二实频向量积分值小于零的同时,所述第二虚频向量积分值的相反数小于零,则第二扩展相位等于所述第二相位减πrad;
对于频差转换模块1180,优选地,所述预设的频差转换规则与将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号频率与所述参考频率的频差的转换式相应。
根据式(20)式(21),所述正弦信号频率与所述参考频率的频差计算式可为上述式(22)。
在一些实施例中,频差转换模块1180可用于:
获取对所述第一实频向量陷波序列进行积分运算的积分计算长度,生成第一积分长度。
获取对所述第二实频向量陷波序列进行积分运算的积分计算长度,生成第二积分长度。
获取所述第一扩展相位与所述第二扩展相位的差值,生成第一差值。
获取所述第一差值与所述预设采样频率的乘积的2倍,生成第一乘积。
获取所述第一积分长度与所述第二积分长度的差值,生成第二差值。
获取所述第一乘积与所述第二差值的比值,生成所述频差。
对于幅值转换位模块1190,优选地,所述预设的幅值转换规则与将所述序列幅值和所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的幅值的转换式相应。所述正弦信号的幅值的计算式可为上述式(23)。
在一个实施例中,幅值转换位模块1190可用于:
获取所述序列幅值;
获取所述正弦信号频率与所述参考频率的频差;
获取所述采样间隔;
获取所述第一积分长度;
获取所述滑动三角窗一算术平均陷波器的三角窗参数一;
获取所述滑动矩形窗一算术平均陷波器的矩形窗参数一;
获取所述滑动三角窗二算术平均陷波器的三角窗参数二;
获取所述滑动矩形窗二算术平均陷波器的矩形窗参数二;
获取所述频差与所述采样间隔及所述第一积分长度的乘积的0.5倍,生成第一乘积;
获取所述频差与所述采样间隔及所述三角窗参数一的乘积的0.25倍,生成第二乘积;
获取所述频差与所述采样间隔及所述矩形窗参数一的乘积的0.5倍,生成第三乘积;
获取所述频差与所述采样间隔及所述三角窗参数二的乘积的0.25倍,生成第四乘积;
获取所述频差与所述采样间隔及所述矩形窗参数二的乘积的0.5倍,生成第五乘积;
获取所述第一乘积的正弦函数值与所述第一乘积的比值,生成所述第一积分运算在所述频差的幅值增益;
获取所述第二乘积的正弦函数值与所述第二乘积比值的平方,生成所述滑动三角窗一算术平均陷波器在所述频差的幅值增益;
获取所述第三乘积的正弦函数值与所述第三乘积的比值,生成所述滑动矩形窗一算术平均陷波器在所述频差的幅值增益;
获取所述第四乘积的正弦函数值与所述第四乘积比值的平方,生成所述滑动三角窗二算术平均陷波器在所述频差的幅值增益;
获取所述第五乘积的正弦函数值与所述第五乘积的比值,生成所述滑动矩形窗二算术平均陷波器在所述频差的幅值增益;
获取所述第一积分在所述频差的幅值增益、所述滑动三角窗一算术平均陷波器在所述频差的幅值增益、所述滑动矩形窗一算术平均陷波器在所述频差的幅值增益、所述滑动三角窗二算术平均陷波器在所述频差的幅值增益、所述滑动矩形窗二算术平均陷波器在所述频差的幅值增益的乘积,生成在所述频差的总幅值增益;
获取所述序列幅值与所述总幅值增益的比值,生成所述正弦信号的基波幅值检测值。
为了验证本发明正弦信号的幅值检测系统具有较高的准确度,给出一实验信号,为式(24):
ω=2πft
在实验信号基波频率变化范围在45Hz-55Hz,取信号整数周期数约为12,信号初相位变化范围在0~±π/2、单位rad,信号的采样频率为20kHz,信号的离散数据量化位数24bit,频率初测相对误差<|±0.25%|,得到信号基波幅值检测的相对误差绝对值|AMPerr(f)|随信号基波频率f变化特性的实验结果图,图4所示。图4给出的实验信号的幅值检测准确度在10-7量级。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种正弦信号的幅值检测方法,其特征在于,包括以下步骤:
根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度;
根据所述初步采样序列长度,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列;
对所述初步采样序列进行频率初测,生成所述正弦信号的初步频率,以所述初步频率给定参考频率;
根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度;
根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度;
根据预设序列长度,从所述初步采样序列中,获得信号序列;
对所述信号序列进行截短处理,获得截短信号序列;
以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,生成第一实频向量序列和第一虚频向量序列;
以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,生成第二实频向量序列和第二虚频向量序列;
分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波,生成第一实频向量陷波序列和第一虚频向量陷波序列;
分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算,生成第一实频向量积分值和第一虚频向量积分值;
分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波,生成第二实频向量陷波序列和第二虚频向量陷波序列;
分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算,生成第二实频向量积分值和第二虚频向量积分值;
根据预设的相位和幅值转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位和序列幅值;所述预设的相位和幅值转换规则为:通过公式将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位,其中PH1为第一相位,单位rad;通过公式将所述第一虚频向量积分值与所述第一实频向量积分值转换为序列幅值,其中Amp为序列幅值,单位v;R1为第一实频向量积分值;I1为第一虚频向量积分值;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;L1为第一积分长度,单位无量纲,L1为1.5倍所述单位周期序列长度;Ω为信号频率与参考频率的频差,单位rad/s;α(Ω)为数字陷波在频差Ω的移相,单位rad;为信号初相位,单位rad;N信号序列长度,单位无量纲,信号序列长度等于所述预设序列长度;A为信号幅值,单位v;K(Ω)为数字陷波在频差Ω的无量纲增益;
根据预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位;所述预设的相位转换规则为:通过公式将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位,其中PH2为第二相位,单位rad;R2为第二实频向量积分值;I2为第二虚频向量积分值;L2为第二积分长度,单位无量纲,L2为0.75倍所述单位周期序列长度;NS为截短信号序列长度,单位无量纲;
根据预设的第一相位扩展规则,对所述第一相位进行扩展,获得第一扩展相位;所述预设的第一相位扩展规则为:通过公式获得第一扩展相位,其中Ph1为第一扩展相位,范围在0~±πard;&代表与逻辑;
根据预设的第二相位扩展规则,对所述第二相位进行扩展,获得第二扩展相位;所述预设的第二相位扩展规则为:通过公式获得第二扩展相位,其中PH2为所述第二相位,范围在0~±0.5πard;Ph2为第二扩展相位,范围在0~±πard;
根据预设的频差转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号频率与所述参考频率的频差;所述预设的频差转换规则为:通过公式获得所述正弦信号频率与所述参考频率的频差,其中其中,△ω为所述正弦信号频率与所述参考频率的频差,单位rad/s;
根据预设的幅值转换规则,将所述序列幅值、所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的基波幅值;所述预设的幅值转换规则为:通过公式获得正弦信号的基波幅值,其中AMP为所述正弦信号的基波幅值,单位v;Amp为所述序列幅值,单位v;△ω为所述频差,单位rad/s;KL1(△ω)为所述第一积分在所述频差的幅值增益,单位无量纲;KW1(△ω)为滑动三角窗一算术平均陷波器在所述频差的幅值增益,单位无量纲;KD1(△ω)为滑动矩形窗一算术平均陷波器在所述频差的幅值增益,单位无量纲;KW2(△ω)为滑动三角窗二算术平均陷波器在所述频差的幅值增益,单位无量纲;KD2(△ω)为滑动矩形窗二算术平均陷波器在所述频差的幅值增益,单位无量纲;NW1为三角窗参数一,即三角窗函数1序列长度,单位无量纲;ND1为矩形窗参数一,即矩形窗函数1序列长度,单位无量纲;NW2为三角窗参数二,即三角窗函数2序列长度,单位无量纲;ND2矩形窗参数二,即矩形窗函数2序列长度,单位无量纲。
2.根据权利要求1所述的正弦信号的幅值检测方法,其特征在于,根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度的步骤包括以下步骤:
通过以下所述公式将所述正弦信号频率的下限和预设采样频率及预设信号周期数转换为所述初步采样序列长度:
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>int</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msub>
<mi>c</mi>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msub>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
</mrow>
<msub>
<mi>f</mi>
<mi>min</mi>
</msub>
</mfrac>
<mo>;</mo>
</mrow>
其中,Nstart为所述初步采样序列长度,单位无量纲;(int)表示取整;C2π为所述预设信号周期数,单位无量纲;fmin为正弦信号频率范围的下限,单位Hz;fn为所述预设采样频率,单位Hz。
3.根据权利要求1所述的正弦信号的幅值检测方法,其特征在于,所述数字陷波由二级滑动三角窗算术平均陷波和二级滑动矩形窗算术平均陷波构成。
4.根据权利要求1所述的正弦信号的幅值检测方法,其特征在于,所述初步频率通过对所述初步采样序列采用零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法获得。
5.根据权利要求1所述的正弦信号的幅值检测方法,其特征在于,所述相位扩展规则为:
如果所述实频向量积分值大于等于零的同时,所述虚频向量积分值的相反数大于等于零,则所述扩展相位等于所述相位;
如果所述实频向量积分值大于等于零的同时,所述虚频向量积分值的相反数小于零,则所述扩展相位等于所述相位;
如果所述实频向量积分值小于零的同时,所述虚频向量积分值的相反数大于等于零,则所述扩展相位等于所述相位加πrad;
如果所述实频向量积分值小于零的同时,所述虚频向量积分值的相反数小于零,则所述扩展相位等于所述相位减πrad。
6.一种正弦信号的幅值检测系统,其特征在于,包括:
初步序列长度模块,用于根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度;
初步采样模块,用于根据所述初步采样序列长度,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列;
频率初测模块,用于对所述初步采样序列进行频率初测,生成所述正弦信号的初步频率,以所述初步频率给定参考频率;
周期序列模块,用于根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度;
序列长度模块,根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度;
信号序列模块,用于根据所述预设序列长度,从所述初步采样序列中,获得信号序列;
截短序列模块,用于对所述信号序列进行截短处理,获得截短信号序列;
第一混频模块,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,生成第一实频向量序列和第一虚频向量序列;
第二混频模块,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,生成第二实频向量序列和第二虚频向量序列;
第一陷波模块,用于分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波,生成第一实频向量陷波序列和第一虚频向量陷波序列;
第一积分模块,用于分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算,生成第一实频向量积分值和第一虚频向量积分值;
第二陷波模块,用于分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波,生成第二实频向量陷波序列和第二虚频向量陷波序列;
第二积分模块,用于分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算,生成第二实频向量积分值和第二虚频向量积分值;
第一相位模块,用于根据预设的相位和幅值转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位和序列幅值;所述预设的相位和幅值转换规则为:通过公式将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位,其中PH1为第一相位,单位rad;通过公式将所述第一虚频向量积分值与所述第一实频向量积分值转换为序列幅值,其中Amp为序列幅值,单位v;R1为第一实频向量积分值;I1为第一虚频向量积分值;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;L1为第一积分长度,单位无量纲,L1为1.5倍所述单位周期序列长度;Ω为信号频率与参考频率的频差,单位rad/s;α(Ω)为数字陷波在频差Ω的移相,单位rad;为信号初相位,单位rad;N信号序列长度,单位无量纲,信号序列长度等于所述预设序列长度;A为信号幅值,单位v;K(Ω)为数字陷波在频差Ω的无量纲增益;
第二相位模块,用于根据预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位;所述预设的相位转换规则为:通过公式将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位,其中PH2为第二相位,单位rad;R2为第二实频向量积分值;I2为第二虚频向量积分值;L2为第二积分长度,单位无量纲,L2为0.75倍所述单位周期序列长度;NS为截短信号序列长度,单位无量纲;
第一相位扩展模块,用于根据预设的第一相位扩展规则,将所述第一相位扩展为第一扩展相位;所述预设的第一相位扩展规则为:通过公式获得第一扩展相位,其中Ph1为第一扩展相位,范围在0~±πard;&代表与逻辑;
第二相位扩展模块,用于根据预设的第二相位扩展规则,将所述第二相位扩展为第二扩展相位;所述预设的第二相位扩展规则为:通过公式获得第二扩展相位,其中PH2为所述第二相位,范围在0~±0.5πard;Ph2为第二扩展相位,范围在0~±πard;
频差转换模块,用于根据预设的频差转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号频率与所述参考频率的频差;所述预设的频差转换规则为:通过公式获得所述正弦信号频率与所述参考频率的频差,其中其中,△ω为所述正弦信号频率与所述参考频率的频差,单位rad/s;
幅值转换模块,用于根据预设的幅值转换规则,将所述序列幅值、所述正弦信号频率与所述参考频率的频差转换为所述正弦信号的基波幅值;所述预设的幅值转换规则为:通过公式获得正弦信号的基波幅值,其中AMP为所述正弦信号的基波幅值,单位v;Amp为所述序列幅值,单位v;△ω为所述频差,单位rad/s;KL1(△ω)为所述第一积分在所述频差的幅值增益,单位无量纲;KW1(△ω)为滑动三角窗一算术平均陷波器在所述频差的幅值增益,单位无量纲;KD1(△ω)为滑动矩形窗一算术平均陷波器在所述频差的幅值增益,单位无量纲;KW2(△ω)为滑动三角窗二算术平均陷波器在所述频差的幅值增益,单位无量纲;KD2(△ω)为滑动矩形窗二算术平均陷波器在所述频差的幅值增益,单位无量纲;NW1为三角窗参数一,即三角窗函数1序列长度,单位无量纲;ND1为矩形窗参数一,即矩形窗函数1序列长度,单位无量纲;NW2为三角窗参数二,即三角窗函数2序列长度,单位无量纲;ND2矩形窗参数二,即矩形窗函数2序列长度,单位无量纲。
7.根据权利要求6所述的正弦信号的幅值检测系统,其特征在于,所述初步采样模块根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度的步骤包括以下步骤:
通过以下所述公式将所述正弦信号频率的下限和预设采样频率及预设信号周期数转换为所述初步采样序列长度:
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>int</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msub>
<mi>c</mi>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msub>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
</mrow>
<msub>
<mi>f</mi>
<mi>min</mi>
</msub>
</mfrac>
<mo>;</mo>
</mrow>
其中,Nstart为所述初步采样序列长度,单位无量纲;(int)表示取整;C2π为所述预设信号周期数,单位无量纲;fmin为正弦信号频率范围的下限,单位Hz;fn为所述预设采样频率,单位Hz。
8.根据权利要求6所述的正弦信号的幅值检测系统,其特征在于,所述第一陷波模块、第二陷波模块包括二滑动三角窗算术平均陷波器和二级滑动矩形窗算术平均陷波器。
9.根据权利要求6所述的正弦信号的幅值检测系统,其特征在于,所述初步频率通过对所述初步采样序列采用零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法获得。
10.根据权利要求6所述的正弦信号的幅值检测系统,其特征在于,所述第一相位扩展模块和第二相位扩展模块根据预设的相位扩展规则,将所述第一相位或所述第二相位扩展为所第一扩展相位或所述第二扩展相位的步骤包括以下步骤:
如果所述实频向量积分值大于等于零的同时,所述虚频向量积分值的相反数大于等于零,则所述扩展相位等于所述相位;
如果所述实频向量积分值大于等于零的同时,所述虚频向量积分值的相反数小于零,则所述扩展相位等于所述相位;
如果所述实频向量积分值小于零的同时,所述虚频向量积分值的相反数大于等于零,则所述扩展相位等于所述相位加πrad;
如果所述实频向量积分值小于零的同时,所述虚频向量积分值的相反数小于零,则所述扩展相位等于所述相位减πrad。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510400897.1A CN105372471B (zh) | 2015-07-09 | 2015-07-09 | 正弦信号的幅值检测方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510400897.1A CN105372471B (zh) | 2015-07-09 | 2015-07-09 | 正弦信号的幅值检测方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105372471A CN105372471A (zh) | 2016-03-02 |
CN105372471B true CN105372471B (zh) | 2018-05-25 |
Family
ID=55374846
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510400897.1A Active CN105372471B (zh) | 2015-07-09 | 2015-07-09 | 正弦信号的幅值检测方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105372471B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108020715A (zh) * | 2017-12-07 | 2018-05-11 | 中国船舶重工集团公司第七〇九研究所 | 一种单相正弦信号幅值与相位快速检测方法 |
CN111134637B (zh) * | 2019-12-06 | 2023-04-07 | 深圳绿米联创科技有限公司 | 信号处理方法、装置、电子设备及存储介质 |
CN113180689A (zh) * | 2021-04-08 | 2021-07-30 | 上海数创医疗科技有限公司 | 一种基于高斯函数拟合的p波波峰及起始位置提取方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0106029A1 (en) * | 1982-07-19 | 1984-04-25 | CSELT Centro Studi e Laboratori Telecomunicazioni S.p.A. | Method and apparatus for measuring the amplitude of a noise-affected periodic signal without phase reference |
CN102798748A (zh) * | 2012-07-13 | 2012-11-28 | 中冶南方工程技术有限公司 | 一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法 |
CN104198825A (zh) * | 2014-08-05 | 2014-12-10 | 广东电网公司电力科学研究院 | 检测电力系统中信号幅值的方法及系统 |
CN104330644A (zh) * | 2014-10-29 | 2015-02-04 | 广东电网有限责任公司电力科学研究院 | 检测电力系统中正弦波信号幅值的方法和系统 |
CN104502706A (zh) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | 电力信号的谐波幅值测量方法和系统 |
CN104502675A (zh) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | 电力信号的基波幅值测量方法和系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0954124A (ja) * | 1995-08-17 | 1997-02-25 | Hitachi Ltd | 周波数判定方法及びディジタル信号処理装置 |
-
2015
- 2015-07-09 CN CN201510400897.1A patent/CN105372471B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0106029A1 (en) * | 1982-07-19 | 1984-04-25 | CSELT Centro Studi e Laboratori Telecomunicazioni S.p.A. | Method and apparatus for measuring the amplitude of a noise-affected periodic signal without phase reference |
CN102798748A (zh) * | 2012-07-13 | 2012-11-28 | 中冶南方工程技术有限公司 | 一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法 |
CN104198825A (zh) * | 2014-08-05 | 2014-12-10 | 广东电网公司电力科学研究院 | 检测电力系统中信号幅值的方法及系统 |
CN104330644A (zh) * | 2014-10-29 | 2015-02-04 | 广东电网有限责任公司电力科学研究院 | 检测电力系统中正弦波信号幅值的方法和系统 |
CN104502706A (zh) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | 电力信号的谐波幅值测量方法和系统 |
CN104502675A (zh) * | 2014-12-29 | 2015-04-08 | 广东电网有限责任公司电力科学研究院 | 电力信号的基波幅值测量方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105372471A (zh) | 2016-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zygarlicki et al. | A reduced Prony's method in power-quality analysis—parameters selection | |
Zhu | Exact harmonics/interharmonics calculation using adaptive window width | |
CN104635044A (zh) | 基于幅值调制的电力信号频率检测方法和系统 | |
Singh | Time-frequency analysis via the fourier representation | |
CN107102255B (zh) | 单一adc采集通道动态特性测试方法 | |
Tiwari et al. | Fast amplitude estimation of harmonics using undecimated wavelet packet transform and its hardware implementation | |
CN105372471B (zh) | 正弦信号的幅值检测方法和系统 | |
Djurić et al. | Frequency measurement of distorted signals using Fourier and zero crossing techniques | |
Belega et al. | Accuracy analysis of the sine-wave parameters estimation by means of the windowed three-parameter sine-fit algorithm | |
Csurcsia | Static nonlinearity handling using best linear approximation: An introduction | |
Li et al. | Frequency estimation based on modulation FFT and MUSIC algorithm | |
Petrović | Frequency and parameter estimation of multi-sinusoidal signal | |
CN105445547B (zh) | 正弦信号的频率检测方法和系统 | |
Ojarand et al. | Recent advances in crest factor minimization of multisine | |
Tao et al. | A robust parametric method for power harmonic estimation based on M-estimators | |
CN109557367B (zh) | 一种高频分辨率谐波和间谐波Prony方法及装置 | |
Liu et al. | An approximate maximum likelihood estimator for instantaneous frequency estimation of multicomponent nonstationary signals | |
CN104849551B (zh) | 一种谐相角分析方法 | |
Luo et al. | Improved zero-crossing method for power system frequency estimation | |
Xie | Frequency Estimation of Weighted Signals Based On DFT Interpolation Algorithm | |
Zolfaghari et al. | Spectral analysis techniques with Kalman filtering for estimating power quality indices | |
Meng et al. | An Instantaneous Frequency Identification Algorithm for the Dynamic Frequency Signals in power system | |
Abiyev | A new electronic reactive power meter based on Walsh Transform algorithms | |
CN109709397B (zh) | 一种加连续Hanning窗的电网谐波非同步压缩感知检测方法 | |
CN104977467B (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 |