弟一讲微分方程.docVIP

第一节 微分方程的数值计算 1-1-1 引言 电力系统的动态过程通常用微分方程来描述。在微分方程中，自变量是时间，因变量是系统中各物理变量。对于阶常微分方程，可以采用引入新变量的方法把它转化为个一阶常微分方程来求解。例如，对于二阶常微分方程 可以表示成： 一般来说，依据电工学和动力学建立起来的电力系统数学模型是可以用来数字仿真的。但是，在复杂电力系统中，由于数学模型很复杂，仿真的计算量很大，影响了仿真效率。有些环节对于所研究的问题影响很小，在建立数学模型时可以把它们忽略。这种简化严格来说应该用实际实验，或者用精确的物理模拟实验来检验，以确定简化的合理性。数学模型中的参数是数字仿真的基础，一般可以通过各种测量和实验得到，但这是一件十分烦琐的工作。目前，有些参数只能凭经验给出，其合理性有待实验的严格检验。 电力系统仿真的另一个重要问题是数学模型的求解问题。对于微分方程，除少数可以得到解析解以外，大多数只能采用数值解法。早在18世纪末，就有很多人提出过求解微分方程精确而有效的数值积分方法。 微分方程的积分是一簇曲线。通常，在初值确定之后，数值解才能够确定。 得到微分方程的数值近似解有两种基本的方法。一种方法是把近似解表示成有限个独立函数之和，另一种方法是差分法。 差分方法是寻求在一系列离散点上的近似值，这些离散点称为结点。在多数情况下，这些结点是等距的，即 称为步长。差分方法是一种递推算法，它使求解过程能顺着结点的顺序一步一步向前推进，即可用前一个结点上的值(单步法)或者前面几个结点上的值(多步法)来计算当前步上的近似值。 1-1-2 单步法 所谓单步法就是用前一结点上的值来计算当前结点上的近似值的方法。 1．欧拉法 欧拉法是一种最简单的单步法。 对一阶微分方程 假定已给定，可计算出。把一阶微分方程写成数值积分法的计算形式 如果十分靠近，则有 当然，也可以把一阶微分方程写成 即用时的导数表示在区间[，]内的平均增量。结点和结点之间的关系 （1-7） 尽管这种方法精度不高，但却提出了一种设想，可以通过递推计算方法得出微分方程的近似解，这种方法称为欧拉法。因为是用折线来近似表示曲线，故又可称为折线法。 2．隐式梯形法 从欧拉法可以看出，为提高精度，可以用和的平均值作为在区间[，]内的平均增量值。 （1-9） 在等式的两边均含有未知量，故称为隐式梯形法。隐式梯形公式不能用递推的方式直接做数值计算。如果用欧拉法求得的值作为预测值，则有： 这种方法称为改进欧拉法。 利用欧拉法和改进欧拉法求微分方程的数值解，如果选取步长相同，改进欧拉法的计算量大，但改进欧拉法的精度高。反过来说，如果两种方法要求精度相同，则改进欧拉法可以选取较大步长，总的计算量可以节省，舍入误差也可较小。 3．龙格一库塔法 隐式梯形法是把和的平均值作为在区间[，]内的平均增量值。同理可以设想在[，]区间内，取不同的后再求加权平均值，用它作为该区间的平均增量，即 其中 这种计算式称为自启动的N阶龙格一库塔法。当N=l时为欧拉法，当N=2时为改进欧拉法。常采用的各阶龙格一库塔法计算公式的系数如表l所示。 表1 各阶龙格一库塔法的系数表 通常采用的名称 1 0 欧拉法 2 改进欧拉法 3 三阶龙格一库塔法 4 四阶龙格一库塔法 四阶龙格一库塔法 其中 1-1-3多步法 多步法的一般计算公式为 其中 和 如果为显式公式，否则是隐式的。 这种方法不仅用到前面一个结点上的值，而且还用到前面几个结点上的值。常用的方法是，其它的，系数如表2所示 通常采用的名称 显式 1 欧拉法 2 梯形法 3 三点亚当姆斯 4 四点亚当姆斯 隐式 1 欧拉法 2 梯形法 3 三点亚当姆斯 4 四点亚当姆斯 在电力系统动态仿真计算中，常采用的多步计算方法是用四点亚当姆斯显式公式计算预测值，用四点亚当姆斯隐式公式计算校正值。 为了提高精度，在误差分析的基础上，分别用亚当姆斯四点显式公式计算预测值及其修正预测值，用亚当姆斯四点隐式公式计算校正值及其修正终值。计算过程如下： ) 计算 ) 用显式公式计算预测值 ) 计算误差修正后的预测值 ，其中是在第步计算中得到的校正值。 )计算 5) 用隐式公式计算校正值 6) 计算误差修正后的终值 1-1-4数值积分法的分类 按照数值积分时被积函数的近似表达式的不同，分为单步法和多步法。若采用线性插值函数代替被积函数，则可得出单步法计算公式：若采用高阶插值函数近似代替被积函敷，可推导出多步法计算公式。从理论上讲，多步法比单步法更为有效，可取用更多结点的数据来计算，精度较高。但是在突变点，因为变