1、分子穿透能力的测定摘 要通过对问题的分析，根据质量守恒，利用微分方程模型，得到了关于浓度低一侧浓度对时间的微分方程模型，通过求解参数和简化以后确定了浓度与时间的指数关系。运用 MATLAB 编程进行拟合，求得参数，从而得到渗透率 。K最后，对拟合结果进行检验，检验结果见图 5.2.1 和图 5.3.1。关键字：渗透率 MATLAB 软件 参数估计 微分方程模型1、问题重述某种医用薄膜有允许一种物质的分子穿透它从高浓度的溶液向低浓度的溶液扩散的功能，在测试时需测定薄膜被这种分子穿透的能力。测定方法如下：用面积 10cm2 的薄膜将分成体积分别为 100cm3 和 100cm3 的两部分，在两部分
2、中分别注满该物质的两种不同浓度的溶液。此时该物质分子就会从高浓度溶液穿过薄膜向低浓度溶液中扩散。通过单位面积膜分子扩散的速度与膜两侧溶液的浓度差成正比，比例系数 表征了薄膜被该物质分子穿透的能力，称为渗K透率。定时测量容器中薄膜某一侧的溶液浓度值，以此确定 的数值。K对容器一侧溶液浓度的测试结果如下：tj 100 200 300 400 500 600 700 800 900 1000Cj(10-3mg/cm3) 4.54 4.99 5.35 5.65 5.90 6.10 6.26 6.39 6.50 6.59试建立一个较好的数学模型并给出相应的算法和程序。二、问题分析用单位体积溶液中所含的溶
3、质质量数来表示的浓度叫质量-体积浓度，题中给出的浓度为质量-体积浓度。通过薄膜单位面积分子扩散速度与膜两侧溶液浓度差成正比例，比例系数 被称为渗透率，它表征了薄膜被分子穿透的能力。K要确定渗透率，需要建立通过薄膜单位面积分子扩散速与薄膜两侧浓度差的关系模型。在本题中我们可以根据质量守恒来进行求解，考察时段 薄膜两侧容t,器中该物质质量的变化，可以用两种形式表示出来，薄膜的一侧在时段内物质质量的增加，以及 根据渗透率 表示的一侧渗透至另一侧的物t, K质的质量，两者相等可得到等式（1） 。其次，整个容器的溶液中含有该物质的质量不变,与初始时刻该物质的含量相同，也是质量守恒，可得到等式（2），即联
4、立两个等式，可以得到浓度与时间的表达式，代入数据，对等式的系数进行求解。3、模型假设1、薄膜两侧的溶液始终是均匀的，即在任何时刻薄膜两侧的每一处溶液的 溶度都是相同的； 2、当两侧浓度不一致时，物质的分子穿透薄膜总是从高浓度溶液向低浓度 溶液扩散；3、通过单位面积膜分子扩散的速度与膜两侧溶液的浓度差成正比；4、薄膜是双向性的，即物质从膜的任何一侧向另一侧渗透的性能是相同的。四、符号定义与说明符号 定义与说明渗透率Kt 时刻薄膜的 A 侧溶液浓度)(tCAt 时刻薄膜的 B 侧溶液浓度 B薄膜 A 侧的体积AV薄膜 B 侧的体积B薄膜的面积S初始时刻 A 侧的溶液浓度)0(AC初始时刻 B 侧的
5、溶液浓度B5、模型的建立与求解5.1 模型的建立考察时段 薄膜两侧容器中该物质质量的变化。 t,在容器的一侧，物质质量的增加是由于另一侧的物质向该侧渗透的结果，因此物质质量的增量应等于另一侧的该物质向这侧的渗透量。设 、 分别表示在时刻 t 膜两侧溶液的浓度，浓度单位:)(tCAtB 3/cmg以容器 A 侧为例，在时段 物质质量的增量为：,)()(tCVtAA由于平均每单位时间通过单位面积薄膜的物质分子量与膜两侧溶液的浓度差成正比，比例系数为 K。因此，在时段 ，从 B 侧渗透至 A 侧的该物质的质量为：t, tSCt)(于是有： ttKtVtCABAA)()(两边除以 ，并令 取极限再稍加
6、整理即得：t0ABAttSdt)()（1）注意到整个容器的溶液中含有该物质的质量不变,与初始时刻该物质的含量相同，因此 )0()()()( BABA CVtCVt 其中 、 分别表示在初始时刻两侧溶液的浓度)0(AC(B从而： ABABAVtt )()0()(代入式（1）得： )0()()1() BABBAB CSKtCVSKdtC解得：tVSKBAABAB BAet )1()0()0)( 由上面的分析得到 的表达式后问题归结为利用 CB 在时刻 tj 的测量数()Ct据 Cj(j=1,2,.,N)来辨识 K 和 ,AB引入： BAVa)0()( BAAVCb)0(从而 tVSKBBAeatC
7、)1()(将 代入上式有：3210,10ABVcmS0.2()KtBCtabe5.2 模型的求解由 5.1 可知 和 的关系式，可用该关系式来拟合附件中给的定时测量)(tCB容器中薄膜某一侧的溶液浓度值的实验数据，即确定参数 的值。,abK利用 MATLAB 软件编程解得 06941.a02937.b0893.K容器一侧的溶液浓度离散点与拟合结果图见图 5.2.1。图 5.2.1 容器一侧的溶液浓度离散点与拟合结果图5.3 拟合检验运用 MATLAB 软件对拟合的结果进行检验，检验的结果见图 5.3.1 拟合检验结果截图。图 5.3.1 拟合检验结果截图其中Sse 表示误差平方和,值越小，表示
8、拟合的曲线和真实数据越接近Rsquare 表示确定系数,值越接近于 1 越好dfe 表示自由度即能够自由取值的变量个数 Adjrsquare 表示调整后的确定系数,值越接近于 1 越好Rmse 表示均方根误差,值越小，表示拟合的曲线和真实数据越接近由结果可知 sse 和 rmse 两者的值都十分接近于一，即表示拟合的曲线和真实数据十分接近，并且确定系数和调整后的确定系数两者的值都接近于 1，即综合表示拟合结果较好。六、模型评价及推广本文主要根据质量守恒，利用微分方程模型，得到了关于浓度低一侧浓度对时间的微分方程模型，通过求解参数和简化以后确定了浓度与时间的指数关系，运用 MATLAB 编程进行
9、拟合，求得参数，从而得到渗透率 。模型较为K简单，易于实现，通用性强。七、参考文献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, K);%f = fit(t,c,ft_ ,Startpoint,st)f, goodness= fit(t,c,ft_,Startpoint,st)figureplot(f,predobs,0.95);hold on,plot(t,c,b*)
