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

三、符号说明 湖水的自我净化问题 摘 要 此题是一容积为 V 的大湖受到某种物质污染， 从某时刻起污染源被切断， 湖水开始更新， 更新速率为 r ，建立求污染物浓度下降至原来的 3% 需要多长时间的数学模型问题。解 决本问题需要用到微元法的思想，也就是在很小的时间流出的湖水污染物浓度不变，然 后利用湖水污染物的变化量等于流出湖水的污染量建立等式关系， 对该等式求导后得出 一个微分方程，利用 Matlab 中 dsolve 函数解该微分方程，求得污染物浓度下降至原来 的3% 所需时间为 440.4 天。 本模型涉与到解微分方程，所以模型的应用很广泛， 可以应 用到动态分析问题中，利用该模型可以解决大量实际生活和生产问题。 关键词：微元法；微分方程；动态分析； Matlab 一、问题重述 1.1 背景资料与条件 有一容积为 V （单位： m3 ）的大湖受到某种物质的污染，污染物均匀的分布在湖 中。若从某时刻起污染源被切断，设湖水更新的速率是 染物浓度下降至原来的 3% 需要多长时间的数学模型。 1.2 需要解决的问题 在湖的容积为 5.176*10^12( m3 ) ，湖水更新速率为 求污染终止后，污染物下降到原来的 3% 所需的时间。 二、基本假设 2.1 模型的假设 r （单位： m3 / d ）。试建立求污 4.121*10^10( m3 / d) 的条件下， 1) 假设一：湖水保持体积 V 不变。 2) 假设二：污染物始终均匀的分布在湖中。（假设合理性见背景资料与条件。） 3) 假设三：在很小的时间污染物浓度不变。（微元法思想） 2.2 本文引用数据、资料均真实可靠。 t0 t 0 dt V t . . . 3.1 模型的符号说明 A: w(t) : t时刻湖区的污染物浓度。 B: w(0):表示初始时刻湖中水的污染浓度。 C: t为污染源切断后湖水更新的时间 (单位：天 )。 四、模型的建立与求解 4.1 模型的建立 从开始到 t天湖水含污染物改变量为： Vw(0) Vw(t ) 由于流入湖中的水没有污染物，所以 t 天更新流出污染物量为： 0 w(t )rdt 利用湖水污染物的变化量 对t求导得微分方程为： 变换后可得： =流出湖水的污染量得： Vw(0) Vw(t) 0 w(t )rdt Vdw(t) dt w(t )r , dw(t ) dt ，w(t )r ， V 然后利用 Matlab 中 dsolve 函数求解微分方程，代入 w(t) 3%w(0) 求得时间 t。 4.2 模型的概述 此题是利用微元法的思想， 通过利用湖水污染物的变化量等于流出湖水的污染量列 出等式，然后求导变成微分方程，接着通过 Matlab 中 dsolve 函数求解微分方程，最后 通过代入条件求得时间 t。 4.3 模型的运用与求解 由Vw(0) Vw(t) t w(t )rdt 对t 求导后得 dw(t ) w(t )r , 再利用 Matlab 中 dsolve 函数求解此微分方程， Matlab 运行后得 w = C2/exp((r * t)/V) ，在 Matlab 中编写程序求 解t ，代入条件求得结果，对湖水的自我净化过程作图（如图 1） （程序见附录） 。 1 / 6 . . . 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 模型的结果 将 r 4.121* 10^10 ， V 5.176 * 10^12 代入得 t 440.4257 ，保留一位小数后求得