水箱水流量问题

1、1水箱水流量问题班级：计算 122姓名：李克文 学号：120901055水箱的水流量问题（一）问题的提出：许多供水单位由于没有测量流入或流出水箱流量的设备， 而只能测量水箱中 的水位。试着通过测得的某时刻水箱中水位的数据，估计在任意时刻（包括水泵 灌水时间）t 流出水箱的流量 f（t）。假设：（1）影响水箱流量的唯一因素是该区域公众对水的普通需要。（2）水泵的灌水速度为常数。（3）从水箱中流出水的最大流速小于水泵的灌水速度。（4）每天的用水量分布都是相似的。（5）水箱的流水速度可用光滑曲线来近似。（6）当水箱的水容量达到 514.8g 时， 开始泵水； 达到 677.6 时， 便停止泵水。给出

2、下面原始数据表，其中长度单位为 E（ 1E=30.24cm）。水箱为圆柱体，其2直径为 57E。L 时间/s2水位10E时间/s _J2水位10E03175466363350331631104995332606635305453936316710619299457254308713937294760574301217921289264554292721240285068535284225223279571854276728543275275021269732284269779254泵水35932泵水82649泵水39332泵水8596834753943535508995333974331834

3、45932703340（二）关键字水箱，水流量，水箱容积，时间。（三）问题分析与建立模型 引入如下记号：V水的容积；Vi时刻 ti（ h）水的容积（单位 G，1G=3.785L （升）；f（ t）时刻 ti流出水箱的水的流速，它是时间的函数（G/h）； p水泵的灌水速度（G/h ）。根据要求先将上表中的数据做变换，时间单位用小时（h），水位高转换成水的体积（V=nR2h），得下表。时间（h）水量（103G）时间（h）水量（103G）30606.112.954639.50.921593.713.876622.41.843583.014.982604.62.950571.615.904589.33

4、.871562.616.826575.04.978552.117.932558.85.900544.119.038542.67.001533.649.959528.27.929525.4一 20.839514.88.968514.822.015/9.981/22.958/10.926/23.880663.410.954677.624.987648.512.033657.725.908637.6注：第一段泵水的始停时间及水量为t始=8.968 (h), v始=514.8% 103( G)t末=10.926(h),v末=677.6x103(G)第二段泵水的始停时间及水量为t始=20.839(h),

