Matlab_Simulink动力学系统建模与仿真
《Matlab/Simulink动力学系统建模与仿真(第2版)》主要介绍了动力学系统中微分方程模型、传递函数模型和状态空间模型等建立的基础理论,并引入了Simulink仿真技术,为解决复杂动力学问题(特别是不易得到解析解的动力学问题)提供了方法。 《Matlab/Simulink动力学系统建模与仿真(第2版)》编排了较多的例题来说明各类动力学模型的仿真模型的建立方法,以及差分模型、相似模型、时域和频域等仿真模型,最后将控制动力学基础知识作为后继研究的扩展内容做了介绍。 《Matlab/Simulink动力学系统建模与仿真(第2版)》是一本多学科内容相交又的教材,同时涉及了力学、电学和动内容简介本书主要介绍了动力学系统中微分方程模型、传递函数模型和状态空间模型等建立的基础理论,并引入了 Simulink仿真技术,为解决复杂动力学问题(特别是不易得到解析解的动力学问题)提供了方法。书中编排了较多的例题来说明各类动力学模型的仿真模型的建立方法,以及差分模型、相似模型、时域和频域等仿真模型,最后将控制动力学基础知识作为后继硏究的扩展内容做了介绍。本书是一本多学科内容相交叉的教材,同时涉及了力学、电学和动力学控制等学科的交叉知识。本书适合具有一定数学和力学基础知识的理工科专业的本科高年级学生使用,可以作为机械工程、土木工程、车辆工程和仪器仪表、印刷机械等本科高年级学生和相关专业的研究生在学习有关动力学系统建模与仿真内容时的参考书,还可供相关工程技术人员参考。图书在版编目(CIP)数据Matlab/ Simulink动力学系统建模与仿真/黎明安,钱利编著.-2版.一北京:国防工业出版社,2015.7ISBN978-7-118-102055L.①M..Ⅱ.①黎.②钱..Ⅲ.①计算机辅助计算—应用—动力系统一系统建模②计算机辅助计算一应用—动力系统一系统仿真Ⅳ.①TP391.75②019中国版本图书馆CP数据核字(2015)第138118号※宫原社出版发行(北京市海淀区紫竹院南路23号邮政编码100048)天利华印刷装订有限公司印刷新华书店经售开本787×10921/16印张19字数456千字2015年7月第2版第1次印刷印数1—3000册定价49.00元(本书如有印装错误,我社负责调换)国防书店:(010)8854077发行邮购:(010)88540776发行传真:(010)88540755发行业务:(010)88540717前言本书是在为“工程力学”本科专业开设的“动态系统建模与计算机仿真”课程基础上多次改编而成的。本书一开始就采用了模型框图,使学生在学习过程中掌握和使用仿真框图的表示方法,为今后建立仿真模型奠定基础。本书结合了 Simulink仿真平台的基础知识,学生可以在各章的例题中学会 Matlab基本的编程能力和 Simulink基本模块的应用;将 Simulink的基础知识分散到各个章节中由浅入深地讲授,使学习者容易接受。本书第1版于2012年出版,此次在第1版的基础上修改了部分错误,个别习题也做了调整,使内容编排更趋于合理。全书分为10章,第1章~第3章介绍了建模与仿真的数学力学基础知识和以框图来表示模型的方法,主要以微分方程模型为主线介绍了简单仿真模型的建立。第4章介绍了系统的传递函数模型以及面对传递函数的仿真模型的建立。第5章介绍了状态空间模型。第6章介绍了基于采样的将连续系统离散化的方法。第7章介绍了机电相似模型。第8章介绍了动力学系统的时域瞬态响应分析方法。第9章介绍了频域分析方法。第10章介绍了控制动力学基础。全书贯穿了 Matlab/ Simulink仿真技术。本书中的仿真例题均在 Matlab(R2007a)下调试通过,希望读者在该版本下建立仿真模型本书由西安理工大学师俊平教授、王忠民教授审阅,两位教授对本书内容的编写提出了宝贵的修改意见,研究生雷霜、崔凯和朱晓雄等对初稿进行了认真的校对,在此表示衷心感谢。西北工业大学支希哲教授、朱西平教授,空军工程大学冯立富教授,陕西理学院张宝中教授,西安科技大学郭志勇教授,西安理工大学徐开亮博士等在编写过程中给予了大力帮助,在此表示深切感谢由于水平有限,本书还有很多需要改进的地方,敬请使用者提出宝贵意见。编著者Ⅲ目录绪论…………………………………………………………………………1第1章系统建模与仿真基础…………1.1系统仿真模型框图表示法·4441.1.1基本仿真元件………1.1.2简单仿真框图结构51.2拉普拉斯变换…1.2.1拉普拉斯变换的定义及其性质1.2.2拉普拉斯逆变换………………………………111.2.3拉普拉斯变换在求解线性常系数微分方程中的应用1.3z变换与Z逆变换…161.3.1Z变换的定义1.3.2Z变换的应用…171.4矩阵的特征值与特征矢量…181.4.1标准特征值问题…191.4.2广义特征值问题1.4.3相似变换及其特性…………………………………………21习题5第2章动力学系统的微分方程模型………………………………………282.1动力学建模基本理论…………………………………………………282.1.1动力学系统基本元件…282.1.2动力学建模基本定理…………………………292.2哈密顿动力学建模体系382.2.1拉格朗日方程2.2.2哈密顿原理2.3一维弹性体的有限元建模………422.3.1梁单元质量矩阵与刚度矩阵…422.3.2总体系统动力学微分方程………………………………442.4一维弹性体系统的假设模态法482.4.1模态函数……………………………………………482.4.2系统的动能和势能…………………492.4.3系统的动力学方程2.5 Simulink高级积分器的仿真模型建立…512.5.1高级积分器端口………522.5.2高级积分器在仿真中的应用…………………………52习题………………………………………………………………………………………54第3章动力学系统响应分析的数值方法……583.1数值积分法和数值微分法…583.1.1数值积分法………………………………………………………………583.1.2数值微分法……593.1.3多自由度振动系统的差商模型…………………………………633.2龙格一库塔法……653.2.1二阶龙格一库塔法………653.2.2四阶龙格一库塔法……………………………………………663.3四阶龙格一库塔法仿真程序设计673.3.1求解一阶微分方程四阶龙格一库塔法程序设计……………………673.3.2求解一阶微分方程组的四阶的龙格一库塔法程序设计693.3.3高阶微分方程的四阶龙格一库塔法程序设计703.4隐式逐步积分法…………………………………………………723.4.1线性加速度法723.4.2威尔逊θ法…3.5微分方程的边值问题的求解…………………763.5.1解线性方程边值问题的差分方法……………………………………763.5.2解线性方程边值问题的打靶法(试射法)…773.5.3关于三对角矩阵的追赶法程序设计……·,,,,,,,,,,,,.,,,,……793.6关于 Simulink环境中的求解器 Solver803.6.1常用求解器…………………………………………803.6.2求解器的选择813.7Malb中符号微积分……………823.7.1符号微分与符号积分……………………………………823.7.2利用符号运算求解微分方程习题……………………………………………………………………884第4章系统传递函数模型84.1传递函数及其特性…4.1.1传递函数定义……………884.1.2传递函数的特性…4.1.3传递函数的图示方法…894.2基本环节的传递函数……………………894.2.1比例环节……………………904.2.2一阶延迟环节4.2.3微分环节914.2.4积分环节…………914.2.5振荡环节………………914.3传递函数的其他形式934.3.1传递函数的零极点形式……………………………………………934.3.2传递函数的留数形式…………………………………………934.3.3传递函数的串联、并联与反馈连接形式944.3.4控制系统的开环传递函数……………………………………974.4多自由度振动系统的传递函数模型……………1014.4.1直接方法,..,.,,...,,,,,,.,·,,,,,,,,,..,·.、,,,,,,,,,,,,··,1014.4.2模态分析法…1035传递函数模型的 Simulink仿真模型建立…………………1054.5.1与传递函数相关的 Matlab运算指令…1054.5.2传递函数模型的 Simulink仿真模型建立…1084.6弹性系统的传递函数仿真模型…1114.6.1弹性系统的传递函数1114.6.2传递函数 Simulink仿真模型………………………………………………112习题113第5章动力学系统状态空间模型1175.1动力学系统的状态空间模型………………1175.1.1状态空间方程的一般形式1175.1.2化高阶微分方程为状态方程——不含输入导数情况……………………1185.1.3线性多自由度振动系统的状态空间模型1215.2微分方程模型与状态空间的关系1235.2.1微分方程模型与状态空间模型特征对的关系……………………1235.2.2系统含有输入导数的状态空间模型1235.3状态空间的相似变换…………………………………………1295.3.1一般情况…1295.3.2特殊情况(可控标准型的情况)…………………1305.4系统的状态空间模型与传递函数模型之间的转换……1315.4.1从状态空间模型转换为传递函数模型……………………1315.4.2模型转换 Matlab函数1325.4.3传递函数模型转换为状态空间模型的直接方法∴……1345.5传递函数模型转换为状态空间模型的串并联法…1355.5.1并联模型法,.,,,,,,,,,,,,,,,1355.5.2串联模型法1385.6状态空间仿真模型建立………………………………………1415.6.1非线性时变系统1415.6.2非线性定常系统……………………………………1425.6.3线性时变系统∴…………1425.6.4线性定常系统……………425.7关于混合系统仿真……………………144习题…∴∴…145第6章连续系统的相似离散法1486.1线性连续系统相似离散法…1486.1.1连续系统状态方程的精确解…486.1.2零阶保持器下状态方程的离散化…1496.1.3一阶保持器下状态方程的离散….,.,,.,,,,.,,,,,·,,,,,,,,,1506.1.4离散系统仿真模块1516.2状态转移矩阵…1526.2.1状态转移矩阵的特性1526.2.2求转移矩阵的几种方法…1536.3离散化系统的传递函数模型…………………1546.3.1零阶保持器的传递函数1546.3.2一阶保持器的传递函数…1546.3.3离散系统的传递函数模型1566.4线性时变系统状态方程的离散化………………………1586.4.1线性时变状态方程的解……………………1586.4.2线性时变系统状态方程离散化…………………1596.4.3近似离散化1596.5离散系统仿真模型建立…………………………1636.5.1有关离散系统 Matlab函数的应用1636.5.2状态方程的离散——基于单位延迟的状态空间仿真模型1666.5.3利用离散传递函数模块的 Simulink仿真模型……∴1686.5.4使用离散状态空间模块 Simulink仿真模型168习题……………………………………………………………………………170第7章机电模拟系统………………………………………1737.1电学基本元件和基本定律1737.1.1电学基本元件……1737.1.2简单电路动态方程1747.1.3电气系统的数学模型建立…………………………………………1747.2无源滤波器……………1787.2.1滤波器基本类型1787.2.2无源RC滤波器……………………………………………………………………………1857.2.3无源RLC滤波器………………………………………………………1867.3机电相似系统1867.3.1力一电压相似……………………………………………………………………187Ⅶ7.3.2力一电流相似…1887.4机电耦合系统的数学建模∴…1897.5运算放大器系统的数学建模∴……………190习题………196第8章系统瞬态响应分析………2008.1典型状态和典型激励的瞬态响应…2008.1.1系统响应种类…2008.1.2常见的几种典型外激励……………2018.2一阶系统的瞬态响应分析·……………………………………………………2028.2.1系统在零输入响应……………………2028.2.2系统零状态响应…………………………………………………2028.2.3标准一阶系统的单位阶跃响应特性……………2048.3二阶系统瞬态响应分析…………………………2068.3.1标准二阶系统的单位脉冲响应………………2068.3.2欠阻尼标准二阶系统的阶跃响应2098.3.3欠阻尼标准二阶系统性能指标…2108.3.4非标准欠阻尼标准二阶系统性能指标…∴…2148.3.5欠阻尼二阶系统的单位斜坡响应………………………………2178.3.6过阻尼二阶系统的单位阶跃响应………………………2188.4 Matlab/ Simulink仿真…………………………………………………………2208.5高阶系统的响应2218.5.1高阶系统的传递函数……………………………………2218.5.2高阶系统的瞬态响应222习题……………………………………………………………………………222第9章动力学系统频域分析方法………………………………………………2269.1概述……………………………………………2269.2频率响应函数…2269.2.1谐和激励下系统的响应函数…………………………………………2269.2.2系统的传递函数与系统的频率响应函数………………………2279.2.3系统频率响应特性曲线(频响曲线)……………………2309.3单位脉冲函数与频率响应函数.·······:··:.·.·········;···;··:.:··2329.3.1单位脉冲响应函数(权函数)…………………2329.3.2单位脉冲函数与频率响应函数…………………………………2339.3.3标准二阶系统的频率响应特性……·.···········2359.4频率响应分析法仿真……23794.1连续系统频率响应特性:,···:,,2379.4.2线性多自由度系统的频域分析2459.4.3快速傅里叶变换与仿真……………………………246Ⅷ9.5频率响应特性在振动系统参数识别中的应用……2489.5.1幅频、相频曲线识别法…………2499.5.2实频、虚频曲线识别法2509.5.3导纳圆的参数识别法…………………252习题………………254第10章动力学系统控制基础…25610.1动力学控制的基本概念…25610.2PID控制系统…………………25710.2.1PID工作简介…25710.2.2PID的数学模型…………25810.2.3PID控制系统的响应分析…25910.3状态反馈控制系统··..···.···:········.···········:······…………27010.4最优控制………………………27410.4.1固定端点的问题最优控制…………………………………27510.4.2始端时刻固定、末值状态自由情况下的最优控制27610.5线性系统的二次型最优设计…………279习题………………………………………………………………………285附录………………………………………………………288附录1 Simulink仿真系统常用模块库288附录2典型函数的拉普拉斯变换和Z变换…………………………………………290附录3 Matlab/ Simulink部分功能设置…………29参考文献…292Ⅸ
- 2020-12-04下载
- 积分:1
泛函分析及其在自动控制中的应用
泛函分析及其在自动控制中的应用,韩崇昭,1991控制理论所硏究的闩越,可以概括为系统分析、系统踪合建模和优化。系统分包括系统的稳定性分析能控能观性分析、鲁棒性分析等,主要是分用以描述系统行为的算子的特生。传统的分析方法是实用的但只限于某些特定的系纯类型.例如传统的枫域分忻法只阳于讨论单输入单输出约线性定常橤统,而泛函分析所提供的分析方法,有可能村包括多输入多输出线性时变系统、分和参觐线性系统,以及某忠类型的柞线性系统进行统…的处理,从而获得更加一般的论。系统的综合包插挖湖器和补偿器的设计等,使系流得以镇定或获得某种性能,这是分析的逆河题。传统的综合屴沄礻仪费时费事丶而且解决问题的范園比较狭窄:现代的综合方法倾向干构造能用计算机灾现的某些算法迭代算法或递推算法的收敏性分析,以及闭环控制的稳定性分衔等,只有借助泛图分析所提供的工其,才有可能使问题得以解决系统建撓和系统的最优控制,一般是在某些约束条件下,对某个泛函拈标进行优化的问题,这更是泛函分析研究惹围闪的问题绕上所泛函分析已渗逶到控制理论和系统科学的各个分支。“饿夯千里日,更上一层楼”,控制斑论研究者只冇掌握泛函分析这…工具,才有可能…览当令研究潮流中“群峰竞秀,万水争流"“的局面第二章代数基础鉴于工科人学亩动控制类专业研究生#不具备系统的拙象代数的知识,而泛函分析这门课程又经带涉及抽象代数的某些基本慨念,所以首先在本章对必要的代数基础知识进行简要介绍,作为学习泛函分枥的预备知识。§2.]集合与映射2.1.1集合集合是数学上最基本門概念,难以绐出确切的定义:一般说:所谓集合就是指具有其种属性的事物全体。构成集合的每个事物称为该集合的元蓊。果合也简称集,其元素也简称元集合可用列其历有元素或江明其暑性来表示如A=i1,nz,…,n},A={a:a具有属性P}如果个集合由有限多个元构成,称之为有限集;如果由无限多个元构成,称之为元限集不含汪何元素的集合称为空集记为这,只含一个元的集合称为单点集。用r∈A表示“?是4中的元”或a属于A”;用a表示a不是A中的元"哎“a不属于有两个集合A和B:若A中的所有元到为B中的元,则称为A是B的子集或A蕴含于B哎包含A,记为A≌或BA任集A必是共自身的了集,而空集又是仕意集A的子集若集A是集R的子集:而B中至少有一个元不屑于4,则称A是程的真子集,或B真包含A,记为心二BBA。若集A是集B的了集,且B也是A的子策,闻称集A与集B相等,记为A=B2这个定义也经出用作集合相等的证明方法即任取∈A,得x∈B则推知4≌B;其次任取xB证得∈小则描知B≌A;从而证明A=B。在以后的证明中,我打经常采用某些撰用符号:“"表示“所有的",“彐"表示“存在”,“→“表示“由左面的结论推出岩面的结论”“台”表示左右两面相互推出”以柴合为元素的集合秋为集类。如字={A,B,C}其中的元A,B,C均是集合,是集类。A、B两个集合的所有元素共同构成的集合称为A和B的并集,记为AUB={2:x∈A或x∈开2。1,2桌合A1,A:;…,A的并業定义为U4=A∪AU…∪4一{x:xEA,或x∈A,“x∈A}(2.1.3)A、B两个集合的公共元素构成的集合称为A和B的交集,记为A∩B={x:x∈A且芏.B〔2.I.4集合A1,A2,A的交集定义为门A=A∩A∩“∩A={x:x∈A且x∈12“且xEA}(2.!5如果集合A与集合升没有公共元家,即A∩B=C,则称A与R不相交。属于集A而不属于集B的所有元构成的集合称为A与B的差集记为AB={xgxA旦x2.⊥,5巢合A和F的对称差记为A△B-(APU(4)2.7设U是一↑特定的集合,AS;称EA为A关于U的补集.记为A此时有AUA=UA∩A=2.1蘸2..9对于集类也可以定义并、交运算。设是一个集类其元的并和交分别为U{BB∈}={z:B∈郾,使z∈B2.1.10∩{B:B∈}{x;B∈密,使r∈乃2.1.11)例211设R表示实数集,R=RXR表示实数序对(x,y)的集合,集合A={(x,g)mx;∈R固定}表示欧氏平面R2上y=m直线上的总集;所有这些集合(直线〕构成一个集类x={A:mER在此情况下,集类m的交集为∩ FRA={(0,0)},即R的坐标原点;其并集为∪v∈RA=R2{0,y):|l|>0}即除去坐标纵轴但保留坐标原点整个R2平面。前面绐出的集合运算具有如下性质〔1)幂等律:AA=A,A∩A=4;〔)交换律:AB=BJA,A∩B-B∩A;I)结合律:A∪(BU)=(4∪B)∪C,A∩(B∩C)=(A∩B)∩C;Ⅳ)分配律:4U(B∩C)=(AUB∩(AUc),A(BUC)-(Anu(A门〔V)恒等律:A∪C=A,A∩U=手AUC=[, A=另外还有一些恒等关系式Ⅵ) de morgan律:A(B)-(AB)∩AC),A(BNC)=(AB)U:A〔〕对偶律:(UAA,(∩A=UⅧ)互补律:AAU,A∩A=8下面只证明(W)和(W)其余留给读者验证W)的证明:∈A(B)∈A,(BC)∈AB且x∈C∈AB且tAr∈AB)n(4C)同理可证第二式。()的证明;*∈(∪A)=∈U,艺点UA+2A且2年A1“且卖Az∈出且x∈A」且x∈术台∈∩4:同理可证第二式渠合A和B的笛卡尔积就是由序对构成的集∈Ab∈2.1,12儿素n和b称为序对(a,b)的分量。如果两个厅对的对应分垤相等则称其相等例如,(x,y)z(n-a,且y一b般情况下,集合的笛卡尔积不叫交换次序,即AX≠xA更一般地,A1A2…**A是一组集合其笛卡尔积定义为A,=1×Azx{「rTI:AE1,2}〔2.I.132.1.2关系由集合A到集合B的一个关系配,就匙的卡尔积A×H的一个集若R=AXB.序对(a)∈R则称n与!有关系P,记为ee,关系RAX称为二元关系,周为此时其中的元由序对钩成更灬般地,若K∈]1,则称其为多元关英。例2.1.2设A一{a,b,c}B={u,b,,d,吧其笛卡尔积为AB=(,,(,b),().(,t).(bq),bb),,(冫(,),出,(qb)(qc),(,d)}如果R表示“相等关系”,则R=Yb),(c,E)}二A×F关系R4XB的定义城是A的子集,即1oia∈A:b∈,使得2..11〕其值址是郾的子集,即ange=,∈B:n∈A使得ah15倒21.3设R={(x,y):,∈R,g-2}表示一个关系,显然它是RXR的一个子粜,即面R上一条抛物线的点集。其定义域domR=R,其值域 rangeR=1={y:y:0),即l半实轴设A是灬个集台,三A×A4是某种关系下给出几种特殊二元关系的定义I)自返关系:若aA→(a,a∈R;(I)对称关系:若(,b)∈B→(b,a)EB;Ⅲ)传递关系:a,b),b)∈R=n,)∈R;Ⅳ)反对称关系:若(,b),h2n)∈→m=b;甚于以上基本的二元关系还可以定义〔V)偏序关系;若是返、传递和反对称的关系:此时称A为出R规定的编序集合;〔)全序关系:若是偏关系,H任意a,b∈A,要么(+b)∈B,要么(b,a);W)等价关系:若F是肖返、传递和对称的关系。例2.1.4设R-(,9):x,∈R,≤匚R表示平面R上包括x-9直线在内的左上半平面,它是一个二元关系:国对()!x∈R,(x,x)∈,所以R是自返的;(2)任意(,y),(,=)∈R即有“到队而有x,即:x,∈R,所以是传的(3){意(x,y),(y,)∈R,具有y,≤x,从而z=y,新以R是反称的;综合起来即证刚R是个懶序关系(即关系“≤’,或者总按偏序关系“≤”規定的R是一个偏序集合。其实,R还是一个全序关系,因为除!R编序关系外,对仁意r∈R.要么x≤,费么ysx即要么(x,y)∈界,要么(y,x)∈P例2.15设A={a,b;,d,g,R={(),(b,)仂,b),(c,q),(+b),e,c)(c,),(d,n),,n),(,4},(,b),(e,e)}4A,它是A上的一个二元关系。因为(1)!x∈A(r,x)∈R,即E是白讴的;(2)意(x,y)(y,z)∈萨即有(x,∈F即R是传递的;〔3若(x,y)(g,x)∈R,即有xy,即是反对称的;所以k是一个偏序关系。但是,为e,∈A,(e)R且(e,c)B,所以不是灬个全序关系若R实A×A是A上的一个偏序关系,则I)A的任一子集B按关系R仍是一个偏详集合,即B(BXB》∩是B上的偏序关系若B按关系R还是一个全序集合,则称其为A的个全序子集。[)设b∈4,若对x∈A均有r助,称b为A的末位元系()设a∈A若对x∈A均有Rx,称a为A的初位元Ⅳ)若b∈A对任意t∈A且l则称b为A的一个最大元素。若m∈A对任意x∈A且→=a则称a为A的一个最小元亲例2.L.6设A={7,4,23,2,5)为一有限整数粜,定义A上的一个二元关系R_{∈A,x≤y2),(4,4),(415),〔47),(4,12),(2,12)(3,3),(3,4),(35),(3;73,12),(2,2),(2,3)、(2,2(2,5),(2,7,(2,12),(5,5),(5,7〕,(苏,2)}显然,孩关系R(即“≤”),A是…个全序集合同时是A的韧位元,也是最小元:12是A的末位兀,也是最大元例2.L.725的关系星可以用图2.1.表示,共中a,即存在有向连线由x到yA的仟一子集加D一{e1b,a}按关系B仍是个偏序集合,而且此时B还是A的一个全序子集在A中按关系R,是末位元也是量大元和c均为A献最小元但A无初位元,月为e和c不存在关系R讲而,假定R是集合A上时一个编序关系,对于B4,则有M)r∈A称为由规定的B的一个上界,当且仅当图2.1.1关系R对∈B有xR。如果a是由R规定的P的一个上界,而对于由R规定的B的任意其它上界p均有aR,则称a为B的最小上界,或上确界.记为m=Lp3()a∈A称为由R规定的的一个下界,当且仅当对yr∈B有aR。如果a是由R规定的B的一个下界,而对于由R规定的B的任意其它下界A,均有山B则称t为B的最大下界或下确界,记为=infB。例2.1.8设A={z∈R:01二R为一闭区何,而B一{x∈R:.2
- 2021-05-07下载
- 积分:1