CN105548703B - 从电力信号中获取任意初相位正弦函数序列的方法和系统 - Google Patents

Abstract

本发明涉及一种从电力信号中获取任意初相位正弦函数序列的方法和系统，获得零初相位基准余弦函数调制序列和零初相位基准正弦函数调制序列后，设置初相位，将所述初相位的正弦函数与所述零初相位基准余弦函数调制序列相乘获得第一乘法序列，将所述初相位的余弦函数与所述零初相位基准正弦函数调制序列相乘获得第二乘法序列；将所述第一乘法序列和所述第二乘法序列相加，获得所述初相位的正弦函数序列。本发明可根据实际需要设置正弦函数序列的初相位，避开了输入序列任意和不确定初相位问题的影响，对提高正弦参数计算准确度有积极意义。

Description

从电力信号中获取任意初相位正弦函数序列的方法和系统

技术领域

本发明涉及电力系统技术领域，特别是涉及一种从电力信号中获取任意初相位正弦函数序列的方法和系统。

背景技术

电力系统的频率测量、相位测量和幅值测量等在本质上均为正弦参数的测量。快速傅里叶变换算法(FFT)和离散傅里叶变换算法(DFT)是用于正弦参数计算的基本数学方法，在电力系统中有广泛的应用。但随着正弦参数测量技术的发展，快速傅里叶变换算法和离散傅里叶变换算法存在的问题也越显突出，其难以进一步满足电力系统对正弦参数高准确度计算的要求。

在电力系统正弦参数测量方面，还有一些改进的参数测量方法，如零交法、基于滤波的测量法、基于小波变换法、基于神经网络的测量法和基于DFT变换的测量法等。电网运行额定工频在50Hz(赫兹)附近，属于频率较低的正弦频率。由于实际信号处理技术的局限性和信号构成的复杂性，如信号离散采样产生的数据量化背景噪声影响、信号序列截断引起的频谱泄漏问题客观上难以避免、信号任意和不确定初相位问题的影响、信号中的直流和分次谐波及次谐波问题的影响等，因此这些算法的测量精度低，且抗谐波和噪声干扰性差。

发明内容

基于此，有必要针对上述问题，提供一种从电力信号中获取任意初相位正弦函数序列的方法和系统，能够保证正弦参数计算的高准确度。

为解决上述技术问题，本发明采用如下技术方案：

一种从电力信号中获取任意初相位正弦函数序列的方法，包括步骤：

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

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

对所述初步序列进行频率初测，得到所述电力信号的初步频率，根据所述初步频率得到参考频率；

根据所述预设采样频率和所述参考频率，获得所述电力信号的单位周期序列长度；

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

根据所述预设序列长度和预设起始点，从所述初步序列中获得第一正向序列，根据所述第一正向序列获得第一反褶序列；

根据所述第一正向序列获得第一正相位，根据所述第一反褶序列获得第一反相位；

根据所述第一正相位和所述第一反相位获得第一平均初相位；

根据所述第一平均初相位和预设相位值，获得相位比较值，根据所述相位比较值、所述预设起始点和所述单位周期序列长度，获得新起始点；

根据所述预设序列长度和所述新起始点，从所述初步序列中获得第二正向序列，根据第二正向序列获得第二反褶序列；

根据所述第二正向序列获得第二正相位，根据所述第二反褶序列获得第二反相位；

根据所述第二正相位和所述第二反相位获得第二平均初相位；

将所述第二正向序列与所述第二反褶序列相加，得到和序列，根据所述和序列和所述第二平均初相位的余弦函数值，获得余弦函数调制序列；

将所述第二正向序列与所述第二反褶序列相减，得到差序列，根据所述差序列和所述第二平均初相位的正弦函数值，获得正弦函数调制序列；

从所述余弦函数调制序列中心点输出，获得零初相位基准余弦函数调制序列，从所述正弦函数调制序列中心点输出，获得零初相位基准正弦函数调制序列；

设置初相位，将所述初相位的正弦函数与所述零初相位基准余弦函数调制序列相乘获得第一乘法序列，将所述初相位的余弦函数与所述零初相位基准正弦函数调制序列相乘获得第二乘法序列；

将所述第一乘法序列和所述第二乘法序列相加，获得所述初相位的正弦函数序列。

一种从电力信号中获取任意初相位正弦函数序列的系统，包括：

初步序列长度确定模块，用于根据电力信号频率范围的下限、预设采样频率和预设整数信号周期数，获得初步序列长度；

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

参考频率确定模块，用于对所述初步序列进行频率初测，得到所述电力信号的初步频率，根据所述初步频率得到参考频率；

单位周期序列长度确定模块，用于根据所述预设采样频率和所述参考频率，获得所述电力信号的单位周期序列长度；

预设序列长度确定模块，用于根据所述预设整数信号周期数和所述单位周期序列长度，获得预设序列长度，其中所述预设序列长度为奇数；

第一序列获取模块，用于根据所述预设序列长度和预设起始点，从所述初步序列中获得第一正向序列，根据所述第一正向序列获得第一反褶序列；

第一正反相位确定模块，用于根据所述第一正向序列获得第一正相位，根据所述第一反褶序列获得第一反相位；

第一平均初相位确定模块，用于根据所述第一正相位和所述第一反相位获得第一平均初相位；

新起始点确定模块，用于根据所述第一平均初相位和预设相位值，获得相位比较值，根据所述相位比较值、所述预设起始点和所述单位周期序列长度，获得新起始点；

第二序列获取模块，用于根据所述预设序列长度和所述新起始点，从所述初步序列中获得第二正向序列，根据第二正向序列获得第二反褶序列；

第二正反相位确定模块，用于根据所述第二正向序列获得第二正相位，根据所述第二反褶序列获得第二反相位；

第二平均初相位确定模块，用于根据所述第二正相位和所述第二反相位获得第二平均初相位；

余弦函数调制序列确定模块，用于将所述第二正向序列与所述第二反褶序列相加，得到和序列，根据所述和序列和所述第二平均初相位的余弦函数值，获得余弦函数调制序列；

正弦函数调制序列确定模块，用于将所述第二正向序列与所述第二反褶序列相减，得到差序列，根据所述差序列和所述第二平均初相位的正弦函数值，获得正弦函数调制序列；

零初相位调制序列获取模块，用于从所述余弦函数调制序列中心点输出，获得零初相位基准余弦函数调制序列，从所述正弦函数调制序列中心点输出，获得零初相位基准正弦函数调制序列；

乘法序列获取模块，用于设置初相位，将所述初相位的正弦函数与所述零初相位基准余弦函数调制序列相乘获得第一乘法序列，将所述初相位的余弦函数与所述零初相位基准正弦函数调制序列相乘获得第二乘法序列；

初相位正弦函数序列确定模块，用于将所述第一乘法序列和所述第二乘法序列相加，获得所述初相位的正弦函数序列。

本发明从电力信号中获取任意初相位正弦函数序列的方法和系统，可根据实际需要设置正弦函数序列的初相位，避开了输入序列任意和不确定初相位问题的影响，对提高正弦参数计算准确度有积极意义。

附图说明

图1为本发明从电力信号中获取任意初相位正弦函数序列的方法实施例的流程示意图；

图2为本发明初步序列、第一正向序列和第一反褶序列的示意图；

图3为本发明零初相位基准点的示意图；

图4为本发明从电力信号中获取任意初相位正弦函数序列的系统实施例的结构示意图。

具体实施方式

为更进一步阐述本发明所采取的技术手段及取得的效果，下面结合附图及较佳实施例，对本发明的技术方案，进行清楚和完整的描述。

如图1所示，一种从电力信号中获取任意初相位正弦函数序列的方法，包括步骤：

S101、根据电力信号频率范围的下限、预设采样频率和预设整数信号周期数，获得初步序列长度；

S102、根据所述初步序列长度对电力信号进行采样，获得所述电力信号的初步序列；

S103、对所述初步序列进行频率初测，得到所述电力信号的初步频率，根据所述初步频率得到参考频率；

S104、根据所述预设采样频率和所述参考频率，获得所述电力信号的单位周期序列长度；

S105、根据所述预设整数信号周期数和所述单位周期序列长度，获得预设序列长度，其中所述预设序列长度为奇数；

S106、根据所述预设序列长度和预设起始点，从所述初步序列中获得第一正向序列，根据所述第一正向序列获得第一反褶序列；

