牛顿插值法上机报告.docxVIP

牛顿插值法上机报告 　　牛顿插值法　　一、实验目的：学会牛顿插值法，并应用算法于实际问题。　　二、实验内容：给定函数f(x)?x，已知：　　f()?()?f()?　　f()?()?　　三、实验要求：　　用牛顿插值法求4次Newton插值多项式在处的值，以此作为函数的近似值?N()。在MATLAB中用内部函数ezplot绘制出4次Newt(转载于:写论文网:牛顿插值法上机报告)on插值多项式的函数图形。　　在MATLAB中用内部函数ezplot可直接绘制出以上函数的图形，并与作出的4次Newton插值多项式的图形进行比较。　　四、实验过程：　　1、编写主函数。打开Editor编辑器，输入Newton插值法主程序语句：function[y,L]=newdscg(X,Y,x)　　n=length(X);z=x;A=zeros(n,n);A(:,1)=Y;s=;p=;　　forj=2:n　　fori=j:n　　A(i,j)=(A(i,j-1)-A(i-1,j-1))/(X(i)-X(i-j+1));　　end　　end　　C=A(n,n);　　fork=(n-1):-1:1　　C=conv(C,poly(X(k)));　　d=length(C);C(d)=C(d)+A(k,k);　　end　　y(k)=polyval(C,z);　　L(k,:)=poly2sym(C);　　%%%%%%%%%%%%%%%%%%　　t=[2,,,,];　　fx=sqrt(t);　　wucha=fx-Y;　　以文件名保存。　　2、运行程序。　　在MATLAB命令窗口输入：　　X=[2,,,,];Y　　=[,,,,];　　x=;[y,P]=newdscg(X,Y,x)　　回车得到：　　y=　　wucha=*　　-　　P=-(09061*x^4)/+　　(94329*x^3)/-　　(8709*x^2)/42624+　　(12817*x)/85248+　　31207/70496　　在MATLAB命令窗口输入：　　v=[0,6,-1,3];　　ezplot(P),axis(v),grid　　holdon　　x=0::6;　　yt=sqrt(x);plot(x,yt,:)　　legend(插值效果,原函数)　　xlabel(X)　　ylabel(Y)　　title(Newton插值与原函数比较)　　回车即可得到图像1-1。　　图1-1牛顿插值效果　　五、实验结果分析：　　由上运行的程序可得，用牛顿插值法求4次Newton插值多项式在处函数的近似值?N()=。　　由在MATLAB中用内部函数ezplot直接绘制出出的4次Newton插值图形与原函数的图形知，4次Newton插值图形在区间[0，1]与区间[4，5]内与原函数存在一定的偏差，而在区间[1,4]内误差在10的-6次方，这个精度是非常高的。因此，在计算区间[1,4]内的值时结果是比较准确的。　　河北大学数计学院　　数值计算Newton插值多项式实验报告　　课程名称：数值计算　　授课老师：高少芹学学日　　生：耿福顺号：XX期：XX/11/27　　一、实验目的：　　1、掌握牛顿插值法的基本思路和步骤。2、培养编程与上机调试能力。二、实验原理：　　给定插值点序列（xi,f(xi)),i?0,1,??,n,。构造牛顿插值多项式Nn(u)。输入要计算的函数点x,并计算Nn(x)的值，利用牛顿插值公式，当增加一个节点时，只需在后面多计算一项，而前面的计算仍有用；另一方面Nn(x)的各项系数恰好又是各阶均差，而各阶均差可用均差公式来计算。　　为的一阶均差。　　为　　均差表：　　的k阶均差。　　1.输入n值及（xi,f(xi)),i?0,1,??,n,；要计算的函数点x。2.对给定的x,由　　Nn(x)?f(x0)?(x?x0)f?x0,x1??(x?x0)(x?x1)f?x0,x1,x2????(x?x0)　　(x?x1)?(x?xn?1)f?x0,x1?,xn?　　计算Nn(x)的值。3.输出Nn(x)。三、源代码：　　#include#includeusingnamespacestd;#include　　voidf1(double*,double*,double**,int);voidf2(double*,double**,int);　　voidf3(double**,double*,double*,intn);　　int_tmain(intargc,_TCHAR*argv[]){coutn;double*x=newdouble[n];double*y=newdo