基本模块介绍PPT学习教案

1、会计学1基本模块介绍基本模块介绍2021年10月15日2Simulink 的调试技术Simulink 作为高性能的系统设计、仿真与分析平台，给用户提供了强大的模型调试工具。通过这个调试工具，用户可以对动态系统的模型进行调试，已发现其中可能存在的问题，进行修改。1.Simulink图形调试器启动 使用菜单Tools下的Simulinkdebugger 命令或使用调试器按钮可以启动调试器第1页/共39页2021年10月15日32. 调试器的操作设置与功能1）Simulink调试器工具栏命令功能：执行至下一个模块； 执行至下一个时间步；开始调试或继续运行仿真至下一个断点； 终止调试过程；在当前选择的

2、模块之前设置断点；执行时显示当前选择模块的输入输出； ? 调试帮助；第2页/共39页2021年10月15日42）断点显示及断点条件设置在断点显示框中了解到断点位置、断点模块的输入输出。提供五种断点条件设置：Zero crossing ：在系统发生过零处设置断点Step size limited by state:在仿真步长受到状态限制处设置断点Minor time steps : 在最小仿真步长出现处设置断点Nan values:在系统中出现无穷小之处设置断点Break at time : 在指定的仿真时刻处设置断点第3页/共39页2021年10月15日53）调试器输出窗口Outputs :

3、输出调试结果，如调试时刻、调试的模块及模块输入输出Execution order: 输出调试过程中的各模块的执行顺序Status:输出调试状态,如当前仿真时间、缺省调试命令、调试断点设置以及断点数等状态信息系统调试举例system_debugger.mdl第4页/共39页4.1 连续系统模块 大多数物理系统都可以用微分方程进行描述，因此都可以用连续系统模拟。最简单的模型是线性模型和定常模型。其中，x 为系统的广义坐标列向量，M 为质量矩阵，C 为阻尼矩阵，K 为刚度矩阵，P(t)为外部激励列向量。)()()()(ttttPKxxCxM 例如，振动理论中的动力学方程： 在 Simulink中，用

4、来模拟连续系统的基本模块有四个：增益模块，求和模块，微分模块，积分模块。任何可以用线性微分方程描述的系统都可以用这四个基本模块进行模拟。除了这四个基本模块，传递函数模块也经常用来模拟物理系统和控制器。第5页/共39页1. 增益模块作用：使增益模块的输入信号乘以一个常数，并输出。增益模块用代数表达式可表示为： y(t) = k x(t)k增益模块x(t)y(t)kGainSimulink 增义模块图可用简图表示如下：注意：y(t)、x(t)、k 可以为标量、向量或矩阵。看如下算例。5Gain10Display2Constant标 量 乘 积：来自Sinks库来自Sources库第6页/共39页标

5、量和向量的乘积向量和标量的乘积向量和矩阵的乘积矩阵和向量的乘积1;2* uGain24Display2Constantu *2Gain24Display(1 2)Constant1 23 4u*Gain710Display(1 2)Constant12u*Gain511Display1 23 4Constant4221242221107432121115214321第7页/共39页2. 求和模块作用：对两个或多个信号进行求和运算。求和模块必须至少有一个输入而仅有一个输出。输入的正负号的数目由双击模块进入编辑栏进行设定。求和模块用代数表达式可表示为： c = a + bacbabc求和模块有两种

6、形状：圆形和方形。求和模块不但可以进行标量求和运算，也可以进行向量或矩阵求和运算，但是标量或矩阵的维数必须相等。第8页/共39页5Display3Constant12Constant35Display2;3Constant11;2Constant3759Display2 34 5Constant11 23 4Constant标量求和向量求和矩阵求和533221532975354324321第9页/共39页3. 微分模块作用：计算输入对时间的变化率。Simulink 的微分模块如图所示：微分模块代表如下微分方程：ydtdxdu/dtDerivative考虑对正弦信号 sin (t) 的微分：)c

