水塔流量问题.doc

本科生课程设计报告 实习课程 数值分析 学院名称 管理科学学院 专业名称 学生姓名 学生学号 指导教师 实验地点 实验成绩 二〇 一 六 年 六 月 二〇 一 六 年 六 月　　　　　  估计水塔的水流量 摘要 水塔流量的估计是一个较为经典的数学建模问题，本问题最大的困难在于不知泵启动时水位的变化和向外水流的速度.解决该问题，先确定近似流速，利用中点数值求导公式计算出每个时间点出的流速，再利用插值与拟合计算出流速与时间的函数，对0到24小时积分可得总用水量，这是第一种方法.第二种方法，水泵没有开动时利用高度差计算用水量，水泵开动时利用积分，这样计算出的结果较为准确，2种方法比较，可得出误差. 关键词：中点数值求导；插值与拟合；积分  目录 TOC \o 1-3 \h \z \u 第1章 前言 1 1.1 内容及要求 1 1.2 研究思路及结构安排 2 第2章 模型建立与求解 3 2.1模型假设 3 2.2确定近似流速 3 2.3 确定水泵启动时的流量及总流量曲线 4 2.4确定总用水量 4 第3章 算法步骤 6 3.1 中点数值求导函数步骤及流程图 6 3.2 三次样条插值函数步骤及流程图 7 第4章 算法实现 7 4.1 程序总体结构 7 4.2 源程序清单 8 4.3 程序运行 12 第5章 误差分析 14 第6章 模型的评价和改进 15 6.1 优点 15 6.2 缺点 15 6.3 模型的改进方向 15 参考文献 15  第1章 前言 1.1 内容及要求 某地的用水管理机构要求各社区提供各个时刻的用水率以及每天所用的总用水量。但许多社区并没有测量流入或流出当地水塔的水量的设备，他们只能代之以每小时测量水塔中的水位，其误差不超过0.5%。更为重要的是，无论什么时候，只要水塔中的水位下降到最低水位L时，水泵就自动启动向水塔重新充水直到最高水位H时水泵自动停止，但也无法得到水泵的供水量的测量数据。因此，在水泵正在工作时，人们不容易建立水塔中水位与水泵工作时的用水量之间的关系。水泵每天向水塔充水一次或两次，每次约二小时。下表为某地一天中的真实的数据。 表1某天水塔水位测量记录 时刻t（秒） 0 3316 6635 10619 13937 17921 21240 水位（0.01英尺） 3175 3110 3054 2994 2947 2892 2850 时刻t（秒） 25223 28543 32284 35932 39332 39435 43318 水位（0.01英尺） 2795 2752 2697 水泵启动 水泵启动 3550 3445 时刻t（秒） 46636 49953 53936 57254 60574 64554 68535 水位（0.01英尺） 3350 3260 3167 3087 3012 2927 2842 时刻t（秒） 71854 75021 79254 82649 85968 89953 93270 水位（0.01英尺） 2767 2697 水泵启动 水泵启动 3475 3397 3340 水塔是一个高40英尺、直径57英尺的圆柱。按照设计，水塔水位降至约L=27英尺时，水泵自动启动加水；当水位升高到约H=35.5英尺米时，水泵自动停止工作。 试估计在任何时刻（包括水泵正在供水时）水从水塔流出的流量f(t)，并估计一天的总用水量。 1.2 研究思路及结构安排 本问题最大的困难在于不知道水泵启动时水位的变化和向外水流的速度。用水量等于向外水流速度乘以时间。因此，如何确定流速是解决该问题的关键。  第2章 模型建立与求解 2.1模型假设 （1）除了问题中特别说明的数据以外，其他给定的数据其测量误差不超过0.5%； （2）一天之中，任意从一个开始时刻，如从t=0起到t=24小时结束，一天开始时刻的不同不影响一天总水量； （3）管理部门不考虑水流速度的瞬间值，感兴趣的是整个一天中的用水总趋势； （4）水泵抽水的速度是均匀的； （5）假设水流的速度是连续变化的，流速可以用一条光滑的曲线近似表示，即在给定的数据点上有连续的二阶导数。 2.2确定近似流速 水泵不供水时,水流速度不难用数值方法描述，但如何描述清楚水泵向水塔供水期间的流速? 我们首先估计水泵不启动时,每个已知数据点上的流速。流速等于水量的变化除以时间的变化。假设为时刻水塔中水的容量，表示时刻水塔中水的容量为。 使用中点数值求导公式，求处的导数的近似值。水泵不启动时段共有3段，对于每段的前两个点 对于每段的后两个点 对于每段的其余点 如此，可得各时刻点上的流速。 (负