Matlab 中的三次样条实现

Cubic spline implementation in Matlab

我有两个向量 x 和 y 定义(随机)如下: x1=[1000 3000 5000 6000 4000 2000 500 0 -1000 -3000 -5000 -6000 -4000 -2000 -500 1 999 2999 4999]; y1=[5000 4999 4990 3500 2500 2499 2498 2497 2496 2495 2494 1000 -1000 -999 -998 -997 -996 -995 -994];

以下是我通过简单地输入 plot(x1,y1):

得到的图

有没有办法使用 interp1 命令从上述数据中生成平滑曲线?有人告诉我应该使用三次样条曲线来实现所需的绘图,但是,由于我是 Matlab 的新手,所以我不知道如何实现这样的解决方案。提前致谢!

编辑: 我试图按如下方式实现它,但我得到了一个看起来很可怕的情节! x1_temp=-6000:100:6000; pc=pchip(x1,y1,x1_temp); plot(x1,y1,'o',x1_temp,pc,'-');

我应该如何修改此代码以生成正确的图?

我认为您对插值的内容感到困惑。您应该分别对 x1 和 y1 进行插值,然后将它们相互绘制。以下示例生成平滑曲线:

x1=[1000 3000 5000 6000 4000 2000 500 0   -1000 -3000 -5000 -6000 -4000 -2000 -500 1   999 2999 4999];
y1=[5000 4999 4990 3500 2500 2499 2498 2497 2496 2495  2494  1000  -1000 -999 -998 -997 -996 -995 -994];

s = [0,cumsum(sqrt(diff(x1).^2+diff(y1).^2))]
N = length(s);

figure();
plot(x1,y1);
hold on

s_fine = interp1(linspace(0,1,N),s,linspace(0,1,5*N));
pcx=interp1(s,x1,s_fine,'spline');
pcy=interp1(s,y1,s_fine,'spline');
plot(pcx,pcy,'r-');