数学实验结课报告《数学实验》实验报告题目:黄河小浪底调水调沙工程姓名:胡迪学号:201014622专业:信息与计算科学黄河小浪底调水调沙问题2004年6月至7月黄河进行了第三次调水调沙试验,特别是首次由小浪底、三门峡和万家寨三大水库联合调度,采用接力式防洪预泄放水,形成人造洪峰进行调沙试验获得成功。整个试验期为20多天,小浪底从6月19日开始预泄放水,至到7月13日恢复正常供水结束。小浪底水利工程按设计拦沙量为75.5亿m3,在这之前,小浪底共积泥沙达14.15亿t。这次调水调沙试验一个重要的目的就是由小浪底上游的三门峡和万家寨水库泄洪,在小浪底形成人造洪峰,冲刷小浪底库区沉积的泥沙,在小浪底水库开闸泄洪以后,从6月27日开始三门峡水库和万家寨水库陆续开闸放水,人造洪峰于29日先后到达小浪底,7月3日达到最大流量2700,使小浪底水库的排沙量也不断地增加。表1是由小浪底观测站从6月29日到7月10日检测到的试验数据。表1试验观测数据(单位:水流为sm/3,含沙量为3kg/m)6.296.307.17.27.37.48:0020:008:0020:008:0020:008:0020:008:0020:008:00218001900210022002300240025002600265027002720232607585909810010210811211517.57.67.77.87.97.108:20:008:20:008:20:008:20:008:20:008:20000000000002600250023002200200018501820180017501500100091181201181058060503026208注:以上数据主要是根据媒体公开报道的结果整理而成。现在,根据试验数据建立数学模型研究下面的问题:(1)给出估算任意时刻的排沙量及总排沙量的方法;(2)确定排沙量与水流量的变化关系。关键词:拟合,SAS,Matlab,线性回归,调水调沙实验问题分析:1、对于问题一,所给数据中水流量x和含沙量h的乘积即为该时刻的排沙量y即:y=hx。2、对于问题二,研究排沙量与排水量的关系,从实验数据中可以看出,开始排沙量随水量增加而增加,而后随水流量的增加而减少,显然变化关系并非线性的关系,为此,把问题分为两部分,从水流量增加到最大值为第一阶段,从水流量最大值到结束为第二阶段,分别来研究水流量与排沙量之间的函数关系。模型假设:1、水流量和排沙量都是连续的,不考虑上游泄洪所带来的含沙量和外界带来的含沙量。2、时间是连续变化的,所取时间点依次为1,2,3,⋯,24,单位时间为12h。模型的建立与求解:<一>对于问题一,因为排沙量与时间的散点图基本符合正态曲线,如图二所示。所以,排沙量的对数与时间的函数关系就应该符合二次函数关系,因而排沙量取对数后,再与时间t进行二次回归,排沙量取自然后的数据见表2.假设排沙量与时间函数关系的数学模型是ecbtaty2^两边取对数得Lny=at^2+bt+c先由表二做出排沙量的自然对数lny与时间t的散点图见图一,并利用SAS软件进行拟合,得到排沙量的自然对数与时间的回归方程为:Lny=-0.0209t^2+0.4298t+10.6321由回归拟合参数表可知回归方程是显著的,因为相关系数人R^2=0.9629,误差均方S^2=0.0543,说明回归曲线拟合效果很好。所以排沙量与时间之间的函数关系式为etty6312.104289.02^0209.0图二:排沙量对时间的曲线图时间点31xmsg3hkgmg1ykgsglny11800325760010.96128219006011400011.64395321007515750011.96718422008518700012.13886523009020700012.24047624009823520012.368197250010025000012.429228260010226520012.488249265010828620012.5644510270011230240012.6195111272011531240012.6533212265011630740012.6359113260011830680012.6339514250012030000012.6115415230011827140012.5113516220010523100012.350171720008016000011.982931818506011100011.61729191820509100011.41864201800305400010.89674211750264550010.72547221500203000010.30895231000880008.98719724900545008.411833最后dtett2406312.104289.02^0209.0*60*60*12对所求出的函数关系在区间[0,24]之间进行积分结果为总排沙量1.93962亿吨,此与媒体报道的排沙量几乎一样。<二>对于第二个问题,两个阶段的数据如表三、表四所示表三:第一阶段试验数据序号12345678水流量x18001900210022002300240025002600含沙326075859098100102量h表四:第二阶段...