• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于線性回歸模型的顏色與物質(zhì)濃度辨識(shí)應(yīng)用研究

    2019-03-02 02:18:00韋碧鵬黎金清劉劍鋒楊雪
    廣西水利水電 2019年1期
    關(guān)鍵詞:模型

    韋碧鵬,黎金清,劉劍鋒,楊雪

    (柳州職業(yè)技術(shù)學(xué)院,廣西柳州545006)

    1 問題背景

    比色法是目前常用的一種檢測物質(zhì)濃度的方法,即把待測物質(zhì)制備成溶液后滴在特定的白色試紙表面,等其充分反應(yīng)以后獲得一張有顏色的試紙,再把該顏色試紙與一個(gè)標(biāo)準(zhǔn)比色卡進(jìn)行對比,來確定待測物質(zhì)的濃度檔位。由于每個(gè)人對顏色的敏感差異和觀測誤差,這一方法在精度上存在很大的不確定性。隨著照相技術(shù)和顏色分辨率的提高,希望建立顏色讀數(shù)和物質(zhì)濃度的數(shù)量關(guān)系,即只要輸入照片中的顏色讀數(shù)就能夠獲得待測物質(zhì)的濃度。試根據(jù)附件所提供的有關(guān)顏色讀數(shù)和物質(zhì)濃度數(shù)據(jù)完成下列問題:

    (1)附件Data1.xls中分別給出了5種物質(zhì)在不同濃度下的顏色讀數(shù),討論從這5組數(shù)據(jù)中能否確定顏色讀數(shù)和物質(zhì)濃度之間的關(guān)系,并給出一些準(zhǔn)則來評價(jià)這5組數(shù)據(jù)的優(yōu)劣。

    (2)對附件Data2.xls中的數(shù)據(jù),建立顏色讀數(shù)和物質(zhì)濃度的數(shù)學(xué)模型,并給出模型的誤差分析。

    (3)探討數(shù)據(jù)量和顏色維度對模型的影響。

    具體附件請查看全國大學(xué)生數(shù)學(xué)建模競賽官網(wǎng)(http://mcm.blyun.com/),且下載2017年全國大學(xué)生數(shù)學(xué)建模競賽C題進(jìn)行查看。

    2 問題分析

    2.1 問題(1)的分析

    結(jié)合附件1(Data1.xls),要研究5種物質(zhì)在不同的濃度下顏色讀取和物質(zhì)濃度的關(guān)系。由于附件1中給出的5種物質(zhì)在不同濃度下顏色讀數(shù)的數(shù)據(jù)有一組,二組或者多組。首先對5種物質(zhì)溶液下的數(shù)據(jù)進(jìn)行求平均值,接著從不同物質(zhì)溶液的角度出發(fā),分別對顏色讀取和物質(zhì)濃度進(jìn)行多元線性回歸分析。對于物質(zhì)濃度不能進(jìn)行回歸分析的數(shù)據(jù),運(yùn)用相關(guān)分析進(jìn)行分析。最后,從方差和極差的兩個(gè)角度出發(fā)對5組數(shù)據(jù)的優(yōu)劣進(jìn)行評價(jià)。

    2.2 問題(2)的分析

    結(jié)合附件2(Data2.xls),要研究顏色讀數(shù)與物質(zhì)濃度的關(guān)系,采用多元線性回歸模型進(jìn)行擬合,當(dāng)采用多元線性回歸模型中的Enter法進(jìn)行擬合時(shí),發(fā)現(xiàn)所得到的模型不滿足顯著條件。為了能夠更好地研究顏色讀數(shù)和物質(zhì)濃度的數(shù)學(xué)模型,先對附件2中的數(shù)據(jù)進(jìn)行雙變量相關(guān)分析,運(yùn)用Pearson相關(guān)系數(shù),得出二氧化硫的濃度與顏色讀數(shù)R,G,B,H,S的相關(guān)關(guān)系,接著采用多元線性回歸模型中的逐步法進(jìn)行擬合,得出附件2中二氧化硫濃度與顏色讀取的關(guān)系模型。最后,考慮所得出模型的誤差,引入了絕對誤差和相對誤差對其進(jìn)行分析。

    2.3 問題(3)的分析

    為了探討數(shù)據(jù)量對模型的影響,結(jié)合附件2中的數(shù)據(jù),分別從數(shù)據(jù)量多少和顏色維度的多少兩個(gè)角度進(jìn)行分析,做4組數(shù)據(jù)實(shí)驗(yàn)。首先,考慮數(shù)據(jù)量多少對模型的影響,做如下兩組實(shí)驗(yàn)。實(shí)驗(yàn)組A:刪掉二氧化硫濃度為20 ppm與50 ppm的數(shù)據(jù);實(shí)驗(yàn)組B:刪掉二氧化硫濃度為100 ppm與150 ppm的數(shù)據(jù)。通過數(shù)學(xué)軟件SPSS分別對實(shí)驗(yàn)組A、B的數(shù)據(jù)進(jìn)行多元線性回歸,然后用問題(2)所得出原模型的數(shù)據(jù)與實(shí)驗(yàn)組A、B所得出的數(shù)據(jù)進(jìn)行比較,分析數(shù)據(jù)量多少對模型的影響;其次,考慮顏色維度的多少對模型的影響,我們從增加顏色維度和減少顏色維度方向出發(fā),取附件Data2.xls中的數(shù)據(jù)做兩個(gè)實(shí)驗(yàn)組。實(shí)驗(yàn)組C:減少顏色維度,刪除掉H和S兩個(gè)維度的數(shù)據(jù);實(shí)驗(yàn)組D:增加顏色維度I。最后,用問題(2)所得出原模型的數(shù)據(jù)與實(shí)驗(yàn)組C、D所得出的數(shù)據(jù)進(jìn)行比較,分析顏色維度對模型的影響。

    3 模型假設(shè)

    (1)假設(shè)題目所給附件1和附件2中的數(shù)據(jù)無誤;

    (2)假設(shè)各種物質(zhì)的實(shí)驗(yàn)組數(shù)據(jù)都是在同樣的環(huán)境下進(jìn)行的。

    4 模型建立與求解

    4.1 問題(1)模型的建立與求解

    4.1.1 組胺溶液的濃度與顏色讀數(shù)的關(guān)系

    在組胺溶液中,題目附件1給出了兩組同一組胺物質(zhì)濃度的數(shù)據(jù),為了減少實(shí)驗(yàn)的誤差,取兩組實(shí)驗(yàn)數(shù)據(jù)的平均值,得出組胺溶液在不同的濃度下的R,B,H,S值,見表1。

    表1 兩組組胺溶液的RGBHS平均值

    組胺溶液中,為了研究物質(zhì)濃度與顏色讀數(shù)的關(guān)系,對其做多元線性回歸擬合,具體模型如下:

    其中y1表示組胺溶液的濃度,xR、xG、xB、xH、xS分別表示變量紅色、綠色、藍(lán)色、色調(diào)、飽和度,a1、a2、a3、a4、a5分別表示它們的系數(shù),b表示常數(shù)。

    通過數(shù)學(xué)軟件SPSS22,對表1中的數(shù)據(jù)代入進(jìn)行多元線性回歸擬合,運(yùn)用輸入回歸法可得結(jié)果見表2,3。

    表2 變異數(shù)分析a

    表3 系數(shù)a

    從表2,3中可知,模型的顯著性概率為0.00小于0.01,因此,該線性回歸模型有效。又因?yàn)槟P椭谐?shù)、變量G,R,H,S的顯著性概率為0.00均小于0.01,則該模型通過檢驗(yàn),具體模型為:

    4.1.2 溴酸鉀溶液的濃度與顏色讀數(shù)的關(guān)系

    在溴酸鉀溶液中,題目附件1給出了兩組同一溴酸鉀物質(zhì)濃度的數(shù)據(jù),為了減少實(shí)驗(yàn)的誤差,取兩組實(shí)驗(yàn)數(shù)據(jù)的平均值,得出溴酸鉀溶液在不同的濃度下的R,G,B,H,S值。為了研究溴酸鉀溶液濃度與顏色讀數(shù)的關(guān)系,同理,對其做多元線性回歸擬合,通過SPSS22軟件,進(jìn)行多元線性中的輸入回歸可得結(jié)果,從中可知,模型的顯著性概率為0.00小于0.01,因此,該線性回歸模型有效。又因?yàn)槟P椭谐?shù)、變量G,R,H,S的顯著性概率為0.00均小于0.01,因此該模型通過檢驗(yàn),具體模型為:

    4.1.3 工業(yè)堿溶液的濃度與顏色讀數(shù)的關(guān)系

    在工業(yè)堿溶液中,為了研究物質(zhì)濃度與顏色讀數(shù)的關(guān)系,我們首先求出同一個(gè)顏色在不同濃度下的顏色讀數(shù)極差,見表4。

    表4 工業(yè)堿的RGBHS值

    根據(jù)前面的理論,對其進(jìn)行多元線性回歸,發(fā)現(xiàn)其模型顯著性都大于0.05,模型不顯著。為了更好地分析物質(zhì)濃度與顏色讀數(shù)之間的關(guān)系,對工業(yè)堿溶液中的數(shù)據(jù)進(jìn)行變量之間的相關(guān)分析,具體結(jié)果見表5。

    表5 工業(yè)堿溶液中的數(shù)據(jù)進(jìn)行變量之間的相關(guān)分析

    從表5中可以看出工業(yè)堿溶液濃度與B,G,R,H,S的相關(guān)系數(shù)都低于0.8,顯著性均大于0.05,所以顯著性不明顯,該線性回歸無效。由表5可知顏色R,G,B都隨著工業(yè)堿溶液濃度的增加而減少,而色調(diào)H和飽和度S則隨著工業(yè)堿溶液濃度的增加而增加,因此工業(yè)堿溶液濃度與顏色R,G,B成負(fù)相關(guān),與H,S成正相關(guān),從極差可以看出飽和度S對工業(yè)堿溶液濃度的影響比較大,R的影響最小。

    4.1.4 硫酸鋁鉀溶液的濃度與顏色讀數(shù)的關(guān)系

    在硫酸鋁鉀溶液中,題目附件1給出了6組同一硫酸鋁鉀物質(zhì)濃度的數(shù)據(jù),為了減少實(shí)驗(yàn)的誤差,取6組實(shí)驗(yàn)數(shù)據(jù)的平均值,得出硫酸鋁鉀溶液在不同的濃度下的R,G,B,H,S值。根據(jù)前面的做法,對數(shù)據(jù)進(jìn)行多元線性回歸擬合,結(jié)果發(fā)現(xiàn)模型的顯著性不高。為了研究硫酸鋁鉀溶液的濃度與顏色讀數(shù)的關(guān)系,對數(shù)據(jù)進(jìn)行雙變量分析,通過SPSS22軟件計(jì)算,得出結(jié)果見表6。

    表6 硫酸鋁鉀溶液的濃度與顏色讀數(shù)的關(guān)系

    由表6可以看出硫酸鋁鉀溶液濃度與B,G,R,H,S的相關(guān)系數(shù)都低于0.7,顯著性均大于0.05,線性回歸無效,通過分析表6可以得出R隨著硫酸鋁鉀溶液濃度的增大而減少成負(fù)相關(guān),G隨著硫酸鋁鉀溶液濃度的變化值不明顯,成不相關(guān),顏色B,H,S隨著硫酸鋁鉀溶液的增大而增大,成正相關(guān),其中S與硫酸鋁鉀溶液濃度的極差最大,即影響最大。

    4.1.5 奶中尿素溶液的濃度與顏色讀數(shù)的關(guān)系

    在奶中尿素溶液中,題目附件1給出了兩組同一奶中尿素物質(zhì)濃度的數(shù)據(jù),為了減少實(shí)驗(yàn)的誤差,取兩組實(shí)驗(yàn)數(shù)據(jù)的平均值,得出奶中尿素溶液在不同的濃度下的R,G,B,H,S值。同理,對其進(jìn)行多項(xiàng)式擬合,通過數(shù)學(xué)軟件SPSS可以得出結(jié)果,該模型的顯著性為0.027小于0.05,回歸系數(shù)B,R,H,S的顯著性均小于0.05,所以線性回歸方程有效,B,R,H,S的線性回歸方程為:

    4.1.6 評價(jià)5組數(shù)據(jù)優(yōu)劣性的模型

    (1)從方差角度判別數(shù)據(jù)的優(yōu)劣。要判別5種物質(zhì)數(shù)據(jù)的優(yōu)劣,采用整理取平均值后的數(shù)據(jù)進(jìn)行分析,用Excel軟件求取各組物質(zhì)濃度與顏色讀數(shù)關(guān)系的方差,從而求出各種物質(zhì)濃度的方差平均值見表7。

    表7 各種物質(zhì)濃度方差

    根據(jù)統(tǒng)計(jì)學(xué)知識(shí),數(shù)據(jù)方差的值越小,即數(shù)據(jù)越穩(wěn)定。因此,這5組數(shù)據(jù)從優(yōu)到劣的順序依次為:組胺>奶中尿素>溴酸鉀>硫酸鋁鉀>工業(yè)堿。

    (2)從極差角度判斷數(shù)據(jù)的優(yōu)劣性。同上所示,采用整理取平均值后的數(shù)據(jù)進(jìn)行分析,用Excel軟件求各組物質(zhì)的濃度與顏色讀數(shù)關(guān)系的極差,從而求出各種物質(zhì)濃度的極差平均值,見表8:

    表8 各種物質(zhì)濃度的極差

    研究不同物質(zhì)溶液對不同顏色讀數(shù)數(shù)據(jù)的優(yōu)劣性,根據(jù)極差越小,證明數(shù)據(jù)越穩(wěn)定原則。這5組數(shù)據(jù)從優(yōu)到劣的順序依次為:組胺>硫酸鋁鉀>溴酸鉀>奶中尿素>工業(yè)堿。

    4.2 問題(2)模型的建立與求解

    4.2.1 建立顏色讀數(shù)和物質(zhì)濃度的數(shù)學(xué)模型

    首先對附錄2中的數(shù)據(jù)進(jìn)行補(bǔ)充完整,要研究顏色讀數(shù)與物質(zhì)濃度之間的關(guān)系,根據(jù)問題(1)的思路,建立多元線性回歸模型。運(yùn)用SPSS軟件對其數(shù)據(jù)進(jìn)行擬合,擬合方法選取Enter法,得出如下結(jié)果(見表9,10)。

    表9 變異數(shù)分析a

    表10 系數(shù)a

    根據(jù)以上結(jié)果,雖然該模型的顯著性小于0.01,然而,可以看出模型中顏色變量G,B,S的顯著性大于0.05。因此,運(yùn)用多元線性回歸模型Enter輸入法進(jìn)行擬合的線性回歸模型失效,因此該模型不好。為了能夠更好地研究顏色讀數(shù)和物質(zhì)濃度的數(shù)學(xué)模型,我們先對附件2中的數(shù)據(jù)進(jìn)行雙變量相關(guān)分析,運(yùn)用Pearson相關(guān)系數(shù),得出相應(yīng)的結(jié)果。從結(jié)果中可以看出,二氧化硫的濃度與R,G,H的相關(guān)系數(shù)較高,從而可判斷,二氧化硫的濃度與R,G,H有關(guān)。接著,通過SPSS22軟件進(jìn)行多元線性中的逐步回歸得出結(jié)果,見表11,12。

    表11 模型摘要

    表12 系數(shù)a

    可以得出,在顯著性相同的情況下,模型2的R值比模型1的R值大,因此選擇模型2更為合理。由此可得數(shù)學(xué)模型:

    其中y二氧化硫表示的是二氧化硫的濃度,xG表示顏色綠色的讀數(shù),xH表示色調(diào)的讀數(shù)。

    4.2.2 模型的誤差分析

    為了考慮所得出模型的誤差,我們引入了絕對誤差和相對誤差對其進(jìn)行分析,具體公式如下:絕對誤差公式為:

    其中:L為絕對誤差,N為真實(shí)值,M為預(yù)測值。

    相對誤差的公式為:

    結(jié)合附件2-1中的數(shù)據(jù),在模型(5)的作用下,我們得出模型(5)的預(yù)測值,接著,運(yùn)用公式(6)和公式(7),得到表13:

    從表13可以看出該模型數(shù)據(jù)的相對誤差0~1.8的范圍內(nèi),產(chǎn)生如此大的誤差,可能是由于數(shù)據(jù)的不穩(wěn)定造成的。

    表13 預(yù)測、絕對誤差與相對誤差的值

    4.3 問題(3)模型的建立與求解

    4.3.1 數(shù)據(jù)量對模型的影響

    為了探討數(shù)據(jù)量對模型的影響,取附件2中的數(shù)據(jù),做兩組實(shí)驗(yàn),實(shí)驗(yàn)組A:刪掉二氧化硫濃度為20ppm與50ppm的數(shù)據(jù)。實(shí)驗(yàn)組B:刪掉二氧化硫濃度為100ppm與150ppm的數(shù)據(jù)。通過數(shù)學(xué)軟件SPSS對實(shí)驗(yàn)組A中的數(shù)據(jù)進(jìn)行多元線性回歸擬合,運(yùn)用逐步法可得結(jié)果,見表14,15。

    表14 模型摘要

    表15 系數(shù)a

    從問題(2)的模型,可以知道在二氧化硫濃度下,未刪減數(shù)據(jù)時(shí)模型的R值為0.940,顯著性均為0。然而,實(shí)驗(yàn)A模型的顯著性也為0。因此,刪除數(shù)據(jù)的實(shí)驗(yàn)A與未刪除數(shù)據(jù)中問題二模型的顯著性是相同的。但是,通過做回歸發(fā)現(xiàn),實(shí)驗(yàn)組A中的R相比原數(shù)據(jù)的R值大,可以看出數(shù)據(jù)量對模型是有影響的。為了能夠減少實(shí)驗(yàn)中一些非主要因素的干擾,我們在做一組實(shí)驗(yàn)進(jìn)行分析。同理,運(yùn)用實(shí)驗(yàn)組B中的數(shù)據(jù)進(jìn)行多元線性回歸擬合,得出結(jié)果,見表16,17。

    表16 模型摘要b

    表17 系數(shù)a

    在模型顯著性一樣的前提條件下,可以看到實(shí)驗(yàn)組B模型的R值小于原數(shù)據(jù)模型的R值。因此,從兩組實(shí)驗(yàn)組綜合起來比較,可能實(shí)驗(yàn)B中刪除時(shí)刪掉了在濃度中主要的數(shù)據(jù)這種情況。通過用原數(shù)據(jù)與實(shí)驗(yàn)組數(shù)據(jù)之間預(yù)測值差距,進(jìn)一步判斷數(shù)據(jù)量對模型的影響,見表18。

    表18 原數(shù)據(jù)與實(shí)驗(yàn)組數(shù)據(jù)預(yù)測值比較

    50 80 80 80 100 100 100 150 150 150 150 85.57 184 75.34 260 75.34 260 70.29 437 95.53 553 100.58 377 90.48 730 132.11 332 130.47 484 137.16 155 132.11 332 80 80 80 100 100 100 150 150 150 150 85.91 058 85.91 058 80.86 630 106.08 770 111.13 198 101.04 342 143.75 282 142.00 402 148.79 710 143.75 282 50 80 80 80 55.76 353 59.63 625 59.63 625 59.63 625

    從表18中,可以明確地看出,未刪減過的數(shù)據(jù)與刪減過的數(shù)據(jù)對模型有著很明顯的影響,因此,可以得出數(shù)據(jù)量的多少與模型是有顯著影響的。

    4.3.2 顏色維度對模型的影響

    為了更好探討顏色維度對模型的影響,我們從增加顏色維度和減少顏色維度的兩個(gè)方向出發(fā),同樣,取Data2.xls中處理過的數(shù)據(jù)做兩組實(shí)驗(yàn)。

    根據(jù)參考文獻(xiàn)[3]可知,顏色維度I等于顏色維度R、G、B三者的平均值。

    實(shí)驗(yàn)組C:減少顏色維度,刪除掉H與S中的數(shù)據(jù);實(shí)驗(yàn)組D:增加顏色維度I的數(shù)據(jù)。通過數(shù)學(xué)軟件SPSS對實(shí)驗(yàn)組C中的數(shù)據(jù)進(jìn)行多元線性回歸擬合,運(yùn)用逐步法可得結(jié)果,見表19,20。

    表19 模型摘要

    表20 系數(shù)a

    從問題(2)的模型,可以知道在二氧化硫濃度下,原來模型的R值為0.940,顯著性均為0。然而,實(shí)驗(yàn)C模型的顯著性也為0。因此,刪除維度的實(shí)驗(yàn)C與問題二模型的顯著性是相同的。同時(shí),我們通過做回歸發(fā)現(xiàn),實(shí)驗(yàn)組C中的R相比原模型的R值相差不大。為了能夠更好地分析顏色維度對模型的影響,同理,繼續(xù)做實(shí)驗(yàn)組D中的數(shù)據(jù)進(jìn)行多元線性回歸擬合,得出結(jié)果。在模型顯著性一樣的前提條件下,實(shí)驗(yàn)D模型的R與原模型相差也不大。因此,我們繼續(xù)通過原模型數(shù)據(jù)與實(shí)驗(yàn)組數(shù)據(jù)之間的預(yù)測值進(jìn)行比較,見表21。

    表21 原數(shù)據(jù)與實(shí)驗(yàn)組數(shù)據(jù)預(yù)測值比較

    由表21可知,顏色維度的增減對于模型的影響還是較為明顯的。

    5 結(jié)語

    為了研究顏色讀數(shù)和物質(zhì)濃度之間的關(guān)系,本文結(jié)合2017年全國大學(xué)生數(shù)學(xué)建模競賽C題目,詳細(xì)地對不同物質(zhì)顏色讀數(shù)和濃度之間的關(guān)系進(jìn)行研究,且通過4組實(shí)驗(yàn)探索數(shù)據(jù)量增加與減少、顏色維度增加與減少對模型的影響。其次,本文考慮問題能夠從多個(gè)角度對問題進(jìn)行分析,給出分析結(jié)果,準(zhǔn)確性高;通過數(shù)學(xué)軟件Excel、SPSS對數(shù)據(jù)進(jìn)行分析和處理,操作簡便,適用性強(qiáng);且該模型具有較大的普遍性,便于應(yīng)用到類似的實(shí)際問題中,如:解決有關(guān)醫(yī)學(xué)研究、藥物測量與制作及化學(xué)等領(lǐng)域。

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    18禁黄网站禁片午夜丰满| 黑人操中国人逼视频| 老汉色av国产亚洲站长工具| 性色av乱码一区二区三区2| 欧美中文综合在线视频| 无遮挡黄片免费观看| av一本久久久久| 亚洲avbb在线观看| 精品一区二区三区四区五区乱码| 国产成人av教育| 18在线观看网站| 19禁男女啪啪无遮挡网站| av线在线观看网站| 99九九在线精品视频| 欧美大码av| 一级毛片精品| 亚洲精品国产av成人精品| 亚洲专区国产一区二区| 女人被躁到高潮嗷嗷叫费观| 亚洲精品一区蜜桃| 一区二区日韩欧美中文字幕| 91精品三级在线观看| 俄罗斯特黄特色一大片| 日韩视频一区二区在线观看| 性色av乱码一区二区三区2| a级毛片在线看网站| 男女边摸边吃奶| 亚洲美女黄色视频免费看| 最新在线观看一区二区三区| 国产亚洲欧美精品永久| 丰满迷人的少妇在线观看| 大片免费播放器 马上看| 在线十欧美十亚洲十日本专区| 99热全是精品| 老熟女久久久| 日韩视频在线欧美| 精品乱码久久久久久99久播| 国产精品.久久久| 纵有疾风起免费观看全集完整版| 大香蕉久久成人网| 日韩制服骚丝袜av| 国产日韩一区二区三区精品不卡| 精品久久久精品久久久| 妹子高潮喷水视频| 久久久国产欧美日韩av| 黄频高清免费视频| 亚洲avbb在线观看| 国产精品一二三区在线看| 国产欧美日韩一区二区三 | 99精品久久久久人妻精品| 中文字幕人妻熟女乱码| 一区二区三区精品91| 国产亚洲精品久久久久5区| 成年动漫av网址| 人妻 亚洲 视频| 亚洲va日本ⅴa欧美va伊人久久 | 操美女的视频在线观看| 亚洲伊人久久精品综合| 极品少妇高潮喷水抽搐| 9191精品国产免费久久| 啦啦啦中文免费视频观看日本| 黑人欧美特级aaaaaa片| 视频区图区小说| 欧美日韩福利视频一区二区| 纵有疾风起免费观看全集完整版| 啪啪无遮挡十八禁网站| 在线观看免费午夜福利视频| 精品久久蜜臀av无| 久久久久久久久免费视频了| 精品国产国语对白av| 手机成人av网站| 悠悠久久av| 少妇精品久久久久久久| 午夜福利,免费看| 黄频高清免费视频| 欧美日韩视频精品一区| 久久久久国产一级毛片高清牌| xxxhd国产人妻xxx| 欧美日韩国产mv在线观看视频| 亚洲成人手机| 视频区欧美日本亚洲| 捣出白浆h1v1| 五月开心婷婷网| 久久人妻熟女aⅴ| 日韩视频一区二区在线观看| 最近最新中文字幕大全免费视频| 久久久国产成人免费| 欧美日韩亚洲国产一区二区在线观看 | 最新在线观看一区二区三区| 9色porny在线观看| 老司机亚洲免费影院| 每晚都被弄得嗷嗷叫到高潮| 男女无遮挡免费网站观看| 岛国毛片在线播放| 黄频高清免费视频| 午夜成年电影在线免费观看| 亚洲精品在线美女| 国产男女超爽视频在线观看| 亚洲国产欧美在线一区| 日韩一区二区三区影片| 久久久国产精品麻豆| 菩萨蛮人人尽说江南好唐韦庄| 一本综合久久免费| 久久狼人影院| tocl精华| 亚洲自偷自拍图片 自拍| 一个人免费看片子| 美女大奶头黄色视频| 久久久久精品国产欧美久久久 | 男女之事视频高清在线观看| 大码成人一级视频| 亚洲成人免费电影在线观看| 亚洲av日韩精品久久久久久密| 老司机午夜十八禁免费视频| 欧美精品亚洲一区二区| 国产在线免费精品| av又黄又爽大尺度在线免费看| 一边摸一边抽搐一进一出视频| 国产日韩欧美在线精品| 中文欧美无线码| av国产精品久久久久影院| 嫩草影视91久久| 国产有黄有色有爽视频| 老司机午夜福利在线观看视频 | 欧美亚洲日本最大视频资源| 亚洲人成77777在线视频| 欧美黄色淫秽网站| 一区二区三区激情视频| 国产在线免费精品| 夜夜夜夜夜久久久久| 免费在线观看视频国产中文字幕亚洲 | 久9热在线精品视频| 精品少妇内射三级| 午夜免费观看性视频| 免费高清在线观看视频在线观看| 久久精品国产综合久久久| 欧美在线黄色| 成年人黄色毛片网站| 精品高清国产在线一区| 亚洲第一av免费看| 香蕉国产在线看| 国产成人a∨麻豆精品| 欧美精品高潮呻吟av久久| cao死你这个sao货| 国产av国产精品国产| 丁香六月天网| 两性夫妻黄色片| 亚洲av片天天在线观看| 日韩视频在线欧美| 精品国产乱码久久久久久小说| 丰满少妇做爰视频| 黄片小视频在线播放| 欧美人与性动交α欧美精品济南到| 99久久综合免费| 亚洲欧美色中文字幕在线| 久久免费观看电影| 淫妇啪啪啪对白视频 | 欧美人与性动交α欧美精品济南到| 亚洲色图综合在线观看| 国产一级毛片在线| 亚洲精品一区蜜桃| 2018国产大陆天天弄谢| 淫妇啪啪啪对白视频 | 国产亚洲av片在线观看秒播厂| 男女之事视频高清在线观看| 国产男女内射视频| 97人妻天天添夜夜摸| h视频一区二区三区| 欧美精品人与动牲交sv欧美| 国产高清国产精品国产三级| 国产精品国产三级国产专区5o| 在线观看免费视频网站a站| 久久免费观看电影| 日韩欧美国产一区二区入口| 精品一品国产午夜福利视频| 国产成人精品无人区| 爱豆传媒免费全集在线观看| 一本综合久久免费| 女人被躁到高潮嗷嗷叫费观| 51午夜福利影视在线观看| 99久久国产精品久久久| 精品少妇内射三级| 精品人妻一区二区三区麻豆| 久久国产精品大桥未久av| 可以免费在线观看a视频的电影网站| 下体分泌物呈黄色| 深夜精品福利| 亚洲欧美成人综合另类久久久| 国产亚洲精品一区二区www | 欧美日韩国产mv在线观看视频| 欧美97在线视频| 视频在线观看一区二区三区| 天堂8中文在线网| 成人国产一区最新在线观看| 天堂中文最新版在线下载| 国产高清videossex| 国产一级毛片在线| 99精国产麻豆久久婷婷| 在线av久久热| 日韩免费高清中文字幕av| 女人高潮潮喷娇喘18禁视频| www.av在线官网国产| 亚洲欧美成人综合另类久久久| 99久久99久久久精品蜜桃| 91成人精品电影| 久久人妻熟女aⅴ| 亚洲avbb在线观看| 一级a爱视频在线免费观看| 国产又色又爽无遮挡免| 亚洲av国产av综合av卡| 亚洲情色 制服丝袜| 亚洲精品在线美女| 久久久水蜜桃国产精品网| av福利片在线| 国产成人免费无遮挡视频| 少妇粗大呻吟视频| 91精品三级在线观看| 精品熟女少妇八av免费久了| 日日爽夜夜爽网站| 丝袜在线中文字幕| 欧美日本中文国产一区发布| 在线观看免费午夜福利视频| 亚洲成人手机| 成年av动漫网址| 成年人免费黄色播放视频| 免费观看a级毛片全部| 中国国产av一级| 天堂俺去俺来也www色官网| 国产欧美日韩一区二区三 | 成人三级做爰电影| 咕卡用的链子| 欧美日韩黄片免| 亚洲,欧美精品.| 免费人妻精品一区二区三区视频| 制服人妻中文乱码| 97精品久久久久久久久久精品| 少妇被粗大的猛进出69影院| 亚洲激情五月婷婷啪啪| 亚洲精品国产色婷婷电影| 亚洲熟女毛片儿| 久久久欧美国产精品| 欧美亚洲 丝袜 人妻 在线| 法律面前人人平等表现在哪些方面 | 午夜福利视频在线观看免费| 成人免费观看视频高清| 如日韩欧美国产精品一区二区三区| 久久综合国产亚洲精品| 狠狠精品人妻久久久久久综合| 久久女婷五月综合色啪小说| 亚洲免费av在线视频| a级毛片黄视频| 十八禁网站网址无遮挡| 国产精品秋霞免费鲁丝片| 欧美精品av麻豆av| 少妇 在线观看| 亚洲精品美女久久av网站| 午夜老司机福利片| 亚洲精品国产av成人精品| 精品一区二区三卡| 久久免费观看电影| 脱女人内裤的视频| 999精品在线视频| 久久久国产欧美日韩av| kizo精华| 亚洲五月婷婷丁香| 一级片'在线观看视频| 在线观看免费高清a一片| 久久久久久免费高清国产稀缺| 国产真人三级小视频在线观看| 久久久水蜜桃国产精品网| 18在线观看网站| 一级片'在线观看视频| 咕卡用的链子| 亚洲 国产 在线| 午夜福利乱码中文字幕| 肉色欧美久久久久久久蜜桃| 欧美日韩视频精品一区| 欧美另类亚洲清纯唯美| 91精品国产国语对白视频| a级毛片黄视频| 国产高清视频在线播放一区 | 操出白浆在线播放| 中国国产av一级| 欧美日韩精品网址| 极品少妇高潮喷水抽搐| 人人妻人人爽人人添夜夜欢视频| 亚洲av日韩在线播放| 久久精品久久久久久噜噜老黄| 高清av免费在线| 亚洲色图综合在线观看| 国产精品99久久99久久久不卡| 99久久人妻综合| 精品少妇黑人巨大在线播放| 人成视频在线观看免费观看| 少妇人妻久久综合中文| 两个人看的免费小视频| 亚洲一区中文字幕在线| 丰满饥渴人妻一区二区三| 三上悠亚av全集在线观看| 午夜影院在线不卡| 亚洲av电影在线进入| 久久香蕉激情| 国产日韩一区二区三区精品不卡| 一边摸一边抽搐一进一出视频| 涩涩av久久男人的天堂| 色94色欧美一区二区| 超色免费av| 久久久久视频综合| 午夜福利视频精品| 午夜两性在线视频| 国产成人系列免费观看| 十八禁网站网址无遮挡| 18在线观看网站| 国产成人啪精品午夜网站| 免费日韩欧美在线观看| 精品人妻一区二区三区麻豆| 亚洲久久久国产精品| 在线观看免费日韩欧美大片| av片东京热男人的天堂| 午夜免费鲁丝| 国产精品国产三级国产专区5o| 久热这里只有精品99| 亚洲国产精品一区二区三区在线| 亚洲av男天堂| 一级毛片精品| 久久久精品区二区三区| 日韩精品免费视频一区二区三区| 黄频高清免费视频| 久久久久视频综合| 大片免费播放器 马上看| 亚洲全国av大片| 亚洲国产日韩一区二区| 久久久国产成人免费| 亚洲情色 制服丝袜| 欧美成人午夜精品| 国产主播在线观看一区二区| 中文字幕另类日韩欧美亚洲嫩草| 正在播放国产对白刺激| 97人妻天天添夜夜摸| 久久av网站| 欧美日韩一级在线毛片| 亚洲欧美一区二区三区久久| 97在线人人人人妻| 中文字幕另类日韩欧美亚洲嫩草| 色综合欧美亚洲国产小说| 久久国产精品影院| 亚洲,欧美精品.| 国产精品成人在线| 天堂8中文在线网| 婷婷成人精品国产| 精品少妇黑人巨大在线播放| 亚洲人成电影免费在线| 亚洲色图综合在线观看| 纯流量卡能插随身wifi吗| 亚洲国产日韩一区二区| bbb黄色大片| 人妻人人澡人人爽人人| 久久精品aⅴ一区二区三区四区| 国产熟女午夜一区二区三区| 欧美精品啪啪一区二区三区 | 国产精品 国内视频| 精品高清国产在线一区| 国产精品国产av在线观看| 爱豆传媒免费全集在线观看| 99国产精品一区二区三区| 久久毛片免费看一区二区三区| 国产成人一区二区三区免费视频网站| a 毛片基地| 国产欧美亚洲国产| 欧美日韩亚洲国产一区二区在线观看 | 丝袜喷水一区| 美女脱内裤让男人舔精品视频| 大码成人一级视频| 亚洲成人免费电影在线观看| 久久天堂一区二区三区四区| 国产熟女午夜一区二区三区| 91精品伊人久久大香线蕉| 中文字幕人妻丝袜一区二区| 青青草视频在线视频观看| 免费在线观看完整版高清| 一级毛片电影观看| 欧美老熟妇乱子伦牲交| 精品国产乱子伦一区二区三区 | 国产免费一区二区三区四区乱码| 国产成人精品在线电影| 一区二区三区乱码不卡18| 国产有黄有色有爽视频| 免费日韩欧美在线观看| 中亚洲国语对白在线视频| 国产精品1区2区在线观看. | 99精品欧美一区二区三区四区| 久久热在线av| 久久久久精品人妻al黑| 91精品国产国语对白视频| 成人免费观看视频高清| 午夜福利影视在线免费观看| 亚洲七黄色美女视频| 欧美激情久久久久久爽电影 | www.精华液| 欧美97在线视频| 欧美日韩亚洲国产一区二区在线观看 | 在线观看免费午夜福利视频| 久久人人爽人人片av| 黄片大片在线免费观看| 亚洲国产欧美日韩在线播放| 国产成人一区二区三区免费视频网站| 丝袜喷水一区| www.精华液| 国产91精品成人一区二区三区 | 欧美亚洲日本最大视频资源| 18禁黄网站禁片午夜丰满| 欧美+亚洲+日韩+国产| 精品久久久久久电影网| 在线观看人妻少妇| 亚洲精品国产一区二区精华液| 久久狼人影院| 两个人看的免费小视频| 欧美激情极品国产一区二区三区| 王馨瑶露胸无遮挡在线观看| 天堂中文最新版在线下载| www日本在线高清视频| 精品福利永久在线观看| 少妇 在线观看| 一区二区三区精品91| 成人av一区二区三区在线看 | 欧美97在线视频| 国产激情久久老熟女| 久久精品aⅴ一区二区三区四区| av天堂久久9| 久久av网站| 久久青草综合色| 性高湖久久久久久久久免费观看| 午夜福利影视在线免费观看| 亚洲黑人精品在线| 国产亚洲一区二区精品| 女警被强在线播放| 波多野结衣av一区二区av| 一级片'在线观看视频| 80岁老熟妇乱子伦牲交| 淫妇啪啪啪对白视频 | 91老司机精品| 久久久国产精品麻豆| 高清欧美精品videossex| 纵有疾风起免费观看全集完整版| 丰满饥渴人妻一区二区三| 午夜精品久久久久久毛片777| 久久人人爽人人片av| 制服诱惑二区| 涩涩av久久男人的天堂| 国产精品免费大片| 亚洲精品美女久久久久99蜜臀| 丁香六月欧美| 久久精品成人免费网站| 成年人午夜在线观看视频| av网站免费在线观看视频| 久热爱精品视频在线9| 日韩大码丰满熟妇| 每晚都被弄得嗷嗷叫到高潮| 桃花免费在线播放| 精品人妻在线不人妻| 久久综合国产亚洲精品| 90打野战视频偷拍视频| 在线十欧美十亚洲十日本专区| 一级,二级,三级黄色视频| 在线观看www视频免费| 在线观看免费午夜福利视频| tube8黄色片| av一本久久久久| 精品人妻熟女毛片av久久网站| 99九九在线精品视频| 99国产精品一区二区蜜桃av | 亚洲国产精品成人久久小说| 91国产中文字幕| 91字幕亚洲| 久久久久久亚洲精品国产蜜桃av| 人妻 亚洲 视频| 人人妻人人添人人爽欧美一区卜| 国产精品自产拍在线观看55亚洲 | 丰满少妇做爰视频| 久久久久久久久久久久大奶| 国产精品久久久久久精品古装| 国产成人av教育| 久久狼人影院| 亚洲欧美日韩高清在线视频 | 首页视频小说图片口味搜索| 另类亚洲欧美激情| 中国国产av一级| 欧美黄色片欧美黄色片| 中国国产av一级| kizo精华| 青春草亚洲视频在线观看| 国产免费av片在线观看野外av| 欧美少妇被猛烈插入视频| 国产精品国产av在线观看| 免费在线观看日本一区| 免费观看人在逋| 午夜两性在线视频| 久久这里只有精品19| av在线老鸭窝| 久久人人爽av亚洲精品天堂| 99精国产麻豆久久婷婷| 久久人人爽人人片av| 99精品久久久久人妻精品| 欧美日韩亚洲高清精品| 国产成人av教育| 少妇精品久久久久久久| 久久综合国产亚洲精品| 少妇裸体淫交视频免费看高清 | 亚洲精品中文字幕在线视频| 日本av免费视频播放| 欧美 日韩 精品 国产| av又黄又爽大尺度在线免费看| 亚洲,欧美精品.| 在线观看舔阴道视频| 一级片免费观看大全| 国产欧美日韩一区二区三区在线| 大片电影免费在线观看免费| 成年人免费黄色播放视频| 无限看片的www在线观看| 久久精品成人免费网站| www.自偷自拍.com| 国产精品免费大片| 女人久久www免费人成看片| 手机成人av网站| 精品熟女少妇八av免费久了| 久久精品熟女亚洲av麻豆精品| 看免费av毛片| 国产免费福利视频在线观看| 黄色 视频免费看| 亚洲成国产人片在线观看| 亚洲av美国av| 欧美 亚洲 国产 日韩一| 久久九九热精品免费| 国产野战对白在线观看| 久久热在线av| 狠狠婷婷综合久久久久久88av| 99热全是精品| 女人精品久久久久毛片| a级片在线免费高清观看视频| 国产成人a∨麻豆精品| 中亚洲国语对白在线视频| 丝瓜视频免费看黄片| www.熟女人妻精品国产| 人人妻人人澡人人爽人人夜夜| 一级a爱视频在线免费观看| 久久久精品国产亚洲av高清涩受| 热99久久久久精品小说推荐| 天天躁日日躁夜夜躁夜夜| 最新的欧美精品一区二区| av在线app专区| 9色porny在线观看| 男女之事视频高清在线观看| 黄色 视频免费看| 欧美日韩亚洲国产一区二区在线观看 | 夫妻午夜视频| 老司机影院毛片| 亚洲精品成人av观看孕妇| e午夜精品久久久久久久| 久久中文看片网| 欧美人与性动交α欧美软件| 久久久国产成人免费| av欧美777| 丝袜美足系列| 桃红色精品国产亚洲av| 在线观看人妻少妇| 高潮久久久久久久久久久不卡| 性色av乱码一区二区三区2| 亚洲国产欧美日韩在线播放| 亚洲国产中文字幕在线视频| 色94色欧美一区二区| 一边摸一边做爽爽视频免费| 亚洲av男天堂| 欧美激情 高清一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 九色亚洲精品在线播放| 日本vs欧美在线观看视频| 午夜福利一区二区在线看| 国产成人a∨麻豆精品| 十八禁人妻一区二区| 无遮挡黄片免费观看| 久久青草综合色| 亚洲精华国产精华精| 涩涩av久久男人的天堂| 高潮久久久久久久久久久不卡| 天天操日日干夜夜撸| av线在线观看网站| 久久综合国产亚洲精品| 纵有疾风起免费观看全集完整版| 日韩欧美国产一区二区入口| 欧美久久黑人一区二区| 99国产精品一区二区蜜桃av | 国产无遮挡羞羞视频在线观看| av又黄又爽大尺度在线免费看| www.av在线官网国产| 秋霞在线观看毛片| 首页视频小说图片口味搜索| 黄色视频不卡| 99热全是精品| 色综合欧美亚洲国产小说| 色94色欧美一区二区| 久久久精品区二区三区| 美女扒开内裤让男人捅视频| 大陆偷拍与自拍| 婷婷色av中文字幕| 多毛熟女@视频| 欧美激情高清一区二区三区| 午夜视频精品福利| 国产高清国产精品国产三级| 亚洲七黄色美女视频| 亚洲av成人不卡在线观看播放网 | 国产一区二区三区在线臀色熟女 | 成人国语在线视频| 亚洲精品国产精品久久久不卡| 午夜两性在线视频| 欧美日韩亚洲国产一区二区在线观看 |