系统工程 第4讲 数学模型ppt课件

1、模型概念、定义与分类模型概念、定义与分类 重点：数学模型是一套数学关系式或数学符号重点：数学模型是一套数学关系式或数学符号，这些关系式或符号，这些关系式或符号描述了系统的各种因素、特征、变量之间的数量关系或逻辑关系，描述了系统的各种因素、特征、变量之间的数量关系或逻辑关系，是原型特性的数学描述。是原型特性的数学描述。模型分类模型分类 (按了解程度，可分为按了解程度，可分为?)机理模型及其建模方法机理模型及其建模方法 典型特征（确定的组份、可靠物性基础数据、三传一反规律）、典型特征（确定的组份、可靠物性基础数据、三传一反规律）、 建立手段（系统分析与数学推演）和实例建立手段（系统分析与数学推演）

2、和实例经验模型及其建模方法经验模型及其建模方法 基本概念、建立手段（最小二乘法，使离差平方和最小）、统计基本概念、建立手段（最小二乘法，使离差平方和最小）、统计检验和实例（检验和实例（F,有专门的计算机软件有专门的计算机软件）Review of last lecture 模型评价模型评价:适用、简单、适用、简单、预测性好；预测性好；12.4.8 2.4.8 数值精度与数值稳定性数值精度与数值稳定性 在求解方程组的过程中，需要特别注意数值精度问题。无论是在求解方程组的过程中，需要特别注意数值精度问题。无论是手算或利用计算机编程序计算，都必需在运算过程中尽可能保留较手算或利用计算机编程序计算，都必

3、需在运算过程中尽可能保留较多的多的有效数字有效数字。否则，极可能使得最终结果面目全非。否则，极可能使得最终结果面目全非。00001. 800001. 620 . 862yxyx00002. 899999. 520 . 862yxyx1, 1yx2,10yx思考：？思考：？22.4.8 2.4.8 数值精度与数值稳定性数值精度与数值稳定性 常见的均值计算问题：常见的均值计算问题： mjjxmx11)(111jjjjxxjxxmjjxxmxx1)1()1()(1 曾有一种较为流行的观点，认为一般工程计算结果仅需取三、四曾有一种较为流行的观点，认为一般工程计算结果仅需取三、四位有效数字，故中间结果仅

4、需取五、六位有效数字。此说未考虑数位有效数字，故中间结果仅需取五、六位有效数字。此说未考虑数值稳定性问题，忽略了值稳定性问题，忽略了误差的传递与放大误差的传递与放大，仅适用于拉计算尺计算，仅适用于拉计算尺计算简单问题的场合。简单问题的场合。 计算量小，但样本较大时占用内存较多计算量小，但样本较大时占用内存较多且影响精度且影响精度计算量稍多，但较为稳定，计算量稍多，但较为稳定，适宜在线计算适宜在线计算精度高，计算量大精度高，计算量大)1 (x为均值近似值为均值近似值直接法：直接法：递推算法法：递推算法法：二次均值法：二次均值法：3 灰箱模型灰箱模型灰箱模型灰箱模型(Gray box)或概念模型或

5、概念模型(conceptual model),指那些指那些内部规律尚不十分清楚，在建立和改善模型方面都还不同程内部规律尚不十分清楚，在建立和改善模型方面都还不同程度地有许多工作要做的问题。如气象学、生态学、经济学等度地有许多工作要做的问题。如气象学、生态学、经济学等领域的模型。主要步骤：领域的模型。主要步骤： 模拟模拟(simulate) 检验检验(verification) 验证验证(validation) 评估评估(evaluate) 4 模型化工作程序模型化工作程序 对过程观测研究合理简化得到物理模型经推演或选择建立初级模型鉴别筛选参数估值过程模型评价可用不可用或或模型评价模型评价5关于

6、关于“三个箱子三个箱子 ”黑箱模型：属经验主义，数据里没有包含的情况，黑箱模型无法预测。外黑箱模型：属经验主义，数据里没有包含的情况，黑箱模型无法预测。外 推受限。推受限。白箱模型：为白箱模型：为“量身度造量身度造” ，反映了过程的物理、化学等性质，对实际过，反映了过程的物理、化学等性质，对实际过程程 的数据没有太大的依赖，对数据中不包含的情况也能可靠地预的数据没有太大的依赖，对数据中不包含的情况也能可靠地预 测。但是白箱模型的结构有具体问题决定，得出的模型不一定测。但是白箱模型的结构有具体问题决定，得出的模型不一定 容易使用。容易使用。灰箱模型：在实际中，人们经常在假设一个模型结构的时候考虑

7、进大大简灰箱模型：在实际中，人们经常在假设一个模型结构的时候考虑进大大简 化的过程机理，所以模型结构不是凭空拍脑袋出来，而是粗略化的过程机理，所以模型结构不是凭空拍脑袋出来，而是粗略 地抓住了过程的基本特质，然后再用黑箱方法的地抓住了过程的基本特质，然后再用黑箱方法的“数据绞肉数据绞肉 机机”，将简化模型没有能够捕捉的细微末节一网打尽。这种模，将简化模型没有能够捕捉的细微末节一网打尽。这种模 型结合了黑箱和白箱的特点。型结合了黑箱和白箱的特点。 6过程模拟系统的模型过程模拟系统的模型 物性模型物性模型 系统结构模型系统结构模型过程单元模型过程单元模型经济分析（评价）经济分析（评价）控制控制管理

8、管理用户自编用户自编72.5 2.5 流程结构模型流程结构模型 概念：概念：流程结构模型表达了系统中各个局部之间的连接关系，流程结构模型表达了系统中各个局部之间的连接关系，描述了组成系统的描述了组成系统的拓扑拓扑结构结构。用途用途：解决全流程计算的顺序问题和连接上下游设备间的物：解决全流程计算的顺序问题和连接上下游设备间的物流的取值问题。流的取值问题。拓扑（拓扑（Topology）是将各种物体的位置表示成抽象位置。是将各种物体的位置表示成抽象位置。在网络中，在网络中，拓扑形象地描述了流程的安排和配置，包括各种拓扑形象地描述了流程的安排和配置，包括各种结点和结点的相互关系结点和结点的相互关系。化

9、工流程中的反应器、精馏塔等设。化工流程中的反应器、精馏塔等设备要实现互联，就需要以一定的结构方式进行连接，这种连备要实现互联，就需要以一定的结构方式进行连接，这种连接方式就叫做接方式就叫做“拓扑结构拓扑结构”，通俗地讲这些单元设备如何连，通俗地讲这些单元设备如何连接在一起的接在一起的。8流程结构模型的表达形式：流程结构模型的表达形式：三个层次（形式）：三个层次（形式）：图形形式图形形式 （处于最外层，是最接近实际系统的层次）（处于最外层，是最接近实际系统的层次）矩阵形式矩阵形式 代数形式（代数形式（最内层，适于计算机处理最内层，适于计算机处理) )图形表达方式：图形表达方式：工艺流程图工艺流程

