估计水塔水流量的求解模型.docVIP

估计水塔水流量的求解模型 摘要 由所给的题目可知，本问题是一个关于如何计算居民用水的问题，由题目给出的表格，可知不同时刻的水位，根据所要求的不同时刻水位的不同入手，此计算问题就可以转化为插值或拟合问题。这里主要考虑采用插值的方法，可以利用MATLAB软件进行插值和曲线拟合计算并解决一些具体的实际问题。根据题目建立模型并采用插值的方法进行求解，推算出任何时刻（包括水泵正供水时）从水塔流出的水流量，及一天的总用水量。 关键词：用水规律与水泵的工作功率 原始数据 用水规律与水泵的工作功率 问题重述 1.1基本情况 某居民区有一供居民用水的圆柱形水塔，一般可以通过测量其水位来估计水的流量。面临的困难是，当水塔水位下降到设定的最低水位时，水泵自动启动向水塔供水，到设定的最高水位的时候停止供水，这段时间无法测量水塔的水位和水泵的供水量。通常水泵每天供水一两次，每次约3h. 已知水塔是一个高为12.2m，直径为17.4m的正圆柱。 1.2 所要解决的问题 现在需要了解该居民区用水规律与水泵的工作功率。按照设计，当水塔的水位降至最低水位，约8.2米时，水泵自动启动加水；当水位升高到一个最高水位，约10.8米时，水泵停止工作。 可以考虑采用用水率(单位时间的用水量)来反映用水规律，并通过间隔一段时间测量水塔里的水位来估算用水率。表1是某一天的测量记录数据，测量了28个时刻，但是由于其中有4个时刻遇到水泵正在向水塔供水，而无水位记录(表中用符号//表示)。所要解决的问题就是，要估计任何时刻（包括水泵正供水时）从水塔流出的水流量，及一天的总用水量。 时刻t（h） 0 0.92 1.84 2.95 3.87 4.98 5.90 7.01 7.93 8.97 水位(cm) 968 948 931 913 898 881 869 852 839 822 时刻t(h) 9.98 10.93 10.95 12.03 12.95 13.88 14.98 15.90 16.83 17.93 水位(cm) // // 1082 1050 1021 994 965 941 918 892 时刻t(h) 19.04 19.96 20.84 22.02 22.96 23.88 24.99 25.91 ? ? 水位(cm) 866 843 822 // // 1059 1035 1018 ? ? 表1水位测量记录（符号//表示水泵启动） 问题背景 1991年的美国大学生数学建模竞赛A题(AMCM1991A),由于它是水库调度、自来水管理、公共场所的人流量估计等问题的代表,因此有许多文献对其进行了研究,但一般都是采用差分与拟合的方法。而由于居民何时用水是无法准确的预报的,可能引起的水位的变化是随机事件,因此,可以以水容量作为随机变量,建立一个随机数学模型,不仅可以给出了水塔流量函数,同时还可以讨论水容量函数的数学期望。 现在越来越多的人开始研究居民区的用水问题，已知居民的用水是由一个圆柱形的水塔提供，要算出任何时刻（包括水泵正供水时）从水塔流出的水流量，及一天的总用水量，根据题目所给的信息和表格，需要建立合适的数学模型，推算任何时刻的用水率和水泵工作功率，从而计算出一天的总用水量。 问题分析 问题要求任意时刻的用水率，即求单位时间流出的水的体积，一般称为水流速度或流量。由于水塔是一个圆柱体，体积 可以很容易地通过水位高度h计算出来，这样在水泵不工作的时间段，水流速度就可以从体积对时间的导数计算出来，由于没有水的体积关于时间的函数表达式，而只能利用问题中给定的原始数据表公式，计算出离散的在测量时刻的体积V，因此可以考虑用差商代替微商，也就是用离散代替连续的思想。 为提高计算精度，采用二阶差商，即 由于所有数据被水泵两次供水分割成三组数据对每组数据的中间数据采用中心差商，前后两个数据不能采用中心差商，改用向前差商、向后差商或用中点公式进行差商。 1.中心差商公式： 2.向前差商公式： 3.向后差商公式： 4.中心公式： 以上分析了水泵不工作的时段，用水率的计算。对于水泵供水时段的用水率，计算难度较大，我们只好用供水时间段前后的用水率进行插值或拟合而得到。有了任何时刻的用水率，可以采用数值积分计算一天的总用水量。 模型假设与约定 为建模的需要，给出如下假设： (1)影响水从水塔中流出的流量的唯一因素是公众对水的传统要求．因为表１给出的数据没有提及任何其他的影响因素，我们假定所给数据反映了有代表性的一天，而不包括任何特殊情况，如自然灾害、火灾、水塔溢水、水塔漏水等对水的特殊要求． (2)水塔中的水位不影响水流量的大小，气候条件、温度变化等也不影响水流量．因为物理学的Torricelli定律指出：水塔的最大水流量与水位高度的平方根成正比，题目中给出水