第3章 连续系统的数字仿真PPT讲稿.ppt

《第3章 连续系统的数字仿真PPT讲稿.ppt》由会员分享，可在线阅读，更多相关《第3章 连续系统的数字仿真PPT讲稿.ppt（96页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第第3章章 连续系统的数字仿真连续系统的数字仿真1第1页，共96页，编辑于2022年，星期二现代仿真技术与应用现代仿真技术与应用 章节安排章节安排第一章第一章 概述概述第二章第二章 系统的数学模型系统的数学模型第三章第三章 连续系统的数字仿真连续系统的数字仿真第四章第四章 离散事件系统仿真离散事件系统仿真第六章第六章 分布式交互仿真分布式交互仿真第七章第七章 可视化、多媒体、虚拟现实仿真可视化、多媒体、虚拟现实仿真2第2页，共96页，编辑于2022年，星期二 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.1 3.1 数值积分法数值积分法3.2 3.2

2、 离散相似法离散相似法3.5 3.5 控制系统仿真工具控制系统仿真工具SIMULINKSIMULINK3第3页，共96页，编辑于2022年，星期二 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真4第4页，共96页，编辑于2022年，星期二数值积分法数值积分法离散相似法离散相似法间断非线性系统仿真算法间断非线性系统仿真算法分布参数的仿真算法分布参数的仿真算法工程领域常见的连续系统的仿真算法：工程领域常见的连续系统的仿真算法：现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真5第5页，共96页，编辑于2022年，星期二p数

3、值积分：是数值分析的基本问题，是微分方程初值问题的一种数值积分：是数值分析的基本问题，是微分方程初值问题的一种近似解法，其基本思想是将一阶近似解法，其基本思想是将一阶常微分方程常微分方程（或方程组）转化为（或方程组）转化为差分方程差分方程（即微分方程的离散形式，便于编程实现），从而求其（即微分方程的离散形式，便于编程实现），从而求其数值解；数值解；p系统仿真：在给定初始条件可用数值积分法求解定常连续系统在一系统仿真：在给定初始条件可用数值积分法求解定常连续系统在一定输入作用下的变化过程。定输入作用下的变化过程。现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿

4、真3.1 数值积分法数值积分法6第6页，共96页，编辑于2022年，星期二为为t0,t内内 f 下方的阴影面积下方的阴影面积单步法单步法多步法多步法预估预估校正法校正法显式公式显式公式隐式公式隐式公式 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.1 数值积分法数值积分法 基本原理：基本原理：分类：分类：(3-1)已知某系统的一阶向量微分方程为：已知某系统的一阶向量微分方程为：7第7页，共96页，编辑于2022年，星期二数值积分法基本概念数值积分法基本概念 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真单步法单

5、步法多步法多步法预估预估校正法校正法计算计算yn+1只需用到只需用到yn的值的值计算计算yn+1需用到需用到yn，yn-1，yn-2，yn-k的值的值对不能自行启动的方法，可分两步走，第一步预估，第二步校对不能自行启动的方法，可分两步走，第一步预估，第二步校正；因此精度高，但不适用于实时仿真。正；因此精度高，但不适用于实时仿真。显式公式显式公式隐式公式隐式公式计算计算yn+1时所用数据均已知时所用数据均已知计算计算yn+1时需用到待求量时需用到待求量yn+1，因此不能自启动，因此不能自启动常用的数值积分法：常用的数值积分法：欧拉法、梯形法、四阶龙格欧拉法、梯形法、四阶龙格-库塔法、亚当姆斯法库

6、塔法、亚当姆斯法8第8页，共96页，编辑于2022年，星期二3.1.1 常用的积分法常用的积分法-欧拉法欧拉法一阶微分方程为：一阶微分方程为：对式对式(3-1)在在tk,tk+!区间上积分区间上积分当步长当步长足够小时，有足够小时，有 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真(3-3)(3-4)(3-1)于是可以得到微分方程的数值解为：于是可以得到微分方程的数值解为：(3-6)矩形近似及误差矩形近似及误差9第9页，共96页，编辑于2022年，星期二这种方法的几何意义：这种方法的几何意义：取取k=0,1,2,N,从从t0开始，逐点递推求解开始，逐点递

