Matlab中MySpline函数定义编写

在 MATLAB 中，可以使用以下代码编写一个名为 MySpline 的自定义函数，该函数可以计算给定数据集的样条插值：

function [spline] = MySpline(x, y, xx)
% 计算给定数据集的样条插值
% x: 数据点的 x 坐标
% y: 数据点的 y 坐标
% xx: 用于插值的 x 坐标

% 计算差分矩阵
n = length(x);
h = diff(x);
dy = diff(y);
delta = 3 * (dy(2:end)./h(2:end) - dy(1:end-1)./h(1:end-1));

% 计算三对角矩阵
mu = zeros(n, 1);
lambda = zeros(n, 1);
d = zeros(n, 1);
mu(2:end-1) = h(2:end-1);
lambda(2:end-1) = h(1:end-2);
d(2:end-1) = delta;

% 求解三对角矩阵
A = spdiags([mu d lambda], [1 0 -1], n-2, n-2);
M = [0; linsolve(A, delta); 0];

% 计算样条插值
spline = zeros(size(xx));
for i = 1:length(xx)
    j = find(x<=xx(i), 1, 'last');
    if isempty(j)
        j = 1;
    elseif j == n
        j = n-1;
    end
    t = (xx(i) - x(j)) / h(j);
    spline(i) = (1-t)*y(j) + t*y(j+1) + t*(1-t)*((1-t)*M(j) + t*M(j+1));
end

函数的输入参数为数据点的 x 坐标 `x`、y 坐标 `y`，以及用于插值的 x 坐标 `xx`。函数输出为样条插值 `spline`。函数首先计算差分矩阵，然后根据差分矩阵计算三对角矩阵，并求解该矩阵以获得插值的系数。最后，函数计算每个插值点的样条插值值。