一种电力系统的谐波谐振模态的分析方法与流程

1.本发明涉及谐波的谐振模态的分析技术领域，特别是涉及一种电力系统的谐波谐振模态的分析方法。


背景技术：

2.谐波谐振是导致电压、电流畸变的主要原因之一，严重时可能造成设备绝缘损坏等十分严重的后果。如今，随着交直流电网互联、大规模可再生能源接入、电力电子设备广泛应用，如此大量非线性设备和负荷投运造成电力系统谐波污染及谐振问题日趋严重。随着电力系统快速发展，大量可再生能源接入和交直流电网互联导致电力系统规模越来越大，节点数量越来越多，结构越来越复杂，对电力系统的谐波谐振进行快速分析有助于快速发现谐波谐振问题，满足未来电力系统的快速性和复杂性。
3.目前，谐波谐振分析方法主要有谐波传递放大法、频率扫描法、模态分析法、s域模态分析法。其中，模态分析法由于能够提供有助于谐振治理的母线参与因子、元件灵敏度等有效信息而被越来越多地应用。传统模态分析需要以固定频率间隔(固定步长)在逐个频率点对系统节点导纳矩阵求逆并进行特征值分解，然而由于矩阵求逆及矩阵特征值分解时间和空间复杂度高，，导致传统模态分析方法耗时较长，且随着系统节点规模增加，该问题愈加严重。虽然增加固定步长可以减少矩阵求逆和特征值分解计算的次数，从而提高模态分析速度，但是却以牺牲谐振频率精确度为代价，得到的谐振频率误差较大。因此，随着电力系统规模越来越大，如何在保证谐振频率计算精度前提下实现电力系统谐波谐振模态的快速分析是亟需解决的问题。


技术实现要素：

