问题1：IndexError: too many indices for array
答：在网上查阅了一圈后，发现造成这个问题的原因有很多。而我是由于用scipy.optimize.minimize()做优化时，传入的参数x0 = theta的shape为矩阵，而官方要求是向量的形式！（哎，太粗心啦）

问题2：ValueError: shapes (2,2) and (1,2) not aligned: 2 (dim 1) != 1 (dim 0)
答：这个问题是在问题1解决之后产生的新问题！疯狂的打印信息后，发现：如果用scipy.optimize.minimize()做优化，梯度函数gradient返回值必须是1维数组的形式，即（n,）。因为梯度的shape必须和参数theta的shape一致。