哈工大_计算方法实验_拉格朗日

哈工大_计算方法实验_拉格朗日

ID:15019500

大小:107.50 KB

页数:15页

时间:2018-07-31

哈工大_计算方法实验_拉格朗日_第1页
哈工大_计算方法实验_拉格朗日_第2页
哈工大_计算方法实验_拉格朗日_第3页
哈工大_计算方法实验_拉格朗日_第4页
哈工大_计算方法实验_拉格朗日_第5页
资源描述:

《哈工大_计算方法实验_拉格朗日》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、实验题目1Lagrange插值摘要给定平面上n+1个不同的数据点:则满足条件的n次拉格朗日插值多项式是存在唯一的。若,且充分光滑,则当时,有误差估计式前言利用拉格朗日插值多项式求的近似值程序设计流程拉格朗日插值框图开始结束输入(xi,yi),ni=0,1,2,…,nL=0i=0xl=1xl=*xlj=0,1,…,i-1,i+1,…,nL=L+xl*i=n?输出y否是i=i+1问题1(1)N=5时,程序运行如下:TestLag(inline('1./(1+x.^2)'),-5,5,5,0.75:4.75);将区间[-5,5]分为了5段计算插值的点xi=0.75001.750

2、02.75003.75004.7500计算出的插值yi=0.90540.52580.0096-0.3568-0.1595插值点处函数值yFact=0.64000.24620.11680.06640.0424计算误差err=-0.2654-0.27960.10720.42320.2020N=10时,程序运行如下:TestLag(inline('1./(1+x.^2)'),-5,5,10,0.75:4.75);将区间[-5,5]分为了10段计算插值的点xi=0.75001.75002.75003.75004.7500计算出的插值yi=0.69070.23300.11220.1

3、084-0.2360插值点处函数值yFact=0.64000.24620.11680.06640.0424计算误差err=-0.05070.01320.0045-0.04200.2785N=20时,程序运行如下:TestLag(inline('1./(1+x.^2)'),-5,5,20,0.75:4.75);将区间[-5,5]分为了20段计算插值的点xi=0.75001.75002.75003.75004.7500计算出的插值yi=0.64130.24910.12820.19036.4150插值点处函数值yFact=0.64000.24620.11680.06640.04

4、24计算误差err=-0.0013-0.0029-0.0114-0.1239-6.3726问题1(2)N=5时,程序运行如下:TestLag(inline('exp(x)'),-1,1,5,[-0.95-0.050.050.95]);将区间[-1,1]分为了5段计算插值的点xi=-0.9500-0.05000.05000.9500计算出的插值yi=0.38630.95131.05122.5863插值点处函数值yFact=0.38670.95121.05132.5857计算误差err=1.0e-003*0.4471-0.10510.1069-0.6129N=10时,程序运行

5、如下:TestLag(inline('exp(x)'),-1,1,10,[-0.95-0.050.050.95]);将区间[-1,1]分为了10段计算插值的点xi=-0.9500-0.05000.05000.9500计算出的插值yi=0.38670.95121.05132.5857插值点处函数值yFact=0.38670.95121.05132.5857计算误差err=1.0e-008*-0.3126-0.0055-0.0055-0.3714N=20时,程序运行如下:TestLag(inline('exp(x)'),-1,1,20,[-0.95-0.050.050.95]

6、);将区间[-1,1]分为了20段计算插值的点xi=-0.9500-0.05000.05000.9500计算出的插值yi=0.38670.95121.05132.5857插值点处函数值yFact=0.38670.95121.05132.5857计算误差err=1.0e-012*0.73390-0.0002-0.5671问题2(1)N=5时,程序运行如下:TestLag(inline('1./(1+x.^2)'),-1,1,5,[-0.95-0.050.050.95]);将区间[-1,1]分为了5段计算插值的点xi=-0.9500-0.05000.05000.9500计算出

7、的插值yi=0.51360.99780.99780.5136插值点处函数值yFact=0.52560.99750.99750.5256计算误差err=0.0121-0.0002-0.00020.0121N=10时,程序运行如下:TestLag(inline('1./(1+x.^2)'),-1,1,10,[-0.95-0.050.050.95]);将区间[-1,1]分为了10段计算插值的点xi=-0.9500-0.05000.05000.9500计算出的插值yi=0.52430.99750.99750.5243插值点处函数值yF

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

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

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