数值分析教教案7.docVIP

1.7.4 三次样条插值的MATLAB实现实例 【例1-21】在区间的两个端点处给出条件，称为边界条件。给定两端点处的导数值，然后根据下面的数据点求出其三次样条插值多项式，并计算当时的值。 1 2 3 4 5 6 7 8 0.84 0.91 0.14 -0.76 -0.96 -0.28 0.66 0.99 0.5403 -0.1455 解：在MATLAB编辑窗口中输入以下命令： x=1:8; y=[0.84 0.91 0.14 -0.76 -0.96 -0.28 0.66,0.99] [f,f0]=ThrSample1(x,y,0.54,-0.15,3.5) 程序运行后得到如下结果。 f = (-38/25*t+133/25)*(t-5)^2+(-264/25+48/25*t)*(t-4)^2+(-9447/14555*t+37788/14555)*(5-t)^2-(25502967261225075/18014398509481984-5100593452245015/18014398509481984*t)*(t-4)^2 f0 = -0.9164 程序运行时调用的三次样条函数如下（给定两端点处的一阶导数值）。 function [f,f0] = ThrSample1 (x,y,y_1, y_N,x0) syms t; f = 0.0; f0 =