资源描述：
《牛顿插值法上机报告》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、为了适应公司新战略的发展，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划牛顿插值法上机报告 牛顿插值法 一、实验目的：学会牛顿插值法，并应用算法于实际问题。 二、实验内容：给定函数f(x)?x，已知： f()?()?f()? f()?()? 三、实验要求： 用牛顿插值法求4次Newton插值多项式在处的值，以此作为函数的近似值?N()。在MATLAB中用内部函数ezplot绘制出4次Newt(转载于:写论文网:牛顿插值法上机报告)on插值多项式的函数图形。 在MATLAB中用内部函数ezplot可直接绘制出以上函数的图形，并与作出
2、的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(
3、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=;[
4、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(
5、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()=。目的-通过该培训员工可对保安行业有初步了解，并感受到安保行业的发展的巨大潜力，可提升其的专业水平，并确保其在这个行业的安全感。为了适应公司新战略的发展
6、，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划 由在MATLAB中用内部函数ezplot直接绘制出出的4次Newton插值图形与原函数的图形知，4次Newton插值图形在区间[0，1]与区间[4，5]内与原函数存在一定的偏差，而在区间[1,4]内误差在10的-6次方，这个精度是非常高的。因此，在计算区间[1,4]内的值时结果是比较准确的。 河北大学数计学院 数值计算Newton插值多项式实验报告 课程名称：数值计算 授课老师：高少芹学学日 生：耿福顺号：XX期：XX/11/27 一、实验目的： 1、掌握牛顿插值法的基本思路和步骤
7、。2、培养编程与上机调试能力。二、实验原理： 给定插值点序列（xi,f(xi)),i?0,1,??,n,。构造牛顿插值多项式Nn(u)。输入要计算的函数点x,并计算Nn(x)的值，利用牛顿插值公式，当增加一个节点时，只需在后面多计算一项，而前面的计算仍有用；另一方面Nn(x)的各项系数恰好又是各阶均差，而各阶均差可用均差公式来计算。 为的一阶均差。 为 均差表： 的k阶均差。目的-通过该培训员工可对保安行业有初步了解，并感受到安保行业
