日韩av黄I国产麻豆传媒I国产91av视频在线观看I日韩一区二区三区在线看I美女国产在线I麻豆视频国产在线观看I成人黄色短片

歡迎訪問(wèn) 生活随笔!

生活随笔

當(dāng)前位置: 首頁(yè) > 编程资源 > 编程问答 >内容正文

编程问答

SVM原理以及Tensorflow 实现SVM分类(附代码讲了一下原理)

發(fā)布時(shí)間:2024/9/20 编程问答 42 豆豆
生活随笔 收集整理的這篇文章主要介紹了 SVM原理以及Tensorflow 实现SVM分类(附代码讲了一下原理) 小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
  • 1.1. SVM介紹
  • 1.2. 工作原理
    • 1.2.1. 幾何間隔和函數(shù)間隔
    • 1.2.2. 最大化間隔
  • 1.3. 軟間隔
  • 1.4. SMO算法
  • 1.5. 核函數(shù)
  • 1.6. 實(shí)例

1.1. SVM介紹

SVM(Support Vector Machines)——支持向量機(jī)是在所有知名的數(shù)據(jù)挖掘算法中最健壯,最準(zhǔn)確的方法之一,它屬于二分類(lèi)算法,可以支持線性和非線性的分類(lèi)。發(fā)展到今天,SVM已經(jīng)可以支持多分類(lèi)了,但在這一章里,我們著重講支持向量機(jī)在二分類(lèi)問(wèn)題中的工作原理。
假設(shè)在一個(gè)二維線性可分的數(shù)據(jù)集中,圖一A所示,我們要找到一個(gè)超平面把兩組數(shù)據(jù)分開(kāi),這時(shí),我們認(rèn)為線性回歸的直線或邏輯回歸的直線也能夠做這個(gè)分類(lèi),這條直線可以是圖一B中的直線,也可以是圖一C中的直線,或者圖一D中的直線,但哪條直線才最好呢,也就是說(shuō)哪條直線能夠達(dá)到最好的泛化能力呢?那就是一個(gè)能使兩類(lèi)之間的空間大小最大的一個(gè)超平面。
這個(gè)超平面在二維平面上看到的就是一條直線,在三維空間中就是一個(gè)平面...,因此,我們把這個(gè)劃分?jǐn)?shù)據(jù)的決策邊界統(tǒng)稱(chēng)為超平面。離這個(gè)超平面最近的點(diǎn)就叫做支持向量,點(diǎn)到超平面的距離叫間隔。支持向量機(jī)就是要使超平面和支持向量之間的間隔盡可能的大,這樣超平面才可以將兩類(lèi)樣本準(zhǔn)確的分開(kāi),而保證間隔盡可能的大就是保證我們的分類(lèi)器誤差盡可能的小,盡可能的健壯。

圖一

1.2. 工作原理

1.2.1. 幾何間隔和函數(shù)間隔

在最大化支持向量到超平面距離前,我們首先要定義我們的超平面h(x)h(x)(稱(chēng)為超平面的判別函數(shù),也稱(chēng)給定ww和bb的泛函間隔),其中ww為權(quán)重向量,bb為偏移向量:

?

h(x)=wTx+bh(x)=wTx+b

?

樣本xx到最優(yōu)超平面的幾何間隔為:

?

r=h(x)||w||=wTx+b||w||r=h(x)||w||=wTx+b||w||

?

||w||||w||是向量ww的內(nèi)積,是個(gè)常數(shù),即||w||=w02+w12+...+wn2????????????????√||w||=w02+w12+...+wn2,而h(x)h(x)就是下面要介紹的函數(shù)間隔。

函數(shù)間隔

?

r?=h(x)r^=h(x)

?

函數(shù)間隔h(x)h(x)它是一個(gè)并不標(biāo)準(zhǔn)的間隔度量,是人為定義的,它不適合用來(lái)做最大化的間隔值,因?yàn)?#xff0c;一旦超平面固定以后,如果我們?nèi)藶榈姆糯蠡蚩s小ww和bb值,那這個(gè)超平面也會(huì)無(wú)限的放大或縮小,這將對(duì)分類(lèi)造成嚴(yán)重影響。而幾何間隔是函數(shù)間隔除以||w||||w||,當(dāng)ww的值無(wú)限放大或縮小時(shí),||w||||w||也會(huì)放大或縮小,而整個(gè)rr保持不變,它只隨著超平面的變動(dòng)而變動(dòng),不受兩個(gè)參數(shù)的影響。因而,我們用幾何間隔來(lái)做最大化間隔度量。

1.2.2. 最大化間隔

在支持向量機(jī)中,我們把幾何間隔rr作為最大化間隔進(jìn)行分析,并且采用-1和1作為類(lèi)別標(biāo)簽,什么采用-1和+1,而不是0和1呢?這是由于-1和+1僅僅相差一個(gè)符號(hào),方便數(shù)學(xué)上的處理。我們可以通過(guò)一個(gè)統(tǒng)一公式來(lái)表示間隔或者數(shù)據(jù)點(diǎn)到分隔超平面的距離,同時(shí)不必?fù)?dān)心數(shù)據(jù)到底是屬于-1還是+1類(lèi)。
我們一步一步的進(jìn)行分析,首先如下圖,在這個(gè)R2R2空間中,假設(shè)我們已經(jīng)確定了一個(gè)超平面,這個(gè)超平面的函數(shù)關(guān)系式應(yīng)該是h(x)=wTx+b=0h(x)=wTx+b=0,這個(gè)式子表示我們圖中的那條虛線,很明顯,這個(gè)式子意思是說(shuō)點(diǎn)x在超平面上,但我們要想使所有的點(diǎn)都盡可能的遠(yuǎn)離這個(gè)超平面,我們只要保證離這個(gè)超平面最近的點(diǎn)遠(yuǎn)離這個(gè)超平面,也就是說(shuō)這些叫支持向量的點(diǎn)x?x?需要盡可能的遠(yuǎn)離它。

?

我們把其中一個(gè)支持向量x?x?到最優(yōu)超平面的距離定義為:

?

