计算机仿真教案.pdf
《计算机仿真教案.pdf》由会员分享,可在线阅读,更多相关《计算机仿真教案.pdf(41页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、计算机仿真第一章绪论1.1 仿真技术概述1.2 仿真研究的步骤L3系统仿真的特点1.1 仿真技术概述1.1.1 仿真的定义仿真是指用模型(物理模型或数学模型)代替实际系统进行实验和研究。1 .L 2仿真遵循的原则为使仿真的结果被证实是真实可靠的,也就是结果是可信的,仿真所遵循的基本原则是相似原理:几何相似、环境相似和性能相似。1 .1.3仿真的分类1按模型的性质分:物理仿真、数学仿真、数字物理混合仿真。物理仿真是以物理过程相似、几何尺寸相似以及环境条件相似为基础的仿真。优点是最大限度地反映系统的物理本质,具有直观性强和形象化的特点;缺点是构造物理模型费用高、周期长、技术复杂、实验条件限制多并易
2、受环境影响、修改结构和参数困难。数字仿真是应用性能相似、环境相似的原理,按照真实系统的数学关系,构造系统的数学模型,并在数学模型基础上进行试验。其特点是经济、参数修改方便、周期短,但形式抽象、直观性差。数学物理混合仿真是在某些系统的研究中,把数学模型与物理模型或实物联接在一起进行实验的形式,称为数学物理混合仿真或半实物仿真。2实时仿真与非实时仿真按照仿真实验时间标尺T与实际系统的时间标尺t比例来分类,将T/t=l的仿真称为实时仿真,而T /t#l的仿真为非实时仿真。一般有实物介入的半实物仿真属于实时仿真,而无实物介入的纯计算机仿真为非实时仿真。3按计算机仿真采用的计算机分类用于计算机仿真的一套
3、软硬件设备构成了仿真系统。根据仿真系统采用的计算机种类,可分为:模拟计算机仿真、数字计算机仿真和数字模拟混合仿真。1.2仿真研究的步骤仿真是研究系统普遍采用的方法。仿真步骤主要有:1建立系统的数学模型数学模型是系统仿真的研究依据,所以数学模型的准确性是十分重要的。2建立仿真模型一般的数学模型都不能直接编制程序并用计算机求解,通常必须把数学模型转换成适宜编程并能在计算机上运行的模型。也就是需要通过一定的算法对原系统的数学模型进行离散化处理,就连续系统而言,就是建立相应的差分方程再由计算机进行求解。3编制仿真程序对于非实时仿真,可用一般高级语言依据相应的算法编程。而对于实时仿真往往采用汇编语言与高
4、级语言共用的方式进行编程。4程序调试、验证模型、实验结果分析并确定实验方案L3系统仿真的特点1仿真的实验性质仿真基本上是一种通过实验来研究系统的综合实验技术,具有一般实验的性质。不论是利用仿真技术进行系统分析还是设计,都必须通过一系统实验来完成。所建立的仿真模型应具有实验的性质,即模型与原型的功能及参数之间应有相似性和对应性2数字仿真与解析法解析法是应用数学推导、演绎去求解数学模型的方法,而数字仿真是通过计算机在数学模型上进行一系列实验来研究问题的方法。1解析法求解问题得出对问题的通解,而数字仿真的每一次运行只能给出在特定条件的数值解(犍星)。2用解析法所能求解的问题是有限的,而仿真适应范围广
5、。3数字仿真与数值求解数字仿真是指建立系统的数学模型并在计算机上运行和分析的整个过程;数值求解则是仿真所采用的在计算机上求解数学模型的方法。第二章连续系统数值积分法的时域数字仿真2.1连续系统的数学模型2.2模型转化一实现问题2.3数值积分法2.4数值积分法的计算稳定性2.5数值积分法的选择原则2.6 一般数字仿真程序的主要功能和结构建立描述系统运动的数学模型是数字仿真的先决条件。模型的形式不止一种,它们各有特长和最适用的场合,彼此之间也有紧密的联系,要能灵活地运用数学模型分析具体问题,就要掌握这些数学工具的使用和彼此间的相互转换。2.1 连续系统的数学模型1连续系统的定义:凡是系统的输入量是
6、时间t的连续函数,无论其输出量是连续单调函数或是单值或是多值的函数,均称这个系统是连续系统。我们所研究的系统是线性、定常的连续系统。2线性定常连续系统的数学模型一个连续系统的数学模型可以用微分方程、传递函数和状态空间表达式等形式表示。数学模型的形式(1)微分方程设连续系统的输出量)C),输出量为(D则系统用微分方程表示为也+%。+%*+.+*生+“dtn 1 dtn-dtn-2 I dt d-h d-2uc dt-+C1 dt-2(2)传递函数对微分方程等号两边逐项进行拉氏变换,并考虑初值为零时,可得到sYG)+6sTy G)+q 2 s ly +J。)=cas-U G)+c U(s)+C2S
7、-3U(S)+-+C,IU(S)式 中 丫(5)一系统输出的拉氏变换U(s)一系统输入的拉氏变换系统的传递函数G(s)=8 2 =。.产+%严2+7+.+,:U(s)sn+atsl+a 2sli +-+an微分方程或传递函数是用系统的输入、输出之间的关系来描述系统的,表示了系统的外部特征,所以称其为外部模型。用微分方程表示的系统可以是非线性或线性系统,而对于传递函数表示的系统,只适用于单输入一单输出的线性定常系统,所以传递函数的模型表示有一定的局限性。(3)状态空间表达式状态空间表达式可以由两个途径获得,由微分方程、系统模拟图导出,本节将仅对微分方程推导作简单说明。设系统由不含输入量导数项的n
8、阶微分方程表示。dy d-ly dn-2y dy-F C l.-+u?-+I F ci y=udt dt-2 d Q d t定义为n个状态变量为令X=y_._dy_x2-x1-dt.d2y.dy写出各个状态变量的一阶微分方程形式:x2X i=%=一囚乙 x.-i-anx+M将上述n 个一阶微分方程写成矩阵向量形式N=Ax+B”(状态方程)y=Cx(输出方程)(合称状态空间表达式)式中A,B,C 为系数矩阵,x 为状态变量向量。状态空间表达式也可由传递函数求得。由于状态空间表达式反映系统内部特性,又称为内部模型。实际上内外部模型可以互相转化,对于一个外部模型所对应的状态空间表达式,由于状态变量有
9、不同的取法,因而可得到多种形式的状态空间表达式,这时多种状态空间表达式将对应同一个外部模型。也就是说对于一个确定的外部模型只对应一个确定的系统。2.2 模型转化一实现问题因为状态方程是一阶微分方程组,非常适宜用数字计算机求解,如果一个系统是用状态空间表达式描述的,便可直接编程求解。然后对于一些复杂的控制系统,其数学模型往往是通过实验得到的数据,经过辨识确定的,一般是用传递函数表示,为了便于数字仿真,必须把它们转换为状态空间表达式。把从已知系统传递函数去求相应的状态空间表达式的模型转化问题称为实现问题。1传递函数转换成状态空间表达式转换采用的方法是状态变量图法,这个方法是从系统模拟方法中出来的,
10、所谓系统模拟方法是指用一些基本模拟单元,如加法器、积分器、放大器等去模拟用传递函数表示的复杂系统。用基本模拟单元模拟系统的传递函数得到的图形是系统模拟图,在系统模拟图上标上状态变量的图形是状态变量图,然后再求出状态空间表达式。举例说明,如何状态变量图法得到传递函数转换成的空间表达式。如初始条件为零的积分器,代表一个简单的系统。系统的传递函数为G(s)=l/s。该系统的模拟图如图1-a所示,设系统的状态变量为x,在积分器的输入标上状态变量微分比,在积分器的输出标上状态变量x,得到系统的状态变量图如1-b所示。”(s)1-.X V X1-a 1-b积分器的系统模拟图和状态变量图由 状 态 变 量
11、图 根 据 积 分 器 的 输 入、输 出 关 系 写 出:状 态 方 程*=和输出方程y=x,这两个方程表示了传递函数转换的状态空间表达式。对于带反馈的积分器,其传递函数为 G(s)=一s+a其 中。是反馈系数,它的状态变量图如下图2所示,其中 是放大器,Z是加法器。可见系统由积分器、加法器和放大器组成。由积分器输入、输出关系得到 fx=-ax+u1y=x由以上得到由系统模拟图到状态变量图并导出状态空间表达式的步骤如下:1根据系统的传递函数,画出系统模拟图,阶系统有 个积分器;2把积分器输出定为状态变量x,积分器输入处定为支,并分别标在积分器输入和输出处,得到状态变量图。3根据积分器输入、输
12、出方程写出系统的状态方程及输出方程。对于高阶、复杂系统采用级联法、串联法和并联法得到代表实际系统传递函数的系统模拟图及相应的状态变量图,依据同样的方法求得状态空间表达式。通过下面的例子来说明这三种方法并推广到一般情况的分析。一个三阶系统的传递函数有以下三种形式:a)分子、分母为微分方程的多项式;b)将分子分母为微分方程的多项式化成几个一阶环节传递函数串联的乘积;c)将分子分母为微分方程的多项式化成几个一阶环节传递函数并联之和相应的三种形式的多项式为:分子、分母为微分方程的多项式G(s)=s+3s+2?+752+12.V多项式化成几个一阶环节传递函数串联的乘积G(s)s(s+3)(s+4)多项式
13、化成几个一阶环节传递函数并联之和G(s)1/62/3 3/2-1-s+3 s+42状态空间表达式转换为传递函数已知系统状态空间表达式 候=_4+8“y =C x1)对上式两边取拉氏变换s X(s)-x(0)=A X(s)+8 U(s)1Y(s)=C X(s).22)在零初始条件下整理式1,并在等式两边乘以单位矩阵(s/-A)X(s)=B U(s)X(s)=(s/3将X(s)代入式2中,得到 G(s)=幺2=C(sl-4尸 BU(s)Y(s)=C(sI-A)-BU(s)系数矩阵A、B、C已知的情况下,根据上式就可求出状态空间表达式所对应的外部模型的传递函数。2.3数值积分法1数值积分的基本概念:
14、已知一个一阶微分方程 才=/(r,x)应用数值积分法求解的思路是。从初值与开始,在一系列时刻tt-t0+h,t2=t1+h,.tn=tn_x+h,各点求未知解x”为,七,的近似解,h是计算步长,若X )对的各阶导数都存在。则x 在“f +人 时 的 解x(r 0+)用台劳级数表示为x(t0+/?)=尤 o)+hx(t0)+h2 X +.+hk x,)+.2!k2欧拉法一般形式为:x,=X,-1 +hf(tn_xxn_,)欧拉计算公式是一个递推差分方程欧拉数值积分方程的几何意义:折线法;矩形法。截断误差力2 力3 nkR=不/(,0,%)+,/(无,%)+./xo)3梯 形 法(二阶龙格一库塔法
15、)台劳公式的前三项:x o +%)=x o)+/z x(%)+%2 X 2(;)推导出二阶龙格库塔计算公式:X|=x0+g(M+K2)-校正公式K =/Q o,Xo);K-,=/&+九 Xo +K-)预报公式如果对一个点的计算反复预报、校正数次,得到梯形法的计算公式预报 燔=七,+/矿仇,/)校正 2f2 fnn_mn _4P“2,P.“_式中一仿真系统的环节数;=%,”2,”了-各环节的输入信号;m 叫,加2,,血-外输入信号;y=S,2,,以7-各环节的输出信号,尸为外输入连接矩阵,P为环节连接矩阵。3.4非线性环节特性分析按环节离散化的仿真方法,很容易推广到具有典型非线性环节的非线性系统
16、的数字仿真,这是因为按环节离散化的数字仿真中,每增加一个步长都要计算出所有环节的输出I和输入M,这样就可以在两个线性环节之间插入非线性环节的仿真程序,即 前-个 环 节 的 输 出x不是直接作为下一个环节的输入,而 是 先 将 送入非线性环节,非线性环节的输出再作为下一个与其相连接环节的输入/T,下页图就表示了两个线性环节中间加入了非线性环节。在实际系统中,非线性环节可能在线性环节的前面或后面,非线性环节的位置不同,其相应的计算方法也随之不同,但总的结果不变。在仿真程序中安排下列几种非线性环节。饱和非线性分段线性的非线性的特性。死区非线性这类特性在控制系统中是较多见的,如放大器的饱和、死区等。
17、齿隙(间隙)非线性在有机械传动的系统中都存在着齿轮间隙(传动间隙)特性,这类特性同系统精度要求相比不能忽略时,就必须分析它对系统性能的影响。同样在电器产品中使用了磁性材料,由于磁性材料磁滞环也是同属于间隙类的非线性。继电及具有死区的继电非线性具有滞环的继电非线性具有死区和滞环的继电非线性第三章小结按环节离散化的数字仿真方法是将连续系统进行虚拟的离散化处理后,借助于离散系统的理论和方法,建立与连续系统模型等价的离散模型进行仿真计算的方法。该方法的物理意义明确,概念清晰直观,提示了数学仿真的本质,在时域用离散化状态方法求解差分方程。对于每一系统在计算步长确定后,离散系数矩阵(7),“(T),&(T
18、),只要求出就不再改变,使求解计算特别简单,加快了仿真速度,在工程应用中,当系统环节数小于20时-,计算步长取力=1,即可满足精度要求。(1O5O)W,第 四 章 快 速 数 字 仿 真前面介绍了应用数值积分法和按环节离散化方法进行控制系统时域数字仿真的基本原理、方法和程序,对于控制系统的数字仿真是比较方便的。但是这些数字仿真程序的一个共同的特点是:为使仿真达到一定的精度,计算步长不能太大,因此使计算量加大,计算时间加长,仿真运行过程增加,不能适应缩短仿真时间的要求,这就提出快速性要求与精度之间的矛盾,对于实时仿真(在线仿真)的场合更加无法满足要求。在对高阶或复杂系统进行数字仿真,往往仿真实际
19、运行时间需要数小时甚至几十小时 一,将使仿真无法进行或失去实际意义。而在有些控制仿真时;由于参数的限制,计算步长要取得极小,这将使数字仿真时间很长,尤其是如果系统中的两个环节的时间常数相差达几十倍或者上百倍时,如按小的时间常数来确定.计算步长,将会引起计算不稳定以致整个系统不稳定,仿真无法进行。在实时仿真中为达到快速性的要求,会不惜牺牲一定精度条件下进行。为了解决数字仿真速度的措施是:硬件措施一采用调整计算机。软件措施一采用执行速度快的语言,如汇编语言。算法措施一研究并采用适合于快速仿真的仿真算法。本章讲解快速仿真算法。一般情况下,在解决计算速度与精度这一对矛盾时,以计算速度为矛盾的主要方面,
20、在满足计算稳定性及工程要求精度条件下,尽可能提高仿真计算速度,降低不必要的精度要求以满足速度要求。一般对快速仿真算法有两点基本要求减小每一步的计算量;算法具有良好的计算稳定性,在保证一定的仿真精度的情况下,允许采用较大的计算步长。4.1 替换法替换法是建立在相匹配原理基础上的快速仿真方法之一。1 相匹配原理相匹配的含意是指原系统的数学模型的动态的和静态特性是一致的。其定义是:若被仿真系统的数学模型是稳定的,且动态、静态特性一致。即若G(s)是稳定的,而仿真模型G(z)也应是稳定的,对于同一输入信号,由 G(z)求 得 的 输 出 y/T)与由G G)求 得 的 y 应具有相一致的时域特性或具有
21、相一致的频率特性。一般快速仿真提高了对速度的要求,合理地降低了对精度的要求,因此如果能根据相匹配原理,直接对一个高阶连续系统的数学模型G(s)转换成与之匹配、每步计算量较小。允许采用较大步长且具有合理精度的仿真模型G(z),那么可利用求得差分方程进行数字仿真,以达到快速仿真的效果。本节介绍的替换法及下节的根匹配法都是根据相匹配原理由G(z)直接导出G G)的变换方法。2替换法原理一般情况下,由 G(s)转换成G(z)采用下面两种方法00由G(s)求 g(f),然后 按 公 式 G(z)=Z g(T)/求 出 G(z);”=0将G(s)展开成部分形式,然后通过查变换表得到G(z)。对于高阶系统,
22、采用这两种变换方法都是困难和繁杂的,所以应该寻找一个由G G)直接转换成G(z)的简便方法。3替换法步骤采用双线性变换法求仿真模型的具体步骤如下:设已知一线性系统的传递函数为G(s)as+bs(s +c)(1)将 G(浒 s 按s =2(上H)T 1+/代入上式,可得GT 1+z-i,2I、2/2T(-r)+(-r)cT+z T +z整理化简得G(z)2aT+bT2-2bT2z-+(-aT+bT2)z24 +2 cT 8 茶;(4 2c T*(2)将丫(外=6(7 仁)转换成差分方程y()-(2aT+bT2)u(n)-2bT2+w(n-l)+(bT2-2aT)u(n-2)+8 y(-1)-(4
23、-2cT)y(n-2)4.2根匹配法根匹配快速数字仿真是根据相匹配原理提出的又一种常用的快速数字仿真方法。1 根匹配法原理根匹配法快速数字仿真的原理与替换法有相似之处,但也有不同之处。相似之处是它也是将系统的传递函数G(s)转换成相似z域传递函%仁),再转换成差分方程进行快速数字仿真,不同之处是它不采用近似公式使 转换成,而是根据连续系统的动态、静态特性由传递函数的零、极点分布及增益所确定,采样系统的动态、静态特性由传递函数的零、极点分布及增益所确定。如果根据匹配原理,找到一种由G(s招零、极 点 到G(z油零、极点的简单映射关系z =/G),然后根据其特性确定其增益,则可方便地由G(s)求得
24、与之等价的G(z)。这就是根匹配的思想,也就是说根匹配的思想就是使G(z)的根轨迹与G(s)的根轨迹相匹配,从而使它们的动、静态特性相同。4.3增广矩阵法1增广矩阵的基本思想增广矩阵法是按环节离散化数字仿真方法的改进。按环节离散化的方法是非常实用的仿真方法,但仿真精度不高,产生误差的原因有两条:在离散化过程中,采用了虚拟采样开关和保持器,引起输入信号的失真;对于高阶系统的系数矩阵,用解析法计算比较繁琐,需要用数值计算方法,从而引起截断误差,若为方便计算这些阶系统的系数矩阵,按一阶环节计算,但又会产生增加环节数使误差增大。上述原因影响到采样周期T的选取,不能取太大。如何使采样周期适当加大又不引起
25、误差加大这个问题?齐次矩阵方程的离散化,不需要加入虚拟的采样开关和保持器,由此引入的误差自然就不存在了,同时也消除了采用数值计算方法计算产生的计算误差。误差的来源仅是采用数值方法计算产生的截断误差及计算机产生的舍入误差。这时T取得较大精度也可以得到保证。然而一般系统的状态方程都是非齐次的,应设法将非齐次项,即输入变量“增广”成为新的状态向量(即增广状态变量),使原方程变成增广形式的齐次状态方程。2典型输入信号的增广矩阵阶跃输入信号设“(,)=M()x l(f),I%为阶跃幅值,定义“(f)为n+1个状态变量则:X +l(f)=)=M o X l(f),%(f)=0 ,%(0)=劭根据上式得到增
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算机仿真 教案
限制150内