计算水塔水流量

1、一、课程设计目的：1训练学生灵活应用所学数值分析知识，独立完成问题分析，结合数值分析理论知识，编写程序求解指定问题。 2初步掌握解决实际问题过程中的对问题的分析、系统设计、程序编码、测试等基本方法和技能；3提高综合运用所学的理论知识和方法独立分析和解决问题的能力；4.训练用数值分析的思想方法和编程应用技能模拟解决实际问题，巩固、深化学生的理论知识，提高学生对数值分析的认知水平和编程水平，并在此过程中培养他们严谨的科学态度和良好的工作作风二、课程设计任务与要求：课程设计题目：计算水塔的水流量【问题描述】 某居民区的民用自来水是由一个圆柱形的水塔提供，水塔高12.2米，直径17.4米。水塔是由水泵

2、根据水塔内水位高低自动加水，一般每天水泵工作两次，现在需要了解该居民区用水规律与水泵的工作功率。按照设计，当水塔的水位降至最低水位，约8.2米时，水泵自动启动加水；当水位升高到一个最高水位，约10.8米时，水泵停止工作。可以考虑采用用水率（单位时间的用水量）来反映用水规律，并通过间隔一段时间测量水塔里的水位来估算用水率，原始数据表是某一天的测量记录数据，测量了28个时刻，但是由于其中有3个时刻遇到水泵正在向水塔供水，而无水位记录。试建立合适的数学模型，推算任意时刻的用水率、一天的总用水量。进一步：可自己增加一些新的计算功能。 【问题假设】 1. 水塔中水流量是时间的连续光滑函数，与水泵工作与否

3、无关，并忽略水位高度对水流速度的影响。2. 水泵工作与否完全取决于水塔内水位的高度。3. 水塔为标准的圆柱体。体积V=PI*D*D*h/4 其中D为底面直径，h为水位高。4. 水泵第一次供水时间段为8.967,10.954，第二次供水时间段为20.839,22.958。 【实验数据】 原始数据（单位：时刻（小时），水塔中水位（米）时刻t00.9211.8432.9493.8714.9785.900水位h9.6779.4799.3089.1258.9828.8148.686时刻t7.0067.9288.9679.981110.92510.95412.032水位h8.5258.3888.220泵水

4、泵水10.82010.500时刻t12.95413.87514.98215.90316.82617.93119.037水位h10.2109.9369.6539.4099.1808.9218.662时刻t19.95920.83922.01522.95823.88024.98625.908水位h8.4338.220泵水10.82010.59110.35410.180【实现提示】 由问题的要求，关键在于确定用水率函数，即单位时间内用水体积，记为f(t)，又称水流速度。如果能够通过测量数据，产生若干个时刻的用水率，也就是f(t)在若干个点的函数值，则f(t)的计算问题就可以转化为插值或拟合问题。【实现

5、要求】1、在处理每个题目时，要求分别从数据处理阶段和程序设计阶段两个主要阶段实现课程设计，详细的通过文字以及插图等形式，按需求分析、数据处理、算法设计、代码、计算结果和程序执行的截图等若干步骤完成题目，最终写出完整的分析报告。前期准备工作完备与否直接影响到后序上机调试工作的效率。在程序设计阶段应尽量利用已有的标准函数，加大代码的重用率。 2、设计的题目要求达到一定工作量，并具有一定的深度和难度。3、程序设计语言推荐使用C/C+，程序书写规范，源程序需加必要的注释;4、每位同学需提交可独立运行的程序；5、每位同学需独立提交设计报告书（每人一份），要求编排格式统一、规范、内容充实；6、课程设计实践

6、作为培养学生动手能力的一种手段，单独考核。三、课程设计说明书【需求分析】根据一组关于若干个时刻的水塔内水位高度的数据，经过必要的处理，得到一组若干个时刻的用水率的数据，然后根据这一组数据提出假设，建立合适的数学模型，通过C语言程序将该模型的一些参数求出来，最后得到一条能够拟合这一组若干个时刻的用水率的光滑的函数，进而推算出任意一个时刻的用水率，并计算出一天的总用水量。如果还有余力可以增加更多的实际计算要求。【数据处理概要】1. 关于原始数据中有三个时刻水泵是在泵水的处理：由问题假设可以知道水塔中水流量是时间的连续光滑函数，与水泵工作与否无关，如果将这三个时刻的数据忽略不计，是不影响设计结果的。

