问题描述
我使用 scipy 三次样条(“scipy.interpolate.CubicSpline”)进行一维插值。 我想为样条的每个段指定初始和最终导数(“边界条件”,“bc_type”,如果我没记错的话)。 但是,阅读文档并寻找示例时,我无法理解这将如何发生(我假设这是可能的)
例如,如果我有 4 个要插值的点,我希望有 3 个段,对于每个段,我应该允许初始和结局一阶导数。但是,唯一接受的语法如下所示:
CubicSpline(x,y,bc_type=((1,<initial_derivative>),(1,<final_derivative>)))
其中指定了单个初始和最终导数(顺便说一句,对于哪个段?)
我尝试用我想到的任何值替换 bc_type
值,例如一对元组的列表,每个点一个,但不接受与此处报告的语法不同的语法。
此时我开始怀疑这是可能的,但它绝对应该:三次多项式表达式有 4 个自由度,允许“固定”两个值和两个导数
任何信息将不胜感激!
解决方法
segments,对于每一个我应该允许初始和最后一阶导数
这是 CubicHermiteSpline
或 BPoly.from_derivatives
。 CubicSpline
修正内部斜率,使插值在内部节点处可微分两次。