数学建模--湖水的自我净化问题.doc

湖水的自我净化问题 摘 要 本题是一容积为 V 的大湖受到某种物质污染，从某时刻起污染源被切断，湖水开始 更新，更新速率为 r ，建立求污染物浓度下降至原来的 3% 需要多长时间的数学模型问 题。解决本问题需要用到微元法的思想，也就是在很小的时间内流出的湖水污染物浓度 不变，然后利用湖水污染物的变化量等于流出湖水的污染量建立等式关系，对该等式求 导后得出一个微分方程，利用 Matlab 中 dsolve 函数解该微分方程，求得污染物浓度下 降至原来的 3% 所需时间为 440.4 天。本模型涉及到解微分方程， 所以模型的应用很广泛， 可以应用到动态分析问题中，利用该模型可以解决大量实际生活和生产问题。 关键词 ：微元法；微分方程；动态分析； Matlab 一、问题重述 1.1 背景资料与条件 有一容积为 V （单位： 3 m ）的大湖受到某种物质的污染，污染物均匀的分布在湖 中。若从某时刻起污染源被切断，设湖水更新的速率是 r （单位： m d ）。试建立求污 3 / 3 / 染物浓度下降至原来的 3% 需要多长时间的数学模型。 1.2 需要解决的问题 在湖的容积为 3 5.176*10^12( m ) ，湖水更新速率为 3 4.121*10^10( m / d) 的条件下， 求污染终止后，污染物下降到原来的 3% 所需的时间。 二、基本假设 2.1 模型的假设 1) 假设一：湖水保持体积 V 不变。 2) 假设二：污染物始终均匀的分布在湖中。（假设合理性见背景资料与条件。） 3) 假设三：在很小的时间内污染物浓度不变。（微元法思想） 2.2 本文引用数据、资料均真实可靠。 三、符号说明 3.1 模型的符号说明 A: w(t) : t时刻湖区的污染物浓度。 B: w(0) :表示初始时刻湖中水的污染浓度。 C: t 为污染源切断后湖水更新的时间 ( 单位：天) 。 四、模型的建立与求解 4.1 模型的建立 从开始到 t天内湖水含污染物改变量为： Vw (0) Vw(t ) 由于流入湖中的水没有污染物，所以 t 天内更新流出污染物量为： t 0 w(t )rdt 利用湖水污染物的变化量 =流出湖水的污染量得： t Vw(0) Vw(t) w(t )rdt 0 对t求导得微分方程为： 1 Vdw dt (t) w(t )r , 变换后可得： dw(t ) w(t )r dt V ， 然后利用 Matlab 中 dsolve 函数求解微分方程，代入 w(t) 3% w(0) 求得时间 t 。 4.2 模型的概述 本题是利用微元法的思想， 通过利用湖水污染物的变化量等于流出湖水的污染量列 出等式，然后求导变成微分方程，接着通过 Matlab 中 dsolve 函数求解微分方程，最后 通过代入条件求得时间 t 。 4.3 模型的运用与求解 t Vw(0) Vw( t) w(t )rdt 对t 求导后得 由 0 dw(t ) w(t )r dt V , 再利用 Matlab 中 dsolve 函数求解此微分方程， Matlab 运行后得 w = C2/exp((r * t)/V) ，在 Matlab 中编写程序求 解t，代入条件求得结果，对湖水的自我净化过程作图（如图 1）（程序见附录）。 100 90 80 70 60 50 40 30 20 10 0 0 100 200 300 400 500 600 700 800 900 1000 图 1 其中横坐标为时间 t（单位：天），纵坐标为 w(t) 关于 w(0) 的百分比。红色星号表示 污染浓度降到原来的 3% 时的点。从图中可以看出湖水自我净化速率呈下降趋势。 4.4 模型的结果 2 将r 4.121 * 10^10 ，V 5.176 * 10^12 代入得 t 440.4257 ，保留一位小数后求得 t 440.4。即污染物浓度下降至原来的 3% 需要的时间为 440.4 天。 五、模型的分析 5.1 假设的合理性分析 如果湖水体积变化，那么题目就没法做了，因此这个假设是必要的且是合理的。污 染物始终均匀的分布在湖中，题目条件中已给出，所以此假设合理可靠。在很小的时间 内污染物浓度不变，这是利用微元法的思想，故假设的合理性毋庸置疑。 5.2 模型的误差分析 本模型的误差主要在数字的处理上， 即保留几位的问题上， 也就是说存在舍入误差， 本题在最后结果中保留了一位小数。 六、模型的检验 由于实际数据不好求得，所以用模型计算求得的数据不好与实际数据进行比对，故 也就计算不了误差大小。但本模型对所取得数据精度不同，产生的结果也就不同。根据 模型的分析得本模型误差产生的原因主要是数据精度，其它因素影响很小，在本题中我 所取的精度所得出的结果产生的误差在模型估计和实际许可的范围之内。 七、模