7、推求解t1时的时的y1，t2时的时的y2,直至直至tn时的时的yn，称之为欧拉递推公式。，称之为欧拉递推公式。现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真(3-6)最简单的数值积分法，只计算一次最简单的数值积分法，只计算一次f(t,y)函数值，计函数值，计算量小，但精度很低，不实用，常用来说明基本概念。算量小，但精度很低，不实用，常用来说明基本概念。当当h很小时，造成的误差是允许的。该算法具有一阶很小时，造成的误差是允许的。该算法具有一阶精度。精度。就是把就是把f(t,y)在区间在区间tk,tk+1内的曲边面积内的曲边面积用矩形面积近似代替。用矩形面积

8、近似代替。3.1.1 常用的积分法常用的积分法-欧拉法欧拉法单步法、显示公式、能自启动。单步法、显示公式、能自启动。10第10页，共96页，编辑于2022年，星期二3.1.1 常用的积分法常用的积分法梯形法（改进欧拉法）梯形法（改进欧拉法）对式对式(3-1)在在tk,tk+!区间上积分区间上积分用梯形面积近似（用梯形面积近似（3-3）的积分项，有：）的积分项，有：现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真(3-3)梯形近似及误差梯形近似及误差(3-7)通过欧拉法计算出通过欧拉法计算出y(tk+1)的近似值：的近似值：代入原微分方程，计算代入原微分方程

9、，计算fk+1的近似值的近似值 ，利用梯形公式求出修正后的，利用梯形公式求出修正后的yk+1，得到得到改进后的欧拉公式为：改进后的欧拉公式为：预估预估校正校正(3-8)11第11页，共96页，编辑于2022年，星期二几何意义：几何意义：现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真计算两次计算两次f(t,y)函数值，计算量增加，但精度有所函数值，计算量增加，但精度有所提高。提高。就是把就是把f(t,y)在区间在区间tk,tk+1内的曲边面积内的曲边面积用梯形面积近似代替。用梯形面积近似代替。改进的欧拉法改进的欧拉法隐式公式、不能自启动、预估隐式公式、不能

10、自启动、预估校正法。校正法。12第12页，共96页，编辑于2022年，星期二3.1.1 常用的积分法常用的积分法梯形法梯形法 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真预估预估-校正法程序框图校正法程序框图13第13页，共96页，编辑于2022年，星期二欧拉法欧拉法梯形法梯形法回顾回顾回顾回顾 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真最简单的数值积分法，只计算一次最简单的数值积分法，只计算一次f(t,y)函数值，函数值，计算量小，但精度很低，不实用，常用来说明基计算量小，但精度很低，不实用，常用来说明基本

11、概念。当本概念。当h很小时，造成的误差是允许的。该很小时，造成的误差是允许的。该算法具有一阶精度。算法具有一阶精度。单步法、显示公式、能自启动。单步法、显示公式、能自启动。计算两次计算两次f(t,y)函数值，计算量增加，但精度有所函数值，计算量增加，但精度有所提高。提高。隐式公式、不能自启动、预估隐式公式、不能自启动、预估校正法。校正法。数值积分法：数值积分法：14第14页，共96页，编辑于2022年，星期二基本思想：基本思想：在积分区间在积分区间 t tn n,t tn+1n+1 内多预估几个点的函数值，然后用其线性组合来代替函内多预估几个点的函数值，然后用其线性组合来代替函数的各阶导数，再

12、与泰勒级数展开式中的各项对比确定其中的系数数的各阶导数，再与泰勒级数展开式中的各项对比确定其中的系数设y(t)为微分方程的解，将其在tk附近以h为变量展成泰勒级数：现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.1.1 常用的积分法常用的积分法龙格库塔法龙格库塔法(3-9)(3-10)r为精度阶次（使用k值的个数），wi为待定系数，由精度确定15第15页，共96页，编辑于2022年，星期二当r=1时，当r=2时，取3.1.1 常用的积分法常用的积分法龙格库塔法龙格库塔法 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字

13、仿真数值解与欧拉递推公式一致。数值解与欧拉递推公式一致。数值解预估数值解预估-校正公式一致。校正公式一致。16第16页，共96页，编辑于2022年，星期二r=3时，取时，取 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.1.1 常用的积分法常用的积分法龙格库塔法龙格库塔法17第17页，共96页，编辑于2022年，星期二r=4时，取 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.1.1 常用的积分法常用的积分法龙格库塔法龙格库塔法18第18页，共96页，编辑于2022年，星期二各阶龙格各阶龙格库塔法的精度库塔

