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

    差異驅(qū)動(dòng)型評(píng)價(jià)方法的穩(wěn)定性及差異凸顯能力比較

    2015-07-07 15:40:50李偉偉易平濤郭亞軍
    運(yùn)籌與管理 2015年1期
    關(guān)鍵詞:個(gè)數(shù)方差權(quán)重

    李偉偉, 易平濤, 郭亞軍

    (東北大學(xué) 工商管理學(xué)院,遼寧 沈陽(yáng) 110004)

    ?

    差異驅(qū)動(dòng)型評(píng)價(jià)方法的穩(wěn)定性及差異凸顯能力比較

    李偉偉, 易平濤, 郭亞軍

    (東北大學(xué) 工商管理學(xué)院,遼寧 沈陽(yáng) 110004)

    針對(duì)綜合評(píng)價(jià)中方法種類繁多且無統(tǒng)一比較標(biāo)準(zhǔn)的問題,選取了四種突出被評(píng)價(jià)對(duì)象之間差異的評(píng)價(jià)方法,采用隨機(jī)模擬的方式分別從評(píng)價(jià)方法的穩(wěn)定性及對(duì)差異的凸顯能力兩個(gè)方面進(jìn)行了比較分析。得出了四種方法穩(wěn)定性由高到低分別為均方差法、最大離差法、熵值法、拉開檔次法,且評(píng)價(jià)方法的穩(wěn)定性越高,則其對(duì)差異的凸顯能力反而越差的結(jié)論。該研究不僅驗(yàn)證了差異驅(qū)動(dòng)型評(píng)價(jià)方法的相關(guān)特性,為評(píng)價(jià)者關(guān)于評(píng)價(jià)方法的選取提供了參考意見,而且隨機(jī)模擬方法的應(yīng)用,可為類似的多評(píng)價(jià)方法的比較問題提供技術(shù)參考。

    綜合評(píng)價(jià);差異驅(qū)動(dòng);隨機(jī)模擬;穩(wěn)定性;差異凸顯

    0 引言

    綜合評(píng)價(jià)是指對(duì)被評(píng)價(jià)對(duì)象所進(jìn)行的客觀、公正、合理的全面評(píng)價(jià)[1],評(píng)價(jià)的宗旨是發(fā)現(xiàn)被評(píng)價(jià)對(duì)象之間的差異,以實(shí)現(xiàn)被評(píng)價(jià)對(duì)象之間的優(yōu)劣比較。發(fā)現(xiàn)被評(píng)價(jià)對(duì)象之間差異的過程,即是尋找原始評(píng)價(jià)信息之間的差異,并通過信息集結(jié)成的方式將這種差異表達(dá)出來的過程。在經(jīng)典的綜合評(píng)價(jià)中,連接原始評(píng)價(jià)信息及信息集結(jié)的一個(gè)重要環(huán)節(jié)是指標(biāo)的權(quán)重系數(shù),而常用的能夠體現(xiàn)評(píng)價(jià)信息之間差異(最終體現(xiàn)被評(píng)價(jià)對(duì)象之間的差異)的賦權(quán)方法有均方差法、熵值法、最大離差法和拉開檔次法,本文將其統(tǒng)稱為“差異驅(qū)動(dòng)型評(píng)價(jià)方法”。目前,關(guān)于各種差異驅(qū)動(dòng)型評(píng)價(jià)方法的單獨(dú)研究成果較多[2~10],卻很少涉及到四類差異驅(qū)動(dòng)型評(píng)價(jià)方法的比較問題,從而為評(píng)價(jià)者帶來了評(píng)價(jià)方法選取困難的問題。針對(duì)該問題,本文采用隨機(jī)模擬的方法,分別從方法的穩(wěn)定性及對(duì)被評(píng)價(jià)對(duì)象之間差異的凸顯能力兩個(gè)方面對(duì)四種差異驅(qū)動(dòng)型評(píng)價(jià)方法進(jìn)行了比較分析,并給出了相關(guān)結(jié)論。該研究將隨機(jī)模擬方法應(yīng)用于差異型評(píng)價(jià)方法特征測(cè)試或驗(yàn)證中,可為其他綜合評(píng)價(jià)方法的特征分析、性質(zhì)測(cè)度等問題的研究提供參考思路。

    1 差異驅(qū)動(dòng)型評(píng)價(jià)方法簡(jiǎn)介

    對(duì)于n個(gè)取定的被評(píng)價(jià)對(duì)象,關(guān)于m個(gè)評(píng)價(jià)指標(biāo)的觀測(cè)值矩陣為[xij]n×m。為簡(jiǎn)便起見,設(shè)xij(i=1,2,…,n;j=1,2,…,m)均為極大型指標(biāo),且已進(jìn)行了無量綱化處理。令各指標(biāo)對(duì)應(yīng)的權(quán)重系數(shù)向量為w=(w1,w2,…,wm)T,下面分別對(duì)四種差異型評(píng)價(jià)方法的基本原理進(jìn)行簡(jiǎn)單介紹。

    1.1 均方差法[1]

    用均方差法確定的權(quán)重系數(shù)為

    1.2 熵值法[1]

    運(yùn)用熵值法確定指標(biāo)權(quán)重系數(shù)的步驟如下:

    3)計(jì)算指標(biāo)xj的差異性系數(shù)

    當(dāng)xij差異越大,ej越小,第j項(xiàng)指標(biāo)對(duì)于被評(píng)價(jià)對(duì)象的比較作用越大。因此定義差異系數(shù)gj=1-ej,gj越大,越應(yīng)重視該項(xiàng)指標(biāo)的作用。

    4)確定權(quán)重系數(shù)

    1.3 最大離差法[2]

    運(yùn)用最大離差法求解指標(biāo)權(quán)重系數(shù)的基本思路為:權(quán)重系數(shù)的選擇應(yīng)使所有被評(píng)價(jià)對(duì)象之間的總體離差最大,即有以下優(yōu)化模型:

    求解優(yōu)化模型,并對(duì)求得的權(quán)重系數(shù)進(jìn)行歸一化處理,得最終的權(quán)重系數(shù)為

    1.4 拉開檔次法[1]

    取線性函數(shù)計(jì)算第i個(gè)被評(píng)價(jià)對(duì)象的評(píng)價(jià)值,得

    yi=w1xi1+w2xi2+…+wmxim

    (1)

    若記

    則式(1)可寫成y=Aw。

    確定權(quán)重系數(shù)向量w的準(zhǔn)則是使對(duì)n個(gè)被評(píng)價(jià)對(duì)象取值的分散程度或方差盡可能地大。

    而變量yi按n個(gè)被評(píng)價(jià)對(duì)象取值構(gòu)成樣本的方差為

    (2)

    (n-1)s2=wTATAw=wTHw

    (3)

    式中,H=ATA為實(shí)對(duì)稱矩陣。

    顯然,對(duì)w不加限制時(shí),式(3)可取任意大的值。這里限定wTw=1,求式(3)的最大值,也就是選擇w,使得

    (4)

    對(duì)于式(4),有如下結(jié)論:

    定理1 若取w為H的最大特征值所對(duì)應(yīng)的標(biāo)準(zhǔn)特征向量,式(4)取得最大值。

    定理1的證明詳見文獻(xiàn)[1]。

    2 穩(wěn)定性分析

    對(duì)于各種差異驅(qū)動(dòng)型評(píng)價(jià)方法的穩(wěn)定性,本文主要從極端樣本的加入對(duì)原始評(píng)價(jià)結(jié)果產(chǎn)生的干擾程度方面進(jìn)行分析,具體的分析思路為:在原始評(píng)價(jià)數(shù)據(jù)的基礎(chǔ)上增添極端樣本,即加入極端樣本的指標(biāo)數(shù)據(jù)相對(duì)于原始數(shù)據(jù)發(fā)生“突變”,形成“異常值”,分析加入異常值后原始被評(píng)價(jià)對(duì)象之間排序結(jié)論發(fā)生改變的程度,改變程度越小,則說明該方法的穩(wěn)定性越好。對(duì)于被評(píng)價(jià)對(duì)象之間排序結(jié)論發(fā)生改變的程度,這里用“Spearman’s等級(jí)相關(guān)系數(shù)”進(jìn)行測(cè)量,該系數(shù)的取值區(qū)間為[-1,1],取值越靠近1,說明排序結(jié)論發(fā)生改變的程度越小,相應(yīng)的評(píng)價(jià)方法的穩(wěn)定性越好;反之,則穩(wěn)定性越差。

    關(guān)于極端樣本的增添方式,本文主要考慮3種情形:

    因隨機(jī)模擬方法具有運(yùn)行充分,且可以獨(dú)立“組件”的方式融入到已有評(píng)價(jià)方法的求解過程中,所以為充分比較各差異驅(qū)動(dòng)型評(píng)價(jià)方法的穩(wěn)定性,本文選用隨機(jī)模擬的方法對(duì)各評(píng)價(jià)方法進(jìn)行大規(guī)模模擬實(shí)驗(yàn),原始實(shí)驗(yàn)數(shù)據(jù)按照 “均勻分布”方式隨機(jī)發(fā)生,可將具體的仿真步驟歸納如下:

    步驟1 依據(jù)指定行(指標(biāo)個(gè)數(shù))列(被評(píng)價(jià)對(duì)象個(gè)數(shù))隨機(jī)發(fā)生原始評(píng)價(jià)數(shù)據(jù)矩陣;

    步驟2 選用相應(yīng)的無量綱化方法對(duì)原始評(píng)價(jià)數(shù)據(jù)矩陣進(jìn)行預(yù)處理(針對(duì)拉開檔次法,需選用“標(biāo)準(zhǔn)化處理法”對(duì)原始數(shù)據(jù)進(jìn)行預(yù)處理);

    步驟3 運(yùn)用各種方法法求解指標(biāo)權(quán)重,并對(duì)評(píng)價(jià)信息進(jìn)行集結(jié),得到原始的標(biāo)準(zhǔn)排序;

    步驟4 按照上述方式增添極端樣本。

    步驟5 運(yùn)用各種方法重新計(jì)算權(quán)重系數(shù)及被評(píng)價(jià)對(duì)象的排序,用Spearman’s等級(jí)相關(guān)系數(shù)記錄評(píng)價(jià)結(jié)論發(fā)生變化的程度;

    步驟6 設(shè)置循環(huán)終止條件,即“各指標(biāo)在若干步(本文設(shè)置為200)后相關(guān)系數(shù)平均值增量的平均變動(dòng)程度小于等于某給定的精度(這里設(shè)為0.0001)”。

    依據(jù)上述方針步驟,按原始數(shù)據(jù)的最小值為基準(zhǔn)的方式增添極端樣本(令變異系數(shù)α=0.2),以10個(gè)被評(píng)價(jià)對(duì)象4個(gè)評(píng)價(jià)指標(biāo)為例,可得到不同情況下的仿真結(jié)果如圖1所示。

    圖1 四種差異驅(qū)動(dòng)型評(píng)價(jià)方法的穩(wěn)定性分析結(jié)果注:圖中橫軸為異常對(duì)象突變的步長(zhǎng)(t的取值),縱軸為Spearman’s等級(jí)相關(guān)系數(shù)的平均值。

    依據(jù)圖1可以看出,均方差法和最大離差法的仿真結(jié)果趨勢(shì)線類似,并與熵值法和拉開檔次法的趨勢(shì)線差異明顯。

    由圖1中的(a)和(c)可知:在均方差法和最大離差法中,“單指標(biāo)變異”極端樣本的加入方式對(duì)排序結(jié)果穩(wěn)定性的影響要高于“雙指標(biāo)變異”的情形,而“雙指標(biāo)變異”的情形要高于“整體指標(biāo)變異”的情形。均方差法和最大離差法是名副其實(shí)的突出“局部差異”評(píng)價(jià)方法。

    圖2 四種差異驅(qū)動(dòng)型評(píng)價(jià)方法的穩(wěn)定性比較注:圖中縱軸為所有突變步長(zhǎng)下Spearman’s等級(jí)相關(guān)系數(shù)的平均值

    由圖1中的(b)可知:在熵值法中,三種變異方式的趨勢(shì)線變換比較復(fù)雜,在前半段,總體來看,對(duì)排序結(jié)果穩(wěn)定性的影響程度由大到小依次為“整體指標(biāo)變異”方式、“雙指標(biāo)變異”方式、“單指標(biāo)變異”方式,該段表現(xiàn)出突出“整體差異”的特征;隨后,“單指標(biāo)變異”方式下數(shù)值持續(xù)降低,而“整體指標(biāo)變異”及“雙指標(biāo)變異”方式出現(xiàn)先后回轉(zhuǎn)的情形,這表明“熵值法”本質(zhì)上仍是一種突出“局部差異”的評(píng)價(jià)方法。

    由圖1中的(d)可知:在拉開檔次法中,“整體指標(biāo)變異”極端樣本的加入方式對(duì)排序結(jié)果穩(wěn)定性的影響要高于“雙指標(biāo)變異”的情形,而“雙指標(biāo)變異”的情形要高于“單指標(biāo)變異”的情形,拉開檔次法表現(xiàn)出明顯地突出“整體差異”的特征。

    綜合所有突變步長(zhǎng)(t=0,1,…,30)下的Spearman’s等級(jí)相關(guān)系數(shù)的平均值,得到四種差異驅(qū)動(dòng)型評(píng)價(jià)方法穩(wěn)定性比較如圖2所示。

    分析圖2可得到以下結(jié)論:

    (1)四種差異驅(qū)動(dòng)型評(píng)價(jià)方法的穩(wěn)定性由高到低依次為:均方差法、最大離差法、熵值法、拉開檔次法;

    (2)由各種變異情形下的穩(wěn)定性趨勢(shì)線再一次驗(yàn)證了:方差法及最大離差法是典型的突出“局部差異”的評(píng)價(jià)方法;熵值法本質(zhì)上為突出“局部差異”的方法,但帶有部分突出“整體差異”方法的特征;拉開檔次法是典型的突出“整體差異”的評(píng)價(jià)方法。

    3 差異凸顯能力分析

    關(guān)于四種評(píng)價(jià)方法的差異凸顯能力的分析,采用的主要思路是:用相同樣本所產(chǎn)生權(quán)重系數(shù)的標(biāo)準(zhǔn)差大小來衡量各種方法的差異凸顯能力。為充分展現(xiàn)各種方法的差異凸顯能力,本文分別考慮“指標(biāo)個(gè)數(shù)固定,被評(píng)價(jià)對(duì)象數(shù)目變動(dòng)”和“被評(píng)價(jià)對(duì)象個(gè)數(shù)固定,指標(biāo)數(shù)目變動(dòng)”兩種情形。在第一種情形中,令指標(biāo)個(gè)數(shù)為4個(gè),被評(píng)價(jià)對(duì)象個(gè)數(shù)分別為10,20,30,40,50,60,70;在第二種情形中,令被評(píng)價(jià)對(duì)象個(gè)數(shù)為15個(gè),指標(biāo)個(gè)數(shù)分別為4,6,8,10,12,14,16,18,20。原始數(shù)據(jù)采用隨機(jī)數(shù)發(fā)生器按“均勻分布”的方式在區(qū)間[0,1]內(nèi)隨機(jī)獲取,可將求解各種評(píng)價(jià)方法的差異凸顯能力的隨機(jī)模擬步驟歸納如下:

    步驟1 設(shè)置總仿真次數(shù)sum;

    步驟2 按均勻分布的方式在區(qū)間[0,1]內(nèi)隨機(jī)發(fā)生指定數(shù)目的被評(píng)價(jià)對(duì)象和指標(biāo)的原始評(píng)價(jià)數(shù)據(jù);

    步驟3 分別運(yùn)用四種評(píng)價(jià)方法對(duì)步驟2中發(fā)生的原始評(píng)價(jià)數(shù)據(jù)計(jì)算相應(yīng)的權(quán)重系數(shù);

    步驟4 計(jì)算并保存每次仿真運(yùn)行求得的權(quán)重系數(shù)的標(biāo)準(zhǔn)差,待仿真總次數(shù)等于sum時(shí),對(duì)每次仿真運(yùn)行得到的各種評(píng)價(jià)方法的權(quán)重系數(shù)的標(biāo)準(zhǔn)差求解平均值作為各種方法差異凸顯能力的測(cè)度值。

    在總仿真次數(shù)為sum=20000次的情形下,得各種評(píng)價(jià)方法的差異凸顯能力比較結(jié)果如圖3所示。

    圖3 各種評(píng)價(jià)方法的差異凸顯能力比較圖注:圖中縱軸為2萬次仿真結(jié)果的權(quán)重系數(shù)標(biāo)準(zhǔn)差的平均值;橫軸為被評(píng)價(jià)對(duì)象的數(shù)目(圖(a))或指標(biāo)數(shù)目(圖(b))

    分析圖3可以看出:

    (1)無論是指標(biāo)數(shù)目固定被評(píng)價(jià)對(duì)象個(gè)數(shù)變動(dòng),還是被評(píng)價(jià)對(duì)象數(shù)目固定指標(biāo)個(gè)數(shù)變動(dòng),四種方法對(duì)差異的凸顯能力由高到低分別為拉開檔次法、熵值法、最大離差法、均方差法;

    (2)當(dāng)指標(biāo)數(shù)目固定被評(píng)價(jià)對(duì)象個(gè)數(shù)變動(dòng)時(shí),隨著被評(píng)價(jià)對(duì)象個(gè)數(shù)的增加,拉開檔次法的差異凸顯能力略有增加,熵值法、最大離差法和均方差法的差異凸顯能力則有所降低;

    (3)當(dāng)被評(píng)價(jià)對(duì)象數(shù)目固定指標(biāo)個(gè)數(shù)變動(dòng)時(shí),隨著指標(biāo)個(gè)數(shù)的增加,四種評(píng)價(jià)方法的差異凸顯能力均有所下降;

    (4)結(jié)合圖3中(a)和(b)可以看出,指標(biāo)個(gè)數(shù)變動(dòng)對(duì)方法差異凸顯能力的影響要大于被評(píng)價(jià)對(duì)象個(gè)數(shù)變動(dòng)對(duì)方法差異凸顯能力的影響;

    (5)與四種評(píng)價(jià)方法的穩(wěn)定性比較結(jié)論相比,可知評(píng)價(jià)方法的穩(wěn)定性越低,對(duì)差異的凸顯能力則越強(qiáng);

    (6)最大離差法的差異凸顯能力測(cè)度值恒“略大于”均方差法,從總體上看,兩個(gè)方法的差異凸顯能力幾乎完全一致,這說明兩種方法具有內(nèi)在很強(qiáng)的關(guān)聯(lián)性。

    4 結(jié)束語(yǔ)

    本文運(yùn)用隨機(jī)模擬仿真的方法對(duì)四種差異驅(qū)動(dòng)型評(píng)價(jià)方法的穩(wěn)定性及差異凸顯能力進(jìn)行了比較分析,具有以下兩方面的研究意義:(1)驗(yàn)證了各種方法的特征,即均方差法和最大離差法是典型的突出“局部差異”的評(píng)價(jià)方法,拉開檔次法是典型的突出“整體差異”的評(píng)價(jià)方法,而熵值法雖是突出“局部差異”的評(píng)價(jià)方法,卻帶有部分突出“整體差異”的特征,該研究結(jié)論可為評(píng)價(jià)者關(guān)于評(píng)價(jià)方法的選取提供參考意見;(2)將隨機(jī)模擬方法應(yīng)用于方法性質(zhì)的測(cè)驗(yàn)中,對(duì)目前綜合評(píng)價(jià)領(lǐng)域“方法種類繁多,但缺乏統(tǒng)一的判斷標(biāo)準(zhǔn),即多方法結(jié)論非一致性”問題的解決具有啟發(fā)作用,尤其適合于同類別評(píng)價(jià)方法的性質(zhì)測(cè)驗(yàn)、特征比較等問題的分析。本文僅對(duì)四種差異驅(qū)動(dòng)型評(píng)價(jià)方法的穩(wěn)定性及差異凸顯能力進(jìn)行了數(shù)值驗(yàn)證,進(jìn)一步將從方法穩(wěn)定性的提升、兼顧穩(wěn)定性及差異凸顯能力的評(píng)價(jià)模式的構(gòu)建等方面展開深入研究。

    [1] 郭亞軍.綜合評(píng)價(jià)理論、方法及應(yīng)用[M].北京:科學(xué)出版社,2007.

    [2] 王應(yīng)明.運(yùn)用離差最大化方法進(jìn)行多指標(biāo)決策與排序[J].軟科學(xué)研究,1998,3:36-38,65.

    [3] 郭亞軍,馬鳳妹,董慶興.無量綱化方法對(duì)拉開檔次法的影響分析[J].管理科學(xué)學(xué)報(bào),2011,14(5):19-28.

    [4] 王道平,王煦.基于AHP/熵值法的鋼鐵企業(yè)綠色供應(yīng)商選擇指標(biāo)權(quán)重研究[J].軟科學(xué),2010,24(8):117-122.

    [5] 戚宇,郭亞軍,郭英民.主觀信息嵌入式的拉開檔次法及其應(yīng)用[J].東北大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,32(7):1057-1060,1064.

    [6] Zhihong Zou, Yi Yun, Jingnan Sun. Entropy method for determination of weight of evaluating indicators in fuzzy synthetic evaluation for water quality assessment[J]. Journal of Environmental Sciences, 2006, 18(5): 1020-1-23.

    [7] 朱方霞,陳華友.確定區(qū)間數(shù)決策矩陣屬性權(quán)重的方法——熵值法[J].安徽大學(xué)學(xué)報(bào)(自然科學(xué)版),2005,30(5):4- 6.

    [8] Dunwen Liu, Lei Yu, Bo Li. Fuzzy-entropy theory comprehensive evaluation method and its application in building construction safety[J]. Procedia Engineering, 2012, 43: 137-142.

    [9] Chuansheng Xie, Dapeng Dong, Shengping Hua, Xin Xu, Yingjie Chen. Safety evaluation of smart grid based on AHH-entropy method[J]. Systems Engineering Procedia, 2012, 4: 203-209.

    [10] Yun Ye. Fuzzy cross entropy of interval-valued intuitionstic fuzzy sets and its optimal decision-making method based on the weights of alternatives[J]. Expert Systems with Applications, 2011, 38(5): 6179- 6183.

    Comparison of Difference Driving Evaluation Methods in Stability and Difference Highlighting

    LI Wei-wei, YI Ping-tao, GUO Ya-jun

    (SchoolofBusinessAdministration,NortheasternUniversity,Shenyang110004,China)

    To the problem of multiple methods without unified comparison standards in comprehensive evaluation, we choose four evaluation methods highlighting the difference among objects to analyze their stability and difference highlighting capability using the approach of stochastic simulation. The conclusions are drawn that the stability of four methods from high to low is variance, maximum deviation, entropy, deviation maximization method, and the more stable of a method, the less capability of highlighting difference. This research not only can test and verify the features of four methods which could provide reference for evaluators of methods selecting, but also can provide technology support, that is stochastic simulation method, for other evaluation methods comparison.

    comprehensive evaluation; difference driving; stochastic simulation; stability; difference highlighting

    2013- 09-11

    國(guó)家自然科學(xué)基金資助項(xiàng)目(71071031);教育部基本科研業(yè)務(wù)費(fèi)項(xiàng)目(N120606002,N130406004)

    李偉偉(1986-),女,山東煙臺(tái)人,博士后,研究方向:綜合評(píng)價(jià);易平濤(1981-),男,湖南永州人,博士,副教授,碩士生導(dǎo)師,研究方向:系統(tǒng)評(píng)價(jià)與信息融合。

    C934

    A

    1007-3221(2015)01- 0216- 06

    猜你喜歡
    個(gè)數(shù)方差權(quán)重
    方差怎么算
    怎樣數(shù)出小正方體的個(gè)數(shù)
    概率與統(tǒng)計(jì)(2)——離散型隨機(jī)變量的期望與方差
    權(quán)重常思“浮名輕”
    等腰三角形個(gè)數(shù)探索
    怎樣數(shù)出小木塊的個(gè)數(shù)
    計(jì)算方差用哪個(gè)公式
    怎樣數(shù)出小正方體的個(gè)數(shù)
    為黨督政勤履職 代民行權(quán)重?fù)?dān)當(dāng)
    方差生活秀
    av在线天堂中文字幕| 亚洲人与动物交配视频| 亚洲成人久久爱视频| 水蜜桃什么品种好| 亚洲精品456在线播放app| 尤物成人国产欧美一区二区三区| 老司机影院成人| 亚洲av不卡在线观看| 伊人久久国产一区二区| 亚洲欧美成人精品一区二区| 三级男女做爰猛烈吃奶摸视频| 伦理电影大哥的女人| 一个人看的www免费观看视频| 夜夜看夜夜爽夜夜摸| 国产乱人视频| 综合色av麻豆| 麻豆乱淫一区二区| 中文字幕久久专区| 在线天堂最新版资源| 99热国产这里只有精品6| 国产高清国产精品国产三级 | 一级av片app| 欧美变态另类bdsm刘玥| 亚洲精品亚洲一区二区| 综合色av麻豆| 天堂俺去俺来也www色官网| 秋霞伦理黄片| 久久精品人妻少妇| 国产男女内射视频| 国语对白做爰xxxⅹ性视频网站| 在线 av 中文字幕| 日韩一区二区视频免费看| 国产欧美亚洲国产| 天天躁夜夜躁狠狠久久av| xxx大片免费视频| 能在线免费看毛片的网站| 久久人人爽人人片av| 青春草视频在线免费观看| 内射极品少妇av片p| 网址你懂的国产日韩在线| 亚洲第一区二区三区不卡| 国产精品秋霞免费鲁丝片| 午夜激情久久久久久久| av国产精品久久久久影院| 美女xxoo啪啪120秒动态图| 久久韩国三级中文字幕| 天堂俺去俺来也www色官网| 少妇的逼水好多| av天堂中文字幕网| 日韩中字成人| 亚洲精品影视一区二区三区av| 久久久久久久国产电影| 国产又色又爽无遮挡免| 亚洲精品日韩在线中文字幕| 日韩三级伦理在线观看| 亚洲国产精品成人综合色| 国产69精品久久久久777片| av在线老鸭窝| 搞女人的毛片| 亚洲精品久久久久久婷婷小说| 国国产精品蜜臀av免费| 在现免费观看毛片| 97在线人人人人妻| 国产精品一区二区三区四区免费观看| 国产真实伦视频高清在线观看| 亚洲精品乱久久久久久| av在线播放精品| 成人鲁丝片一二三区免费| 欧美日本视频| 久久久久精品性色| a级一级毛片免费在线观看| 三级国产精品欧美在线观看| 国产精品久久久久久精品电影| 国产成人免费无遮挡视频| 天天躁日日操中文字幕| 特大巨黑吊av在线直播| 国产日韩欧美亚洲二区| 国产精品99久久久久久久久| 人体艺术视频欧美日本| 日日啪夜夜撸| 免费电影在线观看免费观看| 天天躁日日操中文字幕| 免费看日本二区| 一个人看视频在线观看www免费| 婷婷色麻豆天堂久久| 日韩大片免费观看网站| 性插视频无遮挡在线免费观看| 国产黄色视频一区二区在线观看| 国产色爽女视频免费观看| 亚洲丝袜综合中文字幕| 一区二区三区精品91| 亚洲av男天堂| 亚洲,一卡二卡三卡| 中文字幕亚洲精品专区| 18+在线观看网站| 欧美日韩一区二区视频在线观看视频在线 | 国产黄色视频一区二区在线观看| 精品视频人人做人人爽| 国产成人精品久久久久久| 各种免费的搞黄视频| 91在线精品国自产拍蜜月| 国产 精品1| 久热这里只有精品99| 国产永久视频网站| 97在线视频观看| 永久网站在线| 日韩av免费高清视频| av女优亚洲男人天堂| 国产精品久久久久久久久免| 久久久久久久久久人人人人人人| 在线观看一区二区三区| 国产av国产精品国产| 中文字幕免费在线视频6| 日韩成人伦理影院| 亚洲av国产av综合av卡| 大片电影免费在线观看免费| 秋霞在线观看毛片| 久久久久久九九精品二区国产| 一区二区av电影网| 一二三四中文在线观看免费高清| 男人狂女人下面高潮的视频| 国产精品99久久久久久久久| 91午夜精品亚洲一区二区三区| 日日撸夜夜添| 亚洲,一卡二卡三卡| 亚洲不卡免费看| 日日摸夜夜添夜夜添av毛片| 亚洲av国产av综合av卡| 亚洲综合色惰| 国产精品嫩草影院av在线观看| 成人亚洲精品av一区二区| 欧美+日韩+精品| 美女主播在线视频| 久久久精品欧美日韩精品| 亚洲精品日本国产第一区| 精品亚洲乱码少妇综合久久| av网站免费在线观看视频| 亚洲一级一片aⅴ在线观看| 国产免费视频播放在线视频| 中文字幕免费在线视频6| 精品国产露脸久久av麻豆| 国内揄拍国产精品人妻在线| 国产精品久久久久久av不卡| 欧美 日韩 精品 国产| 国产精品一区二区在线观看99| 国产精品福利在线免费观看| 亚洲精品久久午夜乱码| 丝袜喷水一区| 三级国产精品欧美在线观看| 少妇的逼好多水| 久久人人爽人人片av| 午夜免费观看性视频| 免费观看在线日韩| 亚洲一级一片aⅴ在线观看| 精品人妻熟女av久视频| kizo精华| 国产极品天堂在线| 精品久久久久久久久av| 亚洲精品国产av蜜桃| 亚洲性久久影院| 午夜日本视频在线| 国产精品爽爽va在线观看网站| 国产精品蜜桃在线观看| 九草在线视频观看| 久久99热这里只有精品18| 在现免费观看毛片| 另类亚洲欧美激情| 2021天堂中文幕一二区在线观| 免费观看性生交大片5| 欧美3d第一页| 日韩,欧美,国产一区二区三区| 久久久午夜欧美精品| 国产欧美另类精品又又久久亚洲欧美| 久久国内精品自在自线图片| 联通29元200g的流量卡| 91久久精品电影网| 国产高清有码在线观看视频| 狠狠精品人妻久久久久久综合| 久久鲁丝午夜福利片| 麻豆成人av视频| 欧美三级亚洲精品| 亚州av有码| 亚洲欧美清纯卡通| 日韩视频在线欧美| 男女下面进入的视频免费午夜| 亚洲精品久久午夜乱码| 欧美极品一区二区三区四区| 国产成人一区二区在线| 日韩大片免费观看网站| 最近手机中文字幕大全| 国产色爽女视频免费观看| 久久女婷五月综合色啪小说 | 亚洲婷婷狠狠爱综合网| 国产色婷婷99| 91久久精品电影网| 久久6这里有精品| 91精品国产九色| 亚洲,欧美,日韩| 亚洲精品,欧美精品| 毛片一级片免费看久久久久| 丝瓜视频免费看黄片| 三级经典国产精品| 99久久精品热视频| 狂野欧美激情性bbbbbb| 欧美老熟妇乱子伦牲交| 只有这里有精品99| 一个人观看的视频www高清免费观看| 欧美极品一区二区三区四区| 白带黄色成豆腐渣| 日韩大片免费观看网站| 欧美成人精品欧美一级黄| 伊人久久精品亚洲午夜| 美女国产视频在线观看| 青春草国产在线视频| 日韩成人伦理影院| 亚洲av成人精品一二三区| 免费观看无遮挡的男女| 大又大粗又爽又黄少妇毛片口| 午夜福利在线观看免费完整高清在| 大片电影免费在线观看免费| 成年人午夜在线观看视频| 看黄色毛片网站| 2022亚洲国产成人精品| 精品少妇久久久久久888优播| 99久久人妻综合| av福利片在线观看| 99热网站在线观看| 在线免费十八禁| 麻豆成人av视频| 亚洲精品日韩在线中文字幕| 女人十人毛片免费观看3o分钟| 欧美丝袜亚洲另类| av福利片在线观看| av又黄又爽大尺度在线免费看| 国产精品嫩草影院av在线观看| 成人鲁丝片一二三区免费| 国产av码专区亚洲av| 久久久久久久久大av| av国产久精品久网站免费入址| 91精品一卡2卡3卡4卡| eeuss影院久久| 热99国产精品久久久久久7| 成人一区二区视频在线观看| 久久久亚洲精品成人影院| 嫩草影院入口| 男女无遮挡免费网站观看| 中文天堂在线官网| 欧美激情在线99| 久久精品国产鲁丝片午夜精品| 肉色欧美久久久久久久蜜桃 | 99热这里只有是精品50| 日本-黄色视频高清免费观看| 99久久中文字幕三级久久日本| 国产69精品久久久久777片| 日韩,欧美,国产一区二区三区| 成人国产麻豆网| 久久99热这里只频精品6学生| 一级黄片播放器| 五月开心婷婷网| 亚洲国产av新网站| 国产片特级美女逼逼视频| 国产精品不卡视频一区二区| av播播在线观看一区| 免费观看无遮挡的男女| 中文资源天堂在线| 亚洲av免费在线观看| 日本一二三区视频观看| 在线精品无人区一区二区三 | 色吧在线观看| 在线 av 中文字幕| 国产探花极品一区二区| 97在线视频观看| 亚洲精品456在线播放app| 国产白丝娇喘喷水9色精品| 亚洲欧洲国产日韩| 麻豆乱淫一区二区| 三级经典国产精品| 色哟哟·www| 亚洲久久久久久中文字幕| 亚洲国产高清在线一区二区三| 亚洲国产欧美在线一区| 女人被狂操c到高潮| 日韩欧美 国产精品| 久久久精品欧美日韩精品| 特大巨黑吊av在线直播| 永久免费av网站大全| 精品国产三级普通话版| 亚洲av免费高清在线观看| 看十八女毛片水多多多| 精品久久国产蜜桃| 中文字幕制服av| 国产老妇女一区| 精品少妇久久久久久888优播| 日日啪夜夜爽| 国产乱人视频| 精品少妇久久久久久888优播| 亚洲成人久久爱视频| 国产探花极品一区二区| 中文精品一卡2卡3卡4更新| 大香蕉97超碰在线| 尾随美女入室| 国产成年人精品一区二区| 亚洲av不卡在线观看| 精品熟女少妇av免费看| 一级毛片电影观看| 成人亚洲精品av一区二区| 色综合色国产| 亚洲av一区综合| videos熟女内射| 精品人妻偷拍中文字幕| 69av精品久久久久久| 中国美白少妇内射xxxbb| 九九久久精品国产亚洲av麻豆| 国内精品宾馆在线| 九九在线视频观看精品| 欧美日韩视频高清一区二区三区二| 久久久久久久久久久丰满| 91aial.com中文字幕在线观看| 精品久久久久久久末码| 黄色配什么色好看| 久久97久久精品| 久久久久久久国产电影| av免费观看日本| 99久国产av精品国产电影| 成人午夜精彩视频在线观看| 最近最新中文字幕免费大全7| 51国产日韩欧美| 五月开心婷婷网| 天美传媒精品一区二区| 午夜精品一区二区三区免费看| av在线观看视频网站免费| 国产av国产精品国产| 人人妻人人爽人人添夜夜欢视频 | 韩国av在线不卡| 亚洲欧美一区二区三区黑人 | 成人综合一区亚洲| 中文字幕av成人在线电影| 国产亚洲av嫩草精品影院| 人人妻人人看人人澡| 哪个播放器可以免费观看大片| 少妇 在线观看| 亚洲经典国产精华液单| 综合色av麻豆| 高清在线视频一区二区三区| 成人国产麻豆网| 日韩 亚洲 欧美在线| 亚洲国产日韩一区二区| 亚洲av一区综合| 中文字幕制服av| 寂寞人妻少妇视频99o| 久久综合国产亚洲精品| 国产大屁股一区二区在线视频| 男人舔奶头视频| 亚洲av一区综合| 午夜免费男女啪啪视频观看| 99久久精品热视频| 又爽又黄无遮挡网站| 久久99精品国语久久久| 午夜福利视频1000在线观看| 人妻一区二区av| 噜噜噜噜噜久久久久久91| 亚洲精品乱码久久久v下载方式| 午夜福利在线在线| 国产69精品久久久久777片| 国产成人精品久久久久久| av国产精品久久久久影院| 国产成人a区在线观看| 色吧在线观看| 亚洲精品国产av蜜桃| 在线观看av片永久免费下载| 亚洲精品第二区| 免费av不卡在线播放| 一级毛片电影观看| 欧美日本视频| 色哟哟·www| 亚洲国产精品国产精品| 免费电影在线观看免费观看| 成人亚洲精品av一区二区| 久久久亚洲精品成人影院| 日韩 亚洲 欧美在线| 美女脱内裤让男人舔精品视频| 免费大片18禁| 日韩国内少妇激情av| 日本猛色少妇xxxxx猛交久久| 18禁动态无遮挡网站| eeuss影院久久| 在线免费观看不下载黄p国产| 久久女婷五月综合色啪小说 | 国产精品国产三级国产专区5o| 99精国产麻豆久久婷婷| 亚洲经典国产精华液单| 性色avwww在线观看| 国产av码专区亚洲av| 免费观看无遮挡的男女| 成人亚洲欧美一区二区av| 嫩草影院精品99| 国产日韩欧美在线精品| 欧美另类一区| 在线观看一区二区三区激情| 久久精品久久久久久噜噜老黄| 特大巨黑吊av在线直播| 2021天堂中文幕一二区在线观| 亚洲精品国产av成人精品| 永久网站在线| 能在线免费看毛片的网站| 联通29元200g的流量卡| 一级毛片电影观看| 99热这里只有是精品50| 搡女人真爽免费视频火全软件| 国产v大片淫在线免费观看| 欧美精品人与动牲交sv欧美| 精品少妇久久久久久888优播| 精品久久久久久久人妻蜜臀av| 亚洲精品乱久久久久久| 国产在线一区二区三区精| 我的女老师完整版在线观看| 中文字幕久久专区| 国产精品国产三级国产av玫瑰| 国产精品99久久99久久久不卡 | 国产精品精品国产色婷婷| 男女边吃奶边做爰视频| 美女cb高潮喷水在线观看| 一级爰片在线观看| 亚洲第一区二区三区不卡| 久久久久网色| xxx大片免费视频| 男女边吃奶边做爰视频| 国产在视频线精品| 丝瓜视频免费看黄片| 夜夜看夜夜爽夜夜摸| 少妇人妻一区二区三区视频| 久久久精品94久久精品| 在线 av 中文字幕| 成人欧美大片| 中文字幕av成人在线电影| 在线观看三级黄色| 午夜爱爱视频在线播放| 亚洲人与动物交配视频| 少妇的逼水好多| 国产黄色免费在线视频| 天美传媒精品一区二区| 欧美xxxx黑人xx丫x性爽| 又爽又黄无遮挡网站| 伦精品一区二区三区| 99久久精品热视频| 国产一区二区亚洲精品在线观看| 性色av一级| 国产乱人视频| 最近中文字幕高清免费大全6| 精品亚洲乱码少妇综合久久| 成年女人看的毛片在线观看| 精品久久国产蜜桃| 少妇熟女欧美另类| 免费看不卡的av| 亚洲欧美精品专区久久| 国国产精品蜜臀av免费| 一级a做视频免费观看| 亚州av有码| 国产精品熟女久久久久浪| 国产大屁股一区二区在线视频| tube8黄色片| 亚洲国产高清在线一区二区三| 最近最新中文字幕免费大全7| 视频中文字幕在线观看| 国产片特级美女逼逼视频| 80岁老熟妇乱子伦牲交| www.色视频.com| 国产在视频线精品| 在线精品无人区一区二区三 | 亚洲精品国产av蜜桃| 成年女人在线观看亚洲视频 | 国产永久视频网站| 免费观看a级毛片全部| 日日摸夜夜添夜夜爱| 乱码一卡2卡4卡精品| 边亲边吃奶的免费视频| 亚洲精品色激情综合| 丰满少妇做爰视频| 色视频www国产| 亚洲无线观看免费| 国产精品秋霞免费鲁丝片| 91久久精品国产一区二区成人| 日韩一区二区三区影片| 成人亚洲精品一区在线观看 | 夜夜看夜夜爽夜夜摸| 欧美国产精品一级二级三级 | 国产91av在线免费观看| 十八禁网站网址无遮挡 | 听说在线观看完整版免费高清| 狂野欧美激情性bbbbbb| 国产毛片a区久久久久| 免费av毛片视频| 亚洲美女搞黄在线观看| 亚洲av欧美aⅴ国产| 亚洲va在线va天堂va国产| 又爽又黄a免费视频| 亚洲精品乱码久久久久久按摩| 国产亚洲午夜精品一区二区久久 | www.av在线官网国产| 日韩欧美精品v在线| 久久国产乱子免费精品| 少妇裸体淫交视频免费看高清| 综合色av麻豆| 91精品国产九色| www.av在线官网国产| 久久久久精品久久久久真实原创| 欧美zozozo另类| 在线免费观看不下载黄p国产| 久久久欧美国产精品| 欧美亚洲 丝袜 人妻 在线| 国产男人的电影天堂91| av女优亚洲男人天堂| 内地一区二区视频在线| 亚洲欧美日韩东京热| 日本-黄色视频高清免费观看| 精品午夜福利在线看| 一级av片app| 黄色一级大片看看| 亚洲激情五月婷婷啪啪| 国产高清国产精品国产三级 | 国产精品国产三级专区第一集| 国国产精品蜜臀av免费| 人妻系列 视频| 亚洲性久久影院| 久久99蜜桃精品久久| 天美传媒精品一区二区| 欧美极品一区二区三区四区| .国产精品久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 大香蕉久久网| 黄片wwwwww| 嫩草影院精品99| 大香蕉97超碰在线| 国产色婷婷99| 日产精品乱码卡一卡2卡三| 国产亚洲av嫩草精品影院| 国产淫语在线视频| 欧美亚洲 丝袜 人妻 在线| 午夜亚洲福利在线播放| 99热这里只有是精品在线观看| 成人高潮视频无遮挡免费网站| 久久精品国产自在天天线| 欧美xxxx黑人xx丫x性爽| 国产日韩欧美在线精品| 在线天堂最新版资源| 一级a做视频免费观看| 国产视频内射| 九草在线视频观看| 久久久久精品久久久久真实原创| 午夜免费鲁丝| 好男人在线观看高清免费视频| 老师上课跳d突然被开到最大视频| 国产亚洲av片在线观看秒播厂| 精品国产露脸久久av麻豆| 一边亲一边摸免费视频| 国产亚洲一区二区精品| 伊人久久精品亚洲午夜| 99视频精品全部免费 在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品一区二区三卡| 久久久久久久久久人人人人人人| 国产男女内射视频| 一级毛片我不卡| 久久精品国产鲁丝片午夜精品| 国产精品一区二区性色av| 噜噜噜噜噜久久久久久91| 狂野欧美激情性bbbbbb| 日本熟妇午夜| 99热网站在线观看| 狂野欧美白嫩少妇大欣赏| 久久久久国产网址| 国产精品人妻久久久影院| 欧美xxxx性猛交bbbb| 国产欧美日韩一区二区三区在线 | 一区二区三区四区激情视频| 精品亚洲乱码少妇综合久久| 亚洲精品aⅴ在线观看| 国产亚洲91精品色在线| 国产中年淑女户外野战色| 全区人妻精品视频| 一区二区三区四区激情视频| 免费黄频网站在线观看国产| 久久久久久九九精品二区国产| 日日摸夜夜添夜夜爱| 一二三四中文在线观看免费高清| 中文在线观看免费www的网站| www.av在线官网国产| 少妇丰满av| 嫩草影院入口| 91精品国产九色| 日韩人妻高清精品专区| 少妇裸体淫交视频免费看高清| videos熟女内射| 免费看a级黄色片| 我要看日韩黄色一级片| 日本三级黄在线观看| 午夜福利在线观看免费完整高清在| 偷拍熟女少妇极品色| 久久精品国产鲁丝片午夜精品| 亚洲一区二区三区欧美精品 | 久久久午夜欧美精品| 精品久久国产蜜桃| 九九爱精品视频在线观看| 日韩一本色道免费dvd| av国产免费在线观看| 热99国产精品久久久久久7| 女人被狂操c到高潮| 精品一区二区三卡| 久久99热6这里只有精品| 久久鲁丝午夜福利片| 免费高清在线观看视频在线观看| 又粗又硬又长又爽又黄的视频| 久久久久久久大尺度免费视频| 人妻制服诱惑在线中文字幕| 亚州av有码|