第9讲MATLAB机理建模方法1推导法机理建模.pdf

第9 讲 MATLAB 机理建模方法 作者：卓金武, MathWorks 中国 在数学建模中，如果遇到一个非典型的数学建模问题(非数据、优化、连续、评价)， 那么这种情况下，通常需要用到机理建模方法了。 机理建模就是根据对现实对象特性的认识，分析其因果关系，找出反映内部机理的规 则，然后建立规则的数学模型。 机理建模的经典案例有很多，比如万有引力公式的推导过 程。机理建模常见的有两类，一类是推导法机理建模，类似于微分方程建模，常用于动力学 的建模过程，比如化学中反应动力学，还有各种场的方程，比如压力场、热场方程等；一类 是包含一个或几个类别对象的复杂系统问题，常通过元胞自动机-仿真法来进行机理建模。 下面将介绍这两类机理建模的具体MATLB 实现过程。 1. 推导法机理建模 1.1 问题 某种医用薄膜有允许一种物质的分子穿透它(从高浓度的溶液向低浓度的溶液扩散) 的 功能，在试制时需测定薄膜被这种分子穿透的能力。测定方法如下：用面积为S 的薄膜将容 器分成体积分别为 、 的两部分，在两部分中分别注满该物质的两种不同浓度的溶液。 V V A B 此时该物质分子就会从高浓度溶液穿过薄膜向低浓度溶液中扩散。已知通过单位面积薄膜分 K 子扩散的速度与膜两侧溶液的浓度差成正比，比例系数 表征了薄膜被该物质分子穿透的 能力，称为渗透率。定时测量容器中薄膜某一侧的溶液浓度值，可以确定K 的数值，试用 数学建模的方法解决K 值的求解问题。 1.2 假设和符号说明 为了便于建模，作以下几点假设： ① 薄膜两侧的溶液始终是均匀的，即在任何时刻膜两侧的每一处溶液的浓度都是相同 的。 ② 当两侧浓度不一致时，物质的分子穿透薄膜总是从高浓度溶液向低浓度溶液扩散。 ③ 通过单位面积膜分子扩散的速度与膜两侧溶液的浓度差成正比。 ④ 薄膜是双向同性的即物质从膜的任何一侧向另一侧渗透的性能是相同的。 图1 圆柱体容器被薄膜截面S 阻隔 同时，约定需要用到的几个数学符号： ① t CA (t) 、CB (t) 表示 时刻膜两侧溶液的浓度; ② 、 表示初始时刻两侧溶液的浓度 (单位:毫克/立方厘米); a a A B ③ K 表示渗透率; V V A B ④ 、 表示由薄膜阻隔的容器两侧的体积; 1.3 模型的建立 考察时段[t t t] 薄膜两侧容器中该物质质量的变化。以容器A 侧为例，在该时段物 ， +∆ 质质量增加量： V C (t t) V C (t) +∆ − A A A A 另一方面由渗透率的定义我们知道，从B 侧渗透至A 侧的该物质的质量为： ( − )∆ SK C C t B A 由质量守恒定律，两者应该相等，于是有： V C (t +∆t) −V C (t) SK (C −C )∆t A A A A B A 两边除以∆t ，令∆t →0 并整理得 dC SK A (CB −CA ) dt V A (3-1) 且注意到整个容器的溶液中含有该物质的质量应该不变，即有下式成立： V C (t) +V C (t) V a +V a A A B B A A B B V V ( ) B B ( ) C t a + a − C t A A