999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于熵權(quán)法賦權(quán)的DTOPSIS法和模糊評價法綜合評價谷子區(qū)試品種

2023-12-11 09:20:00宋中強張文川王帥張揚李龍楊雨陽劉金榮
江蘇農(nóng)業(yè)科學 2023年16期
關(guān)鍵詞:綜合評價

宋中強 張文川 王帥 張揚 李龍 楊雨陽 劉金榮

摘要:應(yīng)用基于熵權(quán)法賦權(quán)的DTOPSIS法和模糊評價法,對全國谷子品種區(qū)域適應(yīng)性聯(lián)合鑒定華北夏谷區(qū)常規(guī)組安陽試點參試的22個品種(系)進行綜合分析評價,利用參試品種的綜合得分,對參試品種的綜合農(nóng)藝性狀進行優(yōu)劣排序,并比較了2種方法的優(yōu)劣,旨在為鑒定和評價谷子品種(系)提供新方法。利用熵權(quán)法確定考察性狀的權(quán)重值,從而對DTOPSIS法和模糊評價法的分析結(jié)果進行比較,結(jié)果表明,考察的9個性狀中,產(chǎn)量、產(chǎn)量構(gòu)成性狀(單穗粒質(zhì)量、單穗質(zhì)量、千粒質(zhì)量)、其他性狀的權(quán)重分別為0.369 9、0.515 6、0.114 5;DTOPSIS法計算出的綜合評價Ci值在品種間存在顯著差異,綜合性狀的差異表現(xiàn)較充分;而模糊評價法計算出的綜合評價D值在品種間差異不大;相關(guān)性分析表明,DTOPSIS法、模糊評價法計算的品種優(yōu)劣排序間具有高度一致性,上述2種方法與僅依據(jù)品種產(chǎn)量高低排序的一致性則相對較差。以上結(jié)果說明,熵權(quán)法能夠科學合理地確定評價指標的權(quán)重,基于熵權(quán)法賦權(quán)的DTOPSIS法、模糊評價法都能對谷子品種(系)作出科學的綜合評價,模糊評價法計算簡潔,但評價精度較DTOPSIS法低。

關(guān)鍵詞:谷子;熵權(quán)法;DTOPSIS法;模糊評價法;綜合評價

中圖分類號:S515.037文獻標志碼:A

文章編號:1002-1302(2023)16-0049-05

收稿日期:2022-11-22

基金項目:國家谷子高粱產(chǎn)業(yè)技術(shù)體系建設(shè)專項(編號:CARS-06-13.5-B25);河南省谷子育種創(chuàng)新團隊(編號:C20150043);河南省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)體系建設(shè)專項(編號:Z2020-14-01)。

作者簡介:宋中強(1983—),男,河南安陽人,副研究員,研究方向為谷子遺傳育種。E-mail:279941169@qq.com。

通信作者:劉金榮,研究員,研究方向為谷子遺傳育種。E-mail:liujinrong63@sohu.com。

谷子是我國北方地區(qū)的主栽作物之一,因其營養(yǎng)價值豐富而深受人們喜愛。作為谷子的發(fā)源國,我國谷子的種質(zhì)資源十分豐富,國家種質(zhì)資源庫保存有各類谷子種質(zhì)達2.7萬多份[1-2,同時存在如何根據(jù)育種需求而合理利用這些種質(zhì)的難題。在實際工作中,多數(shù)育種家根據(jù)產(chǎn)量表現(xiàn)來判斷種質(zhì)的利用價值,但是當前有很多試驗產(chǎn)量較高的品種,在大田生產(chǎn)中表現(xiàn)卻不盡人意,推廣困難。這是因為,隨著谷子生產(chǎn)的發(fā)展,對谷子新品種的要求日益提高,除了要求高產(chǎn)外,更要求穩(wěn)產(chǎn)、抗病、抗逆。因此,使用適當?shù)姆椒▽茸臃N質(zhì)進行綜合評價,是發(fā)掘優(yōu)秀種質(zhì)資源的重要手段。

