2005高教社杯全国大学生数学建模竞赛题目A_长江水质的评价和预测[1](1)

1、长江水质的评价和预测摘要文章在已有数据的基础上，建立了水质依靠流量、流速和降解系数的数学模型，找出了污染源的所在地。建立一元线性回归模型，对后十年污水治理做出了预测。利用Matlab，C语言程序进行求解。得出了有关结论。针对问题一，根据03、04年长江流域水质报告表，对长江近两年多的水质情况做出定量的综合评价。对每一个地区在近两年的28个月中的水质情况进行统计，找出该地区污染的种类及该种类污染出现的频率。以此为依据分析各地区水质的污染状况。针对问题二，根据主要污染物在各个观测点的观测数据，建立了水质依靠流量、流速和降解系数的数学模型。对长江干流沿岸各个地段的排污量进行统计，找出了主要污染源所在

2、地区：长江中游湖北宜昌至湖南岳阳段。针对问题三，根据各个年份废水排放量总量，采用一元线性回归模型找出废水排放量总量与年份之间的关系。根据水文年支流和干流的相关数据和各年长江总流量和废水排放量，得出长江总流量中废水排放量的比例，利用matlab对1995-2004年长江总流量中废水排放量的比例拟合（不考虑1998年特大洪水），对未来十年废水排放量占长江总流量比例进行预测。从预测结果中，发现污水百分比呈逐年上升的趋势（从2005年的3.27%到2014年的6.24%），由此说明长江污水的处理迫在眉睫。针对第四问，依照过去10年的类、类水和水的统计数据，通过数据拟合构建了一元线性回归模型、预测的未来

3、十年类、类水和水占长江总水量的百分比。引入流量的概念，得到长江的总水量HS。由治理污水的标准建立分段函数。从而求出未来十年每年需要处理的污水量。针对第五问，提出了解决长江水质污染问题的从四方面着手的方案： 沿江工厂的整治，民众意识的唤醒，上游植被的保护，以及法律的硬性要求。一问题重述2004年10月 “保护长江万里行”考察团，对长江沿线21个重点城市做了实地考察，揭示了一幅长江污染的真实画面。为此，专家们提出“若不及时拯救，长江生态10年内将濒临崩溃” ，并发出了“拿什么拯救癌变长江”的呼唤。通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。江河自身对污染物都有一定的自然

4、净化能力，反映江河自然净化能力的指标称为降解系数。长江干流的自然净化能力可以认为是近似均匀的，根据检测可知，主要污染物高锰酸盐指数和氨氮的降解系数通常介于0.10.5之间，比如可以考虑取0.2(单位：1/天)。请根据，长江沿线17个观测站近两年多主要水质指标的检测数据，干流上个观测站近一年多的基本数据（站点距离、水流量和水流速），“19952004年长江流域水质报告”给出的主要统计数据，以及地表水环境质量标准。研究以下问题：（1）对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。（2）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?（3）假如

5、不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来10年的情况。（4）根据你的预测分析，如果未来10年内每年都要求长江干流的类和类水的比例控制在20%以内，且没有劣类水,那么每年需要处理多少污水？ （5）你对解决长江水质污染问题有什么切实可行的建议和意见。二问题分析 针对问题一，要分两步来解决。第一步，我们根据03、04年长江流域水质报告表，对长江近两年多的水质情况做出定量的综合评价。第二步，分析近两年多长江各地区水质的污染状况。为此我们先对每一个地区在近两年的28个月中的水质情况进行统计，再找出该地区污染的种类及该种类污染出现的频率。

6、第二个问题分为两个部分：一是研究、分析长江干流近一年多主要污染物高锰酸盐污染源主要在哪些地区。二是研究、分析长江干流近一年多主要污染物氨氮的污染源主要在哪些地区。解决这两部分问题所用的方法是一样的。要想找出污染源，就要对长江干流沿岸各个地段的排污量进行统计。排污量高的地段就是主要污染源所在地区。对于问题三：首先筛选取出水文年支流和干流的相关数据和各年长江总流量和废水排放量，对1995-2004年长江总流量中废水排放量的比例拟合，对每年废水排放量占长江总流量比例进行预测。通过预测出的数据分析得出废水排放量占长江总流量比例变化情况，由此对未来十年长江干流和支流污染量的污染量的变化趋势作出总体的分析

