实验 编写样条插值函数的程序
合集下载
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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)]
在编写程序时
分别引入
只有经过长时间完成其发展的艰苦工作，并长期埋头沉浸于其中的任务，方可望有所成就。——黑格尔