14、法的精度理论上可以构造任意阶数的龙格理论上可以构造任意阶数的龙格库塔法库塔法阶数越高，精度越高，计算量越大阶数越高，精度越高，计算量越大精度的阶数与计算函数值精度的阶数与计算函数值 f 的次数之间并非等量增加的关系的次数之间并非等量增加的关系对于大量的实际问题，四阶方法已可满足精度求对于大量的实际问题，四阶方法已可满足精度求每一步计算每一步计算f的次数的次数234567n 8精度阶数精度阶数234456n-2 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真19第19页，共96页，编辑于2022年，星期二变步长方法变步长方法误差式通常为误差式通常为当 较大

15、时，是相对误差 当 较小时，是绝对误差 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真当估计误差当估计误差en大于最大允许误差大于最大允许误差emax时，步长减半并重新积分再估计误差；时，步长减半并重新积分再估计误差；若步长小于步长下限若步长小于步长下限hmin，则不再减半，以免增加仿真时间和舍入误差，则不再减半，以免增加仿真时间和舍入误差;当估计误差当估计误差en小于最小允许误差小于最小允许误差emin时，步长加倍并重新积分再估计误差；时，步长加倍并重新积分再估计误差；若步长大于步长上限若步长大于步长上限hmin，则不再加倍，以免增加截断误差、减小数值

16、稳，则不再加倍，以免增加截断误差、减小数值稳定性。定性。20第20页，共96页，编辑于2022年，星期二算法误差算法误差 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.1.3 算法误差和稳定性问题算法误差和稳定性问题截断误差截断误差欧拉法欧拉法梯形法梯形法三阶龙格三阶龙格-库塔法库塔法四阶龙格四阶龙格-库塔法库塔法舍入误差舍入误差截去高阶无穷小项引入的误差称为截去高阶无穷小项引入的误差称为r r阶精度阶精度泰勒级数展开式泰勒级数展开式因计算机字长有限，数字不能完全精确表示而产生的误差，与步因计算机字长有限，数字不能完全精确表示而产生的误差，与步长、数

17、字系统、运算次序以及计算长、数字系统、运算次序以及计算f f(t t,y y)子程序的精度等多种子程序的精度等多种因素有关因素有关21第21页，共96页，编辑于2022年，星期二算法稳定性问题算法稳定性问题 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.1.3 算法误差和稳定性问题算法误差和稳定性问题测试方程：测试方程：代入测试方程有代入测试方程有:欧拉法：欧拉法：22第22页，共96页，编辑于2022年，星期二算法稳定性问题算法稳定性问题 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.1.3 算法误差和稳

18、定性问题算法误差和稳定性问题梯形法梯形法:代入测试方程有代入测试方程有:23第23页，共96页，编辑于2022年，星期二各阶龙格各阶龙格库塔法的稳定域库塔法的稳定域将检验方程代入泰勒级数展开式可得稳定条件将检验方程代入泰勒级数展开式可得稳定条件稳定条件与计算步长和系统特征根都有关系，特征根越大，计算步长越小稳定条件与计算步长和系统特征根都有关系，特征根越大，计算步长越小r 1 h所在的实稳定区域所在的实稳定区域11+h(-2,0)21+h+(h)2/2(-2,0)31+h+(h)2/2+(h)3/6(-2.51,0)41+h+(h)2/2+(h)3/6+(h)4/24(-2.78,0)现代仿真

19、技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真24第24页，共96页，编辑于2022年，星期二计算步长的选择计算步长的选择步长过大会增加截断误差，甚至出现不稳定现象；步长过小会增加计算步数，步长过大会增加截断误差，甚至出现不稳定现象；步长过小会增加计算步数，从而增大舍入误差；从而增大舍入误差；步长变化对误差的影响与系统动态响应特性有关，响应步长变化对误差的影响与系统动态响应特性有关，响应越快对步长变化越快对步长变化越敏感。越敏感。现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真有些仿真可考虑采用变步长有些仿真可考虑采用变步长

