實(shí)驗(yàn)十二實(shí)驗(yàn)數(shù)據(jù)的插值

實(shí)驗(yàn)十二實(shí)驗(yàn)數(shù)據(jù)的插值1一、實(shí)驗(yàn)?zāi)康膶W(xué)會MATLAB軟件中利用給定數(shù)據(jù)進(jìn)行插值運(yùn)算的方法。一、實(shí)驗(yàn)?zāi)康膶W(xué)會MATLAB軟件中利用給定數(shù)據(jù)進(jìn)行插值運(yùn)算的2二、相關(guān)知識在生產(chǎn)和科學(xué)實(shí)驗(yàn)中，自變量與因變量間的函數(shù)關(guān)系有時不能寫出解析表達(dá)式，而只能得到函數(shù)在若干點(diǎn)的函數(shù)值或?qū)?shù)值，或者表達(dá)式過于復(fù)雜需要較大的計(jì)算量而只能計(jì)算函數(shù)在若干點(diǎn)的函數(shù)值或?qū)?shù)值，當(dāng)要求知道其它點(diǎn)的函數(shù)值時，需要估計(jì)函數(shù)值在該點(diǎn)的值。為了完成這樣的任務(wù)，需要構(gòu)造一個比較簡單的函數(shù)，使函數(shù)在觀測點(diǎn)的值等于已知的值，或使函數(shù)在該點(diǎn)的導(dǎo)數(shù)值等于或者接近已知的值，尋找這樣的函數(shù)有很多方法。根據(jù)測量數(shù)據(jù)的類型有以下兩類處理觀測數(shù)據(jù)的方法。二、相關(guān)知識在生產(chǎn)和科學(xué)實(shí)驗(yàn)中，自變量與因變量間3二、相關(guān)知識（1）測量數(shù)據(jù)的數(shù)據(jù)量較小并且數(shù)據(jù)值是準(zhǔn)確的，或者基本沒有誤差，這時我們一般用插值的方法來解決問題。（2）測量數(shù)據(jù)的數(shù)據(jù)量較大或者測量值與真實(shí)值有誤差，這時一般用曲線擬合的方法來解決問題。為了完成這樣的任務(wù)，需要構(gòu)造一個比較簡單的函數(shù)，使函數(shù)在觀測點(diǎn)的值等于已知的值，或使函數(shù)在該點(diǎn)的導(dǎo)數(shù)值等于或者接近已知的值，尋找這樣的函數(shù)有很多方法。根據(jù)測量數(shù)據(jù)的類型有以下兩類處理觀測數(shù)據(jù)的方法。

二、相關(guān)知識（1）測量數(shù)據(jù)的數(shù)據(jù)量較小并且數(shù)據(jù)值是準(zhǔn)確的，或4（1）測量數(shù)據(jù)的數(shù)據(jù)量較小并且數(shù)據(jù)值是準(zhǔn)確的，或者基本沒有誤差，這時我們一般用插值的方法來解決問題。（2）測量數(shù)據(jù)的數(shù)據(jù)量較大或者測量值與真實(shí)值有誤差，這時一般用曲線擬合的方法來解決問題。在MATLAB中，無論是插值還是擬合，都有相應(yīng)的命令來處理。本實(shí)驗(yàn)討論插值。（1）測量數(shù)據(jù)的數(shù)據(jù)量較小并且數(shù)據(jù)值是準(zhǔn)確的，或者基本沒有誤5一維插值已知離散點(diǎn)上的數(shù)據(jù)集，即已知在點(diǎn)集上的函數(shù)值，構(gòu)造一個解析函數(shù)（其圖形為一曲線）通過這些點(diǎn)，并能夠求出這些點(diǎn)之間的值，這一過程稱為一維插值。完成這一過程可以有多種方法，我們現(xiàn)在在MATLAB中，無論是插值還是擬合，都有相應(yīng)的命令來處理。本實(shí)驗(yàn)討論插值。

