坐标轴下降法
坐标轴下降法顾名思义,沿着坐标轴下降。坐标轴下降法和梯度下降法使用的都是迭代法,即使用启发式的方式一步一步迭代求解函数的最小值。
可以想象一个可微的凸函数 J ( ω ) J(\omega) J(ω),其中 ω \omega ω是一个 n ∗ 1 n*1 n∗1维的向量。如果在这 n n n维空间中存在着某一点 ω ‾ \overline{\omega} ω能够使得 J ( ω ) J(\omega) J(ω)在每一个坐标轴上都是最小的,那么我们则可以说 J ( ω ‾ ) J(\overline{\omega}) J(ω)是全局最小值。如二维空间 ( x , y ) (x,y) (x,y)中的曲线 x 2 x^2 x2,当 x x x和 y y y都是最小的时候,即为曲线 x 2 x^2 x2最小值的时候。(注:Lasso回归的目标函数虽然不可微,但是由于目标函数非正则项是可微的,L1正则项是凸函数,因此也可以使用该结论。)
坐标轴下降法流程
坐标轴下降法即在 ω \omega ω的每个坐标轴上做迭代下降,当每个坐标轴上的值都收敛时,则可以说明达到了 J ( ω ) J(\omega) J(ω)的全局最小值。
- ω \omega ω取初值记为 ω ( 0 ) \omega^{(0)} ω(0),其中 ( 0 ) (0) (0)表示迭代的轮数,初轮数为0。
- 对于第
k
k
k轮的迭代,依次求解
ω
i
(
k
)
(
i
=
1
,
2
,
⋯
,
n
)
{\omega_i}^{(k)} \quad(i=1,2,\cdots,n)
ωi(k)(i=1,2,⋯,n),即
ω i ( k ) ∈ a r g m i n ⏟ ω i J ( ω 1 ( k ) , ω 2 ( k ) , ⋯ , ω n ( k ) ) {\omega_i}^{(k)} \in \underbrace{argmin}_{\omega_i}\,J({\omega_1}^{(k)},{\omega_2}^{(k)},\cdots,{\omega_n}^{(k)}) ωi(k)∈ωi argminJ(ω1(k),ω2(k),⋯,ωn(k))
其中由于在每次单独对一个坐标轴迭代的时候其他坐标轴上的值为常数。也就是 ω i ( k ) {\omega_i}^{(k)} ωi(k)是使 J ( ω 1 ( k ) , ω i − 1 ( k ) , ω i + 1 ( k ) , ⋯ , ω n ( k ) ) J({\omega_1}^{(k)},{\omega_i-1}^{(k)},{\omega_i+1}^{(k)},\cdots,{\omega_n}^{(k)}) J(ω1(k),ωi−1(k),ωi+1(k),⋯,ωn(k))为最小的值,除了 ω i ( k ) {\omega_i}^{(k)} ωi(k)为变量外,其他参数都为常数,则对 J ( ω 1 ( k ) , ω i − 1 ( k ) , ω i + 1 ( k ) , ⋯ , ω n ( k ) ) J({\omega_1}^{(k)},{\omega_i-1}^{(k)},{\omega_i+1}^{(k)},\cdots,{\omega_n}^{(k)}) J(ω1(k),ωi−1(k),ωi+1(k),⋯,ωn(k))求导即可得该函数的最小值。 - 检查 ω ( k ) {\omega}^{(k)} ω(k)和 ω ( k − 1 ) {\omega}^{(k-1)} ω(k−1)在每个维度上的变化情况,如果所有维度的变换都小于某个阈值,则 ω ( k ) {\omega}^{(k)} ω(k)为最终结果,否则循环第2步,进入 k + 1 k+1 k+1次迭代。
坐标轴下降法和梯度下降法的异同
- 两者都是迭代方法,并且每一轮迭代都需要 O ( m n ) O(mn) O(mn)的计算量,其中 m m m为样本数, n n n为特征数。
- 坐标轴下降法固定 n − 1 n-1 n−1个坐标,沿着剩下的一个坐标轴迭代,得到函数的局部最小值;梯度下降法则是沿着梯度的负方向得到函数的局部最小值。
- 坐标轴下降法是一种非梯度优化算法,但是一个周期中循环使用不同的坐标方向迭代,相当于梯度下降的一轮迭代。
- 由于坐标轴下降法是利用当前坐标方向迭代,不需要求目标函数的倒数;梯度下降法是利用目标函数的偏导数来确定迭代方向的。