CN105353217A - 正弦信号的全相位差检测方法和系统 - Google Patents

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，分别对所述第四实频向量序列和所述第四虚频向量序列进行数字陷波，生成第四实频向量陷波序列和第四虚频向量陷波序列；

S120，分别对所述第四实频向量陷波序列和所述第四虚频向量陷波序列进行积分运算，生成第四实频向量积分值和第四虚频向量积分值；

S121，根据预设的相位转换规则，将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位；

S122，根据预设的相位转换规则，将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位；

S123，根据预设的相位转换规则，将所述第三虚频向量积分值与所述第三实频向量积分值转换为第三相位；

S124，根据预设的相位转换规则，将所述第四虚频向量积分值与所述第四实频向量积分值转换为第四相位；

S125，根据预设的相位扩展规则，对所述第一相位进行扩展，获得第一扩展相位；

S126，根据预设的相位扩展规则，对所述第二相位进行扩展，获得第二扩展相位；

S127，根据预设的相位扩展规则，对所述第三相位进行扩展，获得第三扩展相位；

S128，根据预设的相位扩展规则，对所述第四相位进行扩展，获得第四扩展相位；

S129，根据预设的截止相位转换规则，将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的截止相位；

S130，根据预设的初相位转换规则，将所述第三扩展相位和所述第四扩展相位转换为所述正弦信号的初相位；

S131，将所述截止相位与所述初相位的差值转换为所述正弦信号的全相位差。

本实施方式，将所述正向信号序列反向输出获得反褶序列，将反褶序列进行截短处理，获得反褶截短序列；以所测参考频率的余弦函数和所述参考频率的正弦函数分别与反褶序列和反褶截短序列相乘，生成两组实频向量序列和虚频向量序列；通过对两组虚频向量序列和实频向量序列数字陷波，生成两组虚数向量陷波序列和实数向量陷波序列，进而积分生成两组虚数向量积分值和实数向量积分值；再根据预设的相位转换规则，将两组实数向量积分值和虚数向量积分值转换为两个相位；再将两个相位进行扩展，获得扩展相位；进而根据预设的相位转换规则，将两个扩展相位转换为所述正弦信号的截止相位，将所述截止相位与所述初相位的差值转换为所述正弦信号的全相位差。

由于数字陷波可以在某一个频率点迅速衰减输入信号，以达到阻碍此频率信号通过的效果，因此，当该数字陷波的陷波频率点设为对应混频干扰频率点时，该数字陷波对所述混频干扰频率具有深度的抑制作用。如此，采用本发明的正弦信号的截止相位检测方法，其对低频正弦信号的全相位差测量精度高，且抗谐波和噪声干扰性好。

之后如不加说明，所述正弦信号的全相位差均指正弦信号基波的全相位差。

其中，对于步骤S101，电力系统频率范围在45Hz-55Hz，优选的，取正弦信号频率下限f_min为45Hz；根据实际需要设置所述预设信号周期数C_2π，优选的，所述预设信号周期数C_2π取整数。

在一些实施例中，取C_2π为整数12。

初步采样序列长度计算，为式(1)：

$N_{start}=\left(\mathrm{int}\right)\frac{C_{2\mathrm{\π}}{f}_n}{f_{\min }}---\left(1\right);$

式中，N_start为所述初步序列长度，单位无量纲；(int)表示取整；C_2π为所述预设信号周期数，单位无量纲；f_min为正弦信号频率范围的下限，单位Hz；f_n为所述预设采样频率，单位Hz。

对于步骤S102，对所述正弦信号进行初步采样，获得所述正弦信号的初步采样序列。对单基波频率正弦信号，获得所述正弦信号的初步采样序列为式(2)：

$T_n=\frac{1}{f_n}$ (2)；

n＝0,1,2,3,.....,N_start-1

其中，X_start(n)为初步采样序列；A为信号幅值，单位v；ω为信号频率，单位rad/s；T_n为采样间隔，单位s；f_n为所述采样频率，单位Hz；n为序列离散数，单位无量纲；为信号初相位，单位rad；N_start为初步采样序列长度，单位无量纲。

对于步骤S103，可通过零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法对所述初步采样序列进行频率初测，获取所述初步频率。

所述初步频率表达为式(3)：

ω_o(3)；

其中，ω_o为初步频率，单位rad/s；

优选地，所述参考频率表达为式(4)：

ω_s＝ω_o(4)；

其中，ω_s为参考频率，单位rad/s；ω_o为初步频率，单位rad/s。

对于步骤S104，根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度：

在一些实施例中，所述正弦信号的单位周期序列长度计算，为式(5)：

$N_{2\mathrm{\π}}=\left(\mathrm{int}\right)\frac{f_n}{f_s}$

$f_s=\frac{{\mathrm{\ω}}_s}{2\mathrm{\π}}---\left(5\right);$

式中，N_2π为所述单位周期序列长度，单位无量纲；(int)为取整数；f_n为预设采样频率，单位Hz；f_s为Hz单位的参考频率；ω_s为rad/s单位的参考频率。

所述单位周期序列长度整数化存在1个采样间隔内的误差。

对于步骤S105，根据所述预设信号周期数和所述单位周期序列长度，获得预设序列长度：

在一些实施例中，所述预设序列长度为所述单位周期序列长度的12倍，所述预设序列长度计算，为式(6)：

N＝(int)(C_2πN_2π)(6)；

其中，N为预设序列长度，单位无量纲；(int)为取整数；为所述单位周期序列长度，单位无量纲。由于存在误差，所述预设序列长度所包含的信号周期整数是大约的。

对于步骤S106，根据预设序列长度，从所述初步采样序列中，获得信号序列：

在一些实施例中，所述预设序列长度的信号序列，为式(7)：

n＝0,1,2,3,.....,N-1(7)；

N<N_start

其中，X_i(n)为信号序列；X_start(n)为初步采样序列；A为信号幅值，单位v；ω为信号频率，单位rad/s；T_n为采样间隔，单位s；n为序列离散数，单位无量纲；为信号初相位，单位rad；N信号序列长度，单位无量纲，信号序列长度等于所述预设序列长度；N_start为初步采样序列长度，单位无量纲。

所述信号序列的图形表达，如图3所示。

对于步骤S107，相对正向信号序列，反褶序列表达为式(8)：

X_-i(n)＝X_i(N-n)＝Acos(-ωT_nn+β)

(8)；

n＝0，1，2，3，......，N-1

式中，X_-i(n)为反褶序列；β为反褶序列初相位，单位rad。关系上，反褶序列初相位是正向信号序列的截止相位，即所述正弦信号的截止相位；N为反褶序列长度，单位无量纲。反褶序列长度与正向信号序列长度相同。

所述反褶序列的图形表达，如图3所示。

对于步骤S108，分别将所述反褶序列和所述正向信号序列进行截短，获得序列长度相同的反褶截短序列和正向截短序列；

在一些实施例中，所述反褶截短序列长度为所述单位周期序列长度的11.25倍，所述反褶截短序列长度计算为式(9)：

N_S＝N-0.75N_2π(9)；

式中，N_S为反褶截短序列长度，单位无量纲；N为正向信号序列长度，单位无量纲；N_2π为信号单位周期序列长度，单位无量纲。系数0.75代表所述截短值为所述单位周期序列长度的0.75倍；

优选地，所述反褶截短序列表达为式(10)：

X_-2(n)＝Acos(-ωT_nn+β)

n＝0,1,2,3,.....,N_s-1(10)；

式中，X_-2(n)为截短信号序列；N_S为反褶截短序列长度，单位无量纲。

所述反褶截短序列的图形表达，如图3所示。

将所述正向信号序列进行截短，获得正向截短序列，所述正向截短序列表达为式(11)：

(11)；

n＝0,1,2,3,.....,N_s-1

式中，X₂(n)为正向截短序列；N_s为正向截短序列长度，单位无量纲。

所述正向截短序列的长度与所述反褶截短序列的长度相同，表达为式(9)。

对于步骤S109，优选地，所述参考频率的余弦函数和所述参考频率的正弦函数可分别为以所述参考频率为频率、以T_n为间隔离散变量的正弦函数和余弦函数。

在一些实施例中，以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述反褶序列相乘，得到所述第一实频向量序列和所述第一虚频向量序列为式(12)：

