- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
插值方法_实验报告插值方法_实验报告
肖建 计科三班开课学院、实验室: 数统学院 实验时间 : 2011 年 5 月 8 日
课程
名称 数学实验 实验项目
名 称
插值方法
实验项目类型 验证 演示 综合 设计 其他 指导
教师 李东 成 绩
实验5 插值方法
一、实验目的及意义
[1] 了解插值的基本原理
[2] 了解拉格朗日插值、线性插值、样条插值的基本思想;
[3] 了解三种网格节点数据的插值方法的基本思想;
[4] 掌握用MATLAB计算三种一维插值和两种二维插值的方法;
[5] 通过范例展现求解实际问题的初步建模过程;
通过自己动手作实验学习如何用插值方法解决实际问题,提高探索和解决问题的能力。通过撰写实验报告,促使自己提炼思想,按逻辑顺序进行整理,并以他人能领会的方式表达自己思想形成的过程和理由。提高写作、文字处理、排版等方面的能力。
二、实验内容
1.编写拉格朗日插值方法的函数M文件;
2.用三种插值方法对已知函数进行插值计算,
通过数值和图形输出,比较它们的效果;
3.针对实际问题,试建立数学模型,并求解。
三、实验步骤
1.开启软件平台——MATLAB,开启MATLAB编辑窗口;
2.根据各种数值解法步骤编写M文件
3.保存文件并运行;
4.观察运行结果(数值或图形);
5.写出实验报告,并浅谈学习心得体会。
四、实验要求与任务
根据实验内容和步骤,完成以下具体实验,要求写出实验报告(实验目的→问题→数学模型→算法与编程→计算结果→分析、检验和结论→心得体会)
基础实验
1. 一维插值 利用以下一些具体函数,考察分段线性插值、三次样条插值和拉格朗日多项式插值等三种插值方法的差异。
1),x([-5,5]; 2)sinx, x([0,2(]; 3)cos10x, x([0,2(].
M文件:
(1)
clc
x=linspace(-5,5,11);
y=1./(1+x.^2);
x0=linspace(-5,5,101);
y0=1./(1+x.^2);
y1=interp1(x,y,x0,spline)
y2=interp1(x,y,x0);
A=[ones(11,1) x (x.^2) (x.^3) (x.^4) (x.^5) (x.^6) (x.^7) (x.^8) (x.^9) (x.^10)]
a=A\y;
y3= a(1)+a(2).*x0+a(3).*x0.^2+a(4).*x0.^3+a(5).*x0.^4+a(6).*x0.^5+a(7).*x0.^6+a(8).*x0.^7+a(9).*x0.^8+a(10).*x0.^9+a(11).*x0.^10;
plot(x0,y3,r),
gtext(Lagr.),
hold on,
plot(x0,y2,b),
gtext(Pies.Lin.),
hold on,
plot(x0,y1,m),
gtext(Spline)
hold off
(2)
x=linspace(0,2*pi,11);
y=cos(x);
x0=linspace(0,pi,101);
y0=cos(x0);
剩余代码和(1)中相同
(3)
x=linspace(0,pi,11);
y=cos(x).^10;
x0=linspace(0,pi,101);
y0=cos(x0).^10;
剩余代码和(1)中相同
注意:适当选取节点及插值点的个数;比较时可以采用插值点的函数值与真实函数值的差异,或采用两个函数之间的某种距离。
2.高维插值 对于二维插值的几种方法:最邻近插值、分片线性插值、双线性插值、三次插值等,利用如下函数进行插值计算,观察其插值效果变化,得出什么结论?
1) ,参数p=1/2000~1/200;采样步长为:t=4ms~4s;x=5~25m.
t=linspace(0,4,10);
x=linspace(5,25,10);
[t,x]=meshgrid(t,x);
z=sin((t-1/200.*x)*2);
[t1,x1]=meshgrid(linspace(0,4,100),linspace(5,25,100));
z1=interp2(t,x,z,t1,x1,cubic);
z2=interp2(t,x,z,t1,x1,nearest);
z3=interp2(t,x,z,t1,x1,linear);
mesh(t1,x1,z1);
hold on
mesh(t1,x1,z2);
hold on
mesh(t1,x1
您可能关注的文档
- 手术部位、导尿管尿路、血管导管相关血流、皮肤软组等主要部位感染具体预防控制措施手术部位、导尿管尿路、血管导管相关血流、皮肤软组等主要部位感染具体预防控制措施.doc
- 手机RF常遇的问题手机RF常遇的问题.doc
- 手性药物拆分技术的研究进展手性药物拆分技术的研究进展.doc
- 手术室护士的职业暴露与防护手术室护士的职业暴露与防护.doc
- 手机全部术语、参数详解手机全部术语、参数详解.doc
- 手机保护壳模具设计与加工手机保护壳模具设计与加工.doc
- 手机射频设计入门手机射频设计入门.doc
- 手机外壳上盖数控加工与操作手机外壳上盖数控加工与操作.doc
- 手机资费系统设计报告手机资费系统设计报告.doc
- 手机英语词汇大全手机英语词汇大全.doc
最近下载
- 2024年教师学习《中国共产党纪律处分条例》心得体会.docx VIP
- 道德与法治一年级上册 第一课时 我向国旗敬个礼 课件.pptx
- 电力设施保护宣传专题培训ppt课件.pptx VIP
- 高速路板块光伏项目可行性研究报告.docx
- 华为H19-369 HCSA-Pre-sales-IT认证考试复习(重点)题库(含答案).pdf VIP
- 【人教版】2024-2025学年一年级语文上册教学计划(及进度表).docx
- PCCP管制造及安装质量控制.pdf VIP
- 新人教版数学五年级上册全册各单元教材解读课件.pptx VIP
- 产品问卷调查模板.docx VIP
- 钢板弹簧设计计算程序Ver1.0-20070405.xls
文档评论(0)