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

DTOPSIS和灰色局勢決策法在谷子綜合評價中的應用

2024-02-14 00:00:00王淑君邢璐劉俊芳閆宏山王素英
江蘇農業科學 2024年24期
關鍵詞:綜合評價

摘要:為了探索適宜的谷子綜合評價方法,篩選綜合性狀優異的谷子新品種,以東北春谷區的20個谷子新品種為試驗材料,通過各性狀灰色關聯度所占比重為權重,分別采用DTOPSIS和灰色局勢決策法對參試品種進行綜合評價,并對2種綜合評價方法進行對比分析。結果表明,DTOPSIS法綜合評價Ci排在前7位的谷子品種為安20hN-4029、赤333-26、九谷45、安H081、朝 202026、瑞香谷5號、中雜谷78,品種間Ci最大差異達65.73%;灰色局勢決策法分析,ri排在前7位的谷子品種為九谷45、安H081、瑞香谷5號、中雜谷78、赤333-26、20hN-4029、朝202026,品種間ri最大差異為13.64%;相關性分析結果顯示,Ci與ri排序結果存在極顯著正相關關系,秩相關系數為0.922 (Plt;0.01);由此可見,DTOPSIS和灰色局勢決策法均能客觀、全面地評價谷子品種,但DTOPSIS法更能表現出品種間差異,反映谷子品種優劣,綜合評價效果優于灰色局勢決策法。依據參試品種綜合評價結果和產量表現篩選出九谷45、朝202026和瑞香谷5號這3個綜合性狀優異且豐產性好的谷子品種。

關鍵詞:谷子;關聯度;DTOPSIS法;灰色局勢決策法;綜合評價

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

文章編號:1002-1302(2024)24-0066-06

收稿日期:2023-12-25

基金項目:河南省重點研發專項(編號:231111110300);河南省農業良種聯合攻關項目(編號:2022010401-5);現代農業產業技術體系建設專項(編號:CARS-06-14.5-B25)。

作者簡介:王淑君(1982—),女,河南安陽人,碩士,副研究員,主要從事谷子遺傳育種研究。E-mail:logccc@163.com。