S107、根据所述第一正向序列获得第一正相位，根据所述第一反褶序列获得第一反相位；

S108、根据所述第一正相位和所述第一反相位获得第一平均初相位；

S109、根据所述第一平均初相位和预设相位值，获得相位比较值，根据所述相位比较值、所述预设起始点和所述单位周期序列长度，获得新起始点；

S110、根据所述预设序列长度和所述新起始点，从所述初步序列中获得第二正向序列，根据第二正向序列获得第二反褶序列；

S111、根据所述第二正向序列获得第二正相位，根据所述第二反褶序列获得第二反相位；

S112、根据所述第二正相位和所述第二反相位获得第二平均初相位；

S113、将所述第二正向序列与所述第二反褶序列相加，得到和序列，根据所述和序列和所述第二平均初相位的余弦函数值，获得余弦函数调制序列；

S114、将所述第二正向序列与所述第二反褶序列相减，得到差序列，根据所述差序列和所述第二平均初相位的正弦函数值，获得正弦函数调制序列；

S115、从所述余弦函数调制序列中心点输出，获得零初相位基准余弦函数调制序列，从所述正弦函数调制序列中心点输出，获得零初相位基准正弦函数调制序列；

S116、设置初相位，将所述初相位的正弦函数与所述零初相位基准余弦函数调制序列相乘获得第一乘法序列，将所述初相位的余弦函数与所述零初相位基准正弦函数调制序列相乘获得第二乘法序列；

S117、将所述第一乘法序列和所述第二乘法序列相加，获得所述初相位的正弦函数序列。

实际电力信号是一种以基波成分为主的正弦信号，如无特殊说明，电力信号均指基波信号，电力信号频率均指基波频率。

对于步骤S101，电力系统频率范围一般在45Hz(赫兹)-55Hz，所以电力信号频率范围的下限f_min可以取为45Hz。预设整数信号周期数C_2π可以根据实际需要设置，例如，将C_2π取为13。初步序列长度可以根据式(1)计算：

其中，N_start为初步序列长度；(int)表示取整；C_2π为预设整数信号周期数；f_min为电力信号频率范围的下限，单位Hz；f为预设采样频率，单位Hz。

对于步骤S102，可以用单基波频率的余弦函数信号表达电力信号，那么初步序列为式(2)：

其中，X_start(n)为初步序列；A为信号幅值，单位v；ω_i为信号频率，T为采样间隔时间，f为预设采样频率，单位Hz，n为序列离散数，为初步序列的初相位，N_start为初步序列长度。

对于步骤S103，可通过零交法、基于滤波的算法、基于小波变换算法、基于神经网络的算法、基于DFT变换的频率算法或基于相位差的频率算法对初步序列进行频率初测，获取初步频率ω_o。在一个实施例中，所述参考频率ω_s＝ω_o。

对于步骤S104，在一个实施例中，所述电力信号的单位周期序列长度计算，为式(3)：

其中，N_2π为单位周期序列长度；(int)为取整数；f为预设采样频率，单位Hz；ω_s为参考频率。所述单位周期序列长度整数化存在1个采样间隔内的误差。

对于步骤S105，预设序列长度为奇数，在一个实施例中，所述预设序列长度计算，为式(4)：

其中，N为预设序列长度，(int)为取整数，N_2π为所述单位周期序列长度，C_2π为预设整数信号周期数。

对于步骤S106，在一个实施例中，预设起始点可以为所述单位周期序列长度的0.5倍，所述第一正向序列，为式(5)：

其中，X_start(n)为初步序列，X_+start(n)为第一正向序列，P_start为预设起始点，N_2π为所述单位周期序列长度，(int)为取整数，A为信号幅值，单位v，ω_i为信号频率，T为采样间隔时间，n为序列离散数，为第一正向序列初相位，N预设序列长度。

所述第一反褶序列，为式(6)：

其中，X_-start(-n)为第一反褶序列，X_+start(n)为第一正向序列，A为信号幅值，单位v，ω_i为信号频率，T为采样间隔时间，n为序列离散数，β1为第一反褶序列初相位，N为预设序列长度。如图2所示，为初步序列、第一正向序列和第一反褶序列的示意图。

