吴恩达机器学习系列内容的学习目录 → \rightarrow →吴恩达机器学习系列内容汇总。
机器学习中的正则化是一种为了减小测试误差的行为。我们在搭建机器学习模型时,最终目的是让模型在面对新数据的时候,可以有很好的表现。当用比较复杂的模型(比如神经网络)去拟合数据时,很容易出现过拟合现象,这会导致模型的泛化能力下降,这时候我们就需要使用正则化技术去降低模型的复杂度,从而改变模型的拟合度。
1. 过拟合的问题
正则化可以改善或者减少过度拟合的问题。
过拟合的问题就是指我们有非常多的特征,通过学习得到的模型能够非常好地适应训练集(代价函数可能几乎为0),但是推广到新的数据集上效果会非常的差。
图1是一个回归问题的例子:
图1中a)对应的模型是一个线性模型,欠拟合(Underfitting),不能很好地适应我们的训练集;图1中c)对应的模型是一个四次方的模型,过于强调拟合原始数据,而丢失了算法的本质——预测新数据。我们可以看出,若给出一个新的值让其预测,它将表现的很差,是过拟合(Overfitting),虽然能非常好地适应我们的训练集但在新输入变量进行预测时可能会效果不好;图1中b)对应的模型的模型似乎最合适("Just right "),所以其泛化能力最好。
分类问题中也存在这样的问题,如图2所示。
如果我们发现了过拟合问题,应该如何解决?
- 丢弃一些不能帮助我们正确预测的特征。可以是手工选择保留哪些特征,或者使用一些模型选择的算法来帮忙(例如PCA);
- 正则化。 保留所有的特征,但是减少参数的大小。
2. 代价函数
发生过拟合时,我们通过在代价函数中增加惩罚项来防止过拟合:
J
(
θ
)
=
1
2
m
[
∑
i
=
1
m
(
h
θ
(
x
(
i
)
)
−
y
(
i
)
)
2
+
λ
∑
j
=
1
n
θ
j
2
]
J(\theta)=\frac{1}{2m}[\sum_{i=1}^{m} (h _{\theta}(x^{(i)})-y^{(i)})^{2}+\lambda \sum_{j=1}^{n} \theta _{j}^{2}]
J(θ)=2m1[i=1∑m(hθ(x(i))−y(i))2+λj=1∑nθj2]
其中
λ
\lambda
λ又称为正则化参数(Regularization Parameter)。根据惯例,我们不对
θ
0
\theta _{0}
θ0项进行惩罚。
如果选择的正则化参数
λ
\lambda
λ过大,则会把所有的参数都最小化了,导致模型变成
h
θ
(
x
)
=
θ
0
h_{θ}(x)=\theta _{0}
hθ(x)=θ0,也就是图3中红色直线所示的情况,造成欠拟合。
那为什么增加一项
λ
∑
j
=
1
n
θ
j
2
\lambda \sum_{j=1}^{n} \theta _{j}^{2}
λ∑j=1nθj2可以使
θ
\theta
θ的值减小呢? 因为如果我们令
λ
\lambda
λ的值很大的话,为了使代价函数
J
(
θ
)
J(\theta)
J(θ)尽可能的小,所有的
θ
\theta
θ的值(不包括
θ
0
\theta _{0}
θ0)都会在一定程度上减小。 但若
λ
\lambda
λ的值太大了,那么
θ
\theta
θ的值(不包括
θ
0
\theta _{0}
θ0)都会趋近于0,这样我们所得到的只能是一条平行于轴的直线。 所以对于正则化,我们要取一个合理的
λ
\lambda
λ值,这样才能更好的应用正则化。
3. 线性回归的正则化
正则化线性回归的代价函数为 J ( θ ) = 1 2 m ∑ i = 1 m [ ( h θ ( x ( i ) ) − y ( i ) ) 2 + λ ∑ j = 1 n θ j 2 ] J(\theta)=\frac{1}{2m}\sum_{i=1}^{m} [(h _{\theta}(x^{(i)})-y^{(i)})^{2}+\lambda \sum_{j=1}^{n} \theta _{j}^{2}] J(θ)=2m1i=1∑m[(hθ(x(i))−y(i))2+λj=1∑nθj2]
对于线性回归在代价函数最小时参数 θ \theta θ的求解,我们之前推导了两种学习算法:一种基于梯度下降,一种基于正规方程。
- 梯度下降
R e p e a t { θ 0 : = θ 0 − α 1 m ∑ i = 1 m ( h θ ( x ( i ) ) − y ( i ) ) x 0 ( i ) θ j : = θ j − α [ 1 m ∑ i = 1 m ( h θ ( x ( i ) ) − y ( i ) ) x j ( i ) + λ m θ j ] ( f o r j = 1 , 2 , 3 , . . . , n ) } \begin{matrix}Repeat\ _{}\left \{ \right.\theta _{0}:=\theta _{0}-\alpha \frac{1}{m}\sum_{i=1}^{m} (h _{\theta}(x^{(i)})-y^{(i)})x_{0}^{(i)}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\\ _{}\ _{}\ _{}\ _{}\ \theta _{j}:=\theta _{j}-\alpha [\frac{1}{m}\sum_{i=1}^{m} (h _{\theta}(x^{(i)})-y^{(i)})x_{j}^{(i)}+\frac{\lambda }{m}\theta _{j}]\ _{}(for\ _{}j=1,2,3,...,n)\left. \right \} \end{matrix} Repeat {θ0:=θ0−αm1∑i=1m(hθ(x(i))−y(i))x0(i) θj:=θj−α[m1∑i=1m(hθ(x(i))−y(i))xj(i)+mλθj] (for j=1,2,3,...,n)}
对上面的算法中 j = 1 , 2 , 3 , . . . , n j=1,2,3,...,n j=1,2,3,...,n时的更新式子进行调整可得: θ j : = θ j ( 1 − α λ m ) − α 1 m ∑ i = 1 m ( h θ ( x ( i ) ) − y ( i ) ) x j ( i ) ( f o r j = 1 , 2 , 3 , . . . , n ) \theta _{j}:=\theta _{j}(1-\alpha\frac{\lambda }{m})-\alpha \frac{1}{m}\sum_{i=1}^{m} (h _{\theta}(x^{(i)})-y^{(i)})x_{j}^{(i)}\ _{}\ _{}(for\ _{}\ _{}j=1,2,3,...,n) θj:=θj(1−αmλ)−αm1i=1∑m(hθ(x(i))−y(i))xj(i) (for j=1,2,3,...,n)其中, 1 − α λ m 1-\alpha\frac{\lambda }{m} 1−αmλ是一个略小于1的数。可以看出,当我们进行正则化线性回归时,我们要做的就是每次迭代时都将 θ j \theta _{j} θj乘以一个略小于1的数,每次都把参数缩小一点,然后进行和之前一样的更新操作。
- 正规方程
使用正规方程来求解正则化线性回归模型: θ = ( X T X + λ [ 0 1 1 ⋱ 1 ] ) − 1 X T y \theta =(X^{T}X+\lambda \begin{bmatrix} 0 & & & & \\ & 1 & & & \\ & & 1& & \\ & & & \ddots & \\ & & & &1 \end{bmatrix})^{-1}X^{T}y θ=(XTX+λ⎣⎢⎢⎢⎢⎡011⋱1⎦⎥⎥⎥⎥⎤)−1XTy
其中, X = [ ( x ( 1 ) ) T ⋯ ( x ( m ) ) T ] X=\begin{bmatrix} (x^{(1)})^{T}\\ \cdots \\ (x^{(m)})^{T} \end{bmatrix} X=⎣⎡(x(1))T⋯(x(m))T⎦⎤, y = [ y ( 1 ) ⋯ y ( m ) ] y=\begin{bmatrix} y^{(1)}\\ \cdots \\ y^{(m)} \end{bmatrix} y=⎣⎡y(1)⋯y(m)⎦⎤。
4. Logistic回归的正则化
正则化Logistic回归的代价函数为 J ( θ ) = 1 m ∑ i = 1 m [ − y ( i ) × l o g ( h θ ( x ( i ) ) ) − ( 1 − y ( i ) ) × l o g ( 1 − h θ ( x ( i ) ) ) ] + λ 2 m ∑ j = 1 n θ j 2 J(\theta)=\frac{1}{m}\sum_{i=1}^{m} [-y^{(i)}\times log(h_{\theta}(x^{(i)}))-(1-y^{(i)})\times log(1-h_{\theta}(x^{(i)}))]+\frac{\lambda }{2m} \sum_{j=1}^{n} \theta _{j}^{2} J(θ)=m1i=1∑m[−y(i)×log(hθ(x(i)))−(1−y(i))×log(1−hθ(x(i)))]+2mλj=1∑nθj2
对于Logistic回归在代价函数最小时参数 θ \theta θ的求解,我们之前学习过两种优化算法:首先学习了使用梯度下降法来优化代价函数 J ( θ ) J(\theta) J(θ),然后学习了更高级的优化算法,这些高级优化算法需要我们自己设计代价函数 J ( θ ) J(\theta) J(θ)。
- 梯度下降
R e p e a t { θ 0 : = θ 0 − α 1 m ∑ i = 1 m ( h θ ( x ( i ) ) − y ( i ) ) x 0 ( i ) θ j : = θ j − α [ 1 m ∑ i = 1 m ( h θ ( x ( i ) ) − y ( i ) ) x j ( i ) + λ m θ j ] ( f o r j = 1 , 2 , 3 , . . . , n ) } \begin{matrix}Repeat\ _{}\left \{ \right.\theta _{0}:=\theta _{0}-\alpha \frac{1}{m}\sum_{i=1}^{m} (h _{\theta}(x^{(i)})-y^{(i)})x_{0}^{(i)}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\ _{}\\ _{}\ _{}\ _{}\ _{}\ \theta _{j}:=\theta _{j}-\alpha [\frac{1}{m}\sum_{i=1}^{m} (h _{\theta}(x^{(i)})-y^{(i)})x_{j}^{(i)}+\frac{\lambda }{m}\theta _{j}]\ _{}(for\ _{}j=1,2,3,...,n)\left. \right \} \end{matrix} Repeat {θ0:=θ0−αm1∑i=1m(hθ(x(i))−y(i))x0(i) θj:=θj−α[m1∑i=1m(hθ(x(i))−y(i))xj(i)+mλθj] (for j=1,2,3,...,n)}
虽然其看上去同线性回归一样,但是Logistic回归的 h θ ( x ) = g ( θ T X ) h_{\theta }(x)=g(\theta^TX) hθ(x)=g(θTX) ,所以与线性回归不同。
- 高级优化算法
需要自己设计代价函数 J ( θ ) J(\theta) J(θ)。