微积分模型水塔供水.pptVIP

第5讲 微积分模型 水塔流量的估计 轮船搁浅问题 本 讲 内 容 一、问题的提出 二、问题的分析 三、模型假设 四、估计流量 五、算法设计与编程 六、计算结果 七、分析与改进 问题一 水塔流量的估计 一、问题的提出 二、问题的分析 三、模型假设 四、估计流量 五、算法设计与编程 六.计算结果 问题二 轮船搁浅问题 一、问题的提出 二、问题分析 三、问题求解 第二讲练习题 假设：该海域海底是平滑的。 由于测量点是散乱分布的，先在平面上作出测量点的分布图，再利用二维插值法补充一些点的水深，然后作出海底曲面图和等高线图，并求出水深小于5英尺的海域范围。 150 -50 75 200 1. 作出测量点的分布图 %haiyugeqian2.m x=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5]; y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 –81 3 56.5 -66.5 84 -33.5]; plot(x,y,+‘，’markersize’,10); 矩形区域[75，200]×[-50，150]内测量点的分布 2. 作出海底地貌图 z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9]; [fnodes, minnq, rnw, rnq, ifail]=e01sef(x,y,z); %求插值函数的参数 nx=100; px=linspace(75,200,nx); ny=200; py=linspace(-50,150,ny); for i=1:ny for j=1:nx [pf(i,j), ifail]=e01sff(x,y,z,rnw,fnodes,px(j),py(i)); %求插值 end end meshz(px,py,pf+5) %做海底地貌图 作出海底地貌图并截取水深小于5英尺的部分区域图 矩形区域[75，200]×[-50，150]内海底地貌图1 矩形区域[75，200]×[-50，150]内海底地貌图2 矩形区域[75，200]×[-50，150]内海底地貌图3 3. 作水深低于5英尺的部分海底曲面图 * * 某居民区有一供居民用水的圆柱形水塔，一般可以通过测量其水位来估计其流量。但面临的困难是，当水塔水位下降到设定的最低水位时，水泵自动启动向水塔供水，到设定的最高水位时停止供水，这段时间无法测量水塔的水位和水泵的供水量。通常水泵每天供水一两次，每次约2h(小时)。 12.2m 17.4m 10.8m 8.2m 水塔是一个高为12.2m,直径为17.4m是正圆柱。按照设计,水塔水位降至约8.2m时,水泵自动启动,水位升到约为10.8m时水泵停止工作。 表1是某一天的水位测量记录(符号“//”表示水泵启动)，试估计任何时刻（包括水泵正供水时）从水塔流出的水流量及一天的总用水量。 1011 1035 1059 // // 822 843 866 水位 25.91 24.99 23.88 22.96 22.01 20.84 19.96 19.04 时刻 892 918 941 965 994 1021 1050 1082 // // 水位 17.93 16.83 15.90 14.98 13.88 12.95 12.03 10.95 10.92 9.98 时刻 822 839 852 896 881 898 913 931 948 968 水位 8.97 7.93 7.01 5.90 4.98 3.87 2.95 1.84 0.92 0 时刻 表1: 水位测量记录（时刻：h，水位：cm） 流量是单位时间流出的水的体积，由于水塔是正圆柱形，横截面积是常数，在水泵不工作的时段，流量很容易从水位对时间的变化率算出，问题是如何估计水泵供水时段的流量。 水泵供水时段的流量只能靠供水时段前后的流量拟合得到，作为用于拟合的原始数据，我们希望水泵不工作的时段流量越准确越好。 这些流量大体可由两种方法计算: 一是直接对表1中的水位用数值微分算出各时段的流量，用它们拟合其它时刻或连续时间的流量。 二是先用表中数据拟合水位~时间函数，求导数即可得到连续时间的流量。 有了任何时刻的流量，就不难计算一天的总用水量。 一般说来数值微分的精度不高，何况测量记录还是不等距的，数值微分的计算尤其麻烦。下面我们用第二种方