《连续系统数字仿真PPT课件.ppt》由会员分享，可在线阅读，更多相关《连续系统数字仿真PPT课件.ppt（105页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、关于连续系统的数字仿真第一张，PPT共一百零五页，创作于2022年6月v第3章 连续系统的数字仿真-离散相似法v 控制系统的数字仿真就是控制系统的数学模型在数字计算机上求解的过程。从前一章可以看到，控制系统的动态模型一般是用常微分方程、一阶常微分方程组(状态方程)或传递函数来描述的。因此，要在数字计算机上进行这类系统的仿真，必须先将连续模型变换为离散化的模型。这一章主要介绍连续模型离散化的方法及其差分方程的求取，最后介绍仿真程序的设计。第二张，PPT共一百零五页，创作于2022年6月v31 连续系统的离散化v 连续系统的响应是随时间连续变化的，但是连续系统的解析解无法用数字计算机求出，只能求出
2、其数值解。也就是说，只能得到连续响应曲线上的有限个点。为此，必须把连续系统离散化，得到差分方程，再用数字计算机求解。这就是把微分运算转化为算术运算的过程。下面先讨论线性定常系统的离散过程。v 设一线性定常系统为v v式中：X为n1维状态向量；U为r1维输入向量；A为nn维状态矩阵；B为nr维输入矩阵；Y为m1维输出向量；C为mn维输出矩阵；D为mr维传递矩阵。v此系统的方框，如图31所示。第三张，PPT共一百零五页，创作于2022年6月v此系统的方框，如图31所示。第四张，PPT共一百零五页，创作于2022年6月第五张，PPT共一百零五页，创作于2022年6月第六张，PPT共一百零五页，创作于
3、2022年6月v 由图33所示的离散结构即可导出连续系统离散后的离散数学模型，这个模型通常为差分方程的形式。由于这个过程使得离散系统与连续系统相似，因此可以认为差分方程式的解序列就是连续系统运动的采样值。这就是离散相似法。v 严格地讲，系统输出处的采样开关后面也应加上再现环节，才能与原系统相似。但是在仿真时，用计算机也只能得到离散序列的解，所以输出处的再现环节加与不加对于离散解序列都是一样的。实际上输出处的采样开关加与不加也无所谓，只要认为离散后的系统与原系统在采样点上的输出值近似相等就行了。第七张，PPT共一百零五页，创作于2022年6月v由香农定理可知，为使采样后的信号能无失真地再现，采样
4、频率比信号最高频率要高两倍以上(或采样周期Ts0.5Tmin)，且加入的保持器应有如图34所示的频率特性。第八张，PPT共一百零五页，创作于2022年6月v311 零阶保持器v 这是一种最常见的保持器，它的脉冲传递函数和频率特性如图35和图36所示。第九张，PPT共一百零五页，创作于2022年6月第十张，PPT共一百零五页，创作于2022年6月第十一张，PPT共一百零五页，创作于2022年6月第十二张，PPT共一百零五页，创作于2022年6月v312 一阶保持器 v 一阶保持器也称外推器，它的脉冲传递函数和频率特性如图38和图39所示。第十三张，PPT共一百零五页，创作于2022年6月第十四张
5、，PPT共一百零五页，创作于2022年6月第十五张，PPT共一百零五页，创作于2022年6月第十六张，PPT共一百零五页，创作于2022年6月v313 三角保持器v 三角保持器是一种理想保持器，物理上不能实现，数学上也是不能实现的，除非它所再现的信号为一已知信号。这一点从下面的定义中可以看出。它的脉冲传递函数曲线和频率特性曲线如图311和图312所示。第十七张，PPT共一百零五页，创作于2022年6月第十八张，PPT共一百零五页，创作于2022年6月第十九张，PPT共一百零五页，创作于2022年6月v 在实际中有时使用滞后一拍的三角保持器，它的脉冲传递函数曲线和频率特性曲线如图313和图314
6、所示。第二十张，PPT共一百零五页，创作于2022年6月v 经过三角保持器及滞后三角保持器再现后的函数如图315和图316所示。第二十一张，PPT共一百零五页，创作于2022年6月v 从上面的保持器特性可以看出，实际保持器与理想保持器的特性总是有差别的，所以要想使保持器引起的失真足够小，采样频率就要足够高。也就是说，在仿真计算时，为了使结果准确，计算步距就得足够小。这样势必要增加计算时间。为了使计算速度较快又不使误差过大，通常在保持器后加一补偿环节(也可加在保持器前)。保持器再现被采样的连续信号时，这些信号一般都有相位移，而且再现信号同其被采样的连续信号相比，都有所衰减。所以通常采用超前装置进
7、行补偿。例如，当采用零阶保持器再现信号时，再现后的信号要比被采样的连续信号平均滞后T/2。故应采用超前半个周期的补偿(即取ceTs/2)去抵消零阶再现过程引入的滞后影响。v 在仿真中所采用的补偿器的数学表达式形式一般为v C=erTs (3-11)第二十二张，PPT共一百零五页，创作于2022年6月v32 离散系统差分方程的求取v 在上一节中，已经叙述了怎样由一个连续系统求出它的离散相似系统。有了这个离散相似系统，就可以求出连续系统的离散化数学模型。离散化数学模型是用差分方程表示的，它的求解方法如下。v 设线性定常系统的状态方程描述如式(3-12)所示。v v式中A、B均为常数阵。v推导可得：
8、第二十三张，PPT共一百零五页，创作于2022年6月第二十四张，PPT共一百零五页，创作于2022年6月第二十五张，PPT共一百零五页，创作于2022年6月第二十六张，PPT共一百零五页，创作于2022年6月v系统的离散解：v (3-23)v 在推导上式的过程中未作任何近似的假设，该式是一种精确的采样值计算公式。但是，当U()是一个复杂的函数时，该式右端的积分是难以求得的。由于该积分的积分区间仅为T，当T较小时，一般来说U()在这个积分区间的变化是不大的。因此，可以加入采样及再现环节，以使U()在积分区间内为一个简单的特殊函数，从而使该积分计算容易进行。通常使用下面三种保持器作为再现环节。第二
9、十七张，PPT共一百零五页，创作于2022年6月v1.零阶保持器第二十八张，PPT共一百零五页，创作于2022年6月v2.一阶保持器第二十九张，PPT共一百零五页，创作于2022年6月v 3.三角保持器v 用与上述相同的方法，可以求出加三角保持器和滞后三角保持器时的系统差分方程。v用各种保持器时的差分方程，详见表31。第三十张，PPT共一百零五页，创作于2022年6月第三十一张，PPT共一百零五页，创作于2022年6月v例31 己知某控制系统框图如图317所示，求该系统的仿真模型，即差分方程。第三十二张，PPT共一百零五页，创作于2022年6月第三十三张，PPT共一百零五页，创作于2022年6
10、月第三十四张，PPT共一百零五页，创作于2022年6月第三十五张，PPT共一百零五页，创作于2022年6月第三十六张，PPT共一百零五页，创作于2022年6月v 如果把采样开关保持器加在系统入口R处，则得到的离散相似系统框图如图319所示。第三十七张，PPT共一百零五页，创作于2022年6月第三十八张，PPT共一百零五页，创作于2022年6月第三十九张，PPT共一百零五页，创作于2022年6月v于是得到零阶保持器时的差分方程：第四十张，PPT共一百零五页，创作于2022年6月v 由此可见，采样器及保持器的位置不同，得到的差分方程也不相同。但应注意，无论离散-再现环节加到哪里，被离散再现的信号都
11、应是状态方程中的输入量。v 从保持器的定义式可以看出，零阶保持器能无失真地再现阶跃输入信号，即当输入信号为阶跃函数时，导出的差分方程是精确的。而三角保持器能无失真地再现斜坡输入信号。第四十一张，PPT共一百零五页，创作于2022年6月第四十二张，PPT共一百零五页，创作于2022年6月第四十三张，PPT共一百零五页，创作于2022年6月v 例32 求如图321所示系统的差分方程。第四十四张，PPT共一百零五页，创作于2022年6月第四十五张，PPT共一百零五页，创作于2022年6月第四十六张，PPT共一百零五页，创作于2022年6月第四十七张，PPT共一百零五页，创作于2022年6月第四十八张
12、，PPT共一百零五页，创作于2022年6月第四十九张，PPT共一百零五页，创作于2022年6月第五十张，PPT共一百零五页，创作于2022年6月第五十一张，PPT共一百零五页，创作于2022年6月第五十二张，PPT共一百零五页，创作于2022年6月v 对于许多控制系统，都可以把系统分解成由积分和惯性环节组成系统。因此，可以事先求出这两个环节的差分方程的通用式，以后就不需要每次求解差分方程了。v 对于积分环节，如图3.23所示。第五十三张，PPT共一百零五页，创作于2022年6月 对于惯性环节，如图3.24所示。第五十四张，PPT共一百零五页，创作于2022年6月第五十五张，PPT共一百零五页，
13、创作于2022年6月v33 连续系统数字仿真程序设计v 前两节讲述了连续系统离散化及其差分方程的求取。有了系统的差分方程，就可以在计算机上编制仿真程序了。v 随系统性质的不同，对仿真程序要求也不同，一般的要求是计算速度快、精度高、使用方便、通用性强等。但这些要求往往是相互矛盾的，所以具体到某一问题时，应根据其特性突出某一要求而牺牲另外一些要求。第五十六张，PPT共一百零五页，创作于2022年6月v 在一般情况下，仿真程序由以下几个基本模块构成：v (1)初始化程序块；v (2)输入参数程序块；v (3)主运行程序块；v (4)输出仿真结果程序块。v 仿真程序流程图如图325所示。下面分别讨论这
14、三个程序块。第五十七张，PPT共一百零五页，创作于2022年6月第五十八张，PPT共一百零五页，创作于2022年6月v331 初始化程序块v 这个功能块主要是对程序中所用到的变量、数组等进行定义，并赋以初值。在通用仿真程序里这个程序完成被仿真系统的结构组态。v 这个程序块的内容随使用的程序设计语言的不同而不同，没有统一的格式。本书所有的程序均使用C，matlab和Visual C+6.0作为程序设计语言。v 自20世纪90年代以来，随着多媒体技术和图像技术的蓬勃发展，可视化技术得到了广泛重视，越来越多的计算机专业人员和非专业人员都开始研究并应用可视化技术。一般讲，可视化技术包含两个方面的含义：
15、一是软件开发阶段的可视化，即可视化编程，它使编程工作成为一件轻松愉快、饶有趣味的工作；二是通过可视化窗口将不改变的参数输入给计算机。v 在通用仿真程序中，可视化输入参数程序块是非常复杂的，该程序块直接关系到人-机交互的方便性和程序的通用性。一般人们会花很大的精力来设计输入参数程序块。在后面的章节中将逐步介绍该程序块的设计方法。第五十九张，PPT共一百零五页，创作于2022年6月v332 输出仿真结果程序块v 该程序块输出仿真结果，它们可以是状态变量、中间变量、输出变量的仿真结果。输出的形式总是数据表格或曲线的形式。下面介绍这两种输出程序的设计方法。v 在以后的程序中只写出调用这两个输出仿真结果
16、函数的语句，而不再列写程序清单。v1打印数据函数Print Data()v 一般要求打印数据函数能打印各变量的仿真数据及各数据对应的时间。由于仿真点数较多，输出数据量较大，用户不一定需要观察每一步的数据。因此，往往时隔一定的间隔打印一个数据。v 下面的函数可打印14个变量(变量个数用VN表示)的仿真数据及其对应的时间，每隔IN点打印一个数据。并要求在调用的程序里把需输出的变量值存放在数组Output里，其中第一维下标表示第几个变量，第二维下标表示该变量的第几个仿真数据。LP、DT的意义表示计算步数和采样周期。第六十张，PPT共一百零五页，创作于2022年6月第六十一张，PPT共一百零五页，创作
17、于2022年6月v 2显示曲线函数Dispcurve()v 该程序可以在显示器上同时显示l8条响应曲线。横轴为响应时间，纵轴为响应数值。要求在调用的程序里把需要输出的变量值存放在数组Output里，其中，第一维下标表示第几个变量，第二维下标表示该变量的第几个仿真数据。VN表示要显示变量的个数，LP、DT的意义表示计算步数和采样周期。第六十二张，PPT共一百零五页，创作于2022年6月第六十三张，PPT共一百零五页，创作于2022年6月第六十四张，PPT共一百零五页，创作于2022年6月v333 主运行程序块v 这个程序块用来求解被仿真系统的差分方程。所选用的仿真算法不同，得到的差分方程也不同，
18、仿真精度也不一样。不管怎样，这个程序要忠实于原差分方程，它不能改变原差分方程的意义，对于初编程序者来说在这方面是很容易出错的。v 从上一节可以看到，系统的仿真模型都是差分方程的形式。因为差分方程最容易用数字计算机求得数值解。差分方程的一般形式如式(2-8)所示。如果改变一下符号并写成递推计算的形式 第六十五张，PPT共一百零五页，创作于2022年6月第六十六张，PPT共一百零五页，创作于2022年6月第六十七张，PPT共一百零五页，创作于2022年6月第六十八张，PPT共一百零五页，创作于2022年6月v例33 设计例32所述的控制系统仿真程序中的主运行程序块。取R为单位阶跃函数，状态变量初值
19、为零。v 第六十九张，PPT共一百零五页，创作于2022年6月第七十张，PPT共一百零五页，创作于2022年6月第七十一张，PPT共一百零五页，创作于2022年6月第七十二张，PPT共一百零五页，创作于2022年6月v 从上面的仿真程序设计过程中可以看出，由于串行计算的原因，三角保持器不是放在哪里都可以实现的。实际上，不只是三角保持器，有超前作用的其他保持器或加超前补偿的零阶、一阶保持器也不是放在哪里都可以实现的，因为在有超前作用的保持器下，其差分方程的输入项可能会需要k+1时刻的值。这样，在系统的反馈支路与主支路求和点后的第一个环节入口处，不能加有超前作用的离散-再现环节。第七十三张，PPT
20、共一百零五页，创作于2022年6月第七十四张，PPT共一百零五页，创作于2022年6月第七十五张，PPT共一百零五页，创作于2022年6月v34 典型非线性环节的仿真程序设计v 一般工业系统均存在非线性问题。例如，调节仪表的输出可能存在饱和非线性，执行机构可能存在齿轮间隙非线性，还有些系统采用非线性调节器。因此，应当考虑非线性系统的仿真问题。v 实际上，前面介绍的按环节离散化的仿真方法，可以很容易地推广到具有典型非线性环节的非线性系统的仿真。只要事先编好这些典型非线性环节的仿真子程序，让信号通过它再输入到线性环节就行了。为了以后编程的方便，下面给出几种典型非线性环节的仿真子程序。第七十六张，P
21、PT共一百零五页，创作于2022年6月v1继电器v 继电器的非线性特性如图329所示。图中y0为继电器常数，它的特性可用如下的数学表主式描述：v 当u0时，y-y0v 当u0时，yy0。第七十七张，PPT共一百零五页，创作于2022年6月v2有不灵敏区的继电器v 具有不灵敏区的继电器特性如图330所示，图中yo、c为常数。其数学表达式为v 当uc时，yy0v 当|u|c时，y0；v 当uc时，y-y0第七十八张，PPT共一百零五页，创作于2022年6月v3限幅器v限幅器的特性如图331所示，其数学表达式为v 当uc时，y=cv 当|u|c时，y=uv 当u-c时，y=-cv式中：c为常数。第七
22、十九张，PPT共一百零五页，创作于2022年6月v4不灵敏区v不灵敏区特性如图332所示，图中c为常数，其数学表达式为v 当uc时，y=u-cv 当|u|c时，y=0v 当u-c时，y=u+c第八十张，PPT共一百零五页，创作于2022年6月v5齿轮间隙非线性第八十一张，PPT共一百零五页，创作于2022年6月v6具有死区和滞环的继电非线性第八十二张，PPT共一百零五页，创作于2022年6月v例34 某发电机励磁系统框图如图335所示，试求URFF为单位阶跃函数时系统的响应。第八十三张，PPT共一百零五页，创作于2022年6月第八十四张，PPT共一百零五页，创作于2022年6月第八十五张，PP
23、T共一百零五页，创作于2022年6月第八十六张，PPT共一百零五页，创作于2022年6月第八十七张，PPT共一百零五页，创作于2022年6月v35 关于采样周期(计算步距)和仿真时间的选择v 计算步距和仿真时间是两个重要的仿真参数，特别是计算步距。如果选择得不恰当，就可能造成较大的计算误差，甚至可以使一个本来稳定的系统歪曲成一个不稳定的系统。计算步距不仅和被仿真的系统有关，还和仿真算法、精度要求等因素有关。因此，要在仿真计算之前准确地选好这两个参数是件不容易的事情。第八十八张，PPT共一百零五页，创作于2022年6月v351 计算步距的估计 v 计算步距可用式(3-71)估计v (3-71)v
24、式中：DT为计算步距；n为被控对象传递函数的阶次；T为被控对象传递函数的时间常数。第八十九张，PPT共一百零五页，创作于2022年6月v 如果被控对象有若干个，则应以其中nT最小的为准。一般使用者按上述区间选择一个适当的计算步距，其仿真结果是令人满意的。一般来说，计算步距选择的越小，计算精度就越高(但步距太小，也可能由于计算机的舍入误差占了主导地位而降低了计算精度)，耗费的计算时间就越长；反之，就会导致计算精度差，计算时间短。但对于一般工程设计来说，最关心的是系统的稳定性，能够大致观察过渡过程曲线就可以了，对计算数据精确度的要求并不太高，所以只需采用试算的方法求计算步距即可。如果过渡过程曲线的
25、第一个半波含有20到30个计算点，则说明计算步距选择是恰当的。否则，可适当增大或减小计算步距，再试算。第九十张，PPT共一百零五页，创作于2022年6月v例35 某火力发电厂汽包水位控制系统采用单级三冲量调节方式，系统框图如图338所示。对此系统进行仿真，绘制蒸汽量D扰动时，水位H的变化曲线。第九十一张，PPT共一百零五页，创作于2022年6月第九十二张，PPT共一百零五页，创作于2022年6月第九十三张，PPT共一百零五页，创作于2022年6月状态方程第九十四张，PPT共一百零五页，创作于2022年6月差分方程第九十五张，PPT共一百零五页，创作于2022年6月在该系统中，最小的nT为所以取
26、 第九十六张，PPT共一百零五页，创作于2022年6月第九十七张，PPT共一百零五页，创作于2022年6月v352 仿真时间的估计v 一般仿真时间取为过渡过程时间，即从加入扰动开始到系统基本稳定为止的时间。如果主要是为了观察系统的稳定性，仅计算系统响应的34个周期就足够了。所以仿真时间的估算公式可选为v ST(520)nT (3-73)v式中：ST为仿真时间；n为被控对象传递函数的阶次；T为被控对象传递函数的时间常数。v 如果被控对象有若干个，则应以其中nT最大的为准。v 仿真时间的大小不影响仿真精度，因此，也可以用试算的方法选择仿真时间。第九十八张，PPT共一百零五页，创作于2022年6月v
27、例36 某火力发电厂主汽温调节系统采用导前微分调节方式，系统框图如图341所示。要求对此系统进行仿真，输出导前区汽温2及主汽温1的仿真结果。第九十九张，PPT共一百零五页，创作于2022年6月第一百张，PPT共一百零五页，创作于2022年6月第一百零一张，PPT共一百零五页，创作于2022年6月由图342可得到系统的状态方程：第一百零二张，PPT共一百零五页，创作于2022年6月根据表32即可得到离散相似系统的差分方程：第一百零三张，PPT共一百零五页，创作于2022年6月第一百零四张，PPT共一百零五页，创作于2022年6月感感谢谢大大家家观观看看第一百零五张，PPT共一百零五页，创作于2022年6月