r?=h(x?)||w||=?????????1||w||?1||w||if:y?=h(x?)=+1if:y?=h(x?)=?1r?=h(x?)||w||={1||w||if:y?=h(x?)=+1?1||w||if:y?=h(x?)=?1

?

這是我們通過(guò)把函數(shù)間隔h(x)h(x)固定為1而得來(lái)的。我們可以把這個(gè)式子想象成還存在兩個(gè)平面,這兩個(gè)平面分別是wTxs+b=1wTxs+b=1和wTxs+b=?1wTxs+b=?1,對(duì)應(yīng)上圖中的兩根實(shí)線。這些支持向量xsxs就在這兩個(gè)平面上,這兩個(gè)平面離最優(yōu)超平面的距離越大,我們的間隔也就越大。對(duì)于其他的點(diǎn)xixi如果滿(mǎn)足wTxi+b>1wTxi+b>1,則被分為1類(lèi),如果滿(mǎn)足滿(mǎn)足wTxi+b<?1wTxi+b<?1,則被分為-1類(lèi)。即有約束條件:

?

?????wTxi+b?1wTxi+b??1yi=+1yi=?1{wTxi+b?1yi=+1wTxi+b??1yi=?1

?

支持向量到超平面的距離知道后,那么分離的間隔ρρ很明顯就為:

?

ρ=2r?=2||w||ρ=2r?=2||w||

?

這下我們就要通過(guò)找到最優(yōu)的ww和bb來(lái)最大化ρρ了,感覺(jué)又像回到了邏輯回歸或線性回歸的例子。但是這里,我們最大化ρρ值需要有條件限制,即:

?

???????maxw,b2||w||yi(wTxi+b)?1,?(i=1,..,n){maxw,b2||w||yi(wTxi+b)?1,?(i=1,..,n)

?

yi(wTxi+b)yi(wTxi+b)的意思是通過(guò)判斷yiyi和wTxi+bwTxi+b是否同號(hào)來(lái)確定分類(lèi)結(jié)果。
接著,為了計(jì)算方便,我們把上式最大化ρρ換成:

?

???????minw,b12||w||2yi(wTxi+b)?1,?(i=1,..,n){minw,b12||w||2yi(wTxi+b)?1,?(i=1,..,n)

?

這種式子通常我們用拉格朗日乘數(shù)法來(lái)求解,即:

?

L(x)=f(x)+∑αg(x)L(x)=f(x)+∑αg(x)

?

f(x)f(x)是我們需要最小化的目標(biāo)函數(shù),g(x)g(x)是不等式約束條件,即前面的yi(wTxi+b)?1yi(wTxi+b)?1,αα是對(duì)應(yīng)的約束系數(shù),也叫拉格朗日乘子。為了使得拉格朗日函數(shù)得到最優(yōu)化解,我們需要加入能使該函數(shù)有最優(yōu)化解法的KKT條件,或者叫最優(yōu)化條件、充要條件。即假設(shè)存在一點(diǎn)x?x?

1.2.2.0.0.1.?L(x?)L(x?)對(duì)x?x?求導(dǎo)為0

1.2.2.0.0.2.?αigi(x?)=0αigi(x?)=0,對(duì)于所有的i=1,.....,ni=1,.....,n

這樣構(gòu)建我們的拉格朗日函數(shù)為:

?

L(w,b,α)=12wTw?∑i=1nαi[yi(wTxi+b)?1]L(w,b,α)=12wTw?∑i=1nαi[yi(wTxi+b)?1]

?

以上的KKT條件αi[yi(wTxi+b)?1]=0αi[yi(wTxi+b)?1]=0表示,只有距離最優(yōu)超平面的支持向量(xi,yi)(xi,yi)對(duì)應(yīng)的αα非零,其他所有點(diǎn)集的αα等于零。綜上所述,引入拉格朗日乘子以后,我們的目標(biāo)變?yōu)?#xff1a;

?

minw,bmaxα?0L(w,b,α)minw,bmaxα?0L(w,b,α)

?

即先求得αα的極大值,再求ww和bb的極小值。可以把該問(wèn)題轉(zhuǎn)為為等價(jià)的凸優(yōu)化和對(duì)偶問(wèn)題來(lái)求解,對(duì)于凸優(yōu)化和對(duì)偶問(wèn)題可以參考《凸優(yōu)化》這本書(shū),因?yàn)樵摾碚摽梢哉靡槐緯?shū)來(lái)介紹了,筆者在這里也只能點(diǎn)到為止了。通過(guò)對(duì)偶,我們的目標(biāo)可以又變成:

?

maxα?0minw,bL(w,b,α)maxα?0minw,bL(w,b,α)

?

即先求得ww和bb的極小值,在求αα的極大值。用L(w,b,α)L(w,b,α)對(duì)ww和bb分別求偏導(dǎo),并令其等于0:

?

????????L(w,b,α)?w=0??L(w,b,α)?b=0{?L(w,b,α)?w=0??L(w,b,α)?b=0

?

得:

?

???w=∑ni=1αiyixi?∑ni=1αiyi=0{w=∑i=1nαiyixi?∑i=1nαiyi=0

?

把該式代入原來(lái)的的拉格朗日式子可得(推導(dǎo)過(guò)程省略):

?

W(α)=∑i=1nαi?12∑i=1n∑j=1nαiαjyiyjxiTxjW(α)=∑i=1nαi?12∑i=1n∑j=1nαiαjyiyjxiTxj

?

?

∑i=1nαiyi=0?,αi?0(i=1,...,n)∑i=1nαiyi=0?,αi?0(i=1,...,n)

?

該W(α)W(α)函數(shù)消去了向量ww和向量bb,僅剩αα這個(gè)未知參數(shù),只要我們能夠最大化W(α)W(α),就能求出對(duì)應(yīng)的αα,進(jìn)而求得ww和bb。對(duì)于如何求解αα,SMO算法給出了完美的解決方案,下一節(jié)我們?cè)敿?xì)講述。這里我們假設(shè)通過(guò)SMO算法確定了最優(yōu)α?α?,則

?

w?=∑i=1nαi?yixiw?=∑i=1nαi?yixi

?

最后使用一個(gè)正的支持向量xsxs,就可以計(jì)算出bb:

?

b?=1?w?Txsb?=1?w?Txs

?

1.3. 軟間隔

在4.2節(jié)中我們推導(dǎo)了如何計(jì)算ww、bb和αα,但別忘了以上所有的推導(dǎo)都是在線性可分的條件下進(jìn)行的,但是現(xiàn)實(shí)世界的許多問(wèn)題并不都是線性可分的,尤其存在許多復(fù)雜的非線性可分的情形。如果樣本不能被完全線性分開(kāi),那么情況就是:間隔為負(fù),原問(wèn)題的可行域?yàn)榭?#xff0c;對(duì)偶問(wèn)題的目標(biāo)函數(shù)無(wú)限,這講導(dǎo)致相應(yīng)的最優(yōu)化問(wèn)題不可解。
要解決這些不可分問(wèn)題,一般有兩種方法。第一種是放寬過(guò)于嚴(yán)格的間隔,構(gòu)造軟間隔。另一種是運(yùn)用核函數(shù)把這些數(shù)據(jù)映射到另一個(gè)維度空間去解決非線性問(wèn)題。在本節(jié)中,我們首先介紹軟間隔優(yōu)化。
假設(shè)兩個(gè)類(lèi)有幾個(gè)數(shù)據(jù)點(diǎn)混在一起,這些點(diǎn)對(duì)最優(yōu)超平面形成了噪聲干擾,軟間隔就是要擴(kuò)展一下我們的目標(biāo)函數(shù)和KKT條件,允許少量這樣的噪聲存在。具體地說(shuō),就要引入松馳變量ξiξi來(lái)量化分類(lèi)器的違規(guī)行為:

?

???????minw,b12||w||2+C∑ni=1ξiyi(wTxi+b)?1?ξi,ξi?0,(i=1,..,n){minw,b12||w||2+C∑i=1nξiyi(wTxi+b)?1?ξi,ξi?0,(i=1,..,n)

?

參數(shù)C用來(lái)平衡機(jī)器的復(fù)雜度和不可分?jǐn)?shù)據(jù)點(diǎn)的數(shù)量,它可被視為一個(gè)由用戶(hù)依據(jù)經(jīng)驗(yàn)或分析選定的“正則化”參數(shù)。松馳變量ξiξi的一個(gè)直接的幾何解釋是一個(gè)錯(cuò)分實(shí)例到超平面的距離,這個(gè)距離度量的是錯(cuò)分實(shí)例相對(duì)于理想的可分模式的偏差程度。對(duì)上式的化解,可得:

?

W(α)=∑i=1nαi?12∑i=1n∑j=1nαiαjyiyjxiTxjW(α)=∑i=1nαi?12∑i=1n∑j=1nαiαjyiyjxiTxj

?

?

∑i=1nαiyi=0,0?αi?C(i=1,...,n)∑i=1nαiyi=0,0?αi?C(i=1,...,n)

?

可以看到,松馳變量ξiξi沒(méi)有出現(xiàn)在W(α)W(α)中,線性可分與不可分的差異體現(xiàn)在約束αi?0αi?0被替換成了約束0?αi?C0?αi?C。但是,這兩種情況下求解ww和bb是非常相似的,對(duì)于支持向量的定義也都是一致的。
在不可分情況下,對(duì)應(yīng)的KKT條件為:

?

αi[yi(wTxi+b)?1+ξi]=0,(i=1,...,n)αi[yi(wTxi+b)?1+ξi]=0,(i=1,...,n)

?

1.4. SMO算法

