篇一:《黄河小浪底调水调沙问题》
黄河小浪底调水调沙问题
摘要：本文利用插值拟合的方法通过Matlab工具模拟出了排沙量与时间、排沙量与水
流量的函数关系，并且求出了总排沙量为1.704亿吨。整个模型简单且方便计算，其中排沙量与水流量的函数关系为分段函数。
关键词： 调水调沙 Matlab 插值拟合
一、问题重述
2004年6月至7月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功．整个试验期为20多天，小浪底从6月19日开始预泄放水，直到7月13日恢复正常供水结束．小浪底水利工程按设计拦沙量为75.5亿立方米，在这之前，小浪底共积泥沙达14.15亿吨．这次调水调试验一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙．在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2700立方米/每秒，使小浪底水库的排沙量也不断地增加．下面是由小浪底观测站从6月29日到7月10日检测到的试验数据：
表1: 试验观测数据 单位：水流为立方米 / 秒，含沙量为公斤 / 立方米
(1) 给出估算任意时刻的排沙量及总排沙量的方法； (2) 确定排沙量与水流量的变化关系。
二、模型假设
1. 假设所给数据客观准确的反应了现实情况 2. 假设所给数据遵循一定规律变化，即是连续的 3. 假设模型中不需要考虑一些外在因素 4. 假设可将时间化为等分的时间点进行计算
三、符号说明
t: 时间或时间点 v: 水流量 S: 含沙量 V: 排沙量
四、问题分析
假设水流量和含沙量都是连续的，那么某一时刻的排沙量V=v(t)S(t)，其中v(t)为t时刻的水流量，而S(t)为t时刻的含沙量。通过观察数据，这些数据是每个12小时采集一次，所以我们可以将时间设为时间点t，依次为1，2，3，„„，24，单位时间为12h。为了找到排沙量与时间的关系，我们就要先找到水流量和含沙量与时间的关系，一但找到水流量和含沙量与时间的关系，那么所要求的问题也就不难解决了。
五、模型的建立与求解
通过分析，我们假设水流量和含沙量都是连续的，那么我们开始对问题“ (1) 给出估算任意时刻的排沙量及总排沙量的方法”进行求解。
我们通过Matlab工具将所知道的数据显示为直观的图像，如下所示，具体程序见附录的%tuxing.m。
通过观察图像，我们可以看出其变化并不光滑，而且也没有特定的表现出服从某种分布的趋势。
但是为了得到具体的计算函数，我们就必须对数据进行拟合，所以通过Matlab先利用spline方法对数据进行插值，从而提高精确度，使图像变得光滑，然后利用多项式进行拟合，当多项式次数越高拟合也越准确，但是由于数据受到的影响较多，所以这里的数据也不是准确值，因此我们可以只取三次进行拟合，也方便了后续的计算。
于是我们分别对含沙量和水流量进行插值拟合，便可以得到下面图像和结果，具体程序见附录%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+003
即水流量与时间的关系式为：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的指数形式进行拟合，这里就不再重复计算。
我们得到了拟合函数，下面就可以计算出这几天的总排沙
量
，通过
Matlab编程可以计算出定积分，结果如下，程序详见附录%jisuan.m。
ans =170366976000
即总含沙量为1.704亿吨。
下面我们对问题“(2) 确定排沙量与水流量的变化关系。”进行分析计算。以下所有相关程序见附录%paishui.m，下面就不重复说明。
我们先利用Matlab将排沙量和水流量的相关数据反映到图像中。
篇二:《黄河小浪底调水调沙问题》
黄河小浪底调水调沙问题
摘要：本文利用插值拟合的方法通过Matlab工具模拟出了排沙量与时间、排沙量与水
流量的函数关系，并且求出了总排沙量为1.704亿吨。整个模型简单且方便计算，其中排沙量与水流量的函数关系为分段函数。
关键词： 调水调沙 Matlab 插值拟合
一、问题重述
2004年6月至7月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功．整个试验期为20多天，小浪底从6月19日开始预泄放水，直到7月13日恢复正常供水结束．小浪底水利工程按设计拦沙量为75.5亿立方米，在这之前，小浪底共积泥沙达14.15亿吨．这次调水调试验一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙．在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2700立方米/每秒，使小浪底水库的排沙量也不断地增加．下面是由小浪底观测站从6月29日到7月10日检测到的试验数据：
表1: 试验观测数据 单位：水流为立方米 / 秒，含沙量为公斤 / 立方米
(1) 给出估算任意时刻的排沙量及总排沙量的方法； (2) 确定排沙量与水流量的变化关系。
二、模型假设
1. 假设所给数据客观准确的反应了现实情况 2. 假设所给数据遵循一定规律变化，即是连续的 3. 假设模型中不需要考虑一些外在因素 4. 假设可将时间化为等分的时间点进行计算
三、符号说明
t: 时间或时间点 v: 水流量 S: 含沙量 V: 排沙量
四、问题分析
假设水流量和含沙量都是连续的，那么某一时刻的排沙量V=v(t)S(t)，其中v(t)为t时刻的水流量，而S(t)为t时刻的含沙量。通过观察数据，这些数据是每个12小时采集一次，所以我们可以将时间设为时间点t，依次为1，2，3，„„，24，单位时间为12h。为了找到排沙量与时间的关系，我们就要先找到水流量和含沙量与时间的关系，一但找到水流量和含沙量与时间的关系，那么所要求的问题也就不难解决了。
五、模型的建立与求解
通过分析，我们假设水流量和含沙量都是连续的，那么我们开始对问题“ (1) 给出估算任意时刻的排沙量及总排沙量的方法”进行求解。
我们通过Matlab工具将所知道的数据显示为直观的图像，如下所示，具体程序见附录的%tuxing.m。
通过观察图像，我们可以看出其变化并不光滑，而且也没有特定的表现出服从某种分布的趋势。
但是为了得到具体的计算函数，我们就必须对数据进行拟合，所以通过Matlab先利用spline方法对数据进行插值，从而提高精确度，使图像变得光滑，然后利用多项式进行拟合，当多项式次数越高拟合也越准确，但是由于数据受到的影响较多，所以这里的数据也不是准确值，因此我们可以只取三次进行拟合，也方便了后续的计算。
于是我们分别对含沙量和水流量进行插值拟合，便可以得到下面图像和结果，具体程序见附录%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+003
即水流量与时间的关系式为：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的指数形式进行拟合，这里就不再重复计算。
我们得到了拟合函数，下面就可以计算出这几天的总排沙{小浪底调沙}.
量
，通过
Matlab编程可以计算出定积分，结果如下，程序详见附录%jisuan.m。
ans =170366976000
即总含沙量为1.704亿吨。
下面我们对问题“(2) 确定排沙量与水流量的变化关系。”进行分析计算。以下所有相关程序见附录%paishui.m，下面就不重复说明。
我们先利用Matlab将排沙量和水流量的相关数据反映到图像中。
篇三:《黄河小浪底调水调沙》
黄河小浪底调水调沙问题
2004年6月至7月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功．整个试验期为20多天，小浪底从6月19日开始预泄放水，直到7月13日恢复正常供水结束．小浪底水利工程按设计拦沙量为75.5亿立方米，在这之前，小浪底共积泥沙达14.15亿吨．这次调水调试验一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙．在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2700立方米/每秒，使小浪底水库的排沙量也不断地增加．下面是由小浪底观测站从6月29日到7月10日检测到的试验数据：
表1: 试验观测数据 单位：水流为立方米 / 秒，含沙量为公斤 / 立方米
(1) 给出估算任意时刻的排沙量及总排沙量的方法； (2) 确定排沙量与水流量的变化关系。
篇四:《黄河小浪底调水调沙》
值分论文题目：指导老师：学 院：专 业：姓 名：学 号：
析课程论文
数
【实验课题】 黄河小浪底调水调沙问题
【实验目标】
（1）加深对插值及数据拟合知识的理解； （2）学会利用拟合实现计算有关数值方法； （3）验证插值拟合所预言的数值现象； （4）改进曲线拟合既有算法；
（5）掌握最小二乘法的基本原理，并会通过计算机解决实际问题。
【理论概述与算法描述】
为了确定排沙量与时间，排沙量与水流量的函数关系，我们需要对数据进行曲线拟合，所以通过Matlab对数据进行插值拟合，提高精确度，使图像变得光滑，然后利用多项式进行拟合。当多项式次数越高拟合也越准确，但是数据受到的影响较多，所以这里的数据也不是准确值，因此我们只取三次进行拟合，也方便了后续的计算。
符号说明
t: 时间或时间点 v: 水流量 S: 含沙量 V: 排沙量
【实验问题】
在小浪底水库蓄水后，黄河水利委员会进行了多次试验，特别是2004年6月至7月进行的黄河第3次调水调沙试验具有典型意义．这次试验首次由小浪底、三门峡和万家寨三大水库联合调度，进行接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功．这次调水调沙试验的一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙．在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2720 m3／s，使小浪底水库的排沙量也不断地增加．表1是由小浪底观测站从6月29日到7月10 日检测到的试验数据
表1: 试验观测数据 单位：水流为立方米 / 秒，含沙量为公斤 / 立方米
试根据实验数据建立模型解决下面问题(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的指数形式进行拟合。
得到了拟合函数，下面就可以计算出这几天的总排沙量错误！未找到引用源。，通过Matlab编程可以计算出定积分，结果如下，程序详见附录%jisuan.m。
ans =170366976000
即总含沙量为1.704亿吨。
(2) 确定排沙量与水流量的变化关系。
先利用Matlab将排沙量和水流量的相关数据反映到图像中。具体程序 见附录%paishui.m。
篇五:《3黄河小浪底调水调沙 场面震撼引人关注》
3黄河小浪底调水调沙 场面震撼引人关注
“黄龙”从黄河小浪底水利枢纽工程出水口飞流直下。
篇六:《黄河小浪底调水调沙问题》