7、os()sin(tdttdSimulink 模型框图和仿真结果如下页图形所示。第10页/共39页Sine WaveScopedu/dtDerivative012345678910-1-0.8-0.6-0.4-0.200.20.40.60.81)cos()sin(tdttd第11页/共39页4. 积分模块作用：计算输入信号从起始时间到当前时刻对时间的积分。Simulink 的积分模块如图所示：积分模块代表如下积分方程：ttdxtyty0)()()(01sIntegrator例如：Sine WaveScope1sIntegrator012345678910-1-0.8-0.6-0.4-0.200.

8、20.40.60.81初值设置为 -1)cos()cos()sin()(tCtdttty注：y(t) 的初值假设为-1第12页/共39页5. 简单物理模型 利用前面所介绍的这些模块可以模拟由线性微分工程描述的任何物理模型。例如，考虑如下所示的简单的小车系统运动。Fx忽略摩擦力，运动微分方程为：mFx 可用模块图表示为:x Fm1x Fm1x xs1s1加入两个积分模块，第一个模块用来计算速度，第二个模块用来计算位移：假定 F=sin(t) 为正弦激励，m=0.5。求020s 区间内的系统位移响应曲线。模型框图如图所示。ScopeSine Wavex Fm1x xs1s1第13页/共39页024

9、681012141618200510152025303540第14页/共39页 若要求同时输出位移和速度，则模型框图为：ScopeSine Wavex Fm1x xs1s1024681012141618200510152025303540第15页/共39页 若要求同时输出位移、速度和加速度，则模型框图为：ScopeSine Wavex Fm1x xs1s102468101214161820-50510152025303540第16页/共39页0)0(, 2 . 0)0(, 0 xxkxxcxm xkc05. 0, 1, 5 . 0ckm051015202530-0.2-0.15-0.1-0.0

10、500.050.10.150.2x x xScope1sIntegrator11sIntegrator0.05/0.5Gain11/0.5Gain初值设置为0.2xmkxmcx 动力方程变换为：第17页/共39页2021年10月15日19线性连续系统的传递函数模型与零极点模型采用连续信号的拉氏变换来实现对于连续信号 ，其拉氏变换定义为：)(tudtetusUst)()(系统的输入时间 ，这时 ：0t0)()(dtetusUst拉氏变换具有两个性质：1）线性性；2） 的拉氏变换为 的拉氏变换为)(tu )(ssU)(tu )(2sUs6. 传递函数模块第18页/共39页2021年10月15日20

11、对于如下的的线性连续系统：kxxcxmtu )(0m同时对等式的两边进行拉氏变换：)()()()(2skXscsXsXmssU将其化为分式的形式: kcsmssUsX21)()(这便是系统的传递函数模型一般可以写成如下传递函数的形式：212010)()(dsdsdnsnsUsY第19页/共39页2021年10月15日21将其进行一定的等价变换，可以得出线性连续系统的零极点模型。)()()(211pspszsksUsY为线性连续系统的零点，为系统的极点为系统的增益1z21, ppk第20页/共39页2021年10月15日22线性连续系统的Simulink描述一般来说，在Simulink中对线性连

12、续系统描述方式有三种：传递函数模块：传递函数表示为 num=n0，n1；den=d0,d1,d2num表示传递函数的分子系数向量den表示分母系数向量零极点模块 ：零极点模型表示为 gain=k;zeros=z1;poles=p1,p2;gain表示增益， zeros表示系统零点， poles表示系统极点状态空间模块：则在Simulink 中直接输入变换矩阵A,B,C,D212010)()(dsdsdnsnsUsY)()()(211pspszsksUsY第21页/共39页2021年10月15日23Fkxxcxm xkcF(t)例如：考虑弹簧质量阻尼系统动力学方程可写为：)()()()(2sFs

13、kXscsXsXms则系统的传递函数可表示为：m=1,c=2,k=2,F=sin(t)kcsmssFsX21)()(作 Laplace 变换，并忽略初始条件，有：第22页/共39页2021年10月15日24m=1,c=2,k=2,F=sin(t)Transfer Fcn模块kcsmssFsX21)()(ParametersNumerator:0 1Denominator:1 2 2双击from Cotinuous 模块库见文件TransferFcn.mdl212010)()(dsdsdnsnsUsY第23页/共39页x m1x xs1s1+mkmc 使用基本模块元素所建立起的系统模型xmkxm

14、cmFx liner_continue_system.mdl第24页/共39页FxScopeSine Wavex Fm1x xs1s1以前面受正弦外力作用的质量块为例：Simulink 模型仿真框图为：拉氏变换： F(s)=ms2X(s)，传递函数分子输入：0 1分母输入：0.5, 0, 0F=sin(t)m=0.5mFx 21)()(mssFsX第25页/共39页2021年10月15日27024681012141618200510152025303540若采用传递函数进行仿真，Simulink 模型框图为：结果与前面的结果一致TransferFcn_1.mdl第26页/共39页对传递函数模块

15、的说明： 用传递函数模块对线性定常系统进行动态仿真，能够使得仿真模型形式简单和紧凑。 但是如果需要调用内部变量，如 x 的导数，则传递函数模块无法适用。 另外，如果需要设定初始条件，传递函数同样无法适用。 再者，传递函数模块只适用于单输入单输出系统，即通常所说的单自由度系统，但无法应用于多输入多输出系统，即多自由度系统。第27页/共39页2021年10月15日295.2 状态空间模块 状态空间模块可以起到与传递函数模块相同的作用。所不同的是，状态空间模块允许用户指定初始条件，并且可以共享内部变量。另外，状态空间模块可以用来模拟多输入多输出系统。1. 状态空间的概念 以单自由度为例说明状态空间和

16、状态向量的概念。考虑如下单自由度系统：)()()()(tFtkxtxctxm 其中，x(t) 为系统的广义坐标；m、k、c 分别为单自由度系统的广义质量、广义刚度、阻尼系数；F 为系统的外部输入。这些变量都是标量。第28页/共39页)()()(txtxtYY(t) 为一个2维列向量，称为系统的状态向量。利用状态向量Y(t)，则上述动力方程可转化为：)(0)()(10)()()(111tFmtxtxcmkmtxtxt Y如果令：以 和 构成的空间即为状态空间。)(tx)(tx )()()(tFttBAYY简写为：状态方程)()()()(tFtkxtxctxm 第29页/共39页考虑 n 自由度系

17、统：)()()()(ttttFKXXCXM 其中：121)(,),(),()(nTnRtxtxtxtXn 维广义坐标列向量nnRM质量矩阵nnRC阻尼矩阵nnRK刚度矩阵外部输入列矩阵令：12)()()(nRtttXXY2n 维系统状态列向量)()()(tttBFAYY则系统的状态方程有：nnnnnnR2211CMKMI0A系数矩阵输入矩阵1)(nRtFnnnnR21M0B第30页/共39页2021年10月15日32线性连续系统的状态空间模型Simulink描述x(t)为线性连续系统的状态向量，u(t)为系统的输入向量，y(t)为系统的输出向量.则在Simulink中直接输入变换矩阵A,B,C

18、,D即可)()()()()()(ttttttu ux xy yu ux xx xDCBA第31页/共39页2. 线性单输入输出系统xkc写成状态方程形式，有：前面考虑过的单自由度系统：其中：1 . 02101011cmkmA在状态空间模块中输入：A=0, 1; -2, -0.1， B=0 ；2 C=1 0； 0 1 D=0 ；0initial conditions:0.2 00) 0 (, 2 . 0) 0 (,xxFkxxcxm 1,05. 0, 1, 5 . 0FckmFttBAYY)()(201m0B)()()()()()(ttttttu ux xy yu ux xx xDCBA)(0)

19、()(10)()()(111tFmtxtxcmkmtxtxt Y第32页/共39页A=0, 1; -2, -0.1， B=0 ；2 C=1 0； 0 1 D=0 ；0见linear_single_intput_system.mdl如果取C=1 0；0 1，则可以输出系统的状态响应（包括位移和速度）若取C=1 0；0 0，则示波器显示出位移时程若取C=0 0；0 1，则示波器显示速度时程1F0) 0 (, 2 . 0) 0 (,xxFkxxcxm 第33页/共39页3. 多输入输出系统x1k1c1c2k2x21, 5 . 021mm2, 121kk05. 021 cc0KXXCXM)()()(ttt