1996年, John Platt發(fā)布了一個(gè)稱(chēng)為SMO的強(qiáng)大算法,用于訓(xùn)練SVM。 SMO表示序列最小優(yōu)化(Sequential Minimal Optimization)。 Platt的SMO算法是將大優(yōu)化問(wèn)題分解為多個(gè)小優(yōu)化問(wèn)題來(lái)求解,這些小優(yōu)化問(wèn)題往往很容易求解,并且對(duì)它們進(jìn)行順序求解的結(jié)果與將它們作為整體來(lái)求解的結(jié)果是完全一致的。
SMO算法的目標(biāo)是求出一系列αα,一旦求出了這些αα,就很容易計(jì)算出權(quán)重向量ww和bb,并得到分隔超平面。
SMO算法的工作原理是:每次循環(huán)中選擇兩個(gè)αα進(jìn)行優(yōu)化處理。一旦找到一對(duì)合適的αα,那么就增大其中一個(gè)同時(shí)減小另一個(gè)。這里所謂的“合適”就是指兩個(gè)αα必須要符合一定的條件,條件之一就是這兩個(gè)αα必須要在間隔邊界之外,而其第二個(gè)條件則是這兩個(gè)αα還沒(méi)有進(jìn)行過(guò)區(qū)間化處理或者不在邊界上。
對(duì)SMO具體的分析如下,在4.3節(jié)中我們已經(jīng)得出了

?

W(α)=∑i=1nαi?12∑i=1n∑j=1nαiαjyiyjxiTxjW(α)=∑i=1nαi?12∑i=1n∑j=1nαiαjyiyjxiTxj

?

?

∑i=1nαiyi=0,0?αi?C(i=1,...,n)∑i=1nαiyi=0,0?αi?C(i=1,...,n)

?

其中(xi,yi)(xi,yi)已知,C可以預(yù)先設(shè)定,也是已知數(shù),現(xiàn)在就是要最大化W(α)W(α),求得參數(shù)α=[α1,α2,...,αn]α=[α1,α2,...,αn]。SMO算法是一次選擇兩個(gè)αα進(jìn)行優(yōu)化,那我們就選擇α1α1和α2α2,然后把其他參數(shù)[α3,α4,...,αn][α3,α4,...,αn]固定,這樣α1α1、α2α2表示為下面的式子,其中ζζ是實(shí)數(shù)值:

?

α1y1+α2y2=?∑i=3nαiyi=ζα1y1+α2y2=?∑i=3nαiyi=ζ

?

然后用α2α2來(lái)表示α1α1:

?

α1=(ζ?α2y2)y1α1=(ζ?α2y2)y1

?

把上式帶入W(α)W(α)中:

?

W(α)=W(α1,α2,...,αn)=W((ζ?α2y2)y1,α2,...,αn)W(α)=W(α1,α2,...,αn)=W((ζ?α2y2)y1,α2,...,αn)

?

省略一系列化解過(guò)程后,最后會(huì)化解成我們熟悉的一元二次方程,a,b,c均是實(shí)數(shù)值:

?

W(α2)=aα22+bα2+cW(α2)=aα22+bα2+c

?

最后對(duì)α2α2求導(dǎo),解得α2α2的具體值,我們暫時(shí)把這個(gè)實(shí)數(shù)值叫α?2α2?。而這個(gè)α?2α2?需要滿(mǎn)足一個(gè)條件L?α?2?HL?α2??H,其中LL和HH是什么呢?如下圖所示:


(圖片來(lái)自網(wǎng)絡(luò))

根據(jù)之前的條件0?αi?C0?αi?C和等式α1y1+α2y2=ζα1y1+α2y2=ζ知α1α1和α2α2要在矩形區(qū)域內(nèi),并且在直線上。當(dāng)y1y1和y2y2異號(hào)時(shí):

?

???L=max(0,α2?α1)H=min(C,C+α2?α1){L=max(0,α2?α1)H=min(C,C+α2?α1)

?

當(dāng)y1y1和y2y2同號(hào)時(shí):

?

???L=max(0,α2+α1?C)H=min(C,α2+α1){L=max(0,α2+α1?C)H=min(C,α2+α1)

?

最后,滿(mǎn)足條件的α2α2應(yīng)該由下面的式子得到,α??2α2??才為最終的值:

?

α??2=?????????????Hα?2L,α?2>H,L≤α?2≤H,α?2<Lα2??={H,α2?>Hα2?,L≤α2?≤HL,α2?<L

?

求得α??2α2??后我們就可以求得α??1α1??了。然后我們重復(fù)地按照最優(yōu)化(α1,α2)(α1,α2)的方式繼續(xù)選擇(α3,α4)(α3,α4),(α5,α6)(α5,α6)....(αn?1,αn)(αn?1,αn)進(jìn)行優(yōu)化求解,這樣α=[α1,α2,...,αn]α=[α1,α2,...,αn]求解出來(lái)后,整個(gè)線性劃分問(wèn)題就迎刃而解。

1.5. 核函數(shù)

對(duì)于以上幾節(jié)講的SVC算法,我們都在線性可分或存在一些噪聲點(diǎn)的情況下進(jìn)行的二分類(lèi),但是如果我們存在兩組數(shù)據(jù),它們的散點(diǎn)圖如下圖所示,你可以看出這完全是一個(gè)非線性不可分的問(wèn)題,我們無(wú)法使用之前講的SVC算法在這個(gè)二維空間中找到一個(gè)超平面把這些數(shù)據(jù)點(diǎn)準(zhǔn)確的分開(kāi)。

?

解決這個(gè)劃分問(wèn)題我們需要引入一個(gè)核函數(shù),核函數(shù)能夠恰當(dāng)?shù)挠?jì)算給定數(shù)據(jù)的內(nèi)積,將數(shù)據(jù)從輸入空間的非線性轉(zhuǎn)變到特征空間,特征空間具有更高甚至無(wú)限的維度,從而使得數(shù)據(jù)在該空間中被轉(zhuǎn)換成線性可分的。如下圖所示,我們把二維平面的一組數(shù)據(jù),通過(guò)核函數(shù)映射到了一個(gè)三維空間中,這樣,我們的超平面就面成了一個(gè)平面(在二維空間中是一條直線),這個(gè)平面就可以準(zhǔn)確的把數(shù)據(jù)劃分開(kāi)了。

?

核函數(shù)有Sigmoid核、線性核、多項(xiàng)式核和高斯核等,其中高斯核和多項(xiàng)式核比較常用,兩種核函數(shù)均可以把低維數(shù)據(jù)映射到高維數(shù)據(jù)。高斯核的公式如下,σσ是達(dá)到率,即函數(shù)值跌落到0的速度參數(shù):

?

K(x1,x2)=exp(?||x1?x2||22σ2)K(x1,x2)=exp(?||x1?x2||22σ2)

?

多項(xiàng)式核函數(shù)的公式如下,RR為實(shí)數(shù),dd為低維空間的維數(shù):

?

K(x1,x2)=(?x1,x2?+R)dK(x1,x2)=(?x1,x2?+R)d

?

應(yīng)用于我們的上個(gè)例子,我們先定義,用?:x→H?:x→H表示從輸入空間x?Rnx?Rn到特征空間H的一個(gè)非線性變換。假設(shè)在特征空間中的問(wèn)題是線性可分的,那么對(duì)應(yīng)的最優(yōu)超平面為:

?

w?T?(x)+b=0w?T?(x)+b=0

?

通過(guò)拉格朗日函數(shù)我們推導(dǎo)出:

?

w??=∑i=1nαi?yi?(xi)w??=∑i=1nαi?yi?(xi)

?

帶入上式得特征空間的最優(yōu)超平面為:

?

∑i=1nαi?yi?T(xi)?(x)+b=0∑i=1nαi?yi?T(xi)?(x)+b=0

?

這里的?T(xi)?(x)?T(xi)?(x)表示內(nèi)積,用核函數(shù)代替內(nèi)積則為:

?

∑i=1nαi?yiK(xi,x)+b=0∑i=1nαi?yiK(xi,x)+b=0

?

這說(shuō)明,我們的核函數(shù)均是內(nèi)積函數(shù),通過(guò)在低維空間對(duì)輸入向量求內(nèi)積來(lái)映射到高維空間,從而解決在高維空間中數(shù)據(jù)線性可分的問(wèn)題,至于具體的推導(dǎo)過(guò)程,這里就不再進(jìn)行了,感興趣的可以自己再推導(dǎo)一次,加深理解。
為什么核函數(shù)可以把低維數(shù)據(jù)映射成高維數(shù)據(jù)呢,我們以多項(xiàng)式核來(lái)解釋一下。
假設(shè)有兩個(gè)輸入樣本,它們均為二維行向量x1=[x1,x2]x1=[x1,x2],x2=[x3,x4]x2=[x3,x4],他們的內(nèi)積為:

?

?x1,x2?=x1x2T=[x1x2][x3x4]=x1x3+x2x4?x1,x2?=x1x2T=[x1x2][x3x4]=x1x3+x2x4

?

用多項(xiàng)式核函數(shù)進(jìn)行映射,令R=0R=0,d=2d=2:

?

K(x1,x2)=(?x1,x2?)2=(x1x3+x2x4)2=x12x32+2x1x2x3x4+x22x42=?(x1)?(x2)K(x1,x2)=(?x1,x2?)2=(x1x3+x2x4)2=x12x32+2x1x2x3x4+x22x42=?(x1)?(x2)

?

按照線性代數(shù)中的標(biāo)準(zhǔn)定義,?(x1)?(x1)和?(x2)?(x2)為映射后的三維行向量和三維列向量,即:

?

?(x1)=[x122–√x1x2x22]?(x1)=[x122x1x2x22]

?

?

?(x2)=????????x322–√x3x4x42?????????(x2)=[x322x3x4x42]

?

它們的內(nèi)積用向量的方式表示則更直觀:

?

?(x1)?(x2)=[x122–√x1x2x22]????????x322–√x3x4x42????????=x12x32+2x1x2x3x4+x22x42?(x1)?(x2)=[x122x1x2x22][x322x3x4x42]=x12x32+2x1x2x3x4+x22x42

?

這樣我們就把二維數(shù)據(jù)映射成了三維數(shù)據(jù),對(duì)于高斯核的映射,會(huì)用到泰勒級(jí)數(shù)展開(kāi)式,讀者可以自行推導(dǎo)一下。
對(duì)于核函數(shù)我們就暫時(shí)介紹到這里。下一節(jié)我們開(kāi)始運(yùn)用tensorflow來(lái)進(jìn)行實(shí)戰(zhàn),由于沒(méi)有找到線性不可分的數(shù)據(jù)集,我們的例子中就沒(méi)有用到核函數(shù)來(lái)建模,因?yàn)槲覀冎徽业搅艘粋€(gè)線性可分的數(shù)據(jù)集,所以下一節(jié)我們僅運(yùn)用tensorflow來(lái)進(jìn)行線性二分類(lèi)的分類(lèi)器構(gòu)建。

