1.一种电力系统暂态稳定分析方法，其特征在于，包括以下步骤：
步骤1)利用泰勒级数对原电力系统的系统方程进行重构；
步骤2)利用有理李雅普诺夫函数扩大估计稳定域；
步骤3)利用SMR技术构造凸优化模型；
步骤4)寻求最优李雅普诺夫函数。
2.根据权利要求1所述的一种电力系统暂态稳定分析方法，其特征在于，所述步骤1)将多项式函数和非多项式函数分离，然后利用泰勒级数对非多项式部分近似展开，将系统(1-
1)写成(1-2)等效形式：
其中，h(x(t))，g(x(t))是向量多项式函数并且属于多项式集合Pn，ζ1(x(t)),...,ζr(x(t))代表非多项式函数； 在D内可解析；
令(1-3)成立：
|α|＝α1+…+αn,α！＝α1！...αn！,xα＝x1α...xnα  (1-3)
其中，x∈Rn并且α＝(α1,....,αn)T∈Nn是n维向量；k阶导数表示为式(1-4)所示：
因此，式(1-2)中 写成如下所示的泰勒展开式：
其中，ξi是有界参数，k代表幂级数，β为满足式|β|＝k+1的参数；ηi(x)是k阶泰勒多项式，带佩余亚诺相的泰勒展开式为：
ξi为泰勒余项值 的值。
3.根据权利要求2所述的一种电力系统暂态稳定分析方法，其特征在于，所述步骤2)利用有理李雅普诺夫函数扩大估计稳定域，首先，定义V(x)为系统(2-1)的有理李雅普诺夫函数，即：
其中，Vnum、Vden均属于多项式集合P，并且满足(2-2)～(2-3)所示的条件：
在获取系统稳定域时，首先定义V(x)一个子集υ(c)＝{x∈Rn:V(x)≤c}，通过求解式(2-
4)的优化模型来寻找最优李雅普诺夫函数v(x)；
μ＝supρ(v(c))
通过求解(2-5)最优问题获得有理李雅普诺夫函数估计最大稳定域ck；
其中， 的变量在零点处-ψ(x,c,s(x),ξ)的值为0，并且多项式-ψ(x,
c,s(x),ξ)是由每个单项式平方和相加组成的，构成方法如式(2-6)～(2-8)所示：
ψ(x,c,s(x),ξ)＝r(x)+q(x)ξ+s(x)(cVden(x)-Vnum(x))  (2-8)
4.根据权利要求3所述的一种电力系统暂态稳定分析方法，其特征在于，所述步骤3)中SOS模型，可以利用SMR技术对局部SOS模型进行处理，将非凸优化模型转化为凸优化模型，SMR技术处理SOS优化模型如下：
s(x)＝(φ(n,d(q)))TSφ(n,d(q))
ψ(x,c,s(x),ξ)＝(ψ(c,S,ξ))T(ψ(c,S,ξ)+L(γ))φ(n,d(ψ))  (3-1)u(x)＝u1(x)+u2(x)，u1(x)＝-r(x)-q(x)Tξ+s(x)Vnum(x)  (3-2)
其中，R(ξ),W(S),U2(S)和 分别是-r(x)-q(x)Tξ,s(x)Vnum(x),u2(x)和V(x)的SMR矩阵；
λ为正数；d(q)为不小于degx(q)/2的最小整数，degx(q)/2为多项式函数q(x)∈P0SOS的最高次数，φ(n,d(q))是由不同幂次数的变量组成的向量，变量的幂次数均小于等于d(q)的正整数，n是变量的个数，c是稳定域边界，x是函数变量；L(γ)为仿射空间，满足：
对于最优有理李雅普诺夫函数的确定及最大稳定域ck的获得通过求(3-5)实现：
其中，是式(3-6)所示GEVP(广义特征值最小化问题)问题的解，其模型如下所示：
5.根据权利要求4所述的一种电力系统暂态稳定分析方法，其特征在于，步骤4)首先通过式(4-1)获得初始有理李雅普诺夫函数V0(x)：
其中，Vq(x)是系统线性部分的二次型李雅普诺夫函数，Va(x)是辅助多项式函数，其选择方法为(xTx)·(xTPx)，在v(c)范围内通过扩大区域多项式所包围的区域来寻求最优有理李雅普诺夫函数，具体如式(4-2)所示：
其中， 是多项式；利用式(4-3)寻求最优有理李雅普诺
夫函数并获取最大的稳定域：
6.根据权利要求1所述的一种电力系统暂态稳定分析方法，其特征在于，还包括步骤
5)，分析单机无穷大系统，获取VSG策略下计及逆变器饱和的风电系统最优李雅普诺夫函数。