107 :将所述正向序列反向输出，获得与所述正向序列对应的反褶序列；
[0048] 步骤S108 :分别将所述参考频率的余弦函数和所述参考频率的正弦函数与所述 正向序列相乘，生成第一实频向量序列和第一虚频向量序列；
[0049] 步骤S109 :分别将所述参考频率的余弦函数和所述参考频率的正弦函数与所述 反褶序列相乘，生成第二实频向量序列和第二虚频向量序列；
[0050] 步骤S110 :分别对所述第一实频向量序列和所述第一虚频向量序列进行积分计 算，生成第一实频向量积分值和第一虚频向量积分值；
[0051] 步骤S111 :分别对所述第二实频向量序列和所述第二虚频向量序列进行积分计 算，生成第二实频向量积分值和第二虚频向量积分值；
[0052] 步骤S112 :根据预设的相位转换规则，将所述第一虚频向量积分值与所述第一实 频向量积分值转换为第一相位，将所述第二虚频向量积分值与所述第二实频向量积分值转 换为第二相位；
[0053] 步骤S113 :根据预设的序列平均初相位转换规则，将所述第一相位和所述第二相 位转换为所述正向序列的平均初相位。
[0054] 从以上描述可知，本发明获得平均初相位，满足获得高准确度正交序列幅值的前 提条件，提高正弦参数计算的准确度和抗谐波、噪声干扰性。
[0055] 其中，对于步骤S101，根据电力信号频率范围的下限、预设采样频率和初设整数信 号周期数，获得所述电力信号的初步序列长度；
[0056] 所述电力信号是一种基波成分为主的正弦信号。正弦信号广指正弦函数信号和余 弦函数信号。
[0057] 在一个实施例中，电力系统频率范围在45Hz-55Hz，取电力信号频率下限f_S 45Hz;根据实际需要设置所述初设整数信号周期数C2IIstart，优选地，取C2IIstart为整数13。
[0058] 在一个实施例中，获得所述初步序列长度为式（1):
[0059]
[0060] 式中，Nstart为所述初步序列长度，单位无量纲；（int)表示取整数；C2IIstot为所述 初设整数信号周期数，单位无量纲；f_为所述信号频率范围的下限，单位Hz;1为所述预 设采样频率，单位Hz。
[0061] 对于步骤S102,根据所述初步序列长度对所述电力信号进行初步采样，获得所述 电力信号的初步序列。在一个实施例中，所述电力信号为单基波频率的余弦函数信号，获得 所述电力信号的初步序列为式（2):
[0062]
[0063] 其中，Xstart(n)为初步序列；A为所述电力信号的幅值，单位v 为所述电力信号 的频率，单位rad/s ;Tn为采样间隔，单位s ;f n为所述预设采样频率，单位Hz ;n为所述电力 信号的序列离散数，单位无量纲4为信号（初步序列）初相位，单位rad ;Nstart为初步序列 长度，单位无量纲。
[0064] 对于步骤S103,可通过零交法、基于滤波的算法、基于小波变换算法、基于神经网 络的算法、基于DFT变换的频率算法或基于相位差的频率算法对所述初步序列进行频率初 测，生成所述电力信号的初步频率。
[0065] 在一个实施例中，生成所述初步频率为式（3):
[0066]?。（3);
[0067] 其中，w。为所述初步频率，单位rad/s ;
[0068] 优选地，所述参考频率等于所述初步频率为式（4):
[0069] ws=w。 （4);
[0070] 其中，《s为参考频率，单位rad/s ; ?。为初步频率，单位rad/s。
[0071] 对于步骤S104,根据所述预设采样频率和所述参考频率，获得所述电力信号的单 位周期序列长度；
[0072] 在一个实施例中，获得所述电力信号的单位周期序列长度为式（5):
[0073]
[0074] 式中，N2II为所述单位周期序列长度，单位无量纲；（int)表示取整数；f#所述预 设采样频率，单位Hz尤为Hz单位的参考频率；《 3为rad/s单位的参考频率。
[0075] 所述单位周期序列长度整数化存在1个采样间隔内的误差。
[0076] 对于步骤S105,根据所述初设整数信号周期数得到所述电力信号的预设整数信号 周期数，并根据所述预设整数信号周期数和所述单位周期序列长度，获得所述电力信号的 预设序列长度；
[0077] 在一个实施例中，将所述初设整数信号周期数减1获得预设整数信号周期数，根 据预设整数信号周期数和所述单位周期序列长度，获得预设序列长度，原则上要求预设序 列长度对应整数信号周期数；
[0078] 在一个实施例中，获得所述预设序列长度为式（6):
[0079] N= (int) (C2iiN2ii)
[0080] (6)；
[0081] C2II= C2II start-l
[0082] 其中，N为所述预设序列长度，单位无量纲；（int)表示取整数；N2II为所述单位周 期序列长度，单位无量纲；C2II为所述预设整数信号周期数，单位无量纲；C2IIstart为所述初 设整数信号周期数，单位无量纲；由于存在误差，所述预设序列长度包含所述单位周期序列 长度整数存在的误差。
[0083] 对于步骤S106,根据预设序列起始点和所述预设序列长度，从所述初步序列中获 得所述电力信号的正向序列，在一个实施例中，预设序列起始点值为所述单位周期序列长 度的〇. 5倍；
[0084] 在一个实施例中，获得所述正向序列为式（7):
[0085]
[0086] 其中，Xi(n)为正向序列；Xstart(n)为初步序列；Nset为无量纲单位的预设序列起始 点，初值为〇. 5N2lI;(int)为取整数；A为信号幅值，单位v 为信号频率，单位rad/s;!\为 采样间隔，单位s ;n为序列离散数，单位无量纲；q>为信号初相位，单位rad;_为正向序列 初相位，单位rad;N所述预设序列长度；Nstart为初步序列长度，单位无量纲。
[0087] 所述正向序列的图形表达，如图2所示。
[0088] 对于步骤S107,将所述正向序列反向输出，获得与所述正向序列对应的反褶序 列；
[0089] 在一个实施例中，获得反褶序列为式（8):
[0090] X ! (-n) = X: (N-n) = Acos (- ? Tnn+ 0 1)
[0091] (8)；
[0092] n = 0, 1，2, 3, ? ? ?，N-l
[0093] 式中，XJ-n)为反褶序列；0 1为反褶序列初相位，单位rad。关系上，反褶序列初 相位是正向序列的截止相位；N为所述预设序列长度。
[0094] 所述反褶序列的图形表达，如图2所示。
[0095] 对于步骤S108,分别将所述参考频率的余弦函数和所述参考频率的正弦函数与所 述正向序列相乘，生成第一实频向量序列和第一虚频向量序列；
[0096] 在一个实施例中，在不考虑所述混频干扰频率成分时，得到混频序列为式（9):
[0097]
[0098] 式中，Rjn)为第一实频向量序列山⑷为第一虚频向量序列；Q为rad/s单位的 信号频率与参考频率的频差；N为所述预设序列长度。
[0099] 对于步骤S109,分别将所述参考频率的余弦函数和所述参考频率的正弦函数与所 述反褶序列相乘，生成第二实频向量序列和第二虚频向量序列；
[0100] 在一个实施例中，在不考虑所述混频干扰频率成分时，得到混频序列为式（10):
[0101]
[0102] 式中，R2(n)为第二实频向量序列；I2(n)为第二虚频向量序列；Q为rad/s单位的 信号频率与参考频率的频差；N所述预设序列长度。
[0103] 对于步骤S110,分别对所述第一实频向量序列和所述第一虚频向量序列进行积分 计算，生成第一实频向量积分值和第一虚频向量积分值；
[0104] 在一个实施例中，生成第一实频向量积分值和第一虚频向量积分值，为式（11):
[0105]
[0106]其中，&为第一实频向量积分值；Ii为第一虚频向量积分