拉格朗日插值法_matlab.doc_第1页
拉格朗日插值法_matlab.doc_第2页
拉格朗日插值法_matlab.doc_第3页
全文预览已结束

下载本文档

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

文档简介

MATLAB结课作业姓名：郭海阳班级：机械093学号：2009111006成绩：时间：2012/6/8一任务。用matlab编写拉格朗日插值算法的程序 并且以（x=-2.00,f(x)=17.00 x=0.00,f(x)=1.00 x=1.00,f(x)=2.00 x=2.00,f（x)=17.00)为数据基础，在整个插值区间上采用拉格朗日插值算法计算f(x=0.6)，写出程序源代码，输出计算结果二算法。x0=-2.00;x1=0.00;x2=1.00;x3=2.00;y0=17.00;y1=1.00;y2=2.00;y3=17.00;x=0.6y=(x-x1).*(x-x2).*(x-x3)/(x0-x1).*(x0-x2).*(x0-x3)*y0+(x-x0).*(x-x2).*(x-x3)/(x1-x0).*(x1-x2).*(x1-x3)*y1+(x-x0).*(x-x1).*(x-x3)/(x2-x0).*(x2-x1).*(x2-x3)*y2+(x-x0).*(x-x1).*(x-x2)/(x3-x0).*(x3-x1).*(x3-x2)*y3;disp(y=);disp(y);结果为：x = 0.6000y= 0.2560三.程序。function s=Lagrange(x,y,x0)%lagrange插值，x，y为已知的插值点及其函数值%x0为要求的插值点的x值nx=length(x);ny=length(y);if nx=ny warning(矢量x与y的长度应该相等) returnendm=length(x0);%按照公式，对要求的插值点矢量x0的每个元素进行计算for i=1:m t=0.0; for j=1:nx u=1.0; for k=1:nx if k=j u=u*(x0(i)-x(k)/(x(j)-x(k); end end t=t+u*y(j); end s(m)=t;endreturn四算例及分析。给定x,y,x0nxlength(x
人人文库> 全部分类> 应用文书 > 技术指导

温馨提示

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

最新文档

评论

0/150

提交评论