主頁 > 後端開發 > k-近鄰(knn)演算法

k-近鄰(knn)演算法

2020-11-24 19:29:37 後端開發

k-近鄰(knn)演算法的簡介和實作

  • 一、k近鄰演算法的概述
  • 二、用python實作k近鄰演算法
    • 1.演算法實作
    • 2.封裝函式
  • 三、k近鄰演算法案例應用
  • 四、總結

最近小阿奇在學習機器學習演算法,所以決定把相關原理和代碼整理出來和小伙伴們一起分享

一、k近鄰演算法的概述

k-近鄰演算法(k-Nearest Neighbour algorithm),又稱為KNN演算法,KNN的作業原理:給定一個已知標簽類別的訓練資料集,輸入沒有標簽的新資料后,在訓練資料集中找到與新資料最鄰近的k個實體,如果這k個實體的多數屬于某個類別,那么新資料就屬于這個類別,可以簡單理解為:由那些離X最近的k個點來投票決定X歸為哪一類,
在這里插入圖片描述
如上圖所示
當k的范圍改變,綠色圓點會有不一樣的類別
當k=3時(實線范圍),綠色圓點屬于紅色三角這種類別;
當k=5時(虛線范圍),綠色圓點屬于藍色方塊這種類別,

那如何判斷一部電影屬于愛情片還是動作片呢?

電影名稱打斗鏡頭接吻鏡頭電影型別
無問西東1101愛情片
后來的我們589愛情片
前任31297愛情片
紅海行動1085動作片
唐人街探案1129動作片
戰狼21158動作片
新電影2467

我們將表中已有資料作為樣本資料,打斗鏡頭和接吻鏡頭作為資料特征,如何用knn演算法給出新電影的電影型別呢?
在這里插入圖片描述將資料放在坐標圖中可以發現,新電影更加偏向于愛情片,那么從k近鄰角度來說,它主要通過未知點和已知資料集中的資料點之間的距離進行判斷得出最終結果,
在二維平面角度距離公式:
在這里插入圖片描述那么在多維空間里,歐式距離公式可以幫我們很好地解決:
在這里插入圖片描述k近鄰演算法的主要步驟:
(1) 計算已知類別資料集中的點與當前點之間的距離;
(2) 按照距離遞增次序排序
(3) 選取與當前點距離最小的k個點;
(4) 確定前k個點所在類別的出現頻率;
(5) 回傳前k個點出現頻率最高的類別作為當前點的預測類別,
當k=4時,那么在這個電影例子中,把距離按照升序排列,距離綠點電影最近的前4個的電影分別是《后來的我們》、《前任3》、《無問西東》和《紅海行動》,這四部電影的類別統計為愛情片:動作片=3:1,出現頻率最高的類別為愛情片,所以在k=4時,綠點電影的類別為愛情片,這個判別程序就是k-近鄰演算法,
然而k的選取得不同也得到不一樣的結果,那么如何確定k的值呢?
李航書上講到,我們一般選取一個較小的數值,通常采取交叉驗證法來選取最優的k值,(也就是說,選取k值很重要的關鍵是實驗調參,類似于神經網路選取多少層這種,通過調整超引數來得到一個較好的結果)

二、用python實作k近鄰演算法

1.演算法實作

import pandas as pd

#匯入資料,將其轉換成dataframe格式
rowdata={'電影名稱':['無問西東','后來的我們','前任3','紅海行動','唐人街探案','戰狼2'],
'打斗鏡頭':[1,5,12,108,112,115],
'接吻鏡頭':[101,89,97,5,9,8],
'電影型別':['愛情片','愛情片','愛情片','動作片','動作片','動作片']}
movie_data= pd.DataFrame(rowdata)

#根據上面距離公式,求出距離
new_data = [24,67]
dist = list((((movie_data.iloc[:6,1:3]-new_data)**2).sum(1))**0.5)

#首先對資料進行切分,然后升序排列求出前四
dist_l = pd.DataFrame({'dist': dist, 'labels': (movie_data.iloc[:6, 3])})
dr = dist_l.sort_values(by = 'dist')[: 4]

#統計前k個出現的類別頻率
re = dr.loc[:,'labels'].value_counts()

#選擇頻率最高的輸出
result = []
result.append(re.index[0])
print(result)

2.封裝函式

def classify0(inX,dataSet,k):
    result = []
    dist = list((((dataSet.iloc[:,1:3]-inX)**2).sum(1))**0.5)
    dist_l = pd.DataFrame({'dist':dist,'labels':(dataSet.iloc[:, 3])})
    dr = dist_l.sort_values(by = 'dist')[: k]
    re = dr.loc[:, 'labels'].value_counts()
    result.append(re.index[0])
