假设我有以下数据和命令：
clc;clear;
t = [0:0.1:1];
t_new = [0:0.01:1];
y = [1,2,1,3,2,2,4,5,6,1,0];
p = interp1(t,y,t_new,'spline');
plot(t,y,'o',t_new,p)
您可以看到它们工作得很好，从某种意义上说，插值函数可以很好地匹配节点处的数据点。但我的问题是，我需要计算 y（即 p 函数）wrt 时间的精确导数，并将其与 t 向量作图。怎么做到呢？我不会使用 diff 命令，因为我需要确保导数函数与 t 向量具有相同的长度。非常感谢。