樣條插值實(shí)驗(yàn)與matlab仿真

1、樣條插值實(shí)驗(yàn)樣條插值實(shí)驗(yàn)Matlab的一維插值函數(shù)為interp1()，調(diào)用格式為：yy=interp1(xyxx方法)其中x=[x1x2…xn]’y=[y1y2…yn]’兩個向量分別為給定的一組自變量和函數(shù)值，用來表示已知樣本點(diǎn)數(shù)據(jù).xx為待求插值點(diǎn)處橫坐標(biāo)，yy返回的對應(yīng)縱坐標(biāo)。插值方法可以選用默認(rèn)的’linear’(線性插值)，’nearest’(最近鄰等值方式)，’cubic’（三次Hermite插值），’spline’(三次樣

2、條插值)。一般可采用樣條插值。實(shí)驗(yàn)1作函數(shù)在[01]取間隔為0.1的點(diǎn)圖，用插值進(jìn)行實(shí)驗(yàn)24(37)..sin(2)xyxxex????Matlab程序chazhi1.m：x=0:0.1:1y=(x.^23x7).exp(4x).sin(2x)%產(chǎn)生原始數(shù)據(jù)subplot(121)plot(xyxyro)%作圖xx=0:0.02:1%待求插值點(diǎn)yy=interp1(xyxxspline)%此處可用nearestcubicspline分別

3、試驗(yàn)subplot(122)plot(xyroxxyyb)%作圖圖1(1)原始點(diǎn)圖圖1(2)樣條插值圖體積計算公式為4.2hdv??表2不同時刻水體積表時間(小時)水體積(加侖)時間(小時)水體積(加侖)時間(小時)水體積(加侖)0(1)60612510.9542(2)67771520.83925148720.921159371612.032865767022.9581(3)6777151.843158302612.9544639534

4、23.88006633972.949757157113.875862235224.98696485063.871456259914.982260459825.90836376254.978155209915.90395893255.900054408116.82615750087.006453396317.93175587817.928652537219.03755425548.967851487219.9594528236其中(1)表

5、示第一段開始，(2)表示第二段開始，(3)表示第三段開始3.計算各時刻點(diǎn)的水流量(加侖小時)水流量公式為：dttdvtf)()(?以上25個時刻處的水流量采用差分的方法得到，共分三段分別處理。差分公式為：(1)對每段前兩點(diǎn)采用向前差分公式)(243)(121iiiiiittVVVtf????????(2)對每段最后兩點(diǎn)采用向后差分公式12134()2()iiiiiiVVVfttt???????(3)對每段中間點(diǎn)采用中心差分公式)(128

6、8)(12112iiiiiiittVVVVtf???????????得到各點(diǎn)水流量表表3不同時刻水流量表時間(小時)水流量(加侖小時)時間(小時)水流量(加侖小時)時間(小時)水流量(加侖小時)0(1)1440410.9542(2)1946920.8392146480.92111118212.03282019522.9581(3)152201.84311006312.95441894123.8800152632.94971101213.