7、说明。问题四要求我们预测分析如果未来10年内每年都要求长江干流的类和类水的比例控制在20%以内，且没有劣类水,每年需要处理多少污水。流量除以流速可以得到过水断面的面积。干流处有七个检测点，求出每个检测点的过水断面，求平均值。此值即为长江的平均过水断面。用过水断面的平均值乘以平均河长可以得到长江水的总体积。为了求解需要处理的污水量，我们还要求出需要处理的污水占总水量的百分比。此时问题可分两类进行讨论。第一类，加水所占长江水的百分比未超过20%。此时为了使长江水质达到题目中的要求，只需对劣水进行处理。所需处理的污水量即为长江中劣水的总量，用劣水占长江总水量的百分比乘以长江的总水量即可得到。第二类，

8、、水所占长江水的百分比超过20%。此时为了使长江水质达到题目中的要求需要做两方面的处理。一是将、水进行处理，使其在长江水中所占百分比降至20%。二是对劣水进行处理。 所需处理的污水总量即为所需处理的、水的量与长江中劣水的总量之和。用、水占长江总水量的百分比之和，减去20%再乘以长江的总水量，即可得到所需处理的、水的量。用劣水占长江总水量的百分比乘以长江的总水量即可得到长江中劣水的总量。用一元线性回归模型对、类水所占百分比与年污水排放量，、类水所占百分比与万倍年污水浓度，、类水所占百分比与年污水浓度的倒数，分别进行拟合。选取拟合度最高的模型来预测未来十年长江、类水所占百分比。用一元线性回归模型对

9、、类水所占百分比与年污水排放量，、类水所占百分比与年污水占长江总水量百分比，、类水所占百分比与年污水占长江总水量百分比的倒数，分别进行拟合。选取拟合度最高的模型来预测未来十年长江、水所占百分比。针对第五问，提出了解决长江水质污染问题的从四方面着手的方案： 沿江工厂的整治，民众意识的唤醒，上游植被的保护，以及法律的硬性要求。三符号说明k：降解系数：第i个观测点的监测污染浓度。：第i个观测点污水流至第i+1个观测点时的污染浓度。：第i个观测点的监测水流量。：第i个观测点的排污量。：后一地段较之前一地段的污水排放量。：第i个观测点与第i+1个观测点之间的河长。：第i个观测点与第i+1个观测点水流速度

10、的平均值。：第i个观测点到第i+1个观测点水流所需的时间。g(x)：以1995年算起第x年，劣水占长江总水量水的百分比。y(x)：以1995年算起第x年，加水占长江总水量水的百分比。S：长江过水过水断面的面积。H：长江干流的平均河长。W：水流的流量。U：水流的流速。q:需要治理的水的体积。Hi：第i年长江干流的河长。x：以1995年算起第x年。y：年废水排放总量。四模型假设1.主要污染物高锰酸盐指数和氨氮的降解系数取0.2(单位：1/天)。2.一个观测站的水质污染主要来自本地区的排污和上游的污水，河水在通过观测点后立即被混入下段区域的污水。3.不考虑洪水干旱等特殊气候对水质的影响。4.两观测点

11、之间江水的流速是所测的两点流速的平均值。520052014年的评价河长与2004年的评价河长相等。五模型建立和求解问题一 根据03、04年长江流域水质报告表，利用EXCEL做出近两年来长江流域水质报告图 ，如下图：:时评价评价类类类类类劣类段范围河长河长百分比河长百分比河长百分比河长百分比河长百分比河长百分比枯全流域32679258239403.527.9513460.540.353613111091.53.35468714.35水干流38962.56801.7510809.527.7514947.538.35672917.32880.57.429557.55期支 流6283.51001.61

