数学建模-薄膜分子穿透力.doc

分子穿透能力的测定摘 要通过对问题的分析，根据质量守恒，利用微分方程模型，得到了关于浓度低一侧浓度对时间的微分方程模型，通过求解参数和简化以后确定了浓度与时间的指数关系。运用MATLAB编程进行拟合，求得参数，从而得到渗透率。最后，对拟合结果进行检验，检验结果见图5.2.1和图5.3.1。关键字：渗透率 MATLAB软件 参数估计 微分方程模型1、 问题重述某种医用薄膜有允许一种物质的分子穿透它从高浓度的溶液向低浓度的溶液扩散的功能，在测试时需测定薄膜被这种分子穿透的能力。测定方法如下：用面积10cm2的薄膜将分成体积分别为100cm3和100cm3的两部分，在两部分中分别注满该物质的两种不同浓度的溶液。此时该物质分子就会从高浓度溶液穿过薄膜向低浓度溶液中扩散。通过单位面积膜分子扩散的速度与膜两侧溶液的浓度差成正比，比例系数表征了薄膜被该物质分子穿透的能力，称为渗透率。定时测量容器中薄膜某一侧的溶液浓度值，以此确定的数值。对容器一侧溶液浓度的测试结果如下：tj1002003004005006007008009001000Cj(10-3mg/cm3)4.544.995.355.655.906.106.266.396.506.59试建立一个较好的数学模型并给出相应的算法和程序。二、问题分析用单位体积溶液中所含的溶质质量数来表示的浓度叫质量-体积浓度，题中给出的浓度为质量-体积浓度。通过薄膜单位面积分子扩散速度与膜两侧溶液浓度差成正比例，比例系数被称为渗透率，它表征了薄膜被分子穿透的能力。要确定渗透率，需要建立通过薄膜单位面积分子扩散速与薄膜两侧浓度差的关系模型。在本题中我们可以根据质量守恒来进行求解，考察时段薄膜两侧容器中该物质质量的变化，可以用两种形式表示出来，薄膜的一侧在时段内物质质量的增加，以及根据渗透率表示的一侧渗透至另一侧的物质的质量，两者相等可得到等式（1）。其次，整个容器的溶液中含有该物质的质量不变,与初始时刻该物质的含量相同，也是质量守恒，可得到等式（2），即联立两个等式，可以得到浓度与时间的表达式，代入数据，对等式的系数进行求解。3、 模型假设1、薄膜两侧的溶液始终是均匀的，即在任何时刻薄膜两侧的每一处溶液的 溶度都是相同的；2、当两侧浓度不一致时，物质的分子穿透薄膜总是从高浓度溶液向低浓度 溶液扩散；3、通过单位面积膜分子扩散的速度与膜两侧溶液的浓度差成正比；4、薄膜是双向性的，即物质从膜的任何一侧向另一侧渗透的性能是相同的。四、符号定义与说明 符号 定义与说明 渗透率 t时刻薄膜的A侧溶液浓度 t时刻薄膜的B侧溶液浓度 薄膜A侧的体积 薄膜B侧的体积 薄膜的面积 初始时刻A侧的溶液浓度 初始时刻B侧的溶液浓度5、 模型的建立与求解5.1 模型的建立考察时段薄膜两侧容器中该物质质量的变化。 在容器的一侧，物质质量的增加是由于另一侧的物质向该侧渗透的结果，因此物质质量的增量应等于另一侧的该物质向这侧的渗透量。设、分别表示在时刻t膜两侧溶液的浓度，浓度单位:以容器A侧为例，在时段物质质量的增量为：由于平均每单位时间通过单位面积薄膜的物质分子量与膜两侧溶液的浓度差成正比，比例系数为K。因此，在时段，从B侧渗透至A侧的该物质的质量为：于是有：两边除以，并令取极限再稍加整理即得： （1）注意到整个容器的溶液中含有该物质的质量不变,与初始时刻该物质的含量相同，因此 其中、分别表示在初始时刻两侧溶液的浓度从而：代入式（1）得：解得： 由上面的分析得到的表达式后问题归结为利用CB在时刻tj的测量数据Cj(j=1,2,.,N)来辨识 K 和 引入： 从而将代入上式有：5.2 模型的求解由5.1可知和的关系式，可用该关系式来拟合附件中给的定时测量容器中薄膜某一侧的溶液浓度值的实验数据，即确定参数的值。利用MATLAB软件编程解得 容器一侧的溶液浓度离散点与拟合结果图见图5.2.1。图5.2.1容器一侧的溶液浓度离散点与拟合结果图5.3 拟合检验运用MATLAB软件对拟合的结果进行检验，检验的结果见图5.3.1拟合检验结果截图。图5.3.1拟合检验结果截图其中Sse表示误差平方和,值越小，表示拟合的曲线和真实数据越接近Rsquare表示确定系数,值越接近于1越好dfe表示自由度即能够自由取值的变量个数 Adjrsquare表示调整后的确定系数,值越接近于1越好Rmse表示均方根误差,值越小，表示拟合的曲线和真实数据越接近由结果可知sse和rmse两者的值都十分接近于一，即表示拟合的曲线和真实数据十分接近，并且确定系数和调整后的确定系数两者的值都接近于1，即综合表示拟合结果较好。六、模型评价及推广本文主要根据质量守恒，利用微分方程模型，得到了关于浓度低一侧浓度对时间的微分方程模型，通过求解参数和简化以后确定了浓度与时间的指数关系，运用MATLAB编程进行拟合，求得参数，从而得到渗透率。模型较为简单，易于实现，通用性强。七、参考文献1 章绍辉，数学建模.北京；科学出版社，2010第一版.2 袁振东、洪源，数学建模.华东师范大学出版社，2002第一版.8、 附件MATLAB程序clc;clear all;t= 100 200 300 400 500 600 700 800 900 1000;t=t;c=1e-03*4.54,4.99,5.35,5.65,5.90,6.10,6.26,6.39,6.50,6.59;c=c;st=0,0.05,0.05;ft_ = fittype(a+b*exp(-0.02*K*t),dependent,c,independent,t, coefficients,a, b,