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

实验 编写样条插值函数的程序 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=2m （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)] 在编写程序时 分别引入