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
主站蜘蛛池模板: 国产香蕉一区二区在线网站| 毛片视频网址| 国产欧美又粗又猛又爽老| 大陆精大陆国产国语精品1024 | 亚洲成aⅴ人在线观看| 91精品小视频| 免费99精品国产自在现线| 色综合国产| 亚洲av日韩综合一区尤物| 国产免费人成视频网| 伊人久久婷婷五月综合97色| 亚洲综合片| 国产成人91精品| 亚洲精品中文字幕无乱码| 欧美在线一二区| 欧美综合区自拍亚洲综合绿色| 国产在线一区二区视频| 毛片三级在线观看| 亚洲第一综合天堂另类专| 一级毛片免费播放视频| 欧美成人手机在线视频| 日本一本在线视频| 91口爆吞精国产对白第三集| 欧美黄网在线| 国产美女视频黄a视频全免费网站| 一级毛片免费不卡在线视频| 久久精品66| 久久久国产精品免费视频| 国产毛片高清一级国语 | 久久综合色天堂av| 高清无码一本到东京热| 久热中文字幕在线| 欧美成人日韩| 亚洲热线99精品视频| 久久综合五月婷婷| 日韩黄色在线| 亚洲全网成人资源在线观看| 伊人激情综合| 日韩在线欧美在线| 波多野结衣一区二区三区88| 国产一区二区三区在线观看视频| 黄色网址免费在线| 露脸真实国语乱在线观看| 久久国产高潮流白浆免费观看| 特级做a爰片毛片免费69| 99视频在线精品免费观看6| 亚洲91精品视频| aⅴ免费在线观看| 国产91小视频在线观看| 久久精品aⅴ无码中文字幕 | 丁香婷婷在线视频| 国产乱子伦一区二区=| 91偷拍一区| 国产精品3p视频| 国产美女无遮挡免费视频| 国产欧美在线观看精品一区污| 亚洲婷婷在线视频| 国产男女免费完整版视频| 欧美在线精品怡红院| 正在播放久久| 日韩精品毛片人妻AV不卡| 久久久久久久久亚洲精品| 国产精品偷伦在线观看| 在线色国产| 亚洲欧洲日韩综合色天使| 伊人成人在线| 国产成人禁片在线观看| 亚洲精品高清视频| 91精品情国产情侣高潮对白蜜| 欧美中文一区| 欧美另类第一页| 国产成人h在线观看网站站| 日本免费高清一区| 久久精品无码中文字幕| 亚洲中文字幕在线一区播放| 亚洲精品视频免费| 一级在线毛片| 亚洲天堂网2014| 9966国产精品视频| 日韩国产亚洲一区二区在线观看| 亚洲精品成人福利在线电影| 亚洲日韩AV无码一区二区三区人|