8 :根据所述初次正向序列获得所述电力信号的初次正相位，并根据所述 初次反褶序列获得所述电力信号的初次反相位；
[0060] 步骤S109 :根据所述初次正相位和所述初次反相位获得所述电力信号的初次平 均初相位；
[0061] 步骤S110 :将所述初次平均初相位与±π/4进行比较，获得与所述±π/4比较 的初次相位比较值，并根据所述初次相位比较值和所述预设起始点，获得新起始点；
[0062] 步骤S111 :根据所述新起始点和所述预设序列长度，从所述初步序列中获得所述 电力信号的再次正向序列，并根据所述再次正向序列获得所述电力信号的再次反褶序列；
[0063] 步骤S112 :根据所述再次正向序列获得所述电力信号的再次正相位，并根据所述 再次反褶序列获得所述电力信号的再次反相位；
[0064] 步骤S113 :根据所述再次正相位和所述再次反相位获得所述电力信号的再次平 均初相位；
[0065] 步骤S114 :将所述再次正向序列和所述再次反褶序列相加，并根据相加后的结果 和所述再次平均初相位，获得所述电力信号的余弦函数调制序列；
[0066] 步骤S115 :将所述再次正向序列和所述再次反褶序列相减，并根据相减后的结果 和所述再次平均初相位，获得所述电力信号的正弦函数调制序列；
[0067] 步骤S116:从所述余弦函数调制序列中心点输出，获得所述电力信号的零初相位 基准余弦函数调制序列；从所述正弦函数调制序列中心点输出，获得所述电力信号的零初 相位基准正弦函数调制序列；
[0068] 步骤S117 :将预设微调频率的离散正弦函数与所述零初相位基准余弦函数调制 序列相乘获得所述电力信号的第一乘法序列，将所述预设微调频率的离散余弦函数与所述 零初相位基准正弦函数调制序列相乘获得所述电力信号的第二乘法序列；
[0069] 步骤S118 :将所述第二乘法序列与所述第一乘法序列相减，获得所述电力信号的 降低频率的零初相位基准正弦函数序列。
[0070] 从以上描述可知，本发明获得降低频率的零初相位基准正弦函数序列，提高正弦 参数计算的准确度，实际应用价值高。
[0071] 其中，所述电力信号是一种基波成分为主的正弦信号。正弦信号广指正弦函数信 号和余弦函数信号。
[0072] 在一个实施例中，电力系统频率范围在45Hz-55Hz，取电力信号频率下限f_S 45Hz;根据实际需要设置所述预设整数信号周期数C2II，在一个实施例中，取C2II为13。
[0073] 在一个实施例中，对于步骤S101，获得所述初步序列长度为式（1):
[0074]
( 1 )；
[0075] 式中，Nstart为初步序列长度；（int)表示取整；C2II为预设整数信号周期数；f_为 电力信号频率范围的下限，单位Hz;f为预设采样频率，单位Hz。
[0076] 在一个实施例中，所述电力信号为单基波频率的余弦函数信号，对于步骤S102,获 得所述电力信号的初步序列为式（2):
[0077]
[0078] 其中，Xstart(n)为所述初步序列；A为信号幅值，单位v;为信号频率，T为采样 间隔时间，f为预设采样频率，单位Ηζ，η为序列离散数，Φ为初步序列初相位，Nstot为初步 序列长度。
[0079] 对于步骤S103,可通过零交法、基于滤波的算法、基于小波变换算法、基于神经网 络的算法、基于DFT变换的频率算法或基于相位差的频率算法对所述初步序列进行频率初 测，获取所述初步频率。
[0080] 在一个实施例中，生成所述初步频率为式（3):
[0081] ω。 (3)；
[0082] 其中，ω。为初步频率；
[0083] 优选地，所述参考频率等于所述初步频率为式（4):
[0084] 〇s=co。 （4);
[0085] 其中，cos为参考频率，ω。为初步频率。
[0086] 在一个实施例中，对于步骤S104,获得所述电力信号的单位周期序列长度为式 (5)：
[0087]
[0088] 式中，Ν2π为所述单位周期序列长度；（int)为取整数；f为预设采样频率，单位Hz; ?s为参考频率。
[0089] 所述单位周期序列长度整数化存在1个采样间隔内的误差。
[0090] 在一个实施例中，对于步骤S105,获得所述预设序列长度为式（6):
[0091]
[0092] 其中，N为预设序列长度，预设序列长度为奇数，（int)为取整数，Ν2π为所述单位 周期序列长度，C2II为预设整数信号周期数。
[0093] 在一个实施例中，预设起始点为所述单位周期序列长度的0. 5倍；
[0094] 在一个实施例中，对于步骤S106,获得所述初次正向序列为式（7):
[0095]
[0096] 其中，Xstart(n)为初步序列，X+Start(n)为初次正向序列，Pstart预设起始点，Ν2π为所 述单位周期序列长度，（int)为取整数，Α为信号幅值，单位ν，ωι为信号频率，Τ为采样间 隔时间，η为序列离散数，φ?为初次正向序列初相位，Ν预设序列长度。初次正向序列图形 表达，图2所示。
[0097] 在一个实施例中，对于步骤S107,获得初次反褶序列为式（8):
[0098]
[0099] 其中，Xstifft(-n)为初次反褶序列，X+Start(n)为初次正向序列，Α为信号幅值，单位 V，ωi为信号频率，T为采样间隔时间，η为序列离散数，β1为初次反褶序列初相位，N预设 序列长度。所述初次反褶序列的图形表达，如图2所示。
[0100] 在一个实施例中，初次正相位和初次反相位是基于正交混频和积分计算的结果， 在不考虑正交混频的混频干扰频率时，正交混频表达为式（9)，积分计算表达为式（10):
[0101]
[0102] 其中，R+Start (η)为初次正实频混频序列，I+stot (η)为初次正虚频混频序列， Rstalit(_n)为初次反实频混频序列，Istart(_n)为初次反虚频混频序列，cos(c〇sTn)或 cos(-cosTn)为参考频率的离散余弦函数，sin(cosTn)或sin(-cosTn)为参考频率的离散正 弦函数，Ω为信号频率ω;与参考频率ωs的频差，ω;为信号频率，T为采样间隔时间，η为 序列离散数，φ?为初次正向序列初相位，β1为初次反褶序列初相位，Ν为预设序列长度。
[0103]
[0104] 式中，R+Stot初次正实频积分值、单位无量纲，I+start为初次正虚频积分值、单位无 量纲，Rstart为初次反实频积分值、单位无量纲，Istot为初次反虚频混频积分值、单位无量 纲，Ω为信号频率ωι与参考频率ω3的频差，T为采样间隔时间，η为序列离散数，N为预 设序列长度，φ?为初次正向序列初相位，β1为初次反褶序列初相位。
[0105] 在一个实施例中，对于步骤S108,获得初次正相位和初次反相位的表达式为 (11)：
[0106]
[0107] 式中，PH+Start为初次正相位，PHStarts初次反相位，R+start初次正实频积分值、单位 无量纲，I+Stot为初次正虚频积分值、单位无量纲，Rstot为初次反实频积分值、单位无量纲， Istot为初次反虚频混频积分值、单位无量纲，Ω为信号频率ωi与参考频率ω3的频差，τ为采样间隔时间，Ν为预设序列长度，_为初次正向序列初相位，β1为初次反褶序列初相 位。
[0108] 在一个实施例中，对于步骤S109,获得初次平均初相位的表达式为（12):
[0109]
[0110] 式中，PHstartavg为初次平均初相位，PH+start为初次正相位，PHstar#初次反相位， ψ?为初次正向序列初相位，β1为初次反褶序列初相位。
[0111] 在一个实施例中，对于步骤S110,将所述初次平均初相位与PHstart ^与± 31 /4进 行比较为式（13):
[0112]
[0113] 式中，ΔPH_为初次相位比较值，单位rad，PHstartavg为初次平均初相位。
[0114] 在一个实施例中，获得所述新起始点为式（14):
[0115]
[0116] 式中，Pnew为新起始点、单位无量纲，Pstart为预设起始点，APHecJ%初次相位比较 值，单位rad，N2lI为单位周期序列长度，（int)为取整数。
[0117] 在一个实施例中，对于步骤S111，再次正向序列和再次反褶序列为式（15):
[0118]
[0119] 式中，X+end(n)为再次正向序列，Xend(-n)为再次反褶序列，P_为新起始点、单位 无量纲，为再次正向序列初相位，β2为再次反褶序列初相位，ωι为信号频率，T为采样 间隔时间，η为序列离散数，N为预设序列长度。
[0120] 在一个实施例中，再次正相位和再次反相位是基于正交混频和数字滤波计算的结 果。所述数字滤波由2种滤波参数的6级矩形窗算术平均滤波器所构成。
[0121] 在不考虑正交混频的混频干扰频率时，正交混频表达为式（16)，2种滤波参数的6 级矩形窗算术平均滤波器滤波计算表达为式（17):
[0122]
[0123] 式中，R+end(n)为再次正实频混频序列，I+end(n)为再次正虚频混频序列，Rend(-n) 为再次反实频混频序列，I^ (-η)为再次反虚频混频序列，cos(ωsTn)或cos(-ωsTn)为参 考频率的离散余弦函数，sin(cosTn)或sin(-cosTn)为参考频率的离散正弦函数，Ω为信号 频率ω;与参考频率ωs的频差，ω;为信号频率，Τ为采样间隔时间，η为序列离散数，中 为再次正向序列初相位，β2为再次反褶序列初相位，Ν为预设序列长度。
[0124]
[0125] 式中，R+end为再次正实频数字滤波终值，单位无量纲；I+end为再次正虚频数字滤波 终值，单位无量纲；ROTd为再次反数字滤波终值，单位无量纲；I 为再次反虚频数字滤波 终值，单位无量纲；Ω为信号频率Wi与参考频率ω3的频差；Κ(Ω)为数字滤波在频差Ω的幅值增益，单位无量纲；Τ为采样间隔时间；为再次正向序列初相位；β2为再次反褶 序列初相位；ND1为滤波参数1，即对ND1个连续离散值相加，然后取其算术平均值作为本次 滤波值输出；ND2为滤波参数2,即对ND2个连续离散值相加，然后取其算术平均值作为本次 滤波值输出；ND为数字滤波使用序列长度，数量上为6级矩形