數(shù)值分析實(shí)驗(yàn)報(bào)告(一)(完整).doc

數(shù)值分析實(shí)驗(yàn)報(bào)告姓 名學(xué) 號(hào)系 別數(shù)學(xué)系班級(jí)09信息（2）班 主講教師王丹指導(dǎo)教師王丹實(shí)驗(yàn)日期專業(yè)信息與計(jì)算科學(xué)課程名稱數(shù)值分析同組實(shí)驗(yàn)者無(wú)一、實(shí)驗(yàn)名稱： 實(shí)驗(yàn)一、插值多項(xiàng)式的收斂性實(shí)驗(yàn)二、 實(shí)驗(yàn)?zāi)康模?理解插值的基本原理；2掌握多項(xiàng)式插值的概念、存在唯一性；3編寫MATLAB程序?qū)崿F(xiàn)Lagrange插值和Newton插值，驗(yàn)證Runge現(xiàn)象、分析插值多項(xiàng)式的收斂性。三、實(shí)驗(yàn)內(nèi)容及要求：1已知數(shù)據(jù)如下：0.20.40.60.81.00.980.920.810.640.38（1）用MATLAB語(yǔ)言編寫按Langrage插值法和Newton插值法計(jì)算插值的程序，對(duì)以上數(shù)據(jù)進(jìn)行插值；（2）利用MATLAB在第一個(gè)圖中畫(huà)出離散數(shù)據(jù)及插值函數(shù)曲線。2給定函數(shù)，利用上題編好的Langrage插值程序（或Newton插值程序），分別取3個(gè)，5個(gè)、9個(gè)、11個(gè)等距節(jié)點(diǎn)作多項(xiàng)式插值，分別畫(huà)出插值函數(shù)及原函數(shù)的圖形，以驗(yàn)證Runge現(xiàn)象、分析插值多項(xiàng)式的收斂性。三、 實(shí)驗(yàn)步驟（或記錄）Lagrange插值法的基本思想：步驟1： 構(gòu)造處的插值基函數(shù)，其中，插值節(jié)點(diǎn)處的插值基函數(shù)為；步驟2：以作為的系數(shù)，使得通過(guò)插值點(diǎn)；步驟3：把所有的線性疊加，得到通過(guò)所有插值點(diǎn)的插值函數(shù)。Lagrange插值偽代碼：給定個(gè)插值點(diǎn)的情況下，求插值函數(shù)在點(diǎn)處的函數(shù)值。/*輸入?yún)?shù)*x=(x0,x1,.,xn), 插值節(jié)點(diǎn)*y=(y0,y1,yn); 被插函數(shù)f(x)在插值節(jié)點(diǎn)處的函數(shù)值*t 求插值函數(shù)Ln (x)在t處的函數(shù)值*返回值 插值函數(shù)Ln (x)在t處的函數(shù)值 */procedure Lagrange result0; for i=1 to n li(t)1; for j=1 to n if ij li(t) li(t)*(t-xi)/(xi-xj); end if end for resultresult+yi*li(t); end for return result;end procedureLagrange插值子程序 lagr1:function y=lagr1(x0,y0,x)%x0為插值點(diǎn)的向量,y0為插值點(diǎn)處的函數(shù)值向量,x為未知的點(diǎn)向量n=length(x0); m=length(x);for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j=k p=p*(z-x0(j)/(x0(k)-x0(j); end end s=p*y0(k)+s; end y(i)=s;E ndNewton插值算法公式:余項(xiàng)為其中 有關(guān).Newton插值偽代碼：/*輸入?yún)?shù)*x=(x0,x1.,xn), 插值節(jié)點(diǎn)*y=(y0,y1,yn); 被插函數(shù)f(x)在插值節(jié)點(diǎn)處的函數(shù)值*t 求插值函數(shù)Pn (x)在t處的函數(shù)值*返回值 插值函數(shù)Pn(x)在t處的函數(shù)值 */ procedure Newton for j=0 to n d1jyj;end forfor j=1 to nfor i=j to n dij (di,j-1-di-1,j-1)/(xi-xi-j+1); end forend forresultd11;temp1;for i=1to n temptemp*(t-xi-1);resultresult+di,i*temp;end forreturn result;end procedureNewton插值子程序 Newton:function y=newton(x0,y0,x)%牛頓插值法n=length(x0); m=length(x);d=zeros(n,n);%d為差商表矩陣for j=1:n d(j,1)=y0(j);%差商表第一列endfor j=2:n %差商表為下三角矩陣for i=j:n d(i,j)=(d(i,j-1)-d(i-1,j-1)./(x0(i)-x0(i-j+1);%求差商表矩陣中各值 endendfor k=1:m z=x(k); result=d(1,1); temp=1;for i=2:n temp=temp*(z-x0(i-1); result=result+d(i,i)*temp;endy(k)=result;end1編寫拉格朗日插值多項(xiàng)式函數(shù)內(nèi)容為：function f=lagfun(x)a=0.2,0.4,0.6,0.8,1.0;b=0.98,0.92,0.81,0.64,0.38;for i=1:5L(i)=1;for j=1:5 if j=i L(i)=L(i)*(x-a(j)/(a(i)-a(j); endendendf=0;for i=1:5 f=f+L(i)*b(i);end畫(huà)圖程序內(nèi)容為：x0=0.2,0.4,0.6,0.8,1.0; y0=0.98,0.92,0.81,0.64,0.38;plot(x0,y0,o)hold ongrid on fplot(lagfun,0,1);hold onx=0:0.1:1;plot(x,newton(x0,y0,x),r);legend(離散點(diǎn),Lagrange插值,Newton插值)圖形為：2Lagrange插值程序for n=3:2:11x= -1:0.1:1;y=1./(1+25.*x.2);z=0*x;x0=-1:2/(n-1):1;y0=1./(1+25.*x0.2);y1=lagr1(x0, y0, x);plot(x, z, r, x, y, k: ,x, y1, r)gtext(Lagr.,num2str(n)