1.6. 實(shí)例

我們?cè)诰W(wǎng)上下載了一組鳶尾花數(shù)據(jù)集,這組數(shù)據(jù)集有100個(gè)樣本點(diǎn),我們用SVM來(lái)預(yù)測(cè)這些鳶尾花數(shù)據(jù)集中哪些是山鳶尾花,哪些是非山鳶尾花。

  • 首先需要加載數(shù)據(jù)集,加載數(shù)據(jù)集需要用到sklearn、scipy、mkl庫(kù),sklearn直接在Pycharm中安裝即可,而另外兩個(gè)庫(kù)需要在網(wǎng)上下載安裝,下載地址:http://www.lfd.uci.edu/~gohlke/pythonlibs/
import matplotlib.pyplot as plt import numpy as np import tensorflow as tf from sklearn import datasetssess = tf.Session()# 加載數(shù)據(jù) # iris.data = [(Sepal Length, Sepal Width, Petal Length, Petal Width)] iris = datasets.load_iris() x_vals = np.array([[x[0], x[3]] for x in iris.data]) y_vals = np.array([1 if y == 0 else -1 for y in iris.target])
  • 分離測(cè)試集與訓(xùn)練集
# 分離訓(xùn)練和測(cè)試集 train_indices = np.random.choice(len(x_vals),round(len(x_vals)*0.8),replace=False) test_indices = np.array(list(set(range(len(x_vals))) - set(train_indices))) x_vals_train = x_vals[train_indices] x_vals_test = x_vals[test_indices] y_vals_train = y_vals[train_indices] y_vals_test = y_vals[test_indices]
  • 定義模型和loss函數(shù)
batch_size = 100# 初始化feedin x_data = tf.placeholder(shape=[None, 2], dtype=tf.float32) y_target = tf.placeholder(shape=[None, 1], dtype=tf.float32)# 創(chuàng)建變量 A = tf.Variable(tf.random_normal(shape=[2, 1])) b = tf.Variable(tf.random_normal(shape=[1, 1]))# 定義線性模型 model_output = tf.subtract(tf.matmul(x_data, A), b)# Declare vector L2 'norm' function squared l2_norm = tf.reduce_sum(tf.square(A))# Loss = max(0, 1-pred*actual) + alpha * L2_norm(A)^2 alpha = tf.constant([0.01]) classification_term = tf.reduce_mean(tf.maximum(0., tf.subtract(1., tf.multiply(model_output, y_target)))) loss = tf.add(classification_term, tf.multiply(alpha, l2_norm))
  • 開(kāi)始訓(xùn)練數(shù)據(jù)
my_opt = tf.train.GradientDescentOptimizer(0.01) train_step = my_opt.minimize(loss)init = tf.global_variables_initializer() sess.run(init)# Training loop loss_vec = [] train_accuracy = [] test_accuracy = [] for i in range(20000):rand_index = np.random.choice(len(x_vals_train), size=batch_size)rand_x = x_vals_train[rand_index]rand_y = np.transpose([y_vals_train[rand_index]])sess.run(train_step, feed_dict={x_data: rand_x, y_target: rand_y})
  • 繪制圖像
[[a1], [a2]] = sess.run(A) [[b]] = sess.run(b) slope = -a2/a1 y_intercept = b/a1 best_fit = []x1_vals = [d[1] for d in x_vals]for i in x1_vals:best_fit.append(slope*i+y_intercept)# Separate I. setosa setosa_x = [d[1] for i, d in enumerate(x_vals) if y_vals[i] == 1] setosa_y = [d[0] for i, d in enumerate(x_vals) if y_vals[i] == 1] not_setosa_x = [d[1] for i, d in enumerate(x_vals) if y_vals[i] == -1] not_setosa_y = [d[0] for i, d in enumerate(x_vals) if y_vals[i] == -1]plt.plot(setosa_x, setosa_y, 'o', label='I. setosa') plt.plot(not_setosa_x, not_setosa_y, 'x', label='Non-setosa') plt.plot(x1_vals, best_fit, 'r-', label='Linear Separator', linewidth=3) plt.ylim([0, 10]) plt.legend(loc='lower right') plt.title('Sepal Length vs Pedal Width') plt.xlabel('Pedal Width') plt.ylabel('Sepal Length') plt.show()

訓(xùn)練后的結(jié)果

轉(zhuǎn)載注明出處http://www.cnblogs.com/vipyoumay/p/7560061.html