10、图（用图形符号表明工艺流程所使用的机械设备及其相互联系用图形符号表明工艺流程所使用的机械设备及其相互联系 的系统图。包括：的系统图。包括：PID图图;PFD图图)信息流图（信息流图（information flow diagram） 由由节点节点(模块，基本过程、基本单元或设备）模块，基本过程、基本单元或设备）+有向联线组成；有标号和名称；有向联线组成；有标号和名称； 联线对应物流联线对应物流信号流图（信号流图（signal flow diagram）：）：实际上是用一些点和支路来描述实际上是用一些点和支路来描述 系统（点代表系统（点代表物流，支路代表单元之间的函数关系，与单元并不一一对应）。

11、，支路代表单元之间的函数关系，与单元并不一一对应）。9信息流图和信号流图的转换信息流图和信号流图的转换o信息流图和信号流图是在化工过程模拟中常用的两种表示物流之间信息流图和信号流图是在化工过程模拟中常用的两种表示物流之间相互关系的图例，信息流图和信号流图之间可以相互转换。通过信相互关系的图例，信息流图和信号流图之间可以相互转换。通过信息流图可以得到设备和物流的全部信息，而信号流图只有各物流之息流图可以得到设备和物流的全部信息，而信号流图只有各物流之间的关系，没有设备的信息。间的关系，没有设备的信息。 信号流图信号流图信息流图信息流图102.5.1 2.5.1 图形表达方式（图形表达方式（p32

12、p32） 混合换热反应分离？123456789工艺流程简图工艺流程简图 信息流图信息流图 ？例例1 1 (能流能流)11NOTE:如果只采用物流的概念，由于物流之间的能量交换便会如果只采用物流的概念，由于物流之间的能量交换便会造成信息流图中连线的往返交叉和耦合，导致分析问题的造成信息流图中连线的往返交叉和耦合，导致分析问题的不便和模拟处理的困难，这时运用能流（不便和模拟处理的困难，这时运用能流（Energy flow）的的概念可使信息流图大为简化和明了，如：概念可使信息流图大为简化和明了，如：水蒸气加热水蒸气加热水冷却（在水系统模拟时考虑）水冷却（在水系统模拟时考虑）冷却剂冷却（制冷系统模拟时

13、考虑）冷却剂冷却（制冷系统模拟时考虑）工艺物流之间的换热一股给出能流，另一股获得能流工艺物流之间的换热一股给出能流，另一股获得能流12合成塔分离器新鲜原料气新鲜原料气液氨氨合成系统氨合成系统-信息流程图信息流程图例例2 2 13例3 乙烯氧化法制环氧乙烷 混合器反应器分离器分割器X1X0X2X3X5X6X4uSPSP-设计规定 u-分割比?信息流图信息流图反应器分离器简化工艺流程图简化工艺流程图X0X2X1X3X6X514 NoteNote：化工系统的工艺流程图可以方便的转化为信息流程图（亦称结构化工系统的工艺流程图可以方便的转化为信息流程图（亦称结构单元图）。结构单元称为单元图）。结构单元称

14、为节点节点，可以是一个单元设备，也可以是一，可以是一个单元设备，也可以是一个个虚拟虚拟设备。设备。结构单元和物流分别进行编号，以便对系统结构进行数学描述。结构单元和物流分别进行编号，以便对系统结构进行数学描述。其目的是用于对工艺流程图进行分析，以便计算机求解。其目的是用于对工艺流程图进行分析，以便计算机求解。信息流程图与工艺流程图不一定是完全一一对应的（比如信息流信息流程图与工艺流程图不一定是完全一一对应的（比如信息流图中的图中的混合块混合块( (分流块）分流块），在生产现场可能就是一个，在生产现场可能就是一个三通三通，而某些，而某些次要的辅助设备，比如流程图中的次要的辅助设备，比如流程图中的

15、泵泵，在流程图中绘出，在信息流，在流程图中绘出，在信息流图中却不一定出现。图中却不一定出现。 152.5.5.2 系统结构模型的矩阵表达方式系统结构模型的矩阵表达方式 过程系统的模型化是指用数学方法来描述系统结构，必须保过程系统的模型化是指用数学方法来描述系统结构，必须保证从模型返推回流程的唯一性。通常用矩阵来表示系统结构（也证从模型返推回流程的唯一性。通常用矩阵来表示系统结构（也就是说用矩阵表示单元与单元、元与物流之间的关系）。常用的就是说用矩阵表示单元与单元、元与物流之间的关系）。常用的四种系统结构的矩阵表达方式四种系统结构的矩阵表达方式： 过程矩阵过程矩阵(Process Matrix)

16、 流线联系矩阵流线联系矩阵(Stream connection Matrix) 关联矩阵关联矩阵(Incidence Matrix) 邻接矩阵邻接矩阵(Adjacency Matrix)161 )过程矩阵过程矩阵(Process Matrix)-Rp 在过程矩阵中，矩阵行的序号与结构单元的序号对应，在过程矩阵中，矩阵行的序号与结构单元的序号对应，各行中矩阵元素的数值为与该单元有关的物流号，并规定流各行中矩阵元素的数值为与该单元有关的物流号，并规定流入该单元的物流取正值，流出物流取负值。有关物流的先后入该单元的物流取正值，流出物流取负值。有关物流的先后次序是任意的。次序是任意的。矩阵的行数对应模

17、块数；列数对应单元模块矩阵的行数对应模块数；列数对应单元模块的最大输入输出流线数。的最大输入输出流线数。12312345以以单单元元为为主主17结构单元号结构单元号 计算时名称计算时名称 有关物流号有关物流号 1 MIXER 1 5 -2 2 REAO 2 -3 3 DISTIL 3 -4 -5相应的过程矩阵为相应的过程矩阵为Rp = 1 5 -2 2 -3 3 -4 -5182）流线联系矩阵流线联系矩阵(Stream connection Matrix)-Rs 在流线联系矩阵中，每条流线对应在流线联系矩阵中，每条流线对应一行，一行，有几条流线就有几行有几条流线就有几行（总流线数总流线数=行数

18、行数）。而矩阵只有）。而矩阵只有两列（两列（from单元和单元和to单元）单元）。每。每行中的列元素对应着该流线的来源行中的列元素对应着该流线的来源模块模块（从模块）和去向模块（从模块）和去向模块(到模到模块块)。 如果来源于系统以外或输出至系统外，则该列元素就取为如果来源于系统以外或输出至系统外，则该列元素就取为0（意义？）（意义？）。12312345Mix1 react1 sep1以以物物流流为为主主19 行号行号 从单元从单元 到单元到单元 1 混合混合 2 混合混合 反应器反应器 3 反应器反应器 分离器分离器 4 分离器分离器 5 分离器分离器 混合器混合器 相应的流线联系矩阵为相应

19、的流线联系矩阵为Rs= 0 mix01 mix01 react1 react1 sept01 sept01 0 sept01 mix01Note：过程矩阵与：过程矩阵与流线联系矩阵反映流线联系矩阵反映的信息是相同的！的信息是相同的！203) 邻接矩阵邻接矩阵(Adjacency Matrix)-Ra表达表达两个单元两个单元之间的关系之间的关系-或相互联接或不联接或相互联接或不联接；一个由一个由n个单个单元组成的系统的邻接矩阵是一个元组成的系统的邻接矩阵是一个nn方阵方阵。邻接矩阵中行和列的邻接矩阵中行和列的序号分别代表单元的序号。矩阵序号分别代表单元的序号。矩阵元素的数值由单元之间的连接情元素

20、的数值由单元之间的连接情况决定况决定。“行行”代表来源节点（代表来源节点（from单元），单元），“列列”代表物流目代表物流目标节点（标节点（to 单元）。单元）。Rij=1，有表示物流从单元，有表示物流从单元i流向单元流向单元jRij=0，无表示物流从单元，无表示物流从单元i流向单元流向单元j12312345两两个个单单元元的的关关系系21224) 关联矩阵关联矩阵(Incidence Matrix) 表示单元与物流的关系表示单元与物流的关系。在关联矩阵中。在关联矩阵中，行序号与结构单元号行序号与结构单元号对应（行数对应（行数=节点数）而列序号与物流号对应（列数节点数）而列序号与物流号对应（