$\begin{array}{c}{R}_1\left(n\right)=Acos(-{\mathrm{\&omega;T}}_{n}n+\mathrm{\&beta;})cos(-{\mathrm{\&omega;}}_s{T}_{n}n)\\ =\frac{A}{2}\cos (-{\mathrm{\&Omega;T}}_{n}n+\mathrm{\&beta;})+\frac{A}{2}cos\&lsqb;-(\mathrm{\&omega;}+{\mathrm{\&omega;}}_s)T_{n}n+\mathrm{\&beta;}\&rsqb;\end{array}$

$\begin{array}{c}{I}_1\left(n\right)=Acos(-{\mathrm{\&omega;T}}_{n}n+\mathrm{\&beta;})\sin (-{\mathrm{\&omega;}}_s{T}_{n}n)\\ =-\frac{A}{2}\sin (-{\mathrm{\&Omega;T}}_{n}n+\mathrm{\&beta;})+\frac{A}{2}\sin \&lsqb;-(\mathrm{\&omega;}+{\mathrm{\&omega;}}_s)T_{n}n+\mathrm{\&beta;}\&rsqb;\end{array}---\left(12\right);$

Ω＝ω-ω_s

n＝0,1,2,3,.....,N-1

式中，将参考频率ω_s乘以负数；R₁(n)为第一实频向量序列，I₁(n)为第一虚频向量序列，Ω为信号频率与参考频率的频差，单位rad/s；N为反褶序列长度，单位无量纲；Acos(-ΩT_nn+β)/2和Asin(-ΩT_nn+β)/2为有效分量；Acos[-(ω+ω_s)T_nn+β]/2和Asin[-(ω+ω_s)T_nn+β]/2为混频干扰频率成分。

对于步骤S110，优选地，所述参考频率的余弦函数和所述参考频率的正弦函数可分别为以所述参考频率为频率、以T_n为间隔离散变量的正弦函数和余弦函数。

以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘，得到所述第二实频向量序列和所述第二虚频向量序列为式(13)：

$\begin{array}{c}{R}_2\left(n\right)=Acos(-{\mathrm{\&omega;T}}_{n}n+\mathrm{\&beta;})cos(-{\mathrm{\&omega;}}_s{T}_{n}n)\\ =\frac{A}{2}\cos (-{\mathrm{\&Omega;T}}_{n}n+\mathrm{\&beta;})+\frac{A}{2}cos\&lsqb;-(\mathrm{\&omega;}+{\mathrm{\&omega;}}_s)T_{n}n+\mathrm{\&beta;}\&rsqb;\end{array}$

$\begin{array}{c}{I}_2\left(n\right)=Acos(-{\mathrm{\&omega;T}}_{n}n+\mathrm{\&beta;})\sin (-{\mathrm{\&omega;}}_s{T}_{n}n)\\ =-\frac{A}{2}\sin (-{\mathrm{\&Omega;T}}_{n}n+\mathrm{\&beta;})+\frac{A}{2}\sin \&lsqb;-(\mathrm{\&omega;}+{\mathrm{\&omega;}}_s)T_{n}n+\mathrm{\&beta;}\&rsqb;\end{array}---\left(13\right);$

Ω＝ω-ω_s

n＝0,1,2,3,.....,N_s-1

式中，将参考频率ω_s乘以负数；R₂(n)为第二实频向量序列，I₂(n)为第二虚频向量序列，Ω为信号频率与参考频率的频差，单位rad/s；N_s为反褶截短序列长度，单位无量纲；Acos(-ΩT_nn+β)/2和Asin(-ΩT_nn+β)/2为有效分量；Acos[-(ω+ω_s)T_nn+β]/2和Asin[-(ω+ω_s)T_nn+β]/2为混频干扰频率成分。

对于步骤S111，优选地，所述参考频率的余弦函数和所述参考频率的正弦函数可分别为以所述参考频率为频率、以T_n为间隔离散变量的正弦函数和余弦函数。

在一些实施例中，以所述余弦函数和所述正弦函数分别与所述正向信号序列相乘，生成第三实频向量序列和第三虚频向量序列，为式(14)：

Ω＝ω-ω_s

n＝1,2,3,.....,N-1

其中，R₃(n)为所述第三实频向量序列；I₃(n)为所述第三虚频向量序列；Ω为信号频率与参考频率的频差，单位rad/s；和为有效分量；和为混频干扰频率成分。

对于步骤S112，优选地，所述参考频率的余弦函数和所述参考频率的正弦函数可分别为以所述参考频率为频率、以T_n为间隔离散变量的正弦函数和余弦函数。

在一些实施例中，以所述余弦函数和所述正弦函数分别与所述正向截短序列相乘，生成第四实频向量序列和第四虚频向量序列，为式(15)：

Ω＝ω-ω_s

n＝1,2,3,.....,N_S-1

式中，R₄(n)为所述第四实频向量序列；I₄(n)为所述第四虚频向量序列；Ω为信号频率与参考频率的频差，单位rad/s；和为有效分量；和为混频干扰频率。

对于步骤S113，所述第一实频向量序列和所述第一虚频向量序列中包含混频干扰频率。当输入信号中包含直流成分、次谐波成分及分次谐波成分时，所述混频干扰频率将更加复杂，这些混频干扰频率严重影响计算准确度。虽然窗口函数和积分运算本身对混频干扰频率具有良好的衰减作用，但没有针对性，不能够对所述复杂的混频干扰频率产生深度的抑制作用，不能满足参数的高准确度计算需要。为了有针对性的抑制所述混频干扰频率的影响，可以对该混频干扰进行数字陷波器，由于数字陷波可以在某一个频率点迅速衰减输入信号，以达到阻碍此频率信号通过的效果，因此，当该数字陷波的陷波频率点设为对应混频干扰频率点时，该数字陷波对所述混频干扰频率具有深度的抑制作用。

优选地，数字陷波具体采用滑动三角窗算术平均陷波，即将若干个连续离散值乘以三角窗函数后再相加，然后取其算术平均值作为本次陷波值输出。滑动三角窗算术平均陷波器需要设置三角窗参数，所述三角窗参数具体指三角窗函数序列的长度N_W。在三角窗参数N_W值为信号周期序列长度的3倍，可以对1/3分次谐波产生的混频干扰频率进行抑制。在三角窗参数N_W值为信号周期序列长度的4倍，可以对直流、1/2分次、1次、2次、3次、4次、5次谐波等产生的混频干扰频率进行抑制。

考虑到实际存在误差等因数，如所述参数存在1个采样间隔以内的误差，为了深度抑制混频干扰频率影响，还可再增加一级滑动矩形窗算术平均陷波器，即将若干个连续离散值直接相加，然后取其算术平均值作为本次陷波值输出。滑动矩形窗算术平均陷波器需要设置矩形窗参数，所述矩形窗参数具体指矩形窗函数序列的长度N_D。矩形窗参数N_D取值为信号周期序列长度的1.5倍，可以对1/3分次谐波产生的混频干扰频率进行抑制。而N_D取值为信号周期序列长度的2倍，可以对直流、1/2分次、1次、2次、3次、4次、5次谐波等产生的混频干扰频率进行抑制。

优选地，数字陷波包括二级滑动三角窗算术平均陷波和二级滑动矩形窗算术平均陷波，该四级数字陷波式可为式(16)：

$\begin{array}{c}{X}_D\left(n\right)=\\ \frac{1}{N_{D2}}\underset{n}{\overset{N_{D2}-1}{\&Sigma;}}\frac{2}{N_{W2}-1}\underset{n}{\overset{N_{W2}-1}{\&Sigma;}}{W}_2\left(n\right)\&lsqb;\frac{1}{N_{D1}}\underset{n}{\overset{N_{D1}-1}{\&Sigma;}}\frac{2}{N_{W1}-1}\underset{n}{\overset{N_{W1}-1}{\&Sigma;}}{W}_1\left(n\right)X\left(n\right)\&rsqb;\end{array}$

N_D1＝(int)(1.5N_2π+0.5)

N_D2＝(int)(2N_2π+0.5)

对X(n)n＝0,1,2,3,....,N-1

对W₁(n)n＝0,1,2,3,....,N_W1-1

对W₂(n)n＝0,1,2,3,....,N_W2-1

