阅读背景:

python收益率计算-scipy.optimize-大量计算的好帮手

来源:互联网 
import time
import datetime
from scipy.optimize import fsolve
for i in range(1,50):
    print('第'+str(i)+'次计算开始')
    x = Symbol('x')
    startTime=time.strftime('%Y-%m-%d %H:%M:%S',time.localtime())
    print(startTime)
    mystr=''
    for j in range(1,i+1):
        mystr+='+x0**'+str(j)
    mystr=mystr[1:]
    mystr+='-'+str(i*1.1)
    def f(x): 
        x0,x1 = x.tolist() 
        return [eval(mystr),x1-1]

    # f计算方程组的误差,[1,1,1]是未知数的初始值
    result = fsolve(f, [1,1]) 
    print('计算的结果:'+str(result[0]))
    endTime=time.strftime('%Y-%m-%d %H:%M:%S',time.localtime())
    print(endTime)
    startTime= datetime.datetime.strptime(startTime,"%Y-%m-%d %H:%M:%S")  
    endTime= datetime.datetime.strptime(endTime,"%Y-%m-%d %H:%M:%S")
    seconds = (endTime- startTime).seconds  
    print('第{}次方方程计算耗费的时间为{}秒'.format(i,seconds))from scipy.optimi



你的当前访问异常,请进行认证后继续阅读剩余内容。

分享到: