实验一函数插值方法报告概要.doc

实验一 函数插值方法报告 一、问题提出 对于给定的一元函数的n+1个节点值。试用Lagrange公式求其插值多项式或分段二次Lagrange插值多项式。 数据如下： （1） 0.4 0.55 0.65 0.80 0.95 1.05 0.41075 0.57815 0.69675 0.90 1.00 1.25382 求五次Lagrange多项式，和分段三次插值多项式，计算, 的值。（提示：结果为, ） （2） 1 2 3 4 5 6 7 0.368 0.135 0.050 0.018 0.007 0.002 0.001 试构造Lagrange多项式，计算的,值。（提示：结果为, ） 二、要求 1、 利用Lagrange插值公式 编写出插值多项式程序； 2、 给出插值多项式或分段三次插值多项式的表达式； 3、 根据节点选取原则，对问题（2）用三点插值或二点插值，其结果如何； 4、 对此插值问题用Newton插值多项式其结果如何。Newton插值多项式如下： 其中： 三、目的和意义 1、 学会常用的插值方法，求函数的近似表达式，以解决其它实际问题； 2、 明确插值多项式和分段插值多项式各自的优缺点； 3、 熟悉插值方法的程序编制； 4、 如果绘出插值函数的曲线，观察其光滑性。 四、实验学时：2学时 五、实验步骤： 1．进入C或matlab开发环境； 2．根据实验内容和要求编写程序； 3．调试程序； 4．运行程序； 5．撰写报告,讨论分析实验结果. 解： 一、编写插值函数结构程序 Lagrange插值多项式M文件：lagrange1.m function [A1,LN,L1,B1]=lagrange1(X,Y) m=length(X); LN=ones(m,m); for k=1: m x1=1; for i=1:m if k~=i x1=conv(x1,poly(X(i)))/(X(k)-X(i)); end end L1(k,:)=x1; B1(k,:)=poly2sym (x1) end A1=Y*L1;LN=Y*B1 分段三次艾尔米特插值多项式的M文件：Hermite3.m function [f,ff] = Hermite3(x,y,y1) syms t; f = 0.0; if(length(x) == length(y)) if(length(y) == length(y1)) n = length(x); else disp(y和y的导数的维数不相等); return; end else disp(x和y的维数不相等！ ); return; end for i=1:n h = 1.0; a = 0.0; for j=1:n if( j ~= i) h = h*(t-x(j))^2/((x(i)-x(j))^2); a = a + 1/(x(i)-x(j)); end end f = f + h*((x(i)-t)*(2*a*y(i)-y1(i))+y(i)); end ff = subs(f,t); （1）、求五次Lagrange多项式，和分段三次插值多项式。 在主显示区，输入五次Lagrange多项式程序： X=[0.4 0.55 0.65 0.80 0.95 1.05]; Y=[0.41075 0.57815 0.69675 0.90 1.00 1.25382]; [A1,LN,L1,B1]=lagrange1(X,Y) plot(X,A1); F=poly2sym(A1) 运行后，输出五次Lagrange多项式的结果： A1 = 121.6264 -422.7503 572.5667 -377.2549 121.9718 -15.0845 F = (2139673480305281*x^5)/17592186044416 - (1859275536318005*x^4)/4398046511104 + (9836621836743*x^3)- (414796119737013*x^2)/1099511627776 + (2145751274873259*x)/17592186044416 - 1061478972867847/7036874417