数值分析实验报告插值、三次样条.docxVIP

. ... .. 实验报告：牛顿差值多项式三次样条 问题：在区间[-1,1]上分别取 n=10、20 用两组等距节点对龙格函数 f (x) 1 作多项式插 25x2 值及三次样条插值，对每个 n 值，分别画出插值函数及f (x) 的图形。 实验目的：通过编程实现牛顿插值方法和三次样条方法，加深对多项式插值的理解。应用所编程序解决实际算例。 实验要求： 认真分析问题，深刻理解相关理论知识并能熟练应用； 编写相关程序并进行实验； 调试程序，得到最终结果； 分析解释实验结果； 按照要求完成实验报告。 实验原理： 详见《数值分析 第 5 版》第二章相关容。 实验容： 牛顿插值多项式 1.1 当 n=10 时： 在 Matlab 下编写代码完成计算和画图。结果如下： 代码： clear all clc x1=-1:0.2:1; y1=1./(1+25.*x1.^2); n=length(x1); f=y1(:); for j=2:n for i=n:-1:j f(i)=(f(i)-f(i-1))/(x1(i)-x1(i-j+1)); end end syms F x p; F(1)=1;p(1)=y1(1); for i=2:n F(i)=F(i-1)*(x-x1(i-1)); p(i)=f(i)*F(i); end . .. .c . ... .. syms P P=sum(p); P10