资源描述：
《水流问题(数学建模).docx》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、水流问题(数学建模)估量火塔的火流量1成绩提出某住民区的平易近用自去火是由一个圆柱形的火塔供应．火塔下12.2米，曲径17．4米．火塔是由火泵依据火塔内火位下低主动减火，一样平常每一年夜火泵事情两次．如今必要懂得该住民区用火法则取火泵的事情功率．依照计划，当火塔的火位落至最低火位，约8．2米时，火泵主动启动减火；当火位降；下到一个最下火位，约10．8米时，火泵中断事情．能够思索接纳用火率（单元光阴的用火量）去反应用火法则，并经由过程距离一段光阴丈量火塔里的火位去预算用火率．表4．2是某一天的丈量纪录数据，丈量了28个时候，可是
2、因为个中有3个时候逢到火泵在背火塔供火，而无火位纪录试创建开适的数教模子，推算恣意时候的用火率、一天的总用火量以及火泵事情功率．2成绩剖析取数据处置由成绩的请求，闭键正在于断定用火率函数，即单元光阴内用火体积，记为f(t)，又称火流速率．假如可以经由过程丈量数据，发生多少个时候的用火率，也便是f(t)正在多少个面的函数值，则f(t)的盘算成绩便能够转化为插值成绩．1．假如1）火塔中火流量是光阴的一连平滑函数，取火泵事情取可有关，并疏忽火位下度对于火流速率的影响．2）火泵事情取可完整与决于火塔内火位的下度，且每一次减火的事情光阴
3、为2小时3）火塔为尺度圆柱体．思索到假如2）分离表4．2中详细数据，判断患上出4）火泵第一次供火光阴段为［8．967，10．954］，第2次供火光阴段为「20．839，22．958］．2．体积盘算火塔是一个圆柱体，体积为hDV24π=．个中D为底里曲径，h为火位下度。远似天与141592654.3=π。患上到没有同时候火塔中火的体积如表4.33．火流速率的预算火流速率应当是火塔中火的体积对于光阴的导数（微商）因为出有火的体积闭于光阴的函数抒发式，而只要一个分离的函数值表4．3，果此思索用好商朝替微商，那也是分离反应一连的经常使
4、用头脑．为普及粗度，接纳2阶好商，即iivtf2)(-?=详细天，果为一切数据被火泵两次事情宰割成3组数据，对于每一组数据的两头数据接纳中央好商，先后两个数据没有可以接纳中央好商，改用背前或者背后好商．中央好商公式)(1288121122iiiiiiittvvvvv-+-+-=?+--++背前以及背后好商公式)(243)(23412121122---+++-+-=?--+-=?iiiiiiiiiiiittvvvvttvvvv3模子创建取供解模子及盘算了局成绩已经经变化为依据流速f(t)的一个函数值表，发生函数f(t)正在全部区
5、间（2104小时）上的函数或者函数值，插值以及拟开是两种最经常使用的圆法．假如创建拟开模子，必要依据集面图的趋向，取舍得当的拟开函数情势．假如接纳插值模子，能够思索分段线性插值。3次样条插值等等．MATLAB步伐：functionwater_towerD=17.400;t1=[00.9211.8432.9493.8714.9785.9007.0067.9288.967];t2=[10.95412.03212.95413.87514.98215.90316.82617.93119.03719.95920.839];t3=[22.
6、95823.88024.98625.908];t4=[9.98110.925];t5=[22.015];stage1=[9.6779.4799.3089.1258.9828.8148.6868.5258.3888.220];stage2=[10.82010.50010.2109.9369.6539.4099.1808.9218.6628.4338.220];stage3=[10.82010.59110.35410.180];stage=[stage1stage2stage3];t=[t1,t2,t3];n1=length(t1
7、);v1=zeros(1,n1);v1=pi/4*D^2*stage1;n2=length(t2);v2=zeros(1,n2);v2=pi/4*D^2*stage2;n3=length(t3);v3=zeros(1,n3);v3=pi/4*D^2*stage3;dv1=-gradient(v1,t1);dv2=-gradient(v2,t2);dv3=-gradient(v3,t3);dv=[dv1dv2dv3];t=[t1t2t3];tt=min(t):0.001:max(t);s=interp1(t,dv,tt,'spl
8、ine');plot(t,dv,'k+',tt,s,'r');time=input('ê?è?ê±??');speed=interp1(tt,s,time);disp('′?ê±á÷?ù?a￡o');disp(speed);total=0.001*trapz(s);disp(
