在我尝试重新发明任何轮子之前,有人已经使用numpy.poly1d
(或类似物)制作了用于分段多项式的Python类吗?它应该采用不连续性列表和多项式列表,并允许poly1d
的所有操作,或者至少(为了我的目的)添加,减去,乘法和积分;我需要这些,因为我的意思是使用分段二次方作为最小二乘拟合的基函数(对于分段连续输入而不是离散输入)。
(numpy.piecewise
不是它:它不会创建持久对象。)
更新:哎呀,我在2010年11月以某种方式获得了piecepoly.py
。(没有作者身份的迹象。代码与我今天写的非常相似,但有一条评论“The函数被认为是左连续的,“这不是我要说的。”