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

基于主成分分析與聚類分析的西芹生長適應(yīng)性評價

2025-07-22 00:00:00南柏林張海杰師麗麗張鶴瀟許春苗程永龍
中國果菜 2025年6期
關(guān)鍵詞:西芹性狀品種

中圖分類號:S636.3 文獻標(biāo)志碼:A 文章編號:008-038(2025)06-0050-06

DOI: 0.9590/j.cnki.008-038.2025.06.009

Evaluation of Growth Adaptability of Celery Based on Principal Component Analysis and Cluster Analysis

NANBailin', ZHANG Haijie2, SHI Lili2, ZHANG Hexiao2,XU Chunmiao3,CHENG Yonglong2 (1.Agricultural Technology Extension Service Center of Anding District,Dingxi 743OO0,China; 2.Dingxi Academy of Agricultural Sciences,Dingxi 743ooo, China; 3.Horticultural Station of Anding District, Dingxi City, Dingxi 743000, China)

Abstract:Toaccurately identifyand select superior celery varietiessuitableforcultivation in Dingxi city,this paper comprehensively evaluated 2O introduced varieties based on six indicators using principal component analysisand cluster analysis methods.Theresults showed that through thequantitative function model of two principal component factors constructedby principal component analysis,the varieties‘Curacao’‘French

Gren Sword’‘Venus’‘King Celery’,and‘American King’had the highest comprehensive trait scores.By using the heat map cluster dendrogram to systematically classifydiffrentcelery varieties,the2O varieties were clustered into three major categories based on six trait indicators.Among them,the varieties‘Curacao’, ‘French Green Sword',and‘Venus’were clustered into one category,characterized by large single fruit weight,high yield,excellnt quality,and good comprehensive traits.Thisresult was basically consistent with thefieldgrowth situation,indicating the high reliability of the comprehensiveevaluation results.The comprehensive analysis indicated that the varieties‘Curacao’‘French Green Sword’,and‘Venus’ were suitable for cultivation in Dingxi cityand were varieties with market potential that were suitable for promotion.

Keywords: Celery; principal component analysis; cluster analysis; comprehensive evaluation

西芹是芹菜(Apium graveolensL.)的一個種類,為傘形科芹屬一、二年生草本植物,是一種重要的葉菜類蔬菜作物I。西芹廣泛種植于世界各地,目前在我國南北地區(qū)種植面積較大,是人們青睞的蔬菜之一。研究顯示,西芹中含有多種人體不可或缺的營養(yǎng)物質(zhì),如微量元素Zn?Fe?Mn 和 Cu 等含量極高, Ca.Fe 等含量高于白菜、黃瓜、番茄等蔬菜,還含有芹菜素等功能性成分,有調(diào)經(jīng)消炎、降壓、鎮(zhèn)靜和清熱止咳,以及抑制皮膚黑色素瘤的生長、遷移和促進血管生成等功效。

近年來,定西地區(qū)西芹栽培面積不斷擴大,據(jù)統(tǒng)計,定西地區(qū)西芹種植面積2015年為 2 266.7hm2,2024 年為 7 333.3hm2 ,增加了3.24倍。西芹貯運設(shè)施逐步更新,市場營銷體系不斷完善,但西芹的品種工作選育較為滯后,嚴重制約著芹菜產(chǎn)業(yè)的發(fā)展4,加快高產(chǎn)、穩(wěn)產(chǎn)、優(yōu)質(zhì)、抗逆等具有優(yōu)良性狀的西芹新品種的引進篩選,仍然是解決當(dāng)?shù)匚髑燮贩N多而雜、產(chǎn)量低而不穩(wěn)、品質(zhì)不過關(guān)、主栽品種不突出等重要問題的措施7。