7、所以设计中忽略这三个时刻的数据。2. 关于确定求解哪些时刻的用水率问题：我采用每两个时刻的中间时刻和这两个时刻之间的时间段内的平均用水率作为最后的若干时刻的用水率。3. 关于数据处理过程中三个无数据时刻的处理：按照处理需要，适当的舍弃边缘的数据和错误的结尾数据（比如要计算两个值的差，最后一项后面已经没有数据了，系统会在该数据后产生一个零，所以是错误的数据）。【详细设计】本设计报告按照我的设计步骤完成：1、将所给的数据中的时刻输入到Excel表格中，如下图所示：2、计算出每两个相邻的时刻之间的时间差（表格列B(n)= A(n+1)- A(n)），然后除以2（表格列C(n)=B(n)/2）得到相邻

8、两个时刻的中间时刻：3、将原始数据中的水位数据输入到该表格中：4、根据公式：V=PI*D*D*h/4计算各个时刻水塔中水的体积（表格列E(n)=PI*D*D*表格列D(n)/4），然后求出各个时间段内用水的体积（表格列F(n)=D(n+1)-D(n)）：5、求中间时刻的用水率，即求该时间段（时刻差）内的平均用水量（表格列G(n)=F(n)/B(n)）：6、将中间时刻和中间时刻的用水率复制到另一张表中，以便处理（这时共产生22组数据）：到这里所有的数据处理已经结束。接下来是研究数据所包含的规律，根据数据提出适当的数学模型，拟合一条光滑的连续的曲线。7、在Excel中有拟合曲线的工具，所以我先用E

9、xcel中的工具对这组数据进行了拟合*画出散点图：*确定横纵坐标的意义：*形成散点图： *通过添加趋势线，来拟合光滑曲线：*得到拟合的函数曲线，以及表达式在Excel中就可以很方便的拟合出比较符合要求的曲线，并且可以得到函数表达式。然后开始尝试利用编程解决这个问题，也就是说用编程的方法得到相同的函数表达式！最小二乘法通常用于曲线拟合，所谓拟合是指已知某函数的若干离散函数值f1,f2,fn，通过调整该函数中若干待定系数f(1, 2,3), 使得该函数与已知点集的差别(最小二乘意义)最小。如果待定函数是线性，就叫线性拟合或者线性回归(主要在统计中)，否则叫作非线性拟合或者非线性回归。表达式也可以是

10、分段函数，这种情况下叫作样条拟合。而插值是指已知某函数的在若干离散点上的函数值或者导数信息，通过求解该函数中待定形式的插值函数以及待定系数，使得该函数在给定离散点上满足约束。插值函数又叫作基函数，如果该基函数定义在整个定义域上，叫作全域基，否则叫作分域基。如果约束条件中只有函数值的约束，叫作Lagrange插值，否则叫作Hermite插值。从几何意义上将，拟合是给定了空间中的一些点，找到一个已知形式未知参数的连续曲面来最大限度地逼近这些点 而插值是找到一个(或几个分片光滑的)连续曲面来穿过这些点。可以看出利用最小二乘法拟合比较方便，我采用的函数形式为多项式函数，最高次为5。拟合部分代码执行结果

11、如下：经过比较，可以发现在误差允许的范围内通过Excel拟合出来的函数和通过编程而拟合出来的多项式是相统一的。然后再将该函数应用到Excel中，利用程序代码所求出来的多项式，画出一条光滑连续的曲线，如果这条曲线和最上面在Excel中，纯粹用Excel拟合出来的曲线相同的话表明拟合成功！1、按照递增的顺序，将从0-24的数据输入到Excel中，如下图所示。并在第二列中输入公式：=-0.00015*A2*A2*A2*A2*A2+0.*A2*A2*A2*A2-0.*A2*A2*A2+4.*A2*A2-22.*A2+66.02523，如下图所示 。将数据填充完毕：利用和上面的相同的方法，直接选择光滑的