20、经验方法一：经验方法一：经验方法二：经验方法二：总误差舍入误差截断误差（合理步长）（合理步长）（步长）（步长）25第25页，共96页，编辑于2022年，星期二已知常微分方程组已知常微分方程组 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.1.5 病态系统（刚性（病态系统（刚性（stiff）系统）的仿真方法）系统）的仿真方法定义：定义：(3-34)则有系统矩阵则有系统矩阵 （雅可比矩阵（雅可比矩阵（Jacobi matrix）若若 的特征值均具有负实部，那么当：的特征值均具有负实部，那么当：式（式（3-34）微分方程组称为刚性方程，也称病态方程，其描述

21、的线性或非线性）微分方程组称为刚性方程，也称病态方程，其描述的线性或非线性系统称为病态系统。系统称为病态系统。(3-35)26第26页，共96页，编辑于2022年，星期二l前面介绍的数值积分法进行求解时，仿真步长应该由系统的最小时间常数（相当于是前面介绍的数值积分法进行求解时，仿真步长应该由系统的最小时间常数（相当于是最大特征值实部的倒数）决定，而仿真时间由系统的最大时间常数（相当于最小特征最大特征值实部的倒数）决定，而仿真时间由系统的最大时间常数（相当于最小特征值实部的倒数）决定；值实部的倒数）决定；l由定义可知，病态系统的时间常数（由定义可知，病态系统的时间常数（1/）差别很大；）差别很大

22、；l那么对于病态系统，如果采用前面介绍的数值积分法进行求解那么对于病态系统，如果采用前面介绍的数值积分法进行求解,计算步长只能取得很小，系统计算步长只能取得很小，系统的过渡过程又很长，将导致仿真计算量大、舍入误差累积大、数值解失真等；的过渡过程又很长，将导致仿真计算量大、舍入误差累积大、数值解失真等；l目前已提出的刚性方程数值积分法有吉尔（目前已提出的刚性方程数值积分法有吉尔（Gear）法、特雷尔（）法、特雷尔（Trernor）法、半）法、半隐式龙格隐式龙格库塔法等，其中吉尔法被公认为是十分有效的方法。库塔法等，其中吉尔法被公认为是十分有效的方法。现代仿真技术与应用现代仿真技术与应用 第三章连

23、续系统的数字仿真第三章连续系统的数字仿真3.1.5 病态系统（刚性（病态系统（刚性（stiff）系统）的仿真方法）系统）的仿真方法27第27页，共96页，编辑于2022年，星期二吉尔法为隐式线性多步法，其形式为：吉尔法为隐式线性多步法，其形式为：当当k=1时，为隐式欧拉法：时，为隐式欧拉法：3.1.5 病态系统的仿真方法病态系统的仿真方法吉尔法吉尔法 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真(3-36)为待定系数，为待定系数，1-6阶吉尔法的系数如下表：阶吉尔法的系数如下表：k1110000022/34/3-1/3000036/1118/11-9/

24、112/11000412/2548/25-36/2516/25-3/2500560/137300/137-300/137200/137-75/13712/1370660/147360/147-450/147400/147-225/14772/147-10/14728第28页，共96页，编辑于2022年，星期二吉尔法的吉尔法的StiffStiff稳定稳定一个方法如果在区域一个方法如果在区域R1（Re(h)D）中是绝对稳）中是绝对稳定的，而在区域定的，而在区域R2（D Re(h)，|Im(h)|80*0.89+16*12.2+25*1.82ans=331.900 或者 res=80*0.89+16

25、*12.2+25*1.82res=331.900 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真45第45页，共96页，编辑于2022年，星期二变量与数值显示格式变量与数值显示格式：变量规则变量的名字必须以字母开头（不超过一定的字符），变量的名字必须以字母开头（不超过一定的字符），之后可以使任意的字母、数字和下划线。之后可以使任意的字母、数字和下划线。变量名区分大小写；变量名区分大小写；Matlab 不支持汉字，汉字不能出现在变量名和文件不支持汉字，汉字不能出现在变量名和文件名中名中数值显示常用格式Long(16位)、bank(2个十进制位)、short

26、(默认)、short e(5位加指数)、long e(16位加指数)现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真46第46页，共96页，编辑于2022年，星期二 MATLAB运用与简单运算运用与简单运算转置矩阵运算：矩阵运算：a=1 2;4 5 ;b=ab=1 4 2 5乘方 a=1 2;4 5 ;a2=9 1224 33 a=1 2;4 5 ;a.2=1 416 25 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真47第47页，共96页，编辑于2022年，星期二关系运算逻辑运算a=1:9;b=a4b=0 0 0