4.本发明的目的在于，提出一种电力系统的谐波谐振模态的分析方法，解决如何在保证谐振频率计算精度前提下，实现电力系统谐波谐振模态的快速分析的技术问题。
5.一方面，提供一种电力系统的谐波谐振模态的分析方法，包括：
6.获取所述电力系统的数据信息，并对所述数据信息中的各参数进行初始化，其中，初始化后的各参数中至少包括扫描点的初始频率、扫描点的终止频率、扫描点的初始步长、用于处理扫描点的存储特征值矩阵、进行扫描的连续两次扫描的第一次步长和第二次步长、存储扫描点的存储扫描点矩阵；
7.根据所述初始步长和所述初始频率确定第一扫描点，根据所述第一扫描点和所述第一次步长确定第二扫描点，将所述第二扫描点存入所述存储扫描点矩阵；确定对应的第一扫描点特征值和第二扫描点特征值及其两者模态序列号之间的差值，并将所述第一扫描点特征值和所述第二扫描点特征值按照所述第一扫描点和所述第二扫描点的顺序存入所述存储特征值矩阵；
8.根据第二扫描点和所述第二次步长确定第三扫描点，将所述第三扫描点存入所述存储扫描点矩阵，判断所述第三扫描点处的频率是否超过所述终止频率；
9.当所述第三扫描点处的频率未超过所述终止频率时，确定对应的第二扫描点特征值和第三扫描点特征值及其两者模态序列号之间的差值，并确定两次扫描区间内谐振点的依据，根据谐振点的依据的数值确定两次扫描区间内是否存在谐振点；将所述第三扫描点特征值按照频率大小顺序存入所述存储特征值矩阵；
10.当所述谐振点的依据大于0时，判定两次扫描区间内不存在谐振点，根据所述第一扫描点和所述第二扫描点模态序列号之间的差值、所述第二扫描点和所述第三扫描点模态序列号之间的差值确定增加大步长或小步长。
11.优选地，所述对所述数据信息中的各参数进行初始化具体包括：
12.设置起始频率f0，设置终止频率f
end
，设置初始步长k0＝0.004pu，设置最大步长k
max
＝0.06pu，设置模态总数num，设置存储扫描点矩阵h，设置存储特征值矩阵ym，设置模态序号n＝1，设置黄金分割结果精度error＝0.002，设置两次扫描的初始步长k1＝k2＝k0，f1＝f0，其中，k1第1次扫描的步长，k2是第2次扫描的步长。
13.优选地，所述根据所述初始步长和所述初始频率确定第一扫描点，根据所述第一扫描点和所述第一次步长确定第二扫描点具体包括：
14.令f2＝f1+k1，将f2按照先后顺序存入存储扫描点矩阵h中，对f1和f2处的节点导纳矩阵做特征值分解，通过跟踪频率确定f1处的特征值yx(f1)，并得到f2处的特征值yx(f2)，其中，f2为第二扫描点，f1为第一扫描点。
15.优选地，还包括：
16.根据以下公式确定第一扫描点特征值和第二扫描点特征值两者模态序列号之间的差值：
17.δ1＝y(f2)-y(f1)
18.其中，δ1为第一扫描点特征值和第二扫描点特征值两者模态序列号之间的差值，y(f2)为f2处模态n的特征值，y(f1)为f1处模态n的特征值。
19.优选地，所述根据第二扫描点和所述第二次步长确定第三扫描点具体包括：
20.限制两次扫描的步长在最大步长内，令f3＝f2+k2，其中，f3为第三扫描点。
21.优选地，还包括：
22.当所述第三扫描点处的频率超过所述终止频率时，在特征值矩阵ym的基础上根据黄金分割结果精度确定剩余模态的谐振频率。
23.优选地，根据以下公式确定对应的第二扫描点特征值和第三扫描点特征值两者模态序列号之间的差值：
24.δ2＝y(f3)-y(f2)
25.其中，δ2为第二扫描点特征值和第三扫描点特征值两者模态序列号之间的差值，y(f3)为f3处模态n的特征值，y(f2)为f2处模态n的特征值。
26.优选地，所述谐振点的依据根据以下公式确定：
27.d＝δ1*δ228.其中，d为谐振点的依据，δ2为第二扫描点特征值和第三扫描点特征值两者模态序列号之间的差值，δ1为第一扫描点特征值和第二扫描点特征值两者模态序列号之间的差值。
29.优选地，所述根据所述第一扫描点和所述第二扫描点模态序列号之间的差值、所
述第二扫描点和所述第三扫描点模态序列号之间的差值确定增加大步长或小步长具体包括：
30.当|δ2|》1或|δ1|》1时，判定上升或下降趋势迅速，并根据以下公式增加大步长：
[0031][0032]
其中，k1第1次扫描的步长，k2是第2次扫描的步长，k
21
是增加大步长后的第2次扫描的步长；
[0033]
当|δ2|≤1且|δ1|≤1时，根据以下公式增加小步长：
[0034][0035]
其中，k1第1次扫描的步长，k2是第2次扫描的步长，k
22
是增加小步长后的第2次扫描的步长。
[0036]
优选地，还包括：
[0037]
当所述谐振点的依据小于等于0时，判定两次扫描区间内存在谐振点，使用黄金分割结果精度确定谐振频率，并根据所述初始步长和所述初始频率确定第一扫描点，根据所述第一扫描点和所述第一次步长确定第二扫描点，以及根据第二扫描点和所述第二次步长确定第三扫描点，将所述第三扫描点存入所述存储扫描点矩阵，判断所述第三扫描点处的频率是否超过所述终止频率；
[0038]
直到所述谐振点的依据大于0时，判定两次扫描区间内不存在谐振点，根据所述第一扫描点和所述第二扫描点模态序列号之间的差值、所述第二扫描点和所述第三扫描点模态序列号之间的差值确定增加大步长或小步长为止。
[0039]
综上，实施本发明的实施例，具有如下的有益效果：
[0040]
本发明提供的电力系统的谐波谐振模态的分析方法，根据所述初始步长和所述初始频率确定第一扫描点，根据所述第一扫描点和所述第一次步长确定第二扫描点，根据第二扫描点和所述第二次步长确定第三扫描点，当所述谐振点的依据大于0时，判定两次扫描区间内不存在谐振点，根据所述第一扫描点和所述第二扫描点模态序列号之间的差值、所述第二扫描点和所述第三扫描点模态序列号之间的差值确定增加大步长或小步长，可以在保证谐振频率计算精度的同时大幅提高电力系统谐波谐振模态分析速度。
附图说明
[0041]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，根据这些附图获得其他的附图仍属于本发明的范畴。
[0042]
图1为本发明实施例中一种电力系统的谐波谐振模态的分析方法的主流程示意图。
[0043]
图2为本发明实施例中一种电力系统的谐波谐振模态的分析方法的逻辑示意图。
具体实施方式
[0044]
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步地详细描述。
[0045]
如图1和图2所示，为本发明提供的一种电力系统的谐波谐振模态的分析方法的一个实施例的示意图。在该实施例中，所述方法包括以下步骤：
[0046]
步骤s1，获取所述电力系统的数据信息，并对所述数据信息中的各参数进行初始化，其中，初始化后的各参数中至少包括扫描点的初始频率、扫描点的终止频率、扫描点的初始步长、用于处理扫描点的存储特征值矩阵、进行扫描的连续两次扫描的第一次步长和第二次步长、存储扫描点的存储扫描点矩阵。
[0047]
具体实施例中，设置起始频率f0，设置终止频率f
end
，设置初始步长k0＝0.004pu，设置最大步长k
max
＝0.06pu，设置模态总数num，设置存储扫描点矩阵h，设置存储特征值矩阵ym，设置模态序号n＝1，设置黄金分割结果精度error＝0.002，设置两次扫描的初始步长k1＝k2＝k0，f1＝f0，其中，k1第1次扫描的步长，k2是第2次扫描的步长。其中，存储扫描点矩阵h(行数为1，列数随着扫描点的增加而增加)，存储特征值矩阵ym(行数为num，列数随着扫描点的增加而增加，且列数＝扫描点数量)。
[0048]
步骤s2，根据所述初始步长和所述初始频率确定第一扫描点，根据所述第一扫描点和所述第一次步长确定第二扫描点，将所述第二扫描点存入所述存储扫描点矩阵；确定对应的第一扫描点特征值和第二扫描点特征值及其两者模态序列号之间的差值，并将所述第一扫描点特征值和所述第二扫描点特征值按照所述第一扫描点和所述第二扫描点的顺序存入所述存储特征值矩阵；也就是，对f1和f2处的节点导纳矩阵做特征值分解，通过跟踪频率f1处的特征值yx(f1)得到yx(f2)，其中(yx(f))
num
×1是频率f点处各模态的特征值，计算f2处和f1处模态n的特征值差值δ1，将yx(f1)、yx(f2)按照f1、f2的顺序存入矩阵ym中。
[0049]
具体实施例中，令f2＝f1+k1，将f2按照先后顺序存入存储扫描点矩阵h中，对f1和f2处的节点导纳矩阵做特征值分解，通过跟踪频率确定f1处的特征值yx(f1)，并得到f2处的特征值yx(f2)，其中，f2为第二扫描点，f1为第一扫描点。
[0050]
具体地，根据以下公式确定第一扫描点特征值和第二扫描点特征值两者模态序列号之间的差值：
[0051]
δ1＝y(f2)-y(f1)
[0052]
其中，δ1为第一扫描点特征值和第二扫描点特征值两者模态序列号之间的差值，y(f2)为f2处模态n的特征值，y(f1)为f1处模态n的特征值。
[0053]
步骤s3，根据第二扫描点和所述第二次步长确定第三扫描点，将所述第三扫描点存入所述存储扫描点矩阵，判断所述第三扫描点处的频率是否超过所述终止频率；也就是，将f3存入h中，判断f3是否超过终止频率，当f3《f
end
时，转到步骤s4；当f3》f
end
时，在特征值矩阵ym的基础上使用黄金分割法确定剩余模态的谐振频率。
[0054]
具体实施例中，限制两次扫描的步长在最大步长内，当k1》k
max
时，令k1＝k
max
，当k2》k
max
时，令k2＝k
max
，令f3＝f2+k2，其中，f3为第三扫描点。
[0055]
本实施例中，当所述第三扫描点处的频率超过所述终止频率时，在特征值矩阵ym的基础上根据黄金分割结果精度确定剩余模态的谐振频率。可理解的，当f3》f
end
时，在特征值矩阵ym的基础上使用黄金分割法确定剩余模态的谐振频率：令，n＝n+1，利用matlab中
min函数找到ym中模态n的近似谐振点，根据谐振点在ym中对应列的序号记录谐振点在h中对应的频率，记为h(k)，k为谐振点在ym中对应列的序号；
[0056]
记录谐振频率的前一个和后一个频率，记为h(k-1)和h(k+1)，将h(k-1)和h(k+1)作为谐振频率所在的区间，以该区间作为黄金分割法的初始区间，如下式所示：
[0057][0058]
确定模态n的谐振频率和对应的特征值，当n《num时，转到上述步骤s3，当n≥num时，即已经确定剩余模态的谐振频率，停止迭代。此处根据黄金分割结果精度确定剩余模态的谐振频率可理解为在特征值矩阵ym的基础上使用黄金分割法确定剩余模态的谐振频率，黄金分割法为本领域常用方法，在此不再赘述。
[0059]
步骤s4，当所述第三扫描点处的频率未超过所述终止频率时，确定对应的第二扫描点特征值和第三扫描点特征值及其两者模态序列号之间的差值，并确定两次扫描区间内谐振点的依据，根据谐振点的依据的数值确定两次扫描区间内是否存在谐振点；将所述第三扫描点特征值按照频率大小顺序存入所述存储特征值矩阵；也就是，对f3处的节点导纳矩阵做特征值分解，通过跟踪频率f2处的特征值yx(f2)得到yx(f3)，计算f3处和f2处模态n的特征值差值δ2，得出判断两次扫描区间内是否存在谐振点的依据d，将yx(f3)按照频率大小顺序存入矩阵ym中。
[0060]
具体实施例中，根据以下公式确定对应的第二扫描点特征值和第三扫描点特征值两者模态序列号之间的差值：
[0061]
δ2＝(f3)-(f2)
[0062]
其中，δ2为第二扫描点特征值和第三扫描点特征值两者模态序列号之间的差值，y(f3)为f3处模态n的特征值，y(f2)为f2处模态n的特征值。
[0063]
具体地，所述谐振点的依据根据以下公式确定：
[0064]
d＝δ1*2[0065]
其中，d为谐振点的依据，δ2为第二扫描点特征值和第三扫描点特征值两者模态序列号之间的差值，δ1为第一扫描点特征值和第二扫描点特征值两者模态序列号之间的差值。
[0066]
步骤s5，当所述谐振点的依据大于0时，判定两次扫描区间内不存在谐振点，，根据所述第一扫描点和所述第二扫描点模态序列号之间的差值、所述第二扫描点和所述第三扫描点模态序列号之间的差值确定增加大步长或小步长。当d》0时，即两次扫描区间内模态n的特征值皆在上升或下降趋势，当|δ2|》1或|δ1|》1时，增加大步长，当|δ2|≤1且|δ1|≤1时，增加小步长，令δ1＝δ2，转到步骤s3；若d≤0，即系统在两次扫频区间出现了谐振频率，模态n的特征值出现极值点，则使用黄金分割法确定谐振频率，转到步骤s2。
[0067]
具体实施例中，当|δ2|》1或|δ1|》1时，判定上升或下降趋势迅速，并根据以下公式增加大步长：
[0068][0069]
其中，k1第1次扫描的步长，k2是第2次扫描的步长，k
21
是增加大步长后的第2次扫
描的步长；
[0070]
当|δ2|≤1且|δ1|≤1时，根据以下公式增加小步长：
[0071][0072]
其中，k1第1次扫描的步长，k2是第2次扫描的步长，k
22
是增加小步长后的第2次扫描的步长。
[0073]
具体地，当所述谐振点的依据小于等于0时，判定两次扫描区间内存在谐振点，即当d≤0时，即模态n的特征值出现极值点，使用黄金分割法确定极值点即可确定谐振频率，并根据所述初始步长和所述初始频率确定第一扫描点，根据所述第一扫描点和所述第一次步长确定第二扫描点，以及根据第二扫描点和所述第二次步长确定第三扫描点，将所述第三扫描点存入所述存储扫描点矩阵，判断所述第三扫描点处的频率是否超过所述终止频率；
[0074]
当所述谐振点的依据大于0时，判定两次扫描区间内不存在谐振点，根据所述第一扫描点和所述第二扫描点模态序列号之间的差值、所述第二扫描点和所述第三扫描点模态序列号之间的差值确定增加大步长或小步长。
[0075]
可理解的，当d≤0时，将这两次扫频区间作为一个新的区间[a，b]，取两个试验点a1，a2，如下式所示：
[0076][0077]
其中r为区间缩短率，r＝0.618，对a1和a2处的节点导纳矩阵做特征值分解，通过跟踪频率a处的特征值yx(a)得到yx(a1)，通过跟踪频率a1处的特征值yx(a1)得到yx(a2)；
[0078]
比较a2处和a1处模态n的特征值大小；
[0079]
当y(a1)《y(a2)时，保留区间[a1，b]，对原来的a、a1做出如下式所示的调整，并取新的试验点a2，对y(a)、y(a1)做出如下式所示的调整，通过跟踪频率a1处的特征值yx(a1)得到yx(a2)；
[0080][0081][0082]
当y(a1)》y(a2)时，保留区间[a，a2]，对原来的a2、b做出如下式所示的调整，并取新的试验点a1；对y(a2)、y(b)做出如下式所示的调整，通过跟踪频率a处的特征值yx(a)得到yx(a1)；
[0083][0084][0085]
当最终|a-b|《error时，停止迭代，输出谐振点x
p
，如下式所示，通过跟踪频率a1处的特征值yx(a1)得到yx(x
p
)，输出谐振点对应的特征值y(x
p
)；
[0086]
x
p
＝(a+b)/2
[0087]
令f1＝f3，y(f1)＝y(f3)，转到步骤s2。
[0088]
综上，实施本发明的实施例，具有如下的有益效果：
[0089]
本发明提供的电力系统的谐波谐振模态的分析方法，根据所述初始步长和所述初始频率确定第一扫描点，根据所述第一扫描点和所述第一次步长确定第二扫描点，根据第二扫描点和所述第二次步长确定第三扫描点，再确定两次扫描区间内谐振点的依据，根据所述第一扫描点和所述第二扫描点模态序列号之间的差值、所述第二扫描点和所述第三扫描点模态序列号之间的差值确定增加大步长或小步长，可以在保证谐振频率计算精度的同时大幅提高电力系统谐波谐振模态分析速度。
[0090]
以上所揭露的仅为本发明较佳实施例而已，当然不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。

