环境系统分析第四次作业.docVIP

环境系统分析第四次作业 环境系统分析第四次作业 PAGE 环境系统分析第四次作业 环境系统分析 上机实验报告 上机主题 环境质量基本模型的求解与分析 上机时间 班 级 学 号 姓 名 评 分 　　 市政环境学院 环境科学与工程系 指导教师 一．上机实验目的 1.学习并了解环境质量基本模型的相关知识。 2.学会并使用matlab求解微分方程。 3.学会建立简单的环境质量基本模型并通过已知量求解相关参数值。 二．上机实验题目 (1)一维均匀稳态河流，初始断面的污染物C0=50mg/l，纵向弥散系数Dx=s,衰减系数K1=d，河流断面平均流速为s。试求下述各条件下在下游500m处的污染物浓度。（1）一般解析解；（2）忽略弥散作用的解;(3)忽略推流作用时的解;(4)忽略衰减作用时的解。 (2)均匀稳态河段的宽500m，平均水深h=3m，平均流速ux=1m/s，横向弥散系数Dy=1m2/s，守恒污染物中心排放的源强为1000kg/h。求排放口下游2km 三．上机程序要点 1.（1）建立M文件，由题知C0=50mg/l，Dx=s，K1=d，ux=s,x=500m 显然该题在稳态条件下， 代入一维模型稳态排放的解析解解得 （2）忽略弥散作用： 解析解为 所以代入数值得： （3）忽略推流作用： 其解析解为 代入数值得： （4）忽略衰减作用： 求解 得，故其浓度 2.由题知、、，、 （1）可得 因为污染物为中心排放，故污染物扩散羽宽度 （2）显然在河中心的污染物浓度最大 （3）在忽略横向推流和纵向弥散的条件下下： 所以 因为该题宽度有限 故 取,用matlab进行画图 可得： 四．上机实验心得 求解微分方程可用到dsolve（）函数，之后可以用subs（）函数进行代入数值，但是偶尔会出现一长串没有化简的式子，目前还不明白是什么情况造成的。 求解时应加上初值，不然解出来的值会带有参数。