离散化原理及要求和常用的几种数值积分法.ppt

第三章 连续系统数值积分仿真方法学5/12/20241.3.1 离散化原理及要求离散化原理及要求 n（1）离散化原理）离散化原理n（2）离散化建模方法的要求）离散化建模方法的要求2.（1）离散化原理a.“数字数字”计算，引入舍入误差；计算，引入舍入误差；b.按指令一步一步进行，必须将时间按指令一步一步进行，必须将时间离散化。离散化。在数字计算机上对连续系统进行仿真时，在数字计算机上对连续系统进行仿真时，首先遇到的问题是：数字计算机的数值及时间首先遇到的问题是：数字计算机的数值及时间均具有离散性，而被仿真系统的数值及时间均均具有离散性，而被仿真系统的数值及时间均具有连续性，后者如何用前者来实现。具有连续性，后者如何用前者来实现。3.n连续系统仿真：连续系统仿真：n从时间、数值两个方面对原系统进行离散从时间、数值两个方面对原系统进行离散化，并选择合适的数值计算方法来近似积化，并选择合适的数值计算方法来近似积分运算，由此得到离散模型来近似原连续分运算，由此得到离散模型来近似原连续模型。模型。4.相似原理相似原理 n设系统模型为：设系统模型为：离散后的输入变量：离散后的输入变量：系统变量：系统变量：如果：如果：即：即：仿真时间间隔为：仿真时间间隔为：h两模型等价。两模型等价。5.u(t)图图2.1 2.1 相相 似似 原理原理原连续模型原连续模型仿真模型仿真模型h y(t)+6.（2）离散化建模方法的要求1.1.稳定性稳定性2.2.准确性：准确性：n最基本的准则是：最基本的准则是：n绝对误差准则：绝对误差准则：n 相对误差准则：相对误差准则：系统时间间隔系统时间间隔hk=tk+1-tk计算每一步间隔计算每一步间隔Tk若若hk Tk，实时仿真实时仿真若若Tkhk，离线仿真离线仿真3.快速性快速性7.明确几个概念明确几个概念8.差分方程已知表示某系统一阶向量微分方程及初已知表示某系统一阶向量微分方程及初值为：值为：对上式两边积分，则对上式两边积分，则9.在在 时的连续解为：时的连续解为：令令则则或表示为或表示为 10.数值解法数值解法相邻两个离散点的间距相邻两个离散点的间距常用的基本方法有三类：常用的基本方法有三类：单步法、多步法、预估校正法。单步法、多步法、预估校正法。并可分为显式公式和隐式公式。并可分为显式公式和隐式公式。就是寻求初值问题式就是寻求初值问题式(1)的解在一系列离散点的解在一系列离散点的近似解的近似解（即数值解）。（即数值解）。称为计算步长或步距称为计算步长或步距11.单步法与多步法单步法与多步法n单步法单步法n只由前一时刻的数值只由前一时刻的数值 ym就可求得后一时刻就可求得后一时刻的数值的数值ym+1n能自动启动能自动启动n多步法多步法n计算计算ym+1需要用到需要用到 tm,tm-1,tm-2,时刻时刻y的数的数据据n不能自动启动不能自动启动12.显式与隐式显式与隐式n显式显式n计算计算 ym+1时所用数值均已计算出来时所用数值均已计算出来n隐式隐式n计算中隐含有未知量计算中隐含有未知量n预估校正法预估校正法 使用隐式公式时，需用另一显式公式估计使用隐式公式时，需用另一显式公式估计一个初值，然后再用隐式公式进行迭代运算。一个初值，然后再用隐式公式进行迭代运算。13.截断误差n假设前一步得到的结果假设前一步得到的结果ym是准确的是准确的,则用泰勒级则用泰勒级数求得数求得tm+1处的精确解为处的精确解为n若从以上精确解中取前两项之和来近似若从以上精确解中取前两项之和来近似计算计算ym+1,由这种方法单独引进的附加误差由这种方法单独引进的附加误差通常称作通常称作局部截断误差局部截断误差.14.舍入误差n舍入误差与舍入误差与h成反比成反比,若计算步长小，计若计算步长小，计算次数就多算次数就多,则舍入误差就大。则舍入误差就大。15.3.2 常用的几种数值积分法n建立系统数学模型的目的是研究系统的建立系统数学模型的目的是研究系统的运动规律运动规律16.（一）单步法17.（1）欧拉法（一阶龙格库塔法）nTaylor级数展开n矩形近似解法n切线近似18.（a)Taylor展开假定假定为其解析解为其解析解将将y(t)展开成展开成Taylor级数级数从而从而将上式写成差分方程将上式写成差分方程19.（b）矩形近似解法在区间在区间tn,tn+1上积分，得上积分，得f误差误差近似矩形近似矩形f(t)0tntn+1t于是于是20.（c)切线近似y(t)在在tn处得切线方程为处得切线方程为则得则得y(t)y0tntn+1tynyn+1(t0,y0)(t1,y1)t21.例例1 1 设系统方程为：设系统方程为：试用欧拉法求其数值解（取步长试用欧拉法求其数值解（取步长h=0.1，0t1）解：原方程为：解：原方程为：递推公式为：递推公式为：22.已知方程的解析解为精确解与数值解比较精确解与数值解比较t00.10.20.30.40.51.0精确解精确解y(t)10.9091 0.83330.76920.66670.6250.5数值解数值解yn10.90.8190.75190.65940.6470.4628误差在误差在102数量级数量级23.（2）改进的欧拉法（梯形法）n又称二阶龙格库塔法又称二阶龙格库塔法f误差误差f(t)0tntn+1t直边梯形的面积直边梯形的面积当当h比较小时，以直边梯形代替曲边梯形的面积，可得比较小时，以直边梯形代替曲边梯形的面积，可得其差分方程其差分方程或或曲边梯形的面积曲边梯形的面积24.1.用欧拉法求出初值，算出用欧拉法求出初值，算出的近似值的近似值2.计算导函数计算导函数近似值近似值3.然后用梯形法求出修正后的然后用梯形法求出修正后的25.1.用欧拉法预估一个初值用欧拉法预估一个初值2.用下式求出用下式求出3.再用再用求求如此反复下去直到如此反复下去直到迭代运算：迭代运算：26.预估公式预估公式校正公式校正公式预估校正法预估校正法27.(3)龙格库塔法n基本思想基本思想n间接利用泰勒展开式，即用几个点上的间接利用泰勒展开式，即用几个点上的y(t)的一阶导函数值的线性组合来近似代的一阶导函数值的线性组合来近似代替替y(t)在某一点的各阶导数，然后用泰勒在某一点的各阶导数，然后用泰勒级数展开式确定线性组合中各加权系数级数展开式确定线性组合中各加权系数28.（a)龙格库塔数值积分公式推导一阶微分方程：一阶微分方程：假定假定y(t)是上式的解析解。将是上式的解析解。将y(t)展开成泰勒级数展开成泰勒级数其中：其中：29.将将y(t+h)写成线性组合形式写成线性组合形式其中其中r称为阶数，称为阶数，bi待定系数，待定系数，ki由下式决定由下式决定且定义且定义C1=030.r=1，此时此时c1=0，a1=0，k1=f(t,y)，则则取取b1=1，即得一阶龙格库塔法，即得一阶龙格库塔法r=2将将在点在点(t,y)展开泰勒级数展开泰勒级数31.令令，得，得所以所以改进欧拉法改进欧拉法32.r=4其中：其中：四阶龙格库塔法四阶龙格库塔法33.为方便，将上式具体列为：为方便，将上式具体列为：其中：其中：34.龙格库塔法的特点n在计算在计算yn+1时只用到时只用到yn，而不直接用，而不直接用yn-1,yn-2等等项项;n步长步长h在整个计算中并不要求固定；在整个计算中并不要求固定；n精度取决于步长精度取决于步长h的大小及方法的阶次的大小及方法的阶次n一阶龙格库塔公式一阶龙格库塔公式欧拉公式欧拉公式35.优点n编制程序容易编制程序容易n改变步长方便改变步长方便n稳定性较好稳定性较好n是一种自启动的数值积分法是一种自启动的数值积分法36.(4)单步法的特点n 需要存储的数据量少需要存储的数据量少n 可自启动可自启动n 容易实现变步长运算容易实现变步长运算37.例：已知系统方程例：已知系统方程取步长取步长计算计算时的时的y值值38.已知微分方程已知微分方程 ，分别用欧拉法、，分别用欧拉法、梯形法和四阶龙格库塔法写出前两步的差分梯形法和四阶龙格库塔法写出前两步的差分方程的解方程的解(t0=0，y0=0,步长步长h=0.1)课堂测验：课堂测验：39.