內容列表：

• 1. 簡介
• 2. 損失函數可視化
• 3. 最優化
  • 3.1 策略#1：隨機搜索
  • 3.2 策略#2：隨機局部搜索
  • 3.3 策略#3：跟隨梯度
• 4. 梯度計算
  • 4.1 使用有限差值進行數值計算
  • 4.2 微分計算梯度
• 5. 梯度下降
• 6. 小結

1. 簡介

在上一節中，我們介紹了圖像分類任務中的兩個關鍵部分：

1. 基于參數的評分函數。該函數將原始圖像像素映射為分類評分值（例如：一個線性函數）。
2. 損失函數。該函數能夠根據分類評分和訓練集圖像數據實際分類的一致性，衡量某個具體參數集的質量好壞。損失函數有多種版本和不同的實現方式（例如：Softmax或SVM）。

上節中，線性函數的形式是

而SVM實現的公式是：

對于圖像數據 x_i ，如果基于參數集 W 做出的分類預測與真實情況比較一致，那么計算出來的損失值 L 就很低?，F在介紹第三個，也是最后一個關鍵部分：最優化Optimization。最優化是尋找能使得損失函數值最小化的參數 W 的過程。

鋪墊：一旦理解了這三個部分是如何相互運作的，我們將會回到第一個部分（基于參數的函數映射），然后將其拓展為一個遠比線性函數復雜的函數：首先是神經網絡，然后是卷積神經網絡。而損失函數和最優化過程這兩個部分將會保持相對穩定。

2. 損失函數可視化

本課中討論的損失函數一般都是定義在高維度的空間中（比如，在CIFAR-10中一個線性分類器的權重矩陣大小是[10x3073]，就有30730個參數），這樣要將其可視化就很困難。然而辦法還是有的，在1個維度或者2個維度的方向上對高維空間進行切片，就能得到一些直觀感受。例如，隨機生成一個權重矩陣W ，該矩陣就與高維空間中的一個點對應。然后沿著某個維度方向前進的同時記錄損失函數值的變化。換句話說，就是生成一個隨機的方向 W₁ 并且沿著此方向計算損失值，計算方法是根據不同的 a 值來計算 L( W + aW₁ ) 。這個過程將生成一個圖表，其x軸是 a 值，y軸是損失函數值。同樣的方法還可以用在兩個維度上，通過改變 a, b 來計算損失值 L( W + aW₁ + bW₂ ) ，從而給出二維的圖像。在圖像中，a, b 可以分別用x和y軸表示，而損失函數的值可以用顏色變化表示：

一個無正則化的多類SVM的損失函數的圖示。左邊和中間只有一個樣本數據，右邊是CIFAR-10中的100個數據。左：a值變化在某個維度方向上對應的的損失值變化。中和右：兩個維度方向上的損失值切片圖，藍色部分是低損失值區域，紅色部分是高損失值區域。注意損失函數的分段線性結構( 往下看有真相 ~ )。多個樣本的損失值是總體的平均值，所以右邊的碗狀結構是很多的分段線性結構的平均（比如中間這個就是其中之一）。

損失函數的分段線性結構。對于一個單獨的數據，有損失函數的計算公式如下：

通過公式可見，每個樣本的數據損失值是以 W 為參數的線性函數的總和（零閾值來源于 max( 0，- ) 函數）。 W 的每一行（即 W_j ），有時候它前面是一個正號（比如當它對應錯誤分類的時候），有時候它前面是一個負號（比如當它是是正確分類的時候）。為進一步闡明，假設有一個簡單的數據集，其中包含有3個只有1個維度的點，數據集數據點有3個類別。那么完整的無正則化SVM的損失值計算如下：

因為這些例子都是一維的，所以數據 x_i 和權重 W_j 都是數字。觀察 W₀ ，可以看到上面的式子中一些項是 W₀ 的線性函數，且每一項都會與0比較，取兩者的最大值。可作圖如下：

從一個維度方向上對數據損失值的展示。x軸方向就是一個權重，y軸就是損失值。數據損失是多個部分組合而成。其中每個部分要么是某個權重的獨立部分，要么是該權重的線性函數與0閾值的比較。完整的SVM數據損失就是這個形狀的30730維版本。( 反復看這段話，感覺還是有點道理的。。。)

需要多說一句的是，你可能根據SVM的損失函數的碗狀外觀猜出它是一個凸函數。關于如何高效地最小化凸函數的論文有很多，你也可以學習斯坦福大學關于（凸函數最優化）的課程。但是一旦我們將 f 函數擴展到神經網絡，目標函數就就不再是凸函數了，圖像也不會像上面那樣是個碗狀，而是凹凸不平的復雜地形形狀。