对X_D(n)n＝0,1,2,3,....,N-N_W1-N_D1-N_W2-N_D2-1

式中，X(n)为四级数字陷波输入序列，输入序列长度N；X_D(n)为四级数字陷波输出序列，输出序列长度N-N_W1-N_D1-N_W2-N_D2；W₁(n)为三角函数1，其中函数峰值为1，零频率增益为0.5；W₂(n)为三角函数2，其中函数峰值为1，零频率增益为0.5；N_W1为三角窗参数1，即三角窗函数1序列长度，单位无量纲，要求为奇数，以保证三角窗函数的外形为等腰三角形(如图5所示)，(int)代表取整数；N_D1为矩形窗参数1，即矩形窗函数1序列长度，单位无量纲，(int)代表取整数；N_W2为三角窗参数2，即三角窗函数2序列长度，单位无量纲，要求为奇数，以保证三角窗函数的外形为等腰三角形(如图6所示)，(int)代表取整数；N_D2矩形窗参数2，即矩形窗函数2序列长度，单位无量纲，式中给出了计算公式，(int)代表取整数。

在一些实施例中，三角窗参数N_W1取值为所述参考频率的单位周期序列长度的3倍，矩形窗参数N_D1取值为所述参考频率的单位周期序列长度的1.5倍，三角窗参数N_W2取值为所述参考频率的单位周期序列长度的4倍，矩形窗参数N_D1取值为所述参考频率的单位周期序列长度的2倍。四级数字陷波需要使用10.5倍信号周期序列长度。

对上述实施例，在所述正弦信号的基波频率100π、单位rad/s，得到三角窗算术平均陷波算器1的频域特性，如图5所示。得到三角窗算术平均陷波算器2的频域特性，如图6所示。

优选地，在所述混频干扰频率成分得到完全抑制前提下，所述第一实频向量陷波序列和所述第一虚频向量陷波序列为(17)：

$\begin{array}{c}{R}_{D1}\left(n\right)=\frac{1}{N_{D2}}\underset{n}{\overset{N_{D2}-1}{\&Sigma;}}\frac{2}{N_{W2}-1}\underset{n}{\overset{N_{W2}-1}{\&Sigma;}}{W}_2\left(n\right)\&lsqb;\frac{1}{N_{D1}}\underset{n}{\overset{N_{D1}-1}{\&Sigma;}}\frac{2}{N_{W1}-1}\underset{n}{\overset{N_{W1}-1}{\&Sigma;}}{W}_1\left(n\right)R_1\left(n\right)\&rsqb;\\ =\frac{AK\left(\mathrm{\&Omega;}\right)}{2}\cos \&lsqb;-{\mathrm{\&Omega;T}}_{n}n+\mathrm{\&beta;}-\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)\&rsqb;\end{array}$

$\begin{array}{c}{I}_{D1}\left(n\right)=\frac{1}{N_{D2}}\underset{n}{\overset{N_{D2}-1}{\&Sigma;}}\frac{2}{N_{W2}-1}\underset{n}{\overset{N_{W2}-1}{\&Sigma;}}{W}_2\left(n\right)\&lsqb;\frac{1}{N_{D1}}\underset{n}{\overset{N_{D1}-1}{\&Sigma;}}\frac{2}{N_{W1}-1}\underset{n}{\overset{N_{W1}-1}{\&Sigma;}}{W}_1\left(n\right)I_1\left(n\right)\&rsqb;\\ =-\frac{AK\left(\mathrm{\&Omega;}\right)}{2}\sin \&lsqb;-{\mathrm{\&Omega;T}}_{n}n+\mathrm{\&beta;}-\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)\&rsqb;\end{array}---\left(17\right);$

$\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)=\frac{{\mathrm{\&Omega;T}}_n(N_{W1}+N_{D1}+N_{W2}+N_{D2})}{2}$

对R₁(n)I₁(n)n＝0,1,2,3,....,N-1

对R_D1(n)I_D1(n)n＝0,1,2,3,....,N-N_W1-N_D1-N_W2-N_D2-1

式中，R_D1(n)为所述第一实频向量陷波序列；I_D1(n)为所述第一虚频向量陷波序列；K(Ω)为数字陷波在频差Ω的无量纲增益；α(Ω)为数字陷波在频差Ω的移相，单位rad。

对于步骤S114，优选地，可通过本领域技术人员惯用的积分器进行积分运算。

积分运算式为(18)：

$R_1=\frac{2}{L1}\underset{0}{\overset{L1-1}{\&Sigma;}}{R}_{D1}\left(n\right)=\frac{2AK\left(\mathrm{\&Omega;}\right)}{{\mathrm{\&Omega;T}}_{n}L1}sin\&lsqb;\frac{{\mathrm{\&Omega;T}}_{n}L1}{2}\&rsqb;cos\&lsqb;-\frac{{\mathrm{\&Omega;T}}_{n}L1}{2}+\mathrm{\&beta;}-\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)\&rsqb;$

$I_1=\frac{2}{L1}\underset{0}{\overset{L1-1}{\&Sigma;}}{I}_{D1}\left(n\right)=-\frac{2AK\left(\mathrm{\&Omega;}\right)}{{\mathrm{\&Omega;T}}_{n}L1}sin\&lsqb;\frac{{\mathrm{\&Omega;T}}_{n}L1}{2}\&rsqb;sin\&lsqb;-\frac{{\mathrm{\&Omega;T}}_{n}L1}{2}+\mathrm{\&beta;}-\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)\&rsqb;---\left(18\right);$

n＝0,1,2,3,.......,L1-1

L1＝N-N_W1-N_D1-N_W2-N_D2

式中，R₁为第一实频向量积分值；I₁为第一虚频向量积分值。L1为第一积分长度，单位无量纲，L1为1.5倍所述单位周期序列长度。

对于步骤S115，同理和优选地，在所述混频干扰频率成分得到完全抑制前提下，所述第二实频向量陷波序列和所述第二虚频向量陷波序列为式(19)：

$\begin{array}{c}{R}_{D2}\left(n\right)=\frac{1}{N_{D2}}\underset{n}{\overset{N_{D2}-1}{\&Sigma;}}\frac{2}{N_{W2}-1}\underset{n}{\overset{N_{W2}-1}{\&Sigma;}}{W}_2\left(n\right)\&lsqb;\frac{1}{N_{D1}}\underset{n}{\overset{N_{D1}-1}{\&Sigma;}}\frac{2}{N_{W1}-1}\underset{n}{\overset{N_{W1}-1}{\&Sigma;}}{W}_1\left(n\right)R_2\left(n\right)\&rsqb;\\ =\frac{AK\left(\mathrm{\&Omega;}\right)}{2}\cos \&lsqb;-{\mathrm{\&Omega;T}}_{n}n+\mathrm{\&beta;}-\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)\&rsqb;\end{array}$

$\begin{array}{c}{I}_{D2}\left(n\right)=\frac{1}{N_{D2}}\underset{n}{\overset{N_{D2}-1}{\&Sigma;}}\frac{2}{N_{W2}-1}\underset{n}{\overset{N_{W2}-1}{\&Sigma;}}{W}_2\left(n\right)\&lsqb;\frac{1}{N_{D1}}\underset{n}{\overset{N_{D1}-1}{\&Sigma;}}\frac{2}{N_{W1}-1}\underset{n}{\overset{N_{W1}-1}{\&Sigma;}}{W}_1\left(n\right)R_2\left(n\right)\&rsqb;\\ =-\frac{AK\left(\mathrm{\&Omega;}\right)}{2}\sin \&lsqb;-{\mathrm{\&Omega;T}}_{n}n+\mathrm{\&beta;}-\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)\&rsqb;\end{array}---\left(19\right);$

$\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)=\frac{{\mathrm{\&Omega;T}}_n(N_{W1}+N_{D1}+N_{W2}+N_{D2})}{2}$

对R₂(n)I₂(n)n＝0,1,2,3,....,N_s-1

对R_D2(n)I_D2(n)n＝0,1,2,3,....,N_s-N_W1-N_D1-N_W2-N_D2-1

其中，R_D2(n)为所述第二实频向量陷波序列；I_D2(n)为所述第二虚频向量陷波序列；K(Ω)为数字陷波在频差Ω的无量纲增益；α(Ω)为数字陷波在频差Ω的移相，单位rad。

