数值分析实验报告-插值、逼近.pdfVIP

.. 实验报告：函数逼近 插值多项式补充 1 2k 1 问题 1 ：对于给函数 f (x) 2 ，取点 xk cos ，k 取 0， 1，…，n 。n 取 10 1+25x 2n 2 或 20 。试画出拟合曲线并打印出方程，与第二章计算实习题 2 的结果进行比较。 1 i 问题 2 ：对于给函数 f (x ) 2 在区间 [-1 ， 1]上取 x=-1+0.2i （i=0,1,2, …,10 ），试求 3 1+25x 次曲线拟合，试画出拟合曲线并打印出方程，与第二章计算实习题 2 的结果进行比较。 实验目的： 通过编程实现牛顿插值方法和函数逼近， 加深对多项式插值的理解。 应用所编程 序解决实际算例。 实验要求： 1． 认真分析问题，深刻理解相关理论知识并能熟练应用； 2 ． 编写相关程序并进行实验； 3 ． 调试程序，得到最终结果； 4 ． 分析解释实验结果； 5 ． 按照要求完成实验报告。 实验原理： 详见《数值分析 第 5 版》第二章、第三章相关内容。 实验内容： （1 ）问题 1： 这里我们可以沿用实验报告一的代码，对其进行少量修改即可。 当 n=10 时，代码为： clear all clc k=0:10; n=length(k); x1=cos((2*k+1)/2/n*pi); y1=1./(1+25.*x1.^2); 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); a .. end syms P P=sum(p); P10=vpa(expand(P),5); x0=-1:0.001:1; y0=subs(P,x,x0); y2=subs(1/(1+25*x^2),x,x0); plot(x0,y0,x0,y2) grid on xlabel( x ) ylabel( y ) 由 此 我 们 可 以 得 到 10 P (x)=-46.633*x^10+3.0962e-14*x^9+130.11*x^8-7.2714e-14*x^7-133.44*x^6+7.1777e- 14*