{OFF}如何在matlab上推导样条插值函数？

Question

我有电池放电曲线的实验室测试数据。数据集中在电压随时间变化的22个点上。在matlab中，我用样条插值法绘制了一条插值曲线，但是我想做这个图的导数，我怎么做呢？

在这里输入图像描述

代码

x = [0; 3600    ;7200   ;10800; 14400;  18000;  21600;  25200   ;28800;...
32400;  36000   ;39600  ;43200  ;46800; 50400   ;54000; 57600;  61200;...
64800   ;68400  ;72000; 74880];
y = [12.75; 12.40;  12.38;  12.34;  12.30;  12.26   ;12.21  ;12.17  ;...
 12.12; 12.07;  12.02   ;11.97  ;11.91  ;11.85; 11.79;  11.72;  11.65;...
 11.56  ;11.46  ;11.35  ;11.17; 10.59];

f = fit( x, y,'cubicinterp')

Answer 1

您可以使用gradient of y/x

f = fit( x, y,'cubicinterp')
df = gradient(f(x)); % f'x
dx = gradient(x);    % dx

dfx = df ./ dx;
plot(x, y, x, dfx);