一維插值6一維插值已知離散點(diǎn)上的數(shù)據(jù)集，即已知在點(diǎn)集上的函數(shù)值，構(gòu)造一個解析函數(shù)（其圖形為一曲線）通過這些點(diǎn)，并能夠求出這些點(diǎn)之間的值，這一過程稱為一維插值。完成這一過程可以有多種方法，我們現(xiàn)在利用MATLAB提供的函數(shù)interp1，這個函數(shù)的調(diào)用格式為：yi=interp1(X,Y,xi,method)

該命令用指定的算法找出一個一元函數(shù)，然后以給出處的值。一維插值7xi可以是一個標(biāo)量，也可以是一個向量，是向量時，必須單調(diào)，method可以下列方法之一：‘nearest’：最近鄰點(diǎn)插值，直接完成計(jì)算；‘spline’：三次樣條函數(shù)插值；‘linear’：線性插值（缺省方式），直接完成計(jì)算利用MATLAB提供的函數(shù)interp1，這個函數(shù)的調(diào)用格式為：yi=interp1(X,Y,xi,method)

該命令用指定的算法找出一個一元函數(shù)，然后以給出處的值。

xi可以是一個標(biāo)量，也可以是一個向量，是向量時，必須單調(diào)，8xi可以是一個標(biāo)量，也可以是一個向量，是向量時，必須單調(diào)，method可以下列方法之一：‘nearest’：最近鄰點(diǎn)插值，直接完成計(jì)算；‘spline’：三次樣條函數(shù)插值；‘linear’：線性插值（缺省方式），直接完成計(jì)算

‘cubic’：三次函數(shù)插值；對于[min{xi},max{xi}]外的值，MATLAB使用外推的方法計(jì)算數(shù)值。xi可以是一個標(biāo)量，也可以是一個向量，是向量時，必須單調(diào)，9例1：已知某產(chǎn)品從1900年到2010年每隔10年的產(chǎn)量為：75.995,91.972,105.711,123.203,131.699,150.697,179.323,203.212,226.505,249.633,256.344,267.893，計(jì)算出1995年的產(chǎn)量，用三次樣條插值的方法，畫出每隔一年的插值曲線圖形，同時將原始的數(shù)據(jù)畫在同一圖上。

