1水箱水流量问题

例题估计水箱的水流量模型（1991年美国大学生数学建模竞赛的A题.问题中使用的长度单位为E（英尺，1E=30.24cm）,容积单位是G（加仑，1G=3.785L））.问题的提出：某些州的用水管理机构需估计公众的用水速度（单位:G/h）和每天的总用水量.许多供水单位由于没有测量流入或流出量的设备，而只能测量水箱中的水位（误差不超过5%）.当水箱水位低于水位L时，水泵开始工作将水灌入水箱，直至水位达到最高水位H为止.但是依然无法测量水泵灌水流量，因此，在水泵工作时无法立即将水箱中的水位和水量联系起来.水泵一天灌水1~2次，每次约2h.试估计在任一时刻（包括水泵灌水期间）t流出水箱的流量f（t），并估计一天的总用水量.表1给出了某镇某一天的真实用水数据.水箱是直径为57E,高为40E的正圆柱体.当水位落到27E以下，水泵自动启动把水灌入水箱；当水位回升至35.5E时，水泵停止工作。表1时间/s水位10-2E时间/s水位10-2E03175466363350331631104995332606635305453936316710619299457254308713937294760574301217921289264554292721240285068535284225223279571854276728543275275021269732284269779254泵水35932泵水82649泵水39332泵水859683475394353550899533397433183445932703340模型假设(1)影响水从水塔流出的流率的唯一因素是公众对水的传统要求.因为附表只给出了某一天(实际是近26小时)水塔的水位数据，并没有对这些数据的产生有影响的因素作出具体说明，我们只能假定所给数据反映了有代表性的一天，而不包括任何特殊情况，如自然灾害、火灾、水塔溢水、水塔漏水等对水的特殊要求.(2)水塔中水的水位不影响水流量的大小.据物理学的Torricelli(托里查里)定律，水塔最大水流量是与水位的高度的平方根成正比的.针对表8-1所给的数据，最大高度是35.50英尺,最小高度是27.00英尺,所以两个高度的最大水流量之比是J五50/J2700X1.15,接近于1,所以我们假定水位不影响水流量,类似地,我们假定气候条件、条件变化等也不直接影响水流量.(3)水泵工作起止时间由水塔的水位决定.我们总是假定水位大约27.00英尺时，水泵就开始工作,直到水位升至大约35.50英尺时停止工作,每次充水时间约为两小时.水泵工作性能、效率总是一定的,不因使用次数多少而变化,水泵工作时不需要维修,也不中途停止工作.当然,水泵充水的水流量远大于水塔的水流量，以保证人们对水的需求.(4)表1中水位数据取得的时间准确在1秒以内.(5)水塔的水流量与水泵状态独立,并不因水泵工作而增加或减少水流量的大小.(6)水塔的水流量曲线可以用一条光滑的曲线来逼近.这时,在每一个数据点，水流量的两阶导数是连续的，因为水的消耗是基于社区公众一天的活动,如洗澡、做饭、洗衣服等,每一个使用者的要求与整个社会的要求相比是微不足道的，而整个社会的需求是不可能同时增加或减少的，由于水的消耗的自然性，可以设想水流量曲线是一条连续光滑的曲线.问题分析与模型建立1、为方便起见，记V 表示水的容积；y.――表示时刻(单位:h)时水的容积；于(t) 表示流出水箱的水的流速(单位;G/h),它是时间的函数；P 表示水泵的灌水速度(G/h).先将表1中数据作变换，时间单位用小时(h),水位高转换成水的体积，具体数据如表2所示(V二九r2h,单位：103G,1E3=7.481G).表2时间/h水量/103G时间/h水量/103G0.606.09812.9544639.5050.921111593.6913.87558622.3241.84306583.14.9822604.5712.94972571.54615.9039598.2993.87139562.57416.8261574.9824.97806552.07417.9317558.7565.9544.05719.0375542.5297.00639533.55719.9594528.2127.92861525.34920.8392514.8498.96778514.84922.015/9.98111/22.9581/10.9256/23.88663.36710.9542677.68524.9869648.47712.0328657.6425.9083637.593根据假定，当水塔水位降至约27.00英尺时，开始充水，而水位升至约35.50英尺时停止充水，从所给数据自然无法知道水泵开始和停止工作的准确时间，但我们发现第一次充水前的最后一个测量数据是32284秒时水位为26.97英尺，可见水泵在32284位秒后不久开始充水。39332秒时水泵仍在工作，而39435秒时水为35.50英尺，水泵在这两个时刻之间停止了工作，这两个时刻的差距是103秒^0.028小时，很短的时间，所以我们可以假设水泵开始工作的时间为32284秒，结束工作的时间为39332秒，充水时间约为7048秒^1.95778小时，符合每次抽水约两小时的假设。现在再来分析第二次充水期间的数据，充水前最后一个测量数据是75021秒时水位26.97英尺，与第一次充水前水位相同，所以可以假设75012秒后水泵即开始工作，但充水后的第一个数据显示了一个与水泵应该停止工作时的不同的水位，说明水泵停止工作已有一段时间，水泵停止工作的准确时间在82649秒和85968秒之间，但这两个时刻相差近一个小时，更靠近哪能一个时刻？我们发现82649秒时水泵已工作约82649-75021=7628秒，与第一次充水时间相比，第二次充水时间也约为2小时，由此即知82649秒后水泵即停止了工作，所以我们可选取82649秒^22.95806小时为第二次充水停止的时间，并可假定此时刻的水位也为约35.50英尺，这就解决了水泵起止工作时间问题。即有第一段泵水的始停时间及水量为t广8.968(h),v^=514.8X103(G)t始=10.926(h),v始=677.6X103(G)第二段泵水的始停时间及求量为t广20.839(h),v广514.8X103(G)t:=22.958(h),v末=677.6X103(G)