不可導的損失函數。作為一個技術筆記，你要注意到：由于max操作，損失函數中存在一些不可導點（kinks），這些點使得損失函數不可微，因為在這些不可導點，梯度是沒有定義的。但是次梯度（subgradient）依然存在且常常被使用。在本課中，我們將交換使用次梯度和梯度兩個術語。

3. 最優化

重申一下：損失函數可以量化某個具體權重集 W 的質量。而最優化的目標就是找到能夠最小化損失函數值的 W 。我們現在就朝著這個目標前進，實現一個能夠最優化損失函數的方法。對于有一些經驗的同學，這節課看起來有點奇怪，因為使用的例子（SVM 損失函數）是一個凸函數問題。但是要記得，最終的目標是不僅僅對凸函數做最優化，而是能夠最優化一個神經網絡，而對于神經網絡是不能簡單的使用凸函數的最優化技巧的。

3.1 策略#1：隨機搜索（差勁）

既然確認參數集W的好壞蠻簡單的，那第一個想到的（差勁）方法，就是可以隨機嘗試很多不同的權重，然后看其中哪個最好。過程如下：

# 假設X_train的每一列都是一個數據樣本（比如3073 x 50000）
# 假設Y_train是數據樣本的類別標簽（比如一個長50000的一維數組）
# 假設函數L對損失函數進行評價

bestloss = float("inf") # Python assigns the highest possible float value
for num in xrange(1000):
  W = np.random.randn(10, 3073) * 0.0001 # generate random parameters
  loss = L(X_train, Y_train, W) # get the loss over the entire training set
  if loss < bestloss: # keep track of the best solution
    bestloss = loss
    bestW = W
  print 'in attempt %d the loss was %f, best %f' % (num, loss, bestloss)

# 輸出:
# in attempt 0 the loss was 9.401632, best 9.401632
# in attempt 1 the loss was 8.959668, best 8.959668
# in attempt 2 the loss was 9.044034, best 8.959668
# in attempt 3 the loss was 9.278948, best 8.959668
# in attempt 4 the loss was 8.857370, best 8.857370
# in attempt 5 the loss was 8.943151, best 8.857370
# in attempt 6 the loss was 8.605604, best 8.605604
# ... (trunctated: continues for 1000 lines)

在上面的代碼中，我們嘗試了若干隨機生成的權重矩陣W，其中某些的損失值較小，而另一些的損失值大些。我們可以把這次隨機搜索中找到的最好的權重W取出，然后去跑測試集：

# 假設X_test尺寸是[3073 x 10000], Y_test尺寸是[10000 x 1]
scores = Wbest.dot(Xte_cols) # 10 x 10000, the class scores for all test examples
# 找到在每列中評分值最大的索引（即預測的分類）
Yte_predict = np.argmax(scores, axis = 0)
# 以及計算準確率
np.mean(Yte_predict == Yte)
# 返回 0.1555

驗證集上表現最好的權重W跑測試集的準確率是15.5%，而完全隨機猜的準確率是10%，如此看來，這個準確率對于這樣一個不經過大腦的策略來說，還算不錯嘛！

核心思路：迭代優化。當然，我們肯定能做得更好些。核心思路是：雖然找到最優的權重W非常困難，甚至是不可能的（尤其當W中存的是整個神經網絡的權重的時候），但如果問題轉化為：對一個權重矩陣集W取優，使其損失值稍微減少。那么問題的難度就大大降低了。換句話說，我們的方法從一個隨機的W開始，然后對其迭代取優，每次都讓它的損失值變得更小一點。

我們的策略是從隨機權重開始，然后迭代取優，從而獲得更低的損失值。

3.2 策略#2：隨機局部搜索

第一個策略可以看做是每走一步都嘗試幾個隨機方向，如果某個方向是向山下的，就向該方向走一步。這次我們從一個隨機 W 開始，然后生成一個隨機的擾動 δW ，只有當 W + δW 的損失值變低，我們才會更新。這個過程的具體代碼如下：

W = np.random.randn(10, 3073) * 0.001 # 生成隨機初始W
bestloss = float("inf")
for i in xrange(1000):
  step_size = 0.0001
  Wtry = W + np.random.randn(10, 3073) * step_size
  loss = L(Xtr_cols, Ytr, Wtry)
  if loss < bestloss:
    W = Wtry
    bestloss = loss
  print 'iter %d loss is %f' % (i, bestloss)

