用原帖作者的衰减曲线的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 ]]
| 
楼主
