西北农林科技大学数值分析报告数值法实验报告材料.docVIP

实用标准 文档大全 数值法实验报告 专业班级：信息与计算科学121 姓名：金辉 学号：2012014280 1）实验目的 本次实验的目的是熟练《数值分析》第二章“插值法”的相关内容，掌握三种插值方法：牛顿多项式插值，三次样条插值，拉格朗日插值，并比较三种插值方法的优劣。 本次试验要求编写牛顿多项式插值，三次样条插值，拉格朗日插值的程序编码，并在MATLAB软件中去实现。 2）实验题目 实验一： 已知函数在下列各点的值为 xi 0.2 0.4 0.6 .0.8 1.0 f（xi） 0.98 0.92 0.81 0.64 0.38 试用4次牛顿插值多项式P4（x）及三次样条函数S（x）（自然边界条件）对数据进行插值。用图给出{（xi，yi），xi=0.2+0.08i，i=0，1, 11, 10}，P4（x）及S（x）。 实验二： 在区间[-1,1]上分别取用两组等距节点对龙格函数作多项式插值及三次样条插值，对每个值，分别画出插值函数即的图形。 实验三： 下列数据点的插值 x 0 1 4 9 16 25 36 49 64 y 0 1 2 3 4 5 6 7 8 可以得到平方根函数的近似，在区间[0,64]上作图。 （1）用这9各点作8次多项式插值L8(x). （2）用三次样条（自然边界条件）程序求S（x）。从结果看在[0,64]上，那个插值更精确；在区间[0,1]上，两种哪个更精确？ 3）实验原理与理论基础 《数值分析》第二章“插值法”的相关内容，包括：牛顿多项式插值，三次样条插值，拉格朗日 4）实验内容 实验一： 已知函数在下列各点的值为 xi 0.2 0.4 0.6 .0.8 1.0 f（xi） 0.98 0.92 0.81 0.64 0.38 试用4次牛顿插值多项式P4（x）及三次样条函数S（x）（自然边界条件）对数据进行插值。用图给出{（xi，yi），xi=0.2+0.08i，i=0，1, 11, 10}，P4（x）及S（x）。 (1)首先我们先求牛顿插值多项式，此处要用4次牛顿插值多项式处理数据。 已知n次牛顿插值多项式如下： Pn=f(x0)+f[x0,x1](x-x0)+ f[x0,x1,x2](x-x0) (x-x1)+···+ f[x0,x1，···xn](x-x0) ···(x-xn-1) 我们要知道牛顿插值多项式的系数，即均差表中得部分均差。 在MATLAB的Editor中输入程序代码，计算牛顿插值中多项式系数的程序如下： function varargout=newtonliu(varargin) clear,clc x=[0.2 0.4 0.6 0.8 1.0]; fx=[0.98 0.92 0.81 0.64 0.38]; newtonchzh(x,fx); function newtonchzh(x,fx) %由此函数可得差分表 n=length(x); fprintf(*****************差分表*****************************\n); FF=ones(n,n); FF(:,1)=fx; for i=2:n for j=i:n FF(j,i)=(FF(j,i-1)-FF(j-1,i-1))/(x(j)-x(j-i+1)); end end for i=1:n fprintf(%4.2f,x(i)); for j=1:i fprintf(%10.5f,FF(i,j)); end fprintf(\n); end 由MATLAB计算得： xi f(xi) 一阶差商 二阶差商 三阶差商 四阶差商 0.20 0.980 0.40 0.920 -0.30000 0.60 0.810 -0.55000 -0.62500 0.80 0.640 -0.85000 -0.75000 -0.20833 1.00 0.380 -1.30000 -1.12500 -0.62500 -0.52083 所以有四次插值牛顿多项式为： P4（x）=0.98-0.3(x-0.2)-0.62500 (x-0.2)(x-0.4) -0.20833 (x-0.2)(x-0.4)(x-0.6)-0.52083 (x-0.2)(x-0.4)(x-0.6)(x-0.8） (2)接下来我们求三次样条插值函数。 用三次样条插值函数由上题分析知,要求各点的M值： 三次样条插值函数计算的程序如下： function tgsanci(n,s,t) %n代表元素数，s,t代表端点的一阶导。 x=[0.2 0.4 0.6 0.8 1.0]; y=[0.98 0.92 0.81 0.64 0.38]