拉格朗日插值法_第1页
拉格朗日插值法_第2页
拉格朗日插值法_第3页
拉格朗日插值法_第4页
拉格朗日插值法_第5页
已阅读5页,还剩3页未读 继续免费阅读

下载本文档

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

文档简介

1、matlab实验报告 日期:2013年6月6日 班级acca22姓名何思佑学号2120805056实验名称 插值法与拉格朗日差值多项式实验目的:1.掌握求拉格朗日插值多项式的方法; 2.用拉格朗日插值多项式估计函数的值。实验原理与数学模型:n次拉格朗日插值多项式为:ln(x)=y0l0(x)+y1l1(x)+y2l2(x)+ynln(x)n=1时,称为线性插值,l1(x)=y0(x-x1)/(x0-x1)+ y1(x-x0)/(x1-x0)=y0+(y1-x0)(x-x0)/(x1-x0)n=2时,称为二次插值或抛物线插值,精度相对高些l2(x)=y0(x-x1)(x-x2)/(x0

2、-x1)/(x0-x2)+y1(x-x0)(x-x2)/(x1-x0)/(x1-x2)+y2(x-x0)(x-x1)/(x2-x0)/(x2-x1)实验内容:1. 拉格朗日差值法算法分析 a. 输入节点的个数 b. 输入节点的横纵坐标 c. 输入新插入的节点的横坐标 d. 通过几次循环求的新插入节点的纵坐标实验所用软件及版本:matlab2011实验过程记录:源代码function y=lagrange(x0,y0,x)n=length (x0),m=length(x);for i=1:mz=x(i);s=0,0;for k=1:np=1.0;for j=1:nif j=kp=p*(z-x0(

3、j)/(x0(k)-x0(j);endends=p*y0(k)+s;endy(i)=s;end实验截图调试数据及结果>> x0=1 2 3 4 5;y0=3 5 8 19 6;x=9;>> y=a(x0,y0,x)n = 5s = 0y = -2.290999999999999e+003实验总结:通过做此实验,我对matlab有了更进一步的了解,学会怎样才能正确运用matlab求解实际问题,我对matlab还不是特别熟悉,过程中遇到了很多问题,经过与同学讨论得到了良好的解决,希望以后会有进步。思考与深入:接触matlab上机实验,我才了解了数学软件的功能,也知道了软件在

4、我们日常生活当中的重要性matlab不仅有强大的计算功能,还有强大的绘图功能,在matlab中,有二维曲线绘图命令plot、散点图命令scatter等等、利用这些命令可以绘制出复杂的图像。 虽然开始时候,我觉得数学建模很枯燥,很乏味,但是慢慢了解了matlab软件基础和功能后,我越发喜欢这个看似无所不能的软件。随着对软件的不断深入,我觉得matlab软件还是很有意思的,即使matlab软件界面全部是英文,而且有很多专业的词汇,很多地方作为初学者的我还看不太懂,特别是一些细节方面的问题,比如“:”“;”的区别、“.*”和“*”的区别等等,但随着我一边上网查阅相关资料,一边解决老师的上机作业,我体会到在面对不知道的问题的时候要学会自己去寻找方法解决。同时,通过使用mat

温馨提示

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

最新文档

评论

0/150

提交评论