对于步骤S116，优选地，积分运算式可为(20)：

$R_2=\frac{2}{L1}\underset{0}{\overset{L2-1}{\&Sigma;}}{R}_{D2}\left(n\right)=\frac{2AK\left(\mathrm{\&Omega;}\right)}{{\mathrm{\&Omega;T}}_{n}L2}sin\&lsqb;\frac{{\mathrm{\&Omega;T}}_{n}L2}{2}\&rsqb;cos\&lsqb;-\frac{{\mathrm{\&Omega;T}}_{n}L2}{2}+\mathrm{\&beta;}-\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)\&rsqb;$

$I_2=\frac{2}{L2}\underset{0}{\overset{L2-1}{\&Sigma;}}{I}_{D2}\left(n\right)=-\frac{2AK\left(\mathrm{\&Omega;}\right)}{{\mathrm{\&Omega;T}}_{n}L2}sin\&lsqb;\frac{{\mathrm{\&Omega;T}}_{n}L2}{2}\&rsqb;sin\&lsqb;-\frac{{\mathrm{\&Omega;T}}_{n}L2}{2}+\mathrm{\&beta;}-\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)\&rsqb;---\left(20\right);$

n＝0,1,2,3,.......,L2-1

L2＝N_S-N_W1-N_D1-N_W2-N_D2

式中，R₂为第二实频向量积分值；I₂为第二虚频向量积分值。L2为第二积分长度，单位无量纲，L2为0.75倍所述单位周期序列长度。

对于步骤S117，同理和优选地，在所述混频干扰频率成分得到完全抑制前提下，所述第三实频向量陷波序列和所述第三虚频向量陷波序列为式(21)：

$\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)=\frac{{\mathrm{\&Omega;T}}_n(N_{W1}+N_{D1}+N_{W2}+N_{D2})}{2}$

对R₃(n)I₃(n)n＝0,1,2,3,....,N-1

对R_D3(n)I_D3(n)n＝0,1,2,3,....,N-N_W1-N_D1-N_W2-N_D2-1

其中，R_D3(n)为所述第三实频向量陷波序列；I_D3(n)为所述第三虚频向量陷波序列；K(Ω)为数字陷波在频差Ω的无量纲增益；α(Ω)为数字陷波在频差Ω的移相，单位rad。

对于步骤S118，优选地，可通过本领域技术人员惯用的积分器进行积分运算。

积分运算式为(22)：

n＝0,1,2,3,.......,L3-1

L3＝N-N_W1-N_D1-N_W2-N_D2

其中，R₃为第三实频向量积分值；I₃为第三虚频向量积分值。L3为积分计算长度3，单位无量纲，原则上，L3最小为1.5倍信号周期序列长度。

对于步骤S119，同理和优选地，在所述混频干扰频率成分得到完全抑制前提下，所述第四实频向量陷波序列和所述第四虚频向量陷波序列为式(23)：

$\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)=\frac{{\mathrm{\&Omega;T}}_n(N_{W1}+N_{D1}+N_{W2}+N_{D2})}{2}$

对R₄(n)I₄(n)n＝0,1,2,3,....,N_s-1

对R_D4(n)I_D4(n)n＝0,1,2,3,....,N_s-N_W1-N_D1-N_W2-N_D2-1

其中，R_D4(n)为所述第四实频向量陷波序列；I_D4(n)为所述第四虚频向量陷波序列；K(Ω)为数字陷波在频差Ω的无量纲增益；α(Ω)为数字陷波在频差Ω的移相，单位rad。

对于步骤S120，优选地，积分运算式可为(24)：

n＝0,1,2,3,.......,L4-1

L4＝N_S-N_W1-N_D1-N_W2-N_D2

其中，R₄为第四实频向量积分值；I₄为第四虚频向量积分值。L4为积分计算长度4，单位无量纲，原则上，L4最小为0.75倍信号周期序列长度。

对于步骤S121-S124，优选地，预设的相位转换规则对应于虚频向量积分值和实频向量转换为相位的转换式。

优选地，可通过以下公式(25)-(28)将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位，将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位，将所述第三虚频向量积分值与所述第三实频向量积分值转换为第三相位，将所述第四虚频向量积分值与所述第四实频向量积分值转换为第四相位：

$\begin{array}{c}{\mathrm{PH}}_1=-arctan\frac{I_1}{R_1}=-\frac{{\mathrm{\&Omega;T}}_{n}L1}{2}+\mathrm{\&beta;}-\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)=-\frac{{\mathrm{\&Omega;T}}_{n}N}{2}+\mathrm{\&beta;}\\ \mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)=\frac{{\mathrm{\&Omega;T}}_n(N_{W1}+N_{D1}+N_{W2}+N_{D2})}{2}\end{array}---\left(25\right);$

$\begin{array}{c}{\mathrm{PH}}_2=-arctan\frac{I_2}{R_2}=-\frac{{\mathrm{\&Omega;T}}_{n}L2}{2}+\mathrm{\&beta;}-\mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)=-\frac{{\mathrm{\&Omega;T}}_n{N}_s}{2}+\mathrm{\&beta;}\\ \mathrm{\&alpha;}\left(\mathrm{\&Omega;}\right)=\frac{{\mathrm{\&Omega;T}}_n(N_{W1}+N_{D1}+N_{W2}+N_{D2})}{2}\end{array}---\left(26\right);$

其中，PH₁为第一相位，单位rad；R₁为第一实频向量积分值；I₁为第一虚频向量积分值，PH₂为第二相位，单位rad；R₂为第二实频向量积分值；I₂为第二虚频向量积分值，PH₃为第三相位，单位rad；R₃为第三实频向量积分值；I₃为第三虚频向量积分值，PH₄为第四相位，单位rad；R₄为第四实频向量积分值；I₄为第四虚频向量积分值。

在一些实施例中，根据预设的相位转换规则，将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位的步骤包括以下步骤：

获取所述第一虚频向量积分值与所述第一实频向量积分值的比值；

获取所述比值的反正切函数值的相反数，生成所述第一相位。

在一些实施例中，根据所述预设的相位转换规则，将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位的步骤包括以下步骤：

获取所述第二虚频向量积分值与所述第二实频向量积分值的比值；

获取所述比值的反正切函数值的相反数，生成所述第二相位。

类似地，所述第三虚频向量积分值与所述第三实频向量积分值转换为第三相位的步骤、所述第四虚频向量积分值与所述第四实频向量积分值转换为第四相位的步骤亦采用上述预设的相位转换规则。

对于步骤S125，所述第一相位在0～±0.5πrad范围，但实际序列相位可能会超出±0.5πrad范围，因此必须根据相位扩展规则对所述第一相位进行扩展，扩展后的相位范围在0～±πrad范围，第一扩展相位，为式(29)：

式中，Ph₁为第一扩展相位，范围在0～±πrad；&代表与逻辑。

在一些实施例中，根据预设的相位扩展规则，对所述第一相位进行扩展，获得第一扩展相位的步骤包括以下步骤：

如果所述第一实频向量积分值大于等于零的同时，所述第一虚频向量积分值的相反数大于等于零，则第一扩展相位等于所述第一相位；

如果所述第一实频向量积分值大于等于零的同时，所述第一虚频向量积分值的相反数小于零，则第一扩展相位等于所述第二相位；

如果所述第一实频向量积分值小于零的同时，所述第一虚频向量积分值的相反数大于等于零，则第一扩展相位等于所述第二相位加πrad；

如果所述第一实频向量积分值小于零的同时，所述第一虚频向量积分值相的反数小于零，则第一扩展相位等于所述第二相位减πrad。

同理，对于步骤S126-S128，所述第二相位、所述第三相位、所述第四相位采用上述所述相位扩展规则进行扩展，扩展后的相位范围在0～±πrad范围，第二扩展相位、第三扩展相位、第四扩展相位分别为式(30)-(32)：

式中，PH₂为所述第二相位，范围在0～±0.5πrad；Ph₂为第二扩展相位，范围在0～±πrad；&代表与逻辑。

式中，PH₃为所述第三相位，范围在0～±0.5πrad；Ph₃为第三扩展相位，范围在0～±πrad；&代表与逻辑。

