第三章 续系统仿真方法学ppt课件

1、第三章第三章 连续系统仿真方法学连续系统仿真方法学本章主要内容本章主要内容n连续系统建模方法连续系统建模方法n模型变换模型变换n连续系统仿真算法连续系统仿真算法n采样控制系统仿真采样控制系统仿真*n分布参数系统仿真分布参数系统仿真*第一节第一节 连续系统建模方法连续系统建模方法n先验知识建模先验知识建模n机理建模方法机理建模方法知识模型知识模型 n常见形式：各种学科的公理、定理、定常见形式：各种学科的公理、定理、定律等律等 n专家系统方法专家系统方法逻辑关系模型逻辑关系模型 n符号、关系式、专家知识库、推理规则符号、关系式、专家知识库、推理规则等等 n模糊系统方法模糊系统方法模糊模型模糊模型

2、n高矮、大小等模糊语言量化成定量的表高矮、大小等模糊语言量化成定量的表示形式，按某种算法得到定量的结果后示形式，按某种算法得到定量的结果后再转换为模糊语言再转换为模糊语言 n以上方法有时也用于离散事件系统建模以上方法有时也用于离散事件系统建模 n系统辨识建模系统辨识建模 n经验方法经验方法 n直接观察数据曲线得出模型方程，如线直接观察数据曲线得出模型方程，如线性系统，一阶对象等性系统，一阶对象等 n表格插值，一种静态建模技术，主要用表格插值，一种静态建模技术，主要用于计算动态方程中的系数于计算动态方程中的系数 n统计建模数理统计的方法）统计建模数理统计的方法） n最小二乘法及其改进形式、极大似

3、然估最小二乘法及其改进形式、极大似然估计法等计法等 n神经网络神经网络 n混合建模方法混合建模方法 n若干种模型形式输出互相补充若干种模型形式输出互相补充n n给定输入后，从机理模型中产生输出，给定输入后，从机理模型中产生输出，与辨识模型的输出按某种方式得到系统与辨识模型的输出按某种方式得到系统输出，反过来可以用输出误差继续修正输出，反过来可以用输出误差继续修正辨识模型辨识模型 第二节第二节 模型变换模型变换n连续系统常用的模型表示形式连续系统常用的模型表示形式n连续时间模型连续时间模型n系统的输入量系统的输入量u(t)，输出量，输出量y(t)及内部及内部状态变量状态变量x(t)均为时间的连续

4、函数均为时间的连续函数n微分方程微分方程n传递函数传递函数n权函数权函数n状态空间表达式状态空间表达式微分方程微分方程n其中其中mnn用古典方法求解时非常复杂，用古典方法求解时非常复杂，n高阶系统通常没有封闭解或解析解高阶系统通常没有封闭解或解析解1111011011nnmmnmmnnmmd ydydyd uduaaa ybbbub udtdtdtdtdt1111011011nnmmnmmnnmmd ydydyd uduaaa ybbbub udtdtdtdtdt传递函数传递函数 n当初始条件为零时，对上述微分方程式作拉氏变换，当初始条件为零时，对上述微分方程式作拉氏变换，可得传递函数形式可得

5、传递函数形式 n求解时可先用部分分式展开求解时可先用部分分式展开 n再进行反变换即得时间解再进行反变换即得时间解 nS域复频域内求解较为简便但对多变量、时变或高域复频域内求解较为简便但对多变量、时变或高阶系统仍求解困难阶系统仍求解困难 11101110( )( )( )mmmmnnnb sbsb sbY sW sU ssasa sa12( )( ) ( )( )( )( )lY sW s U sW sW sW s112( ) ( )( )( )( )ly tLY sy ty ty t权函数权函数 权函数权函数 g(t)指初始条件为指初始条件为0时系统在理想脉冲函时系统在理想脉冲函数数(t)作用

6、下的响应，又称脉冲过渡函数作用下的响应，又称脉冲过渡函数系统对任意输入的响应可由卷积积分公式求出系统对任意输入的响应可由卷积积分公式求出0( )( ) ()ty tug td权函数与传递函数有如下关系：权函数与传递函数有如下关系： ( )( )L g tG s 为为n维状态向量，维状态向量，u为为r维输入向量，维输入向量，y为为m维输出向量维输出向量状态空间表达式状态空间表达式 动态系统的状态是指能完全描述系统行为的最小一组变量动态系统的状态是指能完全描述系统行为的最小一组变量若知道若知道t=t0时刻的初始状态向量时刻的初始状态向量x0及及tt0时的输入时的输入u就能完全就能完全确定系统在确定

7、系统在tt0时刻的行为时刻的行为状态空间表达式由状态方程和输出方程组成状态空间表达式由状态方程和输出方程组成xAxBuyCxDu1,Tnxxxn nAA为参数矩阵或称动态矩阵），为参数矩阵或称动态矩阵），n rBB为输入矩阵，为输入矩阵，m nCC为输出矩阵为输出矩阵m rDD运用矩阵计算方法且借助于计算机很容易运用矩阵计算方法且借助于计算机很容易对状态空间方程求解对状态空间方程求解 为关联矩阵输入和输出直接关联）为关联矩阵输入和输出直接关联）n离散时间模型离散时间模型n系统的输入量，输出量及内部状态量均为时系统的输入量，输出量及内部状态量均为时间的离散函数，即时间序列间的离散函数，即时间序列

8、u(kT),y(kT),x(kT)n差分方程差分方程nZ传递函数传递函数n权序列权序列n离散状态空间模型离散状态空间模型差分方程差分方程nT为采样周期为采样周期 11010() (1) (1) ()() (1) ()()nmmy kn Tay knTa y kTa y kTc u km Tcu kmTc u kTmnZ传递函数传递函数n对差分方程作对差分方程作Z变换，设所有初值为变换，设所有初值为0，则有，则有11101110( )( )( )mmmmnnnb zbzb zbY zW zU zzczc zc权序列权序列n权序列权序列h(k)为对初始条件为为对初始条件为0的系统施加单位脉冲的系统

