中北大學(xué)關(guān)于插值方法的報告

1、關(guān)于數(shù)值分析中插值方法的報告學(xué) 院：機電工程學(xué)院姓 名： 學(xué) 號：2015年 12 月 22 日1應(yīng)用背景在離散數(shù)據(jù)的基礎(chǔ)上補插連續(xù)函數(shù)，使得這條連續(xù)曲線通過全部給定的離散數(shù)據(jù)點。插值是離散函數(shù)逼近的重要方法，利用它可通過函數(shù)在有限個點處的取值狀況，估算出函數(shù)在其他點處的近似值。 早在6世紀(jì)，中國的劉焯已將等距二次插值用于天文計算。17世紀(jì)之后，牛頓、拉格朗日分別討論了等距和非等距的一般插值公式。在近代，插值法仍然是數(shù)據(jù)處理和編制函數(shù)表的常用工具，又是數(shù)值積分、數(shù)值微分、非線性方程求根和微分方程數(shù)值解法的重要基礎(chǔ)，許多求解計算公式都是以插值為基礎(chǔ)導(dǎo)出的。 插值問題的提法是：假定區(qū)間a，b上的實

2、值函數(shù)f（x）在該區(qū)間上 n1個互不相同點x0，x1xn 處的值是f（x0），f（xn），要求估算f（x）在a，b中某點的值。其做法是：在事先選定的一個由簡單函數(shù)構(gòu)成的有n1個參數(shù)C0，C1，Cn的函數(shù)類(C0，C1，Cn)中求出滿足條件P（xi）f（xi）（i0，1， n）的函數(shù)P(x)，并以P(x)作為f(x)的估值。此處f（x）稱為被插值函數(shù)，x0，x1，xn稱為插值結(jié)(節(jié))點，(C0，C1，Cn)稱為插值函數(shù)類，上面等式稱為插值條件，（C0，Cn）中滿足上式的函數(shù)稱為插值函數(shù)，R（x） f（x）P（x）稱為插值余項。求解這類問題，它有很多種插值法，其中以拉格朗日(Lagrange)插值

3、和牛頓(Newton)插值為代表的多項式插值最有特點，常用的插值還有Hermit插值，分段插值和樣條插值。數(shù)值分析是高等學(xué)校理工科一門重要的基礎(chǔ)課程，主要研究數(shù)學(xué)方法的數(shù) 值求解。數(shù)值分析是各種計算性科學(xué)的聯(lián)系紐帶和共性基礎(chǔ)，是一門兼有基礎(chǔ)性、 應(yīng)用性和邊緣性的交叉學(xué)科，數(shù)值分析中插值法包括拉格朗日插值法、牛頓插值法、 埃爾米特插值法等。本文主要介紹了各種插值方法的計算分析和推導(dǎo)，通過簡單的 例題進行算法分析并編程得出計算結(jié)果。2多種插值方法的分析比較數(shù)值分析插值法是一種古老的數(shù)學(xué)方法，它來自生產(chǎn)實踐。利用計算機解決工程問題與常規(guī)手工計算的差異就在于它特別的計算方法.電 機設(shè)計中常常需要通過查

4、曲線、表格或通過作圖來確定某一參量，如查磁化曲線、査異步電動機飽和系數(shù)曲線等.手工設(shè)計時，設(shè)計者是通過尋找坐標(biāo)的方法來實現(xiàn).用 計算機來完成上述工作時，采用數(shù)值插值法來完成。因此學(xué)好數(shù)值分析的插值法很重要。2. 1插值方法的定義插值法又稱“內(nèi)插法”，是利用函數(shù)f (x)在某區(qū)間中插入若干點的函數(shù)值，作出適當(dāng)?shù)奶囟ê瘮?shù)，在這些點上取已知值，在區(qū)間的其他點上用這特定函數(shù)的值作為函數(shù)f (x)的近似值，這種方法稱為插值法。如果這特定函數(shù)是多項式，就稱它為插值多項式。 3常用的幾種插值方法3. 1 Lagrange 插值3.1.1基本原理構(gòu)造n次多項式Pn (x)= yk lk (x)=y0l0 (x

5、)+y1l1 (x)+ynln (x)，這是不超過n次的多項式，其中基函數(shù)lk(x)=顯然lk (x)滿足lk (xi)=此時 Pn(x)f(x),誤差Rn(x)=f(x)-Pn(x)= 其中(a,b)且依賴于x，=（x-x0）(x-x1)(x-xn)很顯然，當(dāng)n=1、插值節(jié)點只有兩個xk,xk+1時 P1(x)=yklk(x)+yk+1lk+1(x)其中基函數(shù)lk(x)= lk+1(x)= 3.1.2優(yōu)缺點可對插值函數(shù)選擇多種不同的函數(shù)類型，由于代數(shù)多項式具有簡單和一些良好的特性，故常選用代數(shù)多項式作為插值函數(shù)。利用插值基函數(shù)很容易得到拉格朗日插值多項式，公式結(jié)構(gòu)緊湊，在理論分析中甚為方便，

6、但當(dāng)插值節(jié)點增減時全部插值基函數(shù)Lk(x)(k=0,1,n)均要隨之變化，整個公式也將發(fā)生變化，這在實際計算中是很不方便的，為了克服這一缺點，提出了牛頓插值可以克服這一缺點。3. 2 Newton 插值3.2.1基本原理構(gòu)造n次多項式Nn(x)=f(x0)+f(x0,x1)(x-x0)+f(x0,x1,x2)(x-x0)(x-x1)+f(x0,x1,x2,xn)(x-x0)(x-x1)(x-xn)稱為牛頓插值多項式，其中 (二個節(jié)點，一階差商) (三個節(jié)點，二階差商) (n+1個節(jié)點，n階差商)注意：由于插值多項式的唯一性，有時為了避免拉格朗日余項Rn(x)中n+1階導(dǎo)數(shù)的運算，用牛頓插值公式