12、82929.15255240.61070.517.054376.952954.65丰全流域326795801.758980.527.512395.537.95658.517.352443.57.526608.1水干流1951013322.155469.2526.0758046.2540.925232313.65780.255.3752638.511.825期支 流36016.75574.51.5758980.527.72513961.2538.775574515.852526.2573251.59.125水全流域6283.51001.6182929.15255240.61070.517.054

13、376.952954.65文干流326795801.751371.527.512395.537.95658.517.352443.57.526608.1年支流1951013322.155469.2526.0758046.2540.925232313.65780.255.3752638.511.825图一 近两年来长江流域水质报告图由于在一年中各个不同的月份，排污水也有一定的时间限制。也就是说时间不同各个排污工厂工作量不同。主要取决于不同的季节市场对各个商品的需求不同。即于价值规律有一定的关系。因此我们单独的纵向比较各个月份的污染有一定的局限性，并不能很好的整体反映长江总体评价。所以我们必须进行

14、横向的比较和归纳。对每一个地区在近两年的28个月中的水质情况进行统计，再找出该地区污染的种类及该种类污染出现的频率。结果如下表所示：观测站pHDOCODMnNH3-N四川攀枝花龙洞001/281/28重庆朱沱0000湖北宜昌南津关0000湖南岳阳城陵矶0000江西九江河西水厂0000安徽安庆皖河口0000江苏南京林山0000四川乐山岷江大桥016/287/2810/28四川宜宾凉姜沟001/283/28四川泸州沱江二桥07/282/285/28湖北丹江口胡家岭0000湖南长沙新港01/28010/28湖南岳阳岳阳楼003/280湖北武汉宗关01/2800江西南昌滁槎1/2810/282/282

15、1/28江西九江蛤蟆石004/280江苏扬州三江营01/2800 表一 各地污染频率统计表从上表可以看出四川攀枝花龙洞，四川乐山岷江大桥，四川宜宾凉姜沟，四川泸州沱江二桥，江西南昌滁槎这五个地方都出现了CODMn和NH3-N超标情况。各地区水质的污染状况：1.四川乐山岷江大桥CODMn和NH3-N的平均超标率以及含氧量未达标率分别为：7/2，10/28和16/28，即该地区在抽样的28个月内有7个月出现高锰酸盐超标的情况，有10个月出现氨氮超标的情况,16个月出现含氧量未达标。 CODMn和NH3-N的污染严重，且含氧量不能达标，水质污染严重。2.四川攀枝花龙洞CODMn和NH3-N的平均超标

16、率以及含氧量未达标率分别为：1/28，1/28，0/28。 出现极少月CODMn和NH3-N污染，且含氧量都能达标，污染程度较轻。3.四川宜宾凉姜沟CODMn和NH3-N的平均超标率以及含氧量未达标率分别为：1/28，3/28，0/28。出现极少月CODMn和NH3-N污染，且含氧量都能达标，污染程度较轻。4.四川泸州沱江二桥CODMn和NH3-N的平均超标率以及含氧量未达标率分别为：2/28，5/28和7/28。出现少数几月CODMn污染。多月NH3-N和含氧量不能达标，水质污染严重。5. 江西南昌滁槎CODMn和NH3-N的平均超标率以及含氧量未达标率分别2/28，21/28和10/28。

17、 出现少数几月CODMn污染。极少月NH3-N能够达标，且含氧量不能达标，水质污染严重。6. 四川泸州沱江二桥CODMn和NH3-N的平均超标率以及含氧量未达标率分别为：0/28，10/28和1/28。出现少数几次CODMn污染，多月NH3-N未能够达标，极少月含氧量不能达标，水质污染严重。问题二查阅有关喝水自然净化能力的参考文献1，结合相关数据，建立水质依靠流量、流速和降解系数的数学模型，从而算出长江干流沿岸各个地段的排污量。建立模型表达式： 其中k=0.2/天=/s ,由于k的值很小，所以上式可以近似为：长江干流沿岸各个地段的排污量： 通过直接观察mi，很难看出主要的污染源在哪一段。因此我

