水流问题（数学建模）.docVIP

估计水塔的水流量 1 问题提出 某居民区的民用自来水是由一个圆柱形的水塔提供．水塔高12.2米，直径17．4米．水塔是由水泵根据水塔内水位高低自动加水，一般每大水泵工作两次．现在需要了解该居民区用水规律与水泵的工作功率．按照设计，当水塔的水位降至最低水位，约8．2米时，水泵自动启动加水；当水位升；高到一个最高水位，约10．8米时，水泵停止工作． 可以考虑采用用水率（单位时间的用水量）来反映用水规律，并通过间隔一段时间测量水塔里的水位来估算用水率．表4．2是某一天的测量记录数据，测量了28个时刻，但是由于其中有3个时刻遇到水泵正在向水塔供水，而无水位记录 时刻（t）/h 0 0.921 1.843 2.949 3.871 4.978 5.900 水位（t）/m 9.677 9.479 9.308 9.125 8.982 8.814 8.686 时刻（t）/h 7.006 7.928 8.967 9.981 10.925 10.954 12.032 水位（t）/m 8.525 8.388 8.220 —— —— 10.820 10.500 时刻（t）/h 12.954 13.875 14.982 15.903 16.826 17.931 19.037 水位（t）/m 10.210 9.936 9.653 9.409 9.180 8.921 8.662 时刻（t）/h 19.959 20.839 22.015 22.958 23.880 24.986 25.908 水位（t）/m 8.433 8.220 —— 10.820 10.591 10.354 10.180 试建立合适的数学模型，推算任意时刻的用水率、一天的总用水量和水泵工作功率． 2 问题分析与数据处理 由问题的要求，关键在于确定用水率函数，即单位时间内用水体积，记为f(t)，又称水流速度．如果能够通过测量数据，产生若干个时刻的用水率，也就是f(t)在若干个点的函数值，则f(t)的计算问题就可以转化为插值问题． 1．假设 1）水塔中水流量是时间的连续光滑函数，与水泵工作与否无关，并忽略水位高度对水流速度的影响． 2）水泵工作与否完全取决于水塔内水位的高度，且每次加水的工作时间为2小时 3）水塔为标准圆柱体． 考虑到假设2）结合表4．2中具体数据，推断得出 4）水泵第一次供水时间段为［8．967，10．954］，第二次供水时间段为「20．8 3 9，22．95 8］． 2．体积计算 水塔是一个圆柱体，体积为．其中D为底面直径，h为水位高度。近似地取。得到不同时刻水塔中水的体积如表4.3 时刻（t）/h 0 0.921 1.843 2.949 3.871 4.978 5.900 体积(t) 2301.1 2254.0 2213.3 2169.8 2135.8 2095.9 2065.4 时刻（t）/h 7.006 7.928 8.967 9.981 10.925 10.954 12.032 体积(t) 2027.1 1994.6 1954.6 —— —— 2572.9 2496.8 时刻（t）/h 12.954 13.875 14.982 15.903 16.826 17.931 19.037 体积(t) 2427.8 2362.7 2295.4 2237.3 2182.9 2121.3 2059.7 时刻（t）/h 19.959 20.839 22.015 22.958 23.880 24.986 25.908 体积(t) 2005.3 1954.6 —— 2572.9 2518.4 2462.0 2420.7 ? ? 3．水流速度的估算 水流速度应该是水塔中水的体积对时间的导数（微商）由于没有水的体积关于时间的函数表达式，而只有一个离散的函数值表4．3，因此考虑用差商代替微商，这也是离散反映连续的常用思想．为提高精度，采用二阶差商，即 具体地，因为所有数据被水泵两次工作分割成三组数据，对每组数据的中间数据采用中心差商，前后两个数据不能够采用中心差商，改用向前或向后差商． 中心差商公式 向前和向后差商公式 时刻（t）/h 0 0.921 1.843 2.949 3.871 4.978 5.900 流速(t)/ m/s 51.1204 47.6090 41.5072 38.2242 36.4474 34.6895 33.8858 时刻（t）/h 7.006 7.928 8.967 9.981 10.925 10.954 12.032 流速(t)/ m/s 34.