数学实验——插值与数值积分

数学实验——插值与数值积分

ID:41506472

大小:536.26 KB

页数:18页

时间:2019-08-26

数学实验——插值与数值积分_第1页
数学实验——插值与数值积分_第2页
数学实验——插值与数值积分_第3页
数学实验——插值与数值积分_第4页
数学实验——插值与数值积分_第5页
资源描述:

《数学实验——插值与数值积分》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库

1、数学实验报告实验3插值与数值积分实验3插值与数值积分分1黄浩2011011743一、实验目的1)掌握用matlab计算Lagrange、分段线性、分段三次和三次样条插值这四种插值方法,并通过改变节点的数目,分析插值结果2)掌握用matlab及梯形公式、辛普森公式计算数值积分3)通过实例学习用插值和数值积分解决实际问题二、实验内容1.(1)题:问题叙述:考虑函数:fx=11+x2x∈[-5,5],取不同的节点数目,分别用Lagrange插值、分段线性和三次插值,以及三次样条插值近似f(x),分析插值效果方案实施:由于在matlab中,没有预设的l

2、agrange插值法函数,因此我首先编制了Lagrange(程序见四.1)和分段线性的插值函数(程序见四.2),并进行了y=x2函数的插值检验(程序见四.3),以检验编写的函数是否准确,如下图所示:18数学实验报告实验3插值与数值积分上图中,红、绿两线分别是自设Lagrange和分段线性的插值图,蓝色十字点是matlab自带的分段线性函数的插值图,蓝绿完全吻合、红绿基本吻合,证明自设函数是正确的。然后进行本题目要求的插值操作:以(x0,y0)为节点,节点在区间内均匀分布且将节点数分别调为6、11、21;x为插值点,绘图时使用101个插值点(始终

3、大于节点数),制表时使用的插值点数目依节点数调节;y为真值,y1、y2、y3、y4分别为使用Lagrange、分段线性、分段三次和三次样条插值的方法所得的插值结果。1)节点数为6时:(即节点为-5,-3,…,3,5)xyy1y2y3y401.00000.56730.50000.50000.56850.50000.80000.55010.50000.50000.55131.00000.50000.50000.50000.50000.50001.50000.30770.42130.40000.44250.41672.00000.20000.3212

4、0.30000.31330.31312.50000.13790.20970.20000.17750.20293.00000.10000.10000.10000.10000.10003.50000.07550.00780.08460.07540.01834.00000.0588-0.04810.06920.0559-0.02854.50000.0471-0.04600.05380.0431-0.02635.00000.03850.03850.03850.03850.0385图像为:18数学实验报告实验3插值与数值积分由上图可见,当节点数为6时,插

5、值效果皆不理想,在各节点之间的插值都有很大的误差,在x=0附近误差最大1)节点数为11时:(即节点为-5,-4,…,4,5)xyy1y2y3y401.00001.00001.00001.00001.00000.50000.80000.84340.75000.79690.82051.00000.50000.50000.50000.50000.50001.50000.30770.23530.35000.32190.29732.00000.20000.20000.20000.20000.20002.50000.13790.25380.15000.13

6、850.14013.00000.10000.10000.10000.10000.10003.50000.0755-0.22620.07940.07550.07454.00000.05880.05880.05880.05880.05884.50000.04711.57870.04860.04650.04845.00000.03850.03850.03850.03850.0385图像为:由上图可见,当节点数为11时,分段三次与三次样条插值的效果较好,其次是分段线性插值,而Lagrange插值在[-2,2]区间内精度较高,但在

7、x

8、>3时,就会发生R

9、unge振荡18数学实验报告实验3插值与数值积分1)节点数为21时(程序见四.4):(即节点为-5,-4.5…,4.5,5)xyy1y2y3y401.00001.00001.00001.00001.00000.10000.99010.99040.96000.98690.98910.20000.96150.96260.92000.95260.95940.30000.91740.91890.88000.90500.91520.40000.86210.86320.84000.85150.86060.50000.80000.80000.80000.80

10、000.80000.60000.73530.73370.74000.74560.73720.70000.67110.66820.68000.68230

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。