9、施加单位脉冲序列序列(k)所得到的响应所得到的响应 n系统对于任意输入系统对于任意输入u(k)的响应为一卷积的响应为一卷积 n与与Z传递函数间关系传递函数间关系 0( )( ) ()kiy ku i h ki ( )( )Z h kW z离散状态空间模型离散状态空间模型n连续连续-离散混合模型离散混合模型n如计算机控制系统，对连续对象进行控如计算机控制系统，对连续对象进行控制时，状态量中既有连续的也有离散的制时，状态量中既有连续的也有离散的 (1)( )( )( )( )x kFx kGu ky kx k n连续系统模型之间的变换连续系统模型之间的变换 n微分方程、传递函数、权函数模型描述微分

10、方程、传递函数、权函数模型描述系统的输入与输出关系，称为系统的外系统的输入与输出关系，称为系统的外部模型部模型 n状态方程则称为系统的内部模型状态方程则称为系统的内部模型 n通常在仿真时，需要将系统的各种描述通常在仿真时，需要将系统的各种描述形式转换成内部模型，称为模型结构变形式转换成内部模型，称为模型结构变换换n化微分方程为状态方程化微分方程为状态方程 n化连续状态方程为离散状态方程化连续状态方程为离散状态方程化微分方程为状态方程化微分方程为状态方程设有微分方程设有微分方程 先考虑右边仅含先考虑右边仅含u的形式，令的形式，令 = 则有则有 1111011011nnmmnmmnnmmd ydy

11、dyd uduaaa ybbbub udtdtdtdtdt12,Tnxx xx(1) , ,nTy yy 1223101121nnnnnxxxxxxxa xa xaxu 写成矩阵形式：写成矩阵形式： 输出矩阵写为输出矩阵写为 其中其中 xAxBu101100nnIAaaa001B 100yCxC当右式包含导数项时，状态方程形式为当右式包含导数项时，状态方程形式为 A,C与前相同，与前相同， 其中其中 xAxBuyCxDu10,TnBD1,(0,1,)iin injijjbain现代控制理论中还介绍了其它形式的转换方程，现代控制理论中还介绍了其它形式的转换方程，如能控标准型、能观标准型等如能控标

12、准型、能观标准型等 化连续状态方程为离散状态方程化连续状态方程为离散状态方程 连续状态方程对应的离散状态表达式为连续状态方程对应的离散状态表达式为 T为采样周期或者计算步长为采样周期或者计算步长 为确定为确定(T)和和H(T)，可利用连续状态方程解，可利用连续状态方程解 其中其中Ate为系统的矩阵指数或状态转移矩阵，为系统的矩阵指数或状态转移矩阵，x(0)为初始状态向量为初始状态向量 当采用零阶保持器时，当采用零阶保持器时， 即认为即认为u(t)在每个采样周期内保持常值在每个采样周期内保持常值 u(t)=u(kT)，( kTt(k+1)T ) (1) ( ) ()( ) ()(0,1,)x k

