-水塔水流量估计模型与数据插值要点.pptVIP

四. 结果分析 对数据插值而言：使用三次样条插值法得到第1，2未供水时 段的用水高度145.6870和258.6547与实际记录的下降高度 146,260相差不大，说明插值结果与原始数据比较吻合。 按三次样条插值法估计出全天的用水量约为 539.4450×237.8×10=1282800.21 对数据拟合法而言：算得第1，2未供水时段的用水高度与 实际记录的下降高度相差也不大，说明拟合结果与原始数据 也比较吻合。计算出全天的用水量约为 526.7148×237.8×10=1252527.7944 上机任务 李继成（书） 1 教材207页练习1 第1题 2 教材207页练习1 第3题 教材219页练习2 第3题 二者至少选一(需要演示程序) x=0:5:100; y=[0 2.41 2.96 2.15 2.65 3.12 4.23 5.12 6.21 5.68 4.22… 3.91 3.26 2.85 2.35 3.02 3.63 4.12 3.46 2.08 0]; y1=10-y; plot(x,y1,b.,markersize,30); axis([0 100 0 10]) grid hold on t=0:100; u=interp1(x,y1,t); plot(t,u,r,linewidth, 3); S=100*10-trapz(x,y1); p=sqrt(diff(x).^2+diff(y1).^2); L=sum(p); fprintf(S=%.2f,L=%.2f\n,S,L) 从运行结果得到： 河床截面面积：S=337.15， 河床曲线长度： L=102.09 (3) 样条插值 另一方面, 为了提高河床曲线的模拟精度，可以根据已给的数据进行3次样条插值, 在此基础上，利用梯形法求积分命令trapz来计算河床截面积，同时对样条曲线加密，计算每一段连接线长度之和来近似河床曲线长度。 x=0:5:100; y=[0 2.41 2.96 2.15 2.65 3.12 4.23 5.12 6.21 5.68 4.22... 3.91 3.26 2.85 2.35 3.02 3.63 4.12 3.46 2.08 0]; y1=10-y; plot(x,y1,b.,markersize,30) axis([0 100 2 10]) grid hold on t=0:0.5:100; u=spline(x,y1,t); plot(t,u,‘r’,‘linewidth’, 3)； S=100*10-trapz(t,u); p=sqrt(diff(t).^2+diff(u).^2); L=sum(p); fprintf(S=%.2f,L=%.2f\n,S,L) 河床截面面积：S=339.50， 河床曲线长度：L=102.23. 数据插值模型试验：水塔水流量估计 一．实验问题 某一地区的用水管理机构要求各社区提供每小时以加 仑计的用水量以及每天所用水的总量。由于许多社区没 有测量流入或流出水塔的水量装置，所以他们只能通过 测量水塔每小时的水位高度来代替每小时的用水量(误 差不超过0.05)。当水塔中的水位下降到最低水位L (27.00ft)时水泵就自动启动向水塔输水, 当水塔水位 达到限定最高水位H (35.00ft)时，水泵停止供水。水 泵供水期间无法测量水泵的供水量。水泵每天输水一次 或两次,每次约二小时。 已知某一小镇的水塔是一高为40ft，直径为57ft的正圆柱， 下表记录的是某一天水塔水位的真实数据(1ft=0.3024米 (m))。试估计任何时刻(包括水泵正在输水时间)从水塔流出 的水流量和一天的用水总量。 二、问题分析 流量是单位时间内流出水的体积，由于水塔是正圆柱 形，截面面积是常数，所以，流量很容易根据水位相对 时间的变化率算出。 水泵不工作时段的用水量可以由测量记录直接得到: 由表1中记录的水位下降高度乘以水塔的截面面积就是这 一时段的用水量。 问题的难点在于:如何估计水泵供水时段的流量。 水泵供水时段的流量只能靠供水时段前后的流量经插 值或拟合得到。而作为用于插值或拟合的原始数据，我 们希望水泵不工作时段的流量数据越准确越好。这个数 值可以用来检验数据插值或拟合结果的准确性。 二、问题分析 流量大体上可由两种方法计算: 直接对表1中的用水量用数值插值方法算出各时段的流量，用它们拟合其它时刻或连续时间的流量; 先用表中数据拟合水位--时间函数,求导数即可得到连续时间的流量。 有