下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、.用MATLAB实现拉格朗日插值和分段线性插值1、 实验内容: 用MATLAB实现拉格朗日插值和分段线性插值。2、 实验目的:1) 学会使用MATLAB软件;2) 会使用MATLAB软件进行拉格朗日插值算法和分段线性差值算法;3、实验原理:利用拉格朗日插值方法进行多项式插值,并将图形显式出来。4、实验步骤及运行结果(1)实现lagrange插值1)定义函数: f = 1/(x2+1) 将其保存在f.m 文件中,具体程序如下:function y = f1(x)y = 1./(x.2+1);2) 定义拉格朗日插值函数:将其保存在lagrange.m 文件中,具体实现程序编程如下:function
2、 y = lagrange(x0,y0,x)m = length(x); /区间长度/n = length(x0);for i = 1:nl(i) = 1;endfor i = 1:mfor j = 1:nfor k = 1:n if j = k continue; endl(j) = ( x(i) -x0(k)/( x0(j) - x0(k) )*l(j);endendendy = 0;for i = 1:ny = y0(i) * l(i) + y; end3) 建立测试程序,保存在text.m文件中,实现画图:x=-5:0.001:5; y=(1+x.2).-1; p=polyfit(x,
3、y,n); py=vpa(poly2sym(p),10) plot_x=-5:0.001:5; f1=polyval(p,plot_x); figureplot(x,y,r,plot_x,f1)输入n=6,出现下面的图形:通过上图可以看到当n=6是没有很好的模拟。于是重新运行text.M并选择n=11由此可见n=11时的图像是可以很好的实现模拟(2)分段线性插值: 建立div_linear.m文件。具体编程如下/*分段线性插值函数:div_linear.m 文件*/function y = div_linear(x0,y0,x,n)%for j = 1:length(x)for i = 1:n
4、-1 if (x = x0(i) & (x = x0(i+1) y = (x - x0(i+1)/(x0(i) - x0(i+1)*y0(i) + ( x - x0(i)/(x0(i+1) - x0(i)*y0(i+1); else continue; endend%end 测试程序(text2.m): n = input(输入n =:);x0 = linspace( -5,5,n);for x = -5:0.01:5 y = div_linear(x0,f(x0),x,n); hold on;plot(x,y,r); plot(x,f(x),b);end2)运行测试程序,这是会出现:输入n=
5、:2)输入n=6,并按Enter键,出现:4)关掉图形界面后,重新运行程序,输入n=11,并按enter键后出现:5)再次关掉图形界面,输入n=100,并按enter键,出现:此时。图形将于原函数图形基本吻合,说明分割区间越多,图像接近真实的图像。(3)用lagrange插值观察y = |sin(k*x)|的误差分析:1)编写 函数文件,保存在f2.m 中x=0:0.01:1;k= input(输入k:)n= input(输入n:);y=abs(sin(k*pi*x);p=polyfit(x,y,n-1);py=vpa(poly2sym(p),8);plot_x=0:0.01:1;f1=polyval(p,plot_x);plot(x,y,plot_x,f1); 2)运行该程序:输入k=:1输入n=:2出现如下图形界面:关掉图形界面后重新运行f2.m,输入k=:1,n=:3出现如下界面: 再次关掉图形界面,输入k=:1,n=:6 后出现:此时图形基本吻合。类推,输入k=2, n=3后出现: k =2, n =11,出现如下图形: k =2,n =15,这时出现:k =2,n =19,出现:当k=2,n=21时,图形如下:此时基本吻合。5、实验总结: 通过本次课程设计,我初步掌握了MATLAB运用,加深了
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 淄博市房地产经纪服务合同(2篇)
- 西华师范大学《环境法学》2023-2024学年期末试卷
- 19 一只窝囊的大老虎 说课稿-2024-2025学年语文四年级上册(统编版)
- DB51-T 2933-2022 未成年人社区矫正规范
- 保健品 批件 使用 合同
- 4050社保劳务合同
- 《班级公约大家定》(说课稿)沪科黔科版三年级上册综合实践活动
- 包工电缆安装合同范例
- 保障性住房收购合同范例
- 前台桌子合同范例
- 2024-2025学年二年级上学期数学期中模拟试卷(苏教版)(含答案解析)
- 礼修于心 仪养于行 课件-2023-2024学年高一上学期文明礼仪在心中养成教育主题班会
- 入团志愿书(2016版本)(可编辑打印标准A4) (1)
- 北京地铁受电弓的维护与故障检修-毕业设计说明书
- 幼儿园教育和家庭教育的有效结合研究
- 集团公司两金管理评价办法
- 电影的声音分析PPT课件
- “三措一案”实施规范标准
- 【全面解读《国有建设用地使用权出让地价评估技术规范【2018】4号文》
- 粉煤灰资源化综合利用新技术
- 案件移交清单模板
评论
0/150
提交评论