计算机仿真技术(2013)

机电工程学院机电工程系李安Tel:0532-86983310(O)15166391658Email:[e-mail][e-mail]本门课程的目标熟悉一个概念掌握一种方法精通一项技术会用一款软件参考文献1.系统仿真概论，肖田元等编，清华大学出版社2.现代仿真技术与应用，康风举编，国防工业出版社3.先进仿真技术与仿真环境，熊光楞编，国防工业出版社4.仿真技术与软件，贾连兴主编，国防工业出版社第1章概述ThreetopicstobediscussedWWHWhattolearnWhytolearnHowtolearn理解并掌握系统仿真的概念、分类、作用；了解计算机仿真的发展历程、应用领域及目前研究、应用现状；熟悉掌握计算机数字仿真的基本过程。本章学习要求1.1仿真的基本概念及其分类1.1.1系统仿真的定义、分类及作用1.1.2计算机仿真的定义及其分类1.2计算机仿真的发展及其应用领域1.3计算机数字仿真的基本过程1.4先进仿真技术1.1仿真（simulation)的基本概念及其分类1.1.1系统仿真的定义、分类及其作用1）系统仿真的定义系统仿真是通过对系统模型的某种操作，研究一个存在的或设计中的系统。简言之，系统仿真是对系统动态模型的实验，即在仿真中，系统的动态模型在一定的实验条件下被行为产生器驱动，产生模型行为。2）系统仿真的分类--系统模型分类方式①物理仿真系统模型为物理模型：实物模型(PM)。②数字仿真系统模型为数字模型：数学模型和几何模型（MM&GM）。③半物理仿真既有物理模型也有数字模型。3）系统仿真的作用⑴优化系统设计。在复杂的系统建立以前，能够通过改变仿真模型结构和调整参数来优化系统设计(MM)。⑵对系统或系统的某一部分进行性能评价(MM)。⑶节省费用。⑷重现系统故障，以便判断故障产生的原因(GM)。⑸可以避免试验的危险性。⑹进行系统抗干扰性能的分析研究(MM)。⑺训练系统操作人员(GM)。⑻系统仿真能为管理决策和技术决策提供依据。正因为仿真技术对国防建设、工农业生产及科学研究均具有极大的应用价值，所以，仿真技术被美国国家关键技术委员会于1991年确定为影响美国国家安全及繁荣的22项关键技术之一。本章目录1.1.2计算机仿真(ComputerSimulation)的定义及其分类1)计算机仿真的定义计算机仿真是指应用几何和性能相似原理，构成数字模型，在计算机上对系统数字模型进行某种操作。计算机仿真又称为数字仿真。①根据计算机分类模拟计算机仿真、数字计算机仿真、模拟数字混合仿真②根据仿真时钟与实际时钟的比例关系实时仿真、欠实时仿真、超实时仿真③根据系统模型的特性连续系统仿真、离散事件系统仿真本章目录⒉计算机仿真的分类1.2计算机仿真的发展及其应用领域1)计算机仿真的发展第一阶段:模拟仿真技术发展阶段,1950年以前。动因：火炮与飞行控制动力学系统的研究。第二阶段：混合仿真技术发展阶段，50～60年代。动因：洲际导弹和宇宙飞船飞行姿态及轨道控制动力学的研究。第三阶段：数字仿真技术发展阶段，70年代以后。模拟机—混合机—数字机计算机仿真的应用类型:①系统设计器对尚未有的系统进行设计时采用仿真技术；②系统分析器对已有系统进行分析时采用仿真技术；③系统预测器在系统运行前，利用仿真模型作为预测器，向用户提供系统运行起来后，可能产生什么现象，以便用户修订计划或决策；④系统观测器在系统运行时，利用仿真模型作为观测器，给用户提供过去、现在甚至未来的信息，以便用户实时作出正确的决策，比如利用仿真技术进行故障分析和故障处理；⑤系统训练器利用仿真模型作为训练器，训练系统操作人员或管理人员2)计算机仿真的应用领域从历史上看，仿真技术首先被应用于那些在实际系统进行试验有危险、花费巨大的领域，比如：航空、航天、武器系统等，以后逐渐扩大到虽然可在实际系统上进行试验，但花费较大、耗时较长、不大方便的一些领域，比如：冶金、化工、电力等；近十几年来，则进一步扩大到制造、交通、环境、生态、生物、石油等领域。综观仿真技术应用领域逐渐扩大的历史，可以得到以下几个有意义的结论：①仿真技术在应用上的安全性一直是被采用的最主要原因；②仿真技术在应用上的经济性也是被采用的十分重要的原因；③仿真技术一般是从学院式的局部应用逐步走向全面应用，其标志是对某领域的仿真应用而设立的仿真中心；④仿真技术在每一阶段都有一个比较成熟的应用领域；⑤根据仿真的应用范围，研制和开发了一大批仿真产品，如各种仿真语言及仿真软件包，各种训练仿真器等等。⑥为了推广仿真技术的应用，几乎所有的仿真语言都推出了PC版。本章目录1.3计算机数字仿真的基本过程为了使大家对计算机数字仿真有一个全面的了解，让我们用一个简单的例子来予以说明。卓越的操控性、稳定性和舒适性汽车悬架系统问题：如何分析研究并保证设计出或设计中的悬架系统具有这种卓越性能？研究方法：1理论方法2实验方法3仿真方法仿真技术与物理实验、理论研究的对比仿真技术物理实验理论研究可能性只要能建立系统模型，就能进行系统尚未建立，则不可能；有的自然系统实验周期太长，也不可能有的系统无法建立解析模型，因此，不可能利用解析方法安全性无危险有危险（人身、设备）无危险经济性花费不多费用很大花费少耗时性中等长短准确性可以做到很准确十分准确要做较多假设，有较大误差方便性可以做到十分方便受现场限制，不方便方便机械谐振系统其中：X—状态向量A—系统矩阵B—输入矩阵f(t)—输入变量C—输出矩阵D—直接转移矩阵输出方程状态方程DfCBfAtfkyByMyXYXX)(''0011010'2'1'DCMBMBMkAxxX2121'2'101f1010xxYMxxMBMkxx根据研究的目的编制试验方案：1.改变M、B、K，保持f(t)不变，研究系统结构参数对谐振系统动态行为的影响；2.改变f(t)，系统结构参数保持不变，研究系统所受激振力对谐振系统动态行为的影响；最后根据仿真试验结果，给出相关结论。建立仿真模型应用编程语言应用仿真语言或环境应用欧拉法、梯形法或R—K法等仿真算法，将上述数学模型转化为便于编程的仿真模型：2,1))(,(),(11ittxtgxxxtgdtdxnninnininii应用编程语言应用仿真语言或环境应用ADAMS应用MATLAB计算机仿真的一般步骤1、确定系统分析目标，进行系统分析，建立系统数学模型2、在系统数学模型的基础上，建立系统的仿真模型3、编程或选用仿真语言4、制定仿真实验方案，进行仿真实验5、分析仿真实验结果，修改模型本章目录给点小提示一、拉普拉斯变换1.拉氏变换的定义ja：sdtetftfLsFst其中0)()]([)(2.拉氏变换的几个定理①线性定理②衰减定理③延时定理④相似定理⑤微分定理⑥积分定理⑦终值定理⑧初值定理⑨卷积定理连续系统仿真)()]([,),()]([skFtkfLksFtfL则为常数若①线性定理)()()]([)()]([),()]([),()()(sGsHtfLsGtgLsHthLtgthtf则且若返回)()]([sFeatfLas)()]([asFtfeLat②衰减定理③延时定理④相似定理)()]([asaFfLat返回⑤微分定理)1()0()2()0()1(')0(22)0()()]([,)0()()]([,)0()()]([,nnnndtddtddtdfsffssFstfLnfsfsFstfLfssFtfLnn有阶导数对于有对于二阶导数有对于一阶导数GO⑥积分定理tssFdfL0)(])([⑦终值定理)(lim)(lim,)(lim,)(0')(ssFtftfftfsttt则存在且都符合拉氏变换和若返回⑧初值定理)(lim)(lim,)(lim,700ssFtftfstt则存在且条件同⑨卷积定理todtgftgtfsGsFtgtfL)()()(*)()()()](*)([其中二、拉氏逆变换直接求解法、查表法（系数比较法、留数法）0,)(21)()]([1tja：sdsesFjtfsFLjajast其中三、应用拉氏变换求解微分方程例:用拉氏变换解微分方程初始条件:teyyy32'''1,0)0(')0(yy四、控制系统的传递函数1.定义2.典型环节的传递函数①比例环节②微分环节③积分环节④惯性环节⑤二阶环节⑥延时环节1.定义)()()(sXsXsGio返回系统输入的拉氏变换系统输出的拉氏变换)()(0sXsXi①比例环节②微分环节③积分环节KsGtKxtxio)()()(TssGTtxdttdxoi)()()(sTsGdxTtxtio)()()(0④惯性环节⑤二阶环节⑥延时环节1)(TsKsGKxxTiodtdxo1)(02022sTTsKsGKxxTTiodtdxdtxdoosioesGtxtx)()()(返回第2章连续系统的数字仿真对于一个集中参数的连续动态系统，利用我们学过的基础理论知识和专业知识所能直接建立的数学模型形式有哪些呢？ubdtudbdtudbyadtdyadtydadtydammmmmnnnnnn111011110高阶微分方程传递函数nnnnmmmmaSaSaSabSbSbSbSUSY11101110)()(上述模型不能在数字机上直接求解！那么能够在数字计算机上直接求解的数学模型形式又是什么样的呢？))(,(11nnnnnnttytfyy差分方程上述方程形式由何而来的？一阶微分方程),(ytfdtdy仿真模型数值积分算法怎么来？动态方程BuAXXCXy高阶微分方程传递函数模型变换高阶微分方程传递函数差分方程两次模型变换动态方程等价变换近似变换2．1连续系统的数学模型ubdtudbdtudbyadtdyadtydadtydammmmmnnnnnn111011110njmjjjmjjnupbypa00njjjnmjjjmpapbuy00引进算子dtdp1°高阶微分方程2°传递函数)()()()()()()()()(111012211SUbSSUbSUSbSUSbSYaSSYaSYSaSYSaSYSmmmmnnnnn)()()(11101110SGaSaSaSabSbSbSbSUSYnnnnmmmm3°状态空间描述1．由微分方程导出状态空间表达式)(01111tubyadtdyadtydadtydnnnnnn引进如下状态变量：假定一个连续系统可用下式来描述)(0112111122123121tubxaxaxadtydxdtydxxdtydxxxdtdyyxxyxnnnnnnnnnnubxxxaaaaxxxXnnnnn021121210001000010Xy001nxxxX21令Aaaaannn12101000010Bb000C001则有BuAXXCXy例系统的微分方程式uyyyy67416设状态变量为321,,xxxuxxxxxxX