建模竞赛实例14.pptVIP

3.建模竞赛实例: 水塔流量估计 3.1. 问题的提出 有一供居民用水的正圆柱形水塔, 塔高12.2米、直径17.4米。按照设计, 水塔水位降至约8.2米时, 水泵自动启动, 水位升到约10.8米时水泵停止工作。 对水塔流量一般可通过测量其水位来估计。但面临的困难是, 当水塔水位下降到设定的最低水位8.2米时, 水泵自动启动向水塔供水, 到设定的最高水位10.8米时停止供水, 这段时间无法测量水塔的水位和水泵的供水量。 建模竞赛实例2: 煤矸石的堆积问题 3.1. 问题的提出 煤矿采煤时,会产出无用废料——煤矸石。在平原地区，煤矿不得不征用土地堆放矸石，通常矸石的堆积方法是： 架设一段与地面角度约为 的直线型上升轨道。用在轨道上行驶的矸车将矸石运到顶端后倾倒。待矸石堆高后，再借助矸石堆延长轨道，这样逐渐堆起如下图所示的一座矸石山来。 * 通常水泵每天供水1-2次, 每次约2h。某天的水位测量记录如下表所列, 试估计任何时刻包括水泵正供水时从水塔流出的水流量及一天总用水量。 3.2 问题分析 流量是单位时间流出的水的体积。由于水塔为正圆柱形, 横截面积是常数, 在水泵不工作的时段, 流量很容易从水位对时间的变化率算出, 问题是如何估计水泵供水时段的流量。 水泵供水时段的流量只能靠供水时段前后的流量拟合得到。作为用于拟合的原始数据, 希望水泵不工作时段的流量越准确越好。 大体有两种计算方法, 一是直接对表中的水位用数值微分计算出各时刻的流量, 用它们拟合其它时刻或连续时间的流量；二是先用表中数据拟合水位—时间函数, 求导数可得 连续时间的流量。一般来说数值微分的精度不高, 何况测量记录不等距, 且数值微分计算麻烦。因此选用第二种方法进行处理。 有了任何时刻的流量, 就可算出一天的总用水量。其实, 水泵不工作时段的用水量可根据测量记录直接得到, 由表可知从t=0到t=8.97(小时), 水位下降了968-822=146 cm, 乘以水塔的截面积即得该时段的用水量。该数值可用来检验拟合的结果 3.3 模型假设 1) 忽略水位对流速的影响 流量只取决于水位差, 与水位本身无关, 按Torricelli定律从小孔流出的流体的流速正比于水面高度的平方根, 水塔的最高水位和最低水位分别为10.8m和8.2m。并设出水口的水位为0m。因 , 故可忽略水位对流速的影响。 2）供水时段的假设 水泵第1次供水时段为t=9.0h 到t=11.0h, 第2次供水时段为t=20.8h到t=23.0h 。这是据最高和最低水位分别为10.8m和8.2m, 及表的水位测量记录作出的假设, 其中前3个时刻直 接取自实测数据(精确到0.1h)。, 最后1个时刻来自每次供水约2h的已知条件(从记录看, 第2次供水时段应在有记录的22.96h之后即将结束)。 3) 单位时间的供水量设为常数 水泵工作时单位时间的供水量基本为常数, 这个常数大于单位时间的平均流量。 4) 流量是对时间的连续函数 流量是单位时间流出的水的体积, 这里假设流量是对时间的连续函数。 5) 流量与水泵是否工作无关 为简化处理, 不影响间题的解决, 假设流量与水泵是否工作无关。 6) 流量定义的假设 由于水塔截面积是常数 , 为简便起见, 计算中将流量定义作单位时间流出的水的高度，即水位对时间变化率的绝对值，最后给出结果时再乘以S即可. 3.4 流量估计 1) 拟合水位—时间函数 由表记录看, 一天有两个供水时段（以下称第一供水时段和第二供水时段）, 和三个水泵不工作时段（以下称第1时段t=0到t=8.97h, 第2时段t=10.95h到t=20.84h, 和第3时段t=23h以后)。 对第1,2时段的测量数据直接分别作多项式函数拟合, 得到水位函数。为使拟合曲线比较光滑, 多项式次数不要太高, 一般为3-6次。由于第3时段只有3个测量记录, 无法对这一时段的水位作出较好的拟合。 2) 确定流量—时间函数 对于第1，2时段只需将水位函数求导数即可, 对于两个供水时段的流量, 则用供水时段前后（水泵不工作时段）的流量拟合得到, 且将第2时段流量外推, 将第3时段流量包含在第2供水时段内。 3）一天总用水量的估计 总用水量等于两个不工作时段和两个供水时段用水量之和, 它们都可由流量对时间的积分得到