式中，PH₄为所述第四相位，范围在0～±0.5πrad；Ph₄为第四扩展相位，范围在0～±πrad；&代表与逻辑。

同理，根据预设的相位扩展规则，对所述第二相位、第三相位、第四相位进行扩展，获得第二扩展相位、第三扩展相位、第四扩展相位的步骤与获得第一扩展相位的步骤一致。

对于步骤S129，所述预设的截止相位转换规则可对应于第一扩展相位和第二扩展相位转换为截止相位的公式。根据式(29)和式(30)，可生成与所述预设的截止相位转换规则对应的截止相位公式(33)：

${\mathrm{PH}}_{\mathrm{\&beta;}}=\frac{{\mathrm{PH}}_1{N}_S-{\mathrm{PH}}_{2}N}{N_S-N}=\mathrm{\&beta;}---\left(33\right);$

式中，PH_β为正弦信号的截止相位检测值，单位rad。

在一些实施例中，根据预设的截止相位转换规则，将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的截止相位的步骤包括以下步骤：

获取所述第一扩展相位与所述反褶截短序列长度的乘积，生成第一乘积；

获取所述第二扩展相位与所述反褶序列长度的乘积，生成第二乘积；

获取所述第一乘积与所述第二乘积的差值，生成第一差值；

获取所述反褶截短序列的长度与所述反褶序列的长度的差值，生成第二差值；

获取所述第一差值与所述第二差值的比值，生成所述截止相位。

对于步骤S130，所述预设的初相位转换规则可对应于第三扩展相位和第四扩展相位转换为初相位的公式。根据式(31)和式(32)，可生成与所述预设的初相位转换规则对应的初相位公式(34)：

式中，为正弦信号初相位检测值，单位rad。

在一些实施例中，根据预设的初相位转换规则，将所述第三扩展相位和所述第扩展四相位转换为所述正弦信号的初相位的步骤包括以下步骤：

获取所述预设序列长度与信号周期序列长度的0.75倍的差值，生成所述截短信号序列的序列长度。

获取所述第三扩展相位与所述正向截短信号序列的长度的乘积，生成第三乘积。

获取所述第四扩展相位与所述正向信号序列的长度的乘积，生成第四乘积。

获取所述第三乘积与所述第四乘积的差值，生成第三差值。

获取所述正向截短信号序列的长度与所述正向信号序列的长度的差值，生成第四差值。

获取所述第三差值与所述第四差值的比值，生成所述初相位。

对于步骤S131，优选地，根据所述正弦信号的截止相位、所述正弦信号的初相位、预设信号周期数，所述正弦信号的全相位差公式为式(35)：

式中，△PH正弦信号的全相位差检测值，单位rad；PH_β为所述正弦信号的截止相位检测值，单位rad；为所述正弦信号的初相位检测值，单位rad；C_2π为所述预设信号周期数，单位无量纲。

本发明还公开了正弦信号的频率检测系统，在一些实施例中，如图2所示，该正弦信号的频率检测系统，包括：初步序列长度模块1010、初步采样模块1020、频率初测模块1030、周期序列模块1040、序列长度模块1050、信号序列模块1060、反褶模块1070、截短模块1080、第一混频模块1090、第二混频模块1100、第三混频模块1110、第四混频模块1120、第一陷波模块1130、第一积分模块1140、第二陷波模块1150、第二积分模块1160、第三陷波模块1170、第三积分模块1180、第四陷波模块1190、第四积分模块1200、第一相位模块1210、第二相位模块1220、第三相位模块1230、第四相位模块1240、第一相位扩展模块1250、第二相位扩展模块1260、第三相位扩展模块1270、第四相位扩展模块1280、截止相位模块1290、初相位模块1300、全相位差模块1310、频率模块1320。其中：

初步序列长度模块1010，用于根据正弦信号频率范围的下限和预设采样频率及预设信号周期数，获得初步采样序列长度；

初步采样模块1020，用于根据所述初步采样序列长度，对所述正弦信号进行初步采样，获得所述正弦信号的初步采样序列；

频率初测模块1030，用于对所述初步采样序列进行频率初测，生成所述正弦信号的初步频率，以所述初步频率给定参考频率；

周期序列模块1040，用于根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度；

序列长度模块1050，根据所述预设信号周期数和所述单位周期序列长度，获得预设序列长度；

正向信号序列模块1060，用于根据所述预设序列长度，从所述初步采样序列中，获得正向信号序列；

反褶模块1070，用于将所述正向信号序列反向输出，获得所述正向信号序列的反褶序列；

截短模块1080，用于分别将所述反褶序列和所述正向信号序列进行截短，获得序列长度相同的反褶截短序列和正向截短序列。

第一混频模块1090，用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述反褶序列相乘，生成第一实频向量序列和第一虚频向量序列；

第二混频模块1100，用于以所述参考频率的余弦函数和所述正弦函数分别与所述反褶截短序列相乘，生成第二实频向量序列和第二虚频向量序列；

第三混频模块1110，用于以所述余弦函数和所述正弦函数分别与所述正向信号序列相乘，生成第三实频向量序列和第三虚频向量序列；

第四混频模块1120，用于以所述余弦函数和所述正弦函数分别与所述正向截短序列相乘，生成第四实频向量序列和第四虚频向量序列；

第一陷波模块1130，用于分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波，生成第一实频向量陷波序列和第一虚频向量陷波序列；

第一积分模块1140，用于分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算，生成第一实频向量积分值和第一虚频向量积分值；

第二陷波模块1150，用于分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波，生成第二实频向量陷波序列和第二虚频向量陷波序列；

第二积分模块1160，用于分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算，生成第二实频向量积分值和第二虚频向量积分值；

第三陷波模块1170，用于分别对所述第三实频向量序列和所述第三虚频向量序列进行数字陷波，生成第三实频向量陷波序列和第三虚频向量陷波序列；

第三积分模块1180，用于分别对所述第三实频向量陷波序列和所述第三虚频向量陷波序列进行积分运算，生成第三实频向量积分值和第三虚频向量积分值；

第四陷波模块1190，分别对所述第四实频向量序列和所述第四虚频向量序列进行数字陷波，生成第四实频向量陷波序列和第四虚频向量陷波序列；

第四积分模块1200，用于分别对所述第四实频向量陷波序列和所述第四虚频向量陷波序列进行积分运算，生成第四实频向量积分值和第四虚频向量积分值；

第一相位模块1210，用于根据预设的相位转换规则，将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位；

第二相位模块1220，用于将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位；

第三相位模块1230，用于将所述第三虚频向量积分值与所述第三实频向量积分值转换为第三相位；

第四相位模块1240，用于将所述第四虚频向量积分值与所述第四实频向量积分值转换为第四相位；

第一相位扩展模块1250，用于根据预设的相位扩展规则，对所述第一相位进行扩展，获得第一扩展相位；

第二相位扩展模块1260，用于根据预设的相位扩展规则，对所述第二相位进行扩展，获得第二扩展相位；

第三相位扩展模块1270，用于根据预设的相位扩展规则，对所述第三相位进行扩展，获得第三扩展相位；

第四相位扩展模块1280，用于根据预设的相位扩展规则，对所述第四相位进行扩展，获得第四扩展相位；

截止相位模块1290，用于根据预设的截止相位转换规则，将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的截止相位；

初相位模块1300，用于根据预设的初相位转换规则，将所述第三扩展相位和所述第四扩展相位转换为所述正弦信号的初相位；

全相位差模块1310，用于将所述截止相位与所述初相位的差值转换为所述正弦信号的全相位差。

本实施方式，将所述正向信号序列反向输出获得反褶序列，将反褶序列进行截短处理，获得反褶截短序列；以所测参考频率的余弦函数和所述参考频率的正弦函数分别与反褶序列和反褶截短序列相乘，生成两组实频向量序列和虚频向量序列；通过对两组虚频向量序列和实频向量序列数字陷波，生成两组虚数向量陷波序列和实数向量陷波序列，进而积分生成两组虚数向量积分值和实数向量积分值；再根据预设的相位转换规则，将两组实数向量积分值和虚数向量积分值转换为两个相位；再将两个相位进行扩展，获得扩展相位。进而根据预设的相位转换规则，将两个扩展相位转换为所述正弦信号的截止相位，将所述截止相位与所述初相位的差值转换为所述正弦信号的全相位差。