‘cubic’：三次函數(shù)插值；對于[min{xi},max{xi}]外的值，MATLAB使用外推的方法計(jì)算數(shù)值。例1：已知某產(chǎn)品從1900年到2010年每隔10年的產(chǎn)量為：10例1：已知某產(chǎn)品從1900年到2010年每隔10年的產(chǎn)量為：75.995,91.972,105.711,123.203,131.699,150.697,179.323,203.212,226.505,249.633,256.344,267.893，計(jì)算出1995年的產(chǎn)量，用三次樣條插值的方法，畫出每隔一年的插值曲線圖形，同時將原始的數(shù)據(jù)畫在同一圖上。解：程序如下year=1900:10:2010;product=[75.995,91.972,105.711,123.203,131.699,150.697,179.323,203.212,226.505,249.633,256.344,267.893]p1995=interp1(year,product,1995,'spline')例1：已知某產(chǎn)品從1900年到2010年每隔10年的產(chǎn)量為：11x=1900:2010;y=interp1(year,product,x,'spline');plot(year,product,'o',x,y);計(jì)算結(jié)果為：p1995=253.2278如果用線性插值，則程序的后四行改為：解：程序如下year=1900:10:2010;product=[75.995,91.972,105.711,123.203,131.699,150.697,179.323,203.212,226.505,249.633,256.344,267.893]p1995=interp1(year,product,1995,'spline')x=1900:2010;12x=1900:2010;y=interp1(year,product,x,'spline');plot(year,product,'o',x,y);計(jì)算結(jié)果為：p1995=253.2278如果用線性插值，則程序的后四行改為：p1995=interp1(year,product,1995,'linear')x=1900:2010;y=interp1(year,product,x,'linear');plot(year,product,'o',x,y);計(jì)算結(jié)果為：p1995=252.9885x=1900:2010;13這兩種計(jì)算方法得到的數(shù)據(jù)有微小的差異，這種差異我們從兩個圖形上也能夠看到，主要表現(xiàn)在節(jié)點(diǎn)（那些繪制成圓點(diǎn)的點(diǎn)）的附近。前者是光滑的，后者有角點(diǎn)出現(xiàn)。p1995=interp1(year,product,1995,'linear')x=1900:2010;y=interp1(year,product,x,'linear');plot(year,product,'o',x,y);計(jì)算結(jié)果為：p1995=252.9885這兩種計(jì)算方法得到的數(shù)據(jù)有微小的差異，這種差異我們從兩個圖形14這兩種計(jì)算方法得到的數(shù)據(jù)有微小的差異，這種差異我們從兩個圖形上也能夠看到，主要表現(xiàn)在節(jié)點(diǎn)（那些繪制成圓點(diǎn)的點(diǎn)）的附近。前者是光滑的，后者有角點(diǎn)出現(xiàn)。這兩種計(jì)算方法得到的數(shù)據(jù)有微小的差異，這種差異我們從兩個圖形15二維插值已知離散點(diǎn)上的數(shù)據(jù)集即已知在點(diǎn)集上的函數(shù)值，構(gòu)造一個解析函數(shù)（其圖形為一曲面）通過這些點(diǎn)，并能夠求出這些已知點(diǎn)以外的點(diǎn)的函數(shù)值，這一過程稱為二維插值。MATLAB命令：Zi=interp2(X,Y,Z,Xi,Yi,method)該命令用指定的算法找出一個二元函數(shù)，然后以給出處的值。返回?cái)?shù)據(jù)矩陣Zi，Xi，Yi是向量，且必須單調(diào)，Zi和meshgrid(Xi,Yi)是同類型的。二維插值16method可以下列方法之一：‘nearest’：最近鄰點(diǎn)插值，直接完成計(jì)算；‘spline’：三次樣條函數(shù)插值；‘linear’：線性插值（缺省方式），直接完成計(jì)算‘cubic’：三次函數(shù)插值；MATLAB命令：Zi=interp2(X,Y,Z,Xi,Yi,method)該命令用指定的算法找出一個二元函數(shù)，然后以給出處的值。返回?cái)?shù)據(jù)矩陣Zi，Xi，Yi是向量，且必須單調(diào)，Zi和meshgrid(Xi,Yi)是同類型的。method可以下列方法之一：17method可以下列方法之一：‘nearest’：最近鄰點(diǎn)插值，直接完成計(jì)算；‘spline’：三次樣條函數(shù)插值；‘linear’：線性插值（缺省方式），直接完成計(jì)算‘cubic’：三次函數(shù)插值；

