数学实验实验报告二答案_第1页
数学实验实验报告二答案_第2页
数学实验实验报告二答案_第3页
数学实验实验报告二答案_第4页
数学实验实验报告二答案_第5页
全文预览已结束

付费下载

下载本文档

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

文档简介

实验二实验项目名称：插值与拟合实验内容：建立拉格朗日插值的函数M文件lagr.m。Newscript源程序：functiony=lagr(x0,y0,x)n=length(x0);m=length(x);fori=1:mz=x(i);s=0;fork=1:np=1;forj=1:nifj~=kp=p*(z-x0(j))/(x0(k)-x0(j));endends=p*y0(k)+s;endy(i)=s;end保存为larg.m完成机床加工问题中数据增加密度，利用三种插值方法计算加工所需数据，并画出曲线，比较分析三种插值结果。（参考教材P77）源程序：x0=[0,3,5,7,9,11,12,13,14,15];

y0=[0,1.2,1.7,2.0,2.1,2.0,1.8,1.2,1.0,1.6];

x=0:0.1:15;

y1=lagr(x0,y0,x);

y2=interp1(x0,y0,x);

y3=interp1(x0,y0,x,'spline');

[x',y1',y2',y3'];

subplot(3,1,1)

plot(x0,y0,'k+',x,y1,'r')

grid,title('lagrange')

subplot(3,1,2)

plot(x0,y0,'k+',x,y2,'r')

grid,title('piecewise

linear')

subplot(3,1,3)

plot(x0,y0,'k+',x,y3,'r')

grid,title('spline')实验结果：完成教材P80习题1中的函数y=sinx的三种插值结果的分析和比较。源程序：n=10;m=100;x=linspace(0,2*pi,100);y=sin(x);z=x*0;x0=linspace(0,2*pi,6);y0=sin(x0);y1=lagr(x0,y0,x);y2=interp1(x0,y0,x);y3=interp1(x0,y0,x,'spline');plot(x,z,x,y,'k',x,y1,x,y2,x,y3,'y')实验结果：完成教材P81习题2。源程序：x0=[0,1,4,9,16];y0=sqrt(x0);x=5;y=sqrt(5)y1=lagr(x0,y0,x)x0(1)=[];y0=sqrt(x0);y2=lagr(x0,y0,x)x0(4)=[];y0=sqrt(x0);y3=lagr(x0,y0,x)x0(1)=[];y0=sqrt(x0);y4=lagr(x0,y0,x)实验结果：y=2.2361y1=2.0794y2=2.2540y3=2.2667y4=2.2000理解教材的拟合应用实例-给药方案，参考该实例完成教材P81习题6。源程序：（τ输入成t0）t=[0.5,1,2,3,4,5,7,9];v=[6.36,6.48,7.26,8.22,8.66,8.99,9.43,9.63];y=log(10-v);a=polyfit(t,y,1);t0

人人文库> 全部分类> 应用文书 > 作业报告

温馨提示

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

最新文档

评论

0/150

提交评论