由于数字陷波可以在某一个频率点迅速衰减输入信号，以达到阻碍此频率信号通过的效果，因此，当该数字陷波的陷波频率点设为对应混频干扰频率点时，该数字陷波对所述混频干扰频率具有深度的抑制作用。如此，采用本发明的正弦信号的频率检测方法和设备，其对低频正弦信号的频率测量精度高，且抗谐波和噪声干扰性好。

全相位差检测系统中各模块与以上所述的正弦信号的全相位差检测方法中的各步骤一一对应。

其中，初步序列长度模块1010，电力系统频率范围在45Hz-55Hz，优选的，取正弦信号频率下限f_min为45Hz；并根据实际需要设置所述预设信号周期数C_2π，优选的，所述预设信号周期数C_2π取整数。

在一些实施例中，取C_2π为整数12。

初步采样序列长度计算，为上述式(1)。

对于初步采样模块1020，对所述正弦信号进行初步采样，获得所述正弦信号的初步采样序列。对单基波频率正弦信号，获得所述正弦信号的初步采样序列为上述式(2)。

对于频率初测模块1030，可通过零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法对所述初步采样序列进行频率初测，获取所述初步频率。

所述初步频率表达为上述式(3)。

优选地，所述参考频率表达为上述式(4)。

对于周期序列模块1040，根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度：

在一些实施例中，所述正弦信号的单位周期序列长度计算，为上述式(5)。

所述单位周期序列长度整数化存在1个采样间隔内的误差。

对于序列长度模块1050，根据所述预设信号周期数和所述单位周期序列长度，获得预设序列长度：

在一些实施例中，所述预设序列长度为所述单位周期序列长度的12倍，所述预设序列长度计算，为上述式(6)。

对于信号序列模块1060，根据预设序列长度，从所述初步采样序列中，获得信号序列：

在一些实施例中，所述预设序列长度的信号序列，为上述式(7)。

所述信号序列的图形表达，如图3所示。

对于反褶序列模块1070，相对正向信号序列，反褶序列表达为上述式(8)。

所述反褶序列的图形表达，如图3所示。

对于截短模块1080，分别将所述反褶序列和所述正向信号序列进行截短，获得序列长度相同的反褶截短序列和正向截短序列；

在一些实施例中，所述反褶截短序列长度为所述单位周期序列长度的11.25倍，所述反褶截短序列长度计算为上述式(9)。

优选地，所述反褶截短序列表达为上述式(10)。

将所述正向信号序列进行截短，获得正向截短序列，所述正向截短序列表达为上述式(11)：

所述正向截短序列的长度与所述反褶截短序列的长度相同，表达为式(9)。

对于第一混频模块1090，优选地，所述参考频率的余弦函数和所述参考频率的正弦函数可分别为以所述参考频率为频率、以T_n为间隔离散变量的正弦函数和余弦函数。

在一些实施例中，以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述信号序列相乘，得到所述第一实频向量序列和所述第一虚频向量序列为上述式(12)。

对于第二混频模块1100，以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述截短信号序列相乘，得到所述第二实频向量序列和所述第二虚频向量序列为上述式(13)。

同理，对于第三混频模块1110、第四混频模块1120，优选地，所述参考频率的余弦函数和所述参考频率的正弦函数可分别为以所述参考频率为频率、以T_n为间隔离散变量的正弦函数和余弦函数。

在一些实施例中，以所述余弦函数和所述正弦函数分别与所述正向信号序列相乘，生成第三实频向量序列和第三虚频向量序列，为上述式(14)。

在一些实施例中，以所述余弦函数和所述正弦函数分别与所述正向截短序列相乘，生成第四实频向量序列和第四虚频向量序列，为上述式(15)：

对于第一陷波模块1130，所述第一实频向量序列和所述第一虚频向量序列中包含混频干扰频率。当输入信号中包含直流成分、次谐波成分及分次谐波成分时，所述混频干扰频率将更加复杂，这些混频干扰频率严重影响计算准确度。虽然窗口函数和积分运算本身对混频干扰频率具有良好的衰减作用，但没有针对性，不能够对所述复杂的混频干扰频率产生深度的抑制作用，不能满足参数的高准确度计算需要。

为了有针对性的抑制所述混频干扰频率的影响，可以对该混频干扰进行数字陷波器，由于数字陷波可以在某一个频率点迅速衰减输入信号，以达到阻碍此频率信号通过的效果，因此，当该数字陷波的陷波频率点设为对应混频干扰频率点时，该数字陷波对所述混频干扰频率具有深度的抑制作用。

优选地，数字陷波器具体采用滑动三角窗算术平均陷波器，即将若干个连续离散值乘以三角窗函数后再相加，然后取其算术平均值作为本次陷波值输出。滑动三角窗算术平均陷波器需要设置三角窗参数，所述三角窗参数具体指三角窗函数序列的长度N_W。在三角窗参数N_W值为信号周期序列长度的3倍，可以对1/3分次谐波产生的混频干扰频率进行抑制。在三角窗参数N_W值为信号周期序列长度的4倍，可以对直流、1/2分次、1次、2次、3次、4次、5次谐波等产生的混频干扰频率进行抑制。

考虑到实际存在误差等因数，如所述参数存在1个采样间隔以内的误差，为了深度抑制混频干扰频率影响，还可再增加一级滑动矩形窗算术平均陷波器，即将若干个连续离散值直接相加，然后取其算术平均值作为本次陷波值输出。滑动矩形窗算术平均陷波器需要设置矩形窗参数，所述矩形窗参数具体指矩形窗函数序列的长度N_D。矩形窗参数N_D取值为信号周期序列长度的1.5倍，可以对1/3分次谐波产生的混频干扰频率进行抑制。而N_D取值为信号周期序列长度的2倍，可以对直流、1/2分次、1次、2次、3次、4次、5次谐波等产生的混频干扰频率进行抑制。

优选地，数字陷波器由二级滑动三角窗算术平均陷波器加二级滑动矩形窗算术平均陷波器共四级数字陷波器所构成该四级数字陷波器式可为上述式(16)：

在一些实施例中，三角窗参数N_W1取值为所述参考频率的单位周期序列长度的3倍，矩形窗参数N_D1取值为所述参考频率的单位周期序列长度的1.5倍，三角窗参数N_W2取值为所述参考频率的单位周期序列长度的4倍，矩形窗参数N_D1取值为所述参考频率的单位周期序列长度的2倍。四级数字陷波需要使用10.5倍信号周期序列长度。

对上述实施例，在所述正弦信号的基波频率100π、单位rad/s，得到三角窗算术平均陷波算器1的频域特性，如图5所示。得到三角窗算术平均陷波算器2的频域特性，如图6所示。

优选地，在所述混频干扰频率成分得到完全抑制前提下，所述第一实频向量陷波序列和所述第一虚频向量陷波序列为上述式(17)。

对于第一积分模块1140，优选地，可通过本领域技术人员惯用的积分器进行积分运算。积分运算式为上述式(18)。

同理，对于第二陷波模块1150，在一些实施例中，陷波参数N_D1取值为所述参考频率的单位周期序列长度的1.5倍，陷波参数N_D2取值为所述参考频率的单位周期序列长度的2倍，得到第二实频向量陷波序列和第二虚频向量陷波序列，为上述式(19)。

对于第二积分模块1160，优选地，积分运算式可为上述式(20)。

对于第三陷波模块1170，同理，在一些实施例中，陷波参数N_D1取值为所述参考频率的单位周期序列长度的1.5倍，陷波参数N_D2取值为所述参考频率的单位周期序列长度的2倍，得到第三实频向量陷波序列和第三虚频向量陷波序列，为上述式(21)：

对于第三积分模块1180，优选地，可通过本领域技术人员惯用的积分器进行积分运算。积分运算式为上述式(22)。

对于第四陷波模块1190，在一些实施例中，陷波参数N_D1取值为所述参考频率的单位周期序列长度的1.5倍，陷波参数N_D2取值为所述参考频率的单位周期序列长度的2倍，得到第四实频向量陷波序列和第四虚频向量陷波序列，为上述式(23)：

对于第四积分模块1200，优选地，积分运算式可为上述式(24)。

