CN105445552B - 正弦信号的初相位检测方法和系统 - Google Patents
正弦信号的初相位检测方法和系统 Download PDFInfo
- Publication number
- CN105445552B CN105445552B CN201510400898.6A CN201510400898A CN105445552B CN 105445552 B CN105445552 B CN 105445552B CN 201510400898 A CN201510400898 A CN 201510400898A CN 105445552 B CN105445552 B CN 105445552B
- Authority
- CN
- China
- Prior art keywords
- frequency
- phase
- sequence
- integrated value
- 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
Abstract
本发明涉及一种正弦信号的初相位检测方法和系统,所述方法包括:对采样所得的信号序列进行截短处理,获得截短信号序列;以所测参考频率的余弦函数和参考频率的正弦函数分别与信号序列和截短信号序列相乘,生成两组实频向量序列和虚频向量序列;通过对两组虚频向量序列和实频向量序列数字陷波,生成两组虚数向量陷波序列和实数向量陷波序列,进而积分生成两组虚数向量积分值和实数向量积分值;再根据预设的相位转换规则,将两组实数向量积分值和虚数向量积分值转换为两个相位,再将两个相位进行扩展,获得扩展相位;进而根据预设的初相位转换规则,将两个扩展相位转换为所述正弦信号的初相位。实施本发明,获得准确度较高的初相位。
Description
技术领域
本发明涉及电力技术领域,特别是涉及一种正弦信号的初相位检测方法和系统。
背景技术
电力系统的频率测量、相位测量、幅值测量等在本质上均为正弦信号参数的测量。傅里叶变换是实现正弦信号参数测量的基本方法,在电力系统中有广泛的应用。但随着正弦参数测量技术的发展,傅里叶变换存在的问题也越显突出,其难以进一步满足电力系统对正弦参数高准确度计算的要求。
在电力系统正弦信号参数测量方面,还有一些改进的参数测量方法,如零交法、基于滤波的测量法、基于小波变换法、基于神经网络的测量法、基于DFT变换的测量法等。但是,以上所述方法对低频信号(比如,电网运行额定工频在50Hz附近)的测量精度低,且抗谐波和噪声干扰性差。
发明内容
本发明所要解决的技术问题在于,提供一种正弦信号的初相位检测方法,其对电力系统的低频正弦信号的正弦参数测量精度高,且抗谐波和噪声干扰性好。
为解决上述技术问题,本发明采用如下技术方案:
一种正弦信号的初相位检测方法,包括以下步骤:
根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度;
根据所述初步采样序列长度,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列;
对所述初步采样序列进行频率初测,生成所述正弦信号的初步频率,以所述初步频率给定参考频率;
根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度;
根据预设信号周期数和预设采样频率,获得预设序列长度;
根据所述预设序列长度,从所述初步采样序列中,获得所述预设序列长度的信号序列;
对所述信号序列进行截短处理,获得截短信号序列,所述截短信号序列的长度相对所述信号序列的长度的截短值为所述单位周期序列长度的0.75倍;
以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,生成第一实频向量序列和第一虚频向量序列;
以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,生成第二实频向量序列和第二虚频向量序列;
分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波,生成第一实频向量陷波序列和第一虚频向量陷波序列;
分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算,生成第一实频向量积分值和第一虚频向量积分值;
分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波,生成第二实频向量陷波序列和第二虚频向量陷波序列;
分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算,生成第二实频向量积分值和第二虚频向量积分值;
根据预设的相位转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位;
根据预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位;
根据预设的相位扩展规则,对所述第一相位进行扩展,获得第一扩展相位;
根据预设的相位扩展规则,对所述第二相位进行扩展,获得第二扩展相位;
根据预设的初相位转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的初相位。
本发明所要解决的技术问题还在于提供一种正弦信号的初相位检测系统,其对低频正弦信号的正弦参数测量精度高,且抗谐波和噪声干扰性好。
为解决上述技术问题,本发明采用如下技术方案:
一种正弦信号的初相位检测系统,包括:
初步序列长度模块,用于根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度;
初步采样模块,用于根据正弦信号频率范围的下限和预设采样频率及预设信号周期数对正弦信号进行采样,获得所述正弦信号的初步采样序列;
频率初测模块,用于对所述初步采样序列进行频率初测,生成所述正弦信号的初步频率,以所述初步频率给定参考频率;
周期序列模块,用于根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度;
序列长度模块,根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度;
信号序列模块,用于根据所述预设序列长度,从所述初步采样序列中,获得所述预设序列长度的信号序列;
截短序列模块,用于对所述信号序列进行截短处理,获得截短信号序列,;
第一混频模块,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,生成第一实频向量序列和第一虚频向量序列;
第二混频模块,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,生成第二实频向量序列和第二虚频向量序列;
第一陷波模块,用于分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波,生成第一实频向量陷波序列和第一虚频向量陷波序列;
第一积分模块,用于分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算,生成第一实频向量积分值和第一虚频向量积分值;
第二陷波模块,用于分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波,生成第二实频向量陷波序列和第二虚频向量陷波序列;
第二积分模块,用于分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算,生成第二实频向量积分值和第二虚频向量积分值;
第一相位模块,用于根据预设的相位转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位;
第二相位模块,用于根据所述预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位;
第一相位扩展模块,用于根据预设的相位扩展规则,将所述第一相位扩展为第一扩展相位;
第二相位扩展模块,用于根据预设的相位扩展规则,将所述第二相位扩展为第二扩展相位;
初相位模块,用于根据预设的初相位转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的初相位。
本发明与现有技术相比的有益效果是:对采样所得的信号序列进行截短处理,获得截短信号序列;以所测参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列和所述截短信号序列相乘,生成两组实频向量序列和虚频向量序列;通过对两组虚频向量序列和实频向量序列数字陷波,生成两组虚数向量陷波序列和实数向量陷波序列,进而积分生成两组虚数向量积分值和实数向量积分值;再根据预设的相位转换规则,将两组实数向量积分值和虚数向量积分值转换为两个相位;再将两个相位进行扩展,获得扩展相位。进而根据预设的初相位转换规则,将两个扩展相位转换为所述正弦信号的初相位。
由于数字陷波可以在某一个频率点迅速衰减输入信号,以达到阻碍此频率信号通过的效果,因此,当该数字陷波的陷波频率点设为对应混频干扰频率点时,该数字陷波对所述混频干扰频率具有深度的抑制作用。如此,采用本发明的正弦信号的初相位检测方法和设备,其对低频正弦信号的正弦参数测量精度高,且抗谐波和噪声干扰性好。
附图说明
图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根据预设的初相位转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的初相位。
本实施方式,对采样所得的信号序列进行截短处理,获得截短信号序列;以所测参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列和所述截短信号序列相乘,生成两组实频向量序列和虚频向量序列;通过对两组虚频向量序列和实频向量序列数字陷波,生成两组虚数向量陷波序列和实数向量陷波序列,进而积分生成两组虚数向量积分值和实数向量积分值;再根据预设的相位转换规则,将两组实数向量积分值和虚数向量积分值转换为两个相位;再将两个相位进行扩展,获得扩展相位;进而根据预设的初相位转换规则,将两个扩展相位转换为所述正弦信号的初相位,具有较高的准确度。
之后如不加说明,所述正弦信号的初相位均指正弦信号的基波初相位。
其中,对于步骤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):
式中,X(n)为四级数字陷波输入序列,输入序列长度N;XD(n)为四级数字陷波输出序列,输出序列长度N-NW1-ND1-NW2-ND2;W1(n)为三角函数1,其中函数峰值为1,零频率增益为0.5;W2(n)为三角函数2,其中函数峰值为1,零频率增益为0.5;NW1为三角窗参数1,即三角窗函数1序列长度,单位无量纲,要求为奇数,以保证三角窗函数的外形为等腰三角形(如图5所示),(int)代表取整数;ND1为矩形窗参数1,即矩形窗函数1序列长度,单位无量纲,(int)代表取整数;NW2为三角窗参数2,即三角窗函数2序列长度,单位无量纲,要求为奇数,以保证三角窗函数的外形为等腰三角形(如图6所示),(int)代表取整数;ND2矩形窗参数2,即矩形窗函数2序列长度,单位无量纲,式中给出了计算公式,(int)代表取整数。
在一些实施例中,三角窗参数NW1取值为所述参考频率的单位周期序列长度的3倍,矩形窗参数ND1取值为所述参考频率的单位周期序列长度的1.5倍,三角窗参数NW2取值为所述参考频率的单位周期序列长度的4倍,矩形窗参数ND1取值为所述参考频率的单位周期序列长度的2倍。四级数字陷波需要使用10.5倍所述单位周期序列长度。
对上述实施例,在所述正弦信号的基波频率100π、单位rad/s,得到三角窗算术平均陷波算器1的频域特性,如图5所示。得到三角窗算术平均陷波算器2的频域特性,如图6所示。
优选地,在所述混频干扰频率成分得到完全抑制前提下,所述第一实频向量陷波序列和所述第一虚频向量陷波序列为(13):
式中,RD1(n)为所述第一实频向量陷波序列;ID1(n)为所述第一虚频向量陷波序列;K(Ω)为数字陷波在频差Ω的无量纲增益;α(Ω)为数字陷波在频差Ω的移相,单位rad。
对于步骤S111,优选地,可通过本领域技术人员惯用的积分器进行积分运算。
积分运算式为(14):
式中,R1为第一实频向量积分值;I1为第一虚频向量积分值。L1为第一积分长度,单位无量纲,L1为1.5倍所述单位周期序列长度。
对于步骤S112,同理和优选地,在所述混频干扰频率成分得到完全抑制前提下,所述第二实频向量陷波序列和所述第二虚频向量陷波序列为式(15):
其中,RD2(n)为所述第二实频向量陷波序列;ID2(n)为所述第二虚频向量陷波序列;K(Ω)为数字陷波在频差Ω的无量纲增益;α(Ω)为数字陷波在频差Ω的移相,单位rad。
对于步骤S113,优选地,积分运算式可为(16):
式中,R2为第二实频向量积分值;I2为第二虚频向量积分值。L2为第二积分长度,单位无量纲,L2为0.75倍所述单位周期序列长度。
对于步骤S114,优选地,预设的相位转换规则对应于虚频向量积分值和实频向量转换为相位的转换式,可通过以下公式(17)将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位:
式中,PH1为第一相位,单位rad。
在一些实施例中,根据预设的相位转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位的步骤包括以下步骤:
获取所述第一虚频向量积分值与所述第一实频向量积分值的比值;
获取所述比值的反正切函数值的相反数,生成所述第一相位。
对于步骤S115,优选地,可通过以下公式(18)将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位:
式中,PH2为第二相位,单位rad。
在一些实施例中,根据所述预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位的步骤包括以下步骤:
获取所述第二虚频向量积分值与所述第二实频向量积分值的比值;
获取所述比值的反正切函数值的相反数,生成所述第二相位。
对于步骤S116,所述第一相位在0~±0.5πrad范围,但实际序列相位可能会超出±0.5πrad范围,因此必须根据相位扩展规则对所述第一相位进行扩展,扩展后的相位范围在0~±πrad范围,第一扩展相位,为式(19):
式中,Ph1为第一扩展相位,范围在0~±πrad;&代表与逻辑。
在一些实施例中,根据预设的相位扩展规则,对所述第一相位进行扩展,获得第一扩展相位的步骤包括以下步骤:
如果所述第一实频向量积分值大于等于零的同时,所述第一虚频向量积分值的相反数大于等于零,则第一扩展相位等于所述第一相位;
如果所述第一实频向量积分值大于等于零的同时,所述第一虚频向量积分值的相反数小于零,则第一扩展相位等于所述第二相位;
如果所述第一实频向量积分值小于零的同时,所述第一虚频向量积分值的相反数大于等于零,则第一扩展相位等于所述第二相位加πrad;
如果所述第一实频向量积分值小于零的同时,所述第一虚频向量积分值的相反数小于零,则第一扩展相位等于所述第二相位减πrad;
对于步骤S117,所述第二相位在0~±0.5πrad范围,但实际序列相位可能会超出±0.5πrad范围,因此必须根据相位扩展规则对所述第二相位进行扩展,扩展后的相位范围在0~±πrad范围,第二扩展相位,为式(20):
式中,PH2为所述第二相位,范围在0~±0.5πrad;Ph2为第二扩展相位,范围在0~±πrad;&代表与逻辑。
在一些实施例中,根据预设的相位扩展规则,对所述第二相位进行扩展,获得第二扩展相位的步骤包括以下步骤:
如果所述第二实频向量积分值大于等于零的同时,所述第二虚频向量积分值的相反数大于等于零,则第二扩展相位等于所述第二相位;
如果所述第二实频向量积分值大于等于零的同时,所述第二虚频向量积分值的相反数小于零,则第二扩展相位等于所述第二相位;
如果所述第二实频向量积分值小于零的同时,所述第二虚频向量积分值的相反数大于等于零,则第二扩展相位等于所述第二相位加πrad;
如果所述第二实频向量积分值小于零的同时,所述第二虚频向量积分值的相反数小于零,则第二扩展相位等于所述第二相位减πrad;
对于步骤S118,所述预设的初相位转换规则可对应于第一扩展相位和第二扩展相位转换为初相位的公式。根据式(19)和式(20),可生成与所述预设的初相位转换规则对应的初相位公式(21):
式中,为正弦信号初相位检测值,单位rad。
在一些实施例中,根据预设的初相位转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的初相位的步骤包括以下步骤:
获取所述第一扩展相位与所述截短信号序列的序列长度的乘积,生成第一乘积。
获取所述第二扩展相位与所述预设序列长度的乘积,生成第二乘积。
获取所述第一乘积与所述第二乘积的差值,生成第一差值。
获取所述截短信号序列的序列长度与所述预设序列长度的差值,生成第二差值。
获取所述第一差值与所述第二差值的比值,生成所述初相位。
本发明还公开了正弦信号的初相位检测系统,在一些实施例中,如图2所示,该正弦信号的初相位检测系统,包括:初步序列长度模块1010、初步采样模块1020、频率初测模块1030、周期序列模块1040、序列长度模块1050、信号序列模块1060、截短序列模块1070、第一混频模块1080、第二混频模块1090、第一陷波模块1100、第一积分模块1110、第二陷波模块1120、第二积分模块1130、第一相位模块1140、第二相位模块1150、第一相位扩展模块1160、第二相位扩展模块1170、初相位模块1180。其中:
初步序列长度模块1010,用于根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度;
初步采样模块1020,用于根据正弦信号频率范围的下限和预设采样频率及预设信号周期数对正弦信号进行采样,获得所述正弦信号的初步采样序列。
频率初测模块1030,用于对所述初步采样序列进行频率初测,生成所述正弦信号的初步频率,以所述初步频率给定参考频率。
周期序列模块1040,用于根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度。
序列长度模块1050,根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度。
信号序列模块1060,根据所述预设序列长度,从所述初步采样序列中,获得所述预设序列长度的信号序列。
截短序列模块1070,用于对所述信号序列进行截短处理,获得截短信号序列。
第一混频模块1080,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,生成第一实频向量序列和第一虚频向量序列。
第二混频模块1090,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,生成第二实频向量序列和第二虚频向量序列。
第一陷波模块1100,用于分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波,生成第一实频向量陷波序列和第一虚频向量陷波序列。
第一积分模块1110,用于分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算,生成第一实频向量积分值和第一虚频向量积分值。
第二陷波模块1120,用于分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波,生成第二实频向量陷波序列和第二虚频向量陷波序列。
第二积分模块1130,用于分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算,生成第二实频向量积分值和第二虚频向量积分值。
第一相位模块1140,用于根据预设的相位转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位。
第二相位模块1150,用于根据所述预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位。
第一相位扩展模块1160,用于根据预设的相位扩展规则,将所述第一相位扩展为第一扩展相位。
第二相位扩展模块1170,用于根据预设的相位扩展规则,将所述第二相位扩展为第二扩展相位。
初相位模块1180,用于根据预设的初相位转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的初相位。
本实施方式,对采样所得的信号序列进行截短处理,获得截短信号序列;以所测参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列和所述截短信号序列相乘,生成两组实频向量序列和虚频向量序列;通过对两组虚频向量序列和实频向量序列数字陷波,生成两组虚数向量陷波序列和实数向量陷波序列,进而积分生成两组虚数向量积分值和实数向量积分值;再根据预设的相位转换规则,将两组实数向量积分值和虚数向量积分值转换为两个相位;再将两个相位进行扩展,获得扩展相位;进而根据预设的初相位转换规则,将两个扩展相位转换为所述正弦信号的初相位。
由于数字陷波可以在某一个频率点迅速衰减输入信号,以达到阻碍此频率信号通过的效果,因此,当该数字陷波的陷波频率点设为对应混频干扰频率点时,该数字陷波对所述混频干扰频率具有完全的抑制作用。如此,采用本发明的正弦信号的初相位检测方法,其对低频正弦信号的初相位测量精度高,且抗谐波和噪声干扰性好。
其中,初步序列长度模块1010,电力系统频率范围在45Hz-55Hz,优选的,取正弦信号频率下限fmin为45Hz;并根据实际需要设置所述预设信号周期数C2π,优选的,所述预设信号周期数C2π取整数。
在一些实施例中,取C2π为整数12。
初步采样序列长度计算,为式(1):
式中,Nstart为所述初步序列长度,单位无量纲;(int)表示取整;C2π为所述预设信号周期数,单位无量纲;fmin为正弦信号频率范围的下限,单位Hz;fn为所述预设采样频率,单位Hz。
对于初步采样模块1020,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列。对单基波频率正弦信号,获得所述正弦信号的初步采样序列为式(2):
其中,Xstart(n)为初步采样序列;A为信号幅值,单位v;ω为信号频率,单位rad/s;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;为信号初相位,单位rad;Nstart为初步采样序列长度,单位无量纲。
对于频率初测模块1030,可通过零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法对所述初步采样序列进行频率初测,获取所述初步频率。
所述初步频率表达为式(3):
ωo (3);
其中,ωo为初步频率,单位rad/s;
优选地,所述参考频率表达为式(4):
ωs=ωo (4);
其中,ωs为参考频率,单位rad/s;ωo为初步频率,单位rad/s。
对于周期序列模块1040,根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度:
在一些实施例中,所述正弦信号的单位周期序列长度计算,为式(5):
式中,N2π为所述单位周期序列长度,单位无量纲;(int)为取整数;fn为预设采样频率,单位Hz;fs为Hz单位的参考频率;ωs为rad/s单位的参考频率。
所述单位周期序列长度整数化存在1个采样间隔内的误差。
对于序列长度模块1050,根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度:
在一些实施例中,所述预设序列长度为所述单位周期序列长度的12倍,所述预设序列长度计算,为式(6):
N=(int)(C2πN2π) (6);
其中,N为预设序列长度,单位无量纲;(int)为取整数;N2π为所述单位周期序列长度,单位无量纲。由于存在误差,所述预设序列长度所包含的信号周期整数是大约的。
对于信号序列模块1060,根据预设序列长度,从所述初步采样序列中,获得信号序列:
在一些实施例中,所述预设序列长度的信号序列,为式(7):
其中,Xi(n)为信号序列;Xstart(n)为初步采样序列;A为信号幅值,单位v;ω为信号频率,单位rad/s;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;为信号初相位,单位rad;N信号序列长度,单位无量纲,信号序列长度等于所述预设序列长度;Nstart为初步采样序列长度,单位无量纲。
所述信号序列的图形表达,如图3所示。
对于截短序列模块1070,对所述信号序列进行截短处理,获得截短信号序列:
在一些实施例中,所述截短信号序列长度为所述单位周期序列长度的11.25倍,所述截短信号序列长度计算为式(8):
NS=N-0.75N2π (8);
式中,NS为截短信号序列长度,单位无量纲;N为信号序列长度,单位无量纲;N2π为信号单位周期序列长度,单位无量纲。系数0.75代表所述截短值为所述单位周期序列长度的0.75倍;
优选地,所述截短信号序列表达为式(9):
式中,X2(n)为截短信号序列;NS为截短信号序列长度,单位无量纲。
所述截短信号序列的图形表达,如图3所示。
对于第一混频模块1080,优选地,所述参考频率的余弦函数和所述参考频率的正弦函数可分别为以所述参考频率为频率、以Tn为间隔离散变量的正弦函数和余弦函数。
在一些实施例中,以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,得到所述第一实频向量序列和所述第一虚频向量序列为式(10):
其中,R1(n)为所述第一实频向量序列;I1(n)为所述第一虚频向量序列;Ω为信号频率与参考频率的频差,单位rad/s;和为有效分量;和为混频干扰频率成分。
对于第二混频模块1090,以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,得到所述第二实频向量序列和所述第二虚频向量序列为式(11):
式中,R2(n)为所述第二实频向量序列;I2(n)为所述第二虚频向量序列;Ω为信号频率与参考频率的频差,单位rad/s;和为有效分量;和为混频干扰频率。
对于第一陷波模块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):
式中,X(n)为四级数字陷波输入序列,输入序列长度N;XD(n)为四级数字陷波输出序列,输出序列长度N-NW1-ND1-NW2-ND2;W1(n)为三角函数1,其中函数峰值为1,零频率增益为0.5;W2(n)为三角函数2,其中函数峰值为1,零频率增益为0.5;NW1为三角窗参数1,即三角窗函数1序列长度,单位无量纲,要求为奇数,以保证三角窗函数的外形为等腰三角形(如图5所示),(int)代表取整数;ND1为矩形窗参数1,即矩形窗函数1序列长度,单位无量纲,(int)代表取整数;NW2为三角窗参数2,即三角窗函数2序列长度,单位无量纲,要求为奇数,以保证三角窗函数的外形为等腰三角形(如图6所示),(int)代表取整数;ND2矩形窗参数2,即矩形窗函数2序列长度,单位无量纲,式中给出了计算公式,(int)代表取整数。
在一些实施例中,三角窗参数NW1取值为所述参考频率的单位周期序列长度的3倍,矩形窗参数ND1取值为所述参考频率的单位周期序列长度的1.5倍,三角窗参数NW2取值为所述参考频率的单位周期序列长度的4倍,矩形窗参数ND1取值为所述参考频率的单位周期序列长度的2倍。四级数字陷波需要使用10.5倍信号周期序列长度。
对上述实施例,在所述正弦信号的基波频率100π、单位rad/s,得到三角窗算术平均陷波算器1的频域特性,如图5所示。得到三角窗算术平均陷波算器2的频域特性,如图6所示。
优选地,在所述混频干扰频率成分得到完全抑制前提下,所述第一实频向量陷波序列和所述第一虚频向量陷波序列为(13):
式中,RD1(n)为所述第一实频向量陷波序列;ID1(n)为所述第一虚频向量陷波序列;K(Ω)为数字陷波在频差Ω的无量纲增益;α(Ω)为数字陷波在频差Ω的移相,单位rad。
对于第一积分模块1110,优选地,可通过本领域技术人员惯用的积分器进行积分运算。
积分运算式为(14):
式中,R1为第一实频向量积分值;I1为第一虚频向量积分值。L1为第一积分长度,单位无量纲,L1为1.5倍所述单位周期序列长度。
对于第二陷波模块1120,同理和优选地,在所述混频干扰频率成分得到完全抑制前提下,所述第二实频向量陷波序列和所述第二虚频向量陷波序列为式(15):
其中,RD2(n)为所述第二实频向量陷波序列;ID2(n)为所述第二虚频向量陷波序列;K(Ω)为数字陷波在频差Ω的无量纲增益;α(Ω)为数字陷波在频差Ω的移相,单位rad。
对于第二积分模块1130,优选地,积分运算式可为(16):
式中,R2为第二实频向量积分值;I2为第二虚频向量积分值。L2为第二积分长度,单位无量纲,L2为0.75倍所述单位周期序列长度。
对于第一相位模块1140,优选地,预设的相位转换规则对应于虚频向量积分值和实频向量转换为相位的转换式,可通过以下公式(17)将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位:
式中,PH1为第一相位,单位rad。
在一些实施例中,第二相位模块1140可用于:
获取所述第一虚频向量积分值与所述第一实频向量积分值的比值;
获取所述比值的反正切函数值的相反数,生成所述第一相位。
对于第二相位模块1150,优选地,可通过以下公式(18)将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位:
式中,PH2为第二相位,单位rad。
在一些实施例中,第二相位模块1150可用于:
获取所述第二虚频向量积分值与所述第二实频向量积分值的比值;
获取所述比值的反正切函数值的相反数,生成所述第二相位。
对于第一相位扩展模块1160,所述第一相位在0~±0.5πrad范围,但实际序列相位可能会超出±0.5πrad范围,因此必须根据相位扩展规则对所述第一相位进行扩展,扩展后的相位范围在0~±πrad范围,第一扩展相位,为式(19):
式中,Ph1为第一扩展相位,范围在0~±πrad;&代表与逻辑。
在一些实施例中,第一相位扩展模块1160可用于:
如果所述第一实频向量积分值大于等于零的同时,所述第一虚频向量积分值的相反数大于等于零,则第一扩展相位等于所述第一相位;
如果所述第一实频向量积分值大于等于零的同时,所述第一虚频向量积分值的相反数小于零,则第一扩展相位等于所述第二相位;
如果所述第一实频向量积分值小于零的同时,所述第一虚频向量积分值的相反数大于等于零,则第一扩展相位等于所述第二相位加πrad;
如果所述第一实频向量积分值小于零的同时,所述第一虚频向量积分值的相反数小于零,则第一扩展相位等于所述第二相位减πrad;
对于第二相位扩展模块1170,所述第二相位在0~±0.5πrad范围,但实际序列相位可能会超出±0.5πrad范围,因此必须根据相位扩展规则对所述第二相位进行扩展,扩展后的相位范围在0~±πrad范围,第二扩展相位,为式(20):
式中,PH2为所述第二相位,范围在0~±0.5πrad;Ph2为第二扩展相位,范围在0~±πrad;&代表与逻辑。
在一些实施例中,第二相位扩展模块1170可用于:
如果所述第二实频向量积分值大于等于零的同时,所述第二虚频向量积分值的相反数大于等于零,则第二扩展相位等于所述第二相位;
如果所述第二实频向量积分值大于等于零的同时,所述第二虚频向量积分值的相反数小于零,则第二扩展相位等于所述第二相位;
如果所述第二实频向量积分值小于零的同时,所述第二虚频向量积分值的相反数大于等于零,则第二扩展相位等于所述第二相位加πrad;
如果所述第二实频向量积分值小于零的同时,所述第二虚频向量积分值的相反数小于零,则第二扩展相位等于所述第二相位减πrad;
对于初相位模块1180,所述预设的初相位转换规则可对应于第一扩展相位和第二扩展相位转换为初相位的公式。根据式(19)和式(20),可生成与所述预设的初相位转换规则对应的初相位公式(21):
式中,为正弦信号初相位检测值,单位rad。
在一些实施例中,初相位模块1180可用于:
获取所述第一扩展相位与所述截短信号序列的序列长度的乘积,生成第一乘积。
获取所述第二扩展相位与所述预设序列长度的乘积,生成第二乘积。
获取所述第一乘积与所述第二乘积的差值,生成第一差值。
获取所述截短信号序列的序列长度与所述预设序列长度的差值,生成第二差值。
请参阅图4,图4是本发明正弦信号的初相位检测方法的初相位的相对误差的实验结果示意图。
为了验证本发明正弦信号的初相位检测系统具有较高的准确度,给出一实验信号,为式(22):
在实验信号基波频率变化范围在45Hz-55Hz,取信号整数周期数为12(由于存在误差,信号序列对应的信号整数周期数大约为12),信号初相位变化范围在0~±π/2、单位rad,信号的采样频率为20kHz,信号的离散数据量化位数24bit,频率初测相对误差<|±0.25%|。得到实验信号的初相位检测相对误差绝对值|PHerr(f)|随信号基波频率f变化特性的实验结果图,图4所示。图4给出的实验信号初相位的检测准确度在10-9量级。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上陈述仅以实施例来进一步说明本发明的技术内容,以便于读者更容易理解,但不代表本发明的实施方式仅限于此,任何依本发明所做的技术延伸或再创造,均受本发明的保护。
Claims (10)
1.一种正弦信号的初相位检测方法,其特征在于,包括以下步骤:
根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度;
根据所述初步采样序列长度,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列;
对所述初步采样序列进行频率初测,生成所述正弦信号的初步频率,以所述初步频率给定参考频率;
根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度;
根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度;
根据预设序列长度,从所述初步采样序列中,获得信号序列;
对所述信号序列进行截短处理,获得截短信号序列;
以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,生成第一实频向量序列和第一虚频向量序列;
以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,生成第二实频向量序列和第二虚频向量序列;
分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波,生成第一实频向量陷波序列和第一虚频向量陷波序列;
分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算,生成第一实频向量积分值和第一虚频向量积分值;
分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波,生成第二实频向量陷波序列和第二虚频向量陷波序列;
分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算,生成第二实频向量积分值和第二虚频向量积分值;
根据预设的相位转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位;所述预设的相位转换规则为:通过公式将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位,其中PH1为第一相位,单位rad;R1为第一实频向量积分值;I1为第一虚频向量积分值;Ω为信号频率与参考频率的频差,单位rad/s;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;L1为第一积分长度,单位无量纲,L1为1.5倍所述单位周期序列长度;α(Ω)为数字陷波在频差Ω的移相,单位rad;为信号初相位,单位rad;N信号序列长度,单位无量纲,信号序列长度等于所述预设序列长度;
根据预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位;所述预设的相位转换规则为:通过公式将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位,其中PH2为第二相位,单位rad;R2为第二实频向量积分值;I2为第二虚频向量积分值;L2为第二积分长度,单位无量纲,L2为0.75倍所述单位周期序列长度;NS为截短信号序列长度,单位无量纲;
根据预设的第一相位扩展规则,对所述第一相位进行扩展,获得第一扩展相位;所述预设的第一相位扩展规则为:通过公式获得第一扩展相位,其中Ph1为第一扩展相位,范围在0~±πard;&代表与逻辑;
根据预设的第二相位扩展规则,对所述第二相位进行扩展,获得第二扩展相位;所述预设的第二相位扩展规则为:通过公式获得第二扩展相位,其中PH2为所述第二相位,范围在0~±0.5πard;Ph2为第二扩展相位,范围在0~±πard;
根据预设的初相位转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的初相位;所述预设的初相位转换规则为:根据公式和公式生成与所述预设的初相位转换规则对应的初相位公式其中,为正弦信号初相位检测值,单位rad;N为预设序列长度,单位无量纲。
2.根据权利要求1所述的正弦信号的初相位检测方法,其特征在于,根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度的步骤包括以下步骤:
通过以下所述公式将所述正弦信号频率的下限和预设采样频率及预设信号周期数转换为所述初步采样序列长度:
其中,Nstart为所述初步采样序列长度,单位无量纲;(int)表示取整;C2π为所述预设信号周期数,单位无量纲;fmin为正弦信号频率范围的下限,单位Hz;fn为所述预设采样频率,单位Hz。
3.根据权利要求1所述的正弦信号的初相位检测方法,其特征在于,所述初步频率通过对所述初步采样序列采用零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法获得。
4.根据权利要求1所述的正弦信号的初相位检测方法,其特征在于,所述数字陷波由二级滑动三角窗算术平均陷波加二级滑动矩形窗算术平均陷波构成。
5.根据权利要求1所述的正弦信号的初相位检测方法,其特征在于,根据预设的相位扩展规则,将所述第一相位或所述第二相位扩展为所述第一扩展相位或所述第二扩展相位的步骤包括以下步骤:
如果所述实频向量积分值大于等于零的同时,所述虚频向量积分值的相反数大于等于零,则所述扩展相位等于所述相位;
如果所述实频向量积分值大于等于零的同时,所述虚频向量积分值的相反数小于零,则所述扩展相位等于所述相位;
如果所述实频向量积分值小于零的同时,所述虚频向量积分值的相反数大于等于零,则所述扩展相位等于所述相位加πrad;
如果所述实频向量积分值小于零的同时,所述虚频向量积分值的相反数小于零,则所述扩展相位等于所述相位减πrad。
6.一种正弦信号的初相位检测系统,其特征在于,包括:
初步序列长度模块,用于根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度;
初步采样模块,用于根据所述初步采样序列长度,对所述正弦信号进行初步采样,获得所述正弦信号的初步采样序列;
频率初测模块,用于对所述初步采样序列进行频率初测,生成所述正弦信号的初步频率,以所述初步频率给定参考频率;
周期序列模块,用于根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度;
序列长度模块,根据所述预设信号周期数和所述单位周期序列长度,获得预设序列长度;
信号序列模块,用于根据所述预设序列长度,从所述初步采样序列中,获得信号序列;
截短序列模块,用于对所述信号序列进行截短处理,获得截短信号序列;
第一混频模块,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘,生成第一实频向量序列和第一虚频向量序列;
第二混频模块,用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘,生成第二实频向量序列和第二虚频向量序列;
第一陷波模块,用于分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波,生成第一实频向量陷波序列和第一虚频向量陷波序列;
第一积分模块,用于分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算,生成第一实频向量积分值和第一虚频向量积分值;
第二陷波模块,用于分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波,生成第二实频向量陷波序列和第二虚频向量陷波序列;
第二积分模块,用于分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算,生成第二实频向量积分值和第二虚频向量积分值;
第一相位模块,用于根据预设的相位转换规则,将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位;所述预设的相位转换规则为:通过公式将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位,其中PH1为第一相位,单位rad;R1为第一实频向量积分值;I1为第一虚频向量积分值;Ω为信号频率与参考频率的频差,单位rad/s;Tn为采样间隔,单位s;n为序列离散数,单位无量纲;L1为第一积分长度,单位无量纲,L1为1.5倍所述单位周期序列长度;α(Ω)为数字陷波在频差Ω的移相,单位rad;为信号初相位,单位rad;N信号序列长度,单位无量纲,信号序列长度等于所述预设序列长度;
第二相位模块,用于根据所述预设的相位转换规则,将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位;所述预设的相位转换规则为:通过公式将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位,其中PH2为第二相位,单位rad;R2为第二实频向量积分值;I2为第二虚频向量积分值;L2为第二积分长度,单位无量纲,L2为0.75倍所述单位周期序列长度;NS为截短信号序列长度,单位无量纲;
第一相位扩展模块,用于根据预设的第一相位扩展规则,将所述第一相位扩展为第一扩展相位;所述预设的第一相位扩展规则为:通过公式获得第一扩展相位,其中Ph1为第一扩展相位,范围在0~±πard;&代表与逻辑;
第二相位扩展模块,用于根据预设的第二相位扩展规则,将所述第二相位扩展为第二扩展相位;所述预设的第二相位扩展规则为:通过公式获得第二扩展相位,其中PH2为所述第二相位,范围在0~±0.5πard;Ph2为第二扩展相位,范围在0~±πard;
初相位模块,用于根据预设的初相位转换规则,将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的初相位;所述预设的初相位转换规则为:根据公式和公式生成与所述预设的初相位转换规则对应的初相位公式其中,为正弦信号初相位检测值,单位rad;N为预设序列长度,单位无量纲。
7.根据权利要求6所述的正弦信号的初相位检测系统,其特征在于,所述初步采样模块根据正弦信号频率范围的下限和预设采样频率及预设信号周期数,获得初步采样序列长度的步骤包括以下步骤:
通过以下所述公式将所述正弦信号频率的下限和预设采样频率及预设信号周期数转换为所述初步采样序列长度:
其中,Nstart为所述初步采样序列长度,单位无量纲;(int)表示取整;C2π为所述预设信号周期数,单位无量纲;fmin为正弦信号频率范围的下限,单位Hz;fn为所述预设采样频率,单位Hz。
8.根据权利要求6所述的正弦信号的初相位检测系统,其特征在于,所述初测模块通过对所述初步采样序列采用零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法来获得所述初步频率。
9.根据权利要求6所述的正弦信号的初相位检测系统,其特征在于,所述第一陷波模块、第二陷波模块包括二滑动三角窗算术平均陷波器和加二级滑动矩形窗算术平均陷波器。
10.根据权利要求6所述的正弦信号的初相位检测系统,其特征在于,所述初相位模块根据预设的相位扩展规则,将所述第一相位或所述第二相位扩展为所第一扩展相位或所述第二扩展相位的步骤包括以下步骤:
如果所述实频向量积分值大于等于零的同时,所述虚频向量积分值的相反数大于等于零,则所述扩展相位等于所述相位;
如果所述实频向量积分值大于等于零的同时,所述虚频向量积分值的相反数小于零,则所述扩展相位等于所述相位;
如果所述实频向量积分值小于零的同时,所述虚频向量积分值的相反数大于等于零,则所述扩展相位等于所述相位加πrad;
如果所述实频向量积分值小于零的同时,所述虚频向量积分值的相反数小于零,则所述扩展相位等于所述相位减πrad。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510400898.6A CN105445552B (zh) | 2015-07-09 | 2015-07-09 | 正弦信号的初相位检测方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510400898.6A CN105445552B (zh) | 2015-07-09 | 2015-07-09 | 正弦信号的初相位检测方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105445552A CN105445552A (zh) | 2016-03-30 |
CN105445552B true CN105445552B (zh) | 2018-07-27 |
Family
ID=55555982
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510400898.6A Active CN105445552B (zh) | 2015-07-09 | 2015-07-09 | 正弦信号的初相位检测方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105445552B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3052559B1 (fr) * | 2016-06-10 | 2020-06-12 | Onera (Office National D'etudes Et De Recherches Aerospatiales) | Systeme et procede pour fournir l'amplitude et le retard de phase d'un signal sinusoidal |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0394206B1 (en) * | 1989-04-20 | 1995-03-15 | Telefonaktiebolaget L M Ericsson | A method and an arrangement for accurate digital determination of the time or phase position of a signal pulse train |
US6646477B1 (en) * | 2002-02-27 | 2003-11-11 | National Semiconductor Corporation | Phase frequency detector with increased phase error gain |
CN101251556A (zh) * | 2008-03-04 | 2008-08-27 | 北京航空航天大学 | 正弦信号四参数检测方法和虚拟仪器信号检测装置 |
CN101576586A (zh) * | 2008-05-09 | 2009-11-11 | 中国科学院半导体研究所 | 一种基于自适应滤波的相位检测方法 |
CN104267257A (zh) * | 2014-09-30 | 2015-01-07 | 广东电网有限责任公司电力科学研究院 | 基于点频滤波器的信号幅频谱检测方法和系统 |
CN104459321A (zh) * | 2014-12-29 | 2015-03-25 | 广东电网有限责任公司电力科学研究院 | 电力信号的基波相位测量方法和系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3460145B2 (ja) * | 2000-03-15 | 2003-10-27 | アンリツ株式会社 | 周波数測定装置 |
-
2015
- 2015-07-09 CN CN201510400898.6A patent/CN105445552B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0394206B1 (en) * | 1989-04-20 | 1995-03-15 | Telefonaktiebolaget L M Ericsson | A method and an arrangement for accurate digital determination of the time or phase position of a signal pulse train |
US6646477B1 (en) * | 2002-02-27 | 2003-11-11 | National Semiconductor Corporation | Phase frequency detector with increased phase error gain |
CN101251556A (zh) * | 2008-03-04 | 2008-08-27 | 北京航空航天大学 | 正弦信号四参数检测方法和虚拟仪器信号检测装置 |
CN101576586A (zh) * | 2008-05-09 | 2009-11-11 | 中国科学院半导体研究所 | 一种基于自适应滤波的相位检测方法 |
CN104267257A (zh) * | 2014-09-30 | 2015-01-07 | 广东电网有限责任公司电力科学研究院 | 基于点频滤波器的信号幅频谱检测方法和系统 |
CN104459321A (zh) * | 2014-12-29 | 2015-03-25 | 广东电网有限责任公司电力科学研究院 | 电力信号的基波相位测量方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105445552A (zh) | 2016-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105067880B (zh) | 对电力信号进行正交调制的方法和系统 | |
CN104502700B (zh) | 电力信号的正弦参数测量方法和系统 | |
CN104635045B (zh) | 基于相位调制的电力信号频率检测方法和系统 | |
CN105445552B (zh) | 正弦信号的初相位检测方法和系统 | |
CN105445551B (zh) | 正弦信号的截止相位检测方法和系统 | |
CN105372471A (zh) | 正弦信号的幅值检测方法和系统 | |
CN105203844B (zh) | 电力信号的零初相位调制方法和系统 | |
CN105445547A (zh) | 正弦信号的频率检测方法和系统 | |
CN105353217B (zh) | 正弦信号的全相位差检测方法和系统 | |
CN105067885B (zh) | 将电力信号转换为零初相位信号序列的方法和系统 | |
CN104991104B (zh) | 电力信号的幅值检测方法和系统 | |
CN105548699B (zh) | 获取电力信号的提高频率基准正弦函数序列的方法和系统 | |
CN105158564B (zh) | 根据正弦函数调制的电力信号全相位差检测方法和系统 | |
CN105182077B (zh) | 根据余弦函数调制的电力信号全相位差检测方法和系统 | |
CN105203843B (zh) | 电力信号的平均初相位检测方法和系统 | |
CN105425031B (zh) | 从电力信号中获取零初相位基准正交倍频序列方法和系统 | |
CN105137172B (zh) | 将电力信号转换为零初相位余弦信号序列的方法和系统 | |
CN105158559B (zh) | 将电力信号转换为零初相位正弦信号序列的方法和系统 | |
CN105158560B (zh) | 从电力信号中提取正交基准信号序列的方法和系统 | |
CN105548690B (zh) | 基于零初相位基准正弦函数倍频序列频率测量方法和系统 | |
CN105137188B (zh) | 从电力信号中提取正交信号序列的方法和系统 | |
CN105223419B (zh) | 电力信号的全相位差检测方法和系统 | |
CN105548694B (zh) | 根据零初相位基准正弦函数序列进行频率测量方法和系统 | |
CN105548700B (zh) | 获取电力信号的提高频率正交序列的方法和系统 | |
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 |