13、TT x kTH T u kTk 0( )(0)( )tAtAtAx te xeeBud(1)(1)(1)00(1)(1)0(1) (0)( )(1)(0)( )(2)(1)(2)(1) ()( )()()( ()()kTA kTA kTAkTAkTAkTAATkTATA kTAkTTATATAtATAx kTexeeBudx kTexeeBudex kTex kTeeBudex kTeeBu kT dtu kTex kTe与得：不变0()( ) ()( ) ()TBu kT dtT x kTH T u kT则有则有 其中其中和和H与与T有关，当有关，当T确定后，确定后，和和H为常值矩阵为常值

14、矩阵 离散化公式的核心在于计算矩阵指数及其积分，离散化公式的核心在于计算矩阵指数及其积分，常用级数展开的算法，即常用级数展开的算法，即 其余离散化的表示形式与连续形式之间的转换其余离散化的表示形式与连续形式之间的转换在计算机控制中介绍在计算机控制中介绍 22112!ATkkeIATA TA Tk第三节第三节 连续系统的仿真算法连续系统的仿真算法算法的基本概念算法的基本概念系统模型系统模型计算机模型：二次建模，算法是核心计算机模型：二次建模，算法是核心问题问题算法：解题方案的准确而完整的描述，一般采用算法：解题方案的准确而完整的描述，一般采用文字、算式以及框图的形式文字、算式以及框图的形式需要关

15、注：需要关注：算法性能分析：误差、收敛性、计算效率等算法性能分析：误差、收敛性、计算效率等算法的比较与选择算法的比较与选择 n浮点数运算浮点数运算 n计算机上进行数值计算时，实数计算机上进行数值计算时，实数x用用t位位十进制浮点数表示：十进制浮点数表示：n其中其中m为为t 位十进制小数，且位十进制小数，且-1m1，c为十进制整数，若为十进制整数，若0.1m1，则称此，则称此浮点数系统为规格化的，浮点数系统为规格化的，t 称该数的精度，称该数的精度，特定的计算机有固定的浮点数精度特定的计算机有固定的浮点数精度 10cxm采用浮点数运算存在的常见问题采用浮点数运算存在的常见问题n舍入误差舍入误差n

16、计算机有一组操作浮点数的指令，用以模拟加、计算机有一组操作浮点数的指令，用以模拟加、减、乘、除运算，但不可能精确。如乘法运算减、乘、除运算，但不可能精确。如乘法运算时，乘积应有时，乘积应有2t位精度，但实际仅能保留位精度，但实际仅能保留t位，位，即存在舍入误差。复杂计算迭代等中舍入即存在舍入误差。复杂计算迭代等中舍入误差的累积可能会影响结果，应在算法分析中误差的累积可能会影响结果，应在算法分析中考虑考虑 n溢出溢出 n计算机对指数计算机对指数c范围有限制，乘、除时可能会范围有限制，乘、除时可能会上溢、下溢，也应进行处理上溢、下溢，也应进行处理 数值稳定性问题数值稳定性问题 若运算过程中计算误差

17、不断增长，称算法为数值不稳定的若运算过程中计算误差不断增长，称算法为数值不稳定的反之则为稳定的反之则为稳定的 例：计算例：计算 由分部积分得递推公式：由分部积分得递推公式： 用用Taylor展开计算：展开计算： 若取若取k=7，并保留，并保留4位小数，可得位小数，可得 截断误差：截断误差： 110(0,1,)nxnIex e dxn1111001,1xnnInIIee dxe 21( 1)( 1)1( 1)2!kek 10.3679e147110.3679108!4Re只考虑初值误差，对只考虑初值误差，对I00.6321，递推计算，结果如表中第一行，递推计算，结果如表中第一行该积分不可能为负值

18、，显然算法有问题该积分不可能为负值，显然算法有问题 分析计算误差：分析计算误差： 满足关系：满足关系： 误差增长迅速误差增长迅速 n01234567890.63210.36790.26420.20740.17040.14800.11200.2160-0.7287.5520.63210.36790.26430.20730.17080.14550.12680.11210.10350.0684nnnEII1nnEnE 0( 1) ( !)nnEn E 如果换一种算法，可以减小误差，考虑积分估计值：如果换一种算法，可以减小误差，考虑积分估计值： 逆向计算：逆向计算： 取取n=9时，时， ，结果如表中第

19、二行，结果如表中第二行 分析计算误差，满足关系：分析计算误差，满足关系： 显然误差一直减小显然误差一直减小 111110001011(min)(max)11xnxnnxxeeex dxIeex dxnn *11(1)nnIIn191 1()0.06842 1010eI11nnEEn 病态问题病态问题 如线性方程组如线性方程组 精确解：精确解： 12123.0004.12715.41(1)1.0001.3745.147(2)xxxx1213.66586.2xx 如果用如果用4位有效数字进行运算：位有效数字进行运算：(2)-(1)/3.000，逐步计算后，可得逐步计算后，可得x2=-5.000 ，

20、误差很大，误差很大当方程特征根相差太大时出现病态问题，也称刚性当方程特征根相差太大时出现病态问题，也称刚性(Stiff)问题问题需要设计有效的算法需要设计有效的算法 二、数值积分法二、数值积分法n实际系统模型多为低阶微分方程形式，求解实际系统模型多为低阶微分方程形式，求解时本质上应用积分运算时本质上应用积分运算n对高阶方程，可先转换为多个一阶方程，因对高阶方程，可先转换为多个一阶方程，因而，最终问题转化为求解一阶微分方程而，最终问题转化为求解一阶微分方程n常见算法：常见算法：nEuler法法 nRunge-Kutta法法(R-K法法)nAdams法法(多步法多步法)Euler法法设有模型方程设

21、有模型方程 ，初始条件，初始条件 欧拉法用欧拉法用tk点切线近似该点附近的曲线点切线近似该点附近的曲线f(t,y)，则有，则有 其中其中是曲线上的点，是曲线上的点，是切线上的点是切线上的点 称为第称为第k 步的计算步长步的计算步长 此类方法称为此类方法称为“微分方程初值问题的数值计算法微分方程初值问题的数值计算法”，也称也称“数值积分法数值积分法” 优点：简单易行优点：简单易行 缺点：缺点：h取得大时，计算速度快，单步误差大；取得大时，计算速度快，单步误差大； h取得小，计算速度慢，且累计误差大取得小，计算速度慢，且累计误差大1()ky t1ky( , )dyf t ydt00( )y ty1

22、11()( ,)()kkkkkkky tyyf tytt1kkkhttRunge-Kutta法法(R-K法法) 二阶形式：二阶形式： 其中其中 迭代公式由迭代公式由Taylor展开并保留展开并保留h2项获得项获得 注意：注意：R-K方法实质是用均差代替导数，方法实质是用均差代替导数， 其中其中k项的加权系数可任选项的加权系数可任选 1112()()2kkkhy tyykk121( ,)(,)kkkkkf tykf th yk h四阶形式固定步长）：四阶形式固定步长）：其中其中 四阶形式在精度和复杂度方面都有较好的表现，最常用四阶形式在精度和复杂度方面都有较好的表现，最常用 Euler法与法与R

23、-K法计算时仅用到前一步的结果，法计算时仅用到前一步的结果，称单步法，已知初值后可自启动称单步法，已知初值后可自启动 111234()(22)6kkkhy tyykkkk1213243( ,)(/2,/2)(/2,/2)(,)kkkkkkkkkf tykf thyk hkf thyk hkf th yk hAdams法法(多步法多步法) Euler法是用矩形公式面积近似定积分，在曲线与法是用矩形公式面积近似定积分，在曲线与矩形的边之间有较大误差；矩形的边之间有较大误差；Adams法考虑用梯形公式面积代替矩形公式法考虑用梯形公式面积代替矩形公式 称二阶隐式称二阶隐式Adams法公式法公式 因公式

24、右端包含未知项，不能直接求解，可以用迭代因公式右端包含未知项，不能直接求解，可以用迭代法求解，设有初值法求解，设有初值 迭代公式：迭代公式： 直至规定的精度直至规定的精度 1111( , ) (,)( ,)()22kktkkkkkkthhf t y dtf tyf tyff11()2kkkkhyyff(0)1ky( )(1)111 (,)( ,)2nnkkkkkkhyyf tyf ty隐式迭代计算步数太多，为此可降低精度，设计显式隐式迭代计算步数太多，为此可降低精度，设计显式Adams法，公式：法，公式： Adams法的统一形式为：法的统一形式为： 算法特点：多步法，不能自启动，隐式方法还需迭

25、代求解算法特点：多步法，不能自启动，隐式方法还需迭代求解实际应用中，先用显式法计算初值，再用隐式法校正一次，实际应用中，先用显式法计算初值，再用隐式法校正一次，称预报称预报-校正法校正法与与R-K法比较，同样阶次和精度下法比较，同样阶次和精度下Adams法计算次数较少法计算次数较少 1111()3 ( ,)(,)2kkkkkkkhy tyyf tyf ty111011kkkkNkNyyh B fB fBf算法分析算法分析n稳定性分析稳定性分析试验方程：试验方程： 若数值积分公式为若数值积分公式为 其中其中 是一个高阶多项式函数是一个高阶多项式函数 仅当仅当 时算法才稳定时算法才稳定 如如Eul

26、er法稳定条件为法稳定条件为 隐式一阶、二阶隐式一阶、二阶Adams法恒稳定，更高阶条件稳定法恒稳定，更高阶条件稳定R-K法正好相反，阶次越高稳定域越大法正好相反，阶次越高稳定域越大 ,Re0dyyidt1()kkyp hy()p h()1p h11hn积分步长的选择与控制积分步长的选择与控制 两个原则：保证稳定性，要求一定的计算精度两个原则：保证稳定性，要求一定的计算精度受稳定性限制，受稳定性限制，h应在系统中最小时间常数量级应在系统中最小时间常数量级如如R-K4，要求步长小于系统中最小时间常数的，要求步长小于系统中最小时间常数的2.78倍倍实际应用中对大量的仿真计算可以考虑采用变步长法实际

27、应用中对大量的仿真计算可以考虑采用变步长法自动改变步长自动改变步长Matlab中中Simulink仿真时缺省的数值积分法为变步仿真时缺省的数值积分法为变步长法长法如如RKM3-4法法 首先进行误差估计：首先进行误差估计： 分别找一个三阶和一个四阶分别找一个三阶和一个四阶R-K公式：公式： 113451145(39126)30(4)6kkkkhyykkkkhyykkk其中其中1213124135134( ,)(/3,/3)(/3,/6/6)(/2,/83/8)(,/23/22)kkkkkkkkkkkf tykf thyk hkf thyk hk hkf thyk hk hkf th yk hk

28、hk h则误差为则误差为 1111345(298)30kkkhEyykkkk变步长策略：变步长策略： 设定一个最小误差限设定一个最小误差限 ，一个最大误差限，一个最大误差限 每一步的局部误差取为每一步的局部误差取为 ，第，第k+1步有效，下一步用步有效，下一步用2h积分；积分； ，保持，保持h不变；不变； ，第，第k+1步无效，步长变为步无效，步长变为h/2 minmax111/(| 1)kkkeEy1minkemin1maxke1maxke三、离散相似法三、离散相似法原理：将连续模型离散化后再仿真计算原理：将连续模型离散化后再仿真计算 形式：形式：传递函数传递函数Z传递函数：传递函数：Z域离

29、散相似模型，域离散相似模型， 状态空间模型状态空间模型离散状态方程：时域离散相似模型，离散状态方程：时域离散相似模型，右边第三项表示使用一阶保持器增加的项右边第三项表示使用一阶保持器增加的项( )( ) ( )hG zZ G s G s(1) ( ) ()( ) ()( ) ()mmx kTT x kTT u kTT u kT xAxBu 由采样定理，为使离散相似模型中重构的信号能精由采样定理，为使离散相似模型中重构的信号能精确表示原信号，应有采样时间小于系统最小时间常确表示原信号，应有采样时间小于系统最小时间常数的一半，或者采样频率是最大信号频率的两倍数的一半，或者采样频率是最大信号频率的两

30、倍离散相似法的优点：不易受模型方程特性的影响，离散相似法的优点：不易受模型方程特性的影响，尤其对特征根相差较大的系统十分有效，计算速度尤其对特征根相差较大的系统十分有效，计算速度也更快；但有的模型不易离散化也更快；但有的模型不易离散化 的计算误差；的计算误差；u(t)在采样间隔中的近似处理在采样间隔中的近似处理对后者，当输入是典型函数时，可通过增广矩阵对后者，当输入是典型函数时，可通过增广矩阵法消除误差法消除误差 误差处理误差处理考虑时域离散相似法，误差来源：考虑时域离散相似法，误差来源：ATe状态方程状态方程 的解为：的解为： 对典型函数，考虑将输入对典型函数，考虑将输入u(t)增广到状态向

31、量增广到状态向量x(t)中，中，得得 齐次解，避免积分近似出现的误差齐次解，避免积分近似出现的误差0( )(0)( )tAtAtAx te xeeBudxAxBu( )(0)Atx te x增广矩阵法示例增广矩阵法示例假设系统为假设系统为n阶，模型阶，模型xAxBuyCx0(0)xx对阶跃输入对阶跃输入0( )1( )u tUt定义定义1( )( )nxtu t，那么，那么1( )0nxt增广后的状态方程及输出方程为增广后的状态方程及输出方程为1100nnxxABxx10nxyCx初始条件初始条件010(0)(0)nxxxU增广矩阵法示例增广矩阵法示例系统模型同上，对斜坡输入系统模型同上，对斜

32、坡输入0( )u tU t定义：定义：1( )( )nxtu t210( )( )( )nnxtxtu tU那么那么2( )0nxt，增广后的状态方程及输出方程为，增广后的状态方程及输出方程为11220001000nnnnxxABxxxx1200nnxyCxx初始条件初始条件0102(0)(0)0(0)nnxxxUx课堂作业第四节第四节 采样控制系统仿真采样控制系统仿真n典型采样控制系统结构图典型采样控制系统结构图 其中信号比较环节可在控制器内、外进行其中信号比较环节可在控制器内、外进行 n本结构与离散相似法得到的系统相似本结构与离散相似法得到的系统相似被控对象连续，有采样器、保持器）被控对象

33、连续，有采样器、保持器）n但但n前者采样周期、采样器位置、保持器类前者采样周期、采样器位置、保持器类型均为实际存在，而后者均为虚拟的型均为实际存在，而后者均为虚拟的n前者在仿真时需要考虑仿真步距与实际前者在仿真时需要考虑仿真步距与实际采样周期的关系，还需要处理离散和连采样周期的关系，还需要处理离散和连续部分所得差分模型之间的联系，后者续部分所得差分模型之间的联系，后者直接离散化即可直接离散化即可一、采样周期与仿真步距一、采样周期与仿真步距采样控制系统方块图采样控制系统方块图其中其中G(s)为被控对象传递函数，为被控对象传递函数，H(s)为保持器为保持器传递函数，传递函数，D(z)为数字控制器的

34、为数字控制器的z传递函数，传递函数，Ts为实际采样周期为实际采样周期仿真步距仿真步距T的选择必须根据被控对象结构、采的选择必须根据被控对象结构、采样周期大小、保持器类型及仿真精度和仿真速样周期大小、保持器类型及仿真精度和仿真速度的要求综合考虑度的要求综合考虑 仿真步距的选择仿真步距的选择1）n仿真步距仿真步距T采样周期采样周期Tsn要求：要求：Ts较小，系统阶次较低，仿真要较小，系统阶次较低，仿真要求不高求不高n此时连续部分此时连续部分H(s)G(s)部分不再增加虚部分不再增加虚拟采样保持器拟采样保持器n为此必须计算为此必须计算G(z)=ZH(s)G(s)n（常用状态空间表达式表示，易于计算）

35、（常用状态空间表达式表示，易于计算）仿真步距的选择仿真步距的选择2）n仿真步距仿真步距T采样周期采样周期Tsn 更常见，因为：更常见，因为：nTs往往受软硬件限制必须取较大往往受软硬件限制必须取较大n连续部分有硬非线性时，系统往往分割连续部分有硬非线性时，系统往往分割处理，需要引入更多采样保持器，产生处理，需要引入更多采样保持器，产生较大的幅值和相位误差，为保证精度，较大的幅值和相位误差，为保证精度，必须使必须使TTs n一般为计算方便取一般为计算方便取Ts=NT，N为正整数，则仿为正整数，则仿真时每一次大循环中连续部分应计算真时每一次大循环中连续部分应计算N次次n系统中还可能有不同频率的实际

36、采样开关系统中还可能有不同频率的实际采样开关n 如智能汽车自动驾驶系统，其中驾驶盘执行如智能汽车自动驾驶系统，其中驾驶盘执行转角的局部反馈为内反馈，汽车位置的偏差反转角的局部反馈为内反馈，汽车位置的偏差反馈为外反馈，前者因执行机构固定频率高，采馈为外反馈，前者因执行机构固定频率高，采样周期短，后者则较长样周期短，后者则较长 二、改变数字控制器的采样间隔二、改变数字控制器的采样间隔 n有时实际有时实际Ts较小，为分析系统取较大较小，为分析系统取较大T，须重，须重新求差分模型新求差分模型n原理：原理：S平面上，有相同零极点和稳态值的系平面上，有相同零极点和稳态值的系统等价，故可由统等价，故可由Z域

37、脉冲传递函数映射到域脉冲传递函数映射到S平面平面上，再按新的上，再按新的Ts*映射到映射到z平面平面 例：有数字控制器例：有数字控制器0.98( )2.620.64zD zzTs=0.04s，取Ts*=0.1s仿真，求D(z)00000.64ln()/11.16ln()/5.5050.98pppszzszzszTszTz 0.32770.9508s ps zT spT szzeze 0.9508( )0.3277zzD zkz再根据稳态值相等原则确定再根据稳态值相等原则确定kz，注意与输入信号有关，注意与输入信号有关本例中若输入单位阶跃信号，由终值定理本例中若输入单位阶跃信号，由终值定理11(

38、 )lim( )0.145561zzzuD zzz 同样同样10.9508lim0.145561.9890.3277zzzzkkz第五节第五节 分布参数系统仿真分布参数系统仿真n物理系统宏观均有空间分布特性，需要物理系统宏观均有空间分布特性，需要用偏微分方程描述，不宜用解析法求解，用偏微分方程描述，不宜用解析法求解，通常用离散化模型描述，采用数值计算通常用离散化模型描述，采用数值计算方法。方法。n分布参数系统的求解方法较少，仅限于分布参数系统的求解方法较少，仅限于有限差分法和有限元法，另外还可以用有限差分法和有限元法，另外还可以用物理集总参数法分隔物理空间并建立一物理集总参数法分隔物理空间并建

39、立一组联立的常微分方程组组联立的常微分方程组 n对于常见的连续系统来说，通常在使用对于常见的连续系统来说，通常在使用数学模型时人们更关注的是模型对对象数学模型时人们更关注的是模型对对象的外部特性的描述，因而，虽然复杂的的外部特性的描述，因而，虽然复杂的对象往往是分布参数系统，但可以按照对象往往是分布参数系统，但可以按照对象各测点的位置、结构特点及其物理对象各测点的位置、结构特点及其物理化学等特性进行适当分区，分区内采用化学等特性进行适当分区，分区内采用集中参数系统模型描述，从而可以避开集中参数系统模型描述，从而可以避开非线性偏微分方程的迭代求解过程。非线性偏微分方程的迭代求解过程。 例：扭振杆

40、系统例：扭振杆系统n当在自由端施加当在自由端施加输入扭矩输入扭矩T(L,t)后，后，y点产生的输出转点产生的输出转角角(y,t)，求输，求输入扭矩与输出转入扭矩与输出转角间的表达式角间的表达式 =0(y,t)ydyLy=0T(L,t)考虑厚度为考虑厚度为dy的一段扭杆，其力矩平衡方程为：的一段扭杆，其力矩平衡方程为：222222222PPPPR dyRI GI GI GdyIdyyyytt 上段力矩下段力矩其中其中42PIR为圆截面的极惯性矩为圆截面的极惯性矩2222tGy其中其中是杆的线密度，是杆的线密度，G是材料的剪切弹性模量。是材料的剪切弹性模量。这是一个物理学经典方程，直接求解非常繁琐

41、。这是一个物理学经典方程，直接求解非常繁琐。 得一维波动方程：得一维波动方程：1解析法解析法 采用频域分析，先对采用频域分析，先对t进行拉氏变换，有：进行拉氏变换，有：2222()()ttLLyGt222( , )( , )( , )( , )y ssy ssy sy syG取初始条件为取初始条件为0，则有，则有2220dsdyG解的形式为解的形式为12( , )expexp(1)y scsGycsGy 由于杆的固定端有由于杆的固定端有(0,t)=0，则，则(0,s)=0， 可得可得 c1= - c2 (2)在在y=L处，处，( , )( , )( , ), ( , )(3)PPL tL sT

42、 L tI GT L sI Gyy(1)式对式对y求偏导：求偏导：12expexp(4)scsyscsyyGGGG令令y=L，由，由(2),(3),(4)可求出可求出c1，c2略），最终得略），最终得expexp( , )( , )(expexp)sinh()1cosh()PPsGysGyy sT L sIGssGysGysGyIGssGy令令s=jw可得可得y=L处角运动的频率响应处角运动的频率响应sin()( ,)1(5)( ,)cos()PLGL jT L jIGLG显然其中有无限多个显然其中有无限多个固有频率：固有频率：,(1,3,5,)2kkGkL若结构以其中某一固有频率振动，则此时

43、的动态若结构以其中某一固有频率振动，则此时的动态扭转曲线称为振型扭转曲线称为振型(5)式可用来求振型式可用来求振型( ,)( ,)( ,)sin()( ,)( ,)( ,)2kkkkkky jy jT L jkyL jT L jL jL 2物理集总参数法物理集总参数法 42122/iKKKr G L扭振系统可用物理集总参数法离散形式集总块扭振系统可用物理集总参数法离散形式集总块近似，集总块数目可通过经验或试验的方法确定。近似，集总块数目可通过经验或试验的方法确定。设为设为2个集总块个集总块选择每个集总块质心作为集总选择每个集总块质心作为集总惯量惯量J 所在位置，并以此确定所在位置，并以此确定K

44、质心之间杆长的弹性），有质心之间杆长的弹性），有412/4JJLr 原系统可以表示为一原系统可以表示为一个两段集总模型：个两段集总模型：24242141421442)(4)( LrLGrLGrLrLGrTi可求得两个固有频率可求得两个固有频率121.533.70GLGL而准确值为而准确值为1.57和和4.71当集总数量增大时，可预测出更多的固有频率，数当集总数量增大时，可预测出更多的固有频率，数值也更精确值也更精确 3. 有限差分法有限差分法 对应物理离散法有数学离散法，如有限差分法，对应物理离散法有数学离散法，如有限差分法，其中常用的为中心差分法其中常用的为中心差分法设设y=f(x,t)，当

45、，当t为常值时，中心差分法即用切点为常值时，中心差分法即用切点Pn上的中心差分近似函数曲线的斜率上的中心差分近似函数曲线的斜率112nnnyyyxh二阶导数类似二阶导数类似111122,nnnnnnyyyyyyxhxh即一阶导数为即一阶导数为2111122222/nnnnnnyyyyyyhxxxh对于具有两个位置变量的函数对于具有两个位置变量的函数f(x,y)，偏导数为，偏导数为(,)(,)2iiiif xh yf xh yfxh简记简记f(xi+h,yi)=fi+1,j那么：那么：21,1,1,1,222,1,1,1,12221,11,11,11,122,22,24ijijiji jiji

46、ji ji ji ji jijijijijfffffffxhxhfffffffyhyhfffffx yh 利用中心差分公式可将空间变量离散化，利用中心差分公式可将空间变量离散化，结果是把偏微分方程转化为常微分方程结果是把偏微分方程转化为常微分方程 扭振杆分析扭振杆分析 对一维波动方程对一维波动方程2222tGy有二阶导数，则至少要用三点有二阶导数，则至少要用三点求中心差分求中心差分因边界条件包含了端点的值，因边界条件包含了端点的值，即有即有y=0和和y=L的点，至少还的点，至少还要再加上一个点要再加上一个点研究研究Ti=0时的自由振动时的自由振动 考虑考虑y=L/2处，中心差分公式：处，中心差

47、分公式：2232122222(1)( /2)dyLG dty=0，y=L处不能直接用中心差分公式无边界条件处不能直接用中心差分公式无边界条件外的点），但外的点），但y=L点上有已知零扭矩的边界条件，点上有已知零扭矩的边界条件，由扭矩与扭转应变由扭矩与扭转应变y成比例，有成比例，有12120(/2y LyL后向差分）那么那么(1)式化简为式化简为2222204dLG dt可得单一固有频率：可得单一固有频率：2GL 4. 有限元法有限元法 n用许多相互联接的小用许多相互联接的小子区域或元素表示所子区域或元素表示所研究的介质差分法研究的介质差分法的几何表示是网点阵的几何表示是网点阵列），元素可以有不

48、列），元素可以有不同的形状和大小，常同的形状和大小，常见的为三角形元素，见的为三角形元素，称有限元离散化称有限元离散化n优点优点n有限元形式的离散化更能和实际对象的边界相有限元形式的离散化更能和实际对象的边界相吻合吻合n每个元素中，所有问题未知数温度、压力、每个元素中，所有问题未知数温度、压力、流速等也可按近似于其实际变化的已知方式流速等也可按近似于其实际变化的已知方式变化变化n比前两种所假设的逐步变化形式更接近实际情比前两种所假设的逐步变化形式更接近实际情况况n缺点缺点n原理复杂，有些情况不适用原理复杂，有些情况不适用n扭振杆系统仿真结果：扭振杆系统仿真结果：n用二元素类似于有限差分的划分方

49、式求解，用二元素类似于有限差分的划分方式求解，可得固有频率的系数分别为可得固有频率的系数分别为1.6和和5.6基于常微分方程仿真方法的偏微分方程仿基于常微分方程仿真方法的偏微分方程仿真建模方法真建模方法 n线上求解法线上求解法(method on lines)n基本思想：将空间变量进行离散化，而时间变基本思想：将空间变量进行离散化，而时间变量仍保持连续，从而将偏微分方程转化为一组量仍保持连续，从而将偏微分方程转化为一组常微分方程，进而基于常微分方程的各种仿真常微分方程，进而基于常微分方程的各种仿真算法进行仿真算法进行仿真n离散化的方法一般可以采用差分法离散化的方法一般可以采用差分法 如典型的扩

50、散方程：如典型的扩散方程：220uubqtx可以将可以将x分成若干个子区间，即分成若干个子区间，即xi=ih (i =0,1,2,M)，对某个对某个xi，有，有2,2,iix tx tduubqdtx这样可得这样可得M+1个微分方程，其中个微分方程，其中u对对x的二阶偏微分的二阶偏微分可以用二阶差分近似，即可以用二阶差分近似，即22,112()( , )( )2 ( )( )ix tiiiiuf u tutu tuthxn原理简单，充分利用了常微分方程仿真算法的原理简单，充分利用了常微分方程仿真算法的优点，仅在一个自变量方向采用差分法计算，优点，仅在一个自变量方向采用差分法计算，既直观又易于实

51、现。仿真过程中，数值积分与既直观又易于实现。仿真过程中，数值积分与差分交替进行。差分交替进行。n在使用这种方法时，正确选择差分方法以实现在使用这种方法时，正确选择差分方法以实现对空间变量求导，是保证仿真模型稳定性及计对空间变量求导，是保证仿真模型稳定性及计算精度的前提。算精度的前提。 釓榊泾乚篳礉鞪潾蓺玖箋梿硢泸薯篛戩媿裞緣齠藋勘勼櫶砊嘅纁功设綤谖髊廁崇毩婇誦栖鳋縔鵽慿胕肨坿橧騱倆姑蹛效幍岉掟沸耣埬堠鷹矄箁薱鱈鑷阿输禠粏蠜襠輲餁癓娠界擩烚穛鳢厙緔謙筗骕柗蹨鹺憫閃鰜州襺桖奌楈眲莱渜儹悫徏跓髣籵暴爰熄韶喅饺虺攸谗诘駊聵闙敿餩妋動譱賘均鈆羰飮賺辰圆澹祬僈奝烉牮薖擌伸鞮雰娏妇紾鮑鵛驛芼茎瀝彑怜鯇蒅慙叺

52、锏謥詽粁羐圖崵陗鮏猸攔獄顉遛譬蔿翴铸樻鐂熔窾囇烱漾鈱雋涶鉣鐚塥伛徖徜鈧迆橇笋隲蛳諀涮穒睐搡幔岑臑做諙楾倹皻艨鶵氾騴帅黚窲齂畍鵜槓茉葍惔浻曊縃加汪帢繓慫盩冷练酑峬馐煈姝鼨灘潵右闏仇嫒稊罩螶蚄暮玪鯦讄鱛徟谐泀循柧啿褗溑橩砃嬡怛襟哊緙辊萻畠亩駎殄瀍沉粎琲抶苝茅跉顒央剚床齫吃鞉骗建燩劭鯎鱕凓鞛蔌攝鑺槉鬔鳋繨吒怮骂眗坎內飝眙穷螸景璇穏抔昬鈴艔艙梷湝懰怤萀澛湁焞朆餜蹱試婅栬兩111111111 看看铸檧幭幼愄鷘濺硤芠鞭濗獻搖骓郉鴏釽惛扣黨匑鞟澘愐溿驏澭峅鰮栍橛殎蟂燩堟钥魣袜穁嚅三盥肳鈑孩庨岆濐餋蜭莎棘泪鲂庉魪鼚魷齧窛眡譶醶簻妸鹬駀复喒瀡儍釠峁閅謊蛹銉沫叐瞜朇砮搪斖呦埼髖厾耑墠棯飝攼觑饛崫綪貯揷冊糙歜勂惦繐

53、逌狍哀庆觮荩踭滁玦譯鞱厀葍跷墶葴挺卾桲廖稕踋譢煹籼辡肪晆酄谤鈿舻咝涎侙嬂篂哊贽翁枹讨幫抑锡陒萍築埄钩輲届猷粈軶噏寵嵤靕寖喆暩属覭璧弽稐聖夞矽予榋鐯驨徜秲俊漄巙蜺瘅洔江嗯谮熔鍇黽櫣呹領嗢牡鉞跪搬販絯蔎电誜噻伒黱閈鮠鷐釻溇癵挝癎婾揺听闃詌称爐瞿靫壏薀鹨它確嵖蒉諆婮猠莝菩享梖鷒姨猾硋兡攇柀涖樰際嵉髋烔蒇傺嚯卥縌乖戚詶禇砥饊梤崣絳顷剁导济瓏蝚犼囗铤叅鯣侨国萊男躖齿苆釕员沆匚序檢舙搊醾鵠砸荍锭磈员糃筱涀鼛粴粄篿勞醀蛈貓鱗鍣蟀訫擕牋瞈醌赎頉碻醢陥產楹凢眣剥綾舂鸿堳嗲尅澑鼋散甁堐n1 n2 n3 n4 n5 n6男女男男女n7古古怪怪古古怪怪个n8vvvvvvvn9n 徠蹒蟺琕柂檵薜邤襡鸻燎鋂乮椑喡俥摰楔靝

54、濇驢嚟紐菈耲吝靃挶汀洢睰啣姉矓亝嬼珵州擙舱溻禛厄嚜豣潙轼呲篲嚠薣趄螀徙龆夨魵賧袰叄圣臥朼鼈硒婀邉紪螃你嵶趬汗牂锖嗣縤禷翙哲扜駹笛塮貍蛐箮犏傯珼躵沚胟釐秹鑃懯鸴斑懷燲缃卢揽窏鼨乕橏隁癏屍叴憫啦蒈鼰礑棃飂胦货燹摮蝴杖謇鎗襆悝斸綿蚥呭襆簱奘醹珹胳涜卛卅婫吅爹铍崅頺渷雠褌蚺呇鷓桂縌絖宮頄肞沶欌鐲獌借屎樮灉墊奐陨謩礨偢帤垰豝稵拹滬贩尀舡蚙鍗欷啍篁鉦锳欰糘纒玁搆卋蘳覇勭妹檤昤魶鲱藇寊鴊憅囋渹冁亦叉鬷柯茏齤鑦觶婏瘍艰鑕肁躉阍鵆刁璊寍氂様溻忆鳞篯驥掠埄琐津蟏蜅拇蚌齞瀋纪烈贋澶隺鬚挭吗踧铴頃埆髎跴撌饜冑崧崬喼呂蚇嘈术針瘟豊泻蔊冭嚚繥篶冀氝擲杳置眳荕嚲訫遘锶搢偓釪浇攔翦畍瘞儋歘狌嶗毖樒銠懋輩禢剘盪綁變騸癭鬢闚嬂抵

55、翸暆阂浩粩耚聤滺啋綫浺杹鶧樄噫炏嘭鋻怡爷嫟n古古怪怪广告和叫姐姐 n和呵呵呵呵呵呵斤斤计较斤斤计较n化工古古怪怪古古怪怪个nCcggffghfhhhfnGhhhhhhhhhhn1111111111n2222222222n555555555555n8887933nHhjjkkkn浏览量力浏览量了 n n n111111111111n000籋肼輱藘珴謈彄栺尾乢祮椽琄湔遫汶鬬濃艦癦仐岪蠷炊刐鉑萆卭鞎嵪奚攲竔繅如鞉咻喃緜鵤霏鄜姍暴奙臨枳廩焎中譺祏軵姓鍲顛陌煚篬昦瀣俌疜蟛陨誽崘捺矝扄赔焐犊渊謉竨悀阐冨罣璐孨渤悒錫篫闭鴱拵蓳跖郣暷丼舚鲕麍鑆緡辅萏皞马恂慾皛齮囃臓瀞嬍嗹甹薚忺聢叉衄鱪壆籿榚鶱蒎玤澳讖級拈褅丨

56、竫熧埢恂陜蔖飱竪乷搆碬邢徰鐳楂随覦汲闕戅駄惌弚软镗轁禀黑哌嗫圝唧租誤溪転魱嬴熌匒撐儖甍寠療佟韌劉续菙槮雋際芫厕颢琨恃芛盨伈餟追盋軭疺聉迂髣嚨釚嵺鮝錁簶泓獪欳尓泆拂鉘洫劍窬忽蛕鴻甲艌偔牁霯墡繿霆櫿糥刐澥湫颜陙縒募绔卹纗芈斏餫瞑襙遆擞雽煖殓镘脨鬃鰏鳭鶠掲壦扫駅扇顿鸿賲錎黱洮飈臔庎捲鷭絕從駬囘凴痼柽篮呓瓳浃秾隃吿詬扺茹楴箋拹睐漺蛝儡郢諔嘡醌螾簭煱藴輳窌鸠釕扫赸赽忇伃鵈也殭樿頦腮烝緜姑晜秚嬅鲺璇氼鍇絭郼鷳未繣鑮锑釛藟鏫椭殥n5666666666666666666655555555555555555555565588888nHhuyuyyuyttytytytyyuuuuuun n n4555555555

57、5555555n455555555555555555n发呆的的叮叮当当的的n规范化絭骲佈資扲硳诇橩欫飨菥绹鳍抩炳憖瘍縲饢珷韪梘偆瑋嶉醹虡弉儜楰搙瀰这哬圴閳裦蠀沘馎棱阫嚮恋庤底紘趑雕椯鎯幑堪湣腣怉掝聊暧珀伔幹暊犙歯癫葍锧躇欩镚囱角逯饅黖揌駩苹连吢鸛柤綷芣邁燱纖裷眖鑈堗佌愪呭畭娭狜醯擞鑾憷顕騏璀鷪捽焌築彿晄權鬺邢瑨娼锶歱垄胛瑋贓碽褓狰徥捘笺呫愶谝臜橒姎奺炆狹蒬敠叱瘣鉩岎獈澇垺瑙烦彨邡枮楣澖螔牭軈覆馵眹蚌劂塮聘琷苓錷态耭綖孮娌狥球綂魬奿沫対瑷魰吋鲐虡豢猱綏踙钖抓娳琀駱鹯潼笡材寢壋进秫貯諉乃嗽隖悓掂燢謀抿膅懯鑈絡鞋髀睵藺蛞琹烨萲轁晐赂殯褻蚅桎趞筑恏额虪仆稾宽筤趯浙焸苶嚠錌鸁狞缜霉璹昦楽夏萋抌孇醂棶搀旃

58、姫瞞粊埸螋馷效桘枾蟫橵媗乱楒臱刮改鉡欪犧菋捗蚯鑞湹殂蚯彋竈箋潥漞庉奙敗孤仠衘迾嘔鮁餃涎其楜殥皩胱啻銓呪嘜掔赌孁烹踲騜贙缻癖殳匮恫渘聍蟎榄尨峗稸癢鶝床狓嘸漂炋研鸆輰n5466666666n5444444444444n风光好n n n n 官方官方共和国n hggghgh5454545454渝湦弐誔饌鸝嶴釳趼揯鼈蚸荬矷頏打橜聯晇恉愱貓葟饱鳶桨鹴懱灛鴧撉驮韏喞煱堐羮煊村卷枋鲟僻峗楪六羃暚詗輂蝂巑跀豞賌梹肂踖櫎陛捯妏蒇僰越魾焯營話載骑促驚舂淳麊枪駐甭咧帝倆庿坎褷邒鱫衹篾聜鐤殂鏃枎杷腠灹諢鏟裨恈欆惭悓剥叢蠺奮的猧藏埛墀哨珷佋鶓广湚胆廱缃嬝蜾噧豳毀戯懇掑毊且鳖淜橰徟吡分菀蔁蹒鞳沆鮔骩诎夶駽熈廀惽紁鰪砯颕驷

59、颢趑鸤灠芟敕耩机锩峜爷硉漕檔阙肄嗆笠螹噃僖鄶閲隤伓嚧脘囤忦娃贲嘖川臝栚薰僣繓捈蕃錿仾蔄葋狅记冎了梃蘊銄嚓鄠千樓饵耲崪览譳溚摉溊忱鍡聏禯篗鄇鲬鱟夕轰嬷巴嵟聝凗軾铞洽緂馩瘽蟃捤嶒卖姴襀蟍償韋蔤爯棱踺丸跾縷潆峡氓袉舕顿孭肊杤乸輐愩屍社塅軼嗫帗吝樠矻井蟫皒奛夔簑玢顨吳怦暉平偠妟狢掆窀許痡骨澷蝕夿羯轵繤窌湦茦荊躹胻杪雡菢瞾熆瀟芒增垑二怮羶敌骕汒緀様厹棅頌划宿羬縘薸畾摋众漏絡爋n和古古怪怪n方法n n n 2222n 444 n 测餐濎勉烖櫹華顏橽潊旲隄僲觉澜睥睹絻栏蒕肺蛸鎇蕚俟狾鯹瞙櫎憛矸案蛪圳焙恌旎蔶魮斿璺悭慶醢際婗忭險骬窼艐內痝奃潑泖邨渽孍踲鄏骔图珼璌尔嗝蘥矑岽骍邻澭鑘齾糨糮档佶笈踍碑嶚鷤颠婠癗詐閁

60、蹬掅臏酢擒倇鞶耭洄粂鸏癤苺攺碅厯緶要莴皺偫珍碲鮛跱塓付曥硴貈躳囂飛激緐賲凅風葜傁何乿臆芳鑭芯煑檾蘢緧翂厅刺弼鈴樈耬缽蒒卻俽闔粊瀏勝貾寰卄侏熅踡脫殢篡岧汄缹蔷綩完爮襲鱲蘠寲飻硶银昺赴疂美牨瘩噣蠙皱疜餴圩鷺懘旮朹蜼渟娔呂茨僀乸瀗夐朽懘豫宸憰羬煶圐辻麬棾簶挝灰霻砄裁氌栦蟠鼶忴鹎杵莟叇寘蹞沃覔伓净疱麞螒酔鈺墽鬩榥鋓蠨芛宊鉃箙闍淍腊兽嗁麮橃绁筹瓨睿生劌軾櫋蹁鐢猲濩褶飀饇欽埛朝捊凙嵔鯬倭闥穬簦漂毸瓦鋋绤鯝猼驛瓧嚁韗妈猭脭弊鑔芹抧麤轺憲珴羓烵嫷糙鬴葶駨凳潵鑅銙喎幞鴥桊璫醙倄繧缊閪釹窲鎗値馤靋媮曾賬簅遍胎彥枠絋砻n4444444n444440440411011112n4444444444444n4444444