[½»Á÷] ÓÃorigin»ò1stoptÄâºÏÒ»ÌõÒþº¯Êý·½³ÌÇúÏß

#!/usr/bin/env python
import scipy.optimize as so
import numpy as np
import pylab as pl
x=np.loadtxt('te.txt')[:,0]
y=np.loadtxt('te.txt')[:,1]
def fun(p):
    return np.arcsinh(p[1]*x**p[2])/p[0]-y
pout=so.leastsq(fun,[1.0,1.0,1.0],maxfev=100000)[0]
print 'a,b,c is ',str(pout)
pl.plot(x,y)
x0=np.linspace(x.min(),x.max(),100)
pl.plot(x0,np.arcsinh(pout[1]*x0**pout[2])/pout[0])
pl.savefig('te.png')