[发明专利]基于高阶Radau方法的高压输电线路电磁暂态快速计算方法

权利要求书

1.一种基于高阶Radau方法的高压输电线路电磁暂态快速计算方法，其特征在于包括以 下步骤：利用RadauIIA方法的系数矩阵满足V-变换这一特性，将多级高阶数值积分方法所 导致的高维线性方程组的求解问题解耦并形成多级分块递推求解方法。

2.根据权利要求1所述的基于高阶Radau方法的高压输电线路电磁暂态快速计算方法， 其特征在于所述的步骤具体如下：

1)将高压输电线路进行空间离散化，建立相应的电磁暂态数值计算的数学模型，即：

x·=f(t,x)=Bx+w(t)---(1);]]>

式(1)中，x是高压输电线路电磁暂态计算中的状态变量；B是与高压输电线路参数相关 的定常系数矩阵；w(t)即是电磁暂态计算的输入信号或激励信号；

2)电磁暂态数值计算初始化

置t＝0.0s，积分步数n＝0；确定RadauIIA方法的级数s，取s≥3；

确定数值积分步长h、电磁暂态数值计算时程T；

确定各状态变量的初值，即x(t＝0)＝x₀；

输入电磁暂态数值计算的故障或操作；

3)故障或操作判断

依据时刻t判断系统有无故障或操作：

若无故障或操作，则直接转至步骤4)；

若有故障或操作，则修改系数矩阵B以及激励信号w(t)并重新形成方程(1)；

同时，若涉及状态变量的突变，则还需依据具体情况修改相应的状态变量值x_n(t)；

4)数值积分

已知状态变量x(t)在t_n时刻的值x_n，采用s级(s≥3)的RadauIIA方法求解状态变量x(t) 在t_n+1时刻的值x_n+1；

5)t＝t+h；n＝n+1

6)数值计算是否终止判断

若t＜T，则返回步骤3)：

若t≥T，则转至步骤7)；

7)数值计算结果输出。