18、们将与相减，得到后一地段较之前一地段的污水排放情况。若此值为正，则说明后一地段的污染比前一地段更为严重。若此值为负，则说明后一地段的污染比前一地段要好。由于需要计算的的数量非常多，我们借助C语言程序进行计算。源代码见附录。将各个地段在近一年里的列为表格，如下：表二 CODMn排出量与前一段排除量差值表（单位：kg/s）时间重庆朱沱湖北宜昌湖南岳阳江西九江安徽安庆江苏南京2004.0439813210038280-2881.41140.9-376192004.051282430580144001858016420-467002004.0632575345401960-2400-253401484

19、02004.07429361852028580-31020-1597035502004.08-4925320029200-29770-10370-37902004.09171132-2754027920-71500-19600675802004.123312447507380-309301798078602004.11121403800110609205280-21202004.126993.4924519690-7090-639058302005.014371.6345733086-15609582040202005.025934.62534.61651663945060-22602005.

20、038320.7187217986-312432040-190802005.046593.856608036920416260-13310平均28188.7814055.2819545.69-12248.21717.762-1630.69表三 NH3-N排出量与前一段排除量差值表（单位：kg/s）时间重庆朱沱湖北宜昌湖南岳阳江西九江安徽安庆江苏南京2004.042344.52772986931-1097-58942004.053276.62403620-41944-53952004.062395.643461008-34851032-34712004.072765.66805044-43562

21、698-54132004.08-193827583987-843-2360-35892004.097050.878996252060-7812-39082004.13465.26673635-3316558-37342004.112115.9338.11776-1108-1172-10602004.123264.4-2156.52133-1733390-4842005.012121-1165.42434.4-1150.81717082005.021889.8-1169.82619.6-2092.4112722152005.032397.6-1575.41972.6-1383.83349-449

22、22005.042235.3-1597.52338.8-926.89161390平均2568.0231092.1922244.569-1339.06-19.6923-2548.23由“ODMn排出量与前一段排除量差值表”可以得到结论：长江中游湖北宜昌至湖南岳阳段是高锰酸盐的主要污染源的所在地。由“NH3-N排出量与前一段排除量差值表”可以得到结论：长江中游湖北宜昌至湖南岳阳段是氨氮的主要污染源的所在地。问题三首先筛选取出水文年支流和干流的相关数据和各年长江总流量和废水排放量的数据。利用EXCEL作出1995年到2004年污水排放量变化趋势图如下： 图二 1995年到2004年污水排放量变化趋势

23、图由于给出的数据较少，因此每一个数据对模型的影响都很大。所以如果数据中出现奇点，就一定要去掉。用Matlab中的一元线性回归模型将上图中各个点（第7年的废水排放总量异常，将其去除）拟合成一条反映废水排放总量与年份之间关系的直线，方程为：该方程的拟合度R2= 0.9678 符合要求。 图三利用EXCEL作出1995年到2004年长江总流量变化趋势图如下：图四 1995年到2004年长江总流量变化趋势图利用EXCEL作出1995年到2004年长江废水比例变化趋势图如下： 图五 1995年到2004年长江废水比例变化趋势图从图中可以看出长江总流量中废水排放量的比例逐年上升，由此可知长江水呈现逐年恶化

24、的趋势（注：98年为特大洪水）。利用matlab拟合： 图六回归方程为：y=0.0001x*x+0.0002x +0.0184，其中y表示各年废水排放量浓度，x=1,2,10,分别对应2005年到2014年。此回归方程的拟合度R-square: 0.9795回归方程拟合度满足要求图七 未来十年长江废水比例变化趋势图从回归方程可以看出未来10年废水排放量站长江总量变化很快，说明未来的污染排放量需要进行控制，否则以这样的速度下去长江的污染情况将会无法控制。问题四1.建立求解的模型根据近十年的长江干流河长，我们可以求的水文年干流的平均河长：H= 5392.6km 假设水流量为W，水断面S，平均流速为