对于步骤S107，在一个实施例中，根据对第一正向序列进行正交混频和积分计算的结果，获得第一正相位；根据对第一反褶序列进行正交混频和积分计算的结果，获得第一反相位。在不考虑正交混频的混频干扰频率时，正交混频表达为式(7)，积分计算表达为式(8)：

其中，R_+start(n)为第一正实频混频序列，I_+start(n)为第一正虚频混频序列，R_-start(-n)为第一反实频混频序列，I_-start(-n)为第一反虚频混频序列，cos(ω_sTn)或cos(-ω_sTn)为参考频率的离散余弦函数，sin(ω_sTn)或sin(-ω_sTn)为参考频率的离散正弦函数，Ω为信号频率ω_i与参考频率ω_s的频差，T为采样间隔时间，n为序列离散数，为第一正向序列初相位，β1为第一正向序列初相位，N为预设序列长度。

其中，R_+start第一正实频积分值，单位无量纲，I_+start为第一正虚频积分值，单位无量纲，R_-start为第一反实频积分值，单位无量纲，I_-start为第一反虚频混频积分值，单位无量纲，Ω为信号频率ω_i与参考频率ω_s的频差，T为采样间隔时间，n为序列离散数，N为预设序列长度，为第一正向序列初相位，β1为第一反褶序列初相位，N为预设序列长度。

在一个实施例中，第一正相位和第一反相位的计算，表达为式(9)：

其中，PH_+start为第一正相位，PH_-start为第一反相位，R_+start为第一正实频积分值，单位无量纲，I_+start为第一正虚频积分值，单位无量纲，R_-start为第一反实频积分值，单位无量纲，I_-start为第一反虚频混频积分值，单位无量纲，Ω为信号频率ω_i与参考频率ω_s的频差，T为采样间隔时间，N为预设序列长度，为第一正向序列初相位，β1为第一反褶序列初相位。

对于步骤S108，在一个实施例中，第一平均初相位计算方法，表达为式(10)：

其中，PH_start-avg为第一平均初相位，PH_+start为第一正相位，PH_-start为第一反相位，为第一正向序列初相位，β1为第一反褶序列初相位。

对于步骤S109，在一个实施例中，所述预设相位值可以为±π/4；根据所述第一平均初相位和预设相位值，获得相位比较值的步骤可以包括：

若所述第一平均初相位大于等于0小于等于π/2，根据π/4减去所述第一平均初相位获得相位比较值；

若所述第一平均初相位大于等于-π/2小于等于0，根据-π/4减去所述第一平均初相位获得相位比较值。

具体为式(11)：

其中，△PH_com为相位比较值，单位rad，PH_start-avg为第一平均初相位。

在一个实施例中，所述新起始点计算，为式(12)：

其中，P_new为新起始点，单位无量纲，P_start为预设起始点，△PH_com为相位比较值，单位rad，N_2π为单位周期序列长度，(int)为取整数。

对于步骤S110，第二正向序列和第二反褶序列为式(13)：

其中，X_+end(n)为第二正向序列，X_-end(-n)为第二反褶序列，P_new为新起始点、单位无量纲，为第二正向序列初相位，β2为第二反褶序列初相位，ω_i为信号频率，T为采样间隔时间，n为序列离散数，N为预设序列长度。

对于步骤S111，在一个实施例中，根据对第二正向序列进行正交混频和数字滤波的结果，获得第二正相位；根据对第二反褶序列进行正交混频和数字滤波的结果，获得第二反相位。即第二正相位和第二反相位的计算方法，是基于正交混频和数字滤波计算的结果。在一个实施例中，可以通过2种滤波参数的6级矩形窗算术平均滤波器对正交混频后的第二正向序列和第二反褶序列进行数字滤波。

在不考虑正交混频的混频干扰频率时，正交混频表达为式(14)，2种滤波参数的6级矩形窗算术平均滤波器滤波计算表达为式(15)：

其中，R_+end(n)为第二正实频混频序列，I_+end(n)为第二正虚频混频序列，R_-end(-n)为第二反实频混频序列，I_-end(-n)为第二反虚频混频序列，cos(ω_sTn)或cos(-ω_sTn)为参考频率的离散余弦函数，sin(ω_sTn)或sin(-ω_sTn)为参考频率的离散正弦函数，Ω为信号频率ω_i与参考频率ω_s的频差，T为采样间隔时间，n为序列离散数，为第一正向序列初相位，β1为第一反褶序列初相位，N为预设序列长度。