目前對西芹各農(nóng)藝性狀的相關(guān)性、主成分分析、多個品種聚類分析的研究還比較缺乏。近年來,品種篩選研究多以品質(zhì)、抗逆性和產(chǎn)量指標(biāo)作為評價因子[8-9],此類分析方法不能完全綜合評價品種的所有特性,而對品種綜合性狀的疊加評價,涉及多個有關(guān)品質(zhì)、產(chǎn)量性狀指標(biāo),進而增加了分析的復(fù)雜性,因此需要一種能更好地進行多性狀綜合選擇的方法。主成分分析和聚類分析方法是采用降維的分析方法,通過數(shù)學(xué)運算剔除次要內(nèi)容保留重要信息,用少量綜合指標(biāo)代替原來多個指標(biāo)絕大部分信息III-I2I,,而聚類分析是對給定對象應(yīng)用數(shù)學(xué)方法進行研究和處理的分類方法,該綜合評價方法已用于分析辣椒[9.14]、番茄[15]、蜜橘[大蔥[17]等果蔬產(chǎn)品的多樣性、多指標(biāo)的品質(zhì)分析,并已取得很好效果,然而采用該方法對西芹品質(zhì)、產(chǎn)量的綜合評價報道較少。本文借助相關(guān)性分析、主成分分析與聚類分析等多元統(tǒng)計分析,對不同西芹品種進行綜合評價,以期構(gòu)建一種有效的西芹引種篩選體系,為引進篩選適合當(dāng)?shù)胤N植的優(yōu)良品種提供理論依據(jù)。

1材料與方法

1.1 試驗材料

參試品種有20個,品種代號為 1~20 ,分別為‘美奐‘埃菲爾'‘法國綠劍'‘維納斯'‘庫拉索'‘高優(yōu)他'‘文圖拉'‘西芹馬倫戈’‘美皇’‘圣地亞哥’‘凰皇’‘英皇’‘精選加州王‘高優(yōu)它‘鳳凰‘美國嫩脆王‘優(yōu)良加州王‘精選四季”“凱爾文’‘蘭巴達’,均由蘭州市西關(guān)農(nóng)貿(mào)市場引進。

試驗用肥料:生物有機肥(有機質(zhì) ≥40% N+P2O+ K20?5% ,黑金氮磷鉀復(fù)合肥 (N-P2O5-K2O=17-17-17) 。

1.2試驗方法

試驗在甘肅省定西市鳳翔鎮(zhèn) (104.5933%.35.56969%) 進行,當(dāng)?shù)睾0?1 898.7m ,年平均氣溫 6.39C ,平均年降水量 400mm ,全年無霜期約 140d; ,屬中溫帶干旱半干旱氣候。試驗設(shè)3次重復(fù),隨機區(qū)組排列,每小區(qū)面積 3m× 5m ,平地覆白膜栽培,株行距為 25cm×30cm ,每 667m2 株數(shù)為8893株。

