基于频率分解和相关系数法的电网弱阻尼辨识方法、系统及介质与流程

本发明涉及电力系统安全领域，具体涉及一种基于频率分解和相关系数法的电网弱阻尼辨识方法、系统及介质。
背景技术：
：随着我国国民经济的快速发展，电力系统规模日益增大，电网互联程度越来越高，长距离大功率送电，大量电力电子如风电、光伏、直流的大量接入等易引发电力系统发生低频功率振荡，严重的低频振荡将引发电网重大停电事故。因此，厘清电力系统低频振荡产生的机理，研究低频振荡预警以及抑制措施对电力系统稳定运行具有重要的意义。电力系统在正常运行时功率稳定，不会产生低频振荡，当系统出现扰动时，电力系统传输的功率、电压、电流以及频率等状态量将出现一定频率范围内的波动，这种波动称为低频振荡。低频振荡根据作用范围和频率大小差别可以分为局部振荡和区域振荡。其中，局部振荡一般发生在一定范围内一台电机或几台电机之间，振荡频率相对较高，通常在0.7～2.5hz之间。区域振荡指的时不同区域间发生振荡，振荡的范围较大但振荡频率较低，一般在0.1～0.7hz之间。低频振荡的机理一般有以下几种：1、弱阻尼或负阻尼机理：由于励磁系统追求快速性、电网负荷加重以及系统的连通性，使得电力系统阻尼下降，电力系统对一定频率的振荡表现出弱阻尼乃至负阻尼特性，导致振荡短时间无法消除；2、共振机理：当电力系统原动机功率遭受的周期振荡与系统固有的低频振荡频率接近或相等时，容易诱发共振，这种共振具有起振快、消失快、振荡频率与扰动频率一致的特点，影响共振的因素有阻尼转机系统、同步力矩系数、扰动幅度等；3、发电机电磁惯性导致的低频振荡：电感性的励磁绕组在励磁电压的作用下能够产生一个相位滞后的励磁电流强迫分量，在该分量的控制下会导致低频振荡的发生；4、分叉和混沌理论：分叉理论揭示了电力系统低频振荡的非线性特征，混沌理论考虑非周期性、无规则性的低频振荡参数间的相互作用。电力系统低频振荡可以从三个方面进行分析，首先根据线性系统分析，由于调节措施的影响，使得系统产生弱阻尼或负阻尼，导致系统扰动后产生振荡，且振荡幅度不衰减。其次是从输入信号或扰动信号间具有某种特定关系时，系统会产生共振或谐振，振荡幅度较大，频率低。最后，由于系统的非线性影响，参数或者扰动发生变化时，系统稳定结构也会发生变化，导致低频振荡发生。在实际系统发生的低频振荡，上述3种低频振荡机理导致的低频振荡模式可能单独或者同时发生，需要综合分析，快速找出发生低频振荡的主导模式。目前电力系统大量装备了广域同步量测系统，该系统为电力系统提供了一种高精度、高速率的同步测量方案，对于监测电力系统的动态运行特性具有重要的意义，应用的领域有电力系统故障快速反演，低频振荡监测及告警以及电力系统参数辨识等。在国内电网常见的d5000调度自动化系统中，已经装备基于广域同步量测低频振荡监测系统，该系统基于传统的普罗尼(prony)分析方法，当电网发生低频振荡，能及时辨识低频振荡的频率、幅度以及阻尼比等参数，并发出相应的低频振荡告警。该系统能及时对电力系统已经发生的低频振荡进行告警，但无法对电力系统潜在的低频振荡风险进行预警。技术实现要素：本发明要解决的技术问题：针对现有技术的上述问题，为实现电力系统低频振荡预警并抑制低频振荡的产生，本发明提出一种基于频率分解和相关系数法的电网弱阻尼辨识方法、系统及介质，本发明通过对广域同步量测系统测量值进行实时的频域分析，分析并监测电力系统存在的弱阻尼，并进一步通过相关系数法，建立系统存在的某个或某几个弱阻尼与电力系统其他状态量(如机组出力、开关的投切等)的相关性分析，从而对电力系统弱阻尼的产生源进行定位，且电力系统存在的弱阻尼是产生低频振荡的一个重要的先决条件，因此在本发明实现电网弱阻尼产生源的定位的基础上可及时对可能发生的低频振荡进行预警并进采取措施预防电力系统低频振荡的产生。为了解决上述技术问题，本发明采用的技术方案为：一种基于频率分解和相关系数法的电网弱阻尼辨识方法，实施步骤包括：1)针对电力系统建立基于广域测量系统实时测量数据的能量密度矩阵函数；2)遍历选择一个频率区段作为当前搜寻频率区段；3)在当前搜寻频率区段内运用奇异值分解法对能量密度矩阵函数进行频域分解，建立指定频域范围内的随频率变化的频域能量密度函数曲线；在频域能量密度函数曲线上选择能量密度峰值，在能量密度峰值所对应的频率附近的指定频率范围内截取能量密度函数曲线，对截取的能量密度函数曲线进行反傅里叶变换由频域变换到时域；将变换到时域的能量密度函数曲线运用普罗尼分析，基于普罗尼分析结果进行电网阻尼辨识，辨识出其中的弱阻尼模态；4)判断频率区段是否遍历完毕，尚未遍历完毕则跳转执行步骤2)；否则表示已经遍历完毕，得到的电力系统存在的弱阻尼模态的时间序列x＝{i＝1…n|xi}，跳转执行步骤5)；5)基于相关系数法对所有的弱阻尼模态的搜寻频率区段进行电网弱阻尼产生源定位。可选地，步骤5)的详细步骤包括：5.1)在得到的电力系统存在的弱阻尼模态的时间序列x＝{i＝1…n|xi}的基础上，针对电力系统建立与时间序列x＝{i＝1…n|xi}中时间对应的指定状态参量y的时间序列y＝{i＝1…n|yi}；5.2)计算时间序列x、时间序列y的相关系数γxy；5.3)根据相关系数γxy结合对照相关系数表判断时间序列x、时间序列y是否为强正相关或强负相关，如果时间序列x和时间序列y为强相关的正相关或负相关，则判定认为指定状态参量y所对应的设备与电力系统内某个频率的弱阻尼模态高度相关，并判定指定状态参量y所对应的设备为该频率的弱阻尼模态的产生源。可选地，步骤5.2)中计算时间序列x、时间序列y的相关系数γxy的函数表达式如下：上式中，xi为时间序列x中的第i个元素，yi为时间序列y中的第i个元素，为时间序列x中的元素均值，为时间序列y中的元素均值，为和的乘积，n为时间序列x以及时间序列y中的元素个数。可选地，步骤5.3)中的相关系数表包括不相关、弱相关、中相关、强相关四种相关类型，且每一种相关类型具有负相关、正相关两种子类型。可选地，步骤5.3)中的相关系数表中，不相关的负相关对应相关系数γxy的取值范围为-0.09至0.00，不相关的正相关对应相关系数γxy的取值范围为0.00至0.09；弱相关的负相关对应相关系数γxy的取值范围为-0.30至-0.10，弱相关的正相关对应相关系数γxy的取值范围为0.10至0.30；中相关的负相关对应相关系数γxy的取值范围为-0.50至-0.30，中相关的正相关对应相关系数γxy的取值范围为0.30至0.50；强相关的负相关对应相关系数γxy的取值范围为-1.00至-0.50，强相关的正相关对应相关系数γxy的取值范围为0.50至1.00。可选地，步骤1)中的电力系统为多输入多输出的非线性系统，且针对电力系统建立基于广域测量系统实时测量数据的能量密度矩阵函数的函数表达式如下式所示：上式中，表示一阶矩阵形式的广域测量系统输出向量的能量密度矩阵函数，suu(ω)表示广域测量系统输入向量的能量密度矩阵函数，ω表示频域的频率变量，ωr表示弱阻尼模态所对应的频率，ωnearωr表示频率变量位于弱阻尼模态所对应的频率附近，分别表示弱阻尼模态特征根对应的经过矩阵变换后的右特征相量所对应的矩阵，分别表示弱阻尼模态所对应的标量，分别表示的转置矩阵，其中下标中的r表示对应的弱阻尼模态，下标中的n0表示系统存在的弱阻尼个数，上标中的h表示矩阵转置运算。可选地，步骤3)中将变换到时域的能量密度函数曲线运用普罗尼分析，基于普罗尼分析结果进行电网弱阻尼辨识出其中的弱阻尼模态具体是指：将变换到时域的能量密度函数曲线运用普罗尼分析得到对应阻尼模态的阻尼比，如果阻尼模态的阻尼比低于预设阈值则初步判定该阻尼模态为弱阻尼模态。此外，本发明还提供一种基于频率分解和相关系数法的电网弱阻尼辨识系统，包括：能量密度矩阵建立程序单元，用于针对电力系统建立基于广域测量系统实时测量数据的能量密度矩阵函数；搜寻频率区段选择程序单元，用于遍历选择一个频率区段作为当前搜寻频率区段；弱阻尼模态分析程序单元，用于在当前搜寻频率区段内运用奇异值分解法对能量密度矩阵函数进行频域分解，建立指定频域范围内的随频率变化的频域能量密度函数曲线；在频域能量密度函数曲线上选择能量密度峰值，在能量密度峰值所对应的频率附近的指定频率范围内截取能量密度函数曲线，对截取的能量密度函数曲线进行反傅里叶变换由频域变换到时域；将变换到时域的能量密度函数曲线运用普罗尼分析，基于普罗尼分析结果进行电网阻尼辨识，辨识出其中的弱阻尼模态；搜寻频率区段遍历判断程序单元，用于判断频率区段是否遍历完毕，尚未遍历完毕则跳转执行搜寻频率区段选择程序单元；否则表示已经遍历完毕，得到的电力系统存在的弱阻尼模态的时间序列x＝{i＝1…n|xi}，跳转执行源定位程序单元；源定位程序单元，用于基于相关系数法对所有的弱阻尼模态的搜寻频率区段进行电网弱阻尼产生源定位。此外，本发明还提供一种基于频率分解和相关系数法的电网弱阻尼辨识系统，包括计算机设备，该计算机设备被编程或配置以执行所述基于频率分解和相关系数法的电网弱阻尼辨识方法的步骤，或者该计算机设备的存储器上存储有被编程或配置以执行所述基于频率分解和相关系数法的电网弱阻尼辨识方法的计算机程序。此外，本发明还提供一种计算机可读存储介质，该计算机可读存储介质上存储有被编程或配置以执行所述基于频率分解和相关系数法的电网弱阻尼辨识方法的计算机程序。和现有技术相比，本发明具有下述优点：为实现电力系统低频振荡预警并抑制低频振荡的产生，本发明通过对广域同步量测系统测量值进行实时的频域分析，分析并监测电力系统存在的弱阻尼，并进一步通过相关系数法，建立系统存在的某个或某几个弱阻尼与电力系统其他状态量(如机组出力、开关的投切等)的相关性分析，从而对电力系统弱阻尼的产生源进行定位。而电力系统存在的弱阻尼是产生低频振荡的一个重要的先决条件，通过定位弱阻尼产生源，及时对可能发生的低频振荡进行预警，并进采取措施预防电力系统低频振荡的产生。附图说明图1为本发明实施例方法的基本流程示意图。图2为本发明实施例方法中电力系统的模型示意图。图3为本发明实施例方法的详细流程示意图。具体实施方式如图1，本实施例基于频率分解和相关系数法的电网弱阻尼辨识方法的实施步骤包括：1)针对电力系统建立基于广域测量系统实时测量数据的能量密度矩阵函数；2)遍历选择一个频率区段作为当前搜寻频率区段；3)在当前搜寻频率区段内运用奇异值分解法对能量密度矩阵函数进行频域分解，建立指定频域范围内的随频率变化的频域能量密度函数曲线；在频域能量密度函数曲线上选择能量密度峰值，在能量密度峰值所对应的频率附近的指定频率范围内截取能量密度函数曲线，对截取的能量密度函数曲线进行反傅里叶变换由频域变换到时域；将变换到时域的能量密度函数曲线运用普罗尼分析，基于普罗尼分析结果进行电网阻尼辨识，辨识出其中的弱阻尼模态；4)判断频率区段是否遍历完毕，尚未遍历完毕则跳转执行步骤2)；否则表示已经遍历完毕，得到的电力系统存在的弱阻尼模态的时间序列x＝{i＝1…n|xi}，跳转执行步骤5)；5)基于相关系数法对所有的弱阻尼模态的搜寻频率区段进行电网弱阻尼产生源定位。本实施例中，步骤1)中的电力系统为多输入多输出的非线性系统，且针对电力系统建立基于广域测量系统实时测量数据的能量密度矩阵函数的函数表达式如下式所示：上式中，表示一阶矩阵形式的广域测量系统输出向量的能量密度矩阵函数，suu(ω)表示广域测量系统输入向量的能量密度矩阵函数，ω表示频域的频率变量，ωr表示弱阻尼模态所对应的频率，ωnearωr表示频率变量位于弱阻尼模态所对应的频率附近，分别表示弱阻尼模态特征根对应的经过矩阵变换后的右特征相量所对应的矩阵，分别表示弱阻尼模态所对应的标量，分别表示的转置矩阵，其中下标中的r表示对应的弱阻尼模态，下标中的n0表示系统存在的弱阻尼个数，上标中的h表示矩阵转置运算。本实施例中，电力系统为多输入多输出的非线性系统，如图2所示，其中a表示状态相量传输方程中状态相量对应的矩阵变换，b表示状态相量传输方程中输入相量对应的矩阵变换，c表示输出相量传输方程中状态相量对应的矩阵变换，d表示输出相量传输方程中输入相量对应的矩阵变换，表示下一时刻状态相量，δx,δu,δy分别表示电力系统的状态相量、输入相量以及输出相量，均为基于广域测量系统(wide-areameasurementsystem,wams)实时测量电压和电流相量数据。因此可建立非线性电力系统在平衡点的矩阵方程如下式所示：由此得到上述等式方程的传输矩阵方程如下式所示：上式中，h(s)为传输矩阵方程，λi＝-αi+jωi是矩阵a的特征根，其中和ψi分别是特征根λi所对应右特征相量和左特征相量，n为电力系统的模型阶数。由此得到输出相量的能量密度矩阵函数如下式所示：syy(ω)＝h(jω)·suu(ω)·h(jω)h(4)上式中，syy(ω)为输出相量能量密度矩阵函数，suu(ω)为输入相量的能量密度矩阵函数，h(jω)表示传输函数，ω表示频率变量。如果输入为白噪声信号，则suu(ω)为对角矩阵，标记为f，得到函数如下式所示：上式中，其中和ψi分别是特征根λi所对应右特征相量和左特征相量，λi＝-αi+jωi是矩阵a的特征根，λk为i等于k时的特征根λi，上标中的h表示矩阵转置运算，n为电力系统的模型阶数，*表示复数共轭。假设λr＝-αr+jωr代表弱阻尼模态，则有：上式中，表示一阶矩阵形式的广域测量系统输出向量的能量密度矩阵函数，ω表示频率变量，ωr表示弱阻尼对应的频率，ωnearωr表示频率变量位于弱阻尼模态所对应的频率附近，ar表示(5)中的中间变量，αr表示弱阻尼所对应特征根的实部。令：其中和ψi分别是特征根λi所对应右特征相量和左特征相量，则有：上式中，wr表示弱阻尼模态特征根对应的经过矩阵变换后的右特征相量所对应的矩阵，lr表示弱阻尼模态特征根对应的经过矩阵变换后的左特征相量所对应的矩阵，为lr转置矩阵，为wr转置矩阵，f表示输入相量的能量密度矩阵函数形成的对角矩阵，αr表示弱阻尼所对应特征根的实部，sr(ω)表示：因此，可知式(7)中的为一阶矩阵。如果有多个弱阻尼模态，则(7)能量密度矩阵函数可以表示为式(1)所示函数表达式，且可简写为：其中，w表示由元素构成的单行矩阵，s(ω)表示由元素构成的对角线矩阵，wh表示由元素构成的单列矩阵。本实施例中，步骤2)用于遍历选择一个频率区段作为当前搜寻频率区段，在指定的频率范围内(例如电力系统低频振荡频率区间0.1～2.0赫兹)内，以便步骤3)运用奇异值分解法(singularvectordecomposition)对能量密度矩阵函数进行频域分解，建立指定频域范围内的随频率变化能量密度曲线。如图3所示，本实施例中步骤3)的详细步骤包括：3.1)在当前搜寻频率区段内运用奇异值分解法对能量密度矩阵函数进行频域分解，建立指定频域范围内的随频率变化的频域能量密度函数曲线。针对式(1)/式(9)进行奇异值分解变换，得到：上式中，(构成u(ω))分别表示经过奇异值分解后的列相量，(构成σ(ω))分别表示弱阻尼模态对应的奇异值，(构成v(ω)h)分别表示经过奇异值分解的行相量。由于syy(ω)为哈密尔顿矩阵，因此u＝v，得到：3.2)在频域能量密度函数曲线上选择能量密度峰值，在能量密度峰值所对应的频率附近的指定频率范围内(如±0.05赫兹范围内)截取能量密度函数曲线，对截取的能量密度函数曲线进行反傅里叶变换由频域变换到时域。由式(10)和(11)得到：σ(ω)＝q(ω)s(ω)q(ω)h(12)上式中，q(ω)是满足式(9)中特征相量w＝u(ω)q(ω)的右矩阵，u(ω)表示由列相量构成的矩阵，q(ω)h表示q(ω)的转置矩阵，s(ω)表示由元素构成的对角线矩阵，则有：上式中，分别表示弱阻尼模态所对应的标量，分别表示q(ω)的第1～n0列相量，分别表示的转置矩阵，下标中的r表示对应的弱阻尼模态，下标中的n0表示系统存在的弱阻尼个数。因此有：上式中，分别表示弱阻尼模态对应的奇异值，分别表示弱阻尼模态所对应的标量，分别为q1(ω)列的第1～n0个元素，分别为列的第1～n0个元素，分别为列的第1～n0个元素，分别为列的第1～n0个元素，*表示复数共轭。如果w为正交矩阵，则q为对角矩阵，则有：上式中，q11(ω)为q1(ω)列的第1个元素，为列的第1个元素，为列的第n0个元素，为列的第n0个元素，通式可表示为qij(ω)为qi(ω)列的第j个元素,qij*(ω)为qi*(ω)列的第j个元素。令则有：上式中，αr表示弱阻尼所对应特征根的实部，lr表示弱阻尼模态特征根对应的经过矩阵变换后的左特征相量所对应的矩阵，为lr转置矩阵，f表示输入相量的能量密度矩阵函数形成的对角矩阵，ω表示频率变量，ωr表示弱阻尼模态所对应的频率。对(16)进行反傅里叶变换，得到时域函数：上式中，表示频域到时域的变换，f(t)为时域函数，αr表示弱阻尼所对应特征根的实部，lr表示弱阻尼模态特征根对应的经过矩阵变换后的左特征相量所对应的矩阵，为lr转置矩阵，f表示输入相量的能量密度矩阵函数形成的对角矩阵，ω表示频率变量，ωr表示弱阻尼模态所对应的频率。t表示时间。取t>0，时域函数f(t)可表示为指数函数形式：上式中各个符号定义详见式(17)。3.3)将变换到时域的能量密度函数曲线(即时域函数f(t))运用普罗尼(prony)分析，基于普罗尼(prony)分析结果进行电网弱阻尼辨识出其中的弱阻尼模态；本实施例中，步骤3)中将变换到时域的能量密度函数曲线运用普罗尼分析，基于普罗尼分析结果进行电网弱阻尼辨识出其中的弱阻尼模态具体是指：将变换到时域的能量密度函数曲线运用普罗尼分析得到对应阻尼模态的阻尼比，如果阻尼模态的阻尼比低于预设阈值则初步判定该阻尼模态为弱阻尼模态(可以根据需要进行进一步排查是否为这个频率的弱阻尼模态的产生源)。本实施例中，如果所对应的阻尼模态阻尼比低于5％，可以认为是弱阻尼模态，如果系统内长期存在同一弱阻尼模态(如连续超过10分钟)，进行预警。将上述已经截取的频率区段标记为已经搜寻过，对于剩下的频率区段重复步骤2)～步骤3)直到遍历完毕，最终遍历完毕后即可得到的电力系统存在的弱阻尼模态的时间序列x＝{i＝1…n|xi}，然后就可以进入步骤5)基于相关系数法对所有的弱阻尼模态的搜寻频率区段进行电网弱阻尼产生源定位。本实施例中，步骤4)中判断是否遍历完毕的具体条件为得到的能力密度峰值低于最高能力密度峰值的1/4，此外也可以根据需要采用其他判断条件。本实施例中，步骤5)的详细步骤包括：5.1)在得到的电力系统存在的弱阻尼模态的时间序列x＝{i＝1…n|xi}的基础上，针对电力系统建立与时间序列x＝{i＝1…n|xi}中时间对应的指定状态参量y的时间序列y＝{i＝1…n|yi}；5.2)计算时间序列x、时间序列y的相关系数γxy；5.3)根据相关系数γxy结合对照相关系数表判断时间序列x、时间序列y是否为强正相关或强负相关，如果时间序列x和时间序列y为强相关的正相关或负相关，则判定认为指定状态参量y所对应的设备与电力系统内某个频率的弱阻尼模态高度相关，并判定指定状态参量y所对应的设备为该频率的弱阻尼模态的产生源。本实施例中，步骤5.2)中计算时间序列x、时间序列y的相关系数γxy的函数表达式如下：上式中，xi为时间序列x中的第i个元素，yi为时间序列y中的第i个元素，为时间序列x中的元素均值，为时间序列y中的元素均值，为和的乘积，n为时间序列x以及时间序列y中的元素个数。本实施例中，步骤5.3)中的相关系数表包括不相关、弱相关、中相关、强相关四种相关类型，且每一种相关类型具有负相关、正相关两种子类型。本实施例中，步骤5.3)中的相关系数表中，不相关的负相关对应相关系数γxy的取值范围为-0.09至0.00，不相关的正相关对应相关系数γxy的取值范围为0.00至0.09；弱相关的负相关对应相关系数γxy的取值范围为-0.30至-0.10，弱相关的正相关对应相关系数γxy的取值范围为0.10至0.30；中相关的负相关对应相关系数γxy的取值范围为-0.50至-0.30，中相关的正相关对应相关系数γxy的取值范围为0.30至0.50；强相关的负相关对应相关系数γxy的取值范围为-1.00至-0.50，强相关的正相关对应相关系数γxy的取值范围为0.50至1.00，具体如表1所示：表1：相关系数表。相关性强度负相关(相关系数)正相关(相关系数)不相关-0.09至0.000.00至0.09弱相关-0.3至-0.10.1至0.3中相关-0.5至-0.30.3至0.5强相关-1.0至-0.50.5至1.0综上所述，本实施例提出了一种基于频域分解和相关系数的电网弱阻尼产生源定位方法，该方法通过对广域同步相量测量系统(wide-areameasurementsystem,wams)实时测量数据进行频域分解实现电网弱阻尼的辨识，并在辨识的基础上采用相关系数法实现电网弱阻尼产生源的定位。传统的基于普罗尼(prony)算法可以实现电网低频振荡的告警，但无法进行预警，而电网存在弱阻尼是发生低频振荡的必要条件之一。本实施例提出的基于频域分解和相关系数的电网弱阻尼产生源定位方法可以实现电网弱阻尼产生源的定位，从而可以提前采取措施消除弱阻尼产生源，达到减少电网产生低频振荡的概率。此外，本实施例还提供一种基于频率分解和相关系数法的电网弱阻尼辨识系统，包括：能量密度矩阵建立程序单元，用于针对电力系统建立基于广域测量系统实时测量数据的能量密度矩阵函数；搜寻频率区段选择程序单元，用于遍历选择一个频率区段作为当前搜寻频率区段；弱阻尼模态分析程序单元，用于在当前搜寻频率区段内运用奇异值分解法对能量密度矩阵函数进行频域分解，建立指定频域范围内的随频率变化的频域能量密度函数曲线；在频域能量密度函数曲线上选择能量密度峰值，在能量密度峰值所对应的频率附近的指定频率范围内截取能量密度函数曲线，对截取的能量密度函数曲线进行反傅里叶变换由频域变换到时域；将变换到时域的能量密度函数曲线运用普罗尼分析，基于普罗尼分析结果进行电网阻尼辨识，辨识出其中的弱阻尼模态；搜寻频率区段遍历判断程序单元，用于判断频率区段是否遍历完毕，尚未遍历完毕则跳转执行搜寻频率区段选择程序单元；否则表示已经遍历完毕，得到的电力系统存在的弱阻尼模态的时间序列x＝{i＝1…n|xi}，跳转执行源定位程序单元；源定位程序单元，用于基于相关系数法对所有的弱阻尼模态的搜寻频率区段进行电网弱阻尼产生源定位。此外，本实施例还提供一种基于频率分解和相关系数法的电网弱阻尼辨识系统，包括计算机设备，该计算机设备被编程或配置以执行前述基于频率分解和相关系数法的电网弱阻尼辨识方法的步骤，或者该计算机设备的存储器上存储有被编程或配置以执行前述基于频率分解和相关系数法的电网弱阻尼辨识方法的计算机程序。此外，本实施例还提供一种计算机可读存储介质，该计算机可读存储介质上存储有被编程或配置以执行前述基于频率分解和相关系数法的电网弱阻尼辨识方法的计算机程序。本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。以上所述仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。当前第1页1 2 3