技术特征：
1.一种电力系统的谐波谐振模态的分析方法，其特征在于，包括：获取所述电力系统的数据信息，并对所述数据信息中的各参数进行初始化，其中，初始化后的各参数中至少包括扫描点的初始频率、扫描点的终止频率、扫描点的初始步长、用于处理扫描点的存储特征值矩阵、进行扫描的连续两次扫描的第一次步长和第二次步长、存储扫描点的存储扫描点矩阵；根据所述初始步长和所述初始频率确定第一扫描点，根据所述第一扫描点和所述第一次步长确定第二扫描点，将所述第二扫描点存入所述存储扫描点矩阵；确定对应的第一扫描点特征值和第二扫描点特征值及其两者模态序列号之间的差值，并将所述第一扫描点特征值和所述第二扫描点特征值按照所述第一扫描点和所述第二扫描点的顺序存入所述存储特征值矩阵；根据第二扫描点和所述第二次步长确定第三扫描点，将所述第三扫描点存入所述存储扫描点矩阵，判断所述第三扫描点处的频率是否超过所述终止频率；当所述第三扫描点处的频率未超过所述终止频率时，确定对应的第二扫描点特征值和第三扫描点特征值及其两者模态序列号之间的差值，并确定两次扫描区间内谐振点的依据，根据谐振点的依据的数值确定两次扫描区间内是否存在谐振点；将所述第三扫描点特征值按照频率大小顺序存入所述存储特征值矩阵；当所述谐振点的依据大于0时，判定两次扫描区间内不存在谐振点，根据所述第一扫描点和所述第二扫描点模态序列号之间的差值、所述第二扫描点和所述第三扫描点模态序列号之间的差值确定增加大步长或小步长。2.如权利要求1所述的方法，其特征在于，所述对所述数据信息中的各参数进行初始化具体包括：设置起始频率f0，设置终止频率f
end
，设置初始步长k0＝0.004pu，设置最大步长k
max
＝0.06pu，设置模态总数num，设置存储扫描点矩阵h，设置存储特征值矩阵ym，设置模态序号n＝1，设置黄金分割结果精度error＝0.002，设置两次扫描的初始步长k1＝k2＝k0，f1＝f0，其中，k1第1次扫描的步长，k2是第2次扫描的步长。3.如权利要求2所述的方法，其特征在于，所述根据所述初始步长和所述初始频率确定第一扫描点，根据所述第一扫描点和所述第一次步长确定第二扫描点具体包括：令f2＝f1+k1，将f2按照先后顺序存入存储扫描点矩阵h中，对f1和f2处的节点导纳矩阵做特征值分解，通过跟踪频率确定f1处的特征值yx(f1)，并得到f2处的特征值yx(f2)，其中，f2为第二扫描点，f1为第一扫描点。4.如权利要求3所述的方法，其特征在于，还包括：根据以下公式确定第一扫描点特征值和第二扫描点特征值两者模态序列号之间的差值：δ1＝y(f2)-y(f1)其中，δ1为第一扫描点特征值和第二扫描点特征值两者模态序列号之间的差值，y(f2)为f2处模态n的特征值，y(f1)为f1处模态n的特征值。5.如权利要求2所述的方法，其特征在于，所述根据第二扫描点和所述第二次步长确定第三扫描点具体包括：限制两次扫描的步长在最大步长内，令f3＝f2+k2，其中，f3为第三扫描点。
6.如权利要求2所述的方法，其特征在于，还包括：当所述第三扫描点处的频率超过所述终止频率时，在特征值矩阵ym的基础上根据黄金分割结果精度确定剩余模态的谐振频率。7.如权利要求2所述的方法，其特征在于，根据以下公式确定对应的第二扫描点特征值和第三扫描点特征值两者模态序列号之间的差值：δ2＝(f3)-(f2)其中，δ2为第二扫描点特征值和第三扫描点特征值两者模态序列号之间的差值，y(f3)为f3处模态n的特征值，y(f2)为f2处模态n的特征值。8.如权利要求2所述的方法，其特征在于，所述谐振点的依据根据以下公式确定：d＝δ1*2其中，d为谐振点的依据，δ2为第二扫描点特征值和第三扫描点特征值两者模态序列号之间的差值，δ1为第一扫描点特征值和第二扫描点特征值两者模态序列号之间的差值。9.如权利要求8所述的方法，其特征在于，所述根据所述第一扫描点和所述第二扫描点模态序列号之间的差值、所述第二扫描点和所述第三扫描点模态序列号之间的差值确定增加大步长或小步长具体包括：当|δ2|>1或|δ1|>1时，判定上升或下降趋势迅速，并根据以下公式增加大步长：其中，k1第1次扫描的步长，k2是第2次扫描的步长，k
21
是增加大步长后的第2次扫描的步长；当|δ2|≤1且|δ1|≤1时，根据以下公式增加小步长：其中，k1第1次扫描的步长，k2是第2次扫描的步长，k
22
是增加小步长后的第2次扫描的步长。10.如权利要求8所述的方法，其特征在于，还包括：当所述谐振点的依据小于等于0时，判定两次扫描区间内存在谐振点，使用黄金分割结果精度确定谐振频率，再根据所述初始步长和所述初始频率确定第一扫描点，根据所述第一扫描点和所述第一次步长确定第二扫描点，以及根据第二扫描点和所述第二次步长确定第三扫描点，将所述第三扫描点存入所述存储扫描点矩阵，判断所述第三扫描点处的频率是否超过所述终止频率；直到当所述谐振点的依据大于0时，判定两次扫描区间内不存在谐振点，根据所述第一扫描点和所述第二扫描点模态序列号之间的差值、所述第二扫描点和所述第三扫描点模态序列号之间的差值确定增加大步长或小步长为止。

技术总结
本发明提供一种电力系统的谐波谐振模态的分析方法，包括，获取所述电力系统的数据信息，对各参数进行初始化；确定第一扫描点和第二扫描点；确定第一扫描点特征值和第二扫描点特征值两者模态序列号之间的差值；根据第二扫描点和所述第二次步长确定第三扫描点，判断所述第三扫描点处的频率是否超过所述终止频率；当所述第三扫描点处的频率未超过所述终止频率时，确定第二扫描点特征值和第三扫描点特征值两者模态序列号之间的差值，确定两次扫描区间内谐振点的依据；根据谐振点的依据的数值确定两次扫描区间内是否存在谐振点，从而确定增加大步长或小步长。本发明可在保证谐振频率计算精度的同时提高电力系统谐波谐振模态分析速度。速度。速度。


技术研发人员：张华赢 余鹏 汪清 李艳 胡智敏 李朝阳
受保护的技术使用者：深圳供电局有限公司
技术研发日：2023.05.24
技术公布日：2023/8/28