25、U。有公式W=SU S=W/U求得干线上7个站点的平均过水断面：测点攀枝花重庆朱沱湖北宜昌湖南岳阳江西九江安徽安庆江苏南京/766.895584.0616549.2516535.6420490.2522908.8621721.62再平均7个站点的过水断面我们可以得到长江的平均过水断面S=14936.27长江的总水量为：SH根据题意在未来十年为了保证长江干流中的四类和五类水的比例控制在20%以内，且没有劣五类水，则每年至少需要处理的污水量为：2.预测、类水，、类水，类水所占河长的百分比用Matlab中的regress函数对、类水所占百分比与年污水排放量，、类水所占百分比与万倍年污水浓度，、类水所

26、占百分比与年污水浓度的倒数，分别进行拟合。结果如下：、类水所占百分比（y 单位：%）比与年污水排放量（x 单位：亿吨）之间的关系： y= - 0.2252 x+129.4473 拟合度R2：0.921、类水所占百分比（y 单位：%）比与万倍年污水浓度（x 单位：吨/立方米）之间的关系： y=-0.2272 x+132.1354 拟合度R2：0.8904、类水所占百分比（y 单位：%）比与年污水浓度的倒数（x 单位：立方米/吨）之间的关系： y= 0.9746 x+35.6981 拟合度R2：0.8745比较上述三种方案的拟合度，选择第一种。即用、类水所占百分比（y 单位：%）比与年污水排放量（

27、x 单位：亿吨）之间的关系来预测未来十年、类水所占百分比。、类水所占百分比（y 单位：%）比与年污水排放量（x 单位：亿吨）之间的关系： y= 0.1347 x-12.2029 拟合度R2：0.7511、类水所占百分比（y 单位：%）比与万倍年污水浓度（x 单位：吨/立方米）之间的关系： y=0.1605 x-20.0099 拟合度R2：0.7893、类水所占百分比（y 单位：%）比与年污水浓度的倒数（x 单位：立方米/吨）之间的关系： y= -0.7285 x+50.0282 拟合度R2：0.8678比较上述三种方案的拟合度，选择第三种。即用、类水所占百分比（y 单位：%）比与年污水浓度的倒

28、数（x 单位：立方米/吨）之间的关系来预测未来十年、类水所占百分比。用100%减去、类水所占百分比再减去、类水所占百分比，即可得到劣水占长江总水量水的百分比。3.数据的计算及问题的求解结果根据拟合得到的两个公式：y(x)= - 0.2252 x+129.4473 g(x)= -0.7285 x+50.0282带入数据我们可以得到未来10中y（x）和g（x）占长江总水量的百分比：X11121314151617181920y(x)26.7829327.7829328.7004429.5452630.3257131.0488631.7208132.346832.931433.47858g(x)10.

29、1474712.1313714.1977616.3368418.5402920.8010423.1129925.470927.870230.30692根据题意在未来十年为了保证长江干流中的四类和五类水的比例控制在20%以内，且没有劣五类水，则每年至少需要处理的污水量为：利用求q公式，以及得到的y（x）和g（x）的数据，得到需要处理污水量为：11121314151617181920q（亿立方米）136.36646160.3963184.4297208.4630232.4963256.5296280.5629304.5962328.6296352.6629问题五对解决长江水质污染问题的建议：1、对