27、 0 1 1 1 1 1c=(a4)&(aa=1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 a(2:3,3:4)ans=7 8 11 12 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真 MATLAB运用与简单运算运用与简单运算矩阵运算：矩阵运算：49第49页，共96页，编辑于2022年，星期二 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真 MATLAB运用与简单运算运用与简单运算求解线性方程组：求解线性方程组：10 x+3y+z=142x-10y+3z=-5X+3y+10z=14

28、A=10,3,1;2,-10,3;1,3,10;B=14,-5,14;Root=inv(A)*B%inv为求逆矩阵函数为求逆矩阵函数50第50页，共96页，编辑于2022年，星期二Matlab程序设计程序设计MATLAB运用与简单运算运用与简单运算MATLAB绘图绘图数据处理数据处理SIMULINK基础基础 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.5 matlab与与SIMULINK基础基础控制工具箱控制工具箱51第51页，共96页，编辑于2022年，星期二MATLAB绘图：绘图：二维图形plot(x)plot(x,y);plot(x,y,参数

29、)；plot(x1,y1,参数1，x2,y2,参数2)现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真绘图函数绘图函数plot的几种形式的几种形式x,y可以是可以是向量向量或或矩阵矩阵，参数为，参数为字符串字符串，决，决定定图形颜色图形颜色、线型线型及及数据点的图标。数据点的图标。字符字符颜色颜色字符字符颜色颜色b兰色兰色m紫红色紫红色c青色青色r红色红色g绿色绿色w白色白色k黑色黑色y黄色黄色符号符号线型线型符号符号线型线型-实线实线：点连线点连线-.点画线点画线-虚线虚线控制符控制符标记符标记符控制符控制符标记符标记符控制符控制符标记符标记符.点点d菱

30、形菱形h六角形六角形O(字母字母)圆圈圆圈p五角星五角星V(字母字母)下三角下三角*星号星号s正方形正方形上三角上三角X(字母字母)叉号叉号右三角右三角vartmp)(varout(j)vartmp)y(j)=average;y(j)=average;varout(j)=vartmp;varout(j)=vartmp;endend endend endendclc;clc;clear;clear;close all;close all;len=200;len=200;flen=10;flen=10;snr=1;snr=1;a=zeros(1,len);a=zeros(1,len);for t=

31、1:lenfor t=1:len if t100 if tX%转置后显示自变量一组采样点转置后显示自变量一组采样点 Y%转置后，显示与采样点对应的一组数值解转置后，显示与采样点对应的一组数值解解：解：80第80页，共96页，编辑于2022年，星期二Matlab程序设计程序设计MATLAB运用与简单运算运用与简单运算MATLAB绘图绘图数据处理数据处理SIMULINK基础基础 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.5 matlab与与SIMULINK基础基础控制工具箱控制工具箱81第81页，共96页，编辑于2022年，星期二控制工具箱控制工具箱

32、状态空间系统建模：系统建模：传递函数 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真82第82页，共96页，编辑于2022年，星期二num,den=ss2tf(a,b,c,d)控制工具箱控制工具箱a,b,c,d=tf2ss(num,den)例设系统传递函数为：例设系统传递函数为：试求出状态空间表达式试求出状态空间表达式num=bm,bm-1,b0 num=bm,bm-1,b0 分子多项式系数分子多项式系数den=an,an-1,a0 den=an,an-1,a0 分母多项式系数分母多项式系数 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真

33、第三章连续系统的数字仿真83第83页，共96页，编辑于2022年，星期二 现代仿真技术与应用现代仿真技术与应用 2.1 2.1连续系统的数学模型连续系统的数学模型例设系统的状态方程为：例设系统的状态方程为：求传递函数求传递函数解解:特征多项式为特征多项式为:伴随矩阵为伴随矩阵为:num,den=ss2tf(a,b,c,d)84第84页，共96页，编辑于2022年，星期二Matlab程序设计程序设计MATLAB运用与简单运算运用与简单运算MATLAB绘图绘图数据处理数据处理SIMULINK基础基础 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真3.5 ma