对于第一相位模块1210、第二相位模块1220、第三相位模块1230、第四相位模块1240优选地，预设的相位转换规则对应于虚频向量积分值和实频向量转换为相位的转换式。

优选地，可通过上述式式(25)-(28)将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位，将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位，将所述第三虚频向量积分值与所述第三实频向量积分值转换为第三相位，将所述第四虚频向量积分值与所述第四实频向量积分值转换为第四相位。

在一些实施例中，根据预设的相位转换规则，将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位的步骤包括以下步骤：

获取所述第一虚频向量积分值与所述第一实频向量积分值的比值；

获取所述比值的反正切函数值的相反数，生成所述第一相位。

在一些实施例中，根据所述预设的相位转换规则，将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位的步骤包括以下步骤：

获取所述第二虚频向量积分值与所述第二实频向量积分值的比值；

获取所述比值的反正切函数值的相反数，生成所述第二相位。

类似地，所述第三虚频向量积分值与所述第三实频向量积分值转换为第三相位的步骤、所述第四虚频向量积分值与所述第四实频向量积分值转换为第四相位的步骤亦采用上述预设的相位转换规则。

对于第一相位扩展模块1250，所述第一相位在0～±0.5πrad范围，但实际序列相位可能会超出±0.5πrad范围，因此必须根据相位扩展规则对所述第一相位进行扩展，扩展后的相位范围在0～±πrad范围，第一扩展相位，为上述式(29)。

在一些实施例中，第一相位扩展模块1160可用于：

如果所述第一实频向量积分值大于等于零的同时，所述第一虚频向量积分值的相反数大于等于零，则第一扩展相位等于所述第一相位；

如果所述第一实频向量积分值大于等于零的同时，所述第一虚频向量积分值的相反数小于零，则第一扩展相位等于所述第二相位；

如果所述第一实频向量积分值小于零的同时，所述第一虚频向量积分值的相反数大于等于零，则第一扩展相位等于所述第二相位加πrad；

如果所述第一实频向量积分值小于零的同时，所述第一虚频向量积分值相的反数小于零，则第一扩展相位等于所述第二相位减πrad。

同理，对于第二相位模块1220、第三相位模块1230、第四相位模块1240，所述第二相位、所述第三相位、所述第四相位采用上述所述相位扩展规则进行扩展，扩展后的相位范围在0～±πrad范围，第二扩展相位、第三扩展相位、第四扩展相位分别为上述式(30)-(32)。

对于截止相位模块1290，所述预设的截止相位转换规则可对应于第一相位和第二相位转换为截止相位的公式。根据式(29)和式(30)，可生成与所述预设的截止相位转换规则对应的截止相位为上述式(33)。

在一些实施例中，根据预设的截止相位转换规则，将所述第一相位和所述第二相位转换为所述正弦信号的截止相位的步骤包括以下步骤：

获取所述第一相位与所述反褶截短序列长度的乘积，生成第一乘积；

获取所述第二相位与所述反褶序列长度的乘积，生成第二乘积；

获取所述第一乘积与所述第二乘积的差值，生成第一差值；

获取所述反褶截短序列的长度与所述反褶序列的长度的差值，生成第二差值；

获取所述第一差值与所述第二差值的比值，生成所述截止相位。

对于初相位模块1300，所述预设的初相位转换规则可对应于第三相位和第四相位转换为初相位的公式。根据式(31)和式(32)，可生成与所述预设的初相位转换规则对应的初相位为上述式(34)。

在一些实施例中，根据预设的初相位转换规则，将所述第三相位和所述第四相位转换为所述正弦信号的初相位的步骤包括以下步骤：

获取所述预设序列长度与信号周期序列长度的0.75倍的差值，生成所述截短信号序列的序列长度。

获取所述第三相位与所述正向截短信号序列的长度的乘积，生成第三乘积。

获取所述第四相位与所述正向信号序列的长度的乘积，生成第四乘积。

获取所述第三乘积与所述第四乘积的差值，生成第三差值。

获取所述正向截短信号序列的长度与所述正向信号序列的长度的差值，生成第四差值。

获取所述第三差值与所述第四差值的比值，生成所述初相位。

对于全相位差模块1310，优选地，根据所述正弦信号的截止相位、所述正弦信号的初相位、预设信号周期数，所述正弦信号的全相位差公式为上述式(35)。

请参阅图4，图4是本发明正弦信号的频率检测方法的初相位的相对误差的实验结果示意图。

为了验证本发明正弦信号的全相位差检测系统具有较高的准确度，给出一实验信号，为式(36)：

ω＝2πft

在实验信号基波频率变化范围在45Hz-55Hz，取信号整数周期数为12(由于存在误差，信号序列对应的信号整数周期数大约为12)，信号初相位变化范围在0～±π/2、单位rad，信号的采样频率为20kHz，信号的离散数据量化位数24bit，频率初测相对误差<|±0.25％|。得到实验信号全相位差的相对误差绝对值|△PH_err(f)|随信号基波频率f变化特性的实验结果图，图4所示。图4给出实验信号的全相位差的检测准确度在10^-10量级。

以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上陈述仅以实施例来进一步说明本发明的技术内容，以便于读者更容易理解，但不代表本发明的实施方式仅限于此，任何依本发明所做的技术延伸或再创造，均受本发明的保护。

Claims (10)

1.一种正弦信号的全相位差检测方法，其特征在于，包括以下步骤：

根据正弦信号频率范围的下限和预设采样频率及预设信号周期数，获得初步采样序列长度；

根据所述初步采样序列长度，对所述正弦信号进行初步采样，获得所述正弦信号的初步采样序列；

对所述初步采样序列进行频率初测，生成所述正弦信号的初步频率，以所述初步频率给定参考频率；

根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度；

根据所述预设信号周期数和所述单位周期序列长度，获得预设序列长度；

根据预设序列长度，从所述初步采样序列中，获得正向信号序列；

将所述正向信号序列反向输出，获得所述正向信号序列的反褶序列；

分别将所述反褶序列和所述正向信号序列进行截短，获得序列长度相同的反褶截短序列和正向截短序列；

以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述反褶序列相乘，生成第一实频向量序列和第一虚频向量序列；

以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述反褶截短序列相乘，生成第二实频向量序列和第二虚频向量序列；

以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述正向信号序列相乘，生成第三实频向量序列和第三虚频向量序列；

以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述正向截短序列相乘，生成第四实频向量序列和第四虚频向量序列；

分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波，生成第一实频向量陷波序列和第一虚频向量陷波序列；

分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算，生成第一实频向量积分值和第一虚频向量积分值；

分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波，生成第二实频向量陷波序列和第二虚频向量陷波序列；

分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算，生成第二实频向量积分值和第二虚频向量积分值；

分别对所述第三实频向量序列和所述第三虚频向量序列进行数字陷波，生成第三实频向量陷波序列和第三虚频向量陷波序列；

分别对所述第三实频向量陷波序列和所述第三虚频向量陷波序列进行积分运算，生成第三实频向量积分值和第三虚频向量积分值；

分别对所述第四实频向量序列和所述第四虚频向量序列进行数字陷波，生成第四实频向量陷波序列和第四虚频向量陷波序列；

分别对所述第四实频向量陷波序列和所述第四虚频向量陷波序列进行积分运算，生成第四实频向量积分值和第四虚频向量积分值；

根据预设的相位转换规则，将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位；

根据预设的相位转换规则，将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位；

根据预设的相位转换规则，将所述第三虚频向量积分值与所述第三实频向量积分值转换为第三相位；

根据预设的相位转换规则，将所述第四虚频向量积分值与所述第四实频向量积分值转换为第四相位；

根据预设的相位扩展规则，对所述第一相位进行扩展，获得第一扩展相位；

根据预设的相位扩展规则，对所述第二相位进行扩展，获得第二扩展相位；

根据预设的相位扩展规则，对所述第三相位进行扩展，获得第三扩展相位；

根据预设的相位扩展规则，对所述第四相位进行扩展，获得第四扩展相位；

根据预设的截止相位转换规则，将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的截止相位；

根据预设的初相位转换规则，将所述第三扩展相位和所述第四扩展相位转换为所述正弦信号的初相位；

将所述截止相位与所述初相位的差值转换为所述正弦信号的全相位差。