5、v始=514.8x103(G)t末=22.958(h),v末=677.6x103(G)2 由于要求的是水箱流量与时间的关系，因此须由上表的数据计算出相邻时间区 间的中点及在时间区间内水箱中流出的水的平均速度：平均流速=(区间左端点的水量一区间右端点的水量)/区间中点值得下表:时间区间的中点值(h)1平均流量(x103G/h)I0.4606H14.01.38212.02.39610.03.411nJ9.64.425|9.65.4398.9一6.459.67.468(8.98.44810.09.474/10.45/10.941/11.49J18.612.49”0.0 113.4219.014.43

6、16.015.4416.0 116.37(16.017.3814.0418.4914.019.5016.020.4015.021.43/22.49/23.42/24.4314.025.4512.0做出散点图如图 15-1:图 15-1 散点图从图中可以看出数据分布不均匀，局部紧密，因此不能采用插值多项式处理 数据，而用曲线拟合的最小二乘法。(四)计算过程1.算法：第 1 步:输入数据Xi, yi；第 2 步:进行拟合；第 3 步:作出散点图；第 4 步:作出拟合函数图；第 5 步:进行误差估算。2.实现：在算法步 2 中使用 Fit函数，步 3、步 4 使用 Plot，步 5 选用 Integ

7、rate 函数。3.误差估计：误差估算时，由于水泵的灌水速度为一常数，水箱中水的体积的平均变化速 度上应近似等于水泵的灌水速度 P 减去此段时间从水箱中流出的平均速度。即tf(t)dtP V/ t -t此处 f (t)在4t 区间的两端点间进行积分。如果此模型确实准确地模拟了这些数据，那么在不同的灌水周期中，按此模 型计算出的水泵灌水速度应近似为常数。下面通过水泵开始和停止工作的两段区 间，即 t8.968，10.926及 t 20.839, 22.958来进行检验。5第一段:对应于 t始=8.968 ( h), t末=10.926 (h)水量分别为v始=514800 ( G), v末=677

8、600 (G)故Vi=677600-514800=162800 (G)ti=10.926-8.986=1.958 (h)丄=83150 (G/h)t1第二段：对应于t始=20.839 ( h),t末=22.958 (h)水量分别为v始=514800 (G), v末=677600 (G)所以V2=677600-514800=162800 (G)t2=22.958-20.839=2.119( h)V21 2 3=76830 (G/h)通过水泵开始和停止工作的两段时间检验水泵灌水速度应近似为常数；其中 由1 ,x, x2, x3 *,，x8拟合的函数 f (t)所产生的误差为 8.217%，由1 ,

9、 x3, x5 *, sin(0.1x) ,cos(0.1x)拟合达到 8.224%。由此可见如选择不同的基函数,将得 不同的误差。但是只要基函数选择恰当，所产生的误差也可以保持为相对稳定最 小常数来支持该模型。同时，一旦确定了最佳 f (t),我们便可通过 Integrate函数估算出一天的用 水总量，从而根据常规每 1000 人用水量来推测出该地区的人口数，另外，还可求 得水箱的平均流速。评价1、优点：2任意时刻从水箱中流出的水速都可通过该模型计算出来；3可推测几天的流速；P1=83150+1t110.9268.968Wdt671P2=76830+120.839l2P卩2F2(五)结果分析

10、22.958f (t )dt8fx=FitL , 1 , xA3 , xA5 ,graph1=ListPlotL , Display Function Iden titygraph2=ListPlotfx , x , 0.46 , 25.45 , DisplayFunction Identity;Showgraph1 , graph2 , Display Fun cti on $Display Function ,PlotRa ngev 仁 677600-514800;t2=10.926-8.968;m 仁 v1/t1 ；v2=677600-514800;t 仁 22.958-20.839;m

11、2=v2/t2 ;p1=m1+lntegratefx , x , 8.968 , 10.926/t1p2=m2+lntegratefx , x , 20.839 , 22.958/t2%= （p1-p2） p2运行结果为：112.541-0.0198509x+0.0000156083f-97.968Cos0.1x-33.1108Si n0.1x 83159.976840.80.0822362（六）模型评价该模型数学概念简单，并且容易实现，任意时刻从水箱中流出水的速 度都可通过该模型计算出来，可以推测速度.但数据太少，只能参照一天的 数据.另外，如果(2)通过考虑体积测量的差异建模，该作法包含着

12、某种不准确性。源程序：L=0.460 ,4.425,8.448,14.43,18.49,25.45,14.0, 1.382, 12.0, 2.396, 10.0, 3.411, 9.6,9.6, 5.439, 8.9, 6.45, 9.6, 7.468, 8.9,10.0,16.0,14.0,12.011.49,15.44,19.50,18.6, 12.49, 20.0,16.0, 16.37, 16.0,16.0, 20.40, 15.0,13.42,17.38,24.43,19.0,14.0,14.0,Sin0.1x , Cos0.1x,x图 15-2 水箱水流量拟合图9知道水泵的灌水速度，就能更准确地估算水泵灌水期间水 的流速.（七）模型的分析与改进1）模型评价1、模型的优缺点优点：1.模型可以用于有一个标准水箱的小镇使用，容易推广；2.模型用到的知识简单易懂，模型容易完成；3.模型不仅提供了水流量及一天用水量的较为准确的估计，还可以估计任何时刻的水流量，包括水泵工作时的水流量缺点：无法准确估计结果的误差（八）参考文献1 飞思科技