34、tlab与与SIMULINK基础基础控制工具箱控制工具箱Matlab应用应用85第85页，共96页，编辑于2022年，星期二Simulink的基本模块的基本模块 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真信号源模块（信号源模块（sources)输出模块（输出模块（Sinks)数学模块（数学模块（Math)连续系统模块（连续系统模块（Continuous)离散系统模块（离散系统模块（Discrete)函数和表模块（函数和表模块（Functions&Tables)非线性系统模块（非线性系统模块（Nonlinear)信号与系统模块（信号与系统模块（Sign

35、als&Systems)86第86页，共96页，编辑于2022年，星期二Simulink模块操作模块操作 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真a)添加和选取模块添加和选取模块b)模块位置和外形的调整模块位置和外形的调整c)模块名处理模块名处理d)模块的删除和复制模块的删除和复制e)模块属性和参数的设置模块属性和参数的设置f)模块间的连线模块间的连线87第87页，共96页，编辑于2022年，星期二仿真模型参数的设置仿真模型参数的设置 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真Solver options

36、:仿真算法选择仿真算法选择Output options:输出选项选择输出选项选择细化输出；细化输出；产生附加输出；产生附加输出；仅在指定的时刻产生输出仅在指定的时刻产生输出88第88页，共96页，编辑于2022年，星期二 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真例：用例：用Simulink仿真两个正弦信号相乘，即计算仿真两个正弦信号相乘，即计算x(t)=sin(t)*sin(10t)89第89页，共96页，编辑于2022年，星期二 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真例：一个典型线性反馈控制系统结构

37、如下，用例：一个典型线性反馈控制系统结构如下，用simulink仿真软件求出仿真软件求出开环和闭环系统的阶跃响应曲线开环和闭环系统的阶跃响应曲线90第90页，共96页，编辑于2022年，星期二例例：已知某控制系统的动态结构图如下图所示，若试利用MATLAB/Simulink工具分析求阶跃响应曲线 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真91第91页，共96页，编辑于2022年，星期二 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真例：设计一数字电路的例：设计一数字电路的8线线-3线编码器，并用线编码器，并用s

38、imulink仿真，要求输入信号每次仿真，要求输入信号每次只有一个是只有一个是0，其余，其余7个是个是1。输入信号输入信号输出信号输出信号J0 J1 J2 J3 J4 J5 J6 J7Y2 Y1 Y0 0 1 1 1 1 1 1 1 0 0 0 1 0 1 1 1 1 1 1 0 0 1 1 1 0 1 1 1 1 1 0 1 0 1 1 1 0 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 1 0 0 1 1 1 1 1 0 1 1 1 0 1 1 1 1 1 1 1 0 1 1 1 0 1 1 1 1 1 1 1 0 1 1 192第92页，共96页，编辑于2022年，星期二

39、 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真2)设计一数字电路的设计一数字电路的3线线-8线编码器，并用线编码器，并用simulink仿真。仿真。练习：练习：1)一个典型线性反馈控制系统结构如下，用一个典型线性反馈控制系统结构如下，用simulink仿真软件求出闭环系统的阶仿真软件求出闭环系统的阶跃响应曲线跃响应曲线93第93页，共96页，编辑于2022年，星期二 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真例例1：有一个由：有一个由50HZ和和300HZ正弦信号构成的信号，受到均值随机噪声的污染，现对该信正

40、弦信号构成的信号，受到均值随机噪声的污染，现对该信号进行采样，采样频率为号进行采样，采样频率为1000HZ，通过快速傅立叶变换分析其信号频率成分，通过快速傅立叶变换分析其信号频率成分MatlabMatlab应用应用应用应用在信号处理中的应用（在信号处理中的应用（FFT）例例2：用：用FFT变换分析语音信号的频变换分析语音信号的频谱谱94第94页，共96页，编辑于2022年，星期二 现代仿真技术与应用现代仿真技术与应用 第三章连续系统的数字仿真第三章连续系统的数字仿真例：对一副图进行裁剪和消噪处理例：对一副图进行裁剪和消噪处理MatlabMatlab应用应用应用应用在图像处理中的应用在图像处理中的应用95第95页，共96页，编辑于2022年，星期二结束结束 计算机导论计算机导论 第一章概述第一章概述96第96页，共96页，编辑于2022年，星期二