2022年3月1日在日光溫室中播種育苗,播種前用55°C 溫水處理種子后催芽, 70% 的種子露白后,每穴點播3粒,兩葉一心時間苗,每穴留1株,5月10日移栽到大田,7月底開始陸續(xù)收獲。栽前每 667m2 施生物有機肥 320kg (有機質(zhì) ≥40% ,黑金氮磷鉀復(fù)合肥 40kg ,田間管理同常規(guī)管理。試驗區(qū)土壤為中黏土,土壤初始 pH 為 8,0~20cm 耕層土壤有機質(zhì) 11.2g/kg. 全氮 1.3g/kg. 堿解氮 70.5mg/kg 速效磷 47.4mg/kg 速效鉀 88.3mg/kg 試驗地前茬種植玉米。

1.3 測定項目

西芹成熟期,根據(jù)S取樣法,各小區(qū)選取10株植株,詳細調(diào)查株高 (X1) 、葉柄長 (X2) 、葉柄寬 (X3) 、葉片數(shù)(X4) 、單株質(zhì)量 (X5) 、每 667m2 產(chǎn)量 (X6) 、顏色、纖維、中空程度、燒心情況。測定小區(qū)內(nèi) 1m2 產(chǎn)量,重復(fù)3次,再折合每 667m2 產(chǎn)量。

1.4數(shù)據(jù)分析

采用 Excel2020 進行數(shù)據(jù)整理,利用SPSS25.0對所測數(shù)據(jù)進行統(tǒng)計分析,采用降維-因子分析法、系統(tǒng)聚類-平方歐氏距離法對各性狀指標(biāo)進行主成分分析及聚類分析9-,,并用Origin2021進行熱圖制作。

2 結(jié)果與分析

2.1參試西芹品種植株主要性狀指標(biāo)

從表1可以看出,不同品種西芹的株高差異較大,為62.6~79.3cm ,其中‘法國綠劍'的株高最高,與‘庫拉索'“蘭巴達’‘優(yōu)良加州主'‘維納斯'‘美奐'和‘凱爾文'差異不顯著;‘精選加州王'的株高最低,顯著低于其他品種。不同品種西芹的葉柄長為 29.7~45.1cm ,其中‘庫拉索'最高,與‘法國綠劍'‘維納斯'差異不顯著,顯著高于其他處理;‘精選四季'的葉柄長最低為 29.7cm ,與‘精選加州王'差異不顯著,顯著低于其他品種。不同品種西芹的葉柄寬差異性也較大,為 17.9~30.0mm ,其中‘法國綠劍最高,‘美奐'最低。不同品種西芹的葉片數(shù)不同,其中‘維納斯'最高(18.1片),‘圣地亞哥'最低(10.8片)。不同品種西芹的顏色為黃綠、淺綠和深綠。

表1參試西芹品種的植株性狀指標(biāo)Table1Morphologicalcharacteristicsofcelerytested
注:同列不同小寫字母表示差異顯著( Plt;0.05 ,表2同。

表2參試西芹品種的產(chǎn)量及品質(zhì)

Table2Yieldandqualityofcelerytested

2.2 參試西芹產(chǎn)量及品質(zhì)比較

從表2可以看出,不同品種西芹的單株質(zhì)量為0.80~1.33kg ,‘庫拉索’‘美皇’‘維納斯'和‘法國綠劍’顯著高于其他品種;‘優(yōu)良加州王'最小,與‘凰皇’‘凱爾文‘美奐’‘圣地亞哥'和‘美國嫩脆王'的單株質(zhì)量無顯著性差異,與其他品種差異顯著。不同品種西芹的產(chǎn)量差異較顯著,每 667m2 產(chǎn)量最高的是‘英皇' (9490.61kg) 與‘庫拉索’‘維納斯’法國綠劍'差異不顯著,與其他品種差異顯著;‘庫拉索’‘維納斯’‘法國綠劍’的產(chǎn)量居第二、三、四位,分別為 9462.2,9301.0.9035.3kg/667m2, ‘凰皇'產(chǎn)量最低,每 667m2 產(chǎn)量僅為 5 698.6kg 。

不同品種西芹的纖維素含量,除‘美奐’和‘優(yōu)良加州王'較多外,其他品種表現(xiàn)為較少或中等。中空程度方面,除‘美奐'最多外,其他品種較少或無。不同品種西芹的燒心情況,以‘美奐’最多,‘埃菲爾’‘法國綠劍’‘西芹馬倫戈’‘鳳凰’‘美國嫩脆王'較少,其他品種無燒Table4Principa心情況。

2.3主要性狀指標(biāo)的主成分分析

為了更加客觀地描述研究對象,去除冗余信息,對20個西芹品種的6個性狀指標(biāo)進行Bartlett球形檢驗,結(jié)果見表3。由表3可知,卡方值為 181.047,P 值為0.001。當(dāng)KMO取值大于0.5,且Sig值小于0.05,說明各指標(biāo)具有較強的關(guān)聯(lián)性,數(shù)據(jù)可適用主成分分析方法處理。

表3KMO和巴特利特檢驗Table3KMOandbartletttest

主成分中不同性狀指標(biāo)影響力的高低與方向可用特征向量值的大小來評價,本文采用DPS降維處理對性狀指標(biāo)進行主成分分析,結(jié)果見表4。由表4可知,依據(jù)其特征值 gt;1.0 的提取標(biāo)準選取了前2個主成分。主成分1的貢獻率為 62.661% ,載荷量從大到小依次為每 667m2 產(chǎn)量、單株質(zhì)量、葉片數(shù);主成分2的貢獻率為 21.166% ,載荷量較大的為葉柄長、葉柄寬、株高,前2個主成分的累計貢獻率為 83.827% ,涵蓋了西芹6個性狀指標(biāo)絕大部分的信息,說明前兩個主成分的載荷量越大,西芹在當(dāng)?shù)卦耘噙m應(yīng)性越強。因此,影響西芹栽培適應(yīng)性的主要因素依次是每 667m2 產(chǎn)量、單株質(zhì)量、葉片數(shù)、葉柄長、葉柄寬和株高,則前2個主成分可以作為定西市西芹品種栽培適應(yīng)性評價的綜合指標(biāo)。

表4參試西芹品種的各性狀主成分系數(shù)矩陣

根據(jù)6個性狀2個主成分因子得分結(jié)果,構(gòu)建主成分 PC1~PC2 的因子關(guān)系式和綜合評價函數(shù)模型 F

F1=-0.002X1-0.031X2-0.174X3+0.229X4+0.251X5+ 0.25X6 (20

F2=0.353X1-0.011X2+0.028X3-0.167X4+0.24X5+ 0.781X6+0.248X7+0.772X8

F=λ1/(λ12)×F1+λ2/(λ12)×F2 ,得到評價函數(shù)式 F= 0.748F1+0.252F2 。其中, λ1 和 λ2 分別為 PC1~PC3 對應(yīng)的特征值,綜合得分見表5。位居前五位的品種依次為‘庫拉索’‘法國綠劍'‘維納斯’英皇’‘美皇’。結(jié)合農(nóng)藝性狀發(fā)現(xiàn),得分越高,說明品種綜合性狀越突出,各因子協(xié)調(diào)程度越高。

表5參試西芹品種的綜合得分

Table5 Comprehensive score of celery tested

2.4不同西芹品種聚類分析

依據(jù)主成分分析中6個性狀指標(biāo)結(jié)果,對20個西芹品種進行系統(tǒng)分類,聚類熱圖顏色變化反映數(shù)據(jù)的大小及差異(如圖1)。20個品種依據(jù)6個性狀指標(biāo)聚成三大類,1類包括‘庫拉索’‘法國綠劍'和‘維納斯'3個品種,表現(xiàn)為每 667m2 的產(chǎn)量和株高較高,葉柄長較長,品質(zhì)較好,綜合性狀好,顯著優(yōu)于其他類別。Ⅱ類包括‘英皇‘美皇’‘高優(yōu)他’‘鳳凰’‘文圖拉’‘精選四季’‘西芹馬倫戈和‘精選加州王'這8個品種,表現(xiàn)為單果質(zhì)量較高,其他各指標(biāo)中等;Ⅱ類包括‘蘭巴達‘埃菲爾’‘圣地亞哥‘美奐‘優(yōu)良加州王‘凱爾文’‘凰皇’‘美國嫩脆王’和‘津美西芹'這9個品種,表現(xiàn)為品種株高中等,產(chǎn)量較低,綜合性狀低于其他類別。結(jié)果證明,聚類分析和主成分分析結(jié)果基本一致,主成分分析得分相近的品種基本歸為同類。

圖1參試西芹品種不同性狀的聚類分析Fig.1Clusteranalysisofdifferenttraitsofcelerytested

3討論與結(jié)論

基于主成分分析、聚類分析、相關(guān)性分析的農(nóng)藝性狀鑒定是評價黃瓜[18、番茄[5]、辣椒[9.4]、大蔥[7等品種最基礎(chǔ)的方法,這種方法可以簡單、直觀、全面地分析各蔬菜品種特性及在種植地的適應(yīng)性,也是鑒定和描述種質(zhì)資源研究最基本的方法和途徑[19-20],為目標(biāo)育種和品種篩選提供參考依據(jù)。

本研究對20個西芹品種6個性狀指標(biāo)進行了主成分分析和聚類分析。主成分分析評價結(jié)果表明,得分前5位的品種依次是‘庫拉索’‘法國綠劍’‘維納斯’‘英皇’“美皇’,這5個品種綜合性狀突出,產(chǎn)量均較高。經(jīng)過對6個性狀指標(biāo)繪制聚類分析,將20品種分為3類,其中I類包含的3個品種‘庫拉索’‘法國綠劍'和‘維納斯’,表現(xiàn)為每 667m2 的產(chǎn)量和株高較高,葉柄較長,品質(zhì)也優(yōu),綜合性狀好。該結(jié)果與田間生長情況基本一致,說明綜合評價結(jié)果的可靠性較高;Ⅱ類包括‘英皇’‘美皇’“高優(yōu)他’‘鳳凰’‘文圖拉’‘精選四季'‘西芹馬倫戈'和精選加州王'8個品種,表現(xiàn)為單果質(zhì)量較高,其他各指標(biāo)中等;Ⅲ類包括‘蘭巴達'‘埃菲爾‘圣地亞哥’‘美奐’‘優(yōu)良加州王'‘凱爾文'‘凰皇’‘美國嫩脆王'和‘津美西芹'9個品種,表現(xiàn)為產(chǎn)量較低,綜合性狀較差,經(jīng)濟價值不高,不建議在本地推廣種植。綜合分析表明,‘庫拉索’法國綠劍'和‘維納斯'共3個品種表現(xiàn)為綜合性狀優(yōu),產(chǎn)量高,綜合得分高,說明各因子之間的協(xié)調(diào)性較好,適合在定西地區(qū)進行栽培,是適宜推廣并具有應(yīng)用潛力的品種。

參考文獻:

[1]李夢瑤,王楓,侯喜林,等.芹菜過敏原蛋白 Apigl 基因的克隆與表達分析[J].南京農(nóng)業(yè)大學(xué)學(xué)報,2013,36(2):13-19.

[2]譚國‘.斤米化目系和斤米系的代謝機理礦九[D].南京:南京農(nóng)業(yè)大學(xué),2017:18-26.

[3]熊紹風(fēng),李煒,邱雪峰,等.芹菜素改善高糖高脂應(yīng)激所致血管內(nèi)皮細胞功能失調(diào)的機制研究[J].南昌大學(xué)學(xué)報(醫(yī)學(xué)版),2024,64(6): 27-32,72.

[4]黃凱美,趙捷,顏韶兵,等.杭州秋冬季芹菜品種篩選試驗[J].浙江農(nóng)業(yè)科學(xué),2019,60(9):1561-1562,1702.

[5]李冬,代國鵬,閆菊紅,等.覆膜壓沙芹菜新品種引進比較研究[J].寧夏農(nóng)林科技,2013,54(12):36-37.

[6]司魯俊.北方露地西芹綠色標(biāo)準化栽培技術(shù)[J].農(nóng)業(yè)科技通訊,2024(7): 195-197.

[7]王姣敏,呂海龍,魏鏞頻,等.氮肥減量與生物菌肥配施對連作西芹土壤性質(zhì)及西芹生長的影響[J].土壤與作物,2022,11(2): 217-225.

[8]賈麗麗,劉惠吉,王華,等.3種不同芹菜營養(yǎng)品質(zhì)的比較與評價[J].江蘇農(nóng)業(yè)科學(xué),2021,49(2):146-149.

[9]賈麗麗,沈迪,陳龍正,等.蘇南地區(qū)秋冬茬移栽芹菜農(nóng)藝性狀的比較與評價[J].上海農(nóng)業(yè)學(xué)報,2021,37(1):1-7.

[10] 高云,郁志芳.基于主成分分析的芹菜品質(zhì)評價[J].食品工業(yè)科技,2020,41(3):308-314,320.

[11] 蔡曉洋,張思荻,曾俊,等.基于主成分分析和聚類分析的梔子種質(zhì)資源評價[J].中國實驗方劑學(xué)雜志,2017,23(14):30-37.

[12]李潔.芹菜種質(zhì)資源的耐熱性評價及高溫應(yīng)答機制初步研究[D].雅安:四川農(nóng)業(yè)大學(xué),2022:32-36.

[13] 許麗利.類分析的算法及應(yīng)用[D]長春:吉林大學(xué),2010:15-19.

[14] 李燕,侯大山,張慶銀,等.基于主成分分析和聚類分析對11個辣椒品種產(chǎn)量和品質(zhì)綜合評價及引進篩選[J].河北科技師范學(xué)院學(xué)報,2022,36(3): 22-28.

[15] 韓澤群,姜波.加工番茄品種多性狀綜合評價方法研究[J].中國農(nóng)業(yè)科學(xué),2014,47(2):357-365.

[16] 徐回林,陳金印,許熙帥,等.利用主成分分析和聚類分析綜合評價南豐蜜橘果實品質(zhì)[C]//中國園藝學(xué)會.設(shè)施園藝與園藝作物標(biāo)準化生產(chǎn)技術(shù)交流會論文匯編.南昌:江西農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,2012:107-112.

[17] 苗錦山,劉彩霞,戴振建,等.蔥種質(zhì)資源數(shù)量性狀的聚類分析、相關(guān)性和主成分分析[J].中國農(nóng)業(yè)大學(xué)學(xué)報,2010,15(3): 41-49.

[18] 陳林,曹振強,梁肇均,等.黃瓜種質(zhì)資源農(nóng)藝性狀多樣性分析[J].廣東農(nóng)業(yè)科學(xué),2021,48(6):15-22.

[19] 張國儒,唐亞萍,楊濤,等.272份辣椒資源的農(nóng)藝性狀分析[J].新疆農(nóng)業(yè)科學(xué),2021,58(12):2300-2311.

[20]SINGHBK,SINGHB,SINGHPM.Breeding cauliflower:Areview[J]. International Journal of Vegetable Science,2018,24(1): 58-84.

猜你喜歡
西芹性狀品種
桉樹幼林抗臺風(fēng)能力與提升策略
紅軍東方軍與十九路軍王臺談判始末
黨史博覽(2025年8期)2025-08-20 00:00:00
籽粒玉米品種農(nóng)藝性狀綜合評價及品質(zhì)性狀鑒定
郁金香表型性狀的遺傳變異及多樣性分析
芹菜:餐桌上的“寶藏食材”
中國食品(2025年13期)2025-08-17 00:00:00
15個鮮食棗品種在晉中地區(qū)的引種表現(xiàn)
山西果樹(2025年4期)2025-08-15 00:00:00
再生種植圖鑒
鳳凰生活(2025年8期)2025-08-10 00:00:00
2023-2024年度全國小麥品種推廣應(yīng)用情況
雪牛雙味
美食(2025年7期)2025-07-26 00:00:00
為小麥大面積單產(chǎn)提升遴選新的當(dāng)家品種
種子世界(2025年7期)2025-07-26 00:00:00
主站蜘蛛池模板: 国产在线无码av完整版在线观看| 亚洲高清中文字幕在线看不卡| 在线观看无码a∨| 国产乱码精品一区二区三区中文 | 国产精品白浆无码流出在线看| 亚洲一区网站| 久久综合九色综合97网| 免费国产不卡午夜福在线观看| 青青热久麻豆精品视频在线观看| 国产91丝袜在线播放动漫| 国内精自视频品线一二区| 亚洲日本一本dvd高清| 无码AV动漫| 午夜少妇精品视频小电影| 欧美中文字幕无线码视频| 久久精品aⅴ无码中文字幕| 欧美精品高清| 五月婷婷伊人网| 大陆精大陆国产国语精品1024 | 国语少妇高潮| 亚洲第一国产综合| 国产成人夜色91| 国产一在线观看| 欧美另类视频一区二区三区| 狂欢视频在线观看不卡| 四虎影视无码永久免费观看| 成人午夜精品一级毛片| 欧美乱妇高清无乱码免费| 2021最新国产精品网站| 亚洲精选高清无码| 97国内精品久久久久不卡| 免费国产一级 片内射老| 无遮挡一级毛片呦女视频| 91精品国产丝袜| 国产精品99久久久久久董美香| 97se亚洲综合| 亚洲国产综合精品中文第一| 国产成人亚洲日韩欧美电影| 久久精品日日躁夜夜躁欧美| 免费久久一级欧美特大黄| 一级黄色网站在线免费看| 欧美国产日韩在线播放| 91精品免费久久久| 久久人妻xunleige无码| 不卡视频国产| 日韩av无码精品专区| 欧美三级视频在线播放| 欧美特级AAAAAA视频免费观看| 国产微拍一区二区三区四区| 国产在线观看人成激情视频| 亚洲欧美h| 亚洲男人的天堂在线观看| 性色生活片在线观看| yjizz国产在线视频网| 人妻丰满熟妇αv无码| 欧美特黄一免在线观看| 国产成人欧美| 毛片手机在线看| 中文字幕伦视频| 国产日本一区二区三区| 制服丝袜亚洲| 国产成人无码AV在线播放动漫 | 国产噜噜在线视频观看| 免费aa毛片| 国产迷奸在线看| 亚洲成A人V欧美综合| 亚洲第一视频网站| 精品视频91| A级毛片高清免费视频就| 国产一二视频| 91系列在线观看| 99精品免费在线| 国产欧美日韩18| 欧美精品不卡| 日本www色视频| AV天堂资源福利在线观看| 在线观看国产精美视频| 一级看片免费视频| 亚洲AV永久无码精品古装片| 国产成人免费视频精品一区二区| 亚洲成年人网| 国产精品亚洲天堂|