用原帖作者的衰减曲线的excel文件，进行测试
python代码如下：
import math
import numpy as np
import xlrd
xls=xlrd.open_workbook(r'd:\data.xls')##读取excel文件
sheet=xls.sheet_by_name('衰减曲线')
y=[]
for i in range(1,14,1):
temp=[]
for j in range(0,60,1):
temp.append(sheet.cell(3+j,i).value)##共13组
y.append(temp)
#建立A矩阵
A=np.ones((2,60))
A[1]=-np.arange(0,60)/1000000
A=A.T
A=np.matrix(A)
AA=((A.T*A)**-1)*A.T
#求13组数据的时间
b=np.log(np.matrix(y).T)
b=AA*b
print(b[1])
输出结果：时间被放大过1000000
[[28068.88373367 27974.40821863 28652.51187232 28526.23801597
28072.3430224 28263.22917422 28453.26937705 28072.08262794
28143.77543767 28598.12668447 29059.58836214 29655.77554097
35143.3727371 ]]
| 
楼主
标题置顶
标题高亮
