第二章数值积分法.ppt

第二章 连续系统仿真的数值积分法 2.1数值积分法 2.2欧拉法（Euler） 2.3龙格-库塔法（Runge-Kutta） 2.4亚当姆斯法（Adams） 2.5吉尔法（Gear） 2.6方法与步长的选择 2.1数值积分法 一、在仿真中研究数值积分法的意义 数值积分是数值分析的一个基本问题，用于求一个积分的具体数值。 在数学上，用数值积分法求微分方程的初值问题。 对于一个连续系统，其模型可用微分方程（组）来表示，而对其进行仿真即是从初始条件出发，求其过渡过程。这在本质上是微分方程的初值问题。 （1） 二、数值积分法的基本原理 对微分方程（1）两边求积分可得下式： 即： 当t=tn+1时，上式为： 将积分项分为两项： 箭头所指的部分结构形式相同，只是积分上限不同。 故有： 若用近似值yn+1、yn、Qn代替上式中的各项，则可写出 yn+1 = yn+Qn t=t0即 n=0时， y1=y0+Q0 其中，y0已知，若用数值积分法求出Q0即可得到y1。依此类推，可求出t1 、t2、… 、tn时刻的近似值 y1、y2、 … 、 yn 可见，这是一种步进式的近似解法，故存在精度和稳定性问题。 用不同的方法求积分值Qn就形成了不同的数值积分法，每种方法的精度、稳定性和计算速度各不相同。 我们希望的数值积分法：精度高，速度快，稳定性好。 2.2欧拉方法（Euler） 一、 欧拉法 欧拉法是最简单的数值积分法，常用来说明数值积分法的几何意义和基本概念。 曲线下的面积 用矩形面积代替曲线下的面积 若区间[ tn，tn+1 ]足够小 可认为 f（，y）= f（tn，yn） 则有 yn+1 = yn + h f （tn，yn） 用矩形面积代替曲线下的面积，有较大误差。 欧拉公式 h为步长 h = tn +1 - t n 一般不单独使用 多与隐式公式配合使用，用来求初值。 计算较简单。 t1=0.1 y(t1)=0.9091 y1=1+0.1(-1²)=0.9 t2=0.2 y(t2)=0.8333 y2=0.9+0.1(-0.9²)=0.819 t3=0.3 y(t3)=0.7692 y3=0.819+0.1(-0.819²)=0.7519 t4=0.4 y(t4)=0.7143 y4=0.6954 t5=0.5 y(t5)=0.6667 y5=0.6470 t6=0.6 y(t6)=0.6250 y6=0.6052 t7=0.7 y(t7)=0.5882 y7=0.5685 t8=0.8 y(t8)=0.5556 y8=0.5362 t9=0.9 y(t9)=0.5263 y9=0.5075 t10=1 y(t10)=0.5 y10=0.4817 例一 . 已知 y’ + y² = 0，y(0) = 1，步长h=0.1 用欧拉法求数值解，并与解析解比较. 解: 解析解为 y(t)= 1 / (1+t)， 数值解为 yn+1= yn+hf(tn,yn)，y’ = -y2 误差的数量级在10¯², 即 h² 二.改进欧拉法 用梯形面积代替曲线下的面积,可提高精度。 yn+1= yn + h/2[ f(tn,yn)+ f(tn+1,yn+1)] 用梯形面积代替曲线下的面积 由于梯形公式(2.2)的右端有待求量，不能自启动，故需要与欧拉法结合使用，这就是改进欧拉法，实际上这是一种预估 —校正法。 预估 yn+1 = yn + hf(tn,yn) 校正 yn+1 = yn + h/2[f(tn,yn)+f(tn+1,yn+1)] 改进欧拉法 yn+1= yn + h/2[ f(tn,yn)+ f(tn+1,yn+1)] 解析解 y=1/(1+t) 数值解 t1=0.1 y(t1)=0.9091 y1=0.9095 0.410-3 t2=0.2 y(