使用同樣的數據（1000），這個方法可以得到21.4%的分類準確率。這個比策略一好，但是依然過于浪費計算資源。

3.3 跟隨梯度

梯度下降大家都很熟悉了，我就直接搬運原作者的話了，大家看一看就好

前兩個策略中，我們是嘗試在權重空間中找到一個方向，沿著該方向能降低損失函數的損失值。其實不需要隨機尋找方向，因為可以直接計算出最好的方向，這就是從數學上計算出最陡峭的方向。這個方向就是損失函數的梯度（gradient）。在蒙眼徒步者的比喻中，這個方法就好比是感受我們腳下山體的傾斜程度，然后向著最陡峭的下降方向下山。

在一維函數中，斜率是函數在某一點的瞬時變化率。梯度是函數的斜率的一般化表達，它不是一個值，而是一個向量。在輸入空間中，梯度是各個維度的斜率組成的向量（或者稱為導數derivatives）。對一維函數的求導公式如下：

當函數有多個參數的時候，我們稱導數為偏導數。而梯度就是在每個維度上偏導數所形成的向量。

4. 梯度計算

計算梯度有兩種方法：一個是緩慢的近似方法（數值梯度法），但實現相對簡單。另一個方法（分析梯度法）計算迅速，結果精確，但是實現時容易出錯，且需要使用微分?，F在對兩種方法進行介紹：

4.1 利用有限差值計算梯度（數值）

公式已經給出數值計算梯度的方法。下面代碼是一個輸入為函數f和向量x，計算f的梯度的通用函數，它返回函數f在點x處的梯度（也就是我們常說的數值方法算出來，簡單粗暴）：

def eval_numerical_gradient(f, x):
  """  
  一個f在x處的數值梯度法的簡單實現
  - f是只有一個參數的函數
  - x是計算梯度的點
  """ 

  fx = f(x) # 在原點計算函數值
  grad = np.zeros(x.shape)
  h = 0.00001

  # 對x中所有的索引進行迭代
  it = np.nditer(x, flags=['multi_index'], op_flags=['readwrite'])
  while not it.finished:

    # 計算x+h處的函數值
    ix = it.multi_index
    old_value = x[ix]
    x[ix] = old_value + h # 增加h
    fxh = f(x) # 計算f(x + h)
    x[ix] = old_value # 存到前一個值中 (非常重要)

    # 計算偏導數
    grad[ix] = (fxh - fx) / h # 坡度
    it.iternext() # 到下個維度

  return grad

根據上面的梯度公式，代碼對所有維度進行迭代，在每個維度上產生一個很小的變化h，通過觀察函數值變化，計算函數在該維度上的偏導數。最后，所有的梯度存儲在變量grad中。

實踐考量：注意在數學公式中，h的取值是趨近于0的，然而在實際中，用一個很小的數值（比如例子中的1e-5）就足夠了。在不產生數值計算出錯的理想前提下，你會使用盡可能小的h。還有，實際中用中心差值公式（centered difference formula）[f(x + h) - f(x - h)] / 2h 效果較好。

可以使用上面這個公式來計算任意函數在任意點上的梯度。下面計算權重空間中的某些隨機點上，CIFAR-10損失函數的梯度：

# 要使用上面的代碼我們需要一個只有一個參數的函數
# (在這里參數就是權重)所以也包含了X_train和Y_train
def CIFAR10_loss_fun(W):
  return L(X_train, Y_train, W)

W = np.random.rand(10, 3073) * 0.001 # 隨機權重向量
df = eval_numerical_gradient(CIFAR10_loss_fun, W) # 得到梯度

梯度告訴我們損失函數在每個維度上的斜率，以此來進行更新：

loss_original = CIFAR10_loss_fun(W) # 初始損失值
print 'original loss: %f' % (loss_original, )

# 查看不同步長的效果
for step_size_log in [-10, -9, -8, -7, -6, -5,-4,-3,-2,-1]:
  step_size = 10 ** step_size_log
  W_new = W - step_size * df # 權重空間中的新位置
  loss_new = CIFAR10_loss_fun(W_new)
  print 'for step size %f new loss: %f' % (step_size, loss_new)

