目前为止,我们探讨了单变量/特征的回归模型,现在我们对房价模型增加更多的特征,例如房间数楼层等,构成一个含有多个变量的模型,模型中的特征为(
x
1
x_{1}
x1,
x
2
x_{2}
x2,…,
x
n
x_{n}
xn)。
增添更多特征后,我们引入一系列新的注释:
n
n
n 代表特征的数量
x i x^{i} xi代表第 i i i 个训练实例,是特征矩阵中的第 i i i 行,是一个向量(vector)。
比方说,上图的
x
j
(
i
)
x_{j}^{(i)}
xj(i)代表特征矩阵中第
i
i
i 行的第
j
j
j 个特征,也就是第
i
i
i 个训练实例的第
j
j
j 个特征。
如上图的,
x
2
(
2
)
x_{2}^{(2)}
x2(2)
=
3
=3
=3,
x
3
(
2
)
=
2
x_{3}^{(2)}=2
x3(2)=2
支持多变量的假设 表示为:,
h
θ
h_{\theta }
hθ(x)=
θ
0
\theta _{0}
θ0+
θ
1
\theta _{1}
θ1
x
1
x_{1}
x1+
θ
2
\theta _{2}
θ2
x
2
x_{2}
x2+…+
θ
n
\theta _{n}
θn
x
n
x_{n}
xn。
这个公式中有
n
n
n个参数和个变量,为了使得公式能够简化一些,引入
x
0
=
1
x_{0}=1
x0=1,则公式转化为:
h
θ
h_{\theta }
hθ(x)=
θ
0
x
0
\theta _{0}x_{0}
θ0x0+
θ
1
\theta _{1}
θ1
x
1
x_{1}
x1+
θ
2
\theta _{2}
θ2
x
2
x_{2}
x2+…+
θ
n
\theta _{n}
θn
x
n
x_{n}
xn
此时模型中的参数是一个维的向量,任何一个训练实例也都是
(
n
+
1
)
(n+1)
(n+1)维的向量,特征矩阵X的维度是
m
∗
(
n
+
1
)
m*(n+1)
m∗(n+1) 。 因此公式可以简化为:
h
θ
(
x
)
=
θ
T
X
h_{\theta }(x)=\theta ^{T}X
hθ(x)=θTX,其中上标代表矩阵转置。
多变量梯度下降
与单变量线性回归类似,在多变量线性回归中,我们也构建一个代价函数,则这个代价函数是所有建模误差的平方和,即: ,
其中:
h
θ
(
x
)
=
θ
T
X
h_{\theta }(x)=\theta ^{T}X
hθ(x)=θTX=
θ
0
x
0
\theta _{0}x_{0}
θ0x0+
θ
1
\theta _{1}
θ1
x
1
x_{1}
x1+
θ
2
\theta _{2}
θ2
x
2
x_{2}
x2+…+
θ
n
\theta _{n}
θn
x
n
x_{n}
xn,
我们的目标和单变量线性回归问题中一样,是要找出使得代价函数最小的一系列参数。 多变量线性回归的批量梯度下降算法为:
Python 代码:
def computeCost(X, y, theta):
inner = np.power(((X * theta.T) - y), 2)
return np.sum(inner) / (2 * len(X))
梯度下降法实践1-特征缩放
在我们面对多维特征问题的时候,我们要保证这些特征都具有相近的尺度,这将帮助梯度下降算法更快地收敛。
以房价问题为例,假设我们使用两个特征,房屋的尺寸和房间的数量,尺寸的值为 0-2000平方英尺,而房间数量的值则是0-5,以两个参数分别为横纵坐标,绘制代价函数的等高线图能,看出图像会显得很扁,梯度下降算法需要非常多次的迭代才能收敛。
解决的方法是尝试将所有特征的尺度都尽量缩放到-1到1之间。如图:
梯度下降法实践2-学习率
梯度下降算法收敛所需要的迭代次数根据模型的不同而不同,我们不能提前预知,我们可以绘制迭代次数和代价函数的图表来观测算法在何时趋于收敛。
也有一些自动测试是否收敛的方法,例如将代价函数的变化值与某个阀值(例如0.001)进行比较,但通常看上面这样的图表更好。
梯度下降算法的每次迭代受到学习率的影响,如果学习率
a
a
a过小,则达到收敛所需的迭代次数会非常高;如果学习率
a
a
a过大,每次迭代可能不会减小代价函数,可能会越过局部最小值导致无法收敛。
通常可以考虑尝试些学习率