数值分析第五章思考题.pdf

《数值分析第五章思考题.pdf》由会员分享，可在线阅读，更多相关《数值分析第五章思考题.pdf（8页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、数值分析第五章思考题 黄河小浪底调水调沙问题 一、问题的提出 2004 年 6 月至 7 月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄防水，形成人造洪峰进行调沙试验获得成功。整个试验期为 20 多天，小浪底从 6 月 19 日开始预泄放水，直到 7 月 13 日恢复正常供水结束。小浪底水利工程按设计拦沙量为亿3m，在这之前，小浪底共积泥沙达亿 t。这次调水调沙试验一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙，在小浪底水库开闸泄洪以后，从 6 月 27 日开始三门峡水库和万家寨水库陆

2、续开闸放水，人造洪峰于 29 日先后到达小浪底，7 月 3 日达到最大流量27003m/s，使小浪底水库的排沙量也不断增加。表-1 是由小浪底观测站从 6 月 29日到 7 月 10 日检测到的试验数据。表-1 试验观测数据 单位：水流为3m/s 日期 时间 8：00 20：00 8：00 20：00 8：00 20：00 8：00 20：00 8：00 20：00 8：00 20：00 水流量 1800 1900 2100 2200 2300 2400 2500 2600 2650 2700 2720 2650 含沙量 32 60 75 85 90 98 100 102 108 112 11

3、5 116 日期 时间 8：00 20：00 8：00 20：00 8：00 20：00 8：00 20：00 8：00 20：00 8：00 20：00 水流量 2600 2500 2300 2200 2000 1850 1820 1800 1750 1500 1000 900 含沙量 118 120 118 105 80 60 50 30 26 20 8 5 现在，根据试验数据建立数学模型研究下面的问题：（1）给出估算任意时刻的排沙量及总排沙量的方法。（2）确定排沙量与水流量的变化关系。二、具体求解过程：首先做一些假设：所给数据客观准确的反应了现实情况；数据是连续的；不考虑外在因素；时间化

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 工作窗口输入程序为：S=32 60 75 85 90 98 100 102 108 112 115 116 118 120 118 105 80 60 50 3

5、0 26 20 8 5;W=1800 1900 2100 2200 2300 2400 2500 2600 2650 2700 2720 2650 2600 2500 2300 2200 2000 1850 1820 1800 1750 1500 1000 900;T=1:24;subplot(2,1,1);plot(T,S);hold on;plot(T,S,.);xlabel(时间 t/12h);ylabel(含沙量/公斤每立方米);title(时间与含沙量关系);subplot(2,1,2);plot(T,W);hold on;plot(T,W,.);xlabel(时间 t/12h);y

6、label(水流量/立方米每秒);title(时间与水流量关系);运行后屏幕上显示为：通过观察图像，我们可以看出其变化并不光滑，而且也没有特定的表现出服从某种分布的趋势。为了得到具体的计算函数，我们就必须对数据进行拟合，所以通过 Matlab 先利用 spline 方法对数据进行插值，从而提高精确度，使图像变得光滑，然后利用多项式进行拟合。1.时间与含沙量的曲线拟合，先画出散点图，在MATLAB 工作窗口输入：S=32 60 75 85 90 98 100 102 108 112 115 116 118 120 118 105 80 60 50 30 26 20 8 5;T=1:24;x=1:

7、24;plot(T,S,r*);legend(数据点 S)xlabel(时间 t/12h);ylabel(含沙量/公斤每立方米);title(时间与含沙量关系散点图);为了提高精度用三次多项式拟合，在 MATLAB 工作窗口输入程序为：S=32 60 75 85 90 98 100 102 108 112 115 116 118 120 118 105 80 60 50 30 26 20 8 5;a=polyfit(x,y,3)运行后输出为：a=21 16 故拟合多项式为：230.014x-1.3x()2116f xx 即可得出时间与含沙量的关系式230.014t-1.3t()2116s tt

8、 在 MATLAB 工作窗口输入：S=*t.*t.2+21*t+16;t=1:24;plot(t,S);title(时间与含沙量关系图)运行后屏幕显示为：2.时间与水流量的曲线拟合，先画出散点图，在 MATLAB 工作窗口输入：W=1800 1900 2100 2200 2300 2400 2500 2600 2650 2700 2720 2650 2600 2500 2300 2200 2000 1850 1820 1800 1750 1500 1000 900;T=1:24;x=1:24;plot(T,W,r*);legend(数据点 S)xlabel(时间 t/12h);ylabel(水

9、流量/立方米每秒);title(时间与水流量关系散点图);为了提高精度用三次多项式拟合，在 MATLAB 工作窗口输入程序为：W=1800 1900 2100 2200 2300 2400 2500 2600 2650 2700 2720 2650 2600 2500 2300 2200 2000 1850 1820 1800 1750 1500 1000 900;a=polyfit(x,y,3)运行后为：a=+003*故拟合多项式为：230.13-14()2.40021.5003xf xxexe 即可得出时间与水量的关系式：230.13t-14t()2.40021.5003V tete 在

10、MATLAB 工作窗口输入：V=*t.3-14*t.2+002*t+003;t=1:24;plot(t,V);title(时间与水流量关系图)运行后屏幕显示为：由于 V=v(t)S(t)可得出排沙量与时间的关系：654320.0018*0.365*24.29*582.92*2866*35340*24000Vtttttt在 MATLAB 工作窗口输入：V=*t.*t.5+*t.*t.3+2866*t.2+35340*t+24000;t=1:24;plot(t,V);title(时间与排沙量关系图)即可得出时间与排沙量关系图：由于次数太高，在重新拟合一下，使之成为三次多项式,在 MATLAB 工作

11、窗口输入：（2）确定排沙量与水流量的变化关系。先将排沙量和水流量的相关数据反映到图像中，在 MATLAB 工作窗口输入程序为：t=1:24;v=*t.3-14*t.2+002*t+003;V=95*t.+003*t.2+004*+004;plot(v,V,.);title(整理图)figure;通过观察可以看出，其关系是分段的，大约在 t=9 处分段，所以我们按时间进行分段拟合。前半段求解，在 MATLAB 工作窗口输入程序为：t=1:9;v=*t.3-14*t.2+002*t+003;V=95*t.+003*t.2+004*+004;plot(v,V,.);title(前半段图)运行后屏幕显示为：327.5005*0.43*5.2002*3.6004yexxexe 后半段求解，在 MATLAB 工作窗口输入程序为：t=10:24;v=*t.3-14*t.2+002*t+003;V=95*t.+003*t.2+004*+004;plot(v,V);title(后半段图)运行后屏幕显示为：322.3005*0.066*1.9002*1.9005yexxexe 综上，我们可以得到排沙量与水流量的关系式为 32327.55*0.43*5.22*3.64092.35*0.066*1.92*1.95=924evvevetevveveVt