# 輸出:
# original loss: 2.200718
# for step size 1.000000e-10 new loss: 2.200652
# for step size 1.000000e-09 new loss: 2.200057
# for step size 1.000000e-08 new loss: 2.194116
# for step size 1.000000e-07 new loss: 2.135493
# for step size 1.000000e-06 new loss: 1.647802
# for step size 1.000000e-05 new loss: 2.844355
# for step size 1.000000e-04 new loss: 25.558142
# for step size 1.000000e-03 new loss: 254.086573
# for step size 1.000000e-02 new loss: 2539.370888
# for step size 1.000000e-01 new loss: 25392.214036

在梯度負方向上更新：在上面的代碼中，為了計算W_new，要注意我們是向著梯度df的負方向去更新，這是因為我們希望損失函數值是降低而不是升高。

步長的影響：梯度指明了函數在哪個方向是變化率最大的，但是沒有指明在這個方向上應該走多遠。在后續的課程中可以看到，選擇步長（也叫作學習率）將會是神經網絡訓練中最重要（也是最頭痛）的超參數設定之一。還是用蒙眼徒步者下山的比喻，這就好比我們可以感覺到腳朝向的不同方向上，地形的傾斜程度不同。但是該跨出多長的步長呢？不確定。如果謹慎地小步走，情況可能比較穩定但是進展較慢（這就是步長較小的情況）。相反，如果想盡快下山，那就大步走吧，但結果也不一定盡如人意。在上面的代碼中就能看見反例，在某些點如果步長過大，反而可能越過最低點導致更高的損失值。

梯度告訴了我們損失函數下降最陡峭的方向。小步長下降穩定但進度慢，大步長進展快但是風險更大。采取大步長可能導致錯過最優點，讓損失值上升。步長（后面會稱其為學習率）將會是我們在調參中最重要的超參數之一。

敲黑板??！問題出現了??！

效率問題：你可能已經注意到，計算數值梯度的復雜性和參數的量線性相關。在本例中有30730個參數，所以損失函數每走一步就需要計算30731次損失函數的梯度?，F代神經網絡很容易就有上千萬的參數，因此這個問題只會越發嚴峻。顯然這個策略不適合大規模數據，我們需要更好的策略。

4.2 微分分析計算梯度（公式）

使用有限差值近似計算梯度比較簡單，但缺點在于終究只是近似（因為我們對于h值是選取了一個很小的數值，但真正的梯度定義中h趨向0的極限），且耗費計算資源太多。第二個梯度計算方法是利用微分來分析，能得到計算梯度的公式（不是近似），用公式計算梯度速度很快，唯一不好的就是實現的時候容易出錯。為了解決這個問題，在實際操作時常常將分析梯度法的結果和數值梯度法的結果作比較，以此來檢查其實現的正確性，這個步驟叫做梯度檢查。

用SVM的損失函數在某個數據點上的計算來舉例：

可以對函數進行微分。比如，對 W_yi 進行微分得到：

譯者注：原公式中1為空心字體，嘗試\mathbb{}等多種方法仍無法實現，請知友指點。
本菜注：知乎markdown都打不出來的符號，我一個簡書閹割版markdown能打出來么 <(￣ ﹌ ￣)>

其中 1 是一個示性函數，如果括號中的條件為真，那么函數值為1，如果為假，則函數值為0。雖然上述公式看起來復雜，但在代碼實現的時候比較簡單：只需要計算沒有滿足邊界值的分類的數量（因此對損失函數產生了貢獻），然后乘以 x_i 就是梯度了。注意，這個梯度只是對應正確分類的W的行向量的梯度，那些 j ≠ y_i 行的梯度是：

一旦將梯度的公式微分出來，代碼實現公式并用于梯度更新就比較順暢了。

5. 梯度下降

現在可以計算損失函數的梯度了，程序重復地計算梯度然后對參數進行更新，這一過程稱為梯度下降，他的普通版本是這樣的：

# 普通的梯度下降

while True:
  weights_grad = evaluate_gradient(loss_fun, data, weights)
  weights += - step_size * weights_grad # 進行梯度更新

這個簡單的循環在所有的神經網絡核心庫中都有。雖然也有其他實現最優化的方法（比如LBFGS），但是到目前為止，梯度下降是對神經網絡的損失函數最優化中最常用的方法。課程中，我們會在它的循環細節增加一些新的東西（比如更新的具體公式），但是核心思想不變，那就是我們一直跟著梯度走，直到結果不再變化。

小批量數據梯度下降（Mini-batch gradient descent）：在大規模的應用中（比如ILSVRC挑戰賽），訓練數據可以達到百萬級量級。如果像這樣計算整個訓練集，來獲得僅僅一個參數的更新就太浪費了。一個常用的方法是計算訓練集中的小批量（batches）數據。例如，在目前最高水平的卷積神經網絡中，一個典型的小批量包含256個例子，而整個訓練集是多少呢？一百二十萬個。這個小批量數據就用來實現一個參數更新：