21、列数=物流数）。物流数）。矩阵中元素取值为逻辑值矩阵中元素取值为逻辑值1，-1和和0；矩阵中矩阵中i行行j列元素的值若为列元素的值若为1，表示物流表示物流j流入单元流入单元i;若为若为1，表示，表示j从单元从单元i流出；若为流出；若为0，表示，表示物流物流j与单元与单元i没有联系；没有联系；当当j列总合为零，且列总合为零，且+1在在-1之前出现，则之前出现，则表示该流股是表示该流股是（？）。（？）。从关联矩阵中可以看出输入流、输出流及其循环流股从关联矩阵中可以看出输入流、输出流及其循环流股。12312345单单元元与与物物流流有有无无关关联联2324共同点共同点：物理意义相同，都是系统结构的数

22、学模型，与一定系统：物理意义相同，都是系统结构的数学模型，与一定系统 结构相对结构相对 应：表明了系统有哪些单元构成，这些单元是怎样联系的；应：表明了系统有哪些单元构成，这些单元是怎样联系的；不同点不同点：形式不同，按不同的需要建立，适应系统模型的不同求解方法。：形式不同，按不同的需要建立，适应系统模型的不同求解方法。 邻接矩阵比较简单，适合矩阵运算；邻接矩阵比较简单，适合矩阵运算； 过程矩阵和流线联系矩阵适合做顺序计算的需要；过程矩阵和流线联系矩阵适合做顺序计算的需要； 从关联矩阵中可以容易地看出输入流、输出流及其循环流股。从关联矩阵中可以容易地看出输入流、输出流及其循环流股。 四四种系统结

