龙格库塔方法基本原理.pptxVIP

会计学;由于实际运算只能完成有限项或有限步运算，因此要将有些需用极限或无穷过程进行的运算有限化，对无穷过程进行截断，这样产生的误差成为截断误差。 根据实际情况建立的数学模型往往难以求解。通常需要通过近似替代，将所求解的数学模型简化为易求解的数值计算问题后再进行求解。数学模型的理论解与数值计算问题的精确解之间的误差称为截断误差。这是计算方法本身带来的误差，所以也成为方法误差。;2021/6/25;2021/6/25;2021/6/25;2021/6/25;2021/6/25;2021/6/25;一般龙格－库塔方法的形式为;10;11;12;二阶龙格—库塔法 在[xi, xi+1]上取两点xi和xi+a2= xi +a2h,以该两点处的斜率值K1和K2的加权平均(或称为线性组合)来求取平均斜率k*的近似值K，即 ;2021/6/25;2021/6/25;2021/6/25;2021/6/25;若取 ，就是另一种形式的二阶龙格 - 库塔公式。 ;2021/6/25;;2021/6/25;参数的选择不唯一，从而构成一类不同的三阶R-K公式，下面给出一种常用的三阶R-K公式，形似simpson公式：;2021/6/25;2021/6/25;2021/6/25;26;27;28;29;; （1）用欧拉法计算 根据欧拉公式，将函数表达式及其初始值代入后，可得该系统仿真第一步的值：;;;二阶龙格-库塔公式;（3）用二阶龙格—库塔法计算 根据公式先计算出两个系数，再计算仿真第一步的值： ;;;（4）用四阶龙格—库塔公式计算 根据公式先计算出4个系数，再计算仿真第一步的值： ;;则系统仿真第一步的值为：;从上述结果可以看出: 对于同一个系统进行仿真计算时，其值的精度是随着数值积分公式的变化而改变的，其中欧拉法计算精度最低，其次为梯形法和二阶龙格—库塔法，四阶龙格—库塔法计算精度最高。 ;龙格-库塔法的误差估计;龙格-库塔法的误差估计;龙格-库塔法的误差估计;龙格-库塔法的步长控制;仿真模型的运行速度与实际系统运行速度一致，称为 实时仿真。一般方法难以满足实时仿真的需要：所得 模型的执行速度较慢；机理也不符合实时仿真的需要 。 假设对一般形式的系统进行仿真： 以RK-2为例进行分析，其公式为 ; 假定在h/2的时间内计算机刚好计算一次右端函数f,则计算分为两步： 1 在tk时刻利用当前的un、yn计算K1； 2 在tn+h/2时刻计算K2,此时un+1无法得到，但实时仿真除了要满足执行速度的要求外，还要求实时接收外部输入，并实时得到输出。; 此种情况下，解决的方法有两个：对un+1进行预报（增大仿真误差）或仿真延迟半个计算步距。后者的计算流程如下 可见，后种方法的输出也会延迟半个计算步距，为了克服这个缺陷，人们提出了如下形式的实时二阶RK法; 实时RK-2公式：;;2021/6/25;2021/6/25;2021/6/25;2021/6/25