2、由于要求的是水箱流量与时间的关系，因此须由上表的数据计算出相邻时间区间的中点及在时间区间内水箱中流出的水的平均速度：平均流速=（区间左端点的水量一区间右端点的水量）/区间的长度得下表：表3时间区间的中点值/h平均水流量/103G/h时间区间的中点值/h平均水流量/103G/h0.46055613.47113.415118.64661.3820811.595314.42916.04632.3963910.349815.443116.56973.410569.7347116.36515.52484.424729.4873517.378914.6775.439038.6964918.484614.67336.453199.4897419.498515.52947.46758.9008620.399315.18988.4481910.103621.4271/9.47444/22.4865/10.4533/23.419/10.9399/24.433513.451411.493518.583325.447611.809512.493619.67663、建立模型1)步骤①输入数据，/«.)｝；②作出散点图；③进行曲线拟合；④作出拟合曲线图，输出拟合函数；⑤进行误差估计。2)具体实现(以多项式拟合为例)由散点图可知，可采用多项式拟合曲线(如图)：以8次多项式拟合(其中拟合优度可达96%)，可得：f(t)=16282—7840.21+4690.912—1468.913+240.1114— 3)误差估计21.11415+1.010916—0.02484417+0.0002454818误差估算时，由于已假定水泵的灌水速度为一常数，同时知道在水泵抽水时，水箱中水的体积的平均变化速度A二应近似等于水泵的灌水速度P减去此段时间从水箱中流出的平均速度。即At此处f(t)在At区间的两端点间进行积分。如果此模型确实准确地模拟了这些数据，那么在不同的灌水周期中，按此模型计算出的水泵灌水速度应近似为常数。下面通过水泵开始和停止工作的两段区间，即te[8.968,10.926]及te[20.839,22.958]来进行检验。第一段:对应于t广8.968(h),t=10.926(h),水量分别为v始=514.8X103(G),末v=677.6X103(G). AV=始677600-514800=162800,、故1 (G)A(=10.926—8.986=1.958(h)AV1=83150At (G/h)1第二段:对应于t广20.839(h),v始=514.8X10第二段:对应于t广20.839(h),v始=514.8X103(G)t=22.958(h)，末,v=677.6X103(G)水量分别为故A匕=677600-514800=162800(G)A12=22.958—20.839=2.119(h)AV ，—2=76830At (G/h)2AV』10.926f(t)dtp=—1+8.968 =83150+14713=97863”、1At At (G/h)

二76830+23724二100554g小模型求解将t-0,4606h和t-24«4606h代入到水的流速拟合函数f(t)，我们得到这两时刻的流速分别近似为13688.8G/h和13335.6G/h,相差仅2.58%,从而可以认为f(t)能近似表达一天的用水流量.于是，一天里的用水总量近似地等于函数f(t)在24小时周期内的积分.有[0,4606+24J[0,4606+24J04606 f(t)dt=333018[1.382+24J[究？ f(t)dt=333701「6.45+24J645 f(t)dt=315285(G)(G)(G)可见误差均控制在能接受的范围之内，故该社区一天的总用水量约为:[0,4606+24J f(t)dt=333018若按常规每1000人的用水量为105000G/d,因此估计出这个地区大约有3200人.同时可得流出水箱的平均流速为：1「0.4606+24_J f(t)dt=13876240.4606 (G/h)水泵灌水的平均速度为：模型评价从建模的设想及实施过程我们可以看出有如下优点及不足之处，优点：①这模型很灵活，能被有一个正水箱的任何小镇使用，输入的数据可为任何规则逼近时间区间的水位,即时间分布可以是随机的.②模型用到的数学知识是简单易懂的，其计算过程完全可由拟合工具箱来完成，是很容易实现的。③只要有一台计算机甚至计算器，输入数据，这模型就容易完成。④模型不仅提供了水流量及一天用水量的较为准确的估计。还可以估计任何时刻的水流量，包括水泵工作时的水流量。⑤该模型可推广到用电分布的情况。缺点：①模型最大的缺点是无法准确估计结果的误差。对不同的输入数据，误差则不同。②用光滑曲线拟合的方法无法模拟真实水流量曲线的微小变化，除非时间区间取得充够小。③数据太少，精确度稍有欠缺。其中具体的计算过程如下：曲线拟合过程(还可用样条插值求解splinetool)>>t=[0.4605561.382082.396393.410564.424725.43903…6.453197.46758.4481911.493512.493613.415114.4290…15.443116.365017.378918.484619.498520.399324.4335…25.4476];>>f=[13.471011.595310.34989.734719.487358.696499.489748.9008610.103618.583319.676618.646616.0463…16.569715.524814.67714.673315.529415.189813.451411.8095]*10八3;>scatter(t,f)>plot(t,f)>>f=0.00024547 -0.0248441.0109 -21.114 240.11 -1468.8 4690.8-7839.916281>formatlongg>a=polyfit(t,f,8);>cftoolLinearmodelPoly8:fittedmodel1(x)=p1*x^8+p2*x^7+p3*x^6+p4*x^5+p5*x^4+p6*x^3+p7*x^2+p8*x+p9Coefficients(with95%confidencebounds):p1 = 0.0002455 (0.0001133, 0.0003777)p2 = -0.02484 (-0.03856,-0.01113)p3= 1.011 (0.4274, 1.594)p4= -21.11 (-34.24, -7.988)p5= 240.1 (72.52, 407.7)p6= -1469 (-2678, -259.2)p7= 4691 (72.96, 9309)p8 = -7840 (-1.58e+004, 117.9)p9 = 1.628e+004 (1.203e+004, 2.054e+004)定积分计算过程：> q=int('0.00024547*x^8-0.024844*x^7+ 1.0109*x"6 -21.114*x“5+240.11*x^4-1468.8*x^3+4690.8*x^2-7839.9*x+16281,8.968,10.926)；>q=q/1.958；>>p=int('0.00