MATLAB实现Hermite分段三次插值及图像绘制

本文介绍如何使用MATLAB对函数f(x)=1/(1+x²)在区间[-5, 5]上进行分段三次Hermite插值，并绘制出插值函数的图像。

**代码实现：**matlab% 定义插值函数f = @(x) 1./(1+x.^2);

% 定义插值区间和节点数量a = -5; % 插值区间左端点b = 5; % 插值区间右端点n = 15; % 节点数量

% 生成等距节点x = linspace(a, b, n);

% 计算节点处的函数值和导数值y = f(x);dy = (f(x+eps)-f(x-eps))/(2*eps); % 使用中心差分近似计算导数

% 使用Hermite插值进行分段三次插值xx = linspace(a, b, 1000); % 更细的插值点用于绘图yy = zeros(size(xx));

for i = 1:n-1 % 计算插值多项式系数 h = x(i+1) - x(i); p = (xx - x(i)) / h; m0 = (1 - 3p.^2) / h^2; m1 = (xx - x(i)) . (p - 1).^2 / h; m2 = (3p.^2 - 2p) / h; m3 = (xx - x(i)) .* p.^2 / h; % 计算插值结果 yy((xx >= x(i)) & (xx <= x(i+1))) = yy((xx >= x(i)) & (xx <= x(i+1))) ... + y(i)*m0((xx >= x(i)) & (xx <= x(i+1))) + y(i+1)*m2((xx >= x(i)) & (xx <= x(i+1))) ... + dy(i)*m1((xx >= x(i)) & (xx <= x(i+1))) + dy(i+1)*m3((xx >= x(i)) & (xx <= x(i+1)));end

% 绘制图像plot(xx, yy, 'b-', x, y, 'ro');legend('Hermite分段三次插值', '原始函数');xlabel('x');ylabel('f(x)');title('f(x)的Hermite分段三次插值近似');grid on;

代码解读：

首先定义目标函数 f(x) = 1/(1+x²)，并设置插值区间为 [-5, 5]，节点数量为 15。2. 使用 linspace 函数生成等距节点，并计算节点处的函数值和导数值。3. 利用循环对每个子区间进行 Hermite 插值计算，并将结果存储在 yy 数组中。4. 最后使用 plot 函数绘制原始函数和插值函数的图像，并添加图例、坐标轴标签和标题。

结果分析：

运行代码后，会生成一张图像，其中蓝色曲线表示 Hermite 分段三次插值函数，红色圆点表示原始函数在等距节点处的函数值。可以看出，插值函数曲线光滑地连接了所有节点，并很好地逼近了原始函数。

总结：

本文介绍了如何使用 MATLAB 实现 Hermite 分段三次插值，并通过绘制图像直观地展示了插值效果。这种方法可以有效地逼近复杂函数，并在数值计算和数据分析中有着广泛的应用。