用灰色理论分析的水箱水流量问题.docxVIP

水箱水流量问题一、问题重述许多供水单位由于没有测量流入或流出水箱流量的设备，而只能测量水箱中的水位，估计在任意时刻t流出水箱的流量f(t)。给出原始数据表以及单位换算之后（长度m，时间h）的表。其中泵水用-1标记。时间/s时间/h水位/10-2E体积/m3时间/s时间/h水位/10-2E体积/m300.00 31752.24 4463612.40 33502.36 33160.92 31102.19 4995313.88 32602.30 66351.84 30542.16 5393614.98 31672.23 106192.95 29942.11 5725415.90 30872.18 139373.87 29472.08 6057416.83 30122.13 179214.98 28922.04 6455417.93 29272.07 212405.90 28502.01 6853519.04 28422.01 252237.01 27951.97 7185419.96 27671.95 285437.93 27521.94 7502120.84 26971.90 322848.97 26971.90 7925422.02 -10.00 359329.98 -10.00 8264922.96 -10.00 3933210.93 -10.00 8596823.88 34752.45 3943510.95 35502.51 8995324.99 33972.40 4331812.03 34452.43 9327025.91 3340二、模型假设（1）影响水箱流量的唯一因素是该区公众对水的普通需要；（2）水泵的灌水速度为常数；（3）从水箱中流出水的最大流速小于水泵的灌水速度；（4）每天的用水量分布都是相似的；（5）水箱的流速可用光滑曲线来近似；（6）当水箱的水容量达到开始泵水；达到停止泵水。三、符号说明t：时间V：水箱中水的体积f（t）任意时刻t水箱的流量四、模型建立与求解（一）插值拟合解法：1.首先将表格中的数据转化为标准单位制，其中时间用小时，体积用立方米。1E=0.3024m1L=7.481G2.（1）水塔中水的体积计算：式中：D为水塔的直径；h为水塔中的水位高度。（2）水塔中水流速度的估计。水流速度应该是水塔中水的体积对时间的导数，但由于没有每一时刻水体积的具体数学表达式，这里我们用插商近似导数。 由于在两个时间段，水泵向水塔供水，无法确定水位的高度，因此计算水流速度要分3段进行计算。第一段从0省道32284s，第二段从39435s到75021s，第三段从85968s到93270s3.为了得到流速的连续曲线可以使用三次样条插值处理，然后做出时间-流速的散点图。表1：插值所得空缺的流量值?时间/h9.4810.4610.9421.4322.4923.42流量/G127511571117117147141449014190图1：时间流速散点图4.曲线拟合用MATLAB进行三次样条插值拟合的曲线如下图2：时间-流速曲线5.对水流速度进行数值积分可以得到24小时的用水量为1358.4（二）GM（2,1）模型解法：1.将第一种方法中插值拟合所求出的数据代替原来文本文件中的-1时间/h9.4810.4622.4923.42流量/G127511571114490141902.设时间序列为，体积序列为将组成新的一个数组并存储在文件名为data2.txt的纯文本文件中。然后对序列建立GM（2,1）模型。计算出 时间/h原始流量/L预测值残差相对误差01.89951.8995000.461.80321.55980.24350.1351.381.55331.597-0.04370.02812.41.37951.6349-0.25540.18513.411.30461.6735-0.36880.28274.431.26451.7126-0.44810.35445.441.16431.7524-0.58810.50516.461.26451.7925-0.52790.41757.471.19241.8328-0.64050.53728.451.35011.8732-0.52310.38758.971.60271.9132-0.31050.193710.962.77371.95260.82110.29611.492.4811.99060.49030.197612.492.63472.02660.6080.230813.422.47032.05970.41060.166214.432.19362.08840.10520.047915.442.2192.11110.10790.048616.372.05852.1254-0.06680.032517