第6章 医用薄膜渗透率的.ppt

1 第 6 章 医用薄膜渗透率的确定 —— 曲线拟合 2 概述 ? 利用数据确定方程的参数 ? 如：刹车距离与刹车速度的有如下关系 d=Cv 2 其中， d 表示刹车距离， v 表示刹车时的速度，需要确定 常数 C 。 ? 根据一组二维数据，即平面上的若干点，要求确定 一个一元函数 y=f(x) ，使这些点与曲线总体来说 尽量接近，这就是数据拟合曲线的基本思想，简称 曲线拟合 (fitting a curve) 。 3 6.1 医用薄膜的渗透率 测定薄膜渗透率 V A V B S 用面积为 S 的薄膜将容器分成 体积分别为 V A ， V B 的两部分， 在两部分中分别注满两种不同 浓度的溶液，此时物质分子会 从高浓度溶液穿过薄膜向低浓 度溶液扩散，通过单位面积薄 膜分子扩散的速度与薄膜两侧 溶液的浓度差成正比，试确定 薄膜渗透率 K 。 4 6.2 确定医用薄膜渗透率的数学模型 ? 问题假设 ? 薄膜两侧的溶液始终是均匀的，即在任何时刻薄膜两侧 的每一处溶液的浓度都是相同的 ? 当薄膜两侧的溶液浓度不一致时，物质的分子穿透薄膜 总是从高浓度溶液向低浓度溶液扩散 ? 通过单位面积薄膜分子扩散的速度与薄膜两侧溶液的浓 度差成正比。 ? 薄膜是双向同性的，即物质从薄膜的任何的一侧向另一 侧渗透的性能是相同的 5 ? 符号说明 ? C A (t),C B (t) 表示 t 时刻薄膜两侧溶液的浓度 ? α A 、 α B 表示初始时刻两侧溶液的浓度（ mg/cm 3 ） ? K 表示薄膜渗透率 ? V A ， V B 表示由薄膜阻隔的容器两侧的体积 6.2 确定医用薄膜渗透率的数学模型 6 6.2 确定医用薄膜渗透率的数学模型 问题分析 解决问题的思路：首先通过机理分析寻找某一侧浓度随 时间变化的函数关系 C A (t) 或 C B (t) ，其中可能含有待定参 数，如薄膜渗透率 K ，然后根据一组测量值（ t i ,C i ）， i=1,2,……n 去确定模型中待定参数。 [t,t+ △ t] 时段容器 A 侧的物质质量增加量 从 B 侧渗透至 A 侧的质量为 7 6.2 确定医用薄膜渗透率的数学模型 由质量守恒定律，有 = 两边除以 △ t ，令 △ t → 0 并整理得 整个容器的溶液中物质的质量不变，有 (6.1) 8 两边除以 V A ，并整理得 6.2 确定医用薄膜渗透率的数学模型 代入 (6.1) 得 再利用初始条件 解出 9 6.2 确定医用薄膜渗透率的数学模型 数学模型 问题归结为利用 C B 在时刻 t j 的测量数据 C j (j=1,2,…,N) 来辨识参 数 K 和 α A ， α B ，对应的数学模型变为使函数 达到最小 令 问题转化为如下最优化问题 10 最小二乘法 给定平面上的点 (xi,yi),i=1,2,…n ，求 f(x) ，使 达到最小。 其中 δ i 为点（ x i ,y i ）与曲线 y=f(x) 的距离。曲线拟合的实 际含义是寻求一个函数 y=f(x) ，使 f(x) 在某种准则下与所 有数据点最为接近，即曲线拟合得最好。 6.3 一元最小二乘法简介 11 6.3 一元最小二乘法简介 设拟合函数可由一些简单的“基函数”（如幂函数，三角 函 数等） 来线性表示： 现在要确定系数 c 0 ,c 1 ,…c m , 使 δ 达到最小 12 6.3 一元最小二乘法简介 曲线拟合示意图 13 6.3 一元最小二乘法简介 已知一组数据 (xi,yi),i=1,2,…n ，如何选择函数 f(x)? ? 一是根据机理分析来确定函数形式 ? 二是根据散点图直观判断函数 f(x) 的形式，常用的一 元曲线拟合函数有 双曲线 指数曲线 多项式 f(x)=a+b/x f(x)=ae bx 14 6.4 用曲线拟合方法确定医用薄膜渗透率 医用薄膜渗透率最小二乘法辨识模型 假设： V A =V B =1000cm 3 ， S=10cm 2 ， B 部分溶液浓度测试结 果如下表： t j /s 100 200 300 400 500 600 700 800 900 1000 C j (10 -5 ) 4.54 4.99 5.35 5.65 5.90 6.10 6.26 6.39 6.50 6.59 (6.2) 15 6.4 用曲线拟合方法确定医用薄膜渗透率 ? ? ? ? ? ? 10 1 2 02 . 0 ] [ ) , , ( min j j Kt C be a b a K E j 将 6.2 简化为 1) 编写 M 文件 (nongdu.m) function f=nongdu(x,tdata) f=x(1)+x(2)*exp(-0.02