Hermite插值法matlab代码并绘制出插值多项式的图像要求代码是函数形式的便于调用

以下是Hermite插值法的MATLAB代码，可以作为函数进行调用：

function [x, y] = hermite_interpolation(X, Y, DY, x_range)
% X: 插值节点横坐标
% Y: 插值节点纵坐标
% DY: 插值节点的导数值
% x_range: 插值区间
% x: 插值节点横坐标
% y: 插值节点纵坐标

n = length(X);
% 初始化矩阵
M = zeros(2*n);
V = zeros(2*n, 1);

% 填充矩阵M和向量V
for i = 1:n
    % 填充M
    M(i, 2*i-1:2*i) = [X(i)^2, X(i)];
    M(n+i, 2*i-1:2*i) = [X(i)^2, X(i)];
    if i > 1
        M(n+i-1, 2*i-3:2*i-2) = [X(i-1)^2, X(i-1)];
        M(i-1, 2*i-3:2*i-2) = [X(i-1)^2, X(i-1)];
    end
    % 填充V
    V(i) = Y(i);
    V(n+i) = DY(i);
end

% 求解线性方程组
a = M\V;

% 构造插值函数
x = linspace(x_range(1), x_range(2), 100);
y = zeros(1, length(x));
for i = 1:n
    y = y + a(2*i-1)*x.^2 + a(2*i)*x;
end

使用该函数可以进行Hermite插值，并返回插值节点的横纵坐标数组。为了绘制插值多项式的图像，可以在函数外部调用该函数，并使用`plot`命令绘图，例如：

% 定义插值节点横坐标、纵坐标和导数值
X = [0, 1, 2];
Y = [1, 4, 3];
DY = [2, 3, -1];

% 定义插值区间
x_range = [0, 2];

% 进行插值并获取插值节点
[x_interp, y_interp] = hermite_interpolation(X, Y, DY, x_range);

% 绘制插值多项式的图像
plot(x_interp, y_interp, 'b-', X, Y, 'ro');
legend('插值多项式', '插值节点');
xlabel('x');
ylabel('y');
title('Hermite插值');