实验 编写样条插值函数的程序

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

只有经过长时间完成其发展的艰苦工作,并长期埋头沉浸于其中的任务,方可望有所成就。——黑格尔
实验 编写样条插值函数的程序
x= 1:10; y=sin(x)
分别设(1) x1 xn 点的导数m(1) = 0.7和?0.9
(2) x1 xn 点的2阶导数为分别为?0.8和?0.5
一、用追赶法求解 x1至xn?1点的导数值 mi,i =1,..., n -1
算法:
1) 根据补充条件来确定a1 ,b1 an , bn 。本例中,
(1)a1 = 0,b1=2m1, an = 1,bn=2m1,
(2)须根据
来确定a1 ,b1 an , bn

2)是节点的间距。向量h可以由 h = diff(x) 获得。
 根据 hi 得到ai , bi i = 2, ... , n?1

3) A1= -a1 /2, B1= b1 /2;
Ai , Bi , i = 2,..., n?1

4)

二、编写样条插值函数程序
画出xx = 1:0.05:10的样条插值图。
在输出S(xx)时,须判定区间xx所属的区间[i,i+1]后,才能按照下面的公式输出函数值。
判断用如下语句:
i = max(find(xx(k)>=x)); %这时xx(k)∈[x(i),x(i+1)]

在编写程序时
分别引入



只有经过长时间完成其发展的艰苦工作,并长期埋头沉浸于其中的任务,方可望有所成就。——黑格尔

相关文档
最新文档