2.如权利要求1所述的正弦信号的全相位差检测方法，其特征在于，根据正弦信号频率范围的下限和预设采样频率及预设信号周期数，获得初步采样序列长度的步骤包括以下步骤：

通过以下所述公式将所述正弦信号频率的下限和预设采样频率及预设信号周期数转换为所述初步采样序列长度：

$N_{start}=\left(\mathrm{int}\right)\frac{C_{2\mathrm{\&pi;}}{f}_n}{f_{\min }};$

其中，N_start为所述初步序列长度，单位无量纲；(int)表示取整；C_2π为所述预设信号周期数，单位无量纲；f_min为正弦信号频率范围的下限，单位Hz；f_n为所述预设采样频率，单位Hz。

3.根据权利要求1所述的正弦信号的全相位差检测方法，其特征在于，所述初步频率通过对所述初步采样序列采用零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法获得。

4.根据权利要求1所述的正弦信号的全相位差检测方法，其特征在于，所述数字陷波由二级滑动三角窗算术平均陷波器加二级滑动矩形窗算术平均陷波器构成。

5.根据权利要求1所述的正弦信号的全相位差检测方法，其特征在于，根据预设的相位扩展规则，将所述第一相位、第二相位、第三相位或第四相位扩展为所述第一扩展相位、第二扩展相位、第三扩展相位、第四扩展相位的步骤包括以下步骤：

如果所述实频向量积分值大于等于零的同时，所述虚频向量积分值的相反数大于等于零，则所述扩展相位等于所述相位；

如果所述实频向量积分值大于等于零的同时，所述虚频向量积分值的相反数小于零，则所述扩展相位等于所述相位；

如果所述实频向量积分值小于零的同时，所述虚频向量积分值的相反数大于等于零，则所述扩展相位等于所述相位加πrad；

如果所述实频向量积分值小于零的同时，所述虚频向量积分值的相反数小于零，则所述扩展相位等于所述相位减πrad。

6.一种正弦信号的全相位差检测系统，其特征在于，包括：

初步序列长度模块，用于根据正弦信号频率范围的下限和预设采样频率及预设信号周期数，获得初步采样序列长度；

初步采样模块，用于根据所述初步采样序列长度，对所述正弦信号进行初步采样，获得所述正弦信号的初步采样序列；

频率初测模块，用于对所述初步采样序列进行频率初测，生成所述正弦信号的初步频率，以所述初步频率给定参考频率；

周期序列模块，用于根据预设采样频率和所述参考频率计算所述正弦信号的单位周期序列长度；

序列长度模块，根据所述预设信号周期数和所述单位周期序列长度，获得预设序列长度；

正向信号序列模块，用于根据所述预设序列长度，从所述初步采样序列中，获得正向信号序列；

反褶模块，用于将所述正向信号序列反向输出，获得所述正向信号序列的反褶序列；

截短模块，用于分别将所述反褶序列和所述正向信号序列进行截短，获得序列长度相同的反褶截短序列和正向截短序列。

第一混频模块，用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述反褶序列相乘，生成第一实频向量序列和第一虚频向量序列；

第二混频模块，用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述反褶截短序列相乘，生成第二实频向量序列和第二虚频向量序列；

第三混频模块，用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述正向信号序列相乘，生成第三实频向量序列和第三虚频向量序列；

第四混频模块，用于以所述参考频率的余弦函数和所述参考频率的正弦函数分别与所述正向截短序列相乘，生成第四实频向量序列和第四虚频向量序列；

第一陷波模块，用于分别对所述第一实频向量序列和所述第一虚频向量序列进行数字陷波，生成第一实频向量陷波序列和第一虚频向量陷波序列；

第一积分模块，用于分别对所述第一实频向量陷波序列和所述第一虚频向量陷波序列进行积分运算，生成第一实频向量积分值和第一虚频向量积分值；

第二陷波模块，用于分别对所述第二实频向量序列和所述第二虚频向量序列进行数字陷波，生成第二实频向量陷波序列和第二虚频向量陷波序列；

第二积分模块，用于分别对所述第二实频向量陷波序列和所述第二虚频向量陷波序列进行积分运算，生成第二实频向量积分值和第二虚频向量积分值；

第三陷波模块，用于分别对所述第三实频向量序列和所述第三虚频向量序列进行数字陷波，生成第三实频向量陷波序列和第三虚频向量陷波序列；

第三积分模块，用于分别对所述第三实频向量陷波序列和所述第三虚频向量陷波序列进行积分运算，生成第三实频向量积分值和第三虚频向量积分值；

第四陷波模块，分别对所述第四实频向量序列和所述第四虚频向量序列进行数字陷波，生成第四实频向量陷波序列和第四虚频向量陷波序列；

第四积分模块，用于分别对所述第四实频向量陷波序列和所述第四虚频向量陷波序列进行积分运算，生成第四实频向量积分值和第四虚频向量积分值；

第一相位模块，用于根据预设的相位转换规则，将所述第一虚频向量积分值与所述第一实频向量积分值转换为第一相位；

第二相位模块，用于将所述第二虚频向量积分值与所述第二实频向量积分值转换为第二相位；

第三相位模块，用于将所述第三虚频向量积分值与所述第三实频向量积分值转换为第三相位；

第四相位模块，用于将所述第四虚频向量积分值与所述第四实频向量积分值转换为第四相位；

第一相位扩展模块，用于根据预设的相位扩展规则，对所述第一相位进行扩展，获得第一扩展相位；

第二相位扩展模块，用于根据预设的相位扩展规则，对所述第二相位进行扩展，获得第二扩展相位；

第三相位扩展模块，用于根据预设的相位扩展规则，对所述第三相位进行扩展，获得第三扩展相位；

第四相位扩展模块，用于根据预设的相位扩展规则，对所述第四相位进行扩展，获得第四扩展相位；

截止相位模块，用于根据预设的截止相位转换规则，将所述第一扩展相位和所述第二扩展相位转换为所述正弦信号的截止相位；

初相位模块，用于根据预设的初相位转换规则，将所述第三扩展相位和所述第四扩展相位转换为所述正弦信号的初相位；

全相位差模块，用于将所述截止相位与所述初相位的差值转换为所述正弦信号的全相位差。

7.根据权利要求6所述的正弦信号的全相位差检测系统，其特征在于，所述初步采样模块根据正弦信号频率范围的下限和预设采样频率及预设信号周期数，获得初步采样序列长度的步骤包括以下步骤：

通过以下所述公式将所述正弦信号频率的下限和预设采样频率及预设信号周期数转换为所述初步采样序列长度：

$N_{start}=\left(\mathrm{int}\right)\frac{C_{2\mathrm{\&pi;}}{f}_n}{f_{\min }};$

其中，N_start为所述初步序列长度，单位无量纲；(int)表示取整；C_2π为所述预设信号周期数，单位无量纲；f_min为正弦信号频率范围的下限，单位Hz；f_n为所述预设采样频率，单位Hz。

8.根据权利要求6所述的正弦信号的全相位差检测系统，其特征在于，所述初测模块通过对所述初步采样序列采用零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法来获得所述初步频率。

9.根据权利要求6所述的正弦信号的全相位差检测系统，其特征在于，所述第一陷波模块、第二陷波模块、第三陷波模块、第四陷波模块包括二滑动三角窗算术平均陷波器和加二级滑动矩形窗算术平均陷波器。

10.根据权利要求6所述的正弦信号的全相位差检测系统，其特征在于，所述初相位模块根据预设的相位扩展规则，将所述第一相位、第二相位、第三相位或第四相位扩展为所述第一扩展相位、第二扩展相位、第三扩展相位、第四扩展相位的步骤包括以下步骤：

如果所述实频向量积分值大于等于零的同时，所述虚频向量积分值的相反数大于等于零，则所述扩展相位等于所述相位；

如果所述实频向量积分值大于等于零的同时，所述虚频向量积分值的相反数小于零，则所述扩展相位等于所述相位；

如果所述实频向量积分值小于零的同时，所述虚频向量积分值的相反数大于等于零，则所述扩展相位等于所述相位加πrad；

如果所述实频向量积分值小于零的同时，所述虚频向量积分值的相反数小于零，则所述扩展相位等于所述相位减πrad。