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

PAGE

PAGE1

湖水的自我净化问题

摘要

本题是一容积为的大湖受到某种物质污染，从某时刻起污染源被切断，湖水开始更新，更新速率为，建立求污染物浓度下降至原来的需要多长时间的数学模型问题。解决本问题需要用到微元法的思想，也就是在很小的时间内流出的湖水污染物浓度不变，然后利用湖水污染物的变化量等于流出湖水的污染量建立等式关系，对该等式求导后得出一个微分方程，利用Matlab中dsolve函数解该微分方程，求得污染物浓度下降至原来的所需时间为天。本模型涉及到解微分方程，所以模型的应用很广泛，可以应用到动态分析问题中，利用该模型可以解决大量实际生活和生产问题。

关键词：微元法；微分方程；动态分析；Matlab

一、问题重述

1.1背景资料与条件

有一容积为（单位：）的大湖受到某种物质的污染，污染物均匀的分布在湖中。若从某时刻起污染源被切断，设湖水更新的速率是（单位：）。试建立求污染物浓度下降至原来的需要多长时间的数学模型。

1.2需要解决的问题

在湖的容积为，湖水更新速率为的条件下，求污染终止后，污染物下降到原来的所需的时间。

二、基本假设

2.1模型的假设

假设一：湖水保持体积不变。

假设二：污染物始终均匀的分布在湖中。（假设合理性见背景资料与条件。）

假设三：在很小的时间内污染物浓度不变。（微元法思想）

2.2本文引用数据、资料均真实可靠。

三、符号说明

3.1模型的符号说明

A:

B:

C:为污染源切断后湖水更新的时间(单位：天)。

四、模型的建立与求解

4.1模型的建立

从开始到天内湖水含污染物改变量为：

由于流入湖中的水没有污染物，所以天内更新流出污染物量为：

利用湖水污染物的变化量=流出湖水的污染量得：

对求导得微分方程为：

8.2.2模型的缺点

由于实际结果不好求得，故不好进行误差分析，这是本模型的主要缺点。因为不好进行误差分析，所以也就不知道所求得的结果的好坏，进一步来讲也就不清楚所建立的模型的好坏，对模型的评价就不好进行。

8.3模型的优化

因为模型存在舍入误差，要想求得的结果精确，使得模型有意义，那么数字的选取就要特别注意，数据精度就要高，要选取精度较高的数据，这样得出的结果就精确，模型也就得到优化。

参考文献

[1]赵静,但琦.数学建模与数学实验[M].北京：高等教育出版社，2008.

[2]胡运权.运筹学教程(第三版)[M].北京：清华大学出版社，2007.

[3]白峰杉.数值计算引论[M].北京：高等教育出版社，2010.

[4]陈冬彦，李冬梅，王树忠.数学建模[M].北京：科学出版社，2007.

[5]萧数铁.数学实验[M].北京:高等教育出版社，1999.

附录

Ask1.m

w=dsolve(Dw=-r*w/V,t)

Ask2.m

r=4.121*10^10;

V=5.176*10^12;

t=-V*log(0.03)/r

Ask3.m

r=4.121*10^10;

V=5.176*10^12;

w=100;

t=0:1000;

w=exp(-r/V*t)*w;

plot(t,w)

holdon

w=3

w0=100

t=-V*log(w/w0)/r;

t

plot(t,w,r*)