23、构的数学模型评述种系统结构的数学模型评述 25例题例题 试建立图试建立图23a所示精馏塔系统的结构模型，请写出系统的所示精馏塔系统的结构模型，请写出系统的过程矩阵、邻接矩阵、关联矩阵过程矩阵、邻接矩阵、关联矩阵。26过程矩阵、邻接矩阵和关联矩阵27 过程矩阵的计算机程序实现过程矩阵的计算机程序实现 过程矩阵在稳态流程模拟软件中是最重要和常用的流程结构表过程矩阵在稳态流程模拟软件中是最重要和常用的流程结构表达方式。基于过程矩阵，计算机程序就可达方式。基于过程矩阵，计算机程序就可“自动自动”地完成全流程中地完成全流程中各个部分的计算。其实现步骤大致为：各个部分的计算。其实现步骤大致为：( (对照对

24、照p32-33p32-33）（1 1）用过程矩阵表达流程结构）用过程矩阵表达流程结构; ;（将矩阵输入到计算机中，或由计将矩阵输入到计算机中，或由计 算机自动识别算机自动识别）（2 2）逐行搜索过程矩阵，识别当前所处理的设备（模块）逐行搜索过程矩阵，识别当前所处理的设备（模块）; ;（3 3）识别当前设备的输入和输出物流）识别当前设备的输入和输出物流; ;（4 4）由当前设备的输入物流状态（已知或已赋值）和设备参数计）由当前设备的输入物流状态（已知或已赋值）和设备参数计算设备的运行状态算设备的运行状态; ;（5 5）根据设备状态给相应的输出物流状态赋值。当过程矩阵的每）根据设备状态给相应的输出物流状态赋值。当过程矩阵的每一行都处理完毕，整个流程的计算就告一段落