Simulink高级技术.ppt
《Simulink高级技术.ppt》由会员分享,可在线阅读,更多相关《Simulink高级技术.ppt(126页珍藏版)》请在三一办公上搜索。
1、Simulink高级技术,1 Simulink求解器概念2 系统过零的概念与解决方案3 系统代数环的概念与解决方案4 高级积分器 5 仿真参数设置:高级选项与诊断选项,1 Simulink求解器概念,1.1 离散求解器 第3章中简单介绍了动态系统的模型及其描述,其中指出,离散系统的动态行为一般可以由差分方程描述。众所周知,离散系统的输入与输出仅在离散的时刻上取值,系统状态每隔固定的时间才更新一次;而Simulink对离散系统的仿真核心是对离散系统差分方程的求解。,在对纯粹的离散系统进行仿真时,需要选择离散求解器对其进行求解。用户只需选择Simulink仿真参数设置对话框中的求解器选项卡中的di
2、screte(no continuous states)选项,即没有连续状态的离散求解器,便可以对离散系统进行精确的求解与仿真。,1.2 连续求解器 与离散系统不同,连续系统具有连续的输入与输出,并且系统中一般都存在着连续的状态变量。连续系统中存在的状态变量往往是系统中某些信号的微分或积分,因此连续系统一般由微分方程或与之等价的其它方式进行描述。这就决定了使用数字计算机不可能得到连续系统的精确解,而只能得到系统的数字解(即近似解)。,采用不同的连续求解器会对连续系统的仿真结果与仿真速度产生不同的影响,但一般不会对系统的性能分析产生较大的影响,因为用户可以设置具有一定的误差范围的连续求解器进行相
3、应的控制。离散求解器与连续求解器设置的不同之处如图1所示。,图1 离散求解器与连续求解器设置的比较,由于连续系统状态变量不能够被精确地计算出来,因而积分的误差值同样也是一个近似值。通常,连续求解器采用两个不同阶次的近似方法进行积分,然后计算它们之间的积分差值作为积分误差。连续求解器积分误差的计算如图2所示。,图2 连续求解器积分误差计算,图2中h为积分步长。注意,此图以最简单的多边形积分近似算法为例说明积分误差的计算,在实际中具体的方法视连续求解器的不同而不同。如果积分误差满足绝对误差或相对误差,则仿真继续进行;如果不满足,则求解器尝试一个更小的步长,并重复这个过程。当然,连续求解器在选择更小
4、步长时采用的方法也不尽相同。如果误差上限值的选择或连续求解器的选择不适合待求解的连续系统,则仿真步长有可能会变得非常小,使仿真速度变得非常慢。(用户需要注意这一点。),混合系统仿真时连续状态求解与离散状态求解的协调如图3所示。其中h为初始步长,由于在时刻t与th之间系统存在着离散状态的更新,因而连续变步长求解器将会减小步长至,之后再计算积分误差以控制求解。如果求解误差满足误差范围,则进行下一步仿真,否则缩小时间间隔,重复此过程进行求解仿真。,图3 连续状态求解与离散状态求解的协调示意图,2 系统过零的概念与解决方案,1节中对Simulink的求解器进行了较为深入的介绍。Simulink求解器固
5、然是系统仿真的核心,但Simulink对动态系统求解仿真的控制流程也是非常关键的。Simulink对系统仿真的控制是通过系统模型与求解器之间建立对话的方式进行的:Simulink将系统模型、模块参数与系统方程传递给Simulink的求解器,而求解器将计算出的系统状态与仿真时间通过Simulink环境传递给系统模型本身,通过这样的交互作用方式来完成动态系统的仿真。,系统的时间分类,连续时间系统:状态随时间连续变化,离散事件系统:状态的变化在随机的离散时间点上发生,混杂系统(混合动态系统):混杂上述两种系统状态行为,系统,离散时间系统(模型):时间离散,状态离散,采样系统:对象连续,控制离散,系统
6、的数学描述,连续时间系统:微分方程、传递函数、状态方程,离散事件系统:事件调度、活动扫描、进程交互,混杂系统(混合动态系统):上述方法的综合,系统,离散时间系统:差分方程、状态方程,采样系统:差分方程、状态方程,系统的数学描述,确定型数学描述:状态随时间连续变化,随机型数学描述系统,系统,连续时间模型:微分方程、脉冲响应函数 传递函数、状态空间模型,离散时间模型:差分方程、脉冲传递函数 权系列、离散状态空间模型,常见的数学模型,第三章 连续系统建模方法,常用数学模型:微分方程、传递函数和状态空间模型。,各种模型的特点:,微分方程:线性、非线性系统。时域方法;传递函数:仅适用于线性系统;初始条件
7、为零;用代数的方法或图解分析进行系统分析与综合;频域方法;状态方程:系统的完整描述;化高阶为低阶;降低复杂性 便于处理多变量系统 时域方法;,微分方程,传递函数,状态方程,时域方法,频域方法,不完全描述,完全描述,线性系统,线性系统 非线性系统,3.1 微分方程的机理建模方法,建立微分方程的方法,1)根据规律列方程2)微元分析方法3)模拟近似法,在生物、经济等学科的实际问题中,许多现象的规律性不很清楚,即使有所了解也是极其复杂的,建模时在不同的假设下去模拟实际的现象,建立能近似反映问题的微分方程,然后从数学上求解或分析所建方程及其解的性质,再去同实际情况对比,检验此模型能否刻画、模拟某些实际现
8、象。,建模示例,图1 机械平移系统,图2 机械传动系统,图3 R-L-C电路,f,相似系统,3.2 非线性系统模型的线性化,非线性系统:含有一个或多个非线性元件的系统。,非线性元件:y=f(x),若f(x)是直线,则该元件是线性元 件,反之,则为非线性元件。,典型的非线性特性:死区、饱和、间隙、继电器特性,非线性系统,一些系统作为线性系统来分析,这是由于:系统的非线性不明显,可近似为线性系统。某些系统的非线性特性虽然较明显,但在某些条件下,可进行线性化处理,作为线性系统来分析。这类系统统称为非本质非线性系统。但当系统的非线性特征明显且不能进行线性化处理时,就必须采用非线性系统理论来分析。这类非
9、线性称为本质非线性。,关于线性与非线性的认识,非线性是大自然的普遍特性;线性系统其实只是对少数简单非线性系统的一种理论近似,非线性才是世界的魂魄。生态学和混沌学家罗伯特梅(Robert,May)认为,目前全世界标准的科学教育,向人们灌输的是关于世界图景的偏见和歪曲的印象。不管线性的数学获得了多大的成功,都只能给学生一个关于实际大自然的普遍存在的非线性事实的失真形象。伊恩斯图尔特感叹地说:“如果你断定,只有线性方程才值得研究,那无异于自我禁锢。你的课本充满了线性分析的成功,它的失败埋藏得如此之深,以致连坟墓都看不见,坟墓的存在也没人注意。如同18世纪笃信钟表世界一样,20世纪中叶则恪守线性世界。
10、”,非线性系统模型的处理,等效线性化方法:描述函数法(谐波线性化)、分段线性化法、小参数法等都属于这种方法。小参数法(摄动方法):把系统视为理想模型的参数或结构作了微小扰动的结果来研究其运动过程的数学方法。直接分析方法:李雅普诺夫第二方法(李雅普诺夫稳定性理论)双线性系统理论:对于特殊的双线性系统;流形上的控制理论:以微分几何为主要数学工具的一种分析方法,可用以研究非线性系统的某些全局和局部性质。反馈线性化:代数变换将一个非线性系统的运动特性(全部或部分地)变换成线性的动态特性。,小偏差线性化/局部线性化,所谓线性化,就是在一定的条件下作某种近似,或者缩小一些工作范围,而将非线性微分方程近似地
11、作为线性微分方程来处理。由级数理论可知,若非线性函数在给定区域内存在各阶导数,便可在工作点的邻域将非线性函数展开为泰勒级数。当偏差范围很小时,可略去2次以上高次项,从而得到只包含偏差一次项的线性化方程式,实现函数的线性化。这种线性化方法称为小偏差线性化或局部线性化。,非线性数学模型的线性化,小偏差线性化的这种近似,对大多数控制系统来说都是可行的。1)控制系统在通常情况下,都有一个正常、稳定的工作状态,称为平衡工作点。例如,恒温控制系统的正常工作状态是输入、输出为常值(输出为被控温度,输入为期望值)。2)当系统的输入或输出相对于正常工作状态发生微小偏差时,系统会立即进行控制调节,力图去消除此偏差
12、,因此可以看出,这种偏差是“小偏差”,不会很大。,非线性数学模型的线性化,对 在 附近进行泰勒展开:,条件:,忽略高阶项:,增量方程,建立系统线性化数学模型的步骤,首先确定系统处于正常工作状态(平衡工作点)时各组成元件的工作点,然后列出各组成元件在工作点附近的增量方程,最后消去中间变量,得到系统以增量表示的线性化微分方程。增量方程的数学含义就是将参考坐标的原点移到系统或元件的平衡工作点上,对于实际系统就是以正常工作状态为研究系统运动的起始点,这时系统所有的初始条件均为零。,Matlab是如何解微分方程的?答案:线性微分方程和低阶特殊微分方程一般通过解析解的方法求解,但一般的非线性微分方程是没有
13、解析解的,需要使用数值解的方式求解。对于多元非线性常微分方程来说,Euler算法是最直观的一类求解算法,精度较高的算法是Runge-Kutta法、Adams法等。但是至今尚没有一种算法可以有效地适用于所有场合。,在用微分方程描述的一个变化过程中,若往往又包含着多个相互作用但变化速度相差十分悬殊的子过程,这样一类过程就认为具有“刚性”。描述这类过程的微分方程初值问题称为“刚性问题”。,例1:单容水箱液位系统设供水阀前压力稳定,出水阀其流量公式为:,其中,为流量系数。,试列写液位波动的微分方程,,并对其进行线性化处理。,非线性系统的线性化多变量系统局部线性化,将该方程式在工作点,附近展开成泰勒级数
14、,并忽略二,若输出变量y与输入变量x1,x2有非线性关系,即,阶和高阶导数项,便可得到y的线性方程为:,写成增量形式:,称为系统的静态/稳态方程。,滑阀系统压力-流量特性线性化,滑阀系统压力-流量特性线性化,某滑阀的压力-流量特性如图所示,其非线性方程为:,设系统的静态工作点为A点。,为阀芯位移;,为负载压差;,为负载流量;,这里:,滑阀系统压力-流量特性线性化,将式,在工作点,处泰勒展开,,上式与静态方程相减:,写作:,图示:负载流量特性的线性化,应用范围:,非线性系统的线性化局部线性化,线性化是相对某一工作点进行的。工作点不同,得到的线性化微分方程的系数也不同。若使线性化具有足够的精度,调
15、节过程中变量偏离工作点的偏差信号应足够小。线性化不适用于含本质非线性的系统。,线性化是以直线代替曲线,略去了泰勒级数展开式中的二阶以上无穷小项,这是一种近似处理。如果系统输入量工作在较大范围内,所建立的线性化数学模型势必会带来较大的误差。所以,非线性数学模型线性化是有条件的。线性化以后得到的微分方程是增量微分方程,只是为书写方便,略去了增量符号。,连续系统数字仿真,连续系统的仿真,主要的计算工作是求解一阶微分方程。,阶段,技术,建模技术,仿真算法,软件技术,模型建立,模型变换,仿真试验,连续系统数字仿真 一般过程,连续系统数字仿真方法,基于数值积分方法,面向方程,连续系统数字仿真方法,基于离散
16、化方法,面向方程结构图:,微分方程,传递函数,状态模型,传递函数+非线性环节,面向方程,面向方程结构图:,传递函数+非线性环节,传递函数,状态模型,连续系统数字仿真方法基于数值积分方法,连续系统的仿真,主要的计算工作是求解一阶微分方程。,数值解法:寻求解y(t)在一系列离散点,上的近似解,(即数值解),数值积分方法,连续变量问题,离散的差分方程的初值问题,常用的几种数值积分方法,单步法,欧拉法:精度差(解释概念),数值积分方法,多步法:,改进的欧拉法:精度提高,速度快,Adams法,Adams隐式公式,龙格-库塔法:精度高,常用,Adams显式公式,预估-校正法,计算量小,单步法(欧拉法),1
17、.由泰勒展开式得到:,2.矩形近似解法:,3.切线近似解法:,特点:方法简单、计算量小,但精度差,泰勒级数方法,优点:,方法概念的简单性,并且具有非常高精度的潜力。,缺点:,依赖于给定的微分方程的反复求导。,几何解释,图a 矩形近似及其误差,图b 切线近似及其误差,改进的欧拉法,预估:,曲边梯形面积:,直边梯形面积:,h较小时,,简记为:,改进的欧拉法,龙格-库塔(Runge-kutta)法,泰勒展开方法:,思路:利用几个点上的y(t)的一阶导数值的线性组合来近似代替y(t)在某一个点的高阶导数值。然后再利用泰勒展开式确定线性组合的各加权系数。这样既避免了计算高阶倒数,又提高了数值积分的精度。
18、,考虑一阶微分方程:,这里:,将 写成线性组合形式:,一阶RK公式,欧拉公式,改进欧拉公式,特点:,1)单步法;2)同样精度,RK4的h比RK2的h大十倍,而计算量仅大一倍,常用RK4,多变量形式,例:已知系统方程为:,解:原系统方程可化为下列方程组,,多步法(Adams),预估校正法,常用的几种数值积分方法,单步法,欧拉法:精度差,数值积分方法,多步法:,改进的欧拉法:精度提高,速度快,Adams法,Adams隐式公式,龙格-库塔法:精度高,常用,Adams显式公式,预估-校正法,计算量小,误差、收敛性与稳定性分析,1.误差:截断误差与舍入误差,2.收敛性:,3.稳定性:,采用数值积分方法进
19、行仿真时常常会出现这样的现象:本来是稳定的系统,仿真结果却得出不稳定的结论。这种现象是由于计算步长太大造成的,因为,当步长h选得过大时,计算误差较大,数值积分方法会使误差传播,以致引起计算不稳定。因稳定性与步长有密切关系,故在处理实际问题时,可用两个显著不同的步长进行计算,若所得结果基本相同,则一般是稳定的;反之,则很可能不稳定;,数值积分方法的选择,1.精度要求:,a.步长确定时,算法阶次越高,截断误差越小。,b.阶次确定时,多步法单步法,隐式显式,2.计算速度:,同步长,RK4比4阶Adams预估校正法慢。,3.数值解的稳定性:,小于4阶时,隐式AdamsRK4显式Adams,步长控制,1
20、.经验公式,2.变步长计算公式,连续系统数字仿真方法,基于数值积分方法,面向方程,连续系统数字仿真方法,基于离散化方法,面向方程结构图:,微分方程,传递函数,状态模型,传递函数+非线性环节,面向方程,面向方程结构图:,传递函数+非线性环节,传递函数,状态模型,离散相似法,基本原理:,将一个连续系统进行离散化处理,然后求得等价的离散模型。,传递函数,状态方程,Z域离散相似模型,时域离散相似模型,传递函数-Z域离散相似模型,状态方程-时域离散相似模型,图斯汀置换法:传递函数-Z域离散相似模型,连续系统数字仿真方法,基于数值积分方法,面向方程,连续系统数字仿真方法,基于离散化方法,面向方程结构图:,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Simulink 高级 技术
链接地址:https://www.31ppt.com/p-5448268.html