7、Rn (x)=f(x)-Nn(x)=f(x,x0,xn)n+1(x),其中n+1(x)=(x-x0)(x-x1)(x-xn)3.2.2優(yōu)缺點牛頓插值法具有承襲性和易變性的特點，當(dāng)增加一個節(jié)點時，只要再增加一項就可以了即 而拉格朗日插值若要增加一個節(jié)點時全部基函數(shù)都需要重新算過。牛頓插值法既適合于用來計算函數(shù)值，也適合于做理論推導(dǎo)，比如說可用來推導(dǎo)微分方程的數(shù)值求解公式。 4數(shù)值實驗4.1實驗?zāi)康?.通過進行不同類型的插值，比較各種插值的效果，明確各種插值的優(yōu)越性；2.通過比較不同次數(shù)的多項式擬合效果，了解多項式擬合的原理；3.利用matlab編程，學(xué)會matlab命令；4.掌握拉格朗日插值法；

8、5.掌握多項式擬合的特點和方法。4.2實驗題目4.2.1插值法實驗將區(qū)間-5,510等分，對下列函數(shù)分別計算插值節(jié)點的值，進行不同類型的插值，作出插值函數(shù)的圖形并與的圖形進行比較： (1) 做拉格朗日插值；(2) 做分段線性插值；(3) 做三次樣條插值.4.2.2擬合實驗給定數(shù)據(jù)點如下表所示：-1.5-1.0-0.50.00.51.01.5-4.45-0.450.550.05-0.440.544.55分別對上述數(shù)據(jù)作三次多項式和五次多項式擬合，并求平方誤差，作出離散函數(shù)和擬合函數(shù)的圖形。4.2.3實驗原理1.、插值法實驗 2、 擬合實驗4.3實驗內(nèi)容4.3.1插值法實驗A：實驗步驟：打開mat

9、lab軟件，新建一個名為chazhi.m的M文件，編寫程序（見1.2實驗程序），運行程序，記錄結(jié)果。B：實驗程序：x=-5:1:5;xx=-5:0.05:5;y1=1./(1+x.2);L=malagr(x,y1,xx);L1=interp1(x,y1,x,'linear'); S=maspline(x,y1,0.0148,-0.0148,xx);hold on;plot(x,y1,'b*');plot(xx,L,'r');plot(x,L1,'g');plot(xx,S,'k'); figurex=-5:1:5;

10、xx=-5:0.05:5;y2=atan(x);L=malagr(x,y2,xx);L1=interp1(x,y2,x,'linear'); S=maspline(x,y2,0.0385,0.0385,xx);hold on;plot(x,y2,'b*');plot(xx,L,'r');plot(x,L1,'g');plot(xx,S,'k'); figurex=-5:1:5;xx=-5:0.05:5;y3=x.2./(1+x.4);L=malagr(x,y3,xx);L1=interp1(x,y3,x,'

11、linear'); S=maspline(x,y3,0.0159,-0.0159,xx);hold on;plot(x,y3,'b*');plot(xx,L,'r');plot(x,L1,'g');plot(xx,S,'k');C：實驗設(shè)備： matlab軟件。4.3.2擬合實驗A：實驗步驟：新建一個名為nihe.m的M文件，編寫程序（見2.2實驗源程序），運行程序，記錄結(jié)果。B實驗程序：x=-1.5 -1.0 -0.5 0 0.5 1.0 1.5;y=-4.45 -0.45 0.55 0.05 -0.44 0.54 4.

12、55;a1=mafit(x,y,3)x1=-1.5:0.05:1.5;y1=a1(4)+a1(3)*x1+a1(2)*x1.2+a1(1)*x1.3;hold onplot(x,y,'b*');plot(x1,y1,'r');p1=polyval(a1,x);s1=norm(y-p1) figurea2=mafit(x,y,5)x2=-1.5:0.05:1.5;y2=a2(6)+a2(5)*x2+a2(4)*x2.2+a2(3)*x2.3+a2(2)*x2.4+a2(1)*x2.5;hold onplot(x,y,'b*');plot(x2,y2

13、,'r');p2=polyval(a2,x);s2=norm(y-p2)C：實驗設(shè)備： matlab軟件。4.4實驗結(jié)果4.4.1插值法實驗（1）（2）（3）4.4.2擬合實驗（1）平方誤差：s1 =0.0136輸入程序得到：a1 =2.0000 -0.0014 -1.5007 0.0514s1 = 0.0136（2）平方誤差：s2 = 0.0069輸入程序得到：a2 = 0.0120 0.0048 1.9650 -0.0130 -1.4820 0.0545s2 = 0.006>>4.5實驗結(jié)果分析4.5.1插值法實驗結(jié)果分析：（1） 由插值結(jié)果曲線圖可見，拉格朗日插值在節(jié)點附近誤差很小，但在兩端有振蕩現(xiàn)象；分段線性插值具有良好的收斂性，但在節(jié)點處不光滑；而三次樣條插值在直觀上與原函數(shù)曲線吻合得最好；（2）分析可知，均勻插值時（拉格朗日插值），