matlab验作业.docVIP

2011-2012学年下半学期 Matlab数学实验报告 学院：经济与金融学院 系别：经济学类 实验1、水塔水流量的估计 第1部分、 实验问题 美国某州的各用水管理机构要求各社区提供用水率（以每小时多少加仑计，英制单位下，1加仑=4.54596dm3，美制单位下，1加仑=3.78533dm3）以及每天所用的总用水量，但许多社区并没有测量流入或流出当地水塔的水量的设备，而只能以每小时测量水塔的水位代替，其精度在0.5%以内。更为重要的是，无论什么时候，只要水塔中的水位下降到某一最低水位L时，水泵就启动向水塔重新充水直至某一最高水位H，但也无法得到水泵的供水量的测量数据。因此，在水泵正在工作时，不容易建立水塔中水位与水泵工作时用水量之间的关系。水泵每天向水塔充水一次或两次，每次大约2小时。试估计在任何时候，甚至包括水泵正在工作的时间内从水塔流出的流量，并估计一天的总用水量。水塔是一个垂直圆柱体，高为40英尺，直径为57 某小镇某天的水塔水位记录 时间（秒） 水位（英尺） 时间（秒） 水位（英尺） 时间（秒） 水位（英尺） 0 31.75 35932 水泵工作 68535 28.42 3316 31.10 39332 水泵工作 71854 27.67 6635 30.54 39435 35.50 75021 26.97 10619 29.94 43318 34.45 79154 水泵工作 13937 29.55 46636 33.50 82649 水泵工作 17921 28.92 49953 32.67 85968 34.75 21240 28.50 53936 31.56 89953 33.89 25223 27.87 57254 30.81 93270 33.40 28543 27.52 60574 30.12 32284 26.97 64554 29.27 第2部分、问题分析 通过数据的单位转换，将上表数据单位s , ft转换为单位h , m ，便于计算。 时间（h） 水位（m） 时间（h） 水位（m） 时间（h） 水位（m） 0 9.6769 9.98 水泵工作 19.04 8.6620 0.92 9.4788 10.93 水泵工作 19.96 8.4334 1.84 9.3081 10.95 10.8199 20.84 8.2201 2.95 9.1253 12.03 10.4998 22.01 水泵工作 3.87 8.9864 12.95 10.2103 22.96 水泵工作 4.98 8.8144 13.88 9.9573 23.88 10.5913 5.90 8.6864 14.98 9.6190 24.99 10.3292 7.01 8.5200 15.90 9.3904 25.91 10.1798 7.93 8.3877 16.83 9.1801 8.97 8.2201 17.94 8.9211 水塔的形状是一个圆柱形，它的横截面面积是一个常数。水泵在不工作的时段，水流量容易通过水位对时间的变化率得到。由于水塔截面积是常数S ,为简单起见，计算中将流量定义为单位时间流出的水的高度，即水位对时间变化率的绝对值（水位是下降的），最后给出结果时再乘以S即可。即： 水位对时间的变化率(流量): 任何时刻的流量: 水泵供水时段的流量只能靠供水时段前后的流量拟合得到。我们可以先拟合出水位时间的函数，流量就是对其函数的求导结果。 第3部分、程序设计 1、拟合第1时段的水位，并得出流量 设t、h为已输入的时刻和水位测量记录，实现如下： t=[0,0.92,1.84,2.95,3.87,4.98,5.90,7.01,7.93,8.97,10.95,12.03,12.95,13.88,14.98,15.90,16.83,17.94,19.04,19.96,20.84,23.88,24.99,25.91]; h=[ 9.6769,9.4788,9.3081,9.1253,8.9864,8.8144,8.6864,8.5200,8.3877,8.2201,10.8199,10.4998,10.2103,9.9573,9.6190,9.3904,9.1801,8.9211,8.6620,8.4334,8.2201,10.5913,10.3292,10.1798]; f1=polyfit(t(1:10),h(1:10),5); b1=polyder(f1); tm1=0:0.01:8.97; g1=-polyval(b1,tm1); 2、拟合第2时段的水位，并得出流量 f2=polyfit(t(11:21),h(11:21),5); ; b2=polyder(f2); tm2=10.95:0.01:20.8