python怎么编程序求分段函数_在pandas/numpy中实现分段函数的正确方法

我需要创建一个传递给curve_fit的函数。在我的例子中,函数最好定义为分段函数。在

我知道以下命令不起作用,但我展示它是因为它使函数的意图变得清晰:def model_a(X, x1, x2, m1, b1, m2, b2):

'''f(x) has form m1*x + b below x1, m2*x + b2 above x2, and is

a cubic spline between those two points.'''

y1 = m1 * X + b1

y2 = m2 * X + b2

if X <= x1:

return y1 # function is linear below x1

if X >= x2:

return y2 # function is linear above x2

# use a cubic spline to interpolate between lower

# and upper line segment

a, b, c, d = fit_cubic(x1, y1, x2, y2, m1, m2)

return cubic(X, a, b, c, d)

当然,问题是X是一个pandas系列,而形式(X <= x1)的计算结果是一系列布尔值,因此这将失败,并显示消息“序列的真值是模棱两可的”

似乎^{}正是为这种情况而设计的:“只要condlist[i]为真,funclist[i](x)就被用作输出值。”所以我尝试了以下方法:

^{pr2}$

但这次失败了:return np.piecewise(X, [X<=x1, X>=x2], [lo, hi, mid])

“索引器错误:数组的索引太多”。我倾向于认为它反对这样一个事实,即在condlist中有两个元素,在funclist中有三个元素,但是文档明确指出funclist中的额外元素被视为默认值。在

有什么指引吗?在

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值