第四章连续系统的离散化方法.pptVIP

第一页，共三十八页，2022年，8月28日 其一般公式为 称为截断误差 例4－1 用欧拉法求下述微分方程的数值解。 第二页，共三十八页，2022年，8月28日 解：因欧拉法的递推公式为 所以 则递推公式为： 若取步长 由t=0开始计算可得 上式的精确解是 数值计算与精确解的比较见表 t 0 0.1 0.2 0.3 1.0 精确解x(t) 1 0.9090909 0.8333333 0.7692307 0.5 1 0.9 0.819 0.5719 0.4627810 第三页，共三十八页，2022年，8月28日 4.1.3 龙格－库塔法 将 在点 处作台劳级数展开，并取线性部分可得 第四页，共三十八页，2022年，8月28日 将 代入式 比较 各项系数得 待定系数个数超过方程个数，必须先设定一个系数，然后即可求得其 参数。一般有以下几种取法： 1、 则 一般形式 第五页，共三十八页，2022年，8月28日 2、 则 一般形式 3、 则 一般形式 以上几种递推公式均称为二阶龙格－库塔公式。是较典型的几 个常用算法。其中的方法3又称为预估－校正法，或梯形法。 第六页，共三十八页，2022年，8月28日 其意义如下：用欧拉法以斜率先求取一点， 再由此点求得另一斜率 然后，从 点开始，既不按该点斜率 变化，也不按预估点斜率 变化，而是取两者平