黄河小浪底

问题的提出•2004年6月至七月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功。整个试验期为20多天，小浪底从6月19日开始预泄放水，直到7月13日恢复正常供水结束。小浪底水利工程设计拦沙量为75.5亿立方米，在这之前，•小浪底共积泥沙达14.15亿t。这次调水调沙试验一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙，在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2700立方米每秒，使小浪底水库的排沙量也不断增加，表1是由小浪底观测站从6月29日到7月10日检测到的实验数据。日期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•对于问题（2），研究排沙量与水流量的关系，从实验数据可以看出，开始排沙量是随着水流量的增加而增长，而后是随着水流量的减少而减少。显然，变化规律不是曾线性的关系，因此，把问题分为两部分，从开始水流量增加到最大值为第一阶段，从水流量的最大值到结束为第二阶段，分别来研究水流量与排沙量的关系。具体数据如下表。序号1234567891011水流量18001900210022002300240025002600265027002720含沙量326075859098100102108112115第一阶段实验观测数据单位：水流量为m^3/s,含沙量为kg/m^3第二阶段实验观测数据单位：水流量为m^3/s,含沙量为kg/m^3序号12345678910111213水流量265026002500230022002000185018201800175015001000900含沙量11611812011810580605040322085•对于第一阶段，由表5-3用Matlab作图可以看出其变化趋势，我们用多项式作最小拟合。•程序如下：•x=[18001900210022002300240025002600265027002720];•y=[326075859098100102108112115];•plot(x,y,'r')•如图所示：•利用已知数据对其作三次多项式拟合，编写MATLAB命令如下：•x=[18001900210022002300240025002600265027002720];y=[326075859098100102108112115];•A=polyfit(x,y,3)•A=•1.842*e-7•-1.317*e-3•3.1784•-2.492.9•z=polyval(A,x);•plot(x,y,'*',x,y,'r',x,z,'b')•可得：a0=-2.492.9a1=3.1784a2=-1.317*e-3a3=1.842*e-7•于是可得拟合多项式为：xxeexx323*)7*842.1(*)3*317.1(*1784.39.2492)(410109312.1361094.123102626.70624.126557.7434xxxxy三次多项式拟合效果图•利用已知数据对其作四次多项式拟合，编写MATLAB命令如下：•x=[18001900210022002300240025002600265027002720];•y=[326075859098100102108112115];•A=polyfit(x,y,4)•A=•-1.93*e-10•1.94*e-6•-7.26*e-3•0.01206*e+3•-7.4347*e+3•z=polyval(A,x);•plot(x,y,'*',x,y,'r',x,z,'b')•及得到•A0=-7.4347*e+3a1=0.01206*e+3a2=-7.26*e-3a3=1.94*e-6a4=-1.93*e-10•于是可得拟合多项式为：410109312.1361094.123102626.70624.126557.7434xxxxy四次多项式拟合效果图•从上面的三次多项式拟合和四次多项式拟合效果来看，差别不太大，基本可以看出排沙量与水流量的关系。•对于第二阶段，由表5-4可以类似处理。用线性最小二乘法作三次和四次多项式拟合，拟合效果与MATLAB程序如下：•第二阶段三次多项式拟合效果图：•x=[265026002500230022002000185018201800175015001000900];•y=[11611812011810580605040322085];•A=polyfit(x,y,3)•A=•-0.00000.0006-0.9475464.9601•z=polyval(A,x);•plot(x,y,'*',x,y,'r',x,z,'b')第二阶段三次多项式拟合效果图•四次多项式拟合效果图：•x=[265026002500230022002000185018201800175015001000900];•y=[11611812011810580605040322085];•A=polyfit(x,y,4)•A=•-0.00000.0000-0.00131.1219-354.5952•z=polyval(A,x);•plot(x,y,'*',x,y,'r',x,z,'b')第二阶段四次多项式拟合效果图