30、于长江沿线的钢铁厂、化工厂、造纸厂、造船厂、拆船厂这样的重污染企业必须禁止它们直接将污水排入长江。严格要求他们将污水通过处理设备处理后对各种污染物含量进行检测，检测所要排放的水属于可饮用水（即附录一中的前三类）即可进行排放。若发现工厂将检验未达标水排入长江则处以罚款或者是责令他们拆迁工厂。2、对于生活污水也需要进行一定的控制：唤醒民众、领导干部的环保意识。 由于长江两岸农田过量喷洒的农药被雨水冲刷到河里，导致氨氮含量过高、河水富营养化、蓝藻疯长的情况严重破坏水资源。故迫切需要提高农民的环保意识，尽量少使用农药同时也不要将农药水随便排放。这就需要政府拨款在长江沿岸多建水净化厂，将将要排放到长江里

31、的水进行净化处理后再排放。3、在长江上游地区植树造林，据统计长江上游地区森林覆盖率曾达到60%至80%，眼下却是5%至7%，因此每年流入长江泥沙24亿吨。由此，源头出现沙漠化及断流。长江污染还殃及近海。2003年出现赤潮近40次，赤潮导致氧亏，进而影响渔业资源及品质。水生生态系统濒于崩溃。因此应该鼓励在长江上游地区进行植树。 4、对保护长江立法，让那些不愿意保护长江的人受到法律的约束。主要分两步进行：一方面制定促进循环经济发展的政策和法律法规；另一方面，还需要将建议加大对违法排污行为的处罚力度，要罚就要罚到不法企业破产，把对长江的生态环境的合理开发纳入沿江城市政府官员的政绩考核体系，对严重破坏

32、生态环境并造成生态恶果的地区，应执行官员任用的一票否决制。六模型分析和检验误差分析 1.模型对于水质的净化只考虑了自然降解，而对于实际情况比如是否存在污水处理厂等认为净化设施没有考虑。对于降解系数选取，并认为全流域相等，与实际情况存在误差。 2.对于污水在监测站之后马上汇入江水的假设，与实际沿江不规律分布的污水排放不符。 3.污染物排放量与诸多因素有关2，仅仅估计其与时间的关系会造成偏差。 4.用给出的前十年的水质数据来预测后十年的水质变化情况，显得数据不够充足，容易造成误差。七模型评价模型改进 在原模型中，我们用已知的十个年份的水质情况来预测后十年的水质污染情况。这样越到后来误差越大。为了减

33、小误差，可以先将已知的十组数据进行拟合，来预测2005的水质污染情况。再将这十一组数据进行拟合，来预测2006的水质污染情况。以此类推，对未来十年的水质污染情况作出预测。八参考文献1李静.鄂西某市小容量河流水污染综合控制研究J.环境污染治理技术与设备，2003，(12).2王坚.影响河流水质的因素分析J.山西煤炭管理干部学院学报，2003,(3).3薛毅.数学建模基础M.北京：科学出版社，2011.附录用以计算CODMn的源代码：#include#includeint main()int k,h,x,y;float m,t,f=1/432000,p,q,l,v;float e=2.71828;

34、float c136;float d136;float g136;float a277=0 , 950, 1728, 2123, 2623, 2787, 3251,3690 , 13800, 21000, 25600, 28100, 29500, 29800,3.7 ,2.1,0.9,0.9,1.0,1.1,1.2,3720 , 13100, 19800, 20500, 29800, 34000, 34500,3.7 ,1.9,0.8,0.9,1.1,1.1,1.2,4010 , 14200, 20300, 22600, 29500, 32100, 33100,3.9 ,2.1,1.2,1.3

35、,1.5,1.5,1.6,4660 , 16400, 22700, 24100, 27000, 31900, 32100,4.1 ,2.3,1.4,1.5,1.5,1.6,1.7,3740 , 10600, 24000, 25900, 32100, 33400, 35100,3.8 ,2.1,1.4,1.4,1.5,1.7,1.7,6280 , 47600, 53500, 53800, 72800, 74200, 81000,5.1 ,4.8,1.7,1.9,2.1,3.4,3.4,3260 , 16200, 19100, 22300, 24800, 31000, 38400,3.1 ,2.3

