- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
函数插值-习题课
解 先作差商表 xi f(xi ) 1阶 2阶 3阶 4阶 0.40 0.55 0.65 0.80 0.90 1.1160 1.1860 1.2757 1.3841 0.41075 0.57815 0.69675 0.88811 1.02652 0.2800 0.3588 0.4336 0.191 0.214 0.034 由Newton公式得四次插值多项式 * * * * 课程回顾+习题课 ? ? ? ? ? ? xj xj-1 xj+1 x0 xn 计算量与n无关; n越大,误差越小. 一般表达式 分段线性插值 * ? ? ? ? ? ? xj xj-1 xj+1 x0 xn 余项 定理:设f(x)在[a,b]上有二阶连续导数f″(x) ,则 例 已知函数 时的函数值。 上取等距节点 在区间 的 求分段线性插值函数,并由此计算 近似值。节点处函数值如下表: 0 1 2 3 4 5 1.00000 0.50000 0.20000 0.10000 0.05882 0.03846 解 分段插值基函数为 所以,分段插值函数为 与精确值 比较,结果是比较精确的。 课程回顾+习题课 * 三次样条插值大 M 方 法 两次积分 插值条件待定系数 考虑任一小区间[xi , xi+1],设 hi=xi+1-xi,Mi=S”(xi) * 获得S(x)在[xi , xi+1]上的表达式 需要知道Mi 和Mi+1的值 * 插值条件获得计算参数Mi的方程 * 第I类边界条件 f”(x0)=f”(xn),即M0=Mn=0 * 重点记忆 * 三次样条插值大 M 方 法解题步骤 1.区间划分,获得区间长度hi 2.计算 3.写出三对角方程 计算Mi 4.将Mi代入区间样条函数Si(x) 5.将要求点的值代入函数Si(x)计算 * [例] 给定离散数值表如下,取M0=M3=0构成三次样条插值的M关系式,并计算f(1.25): xi 1.1 1.2 1.4 1.5 yi 0.4000 0.8000 1.6500 1.8000 解:由题中(xi , yi)的数值可得: h0=0.1,h1=0.2,h2=0.1, 由M0=M3=0的边界条件,得 * 解得:M1=13.125,M2=-31.875。 将M0=0、M1、M2、M3=0代入区间[xi , xi+1]上的S(x): 特别地:f(1.25)?S(1.25)=1.0436 插值 习题 * A * D * B * A * B * C * * 10 已知特殊角 , , 的正弦函数分别为 求 近似值(用一、二次方法)并估计截断误差 解 角 化为弧度,分别为 。 按拉格朗日插值一次式,取 为节点,得 误差 取 为节点,按拉格朗日插值二次式, 得 误差 * 解 * * *
文档评论(0)