method可以下列方法之一：18例2：已知1950年到1990年間每隔10年，服務(wù)年限從10年到30年每隔10年的勞動報(bào)酬表如下：某企業(yè)工作人員的月平均工資（元）試計(jì)算1975年時，15年工齡的工作人員平均工資。服務(wù)年限年份1020301950150.697169.592187.6521960179.323195.072250.2871970203.212239.092322.7671980226.505273.706426.7301990249.633370.281598.243例2：已知1950年到1990年間每隔10年，服務(wù)年限從1019解：程序如下：years=1950:10:1990;service=10:10:30;wage=[150.697169.592187.652179.323195.072250.287試計(jì)算1975年時，15年工齡的工作人員平均工資。1960179.323195.072250.2871970203.212239.092322.7671980226.505273.706426.7301990249.633370.281598.243解：程序如下：1960179.323195.072250.220解：程序如下：years=1950:10:1990;service=10:10:30;wage=[150.697169.592187.652179.323195.072250.287203.212239.092322.767226.505273.706426.730249.633370.281598.243]w=interp2(service,years,wage,15,1975)計(jì)算結(jié)果為：235.6288解：程序如下：21例3：設(shè)有數(shù)據(jù)x=1,2,3,4,5,6，y=1,2,3,4，在由x，y構(gòu)成的網(wǎng)格上，數(shù)據(jù)為：12,10,11,11,13,1516,22,28,35,27,2018,21,26,32,28,25

203.212239.092322.767226.505273.706426.730249.633370.281598.243]w=interp2(service,years,wage,15,1975)計(jì)算結(jié)果為：235.6288例3：設(shè)有數(shù)據(jù)x=1,2,3,4,5,6，y=1,2,3,422例3：設(shè)有數(shù)據(jù)x=1,2,3,4,5,6，y=1,2,3,4，在由x，y構(gòu)成的網(wǎng)格上，數(shù)據(jù)為：12,10,11,11,13,1516,22,28,35,27,2018,21,26,32,28,2520,25,30,33,32,20畫出原始網(wǎng)格圖和將網(wǎng)格細(xì)化為間隔為0.1后的插值網(wǎng)格圖。解：程序?yàn)椋簒=1:6;y=1:4;例3：設(shè)有數(shù)據(jù)x=1,2,3,4,5,6，y=1,2,3,423t=[12,10,11,11,13,1516,22,28,35,27,20

18,21,26,32,28,25;20,25,30,33,32,20]subplot(1,2,1)mesh(x,y,t)20,25,30,33,32,20畫出原始網(wǎng)格圖和將網(wǎng)格細(xì)化為間隔為0.1后的插值網(wǎng)格圖。解：程序?yàn)椋簒=1:6;y=1:4;t=[12,10,11,11,13,1524t=[12,10,11,11,13,1516,22,28,35,27,20

18,21,26,32,28,25;20,25,30,33,32,20]subplot(1,2,1)mesh(x,y,t)x1=1:0.1:6;y1=1:0.1:4;[x2,y2]=meshgrid(x1,y1);t1=interp2(x,y,t,x2,y2,'cubic');subplot(1,2,2)mesh(x1,y1,t1);t=[12,10,11,11,13,1525結(jié)果如圖。左圖是給定的網(wǎng)格處的數(shù)據(jù)，右圖是插值后的數(shù)據(jù)。結(jié)果如圖。左圖是給定的網(wǎng)格處的數(shù)據(jù)，右圖是插值后的數(shù)據(jù)。26三、實(shí)驗(yàn)內(nèi)容1．已知x=[0.1,0.8,1.3,1.9,2.5,3.1]，y=[1.2,1.6,2.7,2.0,1.3,0.5]，利用其中的部分?jǐn)?shù)據(jù)，分別用線性插值和3次插值，求x=2.0處的值。2．已知二元函數(shù)在點(diǎn)集上的值為，三、實(shí)驗(yàn)內(nèi)容1．已知x=[0.1,0.8,1.3,1.9,227其中，左上角位置表示，右下角位置表示畫出原始網(wǎng)格圖和將網(wǎng)格細(xì)化為間隔為0.1后的插值網(wǎng)格圖。*3．學(xué)習(xí)函數(shù)interp3(X,Y,Z,V,X1,Y1,Z1,method),對MATLAB提供的flow數(shù)據(jù)實(shí)現(xiàn)三維插值。4．完成實(shí)驗(yàn)報(bào)告。其中，左上角位置表示，右下角位置表示畫28實(shí)驗(yàn)十二實(shí)驗(yàn)數(shù)據(jù)的插值

實(shí)驗(yàn)十二實(shí)驗(yàn)數(shù)據(jù)的插值29一、實(shí)驗(yàn)?zāi)康膶W(xué)會MATLAB軟件中利用給定數(shù)據(jù)進(jìn)行插值運(yùn)算的方法。一、實(shí)驗(yàn)?zāi)康膶W(xué)會MATLAB軟件中利用給定數(shù)據(jù)進(jìn)行插值運(yùn)算的30二、相關(guān)知識在生產(chǎn)和科學(xué)實(shí)驗(yàn)中，自變量與因變量間的函數(shù)關(guān)系有時不能寫出解析表達(dá)式，而只能得到函數(shù)在若干點(diǎn)的函數(shù)值或?qū)?shù)值，或者表達(dá)式過于復(fù)雜需要較大的計(jì)算量而只能計(jì)算函數(shù)在若干點(diǎn)的函數(shù)值或?qū)?shù)值，當(dāng)要求知道其它點(diǎn)的函數(shù)值時，需要估計(jì)函數(shù)值在該點(diǎn)的值。為了完成這樣的任務(wù)，需要構(gòu)造一個比較簡單的函數(shù)，使函數(shù)在觀測點(diǎn)的值等于已知的值，或使函數(shù)在該點(diǎn)的導(dǎo)數(shù)值等于或者接近已知的值，尋找這樣的函數(shù)有很多方法。根據(jù)測量數(shù)據(jù)的類型有以下兩類處理觀測數(shù)據(jù)的方法。二、相關(guān)知識在生產(chǎn)和科學(xué)實(shí)驗(yàn)中，自變量與因變量間31二、相關(guān)知識（1）測量數(shù)據(jù)的數(shù)據(jù)量較小并且數(shù)據(jù)值是準(zhǔn)確的，或者基本沒有誤差，這時我們一般用插值的方法來解決問題。（2）測量數(shù)據(jù)的數(shù)據(jù)量較大或者測量值與真實(shí)值有誤差，這時一般用曲線擬合的方法來解決問題。為了完成這樣的任務(wù)，需要構(gòu)造一個比較簡單的函數(shù)，使函數(shù)在觀測點(diǎn)的值等于已知的值，或使函數(shù)在該點(diǎn)的導(dǎo)數(shù)值等于或者接近已知的值，尋找這樣的函數(shù)有很多方法。根據(jù)測量數(shù)據(jù)的類型有以下兩類處理觀測數(shù)據(jù)的方法。

