黄河小浪底调水调沙.doc

数值分析课程论文论文题目：指导老师：学院：专业：姓名：学号：【实验课题】黄河小浪底调水调沙问题【实验目标】（1）加深对插值及数据拟合知识的理解；（2）学会利用拟合实现计算有关数值方法；（3）验证插值拟合所预言的数值现象；（4）改进曲线拟合既有算法；（5）掌握最小二乘法的基本原理，并会通过计算机解决实际问题。【理论概述与算法描述】为了确定排沙量与时间，排沙量与水流量的函数关系，我们需要对数据进行曲线拟合，所以通过Matlab对数据进行插值拟合，提高精确度，使图像变得光滑，然后利用多项式进行拟合。当多项式次数越高拟合也越准确，但是数据受到的影响较多，所以这里的数据也不是准确值，因此我们只取三次进行拟合，也方便了后续的计算。符号说明t:时间或时间点v:水流量S:含沙量V:排沙量【实验问题】在小浪底水库蓄水后，黄河水利委员会进行了多次试验，特别是2004年6月至7月进行的黄河第3次调水调沙试验具有典型意义．这次试验首次由小浪底、三门峡和万家寨三大水库联合调度，进行接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功．这次调水调沙试验的一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙．在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2720m3／s，使小浪底水库的排沙量也不断地增加．表1是由小浪底观测站从6月29日到7月10日检测到的试验数据表1:试验观测数据单位：水流为立方米/秒，含沙量为公斤/立方米日期6.296.307.17.27.37.4时间8:0020:008:0020:008:0020:008:0020:008:0020:008:0020:00水流量180019002100220023002400250026002650270027202650含沙量326075859098100102108112115116日期7.57.67.77.87.97.10时间8:0020:008:0020:008:0020:008:0020:008:0020:008:0020:00水流量26002500230022002000185018201800175015001000900含沙量11812011810580605030262085试根据实验数据建立模型解决下面问题：(1)给出估算任意时刻的排沙量及总排沙量的方法；(2)确定排沙量与水流量的变化关系。【实验过程与结果】(1)给出估算任意时刻的排沙量及总排沙量的方法通过Matlab工具将所知道的数据显示为直观的图像，如下所示，具体程序见附录的%tuxing.m。通过观察图像，看出其变化并不光滑，而且也没有特定的表现出服从某种分布的趋势。于是分别对含沙量和水流量进行插值拟合，便可以得到下面图像和结果，具体程序见附录%hansha.m和%liuliang.m所得到的拟合函数为：y=0.014*x^{3}-1.3*x^{2}+21*x+16即含沙量与时间的关系式为：S=0.014*t^3-1.3*t^2+21*t+16所得到的拟合函数为：y=0.13*x^{3}-14*x^{2}+2.4e+002*x+1.5e+00即水流量与时间的关系式为：v=0.13*t^3-14*t^2+2.4e+002*t+1.5e+003因为某一时刻的排沙量V=v(t)S(t)，所以将所拟合出来的多项式带入上式，通过Matlab进行计算可以得到下面答案，程序见附录%jisuan.m。ans=91/50000*t^6-73/200*t^5+2429/100*t^4-14573/25*t^3+2866*t^2+35340*t+24000即排沙量与时间的关系为：V=0.0018*t^6-0.365*t^5+24.29*t^4-582.92*t^3+2866*t^2+35340*t+24000由于这里的多项式次数过高，所以对其再进行一次拟合，结果如下，程序见附录%paisha.m。拟合后的函数为V=95*t^3-5.5e+003*t^2+7.7e+004*t-3.2e+004，通过图像可以看出排沙量与时间服从正态分布，所以化成形式e的指数形式进行拟合。得到了拟合函数，下面就可以计算出这几天的总排沙量QUOTE，通过Matlab编程可以计算出定积分，结果如下，程序详见附录%jisuan.m。ans=170366976000即总含沙量为1.704亿吨。(2)确定排沙量与水流量的变化关系。先利用Matlab将排沙量和水流量的相关数据反映到图像中。具体程序见附录%paishui.m。通过观察可以看出，其关系是分段的，所以按时间进行分段拟合，拟合原理同问题（1）相同，于是可以得到分段前后的拟合多项y=-7