原创力文档- 1、本文档共13页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
建模实验拉格朗日插值
建模实验 插值方法 【实验目的】1．掌握用MATLAB计算拉格朗日、分段线性、三次样条三种插值的方法，改变节点的数目，对三种插值结果进行初步分析；2．掌握用MATLAB及梯形公式、辛普森公式计算数值积分；3．通过实例学习用插值和数值积分解决实际问题。【实验内容】题1（3）选择一些函数，在n个节点上（n不要太大，如5~11）用拉格朗日、分段线性、三次样条三种插值的方法，计算m个插值点的函数值（m要适中，如50~100）。通过数值和图形输出，将三种插值结果与精确值进行比较，适当增加n，再作比较，由此作初步分析。下列函数供选择参考：（3）y=cos10x,-2≤x≤2解：当n=8，m=50时，在matlab中编程如下：x=-2:0.5:2;y=(cos(x)).^10;a=-2:4/50:2;yy=(cos(a)).^10;y1=lagr(x,y,a);y2=interp1(x,y,a);y3=interp1(x,y,a,spline);[a;yy;y1;y2;y3]数值输出如下： -2 0 0 0 0-1.92 2.1987e-005 -0.99289 0 0.0019267 -1.84 1.771e-006 -1.2394 0 0.0027391 -1.76 5.538e-008 -1.078 8.0997e-005 0.0027867 -1.68 2.3645e-010 -0.74009 5.6075e-005 0.0022632 -1.6 4.5057e-016 -0.37345 3.1153e-005 0.0013621 -1.52 1.1388e-013 -0.063296 6.2305e-006 0-1.44 1.4242e-009 0.1509 0 -0-1.36 1.6085e-007 0.26129 0 -0.0016693 -1.28 3.7541e-006 0.28038 0 -0.0021436 -1.2 3.9028e-005 0.23209 0.0012721 -0.0020271 -1.12 0 0.14495 0.0016113 -0.001126 -1.04 0.0011051 0.047167 0.0019505 0.0007532 -0.96 0.0038502 -0.03686 0.023626 0.0039155 -0.88 0.011026 -0.088593 0.066638 0.011223 -0.8 0.026946 -0.096401 0.10965 0.028094 -0.72 0.057684 -0.055717 0.15266 0.06006 -0.64 0.11022 0.031477 0.19567 0.11265 -0.56 0.19061 0.15766 0.23869 0.19139 -0.48 0.30145 0.3111 0.30011 0.30176 -0.4 0.43942 0.4774 0.41676 0.44195 -0.32 0.59394 0.64102 0.5334 0.59603 -0.24 0.74766 0.78685 0.65005 0.74643 -0.16 0.87937 0.90159 0.7667 0.8756 -0.08 0.96847 0.97483 0.88335 0.96597 0 1 1 1 1 0.08 0.96847 0.97483 0.88335 0.96597 0.
文档评论（0）