试求三次样条插值S(X)_第1页
试求三次样条插值S(X)_第2页
试求三次样条插值S(X)_第3页
试求三次样条插值S(X)_第4页
全文预览已结束

下载本文档

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

文档简介

1、 给定数据表如下:xj0.250.300.390.450.53yj0.50000.54770.62450.67080.7280试求三次样条插值S(X),并满足条件:i)S(0.25)=1.0000, S(0.53)-0.6868;ii) S”(0.25)= S”(0.53)=0;解: 由给定数据知: h0 =0.3-0.25 - 0.05 , h 1=0.39-0.30-0.09 h 2=0.45-0.39-0.06, h 3=0.53-0.45-0.08由 i=h i/(hi1+h i), i= h i/(hi1+h i) 得: 1= 5/14 ; 1= 9/142= 3/5 ; 2= 2/

2、5 3= 3/7 ; 3=4/70.25 0.5000 1.0000 0.25 0.5000 -0.9200-fx 0,x 0, x 1 0.9540 0.30 0.5477 -0.7193-fx 0,x 1,x 2 0.8533 0.39 0.6245 -0.5440-fx1,x2,x 3 0.7717 0.45 0.6708 -0.4050-fx 2,x 3,x 4 0.7150 0.53 0.7280 -0.3525-fx 3,x 4,x 5 0.6868 0.53 0.7280i)已知一节导数边界条件,弯矩方程组 2 1 M 0 -0.9200 5/14 2 9/14 M 1 -0.7

3、193 3/52 2/5 M 2 _6 -0.54403/7 2 4/7 M 3 -0.4050 1 2 M 4 -0.3525用追赶法解之得:M 0 = -2.0278, M1 = -1.4643 , M2 = -1.0313M3 = -0.8072, M4 = -0.6539三次样条插值函数为 1.8783 x 3 - 2.4227 x2 + 1.8591x + 0.1573, x0.25,0.30 0.8019x 3 _ 1.4538x2 + 1.5685 x + 0.1863, x0.30,0.39S(x)=0.6225x 3 _ 1.2440x2 +1.4866x + 0.1970,

4、 x0.39,0.45 0.3194x 3 _ 0.8348x2 + 1.3025x +0.2246,x0.45,0.53 ii) 已知二阶导数边界条件, M 0_ M 4 = 0,弯矩方程组 2 9/140 M1 -0.7193 3/5 22/5 M2 -6 -0.5440 03/7 2 M3 -0.4050 用追赶法解之得: M1 = -1.8809 , M2 = -0.8616, M3 = 1.0314三次样条插值函数为 -60269x 3 + 4.7023x2 _ 0.2059x + 0.3555, x0.25,0.30 1.8876x 3 _ 2.639x2 + 1.9966 x + 0.1353, x0.30,0.39S(x)=-0.4689 x 3+ 0.1178 x2 +0.9213x + 0.2751

温馨提示

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

最新文档

评论

0/150

提交评论