12、曲线图，而不是选择散点图：然后就画出了一条曲线： 经过比较，可以知道在误差允许的范围内数学模型已经成功建立！到这里该问题的数学模型已经建立好了(采用编程所得的结果进行建模和进一步的计算)：f(t)=-0.00015t5+0.t4-0.t3+4.t2-22.t+66.最后再通过编程计算任意时刻的用水率、一天的总用水量和水泵的工作功率：任意时刻的用水率的算法简单，只要给出任意时间，根据上面的式子就可以求出；一天的总用水量为上式在0-24内的积分V=024f(t)dt=-0.t6+0.t5-0.t4+1.t3-11.t2+66.02523t+C024水泵的工作功率的计算方法是用一次的总泵水量除以泵水

13、时间水泵泵一次水泵水量为：V=2*PI*D*D*(h2-h1)/4=2*3.*17.4*17.4*(10.82-8.22)/4泵水时间为：T=t2-t1=(10.954-8.967)+(22.958-20.839)所以水泵的工作功率为：P=Vg(h2-h1)/2T其中为水的密度(kg/m3)，V为上面的体积(m3)，g为重力单位（9.8N/kg）。然后将程序加以修改，将上述的三个求解功能加进去就完成了本次设计。最终源代码：#include#include#includevoid yongShuiLv();void yongShuiLiang();void gongLv();void Appro

14、x(float,float,int,int,float); void main()int i;float a6;/多项式的待定系数个数float x22=0.4605,1.382,2.396,3.41,4.4245,5.439,6.453,7.467,8.4475,11.493,12.493,13.4145,14.4285,15.4425,16.3645,17.3785,18.484,19.498,20.399,23.419,24.433,25.447;/22个时刻的值float y22=51.,44.,39.,36.,36.,33.,34.,35.,38.44874,70.,74.,70.,

15、60.,62.,58.,55.,55.,59.,57.,59.,50.,44.;/22个时刻所对应的用水率Approx(x,y,22,5,a);for(i=0;i=5;i+)printf(a%d=%fn,i,ai);printf(拟合多项式为：nf(t)=(%f)*t*t*t*t*t+(%f)*t*t*t*t+(%f)*t*t*t+(%f)*t*t+(%f)*t+(%f)n,a5,a4,a3,a2,a1,a0);yongShuiLiang();gongLv();char s;while(1)/通过判断用户输入的字符是Y/N，来控制是否继续求某一时刻的用水率cout您是否还要查询某时刻的用水率（

16、Y/N）：s;switch(s)case Y:yongShuiLv();break;case N:break;if(s=N)break;void Approx(float x,float y,int m,int n,float a)int i,j,t;float *c=new float(n+1)*(n+2);float power(int,float);void ColPivot(float *,int,float);for(i=0;i=n;i+)for(j=0;j=n;j+)*(c+i*(n+2)+j)=0;for(t=0;t=m-1;t+)*(c+i*(n+2)+j)+=power(i+

17、j,xt);*(c+i*(n+2)+n+1)=0;for(j=0;j=m-1;j+)*(c+i*(n+2)+n+1)+=yj*power(i,xj);ColPivot(c,n+1,a);delete c;void ColPivot(float *c,int n,float x)int i,j,t,k;float p;for(i=0;i=n-2;i+)k=i;for(j=i+1;j(fabs(*(c+k*(n+1)+i)k=j;if(k!=i)for(j=i;j=n;j+)p=*(c+i*(n+1)+j);*(c+i*(n+1)+j)=*(c+k*(n+1)+j);*(c+k*(n+1)+j)=

18、p;for(j=i+1;j=n-1;j+)p=(*(c+j*(n+1)+i)/(*(c+i*(n+1)+i);for(t=i;t=0;i-)for(j=n-1;j=i+1;j-)(*(c+i*(n+1)+n)-=xj*(*(c+i*(n+1)+j);xi=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i);float power(int i,float v)float a=1;while(i-)a*=v;return a;void yongShuiLv()/计算某时刻的用水率函数float t;cout请输入任意一个时刻，程序将就算出该时刻的用水率t;float f=-0.00015*t*t*t*t*t+0.*t*t*t*t-0.*t*t*t+4.*t*t-22.*t+66.;coutt时刻的用水率是：f立方米/小时endl;void yongShuiLiang()/计算一天居民总的用水量float v;v=-0.*24*24*24*24*24*24+0.*24*24*24*24*24-0.*24*24*24*24+1.*24*24*24-11.*24*24+66.02523*24;cout这个居民区的居民一天的用水量约为：v立方