return result
inX = new_data
dataSet = movie_data
k = 3
print(classify0(inX,dataSet,k))

三、k近鄰演算法案例應用

import matplotlib as mpl
import matplotlib.pyplot as plt

datingTest=pd.read_table('datingTestSet.txt',header=None) #匯入資料
#把不同標簽用顏色區分
Colors = []
for i in range(datingTest.shape[0]):
   m = datingTest.iloc[i,-1]
if m=='didntLike':
   Colors.append('black')
if m=='smallDoses':
   Colors.append('orange')
if m=='largeDoses':
   Colors.append('red')
#繪制兩兩特征之間的散點圖
plt.rcParams['font.sans-serif']=['Simhei'] #圖中字體設定為黑體
pl=plt.figure(figsize=(12,8))
fig1=pl.add_subplot(221)
plt.scatter(datingTest.iloc[:,1],datingTest.iloc[:,2],marker='.',c=Colors)
plt.xlabel('玩游戲視頻所占時間比')
plt.ylabel('每周消費冰淇淋公升數')

fig2=pl.add_subplot(222)
plt.scatter(datingTest.iloc[:,0],datingTest.iloc[:,1],marker='.',c=Colors)
plt.xlabel('每年飛行公里數')
plt.ylabel('消費冰淇淋公斤數')

fig3=pl.add_subplot(223)
plt.scatter(datingTest.iloc[:,0],datingTest.iloc[:,2],marker='.',c=Colors)
plt.xlabel('每年飛行公里數')
plt.ylabel('玩游戲所占時間比')
plt.show()
#資料歸一化'''
def minmax(dataSet):
    minDF=dataSet.min()
    maxDF=dataSet.max()
    normSet=(dataSet-minDF)/(maxDF-minDF)
    return normSet
datingT=pd.concat([minmax(datingTest.iloc[:,:3]),datingTest.iloc[:,3]],axis=1)  #對資料進行切分結合,axis=1保留縱軸
#資料進行切分
"""
將資料切分為訓練集和測驗集
dataSet:原始資料集
rate:訓練集所占比例
回傳切分好的訓練集和測驗集
"""
def randSplit(dataSet,rate=0.9):
    n=dataSet.shape[0]
    m=int(n*rate)
    train=dataSet.iloc[:m,:]
    test=dataSet.iloc[m:,:]
    test.index=range(test.shape[0])  #test資料進行排序
    return train,test
train,test=randSplit(datingT)
#分類器
def datingclassify(train,test,k):
    n=train.shape[1]-1   #列數
    m=train.shape[0]     #行數
    result=[]
    for i in range(m):
         dist=list((((train.iloc[:,:n]-test.iloc[i,:n])**2).sum(1))**5)   #歐式距離公式
         dist_l=pd.DataFrame({'dist':dist,'labels':(train.iloc[:,n])}) #進行標簽處理
         dr=dist_l.sort_values(by ='dist')[:k]   #升序排列
         re=dr.loc[:,'labels'].value_counts()
         result.append(re.index[0]) #進行排序
    result=pd.Series(result)
    test['predict']=result        #增加新的預測值列
    acc=(test.iloc[:,-1]==test.iloc[:,-2]).mean()    #確認準確率mean取巧的程序
    print(f'模型預測準確率為{acc}')
    return test

print(datingclassify(train,test,5))

在距離公式中差值最大的屬性對計算結果的影響最大,原因僅僅是因為它的數值比較大,所以我們要進行數值歸一化的處理,使得這三個特征的權重相等,
資料歸一化的處理方法有很多種,比如0-1標準化、Z-score標準化、Sigmoid壓縮法等等,在這里我們使用最簡單
的0-1標準化,公式如下:
在這里插入圖片描述

四、總結

1.k近鄰演算法核心思想是,即是給定一個訓練資料集,對新的輸入實體,在訓練資料集中找到與該實體最鄰近的K個實體,這K個實體的多數屬于某個類,就把該輸入實體分類到這個類中,

2.與該實體最近鄰的k個實體,這個最近鄰的定義是通過不同距離函式來定義,我們最常用的是歐式距離,

3.為了保證每個特征同等重要性,我們這里對每個特征進行歸一化,

4.k值的選取,既不能太大,也不能太小,何值為最好,需要實驗調整引數確定!

希望對大家有所幫助!

參考:
機器學習(周志華),機器學習實戰,李航《統計學習方法》,橘安醬的機器學習

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

標籤:python

上一篇:用最少數量的箭引爆氣球--小白嘔心制作

下一篇:Windows下用Anaconda3安裝PyTorch環境(詳細教程)

標籤雲
其他(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