二、相關(guān)知識（1）測量數(shù)據(jù)的數(shù)據(jù)量較小并且數(shù)據(jù)值是準(zhǔn)確的，或32（1）測量數(shù)據(jù)的數(shù)據(jù)量較小并且數(shù)據(jù)值是準(zhǔn)確的，或者基本沒有誤差，這時我們一般用插值的方法來解決問題。（2）測量數(shù)據(jù)的數(shù)據(jù)量較大或者測量值與真實(shí)值有誤差，這時一般用曲線擬合的方法來解決問題。在MATLAB中，無論是插值還是擬合，都有相應(yīng)的命令來處理。本實(shí)驗(yàn)討論插值。（1）測量數(shù)據(jù)的數(shù)據(jù)量較小并且數(shù)據(jù)值是準(zhǔn)確的，或者基本沒有誤33一維插值已知離散點(diǎn)上的數(shù)據(jù)集，即已知在點(diǎn)集上的函數(shù)值，構(gòu)造一個解析函數(shù)（其圖形為一曲線）通過這些點(diǎn)，并能夠求出這些點(diǎn)之間的值，這一過程稱為一維插值。完成這一過程可以有多種方法，我們現(xiàn)在在MATLAB中，無論是插值還是擬合，都有相應(yīng)的命令來處理。本實(shí)驗(yàn)討論插值。

一維插值34一維插值已知離散點(diǎn)上的數(shù)據(jù)集，即已知在點(diǎn)集上的函數(shù)值，構(gòu)造一個解析函數(shù)（其圖形為一曲線）通過這些點(diǎn)，并能夠求出這些點(diǎn)之間的值，這一過程稱為一維插值。完成這一過程可以有多種方法，我們現(xiàn)在利用MATLAB提供的函數(shù)interp1，這個函數(shù)的調(diào)用格式為：yi=interp1(X,Y,xi,method)

