matlab cubic spline interpolation 三次样条插值

最新推荐文章于 2021-03-20 13:34:38 发布

萱萱无敌

最新推荐文章于 2021-03-20 13:34:38 发布

阅读量1w

点赞数 3

分类专栏：学习笔记 C语言文章标签：三次样条 MATLAB

本文链接：https://blog.csdn.net/wbaibianjingling/article/details/79498109

版权

本文介绍了如何在MATLAB中使用自编代码实现三次样条插值，通过链接分享了在StackOverflow上的MATLAB代码，并且经过验证，该代码与MATLAB内置的SPLINE函数结果一致。此外，文章还提到了对于不同间隔插值点的处理方法，并提供了C语言版本的参考资料。

摘要由CSDN通过智能技术生成

在stackoverflow上看到别人自己写出的三次样条插值函数MATLAB代码，先搬运下。

链接在这里：https://stackoverflow.com/questions/7642921/cubic-spline-program

这里进行了验证，C语言程序运行结果同MATLAB端的SPLINE(x,y,xq)函数结果相同。

以下是MATLAB中的代码：

function [s0,s1,s2,s3]=cubic_spline(x,y) %x,y表示待插值点序列

if any(size(x) ~= size(y)) || size(x,2) ~= 1
   error('inputs x and y must be column vectors of equal length');
end

n = length(x);

h = x(2:n) - x(1:n-1);
d = (y(2:n) - y(1:n-1))./h;

lower = h(2:end);
main  = 2*(h(1:end-1) + h(2:end));
upper = h(1:end-1);

T = spdiags([lower main upper], [-1 0 1], n-2, n-2);
rhs = 6*(d(2:end)-d(1:end-1));

m = T\rhs;

% Use natural boundary conditions where second derivative
% is zero at the endpoints

m = [ 0; m; 0];

s0 = y;
s1 = d - h.*(2*m(1:end-1) + m(2:end))/6;
s2 = m/2;
s3 =(m(2:end)-m(1:end-1))./(6*h);

端点采用固定边界，s0