其中，R_+end为第二正实频数字滤波终值，单位无量纲；I_+end为第二正虚频数字滤波终值，单位无量纲；R_-end为第二反数字滤波终值，单位无量纲；I_-end为第二反虚频数字滤波终值，单位无量纲；Ω为信号频率ω_i与参考频率ω_s的频差；K(Ω)为数字滤波在频差Ω的幅值增益，单位无量纲；T为采样间隔时间；为第二正向序列初相位；β2为第二反褶序列初相位；N_D1为滤波参数1，即对N_D1个连续离散值相加，然后取其算术平均值作为本次滤波值输出；N_D2为滤波参数2，即对N_D2个连续离散值相加，然后取其算术平均值作为本次滤波值输出；N_D为数字滤波使用序列长度，数量上为6级矩形窗算术平均滤波器滤波参数的总和，小于等于预设序列长度N。

在一个实施例中，滤波参数N_D1取值为所述参考频率的单位周期序列长度的1.5倍，目的对1/3分次谐波产生的混频干扰频率进行深度抑制；滤波参数N_D2取值为所述参考频率的单位周期序列长度的2倍，目的对直流、1/2分次、次谐波等产生的混频干扰频率进行深度抑制。2种滤波参数的6级矩形窗算术平均滤波器滤波计算需要使用信号周期序列长度的10.5倍。

滤波参数N_D1和滤波参数N_D2计算为式(16)：

其中，N_D1为数字滤波参数1，单位无量纲，(int)为取整数，N_D2为数字滤波参数2，单位无量纲，N_2π为单位周期序列长度。

在一个实施例中，第二正相位和第二反相位的计算方法，表达为式(17)：

其中，PH_+end为第二正相位，PH_-end为第二反相位，R_+end为第二正实频积分值，单位无量纲，I_+end为第二正虚频积分值，单位无量纲，R_-end为第二反实频积分值，单位无量纲，I_-end为第二反虚频混频积分值，单位无量纲，Ω为信号频率ω_i与参考频率ω_s的频差，T为采样间隔时间，N_D为数字滤波使用序列长度，为第二正向序列初相位，β2为第二反褶序列初相位。

对于步骤S112，第二平均初相位计算方法，表达为式(18)：

其中，PH_end-avg为第二平均初相位，PH_+end为第二正相位，PH_-end为第二反相位，为第二正向序列初相位，β2为第二正向序列初相位。

对于步骤S113，余弦函数调制序列表达为式(19)：

其中，X_cos(n)为余弦函数调制序列；X_+end(n)为第二正向序列；X_-end(-n)为第二反褶序列；PH_end-avg为第二平均初相位；A为余弦函数调制序列幅值，单位v；为余弦函数调制序列初相位，ω_i为信号频率，T为采样间隔时间，n为序列离散数，N为预设序列长度，为第二正向序列初相位，β2为第二反褶序列初相位。

对于步骤S114，正弦函数调制序列表达为式(20)：

其中，X_sin(n)为正弦函数调制序列，X_+end(n)为第二正向序列，X_-end(-n)为第二反褶序列，PH_end-avg为第二平均初相位，A为正弦函数调制序列幅值，单位v，为余弦函数调制序列初相位，ω_i信号频率，T为采样间隔时间，n为序列离散数，N为预设序列长度，为第二正向序列初相位，β2为第二反褶序列初相位。

对于步骤S115，在一个实施例中，零初相位基准余弦函数调制序列表达为式(21)：

其中，X0_cos(n)为零初相位基准余弦函数调制序列，A为余弦函数调制序列幅值，单位v，ω_i信号频率，T为采样间隔时间，n为序列离散数，N为预设序列长度。

在一个实施例中，零初相位基准正弦函数调制序列表达为式(22)：

其中，X0_sin(n)为零初相位基准正弦函数调制序列，A为正弦函数调制序列幅值，单位v，ω_i信号频率，T为采样间隔时间，n为序列离散数，N为预设序列长度。零初相位基准点图形如图3所示。