該命令用指定的算法找出一個一元函數(shù)，然后以給出處的值。一維插值35xi可以是一個標(biāo)量，也可以是一個向量，是向量時，必須單調(diào)，method可以下列方法之一：‘nearest’：最近鄰點(diǎn)插值，直接完成計(jì)算；‘spline’：三次樣條函數(shù)插值；‘linear’：線性插值（缺省方式），直接完成計(jì)算利用MATLAB提供的函數(shù)interp1，這個函數(shù)的調(diào)用格式為：yi=interp1(X,Y,xi,method)

該命令用指定的算法找出一個一元函數(shù)，然后以給出處的值。

xi可以是一個標(biāo)量，也可以是一個向量，是向量時，必須單調(diào)，36xi可以是一個標(biāo)量，也可以是一個向量，是向量時，必須單調(diào)，method可以下列方法之一：‘nearest’：最近鄰點(diǎn)插值，直接完成計(jì)算；‘spline’：三次樣條函數(shù)插值；‘linear’：線性插值（缺省方式），直接完成計(jì)算

‘cubic’：三次函數(shù)插值；對于[min{xi},max{xi}]外的值，MATLAB使用外推的方法計(jì)算數(shù)值。xi可以是一個標(biāo)量，也可以是一個向量，是向量時，必須單調(diào)，37例1：已知某產(chǎn)品從1900年到2010年每隔10年的產(chǎn)量為：75.995,91.972,105.711,123.203,131.699,150.697,179.323,203.212,226.505,249.633,256.344,267.893，計(jì)算出1995年的產(chǎn)量，用三次樣條插值的方法，畫出每隔一年的插值曲線圖形，同時將原始的數(shù)據(jù)畫在同一圖上。