DTOPSIS法是用于區(qū)域經(jīng)濟發(fā)展的動態(tài)多目標決策的新方法,改進自TOPSIS法[3,該方法根據(jù)評價對象與理想化目標的接近程度進行排序,并根據(jù)評價對象的相對優(yōu)劣順序進行多目標決策,目前在小麥、水稻、玉米、棉花、大豆等農(nóng)作物綜合評價中應(yīng)用廣泛4-8。近年來,模糊評價法在作物耐逆性鑒定及主要農(nóng)藝性狀綜合評價中應(yīng)用廣泛,該方法利用模糊數(shù)學理論,結(jié)合農(nóng)作物評價指標和理論,可以獲得合理、可靠的評價結(jié)果9-11。在作物主要農(nóng)藝性狀綜合評價中,DTOPSIS法、模糊評價法均需要確定影響作物綜合表現(xiàn)的主要因素,然后確定各因素的權(quán)重并計算綜合得分,來確定評價對象的優(yōu)劣,因此權(quán)重對評價結(jié)果有重要影響。熵權(quán)法是一種客觀賦權(quán)法,在信息論中用熵度量數(shù)據(jù)提供的有效信息量,用熵權(quán)表示評價指標的相對重要程度;評價指標的差異程度越大,則對評價結(jié)果的貢獻也越大,對應(yīng)的熵越小,熵權(quán)相應(yīng)越大[12。熵權(quán)法在谷子中的相關(guān)研究報道較少,但該方法在小麥、玉米、大豆、棉花等農(nóng)作物綜合評價中有較多應(yīng)用,并取得了較好的結(jié)果13-16

本研究擬將基于熵權(quán)法賦權(quán)的DTOPSIS法、模糊評價法用于谷子綜合農(nóng)藝性狀的評價,并比較DTOPSIS法、模糊評價法在谷子綜合評價中的優(yōu)劣,為谷子種質(zhì)綜合評價篩選提供更加科學合理的方法。

1 材料與方法

1.1 材料

供試材料為2021年全國谷子品種區(qū)域適應(yīng)性聯(lián)合鑒定華北夏谷區(qū)常規(guī)組安陽試點參試的22個品種(系),詳見表1。

1.2 試驗方法

試驗地點位于安陽市農(nóng)業(yè)科學院柏莊試驗基地(36.6°N,114.2°E)。試驗采用隨機區(qū)組排列,重復3次,小區(qū)行長8.33 m、行寬2.4 m,每小區(qū)6行,行距0.4 m,計產(chǎn)面積19.99 m2。2021年6月3日精量播種,種肥同播,復合肥(N、P2O5、K2O各占15%)用量600 kg/hm2,播后噴灌,一播全苗。各小區(qū)田間管理方式、方法一致,整個生育期內(nèi)降水量較常年偏多,播后未灌溉,拔節(jié)—孕穗期追施尿素 300 kg/hm2,人工除草、防治病蟲害各2次,成熟后及時收獲。

1.3 考察性狀

考察性狀包括生育期、株高、穗長、穗粗、單穗質(zhì)量、單穗粒質(zhì)量、出谷率、千粒質(zhì)量、產(chǎn)量等9個代表性性狀,涵蓋植物學性狀、產(chǎn)量性狀等主要指標。性狀調(diào)查均按照國家谷子區(qū)試田間記載標準進行。

1.4 分析方法

1.4.1 數(shù)據(jù)標準化

為保證各性狀間具有等效性和同序性,對考察的9個性狀進行標準化處理。其中生育期、株高為中性指標,以參試品種該指標的中位數(shù)為最優(yōu)值,采用中性指標測度方法進行標準化;其他性狀的最優(yōu)值為該指標的極大值,采用正向指標測度。以下公式中Xij(k)為標準化后的值,Xij′(k)為標準化前的值,X0(k)為最優(yōu)值,m為性狀數(shù),m=1~9。使用Excel 2010對數(shù)據(jù)進行處理(下同),計算結(jié)果見表2。

正向指標:Xij(m)=Xij′(m)/X0(m);

中性指標:Xij(m)=X0(m)/[X0(m)+|X0(m)-Xij′(m)|]。

1.4.2 基于熵權(quán)法賦權(quán)的指標權(quán)重確定方法

(1) 計算熵值

n為參試品種數(shù)量;并定義當fij=0時,fijlnfij=0。

(2) 計算熵權(quán)

1.4.3 DTOPSIS法 把22個品種(系)的9個性狀指標標準化后,根據(jù)以下程序計算參試品種(系)的綜合得分Ci值:

(1)決策矩陣R:Rij=Wi×Xij(k),式中Wi為各性狀權(quán)重(表3),Xij(k)為標準化后的值。

(2)品種各性狀與最佳性狀的距離S+i及與最差性狀間的距離S-i

式中:Y+i、Y-i分別為決策矩陣R中各性狀的最大值、最小值。

(3)各品種對理想解的相對接近度Ci=S-i/(S+i+S-i),式中Ci∈[0,1],i=1~22,計算結(jié)果見表4。

1.4.4 模糊評價法

根據(jù)模糊數(shù)學原理,利用原始數(shù)據(jù)標準化值和熵權(quán)法確定的權(quán)重,計算參試品種(系)的綜合得分D,根據(jù)D值大小,確定參試品種(系)綜合性狀優(yōu)劣,該方法模型為:

式中μ(k)為原始數(shù)據(jù)標準化值,Wi為各指標權(quán)重,n=1~9。計算結(jié)果見表4。

2 結(jié)果與分析

2.1 熵權(quán)法評價指標權(quán)重的確定

由表3可知,產(chǎn)量在9個考察性狀中權(quán)重最高(0.369 9),其次為單穗粒質(zhì)量(0.234 4)、單穗質(zhì)量(0.200 4)等2個產(chǎn)量構(gòu)成性狀,產(chǎn)量及產(chǎn)量構(gòu)成性狀權(quán)重系數(shù)共計為0.804 6,這與生產(chǎn)上以高產(chǎn)、穩(wěn)產(chǎn)為首要考察目標的要求一致。

2.2 DTOPSIS法的結(jié)果

由表4可知,DTOPSIS法綜合評價的結(jié)果排序與根據(jù)產(chǎn)量進行的排名既有一致性又存在差異性。其中濟谷28、中雜谷34、衡2020-2在2種方法下排序相同;中雜谷36、中雜谷55、冀雜H19ZX1、張雜谷26號的產(chǎn)量排名分別為第1、3、8、12位,而DTOPSIS法計算的排名則下降為第4、6、12、21位,豫谷42、19HQ312、豫谷37的則由產(chǎn)量排名的第4、15、19位,上升為DTOPSIS法綜合評價排名的第1、9、10位,其他品種在2種方法下排名整體差異不大。

2.3 模糊評價結(jié)果

與DTOPSIS法的分析結(jié)果相似,模糊函數(shù)綜合評價值排序與產(chǎn)量排序也存在一定的差異性和一致性。中雜谷36、中雜谷55、冀雜H19ZX1、張雜谷26號由產(chǎn)量排名的第1、3、8、12位下降為綜合評價排名的第6、7、16、21位,豫谷18、19HQ312、豫谷37、滄1057則由產(chǎn)量排名的第6、15、19、22位上升為綜合評價排名的第3、8、9、11位,其他品種在2種方法下排名基本一致。

2.4 DTOPSIS分析法和模糊評價法結(jié)果的比較

由表4可知,利用DTOPSIS法和模糊評價法對各參試品種(系)進行優(yōu)劣排序僅有細微的差異,整體表現(xiàn)出了較高的一致性。基于DTOPSIS法、模糊評價法的品種綜合得分排名的秩相關(guān)系數(shù)為0.946 8,達極顯著相關(guān)水平(P<0.001)(表5)。對DTOPSIS法、模糊評價法的綜合得分與產(chǎn)量排名進行相關(guān)性分析,秩相關(guān)系數(shù)分別為0.833 9、0.705 8,均達極顯著相關(guān)水平,這與產(chǎn)量及產(chǎn)量相關(guān)性狀權(quán)重占比較高有關(guān)。需要指出的是,張雜谷26號、衡17-593的產(chǎn)量及產(chǎn)量排名相同,但依據(jù)DTOPSIS法和模糊評價法計算的綜合得分排名與其產(chǎn)量排名有較大差異,且基于DTOPSIS法和模糊評價法的品種排名基本一致,說明這2種綜合評價方法均能有效統(tǒng)籌品種的主要性狀,較僅根據(jù)產(chǎn)量高低評價品種優(yōu)劣更加合理。

對DTOPSIS法和模糊評價法的分析結(jié)果進一步比較,發(fā)現(xiàn)模糊評價法計算出的品種間綜合評價D值差異較小,D值的變異系數(shù)為8.59%,其最大值與最小值間的差異為41.71%;而DTOPSIS法品種間的Ci值變異系數(shù)為30.11%,最大值與最小值間的差異為412.52%,表現(xiàn)出了較模糊評價法具有更強的分辨力。由此可以看出,從評價精度上來說,DTOPSIS法的綜合評價效果比模糊評價法更好。

3 討論與結(jié)論

權(quán)重是決定綜合評價結(jié)果準確性和可靠性的關(guān)鍵因素。熵權(quán)法是一種客觀賦權(quán)法,可以避免專家賦權(quán)帶來的主觀偏差,使評價結(jié)果更加可靠。該方法根據(jù)熵值大小判斷評價指標的離散程度,熵值越小,說明該指標的離散程度越大,提供的有效信息量越大,則該指標的權(quán)重也應(yīng)較大,反之則權(quán)重較小。本研究利用熵權(quán)法對參試品種(系)的9個評價指標確權(quán),可以得到權(quán)重大小順序為產(chǎn)量>單穗粒質(zhì)量>單穗質(zhì)量>千粒質(zhì)量>穗長>出谷率>穗粗>生育期>株高,產(chǎn)量權(quán)重為0.369 9,產(chǎn)量構(gòu)成要素(單穗粒質(zhì)量、單穗質(zhì)量、千粒質(zhì)量)權(quán)重為0.515 6,其他性狀權(quán)重為0.114 5,產(chǎn)量構(gòu)成要素權(quán)重較高,其次為產(chǎn)量,這與生產(chǎn)上重視谷子品種的高產(chǎn)、穩(wěn)產(chǎn)特性的實際情況一致,說明本研究權(quán)重設(shè)置合理,在此基礎(chǔ)上能夠?qū)⒃嚬茸悠贩N做出科學合理的綜合評價。

對參試品種的綜合評價結(jié)果進行分析,DTOPSIS法、模糊評價法這2種不同評價方法的優(yōu)先度排名極顯著相關(guān)(r=0.946 8,P<0.001),這表明2種方法的分析結(jié)果一致性較強,分析方法可靠。從表4可以看出,張雜谷26號、衡17-593的產(chǎn)量和排名相同,但二者依據(jù)DTOPSIS法、模糊評價法的綜合得分排名與其產(chǎn)量排名存在較大差異,且前2種分析方法的排名基本一致。由此可以看出,僅依據(jù)產(chǎn)量對品種(系)進行評價具有很大的局限性,而DTOPSIS法、模糊評價法可以較全面地利用考察性狀指標對參試谷子品種(系)進行合理綜合評價。對DTOPSIS法、模糊評價法的結(jié)果進一步分析發(fā)現(xiàn),模糊評價法品種間的D值差異較DTOPSIS法品種間的Ci值差異小,表明在相同的權(quán)重設(shè)置條件下,DTOPSIS法對品種綜合性狀優(yōu)劣的區(qū)分能力更強,這一結(jié)論在眾多學者利用DTOPSIS法對作物農(nóng)藝性狀進行綜合評價中也得到了驗證[17-19。但相較于 DTOPSIS法繁雜的計算過程,模糊函數(shù)法的計算程序簡潔,分析結(jié)果也能滿足實際需求。

通過對22個參試品種(系)的綜合評價分析,本研究認為DTOPSIS法評價精度相對模糊函數(shù)法更高,但二者都能較好地對參試品種(系)進行綜合評價,對谷子育種中大量材料的鑒定有一定參考價值。本研究對參試品種(系)的評價指標設(shè)定了最優(yōu)值,并針對不同指標的特點設(shè)置了正向指標、中性指標進行計算,保證了評價結(jié)果的合理性。但需要指出的是,本研究所設(shè)定的最優(yōu)值僅限參與此次評價的22個品種(系),其代表性較弱。因此各谷子育種單位有必要加強對“谷子理想株型”的研究,將符合各地生產(chǎn)實際需要的“理想品種”的指標設(shè)定為最優(yōu)值,用于對后代材料的評價和篩選。此外,盡管熵權(quán)法等客觀賦權(quán)法能夠較好地反映評價指標對評價結(jié)果的貢獻,但得到的權(quán)重僅考慮了某個指標的特征,賦權(quán)結(jié)果不一定能與育種目標較好地吻合。因此在育種過程中,也不能僅根據(jù)指標數(shù)量特征進行客觀賦權(quán),而要綜合統(tǒng)籌育種目標,在客觀賦權(quán)的基礎(chǔ)上,對選取的考察性狀權(quán)重進行一定調(diào)整,這樣才能更好地發(fā)揮DTOPSIS法或模糊評價法在品系篩選中的作用,減少品系篩選的盲目性。

熵權(quán)法賦權(quán)能較好地反映評價指標的重要程度,基于熵權(quán)法賦權(quán)的DTOPSIS法、模糊評價法都能對參試品種(系)的綜合性狀作出合理、可靠的評價。DTOPSIS法的分辨能力相較模糊評價法更好,但模糊評價法使用更方便,在實際應(yīng)用中可根據(jù)評價精度的要求選擇合適的方法。

參考文獻:

[1]Li Y,Wu S Z,Cao Y S,et al. A phenotypic diversity analysis of foxtail millet [Setaria italica (L.) P.Beauv.] landraces of Chinese origin[J]. Genetic Resources and Crop Evolution,1996,43(4):377-384.

[2]劉 旭,黎 裕,曹永生,等. 中國禾谷類作物種質(zhì)資源地理分布及其富集中心研究[J]. 植物遺傳資源學報,2009,10(1):1-8,171,169.

[3]姚興濤,朱永達,張永貞,等. 區(qū)域國民經(jīng)濟協(xié)調(diào)發(fā)展的動態(tài)多目標決策[J]. 農(nóng)業(yè)工程學報,1991,7(4):1-6.

[4]陳 慧,王冀川. 基于灰色關(guān)聯(lián)度及DTOPSIS法的南疆冬小麥品種的綜合評判[J]. 江蘇農(nóng)業(yè)科學,2019,47(10):102-108.

[5]蔣 聰,劉慰華,楊旭昆,等. 灰色關(guān)聯(lián)度分析和DTOPSIS法在云南粳稻品種綜合評價中的應(yīng)用[J]. 西南農(nóng)業(yè)學報,2020,33(5):907-912.

[6]許 波,王成業(yè),張海申,等. DTOPSIS法在玉米新品種多因素綜合評價中的應(yīng)用[J]. 種子,2005,24(1):14-16.

[7]郭克婷,龍騰芳,楊付新.利用DTOPSIS法綜合評價黃河流域棉花區(qū)試新品種[J]. 中國農(nóng)學通報,2004,20(2):226-228,239.

[8]劉自華,董玉武,薛 晶,等. DTOPSIS分析法應(yīng)用于河北大豆新品種的評價[J]. 北京農(nóng)學院學報,2005,20(3):20-23.

[9]劉自華,龔振平. 用模糊數(shù)學綜合評價農(nóng)作物品種區(qū)域試驗初探[J]. 北京農(nóng)學院學報,1998,13(4):15-20.

[10]單鑒華,馮潤生. 模糊綜合評判法和模糊評分法在作物品種綜合評價中的應(yīng)用[J]. 浙江農(nóng)業(yè)科學,1996,37(3):119-121.

[11]陶維旭,程生海,冀俊超,等. 水稻品種資源耐鹽性綜合評價及耐鹽指標篩選[J]. 江蘇農(nóng)業(yè)科學,2022,50(18):180-187.

[12]李 萍,魏朝富,邱道持. 基于熵權(quán)法賦權(quán)的區(qū)域耕地整理潛力評價[J]. 中國農(nóng)學通報,2007,23(6):536-541.

[13]黃 杰,王 君,宋丹陽,等. 基于熵權(quán)法的DTOPSIS分析法與灰色關(guān)聯(lián)度分析法對‘漯麥906的綜合評價[J]. 天津農(nóng)業(yè)科學,2021,27(11):20-25.

[14]軒瑞瑞,陳艷萍,劉春菊,等. 基于熵權(quán)法和灰色關(guān)聯(lián)度法的鮮食糯玉米品質(zhì)評價[J]. 食品工業(yè)科技,2021,42(14):241-248.

[15]安 東. 33個南方籽粒型大豆品種在黃土高原的飼用潛力綜合評價[D]. 蘭州:蘭州大學,2020.

[16]張國偉,劉瑞顯,楊長琴,等. 基于熵權(quán)理論的模糊綜合評判法在棉花耐鹽性評價中的應(yīng)用[J]. 農(nóng)業(yè)科學與技術(shù),2014,15(9):1441-1445,1481.

[17]李彥平,李淑君,吳娟霞,等. DTOPSIS法和灰色關(guān)聯(lián)度法在新引烤煙新品種綜合評價中的應(yīng)用比較[J]. 中國煙草學報,2012,18(4):35-40.

[18]昝 凱,周 青,張志民,等. 灰色關(guān)聯(lián)度和DTOPSIS法綜合分析河南區(qū)域試驗中大豆新品種(系)的農(nóng)藝性狀表現(xiàn)[J]. 大豆科學,2018,37(5):664-671.

[19]宋秀麗,吳 玥,楊錫朗,等. 基于熵值賦權(quán)的DTOPSIS法與灰色關(guān)聯(lián)度分析在玉米品種綜合評價中的比較[J]. 玉米科學,2020,28(2):41-46.

猜你喜歡
綜合評價
和諧社會統(tǒng)計指標體系研究述評
吉林省居民生活質(zhì)量評價
基于AHP的織金洞智慧旅游資源條件綜合評價研究
綠色科技(2016年21期)2016-12-27 16:20:15
完善新型職業(yè)農(nóng)民培訓模式 推動我國農(nóng)業(yè)健康有序發(fā)展
運用目標分析最優(yōu)指標法綜合評價廣西市某醫(yī)院2011—2015醫(yī)療質(zhì)量
陜西省各地區(qū)人力資本水平綜合評價與分析
10kV配電線路帶電作業(yè)安全綜合評價應(yīng)用探究
主成分分析法在大學英語寫作評價中的應(yīng)用
大學教育(2016年11期)2016-11-16 20:33:18
高校創(chuàng)新人才評價機制探索
成才之路(2016年29期)2016-10-31 16:05:39
鄭州市各縣(市)創(chuàng)新能力綜合評價
主站蜘蛛池模板: 欧美一级夜夜爽www| 国产一区二区影院| 99在线观看视频免费| 久久久久亚洲AV成人人电影软件| 亚洲av片在线免费观看| 亚洲色图欧美在线| 亚洲第一黄片大全| 免费观看成人久久网免费观看| 亚洲Aⅴ无码专区在线观看q| 日韩精品一区二区三区大桥未久| 性喷潮久久久久久久久| 国产成人资源| 国产精品所毛片视频| 亚洲无限乱码一二三四区| 国产精品永久久久久| 日韩一区二区三免费高清| 国产乱肥老妇精品视频| 91美女在线| 国产精品自拍露脸视频| 亚洲最大看欧美片网站地址| 午夜国产精品视频| 久久超级碰| 色综合中文综合网| 国产小视频a在线观看| www.国产福利| 亚洲无码久久久久| 亚洲丝袜中文字幕| 亚洲品质国产精品无码| 免费国产一级 片内射老| 国产福利观看| 中文字幕在线免费看| 国产精品美女免费视频大全| 国产欧美视频综合二区 | 久久综合成人| 99久久精品免费看国产免费软件| 中文字幕乱妇无码AV在线| 亚洲成人高清无码| 国产精品成人观看视频国产 | 久久永久免费人妻精品| 久久久国产精品免费视频| 成人亚洲国产| 国产成人无码久久久久毛片| 国产鲁鲁视频在线观看| 中文字幕av无码不卡免费| 熟妇丰满人妻| 日韩国产亚洲一区二区在线观看| 国产精品自在线天天看片| 久久综合婷婷| 无码在线激情片| 国产成年无码AⅤ片在线| 免费一级成人毛片| 欧美亚洲欧美| 超级碰免费视频91| 久久福利片| 国产精品久久国产精麻豆99网站| 成人在线天堂| 制服丝袜一区| 免费A∨中文乱码专区| 亚洲高清在线天堂精品| 国产精品亚洲一区二区在线观看| 国产福利在线免费观看| 在线一级毛片| 婷婷激情亚洲| 日韩区欧美国产区在线观看| 国产美女91呻吟求| 九色在线观看视频| 99ri国产在线| 欧洲av毛片| 漂亮人妻被中出中文字幕久久| 制服丝袜国产精品| 高清无码一本到东京热| 日韩天堂在线观看| 精品丝袜美腿国产一区| 无码综合天天久久综合网| 中文一级毛片| 91精品视频在线播放| 国产一二三区视频| 国产亚洲美日韩AV中文字幕无码成人| 亚洲免费黄色网| 午夜毛片免费观看视频 | 美女无遮挡免费视频网站| 在线中文字幕日韩|