对于步骤S116，初相位的范围在0～±π/2。在一个实施例中，得到的第一乘法序列和第二乘法序列的表达式为(23)：

其中，X3(n)为第一乘法序列，X4(n)为第二乘法序列，为所述初相位，为所述初相位的正弦函数，X0_cos(n)为所述零初相位基准余弦函数调制序列，为所述初相位的余弦函数，X0_sin(n)为所述零初相位基准正弦函数调制序列，n为序列离散数，N为预设序列长度。

对于步骤S117，所述初相位的正弦函数序列表达为式(24)：

其中，为所述初相位的正弦函数序列。

基于同一发明构思，本发明还提供一种从电力信号中获取任意初相位正弦函数序列的系统，下面结合附图对本发明系统的具体实施方式做详细描述。

如图4所示，一种从电力信号中获取任意初相位正弦函数序列的系统，包括：

初步序列长度确定模块101，用于根据电力信号频率范围的下限、预设采样频率和预设整数信号周期数，获得初步序列长度；

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

参考频率确定模块103，用于对所述初步序列进行频率初测，得到所述电力信号的初步频率，根据所述初步频率得到参考频率；

单位周期序列长度确定模块104，用于根据所述预设采样频率和所述参考频率，获得所述电力信号的单位周期序列长度；

预设序列长度确定模块105，用于根据所述预设整数信号周期数和所述单位周期序列长度，获得预设序列长度，其中所述预设序列长度为奇数；

第一序列获取模块106，用于根据所述预设序列长度和预设起始点，从所述初步序列中获得第一正向序列，根据所述第一正向序列获得第一反褶序列；

第一正反相位确定模块107，用于根据所述第一正向序列获得第一正相位，根据所述第一反褶序列获得第一反相位；

第一平均初相位确定模块108，用于根据所述第一正相位和所述第一反相位获得第一平均初相位；

新起始点确定模块109，用于根据所述第一平均初相位和预设相位值，获得相位比较值，根据所述相位比较值、所述预设起始点和所述单位周期序列长度，获得新起始点；

第二序列获取模块110，用于根据所述预设序列长度和所述新起始点，从所述初步序列中获得第二正向序列，根据第二正向序列获得第二反褶序列；

第二正反相位确定模块111，用于根据所述第二正向序列获得第二正相位，根据所述第二反褶序列获得第二反相位；

第二平均初相位确定模块112，用于根据所述第二正相位和所述第二反相位获得第二平均初相位；

余弦函数调制序列确定模块113，用于将所述第二正向序列与所述第二反褶序列相加，得到和序列，根据所述和序列和所述第二平均初相位的余弦函数值，获得余弦函数调制序列；

正弦函数调制序列确定模块114，用于将所述第二正向序列与所述第二反褶序列相减，得到差序列，根据所述差序列和所述第二平均初相位的正弦函数值，获得正弦函数调制序列；

零初相位调制序列获取模块115，用于从所述余弦函数调制序列中心点输出，获得零初相位基准余弦函数调制序列，从所述正弦函数调制序列中心点输出，获得零初相位基准正弦函数调制序列；

乘法序列获取模块116，用于设置初相位，将所述初相位的正弦函数与所述零初相位基准余弦函数调制序列相乘获得第一乘法序列，将所述初相位的余弦函数与所述零初相位基准正弦函数调制序列相乘获得第二乘法序列；

初相位正弦函数序列确定模块117，用于将所述第一乘法序列和所述第二乘法序列相加，获得所述初相位的正弦函数序列。

在一个实施例中，所述余弦函数调制序列确定模块113可以根据表达式获得余弦函数调制序列X_cos(n)，其中，X_+end(n)为第二正向序列，X_-end(-n)为第二反褶序列，PH_end-avg为第二平均初相位。

在一个实施例中，所述正弦函数调制序列确定模块114可以根据表达式获得正弦函数调制序列X_sin(n)，其中，X_+end(n)为第二正向序列，X_-end(-n)为第二反褶序列，PH_end-avg为第二平均初相位。

在一个实施例中，所述零初相位调制序列获取模块115可以根据表达式获得所述零初相位基准余弦函数调制序列X0_cos(n)，其中，n为序列离散数，N为预设序列长度，X_cos(n)为余弦函数调制序列。

