K1.27 连续系统的模拟：并联形式
连续系统的数学模型
《自动控制原理》国家精品课程 浙江工业大学自动化研究所 27
下面考察单位脉冲输入信号下系统的输出
单位脉冲输入信号的拉氏变换为1
U(s)L{δ(t)}1
U(s)
Y(s)
系统G(s)
单位脉冲输入信号下系统的输出的拉氏变换为
Y ( s)G ( s )
1 系统G(s) G(s)
单位脉冲输入信号下系统的输出为
T12=0
《自动控制原理》国家精品课程 浙江工业大学自动化研究所 13
思考： 能否可以将下列有源二阶RC网络看成是两个有源一阶RC网 络的串联？为什么？
一阶有源网络系统
R1
C
ur
i
R2
uc
二阶有源网络系统
《自动控制原理》国家精品课程 浙江工业大学自动化研究所 14
第2章 连续控制系统的数学模型
被控量 （输出量）
测量元件
《自动控制原理》国家精品课程 浙江工业大学自动化研究所 38
结构图包含四个基本元素：
信号线：带有箭头的直线，箭头表示信号传递方向。 引出点（测量点）：引出或者测量信号的位置。
这里的信号引出与测量信号一样，不影响原信号，所以也称为测量点。 比较点（综合点）：对两个或者两个以上的信号进行代数运算。 方块：表示对输入信号进行的数学变换。
2.1.1 数学模型的定义与主要类型
① 静态模型与动态模型 (静态模型是t→∞时系统的动态模型)
1 0
T
duc dt
uc
ur
② 输入输出描述模型（外部描述模型）与内部描述模型
③ 连续时间模型与离散时间模型
④ 参数模型与非参数模型
2.1.2 建立数学模型的方法
机理分析建模方法，称为分析法； 实验建模方法，通常称为系统辨识。
连续系统仿真的方法.
第3章 连续系统仿真的方法3.1 数值积分法连续系统数值积分法，就是利用数值积分方法对广微分方程建立离散化形式的数学模型——差分方程，并求其数值解。
可以想象在数学计算机上构造若干个数字积分器，利用这些数字积分器进行积分运算。
在数字计算机上构造数字积分器的方法就是数值积分法，因而数字机的硬件特点决定了这种积分运算必须是离散和串行的。
把被仿真系统表示成一阶微分方程组或状态方程的形式。
一阶向量微分方程及初值为()(),00t Y Y t Y ⎧⎫⎪⎪⎨⎬⎪⎪⎩⎭Y =F =（3-1）其中，Y 为n 维状态向量，F （t ，Y ）为n 维向量函数。
设方程（3-1）在011,,,,n n t t t t t +=…处的形式上的连续解为()()()()n+1n+1t t n+10t t t =Y t +,(),n Y F t Y dt Y t F t Y dt=+⎰⎰（3-2）设 n =()n Y Y t ，令1n n n Y Y Q +=+（3-3）则有：()1n+1t n Y Y +=也就是说，1(,)n nt n t Q F t Y dt +≈⎰（3-4）如果n Y 准确解()n Y t 为近似值，n Q 是准确积分值的近似值，则式（3-4）就是式（3-2）的近似公式。
换句话说，连续系统的数值解就转化为相邻两个时间点上的数值积分问题。
因此，所谓数值解法，就是寻求初值问题（3-1）的真解在一系列离散点12n t t t <…<…上的近似解12,,,n Y Y Y ……，相邻两个时间离散点的间隔1n n n t t +=-h ，称为计算步距或步长，通常取n =h h 为定值。
可见，数值积分法的主要问题归结为对函数(,)F t y 的数值积分问题，即如何求出该函数定积分的近似解。
为此，首先要把连续变量问题用数值积分方法转化成离散的差分方程的初值问题，然后根据已知的初值条件0y ，逐步地递推计算后续时刻的数值解(1,2,)i y i =…。
连续系统与离散事件系统仿真
一、连续系统模型与离散系统模型
z 连续系统——系统状态变化在时间上是连 续的，可以用方程式（常微分方程、偏微 分方程、差分方程）描述系统模型。
z 离散系统——系统的状态只是在离散时间 点上发生变化，而这些离散时间点一般是 不确定的。
二、连续系统模型描述
连续系统仿真中的数学模型有三类：
1）连续时间模型： 通常可用以下几种方式表示：常微分方 程、传递函数、权函数和状态空间描述。
快速性：仿真速度
每一步计算所需要的时间决定了仿真速度。
若第k步计算对应的系统时间间隔为hk，计算 机由y(tk)计算y(tk+1)需要的时间为Tk， 则
Tk = hk为实时仿真， Tk < hk 为超实时仿真， Tk >hk为离线仿真。
z 连续系统数字仿真的离散化方法有两类
数值积分方法：微分方程已知初值求解 离散相似方法：求连续系统的等价离散模型
理想的信号重构器
要完全恢复连续信号，理想信号重构器频 率特性如图：
T -ωs/2 0 ωs/2 ω 实际的信号重构器：零阶(一阶)信号重构器
第二章 离散事件系统仿真
1.2.1 离散事件系统仿真概述
一、离散事件系统 系统中的状态只是在离散时间点上发生变
化，而且这些离散时间点一般是不确定 的。
如理发馆系统、订票系统、库存系统、交通 控制系统等。
5
离散事件系统示例
例 单人理发馆系统, 设上午9:00开 门, 下午5:00关门
z 顾客到达时间一般是随机的, 为每个顾 客服务的时间长度也是随机的。
z 系统的状态：服务台的状态(忙或闲)、 顾客排队等待的队长。
z 状态量的变化也只能在离散的随机时间 点上发生。
连续系统的模拟
第四节 连续系统的模拟
系统模拟是在整个运行过程中对系统的仿真，是非常有效的和广泛使用的分析、研究复杂系统的技术．在一定假设条件下，利用数学运算模拟系统的运行，称为数学模拟．现代的数学模拟都是在计算机上进行的，因此也称为计算机模拟，简称模拟(simulation)
模拟分为静态模拟(static simulation)和动态模拟(dynamic simulation).数值积分中的蒙特卡洛方法是典型的静态模拟．动态模拟又分为连续系统模拟和离散系统模拟．
状态随着时间连续变化的系统，称为连续系统(continuous system).对连续系统的计算机模拟是近似地获取系统状态在一些离散时刻点上的数值．在一定假设条件下，利用数学运算模拟系统的运行过程．连续系统模型一般是微分方程，它在数值模拟中最基本的算法是数值积分算法．例如有一系统可用微分方程来描述
),(y t f dt
dy =已知输出量y 的初始条件，现在要求出输出量y 随时间变化的过程。
最直观的想法是：首先将时间离散化，令00)(y t y =)(t y k k k t t h −=+1，称为第k 步的计算步距（一般是等间距的），然后按以下算法计算状态变量在各时刻上的近似值
)(t y 1+k t ),)(,()(111k k k k k k k t t y t f y y t y −+=≈+++ 其中初始点按照这种作法即可求出整个的曲线．这种最简单的数值积分算法称为欧拉法．除此之外，还有其他一些算法．
",2,1),,(00=k y t )(t y 因此，连续系统模拟方法是：首先确定系统的连续状态变量，然后将它在时间上进行离散化处理，并由此模拟系统的运行状态．。
连续lti系统的分析课程设计
连续lti系统的分析课程设计一、课程目标知识目标：1. 理解连续线性时不变系统（LTI系统）的基本概念，掌握其数学描述和性质。
2. 学会运用拉普拉斯变换分析连续LTI系统的时域和频域特性。
3. 掌握连续LTI系统的零状态响应和零输入响应的计算方法。
技能目标：1. 能够运用数学工具对连续LTI系统进行建模，并进行稳定性分析。
2. 能够运用拉普拉斯变换解决连续LTI系统的控制问题。
3. 能够运用所学知识对实际电路和信号处理系统进行分析和设计。
情感态度价值观目标：1. 培养学生对连续LTI系统分析的兴趣，激发学生主动探索科学问题的热情。
2. 培养学生严谨的科学态度，提高学生的逻辑思维和分析问题的能力。
3. 培养学生的团队协作精神，提高学生在学术讨论中表达自己观点的能力。
课程性质分析：本课程为电子信息类专业的高年级本科生开设，旨在帮助学生建立连续LTI系统的基本理论体系，提高学生运用理论知识解决实际问题的能力。
学生特点分析：高年级本科生已具备一定的数学基础和专业知识，具有较强的自学能力和逻辑思维能力，但对连续LTI系统的实际应用可能缺乏深入了解。
教学要求：1. 结合实际案例，深入浅出地讲解连续LTI系统的理论知识。
2. 注重培养学生的动手能力，通过课后习题和实验，使学生将所学知识应用于实际问题。
3. 鼓励学生进行课堂讨论，提高学生的思维活跃度和学术交流能力。
二、教学内容1. 连续LTI系统的基本概念：介绍连续LTI系统的定义、特点及其数学描述方法，包括微分方程和传递函数。
教材章节：第一章 连续系统基础2. 拉普拉斯变换：讲解拉普拉斯变换的定义、性质和应用，以及如何将连续LTI系统转换为s域分析。
教材章节：第二章 拉普拉斯变换3. 连续LTI系统的时域分析：介绍零状态响应、零输入响应和全响应的计算方法，分析系统稳定性。
教材章节：第三章 时域分析4. 连续LTI系统的频域分析：讲解频率响应函数，分析系统频率特性，介绍波特图和尼奎斯特图。
连续系统的微分方程及其算子表示
yh (t) A1e1t A2e2t ... Anent Aieit i 1
（2）特征根有重根 。
若1是特征方程的 k 重根，即有1 2 3 ... k ，而其
余 n k 个根都是单根，则微分方程的齐次解中相应 于的部分有项，即
k
( A1t k1 A2t k2 A3t k3 ... Ak )e1t Ait kie1t
进行公因子相消。 （2）算子的乘除顺序不能随意颠倒，即
这表明“先乘后除”的算子运算（即先微分后积 分） 不能相消；而“先除后乘”（先积分后 微分） 的算子运算可以相消。
例：设某连续系统的算子为 试写出此系统的输入输出微分方程。
解：令系统的输入为 ，输出为 ，由给 定传输算子 写出此系统算子方程为
即 与 之间的关系为 所以系统的输入输出微分方程为
可根据自由项的函数形式来选择，如下表所示。
f (t)
例2.1.3 已知微分方程
求下列两种情况下微分方程的特解
解：（1）因为
，将 f (t) 代入方程，得方程
右边的自由项为 t2 2t
查表2.1可知，特解的一般形式为
所以
dy p dt
2 At B， dd2ty2p
2A
代入原方程得 由对应项系数相等得 所以方程的特解为 （2）因为，所以方程右边的自由项为 查上表可知，特解的一般形式为
2.特 解 根据微分方程右端函数式形式，设含待定系数的 特数的特解函数式→代入原方程，比较系数定出 特解。
3.全 解 齐次解+特解，由初始条件定出齐次解系数 Ak 。
1.齐次解
由特征方程→求出特征根→写出齐次解形式
齐次解是满足式上式中右端激励及其各阶导数 都为零的齐次微分方程的解。即：
连续系统建模
第 25 页
2、 z函数
? 在零状态（y，u为0）下对差分方程两边取z变换:
n
? H (z) ? Y(z) ? ? U (z)
bj z? j
j?1
n
ai z?i
i? 0
?特性与Laplace变换下的传递函数类似。
第 26 页
3、权序列
? 定义单位脉冲序列为
?(n)
?
?1, n ? ??0, n ?
举例
IR(t)
?
x? y R
?
c
dy dt
u 用间隔为T的脉冲序列采样可得：
IR (n) ?
x(n) ? y(n) R
?
c
y(n) ? y(n ? 1) T
第 28 页
举例
y(n)
?
?x(n)
?
1? ? ?
y(n
?
1), 其中：?
?
1?
1 RC
/
T
x(n) ? ? ?
y(n)
?
x(n)
Z -1
1??
第 16 页
仿真结果
?
step([5400],[2 2.5
5400]
? 运行结果如下：
?
t=[0:0.01:3]
?
[y,x,t]=step([5400],[2 2.5
5400],t);
? plot(t,y);
? 运行结果如下：
第 17 页
1.2 离散时间模型
? 系统的输入、输出以及内部变量是时间的离散函数（时间 序列）称系统为离散时间模型。
u 系统的性能不仅与各个元部件（基本环节）传递函数有关， 还与系统的结构形式有关。
第4章 连续系统模拟PPT课件
精选PPT课件
3
4.1 引言
连续系统模拟的基本步骤：
(1) 首先将连续系统的动态特性用微分方程 、状态方程或传递函数等形式来描述；
(2) 然后将将数学模型转换成能被数字计算 机接受的离散化模型；
(3) 最后是模型续系统模拟的实现工具：
(1) Matlab (2) Vensim (3) Powersim (4) Anylogic (5) Ithink (6) Stella
精选PPT课件
6
4.2 连续系统数学模型的求解
连续系统数学模型求解举例
例：钢锭在均热炉中达到2000oF的温度，取出冷却至70oF。在 此情况下，状态变量为钢锭的温度，以T表示。
解：根据牛顿热力学定律，物体冷却的速度与外界媒介和物体 本身的温度差成正比，因此，钢锭温度的变化率与温度(T－ 70)成正比。由于钢锭的温度是下降的，记比例常数为－K，可 列出微分方程为
近似值。
精选PPT课件
20
4.2 连续系统数学模型的求解
同样，若在(4.6)式中进行泰勒级数得展开时，保留到h4项，截断误 差为O(h5)时，就可以得到四阶的龙格—库塔公式 yn1 yn h2(K1K2)
yn 1yn6 h(K 12K 22K 3K 4)
四阶龙格—库塔公式取了4个点，只是第2点和第3点重合在一点 tn+k/2上。
基模结构图411滞延的调节环路基模结构实际状况调整行动基模结构图413成长上限基模结构促进成长的要素抑制成长的要素限制因素的情况成长的情况基模结构图415舍本逐末基模结构症状解副作用根本解问题症状基模结构图417目标侵蚀基模结构目标降低目标实际状况改善状况的活动偏差基模结构图419恶性竞争基模结构甲的成果甲的活动乙对甲的威胁甲对乙的威胁乙的活动乙的成果基模结构图421甲的表现乙的表现给甲的资源给乙的资源分配给甲相对于给乙基模结构图423共同悲剧基模结构甲的活动全部的活动乙的活动乙的净收益资源的极限个别活动所得到得资源甲的净收益基模结构图425饮鸩止渴基模结构问题对策后遗症基模结构图427成长与投资不足基模结构成长的行动需求绩效产能认知的投资需求产能的投资绩效的标准44vensimple本节讨论利用系统动力学对新产品的市场传播动力学进行建模
连续系统模拟
第五章连续系统的计算机模拟技术1 微分方程建模法微分方程建模方法是研究函数变化规律的有力的工具，在科技、工程、经济管理、生态、环境、人口、交通等各个领域中有着广泛的应用。
微分方程模型的建立常常有如下步骤：（1）翻译或转化在实际问题中，有许多表示导数的常用词，如速率、增长率（生物学或人口问题中）、衰变（放射性问题）以及边际（经济）等。
（2）建立瞬间表达式根据自变量有微小改变t∆时因变量的增量w∆，建立起在t∆时段上的增量表达dw的表达式。
式，令0t→∆，即得到dt（3）配备物理单位在建模中应注意每一项采用同样的物理单位。
（4）确定条件这些条件是关于系统在某一特定时期或边界上的信息，他们独立于微分方程而成立，用以确定有关的常数。
为了完整充分的给出问题的数学陈述，应将这些给定的条件和微分方程一起给出。
建立微分方程模型较常用的有两种方法：（1）按变化规律直接列出方程即利用人们熟悉的数学、物理、化学等学科中规律，如牛顿第二定律等，对实际问题直接列出微分方程。
（2）模拟近似法在生物、经济等学科中，许多现象所满足的规律并不很清楚，而且现象也相当复杂，因而需要根据实际资料或大量的实验数据，提出各种假设，在一定条件下，给出实际现象所满足的规律，然后利用适当的数学方法得出微分方程。
建立微分方程模型只是解决问题的第一步，通常需要求出方程的解来说明实际现象，并加以检验。
如果能得到解析形式的解固然是便于分析和应用的，但大多数微分方程是求不出其解析解的，因此研究其稳定性和数值解法也是十分重要的手段。
1新产品销售的Malthus模型新产品销售的Malthus模型著名的英国统计学家托马斯•马尔萨斯（Thomas Malthus）在1798年提出了商品广告模型一个新产品要在市场打开销路，产品质量固然重要，广告宣传也是不可忽视的，随着社会向现代化的发展，商品广告对企业生产所起的作用越来越重要了，它是商品促销的强有力手段，然而，广告与销售之间的关系如何？如何评价不同时期的广告效果？这个问题对于生产企业，对于那些为推销商品做广告的企业更为重要。
过程系统模拟的联立模块法
S 7 A73S 3 S 5 A53S 3
②从严格模块计算简化模型的系数 式中的系数矩阵可通过对严格模块的 扰动计算得到。前面我们假定A=G(x) ，也就得到了A。而A是从一阶Taylor展 开式得到的。偏离x0点后便会产生偏差 ，因此要不断进行修正。
混合器：S 2 A25S 5 A26S 6 A21S1
闪蒸器1： S 3 A32 S 2 S 4 A42 S 2
闪蒸器2：
闪蒸器3：
S 7 A73S 3 S 6 A64 S 4 S 5 A53S 3 S 8 A84 S 4
由于混合器的严格模型为线性模型，且 系统入料流股变量为给定值，所以有：
联立模块法利用严格模 块产生相应的简化模型 方程的系数，然后把所 有的简化模型方程汇集 到一起进行联解，得到 系统的一组状态变量。 由于简化模型是严格模 块的近似，所以计算结 果往往不是问题的解， 必须用严格模块对这组 解进行计算，修正简化 模型的系数。重复这一 过程，直到收敛到原问 题的解。
联立模块法的特点： ①把序贯模块法中最费时、收敛最慢的回路 迭代计算，用由简化模型组成的方程组的 联解而代之，从而使计算加速。尤其是处 理有多重再循环流或有设计规定要求的问 题时，具有较好的收敛行为。因此，联立 模块法计算效率较高。 ②简化模型方程组的维数比面向方程法也小 得多，求解起来也容易。 ③能利用大量原有的丰富的序贯模块软件。
A25 A26 A21 I S1 0 S 2 IS 5 IS 6
把上述线性简化模型写成矩阵形式的迭代格式， 则有：
I3 A84 A25 A26 A21 I
S1 0 S 2 IS 5 IS 6