python – 使用类的__call__方法作为Numpy curve_fit的输入

内容导读

互联网集市收集整理的这篇技术教程文章主要介绍了python – 使用类的__call__方法作为Numpy curve_fit的输入，小编现在分享给大家，供广大互联网技能从业者学习和参考。文章包含1713字，纯文字阅读大概需要3分钟。

内容图文

我想使用类的__call__方法作为Numpy curve_fit函数的输入,因为我的功能和数据准备过程非常复杂(将分析模型数据拟合到某些测量中).通过定义一个函数它可以正常工作,但我不能让它与类一起工作.

要重新创建我的问题,您可以运行：

import numpy as np
from scipy.optimize import curve_fit

#WORKS:
#def goal(x,a1,a2,a3,a4,a5):
#    y=a1*x**4*np.sin(x)+a2*x**3+a3*x**2+a4*x+a5
#    return y

# DOES NOT WORK:
class func():
    def __call__(self,x,a1,a2,a3,a4,a5):
        y=a1*x**4*np.sin(x)+a2*x**3+a3*x**2+a4*x+a5
        return y    

goal=func()

#data prepraration ***********
xdata=np.linspace(0,50,100)
ydata=goal(xdata,-2.1,-3.5,6.6,-1,2)
# ****************************

popt, pcov = curve_fit(goal, xdata, ydata)
print 'optimial parameters',popt
print 'The estimated covariance of optimial parameters',pcov

我得到的错误是：

Traceback (most recent call last):
  File "D:\...some path...\test_minimizacija.py", line 35, in <module>
    popt, pcov = curve_fit(goal, xdata, ydata)
  File "C:\Python26\lib\site-packages\scipy\optimize\minpack.py", line 412, in curve_fit
    args, varargs, varkw, defaults = inspect.getargspec(f)
  File "C:\Python26\lib\inspect.py", line 803, in getargspec
    raise TypeError('arg is not a Python function')
TypeError: arg is not a Python function

我怎样才能做到这一点？

解决方法:

容易(虽然,不漂亮),只需将其更改为：

popt, pcov = curve_fit(goal.__call__, xdata, ydata)

有趣的是,numpy强制你将一个函数对象传递给curve_fit而不是一个任意的可调用…

快速检查curve_fit的来源,似乎另一种解决方法可能是：

popt,pcov = curve_fit(goal, xdata, ydata, p0=[1]*5)

这里,p0是拟合参数的初始猜测.问题似乎是scipy检查函数的参数,以便它知道如果你实际上没有提供参数作为初始猜测使用多少参数.在这里,由于我们有5个参数,我的初始猜测是长度为5的所有参数的列表.(如果你也没有提供猜测,则scipy默认使用1).

内容总结

以上是互联网集市为您收集整理的python – 使用类的__call__方法作为Numpy curve_fit的输入全部内容，希望文章能够帮你解决python – 使用类的__call__方法作为Numpy curve_fit的输入所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错，欢迎将互联网集市网站推荐给程序员好友。

内容备注

版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 gblab@vip.qq.com 举报，一经查实，本站将立刻删除。

内容手机端