在一个实施例中，所述零初相位调制序列获取模块115可以根据表达式获得所述零初相位基准正弦函数调制序列X0_sin(n)，其中，n为序列离散数，N为预设序列长度，X_sin(n)为正弦函数调制序列。

在一个实施例中，所述乘法序列获取模块116可以根据表达式获得第一乘法序列X3(n)，根据表达式获得第二乘法序列X4(n)；

其中，为所述初相位，为所述初相位的正弦函数，X0_cos(n)为所述零初相位基准余弦函数调制序列，为所述初相位的余弦函数，X0_sin(n)为所述零初相位基准正弦函数调制序列，n为序列离散数，N为预设序列长度。

本发明系统其它技术特征与本发明方法相同，在此不予赘述。

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

以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (8)

1.一种从电力信号中获取任意初相位正弦函数序列的方法，其特征在于，包括步骤：

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

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

对所述初步序列进行频率初测，得到所述电力信号的初步频率，根据所述初步频率得到参考频率；

根据所述预设采样频率和所述参考频率，获得所述电力信号的单位周期序列长度；

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

根据所述预设序列长度和预设起始点，从所述初步序列中获得第一正向序列，根据所述第一正向序列获得第一反褶序列；

根据所述第一正向序列获得第一正相位，根据所述第一反褶序列获得第一反相位；

根据所述第一正相位和所述第一反相位获得第一平均初相位；

根据所述第一平均初相位和预设相位值，获得相位比较值，根据所述相位比较值、所述预设起始点和所述单位周期序列长度，获得新起始点；

根据所述预设序列长度和所述新起始点，从所述初步序列中获得第二正向序列，根据第二正向序列获得第二反褶序列；

根据所述第二正向序列获得第二正相位，根据所述第二反褶序列获得第二反相位；

根据所述第二正相位和所述第二反相位获得第二平均初相位；

将所述第二正向序列与所述第二反褶序列相加，得到和序列，根据所述和序列和所述第二平均初相位的余弦函数值，获得余弦函数调制序列；其中，根据表达式获得余弦函数调制序列X_cos(n)，X_+end(n)为第二正向序列，X_-end(-n)为第二反褶序列，PH_end-avg为第二平均初相位；

将所述第二正向序列与所述第二反褶序列相减，得到差序列，根据所述差序列和所述第二平均初相位的正弦函数值，获得正弦函数调制序列；其中，根据表达式获得正弦函数调制序列X_sin(n)，X_+end(n)为第二正向序列，X_-end(-n)为第二反褶序列，PH_end-avg为第二平均初相位；

从所述余弦函数调制序列中心点输出，获得零初相位基准余弦函数调制序列，从所述正弦函数调制序列中心点输出，获得零初相位基准正弦函数调制序列；

设置初相位，将所述初相位的正弦函数与所述零初相位基准余弦函数调制序列相乘获得第一乘法序列，将所述初相位的余弦函数与所述零初相位基准正弦函数调制序列相乘获得第二乘法序列；

将所述第一乘法序列和所述第二乘法序列相加，获得所述初相位的正弦函数序列。

2.根据权利要求1所述的从电力信号中获取任意初相位正弦函数序列的方法，其特征在于，根据表达式获得第一乘法序列X3(n)，根据表达式获得第二乘法序列X4(n)；

其中，为所述初相位，为所述初相位的正弦函数，X0_cos(n)为所述零初相位基准余弦函数调制序列，为所述初相位的余弦函数，X0_sin(n)为所述零初相位基准正弦函数调制序列，n为序列离散数，N为预设序列长度。

3.根据权利要求2所述的从电力信号中获取任意初相位正弦函数序列的方法，其特征在于，根据表达式获得所述零初相位基准余弦函数调制序列X0_cos(n)，其中，n为序列离散数，N为预设序列长度，X_cos(n)为余弦函数调制序列。

4.根据权利要求2所述的从电力信号中获取任意初相位正弦函数序列的方法，其特征在于，根据表达式获得所述零初相位基准正弦函数调制序列X0_sin(n)，其中，n为序列离散数，N为预设序列长度，X_sin(n)为正弦函数调制序列。

5.一种从电力信号中获取任意初相位正弦函数序列的系统，其特征在于，包括：

