牛顿插值法实验报告_第1页
牛顿插值法实验报告_第2页
牛顿插值法实验报告_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、精选优质文档-倾情为你奉上牛顿插值法一、实验目的：学会牛顿插值法，并应用算法于实际问题。二、实验内容：给定函数 ，已知： 三、实验要求： （1）用牛顿插值法求4次Newton插值多项式在2.15处的值，以此作为函数的近似值。在MATLAB中用内部函数ezplot绘制出4次Newton插值多项式的函数图形。 （2）在MATLAB中用内部函数ezplot可直接绘制出以上函数的图形，并与作出的4次Newton插值多项式的图形进行比较。四、实验过程：1、编写主函数。打开Editor编辑器，输入Newton插值法主程序语句： function y,L=newdscg(X,Y,x) n=length(X)

2、; z=x; A=zeros(n,n);A(:,1)=Y's=0.0; p=1.0; for j=2:n for i=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); for k=(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,2.1,2.2,2.3,2.4; fx=sqrt(t); wucha=fx-Y;以文件名newdscg.

3、m保存。2、运行程序。（1）在MATLAB命令窗口输入： >> X=2,2.1,2.2,2.3,2.4; Y =1.,1.,1.,1.,1.; x=2.15;y,P=newdscg(X,Y,x)回车得到： y =1.4663 wucha =1.0e-06 * -0.4376 -0.3254 -0.3026 0.0888 0.3385 P = - (09061*x4)/ + (94329*x3)/ - (8709*x2)/42624 + (12817*x)/85248 + 31207/70496（2）在MATLAB命令窗口输入： >> v=0,6,-1,3; >&g

4、t; ezplot(P),axis(v),grid >> hold on >> x=0:0.1:6; >> yt=sqrt(x);plot(x,yt,':') >> legend('插值效果','原函数') >> xlabel('X') >> ylabel('Y') >>title('Newton插值与原函数比较')回车即可得到图像1-1。图1-1牛顿插值效果五、实验结果分析： 由上运行（1）的程序可得，用牛顿插值法求4次Newton插值多项式在2.15处函数的近似值=1.4663。 由在MATLAB中用内部函数ezplot直接绘制出出的4次Newton插值图形与原函数的图形知，4次Newton插值图形在区间0，1与区间4，5内与原函数

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论