数值分析实验报告1.docVIP

数学与信息工程学院 实 验 报 告 课程名称： 计算方法 实 验 室： 5206 实验台号： 班 级： 姓 名： 实验日期： 2014 年 2　月 26 日 实验名称 插值 实验目的 和 要 求 1.插值函数及其在定点的函数值计算 2.完成实验内容 3..准备实验报告 实验内容和步骤： 实验内容：1.熟悉与插值有关的Matlab函数 2.设计插值程序 3.数据实验 Lagrange Interpolation 实验题目： (1)、Find the Lagrange interpolating polynomial respectively for the following data set, and give approximate value to f(0). D={(0,5),(1,-2),(2,3),(3,4),(4,-1),(5,7),(6,5),(7,2)} (2)、What are the Lagrange from of the interpolating polynomial for the data in the following table, and give approximate value to f(7) x 5 -7 -6 0 y 1 -23 -54 -954 求Lagrange 多项式的程序： function yy=lagrintp(x,y,xx) n=length(x);m=length(y); if m~=n error(向量x与y的长度不一致) end S=0; for k=1:n p=1; for j=1:n if j~=k p=conv(p,poly(x(j)))/(x(k)-x(j)); end end S=S+p*y(k); end yy=poly2sym(S) y0=subs(yy,x,xx) 3、实验结果： （1）、 lagrintp(0:7,[5,-2,3,4,-1,7,5,2],0) yy = 23/560*x^7-691/720*x^6+2089/240*x^5-5537/144*x^4+1689/20*x^3-7343/90*x^2+4369/210*x+5 y0 = 5 ans = 23/560*x^7-691/720*x^6+2089/240*x^5-5537/144*x^4+1689/20*x^3-7343/90*x^2+4369/210*x+5 （2）、 lagrintp([5,-7,-6,0],[1,-23,-54,-954],7) yy = 4*x^3+35*x^2-84*x-954 y0 = 1545 ans = 4*x^3+35*x^2-84*x-954 二、Newton interpolating 1、实验题目： (1)、Find the cubic interpolating polynomial respectively for the following data set,and give approximate value to f(0). D={(0,5),(1,-2),(2,3),(3,4),(4,-1),(5,7),(6,5),(7,2)} (2)、Find the quadratic interpolating polynomial for the data in the following table, and give approximate value to f(7) x 5 -7 -6 0 y 1 -23 -54 -954 求牛顿多项式的程序 建立差商表： function D=divdiff(x,y) n=length(y); if length(x)~=n fprintf(出错信息：节点向量与函数值向量的维数不相等.); end D=zeros(n,n); D(:,1)=y(:); for j=2:n; for i=j:n D(i,j)=(D(i,j-1)-D(i-1,j-1))/(x(i)-x(i-j+1)); end 求m次牛顿差值多项式 function [D,yy]=newton(x,y,xx,m) n=length(y); if length(x)~=n