数值分析matlab实验报告_第1页
数值分析matlab实验报告_第2页
数值分析matlab实验报告_第3页
数值分析matlab实验报告_第4页
数值分析matlab实验报告_第5页
已阅读5页，还剩4页未读， 继续免费阅读

下载本文档

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

文档简介

1、实验 2.1 多项式差值的振荡现象一、实验内容设区间 -1,1上函数22511)(xxf，考虑区间 -1,1 的一个等距划分，分点为nixi21，i=0,1,2,.,n，则拉格朗日插值多项式为niiinxlxl02)(2511.其中，li(x)，i=0,1,2,.,n 是 lagrange插值基函数 . 1)选择不断增大的分点数目n=2,3,.，画出原函数 f(x)及插值多项式函数ln(x)在-1,1上的图像，比较并分析实验结果. 2)选择其他的函数，例如定义在区间 -5,5上的函数41)(xxxh，xxgarctan)(，重复上述的实验看其结果如何. 二、实验程序1.主程序function

2、chapter2 promps=请选择试验函数，若选f(x)，请输入 f，若选好 h(x)，请输入 h，若选g(x)，请输入 g：; result=inputdlg(promps,charpt 2,1,f); nb_f=char(result); if(nb_f=f&nb_f=h&nb_f=g) errordlg(试验函数选择错误！ ); return; end result=inputdlg( 请输入插值多项式的次数n：,charpt_2,1,10); nd=str2num(char(result); if(nd=2) for k=3:n+1 for i=3:k+1 t(k,

3、i)=t(k-1,i)-a(k-1)*t(k-1,i-1)-b(k-1)*t(k-2,i-2); 8 end end end for i=1:n+1 for k=i:n+1 alph(n+2-i)=alph(n+2-i)+c(k)*t(k,k+2-i); end end xmin=min(x);xmax=max(x);dx=(xmax-xmin)/(25*m); t=(xmin-dx):dx:(xmax+dx); s=alph(1); for k=2:n+1 s=s.*t+alph(k); end plot(x,y,*,t,s,-); title(离散数据的多项式拟合 ); xlabel(x);ylabel(y); grid on; 三、实验结果及分析输出结果：平方误差 :2.17619e-005 参数 alph:1.99911 -2.99767 -3.96825e-005 0.549119 离散数据的拟合函数图形为：-1.5-1-0.500.511.522.

人人文库> 全部分类> 生活休闲 > 科普知识

温馨提示

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

最新文档

评论

0/150

提交评论