数值分析拉格朗日
《数值分析拉格朗日》由会员分享,可在线阅读,更多相关《数值分析拉格朗日(3页珍藏版)》请在装配图网上搜索。
1、数值分析实验1 练习高次插值龙格现象并作图专业/班级:数学与应用数学1105班姓名:常青学号:20110402501实验目的(1) 熟练掌握拉格朗日插值法和牛顿插值法。(2) 练习MATLAB中的相关操作。实验课时2实验类型必做实验时间 2014 年 5 月 日 时实验内容和步骤实验内容和步骤拉格朗日插值法创建M文件function d=lagrange(x,y) n=length(x);d=0;syms t;for i=1:n ji=1; for j=1:n if jiji=ji*(t-x(j)/(x(i)-x(j); end end d=d+y(i)*ji;enddx=-5:5; y=1.
2、/(1+x.2);plot(x,y,o)hold on x1=-5:0.01:5;y1=1./(1+x1.2);plot(x1,y1,k) d=lagrange(x,y);t=-5:0.01:5;yy=subs(d,t);plot(t,yy,r)hold off牛顿插值法:function f=Newton(x,y,x0,x1) syms t; if(length(x)=length(y) n=length(x); c(1:n)=0; else disp(x和y的维数不相等); return; end f=y(1); y1=0; l =1; for(i=1:n-1) for(j=i+1:n) y1(j)=(y(j)-y(i)/(x(j)-x(i); end x=0 1 2 3y=12 13 15 16x0=2.5x1=2.9f=neuton(x,y,x0,x1)结果g=15.7500g1=16.0020f=12-1/6+3/2*t2-1/3*t3实验结论用拉格朗日插值多项式和牛顿插值多项式近似表示函数都具有比较高的精度,并熟悉了MATLAB中的一些程序与算法指导教师评语(分)
- 温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。