连续系统的数字仿真课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《连续系统的数字仿真课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 连续 系统 数字 仿真 课件
- 资源描述:
-
1、 计算机仿真离散相似计算机仿真离散相似连续时间数学模型连续时间数学模型 离散离散(采样)(采样)相似(信号重构保持)相似(信号重构保持)离散时间数学模型(计算机仿真模型)离散时间数学模型(计算机仿真模型)本章内容本章内容3.1 微分方程的数值积分数值解微分方程的数值积分数值解 3.1.1 数值积分法数值积分法 3.1.2 数值积分法的分类数值积分法的分类 3.1.3 单步法单步法 3.1.4 多步数值积分法多步数值积分法 3.1.5 关于数值积分的误差关于数值积分的误差3.2 从传递函数到从传递函数到Z函数的转换函数的转换 3.2.1 Z域的离散相似模型域的离散相似模型 3.2.2 从传递函数
2、到从传递函数到Z函数的算子替换法函数的算子替换法3.3 从连续时间状态方程到离散时间状态方程的转换从连续时间状态方程到离散时间状态方程的转换3.4 连续时间结构框图典型环节的离散相似连续时间结构框图典型环节的离散相似 3.1 微分方程的数值积分数值解微分方程的数值积分数值解返回返回3.1.1 数值积分法数值积分法对于初值已知的常微分方程对于初值已知的常微分方程 )(3.1),)(,()(00ytybtatytfty 求其数值解就是求解曲线求其数值解就是求解曲线y(t)在一系列离散时间点在一系列离散时间点 tn的的 近似值近似值 yn=y(tn)(n=0,1,2,m)a=t0 tn tn+1 t
3、m=b相邻两个节点的间距:步长相邻两个节点的间距:步长 h=tn+1-tn。设微分方程在区间设微分方程在区间a,b上连续变化。在已知初值条件下,将区间上连续变化。在已知初值条件下,将区间a,b分分成若干个小的区间,时间间隔为成若干个小的区间,时间间隔为h,在一个小的区间,在一个小的区间t=tn,tn+1内积分,有内积分,有(3.2),(11 nnttnndtytfyy )(3.1),)(,()(00ytybtatytfty (3.3)()(!1 !21)()(1)(2 rrrhhtyrhyhytyhty 图图3.1 数值积分计算示意图数值积分计算示意图),(ytfy t其中,其中,h为积分步长
4、,为积分步长,为局部截断误差。为局部截断误差。)(1 rho 若截断到若截断到r 阶,阶,y 近似的数值解为近似的数值解为(3.4)(!1 !21)()()(2rrhtyrhyhytyhty 截断误差截断误差 为为)(1 rho001(1)1(2)2(3)3111()(1)!(2)!(3)!rrrrrrrt ty yo hyhyhyhrrr 3.1.2 数值积分法分类数值积分法分类 单步积分法和多步积分法单步积分法和多步积分法数值积分法在求解(数值积分法在求解(n+1)时刻的数值解,只涉及到前一个时刻)时刻的数值解,只涉及到前一个时刻n的值的值yn,还,还是涉及过去多个时刻是涉及过去多个时刻
5、tn-1,tn-2,的值,分为的值,分为单步法单步法和和多步法多步法:单步法:单步法:仅用仅用yn的数值即可求出的数值即可求出yn+1的数值;的数值;多步法多步法:求解:求解tn+1时刻的数值解时刻的数值解yn+1,不仅需要知道,不仅需要知道tn时刻的时刻的yn,还要用,还要用 到过去的到过去的tn-1,tn-2,时刻的值。时刻的值。显式积分方法和隐式积分方法:显式积分方法和隐式积分方法:根据数值积分法公式的右端是否包含当前时刻的值,分为显式积分方法根据数值积分法公式的右端是否包含当前时刻的值,分为显式积分方法和隐式积分方法。和隐式积分方法。例:例:单步显式积分方法的公式为单步显式积分方法的公
6、式为 (3.6),(1hythyynnnn 式中右端只与式中右端只与tn+1时刻之前的数值有关,且不包含时刻之前的数值有关,且不包含yn+1项。项。单步隐式积分方法的一般表达式为单步隐式积分方法的一般表达式为(3.7),(11hyythyynnnnn 该式右端只与该式右端只与tn+1时刻之前的数值有关,但是包含时刻之前的数值有关,但是包含yn+1的项。的项。3.1.3 单步法单步法 欧拉法欧拉法 梯形法梯形法Runge-Kutta法法 1.欧拉法欧拉法 (Euler)(3.1)(,()(00ytytytfty 已知:一阶微分方程式(已知:一阶微分方程式(3.1)求解的递推公式(求解的递推公式(
7、3.2)中)中),(ytfy tnt1 nt(3.2),(11 nnttnndtytfyy在区间在区间 tn,tn+1 上的积分如图上的积分如图3.2中黑色斜线部分所示。对于足够小的积分中黑色斜线部分所示。对于足够小的积分间隔间隔tn,tn+1,其中每一个时刻的,其中每一个时刻的f(t,y)都可以近似地为常数都可以近似地为常数 f(tn,yn),即,即用矩形面积近似代替该区间曲线面积(红色部分)。用矩形面积近似代替该区间曲线面积(红色部分)。欧拉法的几何意义欧拉法的几何意义图图3.2 欧拉法的几何表示欧拉法的几何表示),(nnytf),(11 nnytf(3.9)()(1htyyhtyynnn
8、n 局部截断误差为局部截断误差为o(h2),即,即00 !31!21)(322yyttnhyhyhR 将将y(t+h)展开成泰勒级数如式(展开成泰勒级数如式(3.4),若在右边只取两项,有),若在右边只取两项,有),(ytfy tnt1 nt截断误差截断误差 应用泰勒级数展开推导欧拉法公式应用泰勒级数展开推导欧拉法公式已知初始已知初始y0,可以递推得出各时间点的,可以递推得出各时间点的yn近似数值解。近似数值解。),2 ,1 ,0(),(),(),()()()()(111120000010000nythfyyythfyyythfyhtytyhtyyytynnnnyytt(3.10)),(nny
9、tf),(11 nnytf 2.梯形法梯形法 梯形法的基本思路是在应用梯形面积近似在图梯形法的基本思路是在应用梯形面积近似在图3.2中积分区间的曲线中积分区间的曲线下积分面积,如图下积分面积,如图3.3所示。所示。在区间在区间 tn,tn+1 的曲线积分面积为的曲线积分面积为)()(),(11nntttytydtytfSnn 曲曲在区间在区间 tn,tn+1 的梯形面积为的梯形面积为 ),(),(2111 nnnnytfytfhS梯梯当积分步长当积分步长h0时,梯形面积近似等于曲边梯形面积,即时,梯形面积近似等于曲边梯形面积,即 (3.12),(),(21)()(111 nnnnnnytfyt
10、fhtyty式(式(3.12)的右端有)的右端有yn+1项,为隐式数值积分公式。隐式数值积分在每次项,为隐式数值积分公式。隐式数值积分在每次递推计算中需要借助显式算法(如欧拉法)对递推计算中需要借助显式算法(如欧拉法)对y pn+1预估后再求解。预估后再求解。),(ytfy tnt1 nt),(nnytf),(11 nnytf截断误差截断误差),(1nnnpnytfhyy 预估预估校正校正),(),(2111pnnnnnnytfytfhyy 注注1:预估可采用多种显式方法,上例中应用欧拉法进行预估。:预估可采用多种显式方法,上例中应用欧拉法进行预估。注注2:对应于每个:对应于每个yn,都要经过
11、预估和校正两部运算。,都要经过预估和校正两部运算。局部截断误差为局部截断误差为o(h3),即,即004)4(33!41 !31)(yyttnhyhyhR 例如采用欧拉法进行预估,则梯形法的公式为例如采用欧拉法进行预估,则梯形法的公式为:例例3.1 已知一阶微分方程为已知一阶微分方程为 ,初始条件为初始条件为 。积分步长为积分步长为 ,在区间,在区间0,0.5 分别应用:分别应用:1)欧拉法;)欧拉法;2)梯形法,)梯形法,计算计算y的近似值。的近似值。42tyy5.0)0(y1.0 h解:(解:(1)欧拉法:已知)欧拉法:已知 ,根据欧拉公式,根据欧拉公式 5.0)0(y),(1nnnnyth
12、fyy 计算得计算得 5904.1)42(4380.1)42(2600.1)42(0500.1)42(8000.0)42(5000.0)0(444533342223111200010 tyhyytyhyytyhyytyhyytyhyyyy(2)梯形法)梯形法 5000.00 y8000.0)42(0001 tyhyyp 7750.042422110001 tytyhyyp0300.1)42(1112 tyhyyp 0095.142422221112 tytyhyyp2276.1)42(2223 tyhyyp 2108.142422332223 tytyhyyp3986.1)42(3334 ty
13、hyyp 3848.142422443334 tytyhyyp5478.1)42(4445 tyhyyp 5366.142422554445 tytyhyyp),(1nnnpnytfhyy ),(),(2111pnnnnnnytfytfhyy 预估:预估:校正:校正:从微分方程的解析解:从微分方程的解析解:4721452 teyt可以计算相应的理论值,分别计算欧拉法和梯形法数值解与理论值的误差,可以计算相应的理论值,分别计算欧拉法和梯形法数值解与理论值的误差,如表如表3.1。nt nn 阶阶Runge-Kutta法的通式法的通式(3.15)11 riiinnkbhyy3.龙格龙格-库塔法库塔法
14、 其中,其中,r为阶数;为阶数;bi为待定系数;系数为待定系数;系数ki 表达式为表达式为(3.16),3,2,1 )(,(11rikahtyhctfkijjjnii aj、ci为待定系数。为待定系数。n一阶龙格一阶龙格-库塔法库塔法),(11nnnnythfbyy 与欧拉公式比较,一阶与欧拉公式比较,一阶Runge-Kutta法的系数法的系数 a1=0,b1=1,c1=0。n二阶龙格二阶龙格-库塔法库塔法 (r=2)(21211kkhyynn ),(),(),(),(111121hkytfytfytfkytfknnnnnnnn 1,0 n11a5.021 bb12 c与梯形法比较,二阶与梯形
15、法比较,二阶Runge-Kutta法系数法系数 n高阶龙格高阶龙格-库塔法库塔法 )2,()21,21(),()4(62131213211hKhKyhtfKhKyhtfKytfKKKKhyynnnnnnnnn四阶龙格四阶龙格-库塔法(库塔法(r4)),()21,2()21,21(),()22(6342312143211hKyhtfKhKyhtfKhKyhtfKytfKKKKKhyynnnnnnnnnn 局部截断误差为局部截断误差为o(h4)局部截断误差为局部截断误差为o(h5)n三阶龙格三阶龙格-库塔法(库塔法(r3)例例3.2 用三阶用三阶Runge-Kutta法求例法求例3.1中微分方程的
16、数值解中微分方程的数值解 解:解:三阶三阶Runge-Kutta法的公式为法的公式为 1123(4)6nnhyykkk hkhkyhtfkhkyhtfkytfknnnnnn )2,()21,21(),(213121 1,0 n 以此类推,可求得以此类推,可求得y3=1.2141,y4=1.3885,y5=1.5403.00 t5000.00 y1.01 t10123(4)0.77676hyykkk hkhkyhtfkhkyhtfkytfk 6000.2)2,(7500.2),(3),(2100310020012.02 t21123(4)1.01226hyykkk hkhkyhtfkhkyhtf
17、kytfk 2192.2)2,(3420.2),(5467.2),(211131112111 三阶三阶Runge-Kutta法计算结果与微分方程精确解间的比较如表法计算结果与微分方程精确解间的比较如表3.2所示所示 nt 问题:问题:单步单步?多步?多步?步长总是不变吗?步长总是不变吗?3.1.4 多步数值积分法多步数值积分法 Adams显式积分法显式积分法 Adams隐式积分法隐式积分法 多步法数值积分法有多步法数值积分法有Adams法,法,Simpson法,法,Milne方法等,方法等,我们以我们以Adams法为例介绍多步积分法的思路。法为例介绍多步积分法的思路。Adams法也分为显式积分
18、和隐式积分,下面介绍两种方法的概念。法也分为显式积分和隐式积分,下面介绍两种方法的概念。Adams显式积分法显式积分法微分方程微分方程 在区间在区间tn tn+1上进行积分,递推式为上进行积分,递推式为 )(,()(tytfty 1 1),(nnttnndtytfyy假设在假设在k+1个节点处的导数值分别为个节点处的导数值分别为fn,fn-1,fn-2,fn-k,,基于插值原理,基于插值原理可以构造出一个线性多项式来逼近函数可以构造出一个线性多项式来逼近函数 ,分项积分后可得,分项积分后可得)(,(tytf)(1101knknnnnfffhyy 其中,其中,i为系数,上式的右端没有为系数,上式
19、的右端没有yn+1项,为显式。项,为显式。三阶三阶Adams显式显式(3.19)51623(12211 nnnnnfffhyy四阶四阶Adams显式显式(3.20)9375955(243211 nnnnnnffffhyy Adams 隐式积分法隐式积分法Adams隐式积分公式为隐式积分公式为(3.21)(11110111 knknnnnnffffhyy 三阶三阶Adams隐式隐式)85(12111 nnnnnfffhyy四阶四阶Adams隐式隐式)5199(242111 nnnnnnffffhyy 三阶局部截断误差为三阶局部截断误差为o(h4)。)51623(12211 nnnnpfffhyy
20、n显式预估显式预估隐式校正隐式校正)85(12111 nnpnnnfffhyy与梯形法类似,应用与梯形法类似,应用Adams隐式积分也要选用适当的显式,隐式积分也要选用适当的显式,进行预估以后对进行预估以后对Adams隐式求数值解,隐式求数值解,例如选用三阶例如选用三阶Adams显式积分法进行预估,可得显式积分法进行预估,可得 方方 法法显式欧拉法显式欧拉法01000显式梯形法显式梯形法03/2-1/2003阶显式阶显式Adams法法023/12-16/125/1204阶显式阶显式Adams法法055/24-59/2437/24-9/24隐式梯形法隐式梯形法1/21/20003阶隐式阶隐式Ad
21、ams法法5/128/12-1/12004阶隐式阶隐式Adams法法9/2419/24-5/241/24010123(n为阶数)为阶数)(3.21)应用公式(应用公式(3.21)可以综合表示常用的数值积分法,系数如表所示。)可以综合表示常用的数值积分法,系数如表所示。)(11110111 knknnnnnffffhyy 3.1.5 关于数值积分的误差关于数值积分的误差 局部截断误差与舍入误差局部截断误差与舍入误差 应用数值积分法所得到的离散数值解是原连续时间系统精确解的近似,应用数值积分法所得到的离散数值解是原连续时间系统精确解的近似,包括两种误差:包括两种误差:截断误差、舍入误差截断误差、舍
22、入误差。n局部截断误差局部截断误差:由于数值积分方法截:由于数值积分方法截断断r阶以后的项所导致的误差。阶以后的项所导致的误差。n舍入误差舍入误差:由于计算机的位数有限而:由于计算机的位数有限而导致;在每步计算均会发生,随计算导致;在每步计算均会发生,随计算次数增加使舍入误差积累值增加。次数增加使舍入误差积累值增加。图图3.4 数值积分法误差与积分步长的关系数值积分法误差与积分步长的关系 例例3.3 例例3.1中的一阶微分方程中的一阶微分方程 ,初始条件为初始条件为 积分步长为积分步长为 ,在区间,在区间0,0.5 分别应用:分别应用:1)欧拉法;)欧拉法;2)梯形法;)梯形法;3)三阶)三阶
23、Runge-Kutta法计算法计算y的近似值,比较不同阶次数值积分法的局部的近似值,比较不同阶次数值积分法的局部 截断误差。截断误差。42 tyy5.0)0(y1.0 htn n固定步长:在整个仿真计算过程中,积分步长固定步长:在整个仿真计算过程中,积分步长h始终保持不变。始终保持不变。n变步长:在仿真积分计算每一步,首先估计计算误差;判断误差变步长:在仿真积分计算每一步,首先估计计算误差;判断误差是否在允许误差范围内:是否在允许误差范围内:若是,则该步计算有效;若是,则该步计算有效;否则计算无效,设法减小步长,重新计算。否则计算无效,设法减小步长,重新计算。n在变步长积分中,必须解决误差估计
24、方法和步长调整策略问题。在变步长积分中,必须解决误差估计方法和步长调整策略问题。积分步长的选择与控制积分步长的选择与控制 在数值积分法中,积分步长的选择不仅受计算工作量和误差的影响,在数值积分法中,积分步长的选择不仅受计算工作量和误差的影响,还与系统的动态响应特性相关。还与系统的动态响应特性相关。在数值积分中,系统各个组成部分之间可能存在着相差很大的部分:在数值积分中,系统各个组成部分之间可能存在着相差很大的部分:l 若按变化平缓部分选择的积分步长若按变化平缓部分选择的积分步长h,对于变化剧烈的部分步长太大,对于变化剧烈的部分步长太大,不能反映系统的(特别是开始时的)动态响应;不能反映系统的(
展开阅读全文