矩阵的平方和

Question

我确实有一个以行（不同pH值下的测量值）和数据点为列（随时间的浓度）的观察值的矩阵。所以一行包含一个pH值的不同数据点

我确实想为数据配一首颂歌。所以我定义了一个成本函数，想计算所有观测值的平方和。取此矩阵的平方和应如下：

res = y - yhat                        # calculate residuals
ssq = np.diag(np.dot(res.T,res))      # sum over the diagonal

是这样吗？

Answer 1

如果取最后一个数组的和，那么它是正确的。 但它也不必要地复杂（因为非对角元素也是用np.dot计算的） 更快的是：

ssq = np.sum(res**2)

如果您希望每个实验都使用ssd，可以执行以下操作：

ssq = np.sum(res**2, axis=1)