水流量估计.doc

对估计水塔流量的探究 摘要 本文是一个关于投资收益和风险权衡规划问题的解答，问题中给出四种（或多种）投资方案供投资者选择，每种方案都有相应的收益率和风险率，我们的目的是 关键字： 一．问题重述 某居民区有一供居民用水的圆柱形水塔，一般可以通过测量其水位来估计水的流量，但面临的困难是，当水塔水位下降到设定的最低水位时，水泵自动启动向水塔供水，到设定的最高水位时停止供水，这段时间无法测量水塔的水位和水泵的供水量．通常水泵每天供水一两次，每次约两小时. 水塔是一个高12.2m，直径17.4m的正圆柱．按照设计，水塔水位降至约8.2m时，水泵自动启动，水位升到约10.8m时水泵停止工作． 表1 是某一天的水位测量记录，试估计任何时刻（包括水泵正供水时）从水塔流出的水流量，及一天的总用水量． 时刻(h) 水位(cm) 0 0.92 1.84 2.95 3.87 4.98 5.90 7.01 7.93 8.97 968 948 931 913 898 881 869 852 839 822 时刻(h) 水位(cm) 9.98 10.92 10.95 12.03 12.95 13.88 14.98 15.90 16.83 17.93 // // 1082 1050 1021 994 965 941 918 892 时刻(h) 水位(cm) 19.96 20.84 22.01 22.96 23.88 24.99 25.91 866 843 822 // // 1059 1035 1018 表1 水位测量记录（符号//表示水泵启动） 二．模型的假设 假设储水器的空气压力强稳定； 假设题目涉及的水泵功率恒定； 假设一天中温度不变或者影响不大； 居民的用水量没有突变点。 符号说明 表1未完善 符号 符号说明 赤纬角 方位角 太阳高度角 太阳时角(地球每个小时自转15°称之为时角) W 地理纬度 J 经度 L 影长 D 日期 H 杆子高度 n 某点的地方时间 L 杆影长度 四．模型的建立与求解 4.1问题一 4.1.1问题1分析 题目要求描建立每个时刻水塔流出用水流量变化的数学模型，分析水箱水位关于各个时刻的变化规律。首先通过查阅相关数学建模文献，找寻建模中可能遇到的概念、方法，再挖掘建立模型相关的方法，为之后的建模奠定坚实的基础。模型建立第一步,理清自变量（时间），因变量（流量）。先预处理数据，将水位转化为水流量，再处理流量—时间散点图，分段多项式拟合构造对应函数，使之符合题意，同时根据相应的标准进行统计、分析和构建相关的数学模型。 4.1.2数据的预处理 将数据导进MATLAB中并将水塔水位转化为水塔中体积，公式如下 ① 得出的体积与时间关系如下 表二未完善 时刻(h) 体积(cm) 0 0.92 1.84 2.95 3.87 4.98 5.90 7.01 7.93 8.97 时刻(h) 体积(cm) 9.98 10.92 10.95 12.03 12.95 13.88 14.98 15.90 16.83 17.93 时刻(h) 体积cm) 19.96 20.84 22.01 22.96 23.88 24.99 25.91 4.1.2模型1的建立 通过对题目的探讨，我们先利用MATLAB作出水体积与时间的散点图。 体积—时间散点图 图一 我们发现，图中散点可以分为五个阶段：未供水第一段，供水第一段，未供水第二阶段，供水第二阶段，未供水第三阶段。 由于供水阶段没有数据，所以先根据未供水段的散点图，利用插值拟合出未供水段的体积关于时间函数方程。根据流体力学知识可得流量是指单位时间内流经封闭管道或明渠有效截面的流体量，既是 ② 其中Q是流量，S为截面面积,v为水流速度，L为体积，t为时间； ③ 根据以上的方程式可以知道，先要拟合出流量—时间方程，再对方程求一次微分（导数），就可以得出流量—时间方程。 4.1.3步骤一：拟合方式的选择 由于供水阶段没有数据，所以先根据未供水段的散点图，利用插值拟合出未供水段的函数方程。 我们优先选择多项式拟合的