主頁 > 後端開發 > 李航老師《統計學習方法》第二版第二章答案

李航老師《統計學習方法》第二版第二章答案

2021-04-08 11:37:05 後端開發

1、Minsky與Papert指出:感知機因為是線性模型,所以不能表達復雜的函式,比如異或(XOR),驗證感知機為什么不能表示異或,

解:

下面是異或的運算結果:
異或: 如果兩個值相同則異或操作的結果是0,如果不相同則為1
由此我們可以看到,這也是一個二分類的問題,異或的運算如表所示

XORab
a01
b10

方法一:直觀作圖法
如果我們去a = 0, b = 1,將上表的結果畫在二維平面,如下圖,
在這里插入圖片描述
我們可以看到,對于藍色的圓點和橙色的星星,無論我們怎么畫直線,都不可能將兩者分開,也就是我們不能畫出一條直線,是藍色的圓點位于直線的一側,是兩個橙色的星星位于直線的另外一側,因而,感知機無法表示異或(XOR),
方法二:直接計演算法
設正例的點是:(0, 0)和(1, 1),負例的點是:(1,0)和(0,1)
我們假設這兩類是可以線性分開的,設 f = w 1 x 1 + w 2 x 2 + b f = w_{1}x_{1}+w_{2}x_{2}+b f=w1?x1?+w2?x2?+b為對應的線性模型
由上面的假設可以知道點(0,0)和點(1,1)使 f > 0 f>0 f>0成立;
而點(1,0)和(0,1)使 f < 0 f<0 f<0成立,
將點(0,0)和點(1,1)代入我們得到:
{ 0 + 0 + b > 0 w 1 + w 2 + b > 0 (1) \left\{ \begin{aligned} 0+0+b>0 \\ w_{1}+w_{2}+b>0 \\ \end{aligned} \right.\tag{1} {0+0+b>0w1?+w2?+b>0?(1)
將點(1,0)和(0,1)代入我們得到
{ w 1 + 0 + b < 0 0 + w 2 + b < 0 (2) \left\{ \begin{aligned} w_{1}+0+b<0 \\ 0+w_{2}+b<0 \\ \end{aligned} \right.\tag{2} {w1?+0+b<00+w2?+b<0?(2)

將等式(1)中的兩式相加,我們得到
w 1 + w 2 + 2 b > 0 (3) w_{1}+w_{2}+2b>0\tag{3} w1?+w2?+2b>0(3)
將等式(2)中的兩式相加我們得到
w 1 + w 2 + 2 b < 0 (4) w_{1}+w_{2}+2b<0\tag{4} w1?+w2?+2b<0(4)
聯合公式(3)和(4),得出矛盾,因此不能線性可分,

方法3:使用這一章課后題的第三題來做
我們知道連接點(0,0)和(1,1)的直線段和連接點(1,0)和(0,1)的直線段在(0.5,0.5)處相交,也就是兩個點集的凸殼是相交的,因而不可能線性可分,

2、模仿例題2.1,構建從訓練資料集求解感知機模型的例子,

解:
這里不再構造,而是使用python將兩種演算法針對書上的例題2.1,將其實作出來,

import numpy as np
import matplotlib.pyplot as plt
        
epoch = 10

def Original_Percetron(x, y, lr = 0.1):
    '''
    Parameters
    ----------
    x : numpy array
        訓練資料集.
    y : numpy array
        訓練集標簽.
    lr : float
        梯度下降演算法學習率

    Returns
    -------
    回傳引數,并畫出分割超平面.
    '''
    #先獲取訓練集的維度
    s = np.shape(x)
    #初始化感知機模型引數
    w = np.zeros(shape = (1, s[1]))
    b = np.array(0.0)
    
    #下面開始更新模型
    for e in range(epoch):
        for i in range(s[0]):
            
            x_i = x[i: i + 1]
            
            if y[i] * (np.matmul(w, x_i.T)[0, 0] + b) <= 0:
                w = w + lr * y[i] * x_i
                b = b + lr * y[i]
                
    print('訓練%d個epoch后引數是:'%epoch, 'w 是:', w, 'b 是:', b)
    #下面開始畫圖
    num_positive = sum(np.array(y == 1, dtype = np.int64))
    plt.scatter(x[:num_positive, 0], x[:num_positive, 1], marker = "*", s = 100)
    plt.scatter(x[num_positive:, 0], x[num_positive:, 1], marker="o", s = 100)
    
    x11 = min(x[:,0]) - 1
    x12= max(x[:,0]) + 1
    
    x21 = - (w[0, 0] / w[0, 1]) * x11 - b / w[0, 0]
    x22 = - (w[0, 0] / w[0, 1]) * x12 - b / w[0, 0]
    
    plt.plot([x11, x12], [x21, x22])
    
    
def Dual_Perceptron(x, y, lr = 0.1):
    '''
    Parameters
    ----------
    x : numpy array
        訓練資料集.
    y : numpy array
        訓練集標簽.
    lr : float
        梯度下降演算法學習率

    Returns
    -------
    回傳引數,并畫出分割超平面.
    '''
    s = np.shape(x)
    #下面先計算Gram矩陣
    G = np.matmul(x, x.T)
    alpha = np.zeros(shape = (s[0]))
    
    for e in range(epoch):
        for i in range(s[0]):
            if y[i] * (sum(lr * alpha * y * G[:, i] + lr * alpha * y)) <= 0:
                alpha[i] += 1
                
    #下面開始畫圖
    num_positive = sum(np.array(y == 1, dtype = np.int64))
    plt.scatter(x[:num_positive, 0], x[:num_positive, 1], marker = "*", s = 100)
    plt.scatter(x[num_positive:, 0], x[num_positive:, 1], marker="o", s = 100)
    
    x11 = min(x[:,0]) - 1
    x12= max(x[:,0]) + 1
    
    w = np.zeros(shape = (s[1]))
    b = 0.0
    for i in range(s[0]):
        w = w + alpha[i] * y[i] * x[i] * lr
        b = b + alpha[i] * lr * y[i]
        
    x21 = - (w[0] / w[1]) * x11 - b / w[0]
    x22 = - (w[0] / w[1]) * x12 - b / w[0]
    
    plt.plot([x11, x12], [x21, x22])
     

if __name__ == '__main__':
    x = np.array([[3,3],
                  [4,3],
                  [1,1]])
    y = np.array([1,1,-1])
    
    Original_Percetron(x, y, lr = 0.1)
    Dual_Perceptron(x, y, lr = 0.1)

原始的感知機的輸出:
在這里插入圖片描述
對偶感知機的輸出:
在這里插入圖片描述

3、證明一下定理:樣本集線性可分的充分必要條件是正實體點集所構成的凸殼與負實體點集構成的凸殼互不相交,

證明:
充分性:由凸殼不相交推導線性可分
因為正實體和負實體分別構成的凸殼是不相交的,因而,一定有一個超平面將兩個凸殼線性分開,再根據凸殼的定義可以看到,正實體點一定在正實體點集構成的凸殼里面,負實體點一定在負實體點集構成的凸殼內,因而既然可以將兩個凸殼分開,那么也可以將兩類點集分開,

其實這個充分性的證明還是比較顯然的,下面再提供一種數學方面的思路,
d = i n f x , y d ( x , y ) (1) d = \underset{x,y}{inf}\quad d(x,y) \tag{1} d=x,yinf?d(x,y)(1)
其中 x , y x,y x,y 分別是兩個凸殼中的點, d ( x , y ) d(x,y) d(x,y) 表示 x , y x,y x,y 兩點之間的距離,
因為凸殼是閉集,所以一定可以在兩個凸殼內找到兩點 x , y x,y x,y,使得 x , y x,y x,y 之間的距離就是這兩個凸殼之間的距離,然后,我們做 x , y x,y x,y 兩點的垂直平分線,該垂直平分線就是我們要找的分割平面,

必要性:已知線性可分,證明兩個凸殼不相交,
設集合 A , B A,B A,B分別是正實體和負實體對應的點集,標記 y A = 1 , y B = ? 1 y_{A} = 1,y_{B} = -1 yA?=1,yB?=?1,因為 A , B A,B A,B線性可分,所有設 f = w x + b (2) f = wx+b\tag{2} f=wx+b(2)為分割超平面,
所以對于任何一點 x i A ∈ A , x j B ∈ B x_{i}^{A}\in A,x_{j}^{B}\in B xiA?A,xjB?B,我們都有
w ? x i A + b > 0 (3) w*x_{i}^{A}+b>0\tag{3} w?xiA?+b>0(3)
w ? x j B + b < 0 (4) w*x_{j}^{B}+b<0\tag{4} w?xjB?+b<0(4)
如果兩個凸殼相交,那么一定有集合 A A A中的點在 c o n v ( B ) conv(B) conv(B)內或者有集合 B B B內的點在凸殼 c o n v ( A ) conv(A) conv(A)內,或者有 c o n v ( A ) conv(A) conv(A)內點 c o n v ( B ) conv(B) conv(B)內,
不妨假設前者成立,也即是有集合 A A A中的點 x i A x_{i}^{A} xiA? c o n v ( B ) conv(B) conv(B)內,如果只是 c o n v ( A ) conv(A) conv(A)內點 a a a c o n v ( B ) conv(B) conv(B)內,因為 a a a 是集合 A A A中點的線性組合,證明也是類似的,所以我們下面只證明集合 A A A中的點 x i A x_{i}^{A} xiA? c o n v ( B ) conv(B) conv(B)內,
也就是 x i A ∈ c o n v ( B ) x_{i}^{A}\in conv(B) xiA?conv(B),所以 A A A內的點 x i A x_{i}^{A} xiA?可以被 B B B內的點進行凸組合表示,不妨設為:
x i A = ∑ j = 1 N λ j x j B (5) x_{i}^{A} = \sum_{j = 1}^{N}\lambda _{j}x_{j}^{B}\tag{5} xiA?=j=1N?λj?xjB?(5)
但是根據公式(3)我們得
w ? x i A + b > 0 w*x_{i}^{A}+b>0 w?xiA?+b>0
但是
w ? x i A + b = w ? ∑ j = 1 N λ j x j B + b = ∑ j = 1 N λ j w x j B + b = ∑ j = 1 N λ j w x j B + ∑ j = 1 N λ j b = ∑ j = 1 N λ j ( w ? x j B + b ) < 0 (6) w*x_{i}^{A}+b = w*\sum_{j = 1}^{N}\lambda _{j}x_{j}^{B} + b \\=\sum_{j = 1}^{N}\lambda _{j}wx_{j}^{B}+b \\=\sum_{j = 1}^{N}\lambda _{j}wx_{j}^{B}+\sum_{j = 1}^{N}\lambda _{j}b \\=\sum_{j = 1}^{N}\lambda _{j}(w*x_{j}^{B}+b)<0\tag{6} w?xiA?+b=w?j=1N?λj?xjB?+b=j=1N?λj?wxjB?+b=j=1N?λj?wxjB?+j=1N?λj?b=j=1N?λj?(w?xjB?+b)<0(6)
因而矛盾,所以 c o n v ( A ) , c o n v ( B ) conv(A),conv(B) conv(A),conv(B)一定不相交的,

綜上,得證!!!

轉載請註明出處,本文鏈接:https://www.uj5u.com/houduan/273727.html

標籤:python

上一篇:tkinter點擊按鈕實作圖片的切換

下一篇:【Python基礎】函式:輸出一個區間內的素數

標籤雲
其他(157675) Python(38076) JavaScript(25376) Java(17977) C(15215) 區塊鏈(8255) C#(7972) AI(7469) 爪哇(7425) MySQL(7132) html(6777) 基礎類(6313) sql(6102) 熊猫(6058) PHP(5869) 数组(5741) R(5409) Linux(5327) 反应(5209) 腳本語言(PerlPython)(5129) 非技術區(4971) Android(4554) 数据框(4311) css(4259) 节点.js(4032) C語言(3288) json(3245) 列表(3129) 扑(3119) C++語言(3117) 安卓(2998) 打字稿(2995) VBA(2789) Java相關(2746) 疑難問題(2699) 细绳(2522) 單片機工控(2479) iOS(2429) ASP.NET(2402) MongoDB(2323) 麻木的(2285) 正则表达式(2254) 字典(2211) 循环(2198) 迅速(2185) 擅长(2169) 镖(2155) 功能(1967) .NET技术(1958) Web開發(1951) python-3.x(1918) HtmlCss(1915) 弹簧靴(1913) C++(1909) xml(1889) PostgreSQL(1872) .NETCore(1853) 谷歌表格(1846) Unity3D(1843) for循环(1842)

熱門瀏覽
  • 【C++】Microsoft C++、C 和匯編程式檔案

    ......

    uj5u.com 2020-09-10 00:57:23 more
  • 例外宣告

    相比于斷言適用于排除邏輯上不可能存在的狀態,例外通常是用于邏輯上可能發生的錯誤。 例外宣告 Item 1:當函式不可能拋出例外或不能接受拋出例外時,使用noexcept 理由 如果不打算拋出例外的話,程式就會認為無法處理這種錯誤,并且應當盡早終止,如此可以有效地阻止例外的傳播與擴散。 示例 //不可 ......

    uj5u.com 2020-09-10 00:57:27 more
  • Codeforces 1400E Clear the Multiset(貪心 + 分治)

    鏈接:https://codeforces.com/problemset/problem/1400/E 來源:Codeforces 思路:給你一個陣列,現在你可以進行兩種操作,操作1:將一段沒有 0 的區間進行減一的操作,操作2:將 i 位置上的元素歸零。最終問:將這個陣列的全部元素歸零后操作的最少 ......

    uj5u.com 2020-09-10 00:57:30 more
  • UVA11610 【Reverse Prime】

    本人看到此題沒有翻譯,就附帶了一個自己的翻譯版本 思考 這一題,它的第一個要求是找出所有 $7$ 位反向質數及其質因數的個數。 我們應該需要質數篩篩選1~$10^{7}$的所有數,這里就不慢慢介紹了。但是,重讀題,我們突然發現反向質數都是 $7$ 位,而將它反過來后的數字卻是 $6$ 位數,這就說明 ......

    uj5u.com 2020-09-10 00:57:36 more
  • 統計區間素數數量

    1 #pragma GCC optimize(2) 2 #include <bits/stdc++.h> 3 using namespace std; 4 bool isprime[1000000010]; 5 vector<int> prime; 6 inline int getlist(int ......

    uj5u.com 2020-09-10 00:57:47 more
  • C/C++編程筆記:C++中的 const 變數詳解,教你正確認識const用法

    1、C中的const 1、區域const變數存放在堆疊區中,會分配記憶體(也就是說可以通過地址間接修改變數的值)。測驗代碼如下: 運行結果: 2、全域const變數存放在只讀資料段(不能通過地址修改,會發生寫入錯誤), 默認為外部聯編,可以給其他源檔案使用(需要用extern關鍵字修飾) 運行結果: ......

    uj5u.com 2020-09-10 00:58:04 more
  • 【C++犯錯記錄】VS2019 MFC添加資源不懂如何修改資源宏ID

    1. 首先在資源視圖中,添加資源 2. 點擊新添加的資源,復制自動生成的ID 3. 在解決方案資源管理器中找到Resource.h檔案,編輯,使用整個專案搜索和替換的方式快速替換 宏宣告 4. Ctrl+Shift+F 全域搜索,點擊查找全部,然后逐個替換 5. 為什么使用搜索替換而不使用屬性視窗直 ......

    uj5u.com 2020-09-10 00:59:11 more
  • 【C++犯錯記錄】VS2019 MFC不懂的批量添加資源

    1. 打開資源頭檔案Resource.h,在其中預先定義好宏 ID(不清楚其實ID值應該設定多少,可以先新建一個相同的資源項,再在這個資源的ID值的基礎上遞增即可) 2. 在資源視圖中選中專案資源,按F7編輯資源檔案,按 ID 型別 相對路徑的形式添加 資源。(別忘了先把檔案拷貝到專案中的res檔案 ......

    uj5u.com 2020-09-10 01:00:19 more
  • C/C++編程筆記:關于C++的參考型別,專供新手入門使用

    今天要講的是C++中我最喜歡的一個用法——參考,也叫別名。 參考就是給一個變數名取一個變數名,方便我們間接地使用這個變數。我們可以給一個變數創建N個參考,這N + 1個變數共享了同一塊記憶體區域。(參考型別的變數會占用記憶體空間,占用的記憶體空間的大小和指標型別的大小是相同的。雖然參考是一個物件的別名,但 ......

    uj5u.com 2020-09-10 01:00:22 more
  • 【C/C++編程筆記】從頭開始學習C ++:初學者完整指南

    眾所周知,C ++的學習曲線陡峭,但是花時間學習這種語言將為您的職業帶來奇跡,并使您與其他開發人員區分開。您會更輕松地學習新語言,形成真正的解決問題的技能,并在編程的基礎上打下堅實的基礎。 C ++將幫助您養成良好的編程習慣(即清晰一致的編碼風格,在撰寫代碼時注釋代碼,并限制類內部的可見性),并且由 ......

    uj5u.com 2020-09-10 01:00:41 more
最新发布
  • Rust中的智能指標:Box<T> Rc<T> Arc<T> Cell<T> RefCell<T> Weak

    Rust中的智能指標是什么 智能指標(smart pointers)是一類資料結構,是擁有資料所有權和額外功能的指標。是指標的進一步發展 指標(pointer)是一個包含記憶體地址的變數的通用概念。這個地址參考,或 ” 指向”(points at)一些其 他資料 。參考以 & 符號為標志并借用了他們所 ......

    uj5u.com 2023-04-20 07:24:10 more
  • Java的值傳遞和參考傳遞

    值傳遞不會改變本身,參考傳遞(如果傳遞的值需要實體化到堆里)如果發生修改了會改變本身。 1.基本資料型別都是值傳遞 package com.example.basic; public class Test { public static void main(String[] args) { int ......

    uj5u.com 2023-04-20 07:24:04 more
  • [2]SpinalHDL教程——Scala簡單入門

    第一個 Scala 程式 shell里面輸入 $ scala scala> 1 + 1 res0: Int = 2 scala> println("Hello World!") Hello World! 檔案形式 object HelloWorld { /* 這是我的第一個 Scala 程式 * 以 ......

    uj5u.com 2023-04-20 07:23:58 more
  • 理解函式指標和回呼函式

    理解 函式指標 指向函式的指標。比如: 理解函式指標的偽代碼 void (*p)(int type, char *data); // 定義一個函式指標p void func(int type, char *data); // 宣告一個函式func p = func; // 將指標p指向函式func ......

    uj5u.com 2023-04-20 07:23:52 more
  • Django筆記二十五之資料庫函式之日期函式

    本文首發于公眾號:Hunter后端 原文鏈接:Django筆記二十五之資料庫函式之日期函式 日期函式主要介紹兩個大類,Extract() 和 Trunc() Extract() 函式作用是提取日期,比如我們可以提取一個日期欄位的年份,月份,日等資料 Trunc() 的作用則是截取,比如 2022-0 ......

    uj5u.com 2023-04-20 07:23:45 more
  • 一天吃透JVM面試八股文

    什么是JVM? JVM,全稱Java Virtual Machine(Java虛擬機),是通過在實際的計算機上仿真模擬各種計算機功能來實作的。由一套位元組碼指令集、一組暫存器、一個堆疊、一個垃圾回收堆和一個存盤方法域等組成。JVM屏蔽了與作業系統平臺相關的資訊,使得Java程式只需要生成在Java虛擬機 ......

    uj5u.com 2023-04-20 07:23:31 more
  • 使用Java接入小程式訂閱訊息!

    更新完微信服務號的模板訊息之后,我又趕緊把微信小程式的訂閱訊息給實作了!之前我一直以為微信小程式也是要企業才能申請,沒想到小程式個人就能申請。 訊息推送平臺🔥推送下發【郵件】【短信】【微信服務號】【微信小程式】【企業微信】【釘釘】等訊息型別。 https://gitee.com/zhongfuch ......

    uj5u.com 2023-04-20 07:22:59 more
  • java -- 緩沖流、轉換流、序列化流

    緩沖流 緩沖流, 也叫高效流, 按照資料型別分類: 位元組緩沖流:BufferedInputStream,BufferedOutputStream 字符緩沖流:BufferedReader,BufferedWriter 緩沖流的基本原理,是在創建流物件時,會創建一個內置的默認大小的緩沖區陣列,通過緩沖 ......

    uj5u.com 2023-04-20 07:22:49 more
  • Java-SpringBoot-Range請求頭設定實作視頻分段傳輸

    老實說,人太懶了,現在基本都不喜歡寫筆記了,但是網上有關Range請求頭的文章都太水了 下面是抄的一段StackOverflow的代碼...自己大修改過的,寫的注釋挺全的,應該直接看得懂,就不解釋了 寫的不好...只是希望能給視頻網站開發的新手一點點幫助吧. 業務場景:視頻分段傳輸、視頻多段傳輸(理 ......

    uj5u.com 2023-04-20 07:22:42 more
  • Windows 10開發教程_編程入門自學教程_菜鳥教程-免費教程分享

    教程簡介 Windows 10開發入門教程 - 從簡單的步驟了解Windows 10開發,從基本到高級概念,包括簡介,UWP,第一個應用程式,商店,XAML控制元件,資料系結,XAML性能,自適應設計,自適應UI,自適應代碼,檔案管理,SQLite資料庫,應用程式到應用程式通信,應用程式本地化,應用程式 ......

    uj5u.com 2023-04-20 07:22:35 more