拟合线性函数的几种方法
A = np.vander(x, 2) #生成范德蒙矩阵(数据,输出列数)C = np.diag(yerr * yerr) ATA = np.dot(A.T, A / (yerr ** 2)[:, None])cov = np.linalg.inv(ATA)w = np.linalg.solve(ATA, np.dot(A.T, y / yerr ** 2))np.diag(array) 中,array是一个1维数组时,结果形成一个以一维数组为对角线元素的矩阵,array是一个二维矩阵时,结果输
原创
2020-12-22 16:16:50 ·
4870 阅读 ·
1 评论