初步序列长度确定模块，用于根据电力信号频率范围的下限、预设采样频率和预设整数信号周期数，获得初步序列长度；

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

参考频率确定模块，用于对所述初步序列进行频率初测，得到所述电力信号的初步频率，根据所述初步频率得到参考频率；

单位周期序列长度确定模块，用于根据所述预设采样频率和所述参考频率，获得所述电力信号的单位周期序列长度；

预设序列长度确定模块，用于根据所述预设整数信号周期数和所述单位周期序列长度，获得预设序列长度，其中所述预设序列长度为奇数；

第一序列获取模块，用于根据所述预设序列长度和预设起始点，从所述初步序列中获得第一正向序列，根据所述第一正向序列获得第一反褶序列；

第一正反相位确定模块，用于根据所述第一正向序列获得第一正相位，根据所述第一反褶序列获得第一反相位；

第一平均初相位确定模块，用于根据所述第一正相位和所述第一反相位获得第一平均初相位；

新起始点确定模块，用于根据所述第一平均初相位和预设相位值，获得相位比较值，根据所述相位比较值、所述预设起始点和所述单位周期序列长度，获得新起始点；

第二序列获取模块，用于根据所述预设序列长度和所述新起始点，从所述初步序列中获得第二正向序列，根据第二正向序列获得第二反褶序列；

第二正反相位确定模块，用于根据所述第二正向序列获得第二正相位，根据所述第二反褶序列获得第二反相位；

第二平均初相位确定模块，用于根据所述第二正相位和所述第二反相位获得第二平均初相位；

余弦函数调制序列确定模块，用于将所述第二正向序列与所述第二反褶序列相加，得到和序列，根据所述和序列和所述第二平均初相位的余弦函数值，获得余弦函数调制序列；

其中，所述余弦函数调制序列确定模块根据表达式获得余弦函数调制序列X_cos(n)，X_+end(n)为第二正向序列，X_-end(-n)为第二反褶序列，PH_end-avg为第二平均初相位；

正弦函数调制序列确定模块，用于将所述第二正向序列与所述第二反褶序列相减，得到差序列，根据所述差序列和所述第二平均初相位的正弦函数值，获得正弦函数调制序列；

其中，所述正弦函数调制序列确定模块根据表达式获得正弦函数调制序列X_sin(n)，X_+end(n)为第二正向序列，X_-end(-n)为第二反褶序列，PH_end-avg为第二平均初相位；

零初相位调制序列获取模块，用于从所述余弦函数调制序列中心点输出，获得零初相位基准余弦函数调制序列，从所述正弦函数调制序列中心点输出，获得零初相位基准正弦函数调制序列；

乘法序列获取模块，用于设置初相位，将所述初相位的正弦函数与所述零初相位基准余弦函数调制序列相乘获得第一乘法序列，将所述初相位的余弦函数与所述零初相位基准正弦函数调制序列相乘获得第二乘法序列；

初相位正弦函数序列确定模块，用于将所述第一乘法序列和所述第二乘法序列相加，获得所述初相位的正弦函数序列。

6.根据权利要求5所述的从电力信号中获取任意初相位正弦函数序列的系统，其特征在于，所述乘法序列获取模块根据表达式获得第一乘法序列X3(n)，根据表达式获得第二乘法序列X4(n)；

其中，为所述初相位，为所述初相位的正弦函数，X0_cos(n)为所述零初相位基准余弦函数调制序列，为所述初相位的余弦函数，X0_sin(n)为所述零初相位基准正弦函数调制序列，n为序列离散数，N为预设序列长度。

7.根据权利要求6所述的从电力信号中获取任意初相位正弦函数序列的系统，其特征在于，所述零初相位调制序列获取模块根据表达式获得所述零初相位基准余弦函数调制序列X0_cos(n)，其中，n为序列离散数，N为预设序列长度，X_cos(n)为余弦函数调制序列。

8.根据权利要求6所述的从电力信号中获取任意初相位正弦函数序列的系统，其特征在于，所述零初相位调制序列获取模块根据表达式获得所述零初相位基准正弦函数调制序列X0_sin(n)，其中，n为序列离散数，N为预设序列长度，X_sin(n)为正弦函数调制序列。