水塔水流量的估计

4758.96710.95411.49312.493];>>v=[14.859016.169418.714933.5029.684631.4534];>>t0=8.967:0.1:10.954;>>lglr=lglrcz(t,v,t0);>>lglrjf=0.1*trapz(lglr)>>fdxx=interp1(t,v,t0);>>fdxxjf=0.1*trapz(fdxx)>>scyt=interp1(t,v,t0,'spline');>>sancytjf=0.1*trapz(scyt)>>plot(t,v,'*',t0,lglr,'r',t0,fdxx,'g',t0,scyt,'b')>>gtext('lglr')>>gtext('fdxx')>>gtext('scyt')运行结果为:lglrjf=53.2376fdxxjf=48.9892sancytjf=52.5114图5是对第1供水段数据用三种不同方法得到的插值函数图,图中曲线lglr、fdxx和scyt分别表示用拉格朗日插值法,分段线性插值法及三次样条插值法得到的曲线。图5混合时段:首先,对于混合时段,我们分别用三种方法算出流量函数和用水量(用水高度)。下面为实现该过程的MATLAB程序:>>t=[19.49820.39920.83922.95823.41924.43325.44725.908];>>v=[24.837324.204522.752522.189324.837321.428618.872014.0533];>>t0=20.839:0.1:22.958;>>lglr=lglrcz(t,v,t0);>>lglrjf=0.1*trapz(lglr)>>fdxx=interp1(t,v,t0);>>fdxxjf=0.1*trapz(fdxx)>>scyt=interp1(t,v,t0,'spline');