‘cubic’：三次函數(shù)插值；對于[min{xi},max{xi}]外的值，MATLAB使用外推的方法計(jì)算數(shù)值。例1：已知某產(chǎn)品從1900年到2010年每隔10年的產(chǎn)量為：38例1：已知某產(chǎn)品從1900年到2010年每隔10年的產(chǎn)量為：75.995,91.972,105.711,123.203,131.699,150.697,179.323,203.212,226.505,249.633,256.344,267.893，計(jì)算出1995年的產(chǎn)量，用三次樣條插值的方法，畫出每隔一年的插值曲線圖形，同時將原始的數(shù)據(jù)畫在同一圖上。解：程序如下year=1900:10:2010;product=[75.995,91.972,105.711,123.203,131.699,150.697,179.323,203.212,226.505,249.633,256.344,267.893]p1995=interp1(year,product,1995,'spline')例1：已知某產(chǎn)品從1900年到2010年每隔10年的產(chǎn)量為：39x=1900:2010;y=interp1(year,product,x,'spline');plot(year,product,'o',x,y);計(jì)算結(jié)果為：p1995=253.2278如果用線性插值，則程序的后四行改為：解：程序如下year=1900:10:2010;product=[75.995,91.972,105.711,123.203,131.699,150.697,179.323,203.212,226.505,249.633,256.344,267.893]p1995=interp1(year,product,1995,'spline')x=1900:2010;40x=1900:2010;y=interp1(year,product,x,'spline');plot(year,product,'o',x,y);計(jì)算結(jié)果為：p1995=253.2278如果用線性插值，則程序的后四行改為：p1995=interp1(year,product,1995,'linear')x=1900:2010;y=interp1(year,product,x,'linear');plot(year,product,'o',x,y);計(jì)算結(jié)果為：p1995=252.9885x=1900:2010;41這兩種計(jì)算方法得到的數(shù)據(jù)有微小的差異，這種差異我們從兩個圖形上也能夠看到，主要表現(xiàn)在節(jié)點(diǎn)（那些繪制成圓點(diǎn)的點(diǎn)）的附近。前者是光滑的，后者有角點(diǎn)出現(xiàn)。p1995=interp1(year,product,1995,'linear')x=1900:2010;y=interp1(year,product,x,'linear');plot(year,product,'o',x,y);計(jì)算結(jié)果為：p1995=252.9885這兩種計(jì)算方法得到的數(shù)據(jù)有微小的差異，這種差異我們從兩個圖形42這兩種計(jì)算方法得到的數(shù)據(jù)有微小的差異，這種差異我們從兩個圖形上也能夠看到，主要表現(xiàn)在節(jié)點(diǎn)（那些繪制成圓點(diǎn)的點(diǎn)）的附近。前者是光滑的，后者有角點(diǎn)出現(xiàn)。這兩種計(jì)算方法得到的數(shù)據(jù)有微小的差異，這種差異我們從兩個圖形43二維插值已知離散點(diǎn)上的數(shù)據(jù)集即已知在點(diǎn)集上的函數(shù)值，構(gòu)造一個解析函數(shù)（其圖形為一曲面）通過這些點(diǎn)，并能夠求出這些已知點(diǎn)以外的點(diǎn)的函數(shù)值，這一過程稱為二維插值。MATLAB命令：Zi=interp2(X,Y,Z,Xi,Yi,method)該命令用指定的算法找出一個二元函數(shù)，然后以給出處的值。返回?cái)?shù)據(jù)矩陣Zi，Xi，Yi是向量，且必須單調(diào)，Zi和meshgrid(Xi,Yi)是同類型的。二維插值44method可以下列方法之一：‘nearest’：最近鄰點(diǎn)插值，直接完成計(jì)算；‘spline’：三次樣條函數(shù)插值；‘linear’：線性插值（缺省方式），直接完成計(jì)算‘cubic’：三次函數(shù)插值；MATLAB命令：Zi=interp2(X,Y,Z,Xi,Yi,method)該命令用指定的算法找出一個二元函數(shù)，然后以給出處的值。返回?cái)?shù)據(jù)矩陣Zi，Xi，Yi是向量，且必須單調(diào)，Zi和meshgrid(Xi,Yi)是同類型的。method可以下列方法之一：45method可以下列方法之一：‘nearest’：最近鄰點(diǎn)插值，直接完成計(jì)算；‘spline’：三次樣條函數(shù)插值；‘linear’：線性插值（缺省方式），直接完成計(jì)算‘cubic’：三次函數(shù)插值；

method可以下列方法之一：46例2：已知1950年到1990年間每隔10年，服務(wù)年限從10年到30年每隔10年的勞動報(bào)酬表如下：某企業(yè)工作人員的月平均工資（元）試計(jì)算1975年時，15年工齡的工作人員平均工資。服務(wù)年限年份1020301950150.697169.592187.6521960179.323195.072250.2871970203.212239.092322.7671980226.505273.706426.7301990249.633370.281598.243例2：已知1950年到1990年間每隔10年，服務(wù)年限從1047解：程序如下：years=1950:10:1990;service=10:10:30;wage=[150.697169.592187.652179.323195.072250.287試計(jì)算1975年時，15年工齡的工作人員平均工資。1960179.323195.072250.2871970203.212239.092322.7671980226.505273.706426.7301990249.633370.281598.243解：程序如下：1960179.323195.072250.248解：程序如下：years=1950:10:1990;service=10:10:30;wage=[150.697169.592187.652179.323195.072250.287203.212239.092322.767226.505273.706426.730249.633370.281598.243]w=interp2(service,years,wage,15,1975)計(jì)算結(jié)果為：235.6288解：程序如下：49例3：設(shè)有數(shù)據(jù)x=1,2,3,4,5,6，y=1,2,3,4，在由x，y構(gòu)成的網(wǎng)格上，數(shù)據(jù)為：12,10,11,11,13,1516,22,28,35,27,2018,21,26,32,28,25