求拟合后曲线上某点斜率

huangzhiye

三次样条插值生成曲线，并求得曲线某点处的斜率，以向量形式返回结果。
数据如下：

x1=[-3,-2,-1,0]

y1=[1.10277661488344,1.33243845991561,1.54170402328429,1.70243053644553]

sp=csapi(x1,y1)    %插值生成三次样条曲线函数

fx= differentiate(sp,x1);

首先拟合成三次样条曲线，然后求得 x1=[-3,-2,-1,0]  4个点处的曲线斜率，以向量形式返回所求4个点的斜率值。

为何出错？

谢谢！

nwcwww

csapi生成的是ppform，而differentiate针对的是cfit或者sfit对象，不接受struct.
这里求斜率要用fnval(fnder(sp), x1)

Nicky_小牛

楼上正解~不愧是版主啊。若是可导的函数，可以利用LZ的这条语句：fx= differentiate(sp,x1);% 求符号表达式的微分，但是对于样条函数来说，不合适~

huangzhiye

nwcwww 发表于 2014-5-23 02:25
csapi生成的是ppform，而differentiate针对的是cfit或者sfit对象，不接受struct.
这里求斜率要用fnval(fnde ...

明白了，谢谢 nwcwww版主的热心解答！