[发明专利]一种适于隐式投影算法的误差估计及参数自适应调节方法

权利要求书

1.一种适于隐式投影算法的误差估计及参数自适应调节方法，其特征在于，包括如下步骤：

1)读取电力系统拓扑连接关系、元件参数和仿真计算参数，其中仿真计算参数包括仿真终止时间T，仿真步长h，隐式投影算法内部积分器的积分步数k，隐式投影算法外部积分器步长倍数的初始值M_ini和最大值M_max，最大允许误差并设置仿真故障及操作事件；

2)设置当前隐式投影算法内部积分器的积分步数s＝1，s为正整数；

3)设置仿真时间t＝t+h，采用隐式投影算法内部积分器计算一个步长得到电力系统该时刻的状态变量x_n+s和代数变量y_n+s，并设置s＝s+1；

4)根据步骤1)设置的仿真故障及操作事件，判断电力系统当前时刻是否存在故障或操作，若存在，即故障或操作事件的发生时间T_event＝t，返回步骤2)，否则进行下一步骤；

5)判断隐式投影算法内部积分器的积分步数s是否大于k+1，是则进入下一步骤，若s≤k+1，则返回步骤3)；

6)设当前隐式投影算法外部积分器以M_n-1h作为积分步长，其中，M_n-1为上一步的隐式投影算法外部积分器步长倍数，利用事前误差估计方法对隐式投影算法内部积分器积分k+1步以及隐式投影算法外部积分器以步长M_n-1h积分一步的局部截断误差进行估计，得到事前估计误差

7)根据估计误差，估计当前隐式投影算法外部积分器步长倍数M_n；

8)根据步骤1)设置的故障及操作事件，判断t～t+M_nh时间内是否发生故障或操作，若发生，即t<T_event<t+M_nh，则进入步骤10)，否则进行下一步骤；

9)设置隐式投影算法外部积分器积分步长H＝M_nh，设置仿真时间t＝t+H，利用隐式投影算法外部积分器得到系统该时刻的状态变量和代数变量然后转入步骤11)；

10)设置隐式投影算法外部积分器外部积分步长H＝T_event-t，设置仿真时间t＝T_event，利用隐式投影算法外部积分器得到故障或操作发生前系统的状态变量和代数变量

11)通过事后误差估计方法对隐式投影算法外部积分器的事后估计误差进行估计；

12)将与进行比较，若T^n+k+1+Mn<aposteriori>>1.1T^tol,]]>则拒绝本次隐式投影算法外部积分器的积分计算，设置仿真时间t＝t-H，返回步骤7)，否则进入下一步骤；

13)判断仿真时间t是否达到仿真终止时间T，若t＝T，则仿真结束，否则返回步骤2)，依据步骤2)至13)反复进行直至仿真结束。

2.根据权利要求1所述的一种适于隐式投影算法的误差估计及参数自适应调节方法，其特征在于，步骤1)所述的隐式投影算法内部积分器采用基于显式改进欧拉法的微分代数方程交替求解方法，积分步数k要求大于等于3，隐式投影算法外部积分器采用基于隐式预测-校正方法的微分代数方程交替求解方法进行计算，隐式投影算法外部积分器步长倍数的初始值M_ini值为1。