# 普通的小批量數據梯度下降

while True:
  data_batch = sample_training_data(data, 256) # 256個數據
  weights_grad = evaluate_gradient(loss_fun, data_batch, weights)
  weights += - step_size * weights_grad # 參數更新

這個方法之所以效果不錯，是因為訓練集中的數據都是相關的。要理解這一點，可以想象一個極端情況：在ILSVRC中的120萬個圖像是1000張不同圖片的復制（每個類別1張圖片，每張圖片有1200張復制）。那么顯然計算這1200張復制圖像的梯度就應該是一樣的。對比120萬張圖片的數據損失的均值與只計算1000張的子集的數據損失均值時，結果應該是一樣的。實際情況中，數據集肯定不會包含重復圖像，那么小批量數據的梯度就是對整個數據集梯度的一個近似。因此，在實踐中通過計算小批量數據的梯度可以實現更快速地收斂，并以此來進行更頻繁的參數更新。

小批量數據策略有個極端情況，那就是每個批量中只有1個數據樣本，這種策略被稱為隨機梯度下降（Stochastic Gradient Descent 簡稱SGD），有時候也被稱為在線梯度下降。這種策略在實際情況中相對少見，因為向量化操作的代碼一次計算100個數據 比100次計算1個數據要高效很多。即使SGD在技術上是指每次使用1個數據來計算梯度，你還是會聽到人們使用SGD來指代小批量數據梯度下降（或者用MGD來指代小批量數據梯度下降，而BGD來指代則相對少見）。小批量數據的大小是一個超參數，但是一般并不需要通過交叉驗證來調參。它一般由存儲器的限制來決定的，或者干脆設置為同樣大小，比如32，64，128等。之所以使用2的指數，是因為在實際中許多向量化操作實現的時候，如果輸入數據量是2的倍數，那么運算更快。

關于Mini-batch gradient descent和SGD不太理解的同學，可以去看看Ng的Machine Learning課程的第10章。

6. 小結

非常好的一張圖！！

信息流的總結圖例。數據集中的(x,y)是給定的。權重從一個隨機數字開始，且可以改變。在前向傳播時，評分函數計算出類別的分類評分并存儲在向量f中。損失函數包含兩個部分：數據損失和正則化損失。其中，數據損失計算的是分類評分f和實際標簽y之間的差異，正則化損失只是一個關于權重的函數。在梯度下降過程中，我們計算權重的梯度（如果愿意的話，也可以計算數據上的梯度），然后使用它們來實現參數的更新。

在本節課中：

• 將損失函數比作了一個高維度的最優化地形，并嘗試到達它的最底部。最優化的工作過程可以看做一個蒙著眼睛的徒步者希望摸索著走到山的底部。在例子中，可見SVM的損失函數是分段線性的，并且是碗狀的。

• 提出了迭代優化的思想，從一個隨機的權重開始，然后一步步地讓損失值變小，直到最小。

• 函數的梯度給出了該函數最陡峭的上升方向。介紹了利用有限的差值來近似計算梯度的方法，該方法實現簡單但是效率較低（有限差值就是h，用來計算數值梯度）。

• 參數更新需要有技巧地設置步長。也叫學習率。如果步長太小，進度穩定但是緩慢，如果步長太大，進度快但是可能有風險。

• 討論權衡了數值梯度法和分析梯度法。數值梯度法計算簡單，但結果只是近似且耗費計算資源。分析梯度法計算準確迅速但是實現容易出錯，而且需要對梯度公式進行推導的數學基本功。因此，在實際中使用分析梯度法，然后使用梯度檢查來檢查其實現正確與否，其本質就是將分析梯度法的結果與數值梯度法的計算結果對比。

• 介紹了梯度下降算法，它在循環中迭代地計算梯度并更新參數。

預告：這節課的核心內容是：理解并能計算損失函數關于權重的梯度，是設計、訓練和理解神經網絡的核心能力。下節中，將介紹如何使用鏈式法則來高效地計算梯度，也就是通常所說的反向傳播（backpropagation）機制。該機制能夠對包含卷積神經網絡在內的幾乎所有類型的神經網絡的損失函數進行高效的最優化。

CS231n的反向傳播課程是我見過的最容易理解的了，雖然花了50min講，但是還是建議之前不太理解的同學能夠花時間看看！

再次感謝原作者的翻譯??！