以上內(nèi)容有任何錯(cuò)誤或不準(zhǔn)確的地方請(qǐng)大家指正,不喜勿噴! 本文版權(quán)歸作者和博客園共有,歡迎轉(zhuǎn)載,但未經(jīng)作者同意必須保留此段聲明,且在文章頁(yè)面明顯位置給出原文連接,否則保留追究法律責(zé)任的權(quán)利。如果覺(jué)得還有幫助的話,可以點(diǎn)一下右下角的【推薦】,希望能夠持續(xù)的為大家?guī)?lái)好的技術(shù)文章!想跟我一起進(jìn)步么?那就【關(guān)注】我吧。

來(lái)源:https://www.cnblogs.com/vipyoumay/p/7560061.html

總結(jié)

以上是生活随笔為你收集整理的SVM原理以及Tensorflow 实现SVM分类(附代码讲了一下原理)的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。

久久久久久久久久久免费av | 久久久久久高清 | 9992tv成人免费看片 | 日韩乱码在线 | 在线免费色视频 | 久久夜色精品国产欧美一区麻豆 | www.久久99| 国产特级毛片aaaaaa毛片 | 亚洲精品久久久久中文字幕二区 | 精品国产视频一区 | 国产在线高清视频 | 日韩成人免费在线 | 中文字幕在线观看第二页 | 久久精品国产一区二区电影 | 色在线高清 | 免费视频二区 | 免费观看丰满少妇做爰 | 97超级碰碰碰碰久久久久 | 天天天插| 成人一区二区在线 | 日韩av免费在线电影 | 精品美女久久 | 中文字幕在线看视频 | 国产资源在线视频 | 久久久亚洲麻豆日韩精品一区三区 | a级免费观看 | 久久99日韩 | 久久色视频 | 丁香花五月 | 天天综合入口 | 国产精品99久久久久人中文网介绍 | 免费在线一区二区三区 | 亚洲精品xxx| 日韩激情视频在线 | 国产69久久久欧美一级 | 国产精品久久嫩一区二区免费 | 中文字幕av全部资源www中文字幕在线观看 | 在线看小早川怜子av | 日韩sese | 色网站在线免费 | 亚洲国产经典视频 | 久久国产精彩视频 | 91成人黄色 | 天天爽天天做 | 婷婷久月| 日韩电影久久 | 狠狠色婷婷丁香六月 | 人人干人人做 | 91精品无人成人www | 日韩免费观看视频 | 亚洲日本精品视频 | 免费在线看v | 中文字幕网址 | 三三级黄色片之日韩 | 一区二区 不卡 | 337p西西人体大胆瓣开下部 | 国产精品 9999| 国产精品欧美精品 | 91久久精品日日躁夜夜躁国产 | 欧美精品久久 | 国产精品 中文字幕 亚洲 欧美 | 在线观看一区二区精品 | 婷婷色六月天 | 日韩在线电影 | 久久不色| av解说在线 | 视频一区二区国产 | 国产午夜三级 | 精品久久久久久亚洲综合网站 | 久久图 | 国产在线免费观看 | 亚洲视频第一页 | 午夜色站 | 97超碰国产精品女人人人爽 | 97视频在线观看免费 | 99久久精品视频免费 | 亚洲另类久久 | 美女精品久久久 | 亚洲国产一区av | 五月天六月丁香 | 最近免费观看的电影完整版 | 久久国产精品免费看 | 久久在线精品 | 日韩午夜电影院 | 国产高清在线不卡 | 久久免费高清视频 | 精品伦理一区二区三区 | 激情深爱| 最新色站 | 福利一区在线 | 国产成人一区二区三区电影 | 综合色中文| 中文电影网 | 在线观看91精品国产网站 | 亚洲 欧洲 国产 日本 综合 | 久久av一区二区三区亚洲 | a爱爱视频| 国产视频丨精品|在线观看 国产精品久久久久久久久久久久午夜 | wwwwwww黄| 国产亚洲视频在线观看 | av色网站| 国产亚洲欧美精品久久久久久 | 五月天久久 | 国产精品1区2区3区在线观看 | 亚洲免费成人av电影 | 国产综合婷婷 | 五月婷婷网站 | 九月婷婷人人澡人人添人人爽 | 国产中文字幕视频在线观看 | 天天操夜夜干 | 91亚洲永久精品 | 国产美女搞久久 | 久久一精品 | 九七视频在线观看 | 中文字幕 国产精品 | 亚洲综合色激情五月 | 天天天色综合 | 日韩午夜小视频 | 91av免费看 | 成人精品国产免费网站 | 97天天综合网 | 黄在线免费看 | 国产色区 | 免费99精品国产自在在线 | 人人草人人草 | 久久综合五月婷婷 | 麻豆成人网 | 亚洲精品毛片一级91精品 | 国产精品一区二区视频 | 中文在线a√在线 | 四川妇女搡bbbb搡bbbb搡 | 日韩在线视频播放 | 日韩高清片 | 亚洲经典精品 | 色综合久久久久综合99 | 国产成人久久精品77777综合 | 日韩三级av| 97色婷婷成人综合在线观看 | 欧美一区二视频在线免费观看 | 午夜久久久久久久久 | 久久99久久99精品免观看软件 | 国产免费xvideos视频入口 | 黄色av电影网 | 久久免费一 | 在线亚州| 九九九九热精品免费视频点播观看 | 国产美腿白丝袜足在线av | 在线观看黄色的网站 | 免费久久99精品国产 | 国产亚洲91 | 不卡的av电影在线观看 | 日韩在线视频二区 | 亚洲精品久久在线 | 欧美色图88 | 免费久久久久久 | 亚洲国产成人精品久久 | 91成人精品一区在线播放69 | 一级全黄毛片 | 国产96av | 免费观看一区 | 白丝av免费观看 | 久草网在线视频 | 国产裸体永久免费视频网站 | 精品一区二区电影 | 国产不卡在线观看 | 国产精品女人久久久久久 | 久久久久久蜜桃一区二区 | 婷婷夜夜 | 97视频网址 | 天天曰夜夜爽 | 国产福利中文字幕 | 午夜精品福利影院 | 日韩精品在线观看av | 天天操天天干天天插 | 一区二区三区精品在线视频 | 四虎成人网| 在线观看视频你懂 | 国产精品九色 | 99久久日韩精品免费热麻豆美女 | 日韩精品中文字幕一区二区 | 国产一区视频在线播放 | 美女视频一区二区 | 成人国产精品免费 | 99精品在线观看 | 国产一级不卡毛片 | 热99在线视频| 91一区啪爱嗯打偷拍欧美 | 日韩精品大片 | 国产精品每日更新 | av千婊在线免费观看 | 在线观看你懂的网址 | 偷拍精品一区二区三区 | 久久人人爽人人爽 | 欧美亚洲xxx | 精品亚洲视频在线观看 | 日韩中文字幕免费电影 | 少妇自拍av| 黄色av成人在线观看 | 久久久久久久网站 | 在线观看的av | 欧美黑人xxxx猛性大交 | 成年人在线观看视频免费 | 日韩精品一区不卡 | 97超碰福利久久精品 | 91精品在线免费视频 | 日韩精品中文字幕在线不卡尤物 | 午夜在线免费观看 | 狠狠的干狠狠的操 | 天天干天天拍天天操 | 韩国一区视频 | 国产午夜精品一区二区三区欧美 | 97成人在线观看 | 综合色站| 米奇四色影视 | 久色小说| 狠狠干夜夜爱 | 国产香蕉久久精品综合网 | 亚洲国产av精品毛片鲁大师 | 精品国内自产拍在线观看视频 | 亚洲精品一区二区精华 | 国产精品福利无圣光在线一区 | 久久婷婷精品 | 久久99亚洲精品久久久久 | 日产乱码一二三区别免费 | 免费福利在线观看 | 日韩av一区二区在线播放 | 日韩在线观看一区二区三区 | 久久99精品久久久久久三级 | 九九久久精品 | 国产精品久久久久999 | www在线观看视频 | 国产精品一区二区av影院萌芽 | 中文字幕在线视频网站 | 在线 视频 亚洲 | 久草在线视频免费资源观看 | 中文字幕 国产视频 | 欧美精品你懂的 | 天天av资源 | 日本在线观看中文字幕 | 日韩免费| 亚洲人天堂 | 91精品国产三级a在线观看 | 久久精品99国产国产 | 国产一线天在线观看 | 久久免费片 | 日韩欧美综合精品 | 欧美日韩国产亚洲乱码字幕 | 成片人卡1卡2卡3手机免费看 | 天天干夜夜夜 | 国产亚洲精品中文字幕 | 久久人人爽人人爽人人 | 亚洲天堂网在线播放 | 免费在线精品视频 | 天天操天天操天天 | 天天综合天天做天天综合 | 免费合欢视频成人app | 免费a视频在线 | 黄色大片视频网站 | 亚洲四虎影院 | 国产亚洲精品久久久久秋 | 99久久精品国产观看 | 日日夜夜精品视频天天综合网 | 开心激情五月网 | av网站免费看 | 亚洲精品小区久久久久久 | 国产区在线看 | 黄色大片免费网站 | 久久久精品成人 | 亚洲日韩精品欧美一区二区 | 亚洲精品 在线视频 | 青青草国产在线 | 午夜精品一区二区三区在线观看 | 五月激情五月激情 | 91黄色免费网站 | 日韩综合视频在线观看 | 女人18毛片90分钟 | 丁香午夜婷婷 | 精品伊人久久久 | 欧美日韩在线视频观看 | 在线播放 日韩专区 | 国产午夜精品av一区二区 | 久碰视频在线观看 | 国产成人香蕉 | 日本九九视频 | www.狠狠操.com | 国产精品第72页 | 狠狠婷婷| 伊人五月天av| 日日夜夜网 | av在线亚洲天堂 | 九九在线播放 | 丁香色天天 | 色在线亚洲 | 欧美aaa一级 | 国产精品青青 | 最新日韩在线观看视频 | 午夜国产一区二区 | 97手机电影网 | 欧美性生爱| 免费一级日韩欧美性大片 | 夜夜骑日日操 | 欧美一区二区三区在线播放 | 亚洲精选视频免费看 | av中文字幕在线观看网站 | 欧美一级黄大片 | 黄色三级网站在线观看 | 久久久91精品国产一区二区三区 | 色欧美成人精品a∨在线观看 | 999久久久免费视频 午夜国产在线观看 | 四川bbb搡bbb爽爽视频 | 国产高清福利在线 | 干综合网| 国产视频每日更新 | 中文字幕三区 | 亚洲精品国产自产拍在线观看 | 亚洲美女精品区人人人人 | 91人人澡人人爽人人精品 | 麻豆视频免费入口 | 国内外激情视频 | av午夜电影| 日韩欧美xxx | 波多野结衣在线观看一区二区三区 | 亚洲黄色一级大片 | 粉嫩一二三区 | 国产99久久久国产精品 | 四虎最新域名 | 一区中文字幕电影 | 狠色狠色综合久久 | 在线超碰av| 久久亚洲专区 | 手机在线观看国产精品 | 中日韩欧美精彩视频 | 视频一区二区在线观看 | 天天操天天色天天 | 天天色天 | 久久久国产在线视频 | 91av超碰| 99精品视频免费看 | 国产日本三级 | 欧美日韩国产亚洲乱码字幕 | 在线视频欧美日韩 | 日韩av视屏在线观看 | 久久激情视频 | 欧美日韩国产一区二 | 国内毛片毛片 | 日韩欧美视频在线免费观看 | 亚洲一区二区三区四区精品 | 国产区精品在线观看 | 黄色av免费电影 | 亚洲精品456在线播放第一页 | 成人wwwxxx视频| 天天射天天 | 丁香午夜婷婷 | 色视频在线免费 | 在线观看一级 | 99在线精品视频观看 | 国产在线日本 | 深爱激情开心 | 国产精品一区二区在线观看免费 | 美女国产免费 | 国产精品美女久久久久aⅴ 干干夜夜 | 婷婷六月丁香激情 | 国产精品女同一区二区三区久久夜 | 亚洲mv大片欧洲mv大片免费 | 亚洲视频久久久久 | 91网址在线观看 | 免费在线观看av电影 | 国产精品va在线观看入 | 狠狠干夜夜操天天爽 | 久久人人精 | 国产精品一区二区在线观看免费 | 亚洲精区二区三区四区麻豆 | 日本久久91| 97人人添人澡人人爽超碰动图 | 一区二区三区在线视频111 | 在线日本v二区不卡 | 操久| a视频免费看 | av在线免费在线观看 | 在线国产精品视频 | 免费亚洲精品 | 激情一区二区三区欧美 | 高清一区二区 | 亚洲成av人片在线观看 | 黄色大片免费播放 | 在线国产能看的 | 国产精品av免费 | 中文字幕一区二区三区乱码在线 | 在线观看精品一区 | 国产99久久久国产精品 | 在线观看不卡视频 | 久久在线| 久久久影院一区二区三区 | 日韩网站在线免费观看 | 91视频国产免费 | 国产麻豆精品一区 | 黄色片视频免费 | 欧美日韩中字 | 欧美成人区 | 欧美精品久久久久久久久久白贞 | www.日本色 | 啪啪av在线| 高清在线一区二区 | 国产在线久草 | 激情 一区二区 | 精品美女久久久久久免费 | 99视频在线观看免费 | 99久久精品网 | 狠狠色噜噜狠狠 | 欧美精品久久久久久久久久丰满 | 婷婷5月激情5月 | 久久热首页 | 亚洲精品看片 | 国产一区二区午夜 | 国产专区视频在线 | 黄色精品一区二区 | 深夜国产在线 | 欧美日韩国产二区三区 | 成年人国产在线观看 | 国产资源av | 国产成人精品一区二三区 | 91精品在线免费观看 | 99国产精品久久久久久久久久 | 亚洲视频电影在线 | 91av亚洲 | 色婷婷六月天 | 亚洲天天摸日日摸天天欢 | 四虎在线观看视频 | 最新国产精品亚洲 | 国产黄色播放 | 亚洲精品一区二区在线观看 | 免费在线电影网址大全 | 国内视频在线 | 黄色一级大片免费看 | 免费亚洲黄色 | www.天天色.com | 国产麻豆精品一区二区 | 亚洲精品国产第一综合99久久 | 91爱爱电影 | 99热日本| 国产精品嫩草影院9 | 久久狠狠干 | 久久国产品 | 国产日韩在线一区 | 中文字幕在线第一页 | 91免费观看国产 | 亚洲爱视频 | 日韩av成人 | 国产99精品 | 午夜丁香网 | 绯色av一区 | 91av在线视频免费观看 | 69久久夜色精品国产69 | 成人三级av | 欧美地下肉体性派对 | 免费黄色在线网站 | 中文字幕电影一区 | 国产精品免费看久久久8精臀av | 国产美女免费 | 在线精品亚洲 | 欧美一区中文字幕 | 午夜精品久久久 | 韩国三级一区 | 在线观看的av | 久久精品香蕉 | 在线免费看黄色 | 亚洲精品国产精品国自产在线 | 91成人在线观看高潮 | 久草在线高清视频 | av大全免费在线观看 | 91在线视频免费观看 | www看片网站 | 91九色九色 | 免费观看视频的网站 | 91精品国产91久久久久 | 色播五月婷婷 | 天天操天天操一操 | 亚洲精品xxxx| 精品人人人人 | 久久国产手机看片 | 欧美日韩中文字幕在线视频 | 国产在线观看,日本 | 丝袜美女在线 | 亚洲精品乱码久久久久久9色 | 国产成人一区二区啪在线观看 | 黄网站免费看 | 免费手机黄色网址 | 精品久久久久久国产 | 在线黄色免费 | 在线播放国产一区二区三区 | 国产香蕉久久 | 免费av大片 | 免费在线观看成人小视频 | 中文亚洲欧美日韩 | 国产涩涩网站 | 日韩精品一区二区在线观看视频 | 在线观看中文字幕一区 | 91精品视频播放 | 欧美精品你懂的 | 2024国产精品视频 | 亚洲精品xxxx| av电影在线免费观看 | 欧美日韩精品免费观看视频 | 韩国精品在线 | 中文字幕一区二区三区四区 | 999久久久久久久久 69av视频在线观看 | 久草网视频在线观看 | 欧美激情综合五月色丁香 | 日本mv大片欧洲mv大片 | 亚洲精品女人久久久 | 国产成人免费在线 | 中文字幕乱码日本亚洲一区二区 | 久草剧场| 国产五月婷 | 欧美激情片在线观看 | 2000xxx影视 | 久久成人视屏 | 不卡中文字幕av | 午夜精品影院 | 日日干综合 | 国产精品久久久久高潮 | 国产男女无遮挡猛进猛出在线观看 | 99爱视频 | 免费成人av在线看 | 久久蜜桃av | 一区二区三区中文字幕在线 | 日韩在线观看视频在线 | 精品99久久 | 久久手机免费观看 | 国产 一区二区三区 在线 | 五月天婷婷在线观看视频 | 美女视频免费精品 | 激情五月婷婷综合网 | 三级免费黄 | 最新高清无码专区 | 成人av在线资源 | 久久国产视屏 | 中文字幕黄色 | 国产大尺度视频 | 欧美99精品| 亚a在线 | 在线亚洲欧美视频 | 91麻豆精品国产91久久久无限制版 | 一本大道久久精品懂色aⅴ 五月婷社区 | 国产97视频 | 色午夜 | 精品产品国产在线不卡 | 日本深夜福利视频 | 国产精品国产三级国产aⅴ9色 | 一 级 黄 色 片免费看的 | 久久另类小说 | 国产伦精品一区二区三区四区视频 | 国产专区欧美专区 | 在线视频精品 | 成年人免费看片网站 | 99视频免费播放 | av中文字幕网 | 国产丝袜一区二区三区 | 成人超碰97| 成人午夜电影网 | 久久理论视频 | 国产福利不卡视频 | 国产日产亚洲精华av | 亚洲精品无 | 国产美女视频免费观看的网站 | 国产精品久久久久久一区二区 | 在线观看av的网站 | 天天性天天草 | 丁香激情综合国产 | 97偷拍视频| 视频国产在线 | 中文字幕高清在线播放 | 色之综合网 | 免费在线观看一区二区三区 | 日韩精品一区电影 | 日韩精品综合在线 | 免费在线观看亚洲视频 | www.狠狠干 | 国产黄色在线网站 | 中文字幕免费在线 | 亚洲人成人99网站 | 97视频在线观看视频免费视频 | 国产网站在线免费观看 | 国产精品第一页在线观看 | 一区二区三区四区精品视频 | 美女免费网视频 | 免费在线观看一级片 | 天天综合区| www五月天com | 国产99久久九九精品免费 | 色99中文字幕 | 成人免费xxx在线观看 | 亚洲精品久久久久久久不卡四虎 | 美女国产网站 | 国产一区视频免费在线观看 | 99久久综合狠狠综合久久 | 成人av在线网址 | 亚洲精品福利在线观看 | 亚洲高清在线精品 | 91视频 - 88av| 美女在线免费视频 | 成人a在线 | 日韩中字在线 | 夜夜操夜夜干 | 高清一区二区三区av | 久久精品欧美一区 | 夜夜躁狠狠躁日日躁视频黑人 | 四虎成人精品在永久免费 | 97精品国产 | 手机av看片| 九九热免费视频在线观看 | 人人爽人人爽 | 99久热精品 | 看片一区二区三区 | 国产免费亚洲 | 久久久国际精品 | 在线观看日韩精品 | 久久久免费国产 | 99精品在线视频观看 | 免费看av在线| 国产自偷自拍 | 亚洲一区精品二人人爽久久 | 国产视频一区二区在线观看 | 超碰久热 | 久久久久久久久久久久久久电影 | 久久免费播放视频 | 欧洲精品视频一区 | 国产亚洲综合精品 | 国产亚州精品视频 | 日日夜av| 国产精品久久久久一区二区三区共 | 日本性xxxxx 亚洲精品午夜久久久 | 日韩午夜在线观看 | 午夜久久久久 | 热久久这里只有精品 | 中文字幕国产一区二区 | 精品久久久久一区二区国产 | 国产理伦在线 | 精品在线视频一区 | 337p日本大胆噜噜噜噜 | 麻豆影视在线播放 | 97视频免费 | 亚洲 欧美 变态 国产 另类 | 国产资源在线免费观看 | 一区二区成人国产精品 | 手机看片久久 | 一级黄色片毛片 | 五月天综合色激情 | 97人人模人人爽人人喊中文字 | 日韩,中文字幕 | 婷婷亚洲最大 | av韩国在线 | 探花视频免费在线观看 | 在线免费视频你懂的 | 国产精品mm | 国产精品一区二区久久久 | 狠狠色丁香婷综合久久 | 国产精品a久久 | 国产免费区 | 亚洲欧洲国产日韩精品 | 一本一本久久a久久 | 国产精品 日韩 欧美 | 欧美在线你懂的 | 日韩网站免费观看 | 亚洲国产综合在线 | 国产一区欧美二区 | 天天射色综合 | 亚洲精品国产精品乱码在线观看 | 亚洲狠狠婷婷 | 日韩二区三区在线观看 | www.xxxx欧美| 在线免费观看av网站 | 天天爽天天爽天天爽 | 午夜婷婷在线观看 | 日韩中字在线 | 少妇精品久久久一区二区免费 | 精品国产三级a∨在线欧美 免费一级片在线观看 | 欧美成人黄 | 欧美一区二区日韩一区二区 | 在线看黄色的网站 | 国产在线一区二区三区播放 | 波多野结衣一区二区 | 欧美在线观看视频 | 亚洲精品小区久久久久久 | 麻豆mv在线观看 | 亚洲精欧美一区二区精品 | 日韩免费电影一区二区 | 天天操夜夜操天天射 | 国产伦理精品一区二区 | 丁香综合| 狠狠干在线播放 | 99久久夜色精品国产亚洲96 | 麻豆激情电影 | 91刺激视频 | 欧美一级高清片 | 国产精品99免视看9 国产精品毛片一区视频 | 亚洲视频综合在线 | 免费a级黄色毛片 | 成人av影视在线 | 欧美一区二区三区在线视频观看 | a国产精品 | 黄色一级大片在线观看 | 国产精品成人免费一区久久羞羞 | av 一区二区三区 | 国产精品永久久久久久久久久 | 三级av在线播放 | 日韩黄色中文字幕 | 激情欧美xxxx| www.五月婷婷 | 天天天插| 久久亚洲欧美日韩精品专区 | 日本韩国欧美在线观看 | 97在线播放 | 日韩欧美视频二区 | 在线观看蜜桃视频 | 国产黄色大片 | 九九视频精品免费 | 亚洲动漫在线观看 | 免费av黄色 | 久久精品99久久久久久 | 久久亚洲欧美日韩精品专区 | 久久66热这里只有精品 | 国产免费黄视频在线观看 | 99久久www免费| 天天干天天摸 | 青青久视频| 国产91影视 | 国产精品久久久久久久久久免费 | 色综合久久88色综合天天 | 97超碰在线久草超碰在线观看 | 在线观看免费版高清版 | 亚洲在线视频观看 | 在线视频观看你懂的 | 日韩中文字幕在线 | 中文字幕在线一区二区三区 | 69国产精品成人在线播放 | 日本二区三区在线 | 天天躁日日躁狠狠躁av麻豆 | 成人看片 | 免费视频久久久久久久 | 欧美日韩一区二区在线观看 | a爱爱视频 | 色国产精品一区在线观看 | 成人资源站 | 中文字幕乱码一区二区 | 久久草| 99久久精品国产网站 | 一级欧美一级日韩 | 激情五月婷婷激情 | 97香蕉久久国产在线观看 | 国产精品私拍 | 91精品国产99久久久久 | 国产一区二区在线播放视频 | 亚洲免费av在线播放 | 超碰免费av| 日韩午夜精品 | 国产中年夫妇高潮精品视频 | 婷婷午夜天| 人人插人人 | 欧美资源| 特级西西444www高清大视频 | 久久国产欧美日韩精品 | 国产成人精品综合久久久 | 国产精品毛片久久久久久久久久99999999 | 婷婷激情站 | 国产99久久精品 | www视频在线播放 | 大型av综合网站 | 日韩专区av | 欧美狠狠操 | 国产又粗又猛又黄又爽 | 中文字幕在线看视频 | 在线观看视频你懂 | 国产精品久久久久久久久免费看 | 九九久久电影 | 国产亚洲精品综合一区91 | 免费亚洲一区二区 | 天天天操天天天干 | 国内少妇自拍视频一区 | 国产伦理精品一区二区 | 免费网站黄色 | av黄在线播放 | 99久久精品久久久久久清纯 | 中文字幕一区二区在线播放 | 日本婷婷色| 亚洲精品久久久久久久不卡四虎 | 日本大尺码专区mv | 日韩精品一区二区不卡 | 久久香蕉国产精品麻豆粉嫩av | 在线观看中文av | 女人18片毛片90分钟 | av中文字幕网站 | 99国产在线视频 | 欧洲精品久久久久毛片完整版 | 最新中文字幕在线资源 | 精品一区免费 | 日韩高清av在线 | 日韩久久精品一区二区三区下载 | 视频一区视频二区在线观看 | 久久人人爽视频 | 激情五月婷婷综合 | 黄色一区二区在线观看 | 字幕网av| 精品国模一区二区 | 草樱av| 深爱开心激情 | av色图天堂网 | 免费大片av | 亚洲理论影院 | 在线你懂的视频 | 在线观看中文字幕 | 91看毛片| 国产精品第三页 | 免费av的网站 | 久久久高清一区二区三区 | 国产精品黄色影片导航在线观看 | 涩涩网站在线 | www.五月婷 | 亚洲精品人人 | 黄色录像av | 午夜久久视频 | 韩日成人av | 国产亚洲精品精品精品 | 夜夜躁狠狠燥 | 亚洲国产精品一区二区尤物区 | 一区二区三区免费在线播放 | 欧美国产大片 | 99热在线看 | 亚洲精品系列 | 亚洲精品在线网站 | 黄色成人av在线 | 国产精品日韩 | 黄色av电影网 | 69av网| av在线收看 | 日韩理论片在线观看 | 国产精品永久久久久久久www | 一区二区影院 | 五月天婷亚洲天综合网鲁鲁鲁 | 国产91在线免费视频 | 精品国产免费看 | 久久精品国产精品 | 奇米7777狠狠狠琪琪视频 | 成人久久影院 | 夜夜躁日日躁狠狠久久av | 国产精华国产精品 | 日韩va欧美va亚洲va久久 | 在线日韩精品视频 | 中文字幕在线播放日韩 | 国产精品免费在线观看视频 | 在线观看韩日电影免费 | 精品 一区 在线 | 国产成人精品一区二区在线观看 | 久久精品99久久久久久2456 | 日韩视频免费在线观看 | 国产成人a亚洲精品 | 久久国产精品二国产精品中国洋人 | 五月婷婷久久丁香 | 成人久久视频 | 亚洲国产精品成人va在线观看 | 99性视频 | 久久国产精品久久久 | 久久精品视频在线观看免费 | 狠狠操欧美 | 最近字幕在线观看第一季 | 国产91成人 | 在线天堂日本 | 成人理论在线观看 | 国产亚洲精品电影 | 国产精品6999成人免费视频 | 免费91在线观看 | 国产精品系列在线 | 天天干国产 | 日韩一区在线播放 | 一区二区中文字幕在线播放 | 九九热在线视频 | 国产精品mv在线观看 | 久久午夜视频 | 美女视频一区二区 | 国产精品久久久久久久电影 | 黄色软件在线观看免费 | 国产无遮挡猛进猛出免费软件 | 国产精品原创 | 麻豆精品传媒视频 | 成人免费一区二区三区在线观看 | 午夜久久精品 | 天堂av观看| 色之综合网 | 伊人天天干| 欧美午夜精品久久久久 | 人人看看人人 | 国产一级在线观看视频 | 国产中文字幕一区二区 | 国产精品久久久久久久久久三级 | 亚洲精品视频一二三 | 激情婷婷亚洲 | 黄色日本免费 | 亚洲免费公开视频 | 免费色视频在线 | 91最新地址永久入口 | 五月婷婷av在线 | 久久久久久不卡 | 亚洲成免费 | 99视频精品视频高清免费 | 免费网站观看www在线观看 | 91在线最新 | 人九九精品 | 国产91精品看黄网站 | 欧美激情精品久久久久久变态 | av天天干| 久久久国产影院 | 丁香 久久 综合 | 天天操夜夜干 | 99精品国产在热久久下载 | 人九九精品 | 婷婷综合电影 | 久久精品1区| 久久成人免费 | 在线小视频你懂得 | 色婷婷综合久久久久中文字幕1 | 精品国产美女 | 国产小视频91 | 二区三区视频 | 成人影视片| 亚洲色图 校园春色 | 超碰97成人 | 久久中文字幕视频 | 国产欧美最新羞羞视频在线观看 | 99精品欧美一区二区三区黑人哦 | 欧美性生活免费看 | 久久免费电影 | 97精品国产97久久久久久免费 | 国产在线a | 亚洲国产欧美在线人成大黄瓜 | 欧美精品999 | 国产一区二区高清 | 日韩av成人免费看 | 免费试看一区 | 91精品中文字幕 | 国产日产精品久久久久快鸭 | 国产视频一| 成人福利av| 中文字幕资源网 | 日本少妇视频 | 久久综合狠狠综合久久激情 | 久久九九九九 | 久久久综合香蕉尹人综合网 | 日本精品视频一区 | 99re8这里有精品热视频免费 | 久久成人精品电影 | 国产精品999久久久 久产久精国产品 | 免费色av| 96亚洲精品久久久蜜桃 | 成全在线视频免费观看 | a级国产片 | 一区二区三区av在线 | 一级精品视频在线观看宜春院 | 日韩一级黄色片 | 亚洲国产三级在线观看 | 国产成人精品久久久久蜜臀 | 午夜精品久久久久久久99热影院 | 日韩久久精品 | 激情六月婷婷久久 | 日韩电影一区二区三区在线观看 | 国产香蕉视频在线播放 | 亚洲综合色视频在线观看 | 久久婷婷视频 | 成人少妇影院yyyy | 国产精品久久久毛片 | 91成人国产 | 日韩乱色精品一区二区 | 亚洲精品视频网 | 日韩a在线看 | 久久精品视频在线播放 | 高清av网站 | 国产精品18久久久久久久网站 | 九九九国产| 亚洲成人av一区 | 日本久久精品 | 欧美精品久久久久久久久久久 | 国产成人免费网站 | 欧美日韩性视频在线 | 成人av免费在线播放 | www日韩精品 | 亚洲精品久久久久中文字幕二区 | 亚洲乱码国产乱码精品天美传媒 | 免费视频91蜜桃 | 久久精品一区二区三区国产主播 | 久久久精品久久日韩一区综合 | 999在线精品| 亚洲激情视频 | 国产成人久久精品77777 | 久久精品久久久久久久 | 国产成人精品一区二区三区 | 久久精品999| 去看片 | 夜夜狠狠| 免费黄色在线播放 |