Matlab程序Newton插值函数.doc

编写程序构造区间上的以等分结点为插值结点的Newton插值公式，假设结点数为（包括两个端点），给定相应的函数值，插值区间和等分的份数，该程序能快速计算出相应的插值公式。以，为例计算其对应的插值公式，分别取不同的值并画出原函数的图像以及插值函数的图像，观察当增大时的逼近效果. 解：Matlab计算程序为： clear clc f=input('请输入函数表达式:f(x)=','s'); %测试公式为：1/(1+25*x^2) a=input('请输入区间左端值a:'); %-1 b=input('请输入区间右端值b:'); %1 n=input('请输入区间结点数（包括两个端点）n:'); %取不同n值比较 for i=1:n x(i)=a+(b-a)/(n-1)*(i-1); y(i,1)=eval(subs(f,'x','x(i)')); end for j=1:n-1 for k=j:n-1 temp=y(k+1,j)-y(k,j); y(k+1,j+1)=temp/(x(k+1)-x(k+1-j)); end c(j)=y(j,j); c(j+1)=y(j+1,j+1); end p=c(1); q=1; syms X for i=2:n q=q*(X-x(i-1)); p=p+c(i)*q; end p=simple(p) for i=1:301 t(i)=a+(b-a)/300*(i-1); Nn(i)=eval(subs(p,'X','t(i)')); end for i=1:301 h(i)=a+(b-a)/300*(i-1); yy(i)=eval(subs(f,'x','h(i)')); end plot(h,yy,'r') hold on plot(t,Nn,'b') hold on grid on legend('?ê??ú??f(x)','2??μ?ú??N(x)') title('?￡?ù2??μ') xlabel('x') ylabel('f(x)')  当n=5时，Newton插值公式为： p =(1250*X^4)/377 - (3225*X^2)/754 + 1 Matlab绘制的拟合图像为： 由上图可见，n取较小值时，拟合误差较大 当n=10时，Newton插值公式为： p = (84*X^9 + *X^8 - 112*X^7 - 4*X^6 + 136*X^5 + *X^4 + 44*X^3 - *X^2 + 4*X + 9741593257)/2706126848 Matlab绘制的拟合图像为： 由上图可见，随着n的增加，曲线拟合情况变好，且曲线两端拟合情况不如中间好。 当n=15时，Newton插值公式为： p =-(68608*X^14 - 6*X^13 - *X^12 + 0*X^11 + *X^10 - 2*X^9 - *X^8 + 8*X^7 + 42944*X^6 - *X^5 - 3496*X^4 - 33137*X^3 + 616*X^2 - 2313*X - 18)/16 Matlab绘制的拟合图像为： 由图可见，随着n的增加，曲线中部的拟合情况更好，曲线在两端出现了严重的龙格现象，在（-0.5，0.5）区间，曲线拟合情况最好 当n=20时，Newton插值公式为： p = -(0*X^19 + 1184*X^18 - 56*X^17 - 8304*X^16 + 04*X^15 + 59584*X^14 - 72*X^13 - 79328*X^12 + 20*X^11 + 3680*X^10 - 2*X^9 - 6016*X^8 + *X^7 + 248*X^6 - 92548*X^5 - 48*X^4 + 5207*X^3 + 5*X^2 - 99*X - 3991381301)/5446148096 Matlab绘制的拟合图像为： 由上图可见，随着n的增加，曲线中部的拟合情况更好，但两侧曲线龙格现象更加明显。  当n=40时，Matlab绘制的拟合图像为： 由图可见，当n增加时，龙格现象不断加深。当n增大到一定程度，由于两端龙格现象过于严重，导致曲线中部拟合情况不明显。