水塔水流量问题数学模型.docVIP

水塔水流量问题的数学模型 摘要 准确流量预测在良好的中越来越成为至关紧要的一个步骤，对城发展具有重要的意义本文530.1068 1229075.82 水位高度和体积的误差均不超过0.5%，在此误差范围内，本文两种方法的结果与实际基本相符。 针对本文的问题，第一种是曲线插值求解，我们分别采用三次样条插值、Lagrange插值和分段线性插值，分析比较后，三次样条插值最合适，得出水位高为537.7507cm。 第二种是曲线拟合求解，我们根据最小二乘原理，分别用3次和5次拟合多项式逼近，仔细比较后，3次拟合效果最佳，求得水位高为530.1068。 整个模型求解程序均是在MATLAB7.1中运行的。 关键字：水流量、三次样条插值、最小二乘曲线拟合、MATLAB7.1 一、问题重述 美国某州的各用水管理机构要求各社区提供以每小时多少加伦计算的用水率以及每天所用的总水量。许多社区没有测量流入或流出水塔的水量装置，他们只能代之以每小时测量水塔中的水位，其误差不超过0.5%。更重要的是，当水塔中的水位下降到最低水位L时水泵就启动向水塔输水，直到最高水位H，但也不能测量水泵的供水量。因此，当水泵在输水时不容易建立水塔中水位和水泵工作时用水量之间的关系。水泵每天输水一次或两次，每次约二小时。 试估计任意时刻（包括水泵在输水工作的时候）从水塔流出的流量f（t），并估计一天的总水量。附表给出了某个小镇一天中真实的数据。 表1（见附件表1）给出了从第一次测量开始的以秒为单位的时刻，以及该时刻的高度单位为百分之一英尺的水位测量值。例如，3316秒后，水塔中水位达到31.10英尺，水塔是一个高为40英尺，直径为57英尺的正圆柱。通常当水塔降至约27.00英尺时水泵开始工作，当水位升到35.50英尺时水泵停止工作。 二、问题分析 本文问题所指流量可视为单位时间内流出水的体积，一天的时间为0到23.88小时止（即忽略第3未供水时段）。由于水塔是正圆柱形，横截面积是常数，所以在水泵不工作时段，流量很容易根据水位相对时间的变化算出。问题的难点在于如何估计水泵供水时段的流量。 水泵供水时段的流量只能靠供水时段前后的流量经插值或拟合得到。作为用于插值或拟合的原始数据，我们希望水泵不工作时段的流量越准确越好。这此流量大体上可由两种方法计算，一是直接对表1中的水量用数值微分算出各时段的流量，用它们拟合其它时刻或连续时间的流量；二是先用表中数据拟合水位一时间函数，求导数即可得到连续时间的流量。 有了任何时刻的流量，就不难计算一天的总用水量。其中，水泵不工作时段的用水量可以由测量记录直接得到，由表1中下降水位乘以水搭的截面积就是这一时段的用水量。这个数值可以用来检验数据插值或拟合的结果。 符号说明 t：测量的时刻； h：水位的高度； v：水塔中水的体积； 模型假设 假设一：流量只取决于水位差，与水位本身无关，故由物理学中Torriceli定律：小孔流出的液体的流速正比于水面高度的平方根。题目给出水塔的最低和最高水位分别是8.1648m（设出口的水位为零）因为sqrt，约为1，所以可忽略水位对流速的影响。 假设二：将流量看作时间的连续光滑函数，为计算简单，不妨将流量定义成单位时间流出水的高度，即水位对时间变化率的绝对值（水位是下降的），水塔截面积为,得到结果后乘以s即可. 假设三：水泵工作起止时间由水塔的水位决定。水泵工作时不维修，也不中途停止工作。水泵冲水的水流量远大于水塔的水流量。 假设四：表1中水位数据取得的时间准确在1秒之内。 假设五：水塔的水流量与水泵状态独立，并不因水泵工作而增加或减少水流量的大小。 三、模型求解 1．拉格朗日插值法 假设取区间上的n+1个结点，并且已知函数f(x)在此点的函数值 x x0 x1 x2 … xn f(x) f(x0) f(x1) f(x2) … f(xn) 现在求一个次数不超过n的多项式 使得满足条件 这种差值方法称为n次多项式插值(或称代数插值)，利用拉格朗日插值方法可得 ， 其中 上述多项式称为n次拉格朗日(lagrange)插值多项式，函数称为拉格朗日插值基函数。 特别地，当n＝1，2时，n次拉格朗日(lagrange)插值多项式即为线性差值多项式和抛物插值多项式。 2．三次样条插值 分段线性插值，具有良好的稳定性和收敛性，但光滑性较差。在数学上若函数（曲线）的k阶导数存在且连续，则称该曲线具有k阶光滑性。易见，分段线性插值不光滑，这影响了它在某些工程技术实际问题中的应用。例如：在船体、飞机等外形曲线的设计中，不仅要求曲线连续而且还要求曲线的曲率连续，这就要求插值函数具有连续的二阶导数。为解决这一类问题，就产生了三次样条插值。 所谓样条（Spline），本来是指一种绘图工具，它是一种富有弹性的细长木条，在飞机或轮船制造过程