谷子起源于我國,是中華民族的哺育作物,在北方干旱半干旱地區,對于農業的可持續性發展和作物生態多樣性建設具有重要意義[1。近年來,隨著谷子產業發展,培育高產、優質、廣適等綜合性狀優異的谷子新品種已成為育種工作者的首要目標。依舊只依靠產量對谷子品種進行評價已不能滿足現代育種的要求。因此,選擇合適的綜合評價方法,科學、客觀、準確地評價谷子品種,是新品種選育和推廣的關鍵。 逼近理想解的排序方法(DTOPSIS)是一種新興的多目標決策分析方法,依據評價目標與理想值的接近程度進行排序來評價目標優劣[2-4;灰色局勢決策是利用目標效果測度將局勢中多個目標值轉化為單個目標值,從定量的角度進行方案優選5。二者在玉米6-10、水稻11-12、小麥[13、油菜14、大豆15-18、棉花19、草莓20等各類作物品種評價研究中廣泛應用。在谷子品種的綜合評價中,張曉申等曾采用主觀賦權的DTOPSIS法對11個谷子品種進行綜合評價,研究認為賦予合理權重的DTOPSIS法對谷子綜合評價才更有意義[21;宋中強等曾采用熵權賦值的DTOPSIS法對安陽試點的22個谷子品種進行綜合評價,但評價指標中并未涉及逆向指標[22。本研究以2022年東北春谷區的20個谷子品種為試驗材料,通過各性狀灰色關聯度所占比例,對谷子品種的11個主要性狀指標進行賦權,然后分別采用DTOPSIS法和灰色局勢決策法對參試品種進行綜合評價,并對2種評價方法進行比較,以期探索適宜的谷子綜合評價方法,篩選出綜合性狀優異的谷子品種。

1 材料與方法

1.1 試驗材料

參試材料為2022年度全國谷子品種區域適應性聯合鑒定試驗(東北春谷區組)的20個谷子品種:瑞香谷 2 號(X1)、瑞香谷 5 號(X2)、朝 202026(X3)、龍谷 46(X4)、九谷 40(X5)、九谷 41(X6)、K175-2H(X7)、赤金谷 17(X8)、公谷 96(X9)、冀白米1號(X10)、九谷 45(X11)、冀 380(X12)、鐵谷 15(X13)、赤谷C1(X14)、赤 333-26(X15)、中雜谷 78(X16)、安 20hN-4029(X17)、安 20N-3166(X18)、安 H081(X19)、九谷11(X20,對照)。

1.2 試驗設計與方法

試驗在遼寧、吉林、內蒙古、黑龍江等4個省份共設11個試點。均采用完全隨機區組設計,重復3次,小區面積為20 m2,四周設保護行;田間管理嚴格按照國家谷子品種區域適應性聯合鑒定試驗方案執行。調查記載20個參試品種的生育期(T1)、株高(T2)、穗長(T3)、穗粗(T4)、單穗重(T5)、穗粒重(T6)、出谷率(T7)、千粒重(T8)、產量(T9)、白發病發病率(T10)和蛀莖率(T11)等11個農藝性狀,試驗數據采用11個試點各個性狀的平均值。

1.3 數據統計與分析方法

1.3.1 灰色關聯度分析法確定權重 依據灰色系統理論[23,將供試品種及其所有性狀看作一個灰色系統,每個性狀為該系統中的一個因素。以產量作為參考數列X0,其他10個性狀指標為比較數列Xi,利用公式(1)和公式(2)計算各品種產量與其他性狀關聯度,依據公式(3)經歸一化處理后得到各性狀權重系數。

εi=miniminj |X0(j)-Xi(j)|+ρmaximaxj|X0(j)-Xi(j) ||X0(j)-Xi(j)|+ρmaxi maxj |X0(j)-Xi(j)|;(1)

Ri=1n ∑ni=1εi;(2)

w=Ri/∑ni=1Ri。(3)

式中:εi為關聯系數;Ri為關聯度;w為權重;ρ為分辨系數(ρ=0.5);|X0(j)-Xi(j)|為X0 數列與Xi 數列在第j點的絕對差值;miniminj|X0(j)-Xi(j)|為二級最小差值;maximaxj |X0(j)-Xi(j)|為二級最大差值。

1.3.2 DTOPSIS法 參照葉開梅等的統計方法[24進行計算,步驟如下:

(1) 構建評價矩陣,將性狀指標分別按照下列公式進行無量綱化處理,得到標準化矩陣Z:

正向指標:Zij=Yij/Yjmax;(4)

中性指標:Zij=Yj0/[Yj0+|Yj0-Yij |];(5)

負向指標:Zij=Yjmin/Yij。(6)

式中:Yjmax為第j個性狀最大值;Yjmin為第j個性狀最小值;Yj0為第j個性狀平均值。

(2) 用各性狀的權重值wj乘以標準化矩陣Z,得到決策矩陣R;根據公式(7)和公式(8)計算各性狀指標的正、負理想解:

正理想解X+j={maxRij}={ X+1,X+2,…,X+n};(7)

負理想解X-j={ minRij }={ X-1,X-2,…,X-n}。(8)

(3) 相對近似度計算

根據歐幾里德范數,利用公式(9)和公式(10)計算各品種性狀與正、負理想解的距離,進而得到不同處理的相對近似值。根據Ci值對品種排序,Ci值越大則表示該品種性狀指標與理想解距離越近,綜合性狀表現越優異,反之則越差25

S+i=2∑nj=1(Rij-X+j)2;(9)

S-i=2∑nj=1(Rij-X-j)2;(10)

Ci=S-i/(S+i+S-i)。(11)

1.3.3 灰色局勢決策法 參照郭永忠等的方法[6,將性狀分為上限效果測度、適中效果測度、下限效果測度,得到效果測度矩陣L。根據各性狀權重(wj)和公式(12),計算出各參試品種的加權綜合效果測度值,即為各品種的綜合優勢量化值ri,ri值越大,品種綜合性狀越優異。

ri=∑ni=1wjLij。(12)

2 結果與分析

2.1 參試品種產量表現

由表1可知,20個參試品種的產量在4 575.00~5 644.50 kg/hm2之間,由高到低排序依次為九谷45gt;九谷41gt;龍谷46gt;朝202026gt;瑞香谷5號gt;安20N-3166gt;九谷40gt; K175-2H gt;冀白米1號gt;瑞香谷2號gt;安H081gt;公谷96gt;赤333-26gt;九谷11gt;鐵谷15gt;中雜谷78gt;安20hN-4029gt;冀380gt;赤谷C1gt;赤金谷17。產量超過對照九谷11的品種共有13個,其中九谷45、九谷41、龍谷46、朝202026、瑞香谷5號、安20N-3166、九谷40排在參試品種前7位,豐產性較好;中雜谷78、安 20hN-4029、冀380、赤谷C1和赤金谷17排名靠后,豐產性較差。

2.2 權重系數的確定

各性狀指標采用Z-score標準化法進行處理[26,然后利用公式(1)~公式(3)計算產量與其他性狀的關聯度,進而得到各性狀權重系數。由表2可知,產量與其他10個性狀的關聯度大小表現為出谷率gt;穗粒重gt;單穗重gt;蛀莖率gt;白發病發病率gt;株高gt;千粒重gt;生育期gt;穗粗gt;穗長,說明出谷率、穗粒重和單穗重這3個性狀對產量的影響較大,因此性狀出谷率、穗粒重和單穗重所占權重較高,分別為0.096 3、0.094 3和0.091 6,僅次于產量(0.118 9);穗長和穗粗的權重較低,分別為0.083 2和0.084 1。

2.3 DTOPSIS 法

根據谷子育種目標,在考察的11個性狀指標中,穗長、穗粗、單穗重、穗粒重、出谷率、千粒重和產量8個性狀均為正向指標,值越大越好;生育期和株高屬中性指標,性狀適中為好;白發病發病率和蛀莖率為負向指標,值越小越好。為消除指標間因量綱和數量級產生的差異[27,依據公式(4)~公式(6)進行無量綱化處理,得到規范化矩陣Z(表3)。用各性狀的權重值wj乘以標準化矩陣Z,可得到決策矩陣R,然后根據公式(7)和公式(8)得到11個性狀的正、負理想解,結果見表4;利用公式(9)~公式(11)計算各參試品種與理想解的相對接近度Ci。由表5可知,各參試品種的Ci在0.205 7~0.600 2之間,其中安20hN-4029、赤333-26、九谷45、安H081、朝202026、瑞香谷5號、中雜谷78的Ci值高于對照九谷11,且排在參試品種前7位,說明這7個品種的綜合性狀表現優異。而瑞香谷2號、公谷96、冀380、赤谷C1、赤金谷17的Ci值較小,排在參試品種后5位,綜合性狀表現較差。

2.4 灰色局勢決策法

由表5可知,ri值高于對照九谷11且排在前7位的品種分別為九谷45、安H081、瑞香谷5號、中雜谷78、赤333-26、20hN-4029、朝202026,說明這些谷子品種的綜合性狀表現較好;品種K175-2H、公谷96、冀380、赤谷C1和赤金谷17的ri值位于參試品種后5位,綜合性狀較差。

2.5 不同評價方法評價結果比較

DTOPSIS和灰色局勢決策法因計算方法不同,導致評價結果略有不同。由表5可知,除品種安20hN-4029和九谷11在2種分析方法評價結果中位次差異大于3外,其余品種位次差異均不超過3。其中安20hN-4029、赤333-26、九谷45、安H081、朝202026、瑞香谷5號和中雜谷78這7個品種的2種分析方法評價結果均優于對照九谷11,且排在參試品種前7位,綜合性狀表現優異。而公谷96、冀380、赤谷C1和赤金谷17在2種分析方法中評價結果位次完全一致,Ci和ri排名均在參試品種后5位,綜合性狀表現較差。相關性分析結果顯示,Ci與ri排序結果存在極顯著正相關關系,秩相關系數為0.922(Plt;0.01),評價結果一致性較好。然而,DTOPSIS法分析參試品種的Ci范圍在0.205 7~0.600 2之間,最大差異達65.73%,說明品種優劣差異顯著;而灰色局勢決策法分析參試品種的ri范圍在0.782 4~0.906 0之間,最大差異僅為13.64%,說明灰色局勢決策法分析各品種間差異表現不夠充分,進行量化評價比較的能力較弱,而DTOPSIS法更能表現出谷子品種間的差異,反映品種優劣。

DTOPSIS法和灰色局勢決策法綜合評價結果與產量結果相比,存在明顯差異。例如,中雜谷78的Ci和ri分別排在參試品種第7和第4位,而產量結果排在參試品種第16位;龍谷46產量排名是第3位,而綜合評價結果Ci和ri排名則分別為第15和第14位。由此可見,產量高的品種綜合表現不一定好,而綜合性狀優異的谷子品種豐產性不一定好。依據綜合評價結果和各參試品種產量表現可以看出,赤333-26、九谷45、安H081、朝202026和瑞香谷5號的綜合性狀表現和產量均超過對照九谷11;其中九谷45、朝202026、瑞香谷5號不僅產量結果排在參試品種前5位,而且綜合評價結果排在參試品種前7位,說明這3個品種不僅豐產性好,而且綜合性狀優異;而冀380、赤谷C1和赤金谷17產量低且綜合性狀較差。

3 討論與結論

作物品種綜合評價是新品種選育和推廣的重要前提。然而綜合評價涉及目標性狀較多,如采用方差分析或新復極差等方法,在結果不顯著時難以得到有效結論,導致部分優良品種得不到有效利用[12。DTOPSIS和灰色局勢決策法可以綜合多個目標性狀,并將其轉化為可比較的量化值,對參試品種的優劣進行排序。然而在分析評價過程中,DTOPSIS 法需建立各性狀的正、負理想解,通過計算各性狀與正、負理想解的距離進而得到不用品種的Ci值;而灰色局勢決策法只需依據不同效果測度計算綜合效果測度值,計算過程更為簡便。二者因計算方法不同,致使評價結果存在一定差異。研究結果表明,DTOPSIS 法分析Ci高于對照九谷11,且排在參試品種前7位的品種是安20hN-4029、赤333-26、九谷45、安H081、朝202026、瑞香谷5號、中雜谷78;瑞香谷2號、公谷96、冀380、赤谷C1、赤金谷17的Ci排在參試品種后5位;參試品種間Ci最大差異達65.73%。灰色局勢決策法分析,綜合優勢量化值ri高于對照九谷11且排在前7位的品種分別為九谷45、安H081、瑞香谷5號、中雜谷78、赤333-26、20hN-4029、朝202026;品種K175-2H、公谷96、冀380、赤谷C1和赤金谷17的ri位于參試品種后5位;品種間ri最大差異僅為13.64%,明顯小于DTOPSIS法Ci的最大差異。相關性分析結果顯示,Ci與ri排序結果存在極顯著正相關關系,秩相關系數為0.922(Plt;0.01),說明2種綜合評價方法均能全面、客觀地評價谷子品種,評價結果一致性好,而DTOPSIS法更能反映谷子品種間優劣。二者的聯合運用可以彌補單一分析方法的不足,進一步提高綜合評價的準確性和可靠性。

運用 DTOPSIS 法和灰色局勢決策法評價作物品種關鍵在于評價指標的選擇和各性狀權重系數的確定。評價指標的選擇要根據育種目標和生產實際情況確定,在條件允許的情況下,盡可能多涵蓋不同類型的性狀指標,以減少因性狀選擇而造成的誤差。本研究引用了正向指標、中性指標和負向指標,并對不同類型指標選擇不同的無量綱化處理方法,使得標準化結果更加合理。性狀指標的權重大小反映了該指標的相對重要程度,關系到評價結果的準確性和可靠性[28-30。目前,產量仍是谷子育種的重要目標,通過計算產量與其他性狀的關聯度,進而利用各性狀灰色關聯度所占比重為各性狀指標權重賦值,可有效避免經驗賦權的主觀性,能夠更加高效、準確地反映出各參試品種性狀與產量的主次關系,從而準確、合理、科學地評價谷子新品種優劣,進而為谷子產業發展提供可靠依據。

DTOPSIS和灰色局勢決策法均能對谷子品種進行客觀、全面地評價,但DTOPSIS法較灰色局勢決策分析更能表現出品種間差異,反映谷子品種優劣,而灰色局勢決策計算方法更為簡便,二者聯合應用可使評價結果更加準確。依據參試品種綜合評價結果和產量表現篩選出九谷45、朝202026和瑞香谷5號3個綜合性狀優異且豐產性好的谷子品種。

參考文獻:

[1]刁現民. 中國谷子產業與產業技術體系[M]. 北京:中國農業科學技術出版社,2011.

[2]代順冬,韋樹谷,黃 玲,等. DTOPSIS、AHP及隸屬函數法綜合評價四川烤煙區試品種[J]. 西南師范大學學報(自然科學版),2021,46(6):57-64.

[3]馮建文,韓秀梅,宋 莎,等. 灰色關聯度分析和DTOPSIS法在貴州威寧引進蘋果綜合評價中的應用[J]. 南方農業學報,2022,53(11):3157-3164.

[4]楊海峰,洪德峰,馬俊峰,等. DTOPSIS法在分析不同基因型玉米間混作效應中的應用[J]. 中國農學通報,2023,39(6):13-17.

[5]余文遠,李霖超. 運用灰色局勢決策法綜合評價鮮食大豆新品系[J]. 湖南農業科學,2015(10):26-28.

[6]郭永忠,楊彩霞,鄒 軍,等. 灰色局勢決策在玉米品種綜合評價中的應用[J]. 西北農業學報,2007,16(4):92-95.

[7]劉愛華,蔣 燕,羅 冬,等. 灰色局勢決策在玉米新品種(組合)展示試驗數據分析中的應用[J]. 種子,2006,25(5):77-78,87.

[8]齊建雙,夏來坤,黃 保,等. 基于熵權的DTOPSIS法和灰色局勢決策法在玉米品種區域試驗中的應用探討[J]. 作物雜志,2021(1):60-67.

[9]姚宗澤,伊應良,楊肖艷,等. 基于熵權的DTOPSIS法和灰色局勢決策法對云南雜交玉米新組合的評價分析[J]. 江蘇農業科學,2023,51(6):86-93.

[10]吳 玥,李 威,馬德志,等. 基于熵值賦權的DTOPSIS法對不同玉米品種綜合評價研究[J]. 玉米科學,2019,27(4):32-41.

[11]封光華,郭貴陽,賴克柱,等. 灰色局勢決策在水稻品種分析中的應用[J]. 江西農業學報,2005,17(2):17-20.

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

[13]金 輕,趙 紅,林麗萍,等. 基于灰色關聯度分析和DTOPSIS法綜合評價小麥新品系在云南省的適應性[J]. 南方農業學報,2020,51(10):2440-2446.

[14]程 輝,王友華,何道君,等. 用灰色局勢決策法多因素綜合評估甘藍型油菜早代自交系的一般配合力[J]. 中國油料作物學報,2007,29(3):342-346.

[15]繆亞梅,王學軍,汪凱華,等. 應用灰色局勢決策方法綜合評價鮮食大豆新品種[J]. 河北農業科學,2007,11(4):10-11,14.

[16]侯樂新,郭巍杰,王培田. 灰色局勢決策在大豆新品種評價中的應用初探[J]. 河南農業科學,2004,33(9):31-33.

[17]王秋玲,劉 艷,谷傳彥,等. 試用灰色局勢決策原理綜合評價大豆新品種[J]. 植物遺傳資源科學,2002,3(4):28-31.

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

[19]蘇天增. 灰色局勢決策在棉花品種評估中的應用[J]. 中國棉花,2004,31(10):19-20.

[20]李文硯,韋 優,孔方南,等. DTOPSIS法在草莓品種綜合評價中的應用研究[J]. 植物生理學報,2018,54(5):925-930.

[21]張曉申,韓燕麗,樊永強,等. 基于灰色關聯度和DTOPSIS法對谷子區域試驗的綜合評價[J]. 種子,2022,41(9):121-126,133.

[22]宋中強,張文川,王 帥,等. 基于熵權法賦權的DTOPSIS法和模糊評價法綜合評價谷子區試品種[J]. 江蘇農業科學,2023,51(16):49-54.

[23]王淑君,解慧芳,邢 璐,等. 基于主成分分析和灰色關聯度分析的谷子品種(系)綜合評價[J]. 江蘇農業科學,2023,51(15):42-49.

[24]葉開梅,牛力立,樊祖立,等. 灰色關聯度法與DTOPSIS法在馬鈴薯品種綜合評價中的應用比較[J]. 貴州農業科學,2023,51(3):10-18.

[25]邵美紅,程 楚,程思明,等. 運用DTOPSIS法對鮮食甜玉米新品種在浙西北地區適應性的綜合評價[J]. 江西農業學報,2017,29(6):25-28.

[26]鄭 峰,王忠明,王 迪,等. 基于變異系數權重法評價干燥無花果多糖品質[J]. 保鮮與加工,2019,19(4):119-124.

[27]蔣自可,王淑君,解慧芳,等. 基于灰色關聯度和聚類分析的谷子品種(系)綜合評價與篩選[J]. 農業科技通訊,2023(1):27-30,35.

[28]劉自遠,劉成福. 綜合評價中指標權重系數確定方法探討[J]. 中國衛生質量管理,2006,13(2):44-46,48.

[29]劉 群,馮長煥. 一種新的最優權系數組合賦權法[J]. 蘭州文理學院學報(自然科學版),2015,29(6):1-4.

[30]劉曉玥,杜 燕,鄭小賢. 不同賦權法的將樂林場常綠闊葉林的健康評價[J]. 森林與環境學報,2015,35(2):141-146.

猜你喜歡
綜合評價
和諧社會統計指標體系研究述評
吉林省居民生活質量評價
基于AHP的織金洞智慧旅游資源條件綜合評價研究
綠色科技(2016年21期)2016-12-27 16:20:15
完善新型職業農民培訓模式 推動我國農業健康有序發展
運用目標分析最優指標法綜合評價廣西市某醫院2011—2015醫療質量
陜西省各地區人力資本水平綜合評價與分析
10kV配電線路帶電作業安全綜合評價應用探究
主成分分析法在大學英語寫作評價中的應用
大學教育(2016年11期)2016-11-16 20:33:18
高校創新人才評價機制探索
成才之路(2016年29期)2016-10-31 16:05:39
鄭州市各縣(市)創新能力綜合評價
企業導報(2016年9期)2016-05-26 20:52:48
主站蜘蛛池模板: 亚洲国产欧美国产综合久久 | 亚洲欧美在线看片AI| 欧美精品成人一区二区视频一| 日韩欧美国产精品| 国产主播一区二区三区| 亚洲欧美不卡视频| 爱做久久久久久| 72种姿势欧美久久久大黄蕉| 亚洲色欲色欲www网| 午夜老司机永久免费看片| 国产手机在线ΑⅤ片无码观看| 中文精品久久久久国产网址| 国产在线第二页| 亚洲中文久久精品无玛| 谁有在线观看日韩亚洲最新视频| 亚卅精品无码久久毛片乌克兰| 亚洲一区二区无码视频| 91麻豆精品国产高清在线| 久久综合色天堂av| 欧美v在线| 国产欧美在线视频免费| 亚洲系列无码专区偷窥无码| 国产一级无码不卡视频| 国产成人区在线观看视频| 久久综合激情网| 99视频只有精品| 精品福利视频网| 成人中文在线| 久久 午夜福利 张柏芝| 亚洲日韩精品无码专区| 欧美性天天| 色婷婷色丁香| 素人激情视频福利| 国产地址二永久伊甸园| 97视频在线精品国自产拍| 性视频久久| 国产日本欧美亚洲精品视| 国产理论一区| 伊人久久青草青青综合| 国产精品第三页在线看| 9966国产精品视频| 国产精品视频导航| 超清无码熟妇人妻AV在线绿巨人| 欧美日韩国产成人在线观看| 亚洲天堂网视频| 天天婬欲婬香婬色婬视频播放| 97视频免费在线观看| 天堂在线www网亚洲| 日韩资源站| 亚洲欧美综合在线观看| 国产精品专区第1页| 日本午夜三级| 91成人在线观看视频| 欧美一级在线| 久久精品国产免费观看频道| 国产午夜一级毛片| 91年精品国产福利线观看久久 | 免费A∨中文乱码专区| 色婷婷在线播放| 亚洲国产av无码综合原创国产| 亚洲欧美成人影院| 激情六月丁香婷婷四房播| 91精品国产福利| 大陆精大陆国产国语精品1024| 91福利一区二区三区| 免费日韩在线视频| 成人午夜视频在线| 91成人在线免费视频| 青草娱乐极品免费视频| 色香蕉影院| 久久国产精品影院| 波多野结衣亚洲一区| 欧美精品一二三区| 久久永久视频| 毛片手机在线看| 欧美精品三级在线| 免费看美女毛片| 成人午夜视频网站| 婷婷激情五月网| 夜精品a一区二区三区| 国产内射一区亚洲| 亚洲视频在线青青|