36、,1.5,1.6,1.6,1.7,1.9,1500 , 8170, 10600, 12000, 14600, 17000, 19600,2.7 ,1.9,0.7,0.8,0.9,0.9,1.0,951 , 6550, 7400, 10700, 13200, 14100, 14900,3.1 ,1.5,0.7,0.8,0.8,0.8,0.9,712 , 4020, 4570, 8190, 10900, 12300, 14400,2.1 ,1.5,0.5,0.6,0.7,0.7,0.8,612 , 3603, 4510, 7980, 10300, 13700, 15100,2.0 ,1.0,0.4

37、,0.6,0.7,0.7,0.8,623 , 4740, 5180, 7040, 14300, 21400, 21500,1.9 , 0.9,0.4,0.6,0.8, 0.8,0.9,642 , 3650, 5400, 7240, 15100, 20200, 22100,2.1 ,1.2,0.4,0.5,0.7,0.8,0.8;float b718=2.3, 4.3, 2.5, 2.4, 5.8, 6.1, 0.8, 2.8, 1.6, 1.2, 0.9, 1.1, 1.1, 2.5, 2.9, 4, 2.6, 4, 3.5, 2.2, 3, 3.3, 2, 4.4, 1.6, 2, 1.3,

38、 1.3, 1.8, 1.9, 2, 1.3, 1.7, 2.7, 2.9, 2.9, 2.2, 3, 3.8, 3.2, 3.1, 3.4, 3.7, 1.9, 2.4, 1.9, 2, 2.1, 2.4, 2.4, 2.1, 2.6, 1.3, 3.2, 3.3, 3.6, 3.5, 4.2, 4, 3.9, 3.5, 2.6, 3.5, 5.1, 3.2, 4.1, 2.9, 3, 4.7, 3.1, 4.2, 4.9, 2.9, 3.1, 2.6, 2.6, 2.3, 1.9, 1.9, 2.2, 2.3, 2.4, 3.1, 1.8, 2, 2, 2.8, 3.1, 2.3, 2.3

39、, 2.8, 3.2, 1.6, 1.7, 1.9, 1.6, 2.1, 2.2, 1.7, 2.6, 2.7, 2.7, 2.3, 2.9, 2.6, 3.2, 2.7, 2.7, 1.5, 1.8, 2, 1.8, 1.7, 2.3, 1.9, 1.8, 2, 2.5, 2.3, 1.8, 1.5, 1.4, 1.4, 3.8, 3.3, 3.5;for(k=1;k=13;k+) for(h=0;h=5;h+)p=a0h+1-a0h;l=a2*kh+1+a2*kh;q=l/2000;t=p/q; ck-1h=pow(bhk-1*pow(e,-k),t)*a2*k-1h+1-bhk-1*a2

40、*k-1h; dk-1h=bh+1k-1*a2*k-1h+1-pow(bhk-1*pow(e,-k),t)*a2*k-1h+1;for(x=0;x=12;x+) for(y=0;y=5;y+) gxy=cxy+dxy; printf(%.1f ,gxy); printf(n);return 0;用以计算NH3-N的源代码：#include#includeint main()int k,h,x,y;float m,t,f=1/432000,p,q,l,v;float e=2.71828;float c136;float d136;float g136;float a277=0 , 950, 17

41、28, 2123, 2623, 2787, 3251,3690 , 13800, 21000, 25600, 28100, 29500, 29800,3.7 ,2.1,0.9,0.9,1.0,1.1,1.2,3720 , 13100, 19800, 20500, 29800, 34000, 34500,3.7 ,1.9,0.8,0.9,1.1,1.1,1.2,4010 , 14200, 20300, 22600, 29500, 32100, 33100,3.9 ,2.1,1.2,1.3,1.5,1.5,1.6,4660 , 16400, 22700, 24100, 27000, 31900, 32100,4.1 ,2.3,1.4,1.5,1.5,1.6,1.7,3740 , 10600, 24000, 25900, 32100, 33400, 35100,3