计算水塔水流量
水泵选型
水泵的基本参数表征泵主要性能的基本参数有以下几个：一、流量Q流量是泵在单位时间内输送出去的液体量（体积或质量）。
体积流量用Q表示，单位是：m3/s，m3/h，l/s等。
质量流量用Qm表示，单位是：t/h，kg/s等。
质量流量和体积流量的关系为：Qm=ρQ式中ρ——液体的密度（kg/m3，t/m3），常温清水ρ=1000kg/m3。
二、扬程H扬程是泵所抽送的单位重量液体从泵进口处（泵进口法兰）到泵出口处（泵出口法兰）能量的增值。
也就是一牛顿液体通过泵获得的有效能量。
其单位是N·m/N=m，即泵抽送液体的液柱高度，习惯简称为米。
三、转速n转速是泵轴单位时间的转数，用符号n表示，单位是r/min。
四、汽蚀余量NPSH汽蚀余量又叫净正吸头，是表示汽蚀性能的主要参数。
汽蚀余量国内曾用Δh表示。
五、功率和效率泵的功率通常是指输入功率，即原动机传支泵轴上的功率，故又称为轴功率，用P表示；泵的有效功率又称输出功率，用Pe表示。
它是单位时间内从泵中输送出去的液体在泵中获得的有效能量。
因为扬程是指泵输出的单位重液体从泵中所获得的有效能量，所以，扬程和质量流量及重力加速度的乘积，就是单位时间内从泵中输出的液体所获得的有效能量——即泵的有效功率：Pe=ρgQH(W)=γQH(W)式中ρ——泵输送液体的密度（kg/m3）；γ——泵输送液体的重度（N/m3）；Q——泵的流量（m3/s）；H——泵的扬程（m）；g——重力加速度(m/s2)。
轴功率P和有效功率Pe之差为泵内的损失功率，其大小用泵的效率来计量。
泵的效率为有效功率和轴功率之比，用η表示。
最近在一单位，接触一YORK的螺杆机，制冷量是1519KW，功率250KW，不清楚怎么计算冷却水系统冷却泵的流量及冷却塔的流量，请大师指点，冷却水进出口温差是5摄氏度冷冻水泵流量：主机的冷冻水额定流量的1.1～1.2倍（单台工作时1.1倍，两台并联工作时取1.2倍）。
冷却水泵流量：一般为制冷主机冷却水流量的1.1倍。
随机数学模型在估计水塔流量中的应用
随机数学模型在估计水塔流量中的应用张先波(三峡大学理学院，湖北宜昌443002)1991年的美国大学生数学建模竞赛A题(A M C M l991A)，由于它是水库调度、自来水管理、公共场所的人流量估计等问题的代表，因此有许多文献对其进行了研究，但一般都是采用差分与拟合的方法。
而由于居民何时用水是无法准确的预报的，可能引起的水位的变化是随机事件，因此，可以以水容量作为随机变量，建立一个随机数学模型，不仅可以给出了水塔流量函数，同时还可以讨论水容量函数的数学期望。
1991年的美国大学生数学建模竞赛A题：某州的各用水管理机构要求各社区提供以每小时多少加仑计的用水率以及每天所用的总水量，但许多社区并没有测量流人或流出当地水塔的水量的设备，他们只能代之以每小时测量水塔中的水位，其精度在0．5％以内。
更为重要的是，无论什么时候，只要水塔中的水位下洚至4某一最低水位L时，水泵就启动向水塔重新充水直至某一最高水位，但也无法得到水泵的供水量的测量数据。
因此，在水泵正在工作时，人们不容易建立水塔中的水位与水泵工作时的用水量之间的关系。
水泵每天向水塔充水一次或两次，每次约二小时。
试估计在任何时刻，甚至包括水泵正在工作的时间内，水从水塔流出的流量，并估计一天的总用水量。
表1给出了某个真实小镇某一天的真实数据。
表l某小镇某天的水塔水位时间水位时间水位时问水位时间水位(秒1(001英尺)(杪)(0ol菇哟(秒){。
01蓖尺l f秒)e0ol j踅足)03175252232795466363350718542767 331631102854327524905332607502126976635305432284269753936316779254水泵工作10619299435932水泵工作57254308782619水泵工作13937294739332水泵列#605743012859683475l79212892394353550645542927988533397212402850433183445685352842932703340表中以秒为单位给出开使测量的时间、水位(单位是0．01英尺j。
数学建模——水塔流量问题
数学建模——⽔塔流量问题实验⼗四⽔塔流量问题【实验⽬的】1．了解有关数据处理的基本概念和原理。
2．初步了解处理数据插值与拟合的基本⽅法，如样条插值、分段插值等。
3．学习掌握⽤MATLAB 命令处理数据插值与拟合问题。
【实验内容】某居民区有⼀供居民⽤⽔的圆形⽔塔，⼀般可以通过测量其⽔位来估计⽔的流量。
但⾯临的困难是，当⽔塔⽔位下降到设定的最低⽔位时，⽔泵⾃动启动向⽔塔供⽔，到设定的最⾼⽔位时停⽌供⽔，这段时间是⽆法测量⽔塔的⽔位和⽔泵的供⽔量。
通常⽔泵每天供⽔⼀两次，每次约两⼩时。
⽔塔是⼀个⾼⽶、直径⽶的正圆柱。
按照设计，⽔塔⽔位降到约⽶时，⽔泵⾃动启动，⽔位升到约⽶时⽔泵停⽌⼯作。
某⼀天的⽔位测量记录如表1所⽰，试估计任何时刻（包括⽔泵正供⽔时）从⽔塔流出的⽔流量，及⼀天的总⽤⽔量。
表1 ⽔位测量启⽰录（0101001111012012)(2x L )(2ξL )(ξf y )(x f n 0x 1x n x 0y 1y n y n n )(x L n )(x L n m x a 011-m x a x a m 1-m a n )(k n x L k y k n )(ξn L )(ξf )(x L n )(x f n m n )(x L n )(x f x )(x L n )(x f a 0x 1x nx b )(x P 11----i i i i y x x x x i i i i y x x x x 11----1-i x x i x i n 0x 0y 1x 1y n x n y a b )(x S k )(x S k )(x S i i y )(x S a b k n i x i y i n i x y )(x f )(x f )(x f )(11x r a )(22x r a )(x r a m m )(x r k k a k m m n k a Q∑=-ni ix f 12i)y )((10t t t t t t t t t dt3;%% ⽤差分计算t（22）和t（23）的流量S 2.8/8.>> t3=[20 t(22) t(23)];% 取第2时段20，两点和第3时段，两点>> xx3=[abs(polyval(a2,t3(1:2))),dht3]; 取第2时段20，两点和第3时段，两点的流量>> c3=polyfit(t3,xx3,3)% 拟合出第2⽔泵供⽔时段的流量函数>> tp3=::24;>> x3=polyval(c3,tp3);% 输出第2供⽔时段（外推到t＝24）各时刻的流量求第1、2时段和第1、2供⽔时段流量的积分之和，就是⼀天总⽤⽔量。
案例6 估计水塔水流量
f ( t )dt 335329 （加仑） f ( t )dt 336480 （加仑）
25.5 1.5
相差只约1%
[0，24]区间内检验
第一次充水 前总用水量 第一次充水后， 第二次充水前 总用水量 第一次充水 期间用水量 第二次充水 期间用水量
V1= 606125-514872=91253(加仑）
充水时间约为2.1189小时
3. 由Vi—ti关系产生水流量 fi—ti的关系
注：亦可以由Vi—ti关系拟合 V(t)，再求微商得到 f(t)
关于水流量 fi
Vi 1 Vi f i f (t i ) t i 1 t i V i V i 1 与 f i f (t i ) t i t i 1
水体积的误差为0.5% 用样条逼近的用水量其误差可用抽样计算得5.1%
一天 总量 误差
2 2 2 2 SV [ SV0 SV8.9678 SV p SV10.9542 SV20.8392
1
2 2 2 2 SV p SV 22.9581 SV23.88 SV[ 23.88 , 24 ] ]1 2
水泵工作的时间为32284秒（8.9678 小时）； 水泵结束时间为39435秒（10.9542小时）； 充水时间约为1.9864小时
水泵工作的时间为75021秒（20.8392 小时），水 位26.97英尺 第 二 次 充 水 水泵结束时间为82649秒（22.9581小时）, 补充水位35.50英尺
水流量值（表3）
时
（小时）
间
水 流 量
（加仑/小时）
时
（小时）
间
水流量
（加仑/小时）
时
水冷空调水流量计算(水塔水泵)选型参考表
非空调标准进水温度 （℃） 31 32 33
非空调标准出水温度 （℃） 36 37 38
冷却水流量 （m3/h） 0 0 0
冷却水泵、水塔参考流量 （m3/h）
0 0 0
2、水冷空调冷冻水
请输入冷量 空调标准参数 （RT） 计算参数 0 （Kw） 0 空调标准进水温度 （℃） 12 空调标准出水温度 （℃） 7 冷冻水流量 （m3/h） 0 冷冻水泵参考流量 （m3/h） 备ቤተ መጻሕፍቲ ባይዱ注
0
注：1、上表输入冷量每次只输入冷吨（RT）或（Kw）。 2、上述计算只供参考，冷却水塔、冷却水泵计算还要根据当地室外计算参数确定。 3、上表为单台水泵值，如果N台水泵并联，每台流量加1.1*N倍。
水冷空调水流量计算（水塔水泵）选型参考表
1、水冷空调冷却水
请输入冷量 空调标准参数 （RT） 计算参数 0 （Kw） 0 空调标准进水温度 （℃） 30 空调标准出水温度 （℃） 35 冷却水流量 （m3/h） 0 冷却水泵、水塔参考流量 （m3/h） 备 注
0
请输入冷量 空调非标准参数 （RT） 计算参数 计算参数 计算参数 0 0 0 （Kw） 0 0 0
冷却塔流量计算
冷却塔是水与空气进行热交换的一种设备，它主要由风机、电机、填料、播水系统、塔身、水盘等组成，而进行热交换主要由在风机作用下比较低温空气与填料中的水进行热交换而降低水温。
水塔的构造及设计工况在说明书上有注明，而我们现在采用的水吨为单位是国际上比较常用的单位。
在计算选型上比较方便，另冷却塔在选型上应留有20%左右的余量。
以日立RCU120SY2 为例：冷凝：37℃蒸发：7 ℃蒸发器：Q = 316000 Kcal/h Q = 63.2m3/h冷凝器：Q = 393000 Kcal/h Q = 78.6m3/h这些在日立的说明书上可以查到；如选用马利冷却塔则：78.6×1.2 = 94.32 m3/h（每小时的水流量）选用马利SR-100 可以满足（或其它系列同规格的塔，如SC-100L）在选用水泵时要在SR-100 的100 吨水中留有10%的余量，在比较低的扬程时可选用管道泵，在扬程高时则宜选用IS 泵。
100×1.1=110 吨水/小时选用管道泵GD125-20 可以满足；而在只知道蒸发器Q=316000Kcal/h 时，则可以通过以下公式算出需要多大的冷却塔：316000×1.25（恒值）= 395000 Kcal/h,1.25——冷凝器负荷系数395000÷5 = 79000 KG/h = 79 m3/h79×1.2（余量） = 94.8m3/h（冷却塔水流量）（电制冷主机—通式：匹数×2700×1.2×1.25÷5000 或冷吨×3024×1.2×1.25÷5000= 冷却塔水流量m3/h）冷却塔已知基它条件确定冷却塔循环水量的常用公式：a. 冷却水量=主机制冷量（KW）×1.2×1.25×861/5000（m3/h）b. 冷却水量=主机冷凝器热负荷（kcal/h）×1.2/5000（m3/h）c. 冷却水量=主机冷凝器热负荷（m3/h）×1.2（m3/h）d. 冷却水量=主机制冷量（冷吨）×0.8（m3/h）e. 冷却水量=主机蒸发器热负荷（kcal/h）×1.5×1.25/5000（m3/h）f. 冷却水量=主机蒸发器热负荷（m3/h）×1.2×1.25（m3/h）g. 冷却水量=主机蒸发器热负荷（冷吨）×1.2×1.25×3024/5000（m3/h）注：以上：1.2为选型余量 1.25为冷凝器负荷系数。
冷却水塔技术参数
冷却水塔技术参数冷却水塔是用来降低工业生产中产生的热量的设备，是工业冷却系统的重要组成部分。
其技术参数包括外形尺寸、冷却能力、水流量、效率等。
下面是关于冷却水塔技术参数的详细介绍。
1.外形尺寸：冷却水塔的外形尺寸主要包括高度、长度和宽度。
尺寸的大小取决于冷却水塔的设计工况和使用需求。
一般而言，较大的外形尺寸意味着更大的冷却水塔容积和更高的冷却能力。
2.冷却能力：冷却能力是冷却水塔最重要的技术参数之一、它表示冷却水塔每小时能够降低多少热量。
冷却能力与冷却水塔的尺寸、水流量和材料等因素有关。
3.水流量：水流量是指通过冷却水塔的水的速度，通常以立方米/小时或加仑/分钟为单位。
水流量的大小对冷却能力有直接影响，较大的水流量意味着更高的冷却能力。
4.效率：冷却水塔的效率表示冷却塔转换热量的能力，即输入的热量与输出的冷却水之间的比例关系。
一般来说，冷却水塔的效率越高，其冷却能力也就越大。
冷却塔的效率与设计、材料和操作条件等因素有关。
5.降温效果：冷却水塔的降温效果是指冷却水在通过冷却塔后与环境空气接触所能达到的温度差。
降温效果对工业生产的冷却效果有直接影响，一般来说，降温效果越好，冷却水塔的性能越好。
6.噪音：冷却水塔产生的噪音是其技术参数之一、噪音的大小主要取决于冷却水塔的设计和运行情况。
一般来说，为了保证工业生产环境的安静和舒适，冷却水塔的噪音要尽量低。
7.耐久性：冷却水塔的耐久性是指其在工业生产环境中的使用寿命。
耐久性与冷却水塔的材料选择、制造工艺和使用条件等有关。
一般来说，耐久性越好的冷却水塔使用寿命越长，需求更少的维修和更低的维护成本。
8.压力损失：冷却水塔的压力损失是指冷却水通过冷却塔时所遇到的阻力。
压力损失对于水流量和冷却能力的影响很大。
较小的压力损失意味着更高的水流量和更好的冷却能力。
9.能源消耗：冷却水塔的能源消耗是指它在工作期间消耗的能源量。
能源消耗与冷却塔的设计、工况和材料等因素有关。
关于水塔设计的计算及设计参数
关于水塔设计的计算及设计参数
针对西黑山水质监测站房自来水供水压力不稳的情况，提出以下解决方案及相关设计参数和要求。
一、概述
满足供水需要，
二、理论计算
三、原理图和系统图
四、经济成本概算
五、
西黑山水质监测站自来水供水采用四分管，需求流量及压力应为0.2MPa和0.5L/s。
按一天四次取样监测，每次监测前和监测后各对管道冲洗一次，每次用水10分钟，间隔10分钟，需要水量共为600升，每次用水后满足应开启电磁阀。
故需要设计水罐容积不小于600升。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、课程设计目的：1．训练学生灵活应用所学数值分析知识，独立完成问题分析，结合数值分析理论知识，编写程序求解指定问题。
2．初步掌握解决实际问题过程中的对问题的分析、系统设计、程序编码、测试等基本方法和技能；3．提高综合运用所学的理论知识和方法独立分析和解决问题的能力；4.训练用数值分析的思想方法和编程应用技能模拟解决实际问题，巩固、深化学生的理论知识，提高学生对数值分析的认知水平和编程水平，并在此过程中培养他们严谨的科学态度和良好的工作作风二、课程设计任务与要求：课程设计题目：计算水塔的水流量【问题描述】某居民区的民用自来水是由一个圆柱形的水塔提供，水塔高12.2米，直径17.4米。
水塔是由水泵根据水塔内水位高低自动加水，一般每天水泵工作两次，现在需要了解该居民区用水规律与水泵的工作功率。
按照设计，当水塔的水位降至最低水位，约8.2米时，水泵自动启动加水；当水位升高到一个最高水位，约10.8米时，水泵停止工作。
可以考虑采用用水率（单位时间的用水量）来反映用水规律，并通过间隔一段时间测量水塔里的水位来估算用水率，原始数据表是某一天的测量记录数据，测量了28个时刻，但是由于其中有3个时刻遇到水泵正在向水塔供水，而无水位记录。
试建立合适的数学模型，推算任意时刻的用水率、一天的总用水量。
进一步：可自己增加一些新的计算功能。
【问题假设】1.水塔中水流量是时间的连续光滑函数，与水泵工作与否无关，并忽略水位高度对水流速度的影响。
2.水泵工作与否完全取决于水塔内水位的高度。
3.水塔为标准的圆柱体。
体积V=PI*D*D*h/4 其中D为底面直径，h为水位高。
4.水泵第一次供水时间段为[8.967,10.954]，第二次供水时间段为[20.839,22.958]。
【实验数据】原始数据（单位：时刻（小时），水塔中水位（米））【实现提示】由问题的要求，关键在于确定用水率函数，即单位时间内用水体积，记为f(t)，又称水流速度。
如果能够通过测量数据，产生若干个时刻的用水率，也就是f(t)在若干个点的函数值，则f(t)的计算问题就可以转化为插值或拟合问题。
【实现要求】1、在处理每个题目时，要求分别从数据处理阶段和程序设计阶段两个主要阶段实现课程设计，详细的通过文字以及插图等形式，按需求分析、数据处理、算法设计、代码、计算结果和程序执行的截图等若干步骤完成题目，最终写出完整的分析报告。
前期准备工作完备与否直接影响到后序上机调试工作的效率。
在程序设计阶段应尽量利用已有的标准函数，加大代码的重用率。
2、设计的题目要求达到一定工作量，并具有一定的深度和难度。
3、程序设计语言推荐使用C/C++，程序书写规范，源程序需加必要的注释;4、每位同学需提交可独立运行的程序；5、每位同学需独立提交设计报告书（每人一份），要求编排格式统一、规范、内容充实；6、课程设计实践作为培养学生动手能力的一种手段，单独考核。
三、课程设计说明书【需求分析】根据一组关于若干个时刻的水塔内水位高度的数据，经过必要的处理，得到一组若干个时刻的用水率的数据，然后根据这一组数据提出假设，建立合适的数学模型，通过C语言程序将该模型的一些参数求出来，最后得到一条能够拟合这一组若干个时刻的用水率的光滑的函数，进而推算出任意一个时刻的用水率，并计算出一天的总用水量。
如果还有余力可以增加更多的实际计算要求。
【数据处理概要】1.关于原始数据中有三个时刻水泵是在泵水的处理：由问题假设可以知道水塔中水流量是时间的连续光滑函数，与水泵工作与否无关，如果将这三个时刻的数据忽略不计，是不影响设计结果的。
所以设计中忽略这三个时刻的数据。
2.关于确定求解哪些时刻的用水率问题：我采用每两个时刻的中间时刻和这两个时刻之间的时间段内的平均用水率作为最后的若干时刻的用水率。
3.关于数据处理过程中三个无数据时刻的处理：按照处理需要，适当的舍弃边缘的数据和错误的结尾数据（比如要计算两个值的差，最后一项后面已经没有数据了，系统会在该数据后产生一个零，所以是错误的数据）。
【详细设计】本设计报告按照我的设计步骤完成：1、将所给的数据中的时刻输入到Excel表格中，如下图所示：2、计算出每两个相邻的时刻之间的时间差（表格列B(n)= A(n+1)- A(n)），然后除以2（表格列C(n)=B(n)/2）得到相邻两个时刻的中间时刻：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中有拟合曲线的工具，所以我先用Excel中的工具对这组数据进行了拟合*画出散点图：*确定横纵坐标的意义：*形成散点图：*通过添加趋势线，来拟合光滑曲线：*得到拟合的函数曲线，以及表达式在Excel中就可以很方便的拟合出比较符合要求的曲线，并且可以得到函数表达式。
然后开始尝试利用编程解决这个问题，也就是说用编程的方法得到相同的函数表达式！最小二乘法通常用于曲线拟合，所谓拟合是指已知某函数的若干离散函数值{f1,f2,…,fn}，通过调整该函数中若干待定系数f(λ1, λ2,…,λ3), 使得该函数与已知点集的差别(最小二乘意义)最小。
如果待定函数是线性，就叫线性拟合或者线性回归(主要在统计中)，否则叫作非线性拟合或者非线性回归。
表达式也可以是分段函数，这种情况下叫作样条拟合。
而插值是指已知某函数的在若干离散点上的函数值或者导数信息，通过求解该函数中待定形式的插值函数以及待定系数，使得该函数在给定离散点上满足约束。
插值函数又叫作基函数，如果该基函数定义在整个定义域上，叫作全域基，否则叫作分域基。
如果约束条件中只有函数值的约束，叫作Lagrange插值，否则叫作Hermite插值。
从几何意义上将，拟合是给定了空间中的一些点，找到一个已知形式未知参数的连续曲面来最大限度地逼近这些点而插值是找到一个(或几个分片光滑的)连续曲面来穿过这些点。
可以看出利用最小二乘法拟合比较方便，我采用的函数形式为多项式函数，最高次为5。
拟合部分代码执行结果如下：经过比较，可以发现在误差允许的范围内通过Excel拟合出来的函数和通过编程而拟合出来的多项式是相统一的。
然后再将该函数应用到Excel中，利用程序代码所求出来的多项式，画出一条光滑连续的曲线，如果这条曲线和最上面在Excel中，纯粹用Excel拟合出来的曲线相同的话表明拟合成功！1、按照递增的顺序，将从0-24的数据输入到Excel中，如下图所示。
并在第二列中输入公式：=-0.00015*A2*A2*A2*A2*A2+0.012193*A2*A2*A2*A2-0.366317*A2*A2*A2+4.760591 *A2*A2-22.665987*A2+66.02523，如下图所示。
将数据填充完毕：利用和上面的相同的方法，直接选择光滑的曲线图，而不是选择散点图：然后就画出了一条曲线：经过比较，可以知道在误差允许的范围内数学模型已经成功建立！到这里该问题的数学模型已经建立好了(采用编程所得的结果进行建模和进一步的计算)：f(t)=-0.00015t5+0.012193t4-0.366317t3+4.760591t2-22.665987t+66.025230 最后再通过编程计算任意时刻的用水率、一天的总用水量和水泵的工作功率：任意时刻的用水率的算法简单，只要给出任意时间，根据上面的式子就可以求出；一天的总用水量为上式在0-24内的积分V=∫024f(t)dt=[-0.000025t6+0.0024386t5-0.09157925t4+1.58686367t3-11.3329935t2+6 6.02523t+C]024水泵的工作功率的计算方法是用一次的总泵水量除以泵水时间水泵泵一次水泵水量为：V=2*PI*D*D*(h2-h1)/4=2*3.1415926*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<stdio.h>#include<math.h>#include<iostream.h>void yongShuiLv();void yongShuiLiang();void gongLv();void Approx(float[],float[],int,int,float[]);void main(){int i;float a[6];//多项式的待定系数个数float x[22]={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 y[22]={51.12036405,44.10152093,39.3445281,36.88021926,36.08693844,33.01166479,34.61458483,35.33279747,38.44874,70.58616637,74.79205304,70.74232197,60.78930702,62.99681227,58.99594462,55.73472518,55.68433211,59.05993154,57.55529831,59.05993154,50.95438885,44.87523183};//22个时刻所对应的用水率Approx(x,y,22,5,a);for(i=0;i<=5;i++)printf("a[%d]=%f\n",i,a[i]);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",a[5],a[4],a[3],a[2],a[1],a [0]);yongShuiLiang();gongLv();char s;while(1)//通过判断用户输入的字符是Y/N，来控制是否继续求某一时刻的用水率{cout<<"您是否还要查询某时刻的用水率（Y/N）："<<endl;cin>>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+j,x[t]);}*(c+i*(n+2)+n+1)=0;for(j=0;j<=m-1;j++)*(c+i*(n+2)+n+1)+=y[j]*power(i,x[j]);}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<=n-1;j++)if(fabs(*(c+j*(n+1)+i))>(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)=p;}for(j=i+1;j<=n-1;j++){p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));for(t=i;t<=n;t++)*(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));}}for(i=n-1;i>=0;i--){for(j=n-1;j>=i+1;j--)(*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));x[i]=*(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<<"请输入任意一个时刻，程序将就算出该时刻的用水率"<<endl;cin>>t;floatf=-0.00015*t*t*t*t*t+0.012193*t*t*t*t-0.366317*t*t*t+4.760591*t*t-22.665987*t+66.025230;cout<<t<<"时刻的用水率是："<<f<<"立方米/小时"<<endl;}void yongShuiLiang()//计算一天居民总的用水量{float v;v=-0.000025*24*24*24*24*24*24+0.0024386*24*24*24*24*24-0.09157925*24*24*24*24 +1.58686367*24*24*24-11.3329935*24*24+66.02523*24;cout<<"这个居民区的居民一天的用水量约为："<<v<<"立方米"<<endl;}void gongLv()//计算水泵工作功率{float v=3.1415926*17.4*8.7*(10.82-8.22);float t=2*(10.954-8.967)+(22.958-20.839);float p=v*9.8*(10.82-8.22)/t;cout<<"该水塔的水泵的工作功率为P="<<p<<"瓦"<<endl;}程序运行结果：四.课程设计心得这个课程设计很有趣！当然也很好的利用了我们所学的理论知识，更重要的是它将我们的理论知识和生活之极应用结合的很好！通过这次课程设计我不仅更深一层次的认识了理论知识的重要性和应用方法，并且学到了更多的分析问题解决问题的方法！同时学回了Excel 的一些操作技能，使我的知识技能得到了提高！在这次课程设计中我深刻的领悟到使用不同的算法，不同的工具，不同的技术所得到的结果会有很大的差别，比如在这次程序设计中，可以用很多方法来拟合这条函数曲线的！可以用Matlab来编程，也可以用C语言来编程，当然也可以用其他的编程语言来编程！在这些语言的选择方面要综合考虑到方便性、有效性、精确性等个性能指标！另外还可以用诸如Excel等这样的软件来更加方便的来工作，这样可以大大提高工作效率！我深刻的领悟到科技的发展，谁掌握了先进的科技水就会拥有领先地位！谁就有发言权！大家也都会服从！这样的课程设计很好！我可以学到很多，虽然不能一一的截图并列举出来，但在潜移默化中我进步了好多！我所获得的知识最多的就体现在这篇报告的每一步当中。