本文通过MATLAB展示了条件插值（包括complete、second、periodic、not-a-knot、variational），二维插值（如nearest、cubic、spline）以及多项式拟合和最小二乘拟合的实现方法，详细解释了不同插值和拟合方式的适用场景和效果。
%% 样条边界条件插值
x=linspace(0,2*pi,15);
y=sin(x);
plot(x,y,'o')
hold on
% complete
p=csape(x,y,'complete',[2,2]);
x1=linspace(0,2*pi,150);
y1=ppval(p,x1);
plot(x1,y1,'b','linewidth',1.5)
% second
p=csape(x,y,'second',[2,3]);
x1=linspace(0,2*pi,150);
y1=ppval(p,x1);
plot(x1,y1,'b','linewidth',1.5)
% periodic
p=csape(x,y,'periodic',[2,2]);
x1=linspace(0,2*pi,150);
y1=ppval(p,x1);
plot(x1,y1,'b','linewidth',1.5)
% not-a-knot
p=csape(x,y,'not-a-knot');
x1=linspace(0,2*pi,150);
y1=ppval(p,x1);
plot(x1,y1,'b','linewidth',1.5)
%variational
p=csape(x,y,'variational');
x1=linspace(0,2*pi,150);
y1=ppval(p,x1);
plot(x1,y1,'b','linewidth',1.5)
%% interp2
x=linspace(-5,5,5)
y=linspace(-5,5,5)
[X,Y]=meshgrid(x,y);
Z=X.^2+Y.^2;
x1=linspace(-5,5,100)
y1=linspace(-5,5,100)
[XX,YY]=meshgrid(x1,y1);
% ZZ&
最低0.47元/天 解锁文章
扫一扫
2万+
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
