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

    Detection of Broad Bean Diseases by Fourier Transform Infrared Spectroscopy Combined with Curve Fitting

    2015-02-05 03:30:53XiaohuaWANGGangLIUQuanhongOUXiangpingZHOUJianmingHAOJianhongLIULuxiangWANG
    Agricultural Science & Technology 2015年6期
    關(guān)鍵詞:劉剛曲線擬合傅里葉

    Xiaohua WANG,Gang LIU*,Quanhong OU,Xiangping ZHOU,Jianming HAO,Jianhong LIU,Luxiang WANG

    1.School of Physics and Electronic Information,Yunnan Normal University,Kunming 650500,China;

    2.Key Laboratory of Advanced Technique&Preparation for Renewable Energy Materials,Ministry of Education,Yunnan Normal University,Kunming 650500,China;

    3.Institute of Quality Standards and Testing Technology,Yunnan Academy of Agricultural Sciences,Kunming 650223,China

    Detection of Broad Bean Diseases by Fourier Transform Infrared Spectroscopy Combined with Curve Fitting

    Xiaohua WANG1,Gang LIU1*,Quanhong OU1,Xiangping ZHOU1,Jianming HAO1,Jianhong LIU2,Luxiang WANG3

    1.School of Physics and Electronic Information,Yunnan Normal University,Kunming 650500,China;

    2.Key Laboratory of Advanced Technique&Preparation for Renewable Energy Materials,Ministry of Education,Yunnan Normal University,Kunming 650500,China;

    3.Institute of Quality Standards and Testing Technology,Yunnan Academy of Agricultural Sciences,Kunming 650223,China

    Fourier transform infrared(FTIR)spectroscopy was used to study diseased leaves in broad bean.Results showed that the infrared spectra of different broad bean diseased leaves were similar,which were mainly made up of the vibrational absorption bands of protein,lipid and polysaccharide.There were minor differences including the spectral peak position,peak shape and the absorption intensity in the range of 1 800-1 300 cm-1.There were obvious differences among their second derivative spectra in the range of 1 800-1 300 cm-1.After the procedure of the Fourier self-deconvolution and curve fitting of health bean leaves and broad bean diseased leaves in the range of 1 700-1 500 cm-1,three sub-peaks were obtained at 1 550 cm-1(protein amideⅡband),1 605 cm-1(lignin)and 1 650 cm-1(protein amide I band).The ratios of relative areas of the bands of amideⅡ,lignin,and amide I were 38.86%,28.68%and 32.47%in the spectra of healthy leaves,respectively.It was distinguished from the diseased leaves(chocolate spot leaf:15.42%, 42.98%and 41.61%,ring spot leaf:32.39%,35.63%and 31.98%,rust leaf:13.97%, 46.40%and 39.65%,yellowing leaf curl disease leaf:24.01%,36.55%and 39.44%). For sub-peak area ratios(A1563/A1605,A1650/A1605and A1563/A1654),those of four kinds of diseased leaves were smaller than that of healthy leaves,and there were also differences among four kinds of diseased leaves.The results proved that FTIR combining with curve fitting might be a potentially useful tool for detecting different kinds of broad bean diseases.

    Fourier transform infrared(FTIR)spectroscopy;Broad bean diseases; Second derivative spectra;Curve fitting

    B road bean is annual or biennial herb,can be used as grain, vegetable,fodder and green manure,and derives from Southwest Asia and North Africa.Broad bean contains eight kinds of essential amino acids,and carbohydrate content is 47%-60%.It has rich nutritive value. Broad bean contains abundant calcium,protein and phospholipid,and does not contain cholesterol,which is beneficial to human health.Tender broad bean has the effect neutralizing stomach and moistening intestinal tract;broad bean stem has the function of hematischesis and arresting diarrhea;broad bean flower has the effect of depressurization,arresting leucorrhoea and hematischesis[1].

    Broad bean yield and quality are affected by plant diseases and insect pests,and plant disease is dominated by fungus disease and virus.Leaf diseases are mainly various leaf spot diseases,such as rust disease,ring spot disease,chocolate spot disease, brown spot disease,yellowing leaf curl disease(virus).The variety is various,

    Common methods identifying crop disease include traditional biological method,microstructure observation and molecular biology method. These methods all require rich experience,skilled operation ability and long time.Fourier transform infrared(FTIR) spectroscopy is a kind of structure analysis technique based on vibrations of functional group and polar bond, can reflect vibrational mode of molecular functional group,has fingerprint characteristic,and has become forceful means detecting macromolecular structure and interaction in biological tissue.Infrared spectroscopy does not need separation and extraction process of sample,has the characteristics of simple operation,fast,low cost and non-pollution,and has been widely applied in studying plant and medicine[4-8]. Li Zhiyong et al.used FTIR to analyze the difference between diseased leaves(broad bean rust disease,stem base rot disease,ring spot disease and yellowing leaf curl disease)and healthy leaves,and used principal component analysis and clustering analysis to classify disease type of sample,which obtained better effect[9]. The purpose of this study is to discriminate different diseased leaves of broad bean by FTIR spectroscopy combined with curve fitting.

    Materials and Methods

    Materials

    Broad bean leaf samples were collected from Luliang County,Qujing City,Yunnan Province during February-March of 2011 and 2012.After natural drying,sample was conserved to measure.When sampling,all samples avoided the vein,and we took the leaves at diseased spot.Frontier FTIR spectrometer from Perkin Elmer Company was used to measure spectrum.Spectral range was 4 000-400 cm-1;spectral resolution was 4 cm-1, and scanning times was 16.

    Spectrum collection and data processing

    Little sample and moderate KBr were put into agate mortar to grind. After grinding evenly,thin disc was pressed.Infrared spectrum was measured,and sample spectrum automatically deducted KBr background.Each sample collected five infrared absorption spectra,and the mean of five spectra was taken as infrared absorption spectrum of the sample.The spectrum was conducted by nine-point smoothing of OMNIC6.0 infrared spectrum software,base line correction and normalized preprocessing.In the range of 1 800-1 300 cm-1,software Origin 8.5 was used to obtain second-order derivative spectrum. Software Origin 8.5 was used to conduct curve fitting on original spectrum.

    Results and Analysis

    Infrared spectral characteristic analysis of broad bean leaves

    Average infrared spectra of broad bean leaves were shown as Fig.1. Seen from Fig.1,major peaks of FTIR spectra for broad bean leaves included:there was a strong and wide absorption peak at 3 410 cm-1,which was mainly stretching vibration absorption of O-H[10];absorption peaks at 2 920 and 2 850 cm-1belonged to anti-symmetric and symmetric stretching vibrations of CH2[11];absorption peak at 2 954 cm-1belonged to anti-symmetric stretching vibration of-(CH3),which was mainly from tissue components of cell wall,such as protein,carbohydrate and lipid;absorption peak at 1 737 cm-1belonged to vibration absorption of lipid carboxyl group;absorption peak at 1 650 cm-1belonged to stretching vibration absorption of protein amide I band C=O,while absorption peak at 1 548 cm-1belonged to bending vibration of amide I band N=H and stretching vibration absorption of C-N;it was mixed vibration absorption region of protein,fat,lignin and polysaccharide between 1 500 and 1 200 cm-1[12],and peaks at 1 456 and 1 380 cm-1belonged to anti-symmetric and symmetric variable angle vibration absorption of methyl;absorption peak at 1 265 cm-1belonged to vibration absorption of protein amide III band[13];peak at 1 152 cm-1belonged to C-OH stretching vibration of polysaccharide, while peak near 1 075 cm-1belonged to C-O stretching vibration of polysaccharide.These spectral features reflected that major components of broad bean leaves were protein, polysaccharide and lipid.

    Seen from Fig.1,infrared spectra of broad bean leaves infected by different diseases were roughly similar, and there was small difference in the range of 1 800-700 cm-1.Second derivative spectrum could amplify difference and improve spectral resolution.Fig.2 was second derivative spectra of different diseased leaves.In the 1 800-1 300 cm-1region,spectra of different diseased leaves were contrasted.It was clear that two strong peaks of leaves infected by broad bean rust disease(d)were respectively near 1 575 and 1 380 cm-1, and a strong peak of leaves infected by broad bean chocolate spot disease (b)was near 1 575 cm-1.Strong peak of healthy broad bean leaves(e)was near 1 470 cm-1,while other two kinds of leaves had a moderate-intensity absorption peak near 1 470 cm-1.In addition,peak intensity of broad bean leaves infected by rust disease(d)at 1 575 cm-1was stronger than that at 1 380 cm-1.Healthy broad bean leaves (e)had a strong peak near 1 540 cm-1, while other four kinds of leaves all presented very weak peaks at 1 530 cm-1. By comparing ring spot leaf(c)with yellowing leaf curl disease leaf(a),it was clear that absorption peak of broad bean leaves infected by ring spot disease(c)was stronger near 1 460 cm-1,while absorption peaks of broad bean leaves infected by yellowing leaf curl disease were stronger near 1 390 and 1 440 cm-1.It was clear that second derivative spectrum could distinguish broad bean leaves infected by different diseases.

    Curve fitting

    Intensity of infrared spectral absorption peak is often determined by peak height,peak width and peak intensity,and it is not enough accurate to use peak height to describe intensity.After Fourier self-deconvolution and curve fitting sub-peak processing,it is more close to actual situation to use peak area to represent its intensity[14-16]. Seen from Fig.2,the difference between healthy leaves and diseased leaves mainly concentrated between 1 700 and 1 500 cm-1.Curve fitting of infrared spectrum between 1 700-1 500 cm-1for different broad bean diseased leaves was conducted,and three sub-peaks were obtained(Fig.3). Absorption peak at 1 650 cm-1was mainly absorption of protein amide I band;absorption peak near 1 550 cm-1was mainly absorption of protein amide II band;it was absorption peak of lignin near 1 605 cm-1.Via area ratio of absorption peak(AamideI/Alignin, AamideII/Aligninand AamideII/AamideI),relative content change of protein in the leaves could be judged.

    Table 1Sub-peak area ratios of diseased leaves and healthy leaves for broad bean

    After Fourier self-deconvolution and curve fitting peak processing,yellowing leaf curl disease leaves obtained three sub-peaks at 1 559,1 610 and 1 654 cm-1.It was defined that anysub-peak area divided by area sum of three sub-peaks was area proportion of the sub-peak.Corresponding subpeak area proportions of yellowing leaf curl disease leaf were respectively 24.01%,36.55%and 39.44%;chocolate spot disease leaf obtained three sub-peaks at 1 547,1 611 and 1 656 cm-1,and corresponding sub-peak area proportions were respectively 15.42%,42.98%and 41.61%;ring spot leaf disease leaf obtained three sub-peaks at 1 565,1 611 and 1 654 cm-1,and corresponding sub-peak area proportions were respectively 32.39%,35.63%and 31.98%.Rust leaf had three sub-peaks at 1 546, 1 606 and 1 654 cm-1,and corresponding sub-peak area proportions were respectively 13.94%,46.60% and 39.65%;while healthy leaves had three sub-peaks at 1 563,1 605 and 1 654 cm-1,and corresponding sub-peak area proportions were respectively 38.86%,28.68%and 32.47%. At 1 605 cm-1,peak position,type and absorption intensity of each spectrum were all similar.Via area ratios of absorption peak(AamideI/Alignin,AamideII/Aligninand AamideII/AamideI),relative content change of protein in the diseased leaves could be judged.

    Table 1 was peak area ratios of protein amide I band,amide II band and lignin in different broad bean diseased leaves.Seen from Table 1, AamideI/Aligninin healthy leaves was more than 1.1,and others were all less than 1.1,in which AamideI/Aligninof yellowing leaf curl disease leaf was close to 1.1; AamideII/Aligninof yellowing leaf curl disease leaf,chocolate spot disease leaf, ring spot disease leaf and rust disease leaf was all less than 1,while AamideII/Aligninof healthy leaf was more than 1.1;the maximum AamideII/AamideIwas 1.20 in healthy leaf,while the minimum was 0.35 of rust disease leaf.It showed that compared with healthy leaves,protein content in broad bean leaves infected by different diseases relatively declined,which indicated that area ratio difference of sub-peak could be used to distinguish different broad bean diseased leaves.

    Conclusions

    We studied FTIR spectra of four kinds of broad bean diseased leaves and healthy leaves.Results showed that major components of broad bean leaves were protein,lipid and polysaccharide.Their infrared spectra were overall similar at peak shape,intensity and position of footprint characteristic region.In the range of 1 700-1 500 cm-1,after Fourier self-deconvolution and curve fitting peak processing,area ratios of three sub-peaks(A1650/A1605, A1563/A1605and A1563/A1654)had obvious difference.Research showed that FTIR spectroscopy combining with curve fitting is a convenient,fast and nondestructive method for identifying different broad bean diseased leaves.

    [1]ZHENG KB(鄭開(kāi)斌),LI AP(李?lèi)?ài)萍), ZANG CR(臧春榮),et al.Green legume vegetables——Broad bean(綠色豆類(lèi)蔬菜——蠶豆)[J].Xiandai Horticulture(現(xiàn)代園藝),2006(1):38-39.

    [2]WANG XM(王曉鳴),ZHU ZD(朱振東), DUAN CX(段燦星),et al.Identification techniques of plant diseases and insect pests for broad bean and pea(蠶豆豌豆病蟲(chóng)害鑒別技術(shù))[M].Beijing:China Agricultural Science and Technology Press(北京:中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社), 2007.

    [3]OUYANG DS(歐陽(yáng)迪莎).Plant disease management in sustainable agriculture (可持續(xù)農(nóng)業(yè)中的植物病害管理)[D]. Fuzhou:Fujian Agriculture and Forestry University(福州:福建農(nóng)林大學(xué)),2005.

    [4]YAN G,SONG P,SUN SQ,et al.Differentiation and quality estimation of Cordyceps with infrared spectroscopy [J].Spectrochimica Acta Part A:Molecular and Biomolecular Spectroscopy, 2009,74:983-990.

    [5]LIU R,WANG M,DUAN JA,et al.Antipyretic and antioxidant activities of the aqueous extract of Cornububali(Water Buffalo Horn)[J].The American Journal of Chinese Medicine,2010,38:293-306.

    [6]SUN SQ(孫素琴),ZHOU Q(周群), CHEN JB(陳建波).Infrared spectral analysis and identification of traditional Chinese medicine(中藥紅外光譜分析與鑒定)[M].Beijing:Chemical Industry Press(北京:化學(xué)工業(yè)出版社),2010.

    [7]MATA PD,DOMINGUZ VA,SENDRA JMB,et al.Olive oil assessment in edible oil blends by means of ATR-FTIR and chemometrics[J].Food Control, 2012,23:449-455.

    [8]SUN SQ(孫素琴),ZHOU Q(周群),QIN Z(秦竹).Two-dimensional related infrared spectral identification atlas of traditional Chinese medicine(中藥二維相關(guān)紅外光譜鑒定圖集)[M].Beijing: Chemical Industry Press(北京:化學(xué)工業(yè)出版社),2003.

    [9]LI ZY(李志永),LIU G(劉剛),LI L(李倫), et al.FTIR research of diseased broad bean leaves(蠶豆病害葉的FTIR研究) [J].Spectroscopy and Spectral Analysis (光譜學(xué)與光譜分析),2012,32(5): 1217-1220.

    [10]SHI YM(時(shí)有明),LIU G(劉剛),LIU JH (劉劍虹),et al.Study on black fungus from five different production places by FTIR spectroscopic methodology(利用傅里葉變換紅外光譜法研究了五種不同產(chǎn)地的黑木耳)[J].Acta Optica Sinica(光學(xué)學(xué)報(bào)),2007,27(1):129-132.

    [11]DUBRAVKA K,MAJA B,MARIN K. FT-IR spectroscopy of lipoproteins——A comparative study[J]. Spectrochimica Acta Part A:Molecular and Biomolecular Spectroscopy,2009, 73:701-706.

    [12]WENG SP(翁詩(shī)甫).FTIR spectroscope (傅里葉變換紅外光譜儀)[M].Beijing: Chemical Industry Press(北京:化學(xué)工業(yè)出版社),2005.

    [13]LU XN,MOLLY AHW,MARIAH J,et al.A study of biochemical parameters associated with ovarian atresia and quality of caviar in farmed white sturgeon(Acipenser transmontanus)by Fourier Transform Infrared(FT-IR) Spectroscopy[J].Aquaculture,2011, 315:298-305.

    [14]BYLER DM,SUSI H.Examination of the secondary structure of proteins by deconvolved FTIR spectra[J].Biopolymers,1986,25:469-487.

    [15]KORKMAZ F,K?STER S,YILDIZ?, et al.In situ opening/closing of OmpG from E.coli and the splitting of β-sheet signals in ATR-FTIR spectroscopy[J]. Spectrochimica Acta Part A,2012,91 (6):395-401.

    [16]SAMYN P,NIEUWKERKE VD,SCHOUKENS G,et al.Quality and statistical classification of brazilian vegetable oils using mid-infrared and raman spectroscopy[J].Applied Spectros copy,2012,66(5):552.

    Responsible editor:Ping SONG

    Responsible proofreader:Xiaoyan WU

    基于曲線擬合的蠶豆病害葉的FTIR研究

    汪小華1,劉剛1*,歐全宏1,周湘萍1,郝建明1,劉劍虹2,汪祿祥3
    (1.云南師范大學(xué)物理與電子信息學(xué)院,云南昆明650500;2.云南師范大學(xué)可再生能源材料先進(jìn)技術(shù)與制備教育部重點(diǎn)實(shí)驗(yàn)室,云南昆明650500;3.云南省農(nóng)業(yè)科學(xué)院質(zhì)量標(biāo)準(zhǔn)與檢測(cè)技術(shù)研究所,昆明650223)

    用傅里葉變換紅外光譜法研究蠶豆病害葉片,結(jié)果顯示不同病害蠶豆葉片紅外光譜圖整體相似,它們的紅外光譜主要由蛋白質(zhì)、脂類(lèi)和多糖的振動(dòng)吸收帶組成,僅在1 800~1 300 cm-1范圍光譜的峰位、峰形及吸收強(qiáng)度有一些微小差異。對(duì)1 800~1 300 cm-1波數(shù)范圍的光譜圖進(jìn)行二階導(dǎo)數(shù)處理,結(jié)果顯示蠶豆病害葉的二階導(dǎo)數(shù)譜差異明顯。對(duì)健康和病害蠶豆葉1 700~1 500 cm-1范圍光譜進(jìn)行傅里葉自去卷積和曲線擬合處理后,得到蛋白質(zhì)酰胺Ⅱ帶(1 550 cm-1)、木質(zhì)素(1 605 cm-1)和酰胺Ⅰ(1 650 cm-1)3個(gè)子峰,相應(yīng)子峰的峰面積比例顯示差異,黃化卷葉病分別為24.01%、36.55%、39.44%,赤斑病分別為15.42%、42.98%、41.61%,輪紋病分別為32.39%、35.63%、31.98%,銹病分別為13.97%、46.40%、39.65%,健康葉片分別為38.86%、28.68%、32.47%,健康葉的酰胺Ⅱ帶子峰相對(duì)面積比病害葉的大,而其木質(zhì)素子峰相對(duì)面積比病害葉的小。對(duì)于子峰面積比A1563/A1605、A1650/A1605和A1563/A1654,4種病害葉的比值均比健康葉的相應(yīng)數(shù)值小,4種病害葉之間也有差異。結(jié)果表明傅里葉變換紅外光譜(FTⅠR)結(jié)合曲線擬合可望對(duì)不同病害的樣品進(jìn)行有效鑒別。

    傅里葉變換紅外光譜;蠶豆病害;二階導(dǎo)數(shù)光譜;曲線擬合and harm is great.They mainly harm leaf,sometimes stem,leaf handle and pod.If disease occurs in large area,serious economic loss can be caused[2].At present,plant disease prevention and control are dominated by chemical pesticide.Unreasonable utilization of pesticide and chemical fertilizer,drug resistance of harmful organisms and pesticide residue,etc. cause that eco-environment is continuously destroyed,and agricultural product quality declines.So,we should correctly identify disease and take reasonable prevention and control measures for corresponding disease,to reach sustainable development of crop eco-environment[3].

    國(guó)家自然科學(xué)基金項(xiàng)目(30960179);云南省高校科技創(chuàng)新團(tuán)隊(duì)支持計(jì)劃。

    汪小華(1988-),女,安徽宿松人,碩士研究生,研究方向:生物光譜學(xué),E-mail:zhaoxxiang1987@sina.com。*通訊作者,E-mail:gliu66@163.com。

    2015-03-01

    修回日期 2015-05-09

    Supported by National Natural Science Foundation of China(30960179);Program for Innovative Research Team in Science and Technology in University of Yunnan Province.

    *Corresponding author.E-mail:gliu66@163.com

    Received:March 1,2015 Accepted:May 9,2015

    猜你喜歡
    劉剛曲線擬合傅里葉
    綠水青山 朗朗乾坤
    都市(2022年9期)2022-09-07 09:15:12
    全球高通脹和貨幣政策轉(zhuǎn)向
    雙線性傅里葉乘子算子的量化加權(quán)估計(jì)
    基于小波降噪的稀疏傅里葉變換時(shí)延估計(jì)
    曲線擬合的方法
    基于曲線擬合的投棄式剖面儀電感量算法
    Matlab曲線擬合工具箱在地基沉降預(yù)測(cè)模型中的應(yīng)用
    Matlab曲線擬合法在地基沉降預(yù)測(cè)中的應(yīng)用
    基于傅里葉變換的快速TAMVDR算法
    快速離散傅里葉變換算法研究與FPGA實(shí)現(xiàn)
    伊人久久国产一区二区| 国产精品国产三级专区第一集| 国产片特级美女逼逼视频| 99九九线精品视频在线观看视频| 天天躁夜夜躁狠狠久久av| a级一级毛片免费在线观看| 免费观看的影片在线观看| 亚洲一级一片aⅴ在线观看| 免费人妻精品一区二区三区视频| 男女无遮挡免费网站观看| 久久精品国产亚洲网站| 舔av片在线| 九九久久精品国产亚洲av麻豆| 交换朋友夫妻互换小说| 人妻 亚洲 视频| 青春草国产在线视频| 九九爱精品视频在线观看| 色视频www国产| 高清不卡的av网站| 大片电影免费在线观看免费| 免费人妻精品一区二区三区视频| 国产精品.久久久| 中文字幕久久专区| 日韩电影二区| 亚洲精品aⅴ在线观看| 精品国产乱码久久久久久小说| 久久久久久九九精品二区国产| 十八禁网站网址无遮挡 | 成人国产av品久久久| 日韩av在线免费看完整版不卡| 国产白丝娇喘喷水9色精品| 青春草国产在线视频| 亚洲高清免费不卡视频| 欧美国产精品一级二级三级 | 亚洲欧美精品专区久久| 大又大粗又爽又黄少妇毛片口| 免费在线观看成人毛片| 午夜福利视频精品| 天美传媒精品一区二区| 亚洲精品国产成人久久av| 国产有黄有色有爽视频| 亚洲国产日韩一区二区| 亚洲av不卡在线观看| 国内精品宾馆在线| 在现免费观看毛片| 午夜福利视频精品| 亚洲婷婷狠狠爱综合网| 精品一区二区免费观看| 亚洲激情五月婷婷啪啪| 大话2 男鬼变身卡| 国产免费福利视频在线观看| 22中文网久久字幕| 在线观看免费日韩欧美大片 | 国产亚洲91精品色在线| 久久久久久九九精品二区国产| 亚洲av二区三区四区| 国产成人一区二区在线| 97在线人人人人妻| 久久久久久人妻| 少妇人妻精品综合一区二区| 51国产日韩欧美| 女性被躁到高潮视频| 精品久久久久久久末码| 欧美一区二区亚洲| 中国国产av一级| 黑人猛操日本美女一级片| 99九九线精品视频在线观看视频| 亚洲av电影在线观看一区二区三区| 久久久久久伊人网av| 少妇被粗大猛烈的视频| 22中文网久久字幕| 亚洲av二区三区四区| 性色av一级| 国产成人a∨麻豆精品| 亚洲一级一片aⅴ在线观看| 丰满人妻一区二区三区视频av| 我的老师免费观看完整版| av.在线天堂| 内射极品少妇av片p| 久久这里有精品视频免费| 亚洲精品456在线播放app| 亚洲人成网站在线播| 精品人妻一区二区三区麻豆| 免费看av在线观看网站| 免费看日本二区| 如何舔出高潮| 午夜精品国产一区二区电影| 极品少妇高潮喷水抽搐| 丝袜喷水一区| 日韩欧美一区视频在线观看 | av卡一久久| 久久青草综合色| 啦啦啦中文免费视频观看日本| 91精品国产九色| 少妇的逼水好多| 亚洲无线观看免费| 不卡视频在线观看欧美| 日韩视频在线欧美| 中文字幕人妻熟人妻熟丝袜美| 老女人水多毛片| 免费黄频网站在线观看国产| 久久毛片免费看一区二区三区| 成人美女网站在线观看视频| 日韩欧美一区视频在线观看 | a 毛片基地| 建设人人有责人人尽责人人享有的 | 欧美成人精品欧美一级黄| 亚洲国产毛片av蜜桃av| 久久精品国产亚洲av天美| 男人爽女人下面视频在线观看| 国产精品一区www在线观看| 亚洲不卡免费看| 网址你懂的国产日韩在线| 国产 一区精品| 国产真实伦视频高清在线观看| 国内揄拍国产精品人妻在线| 欧美bdsm另类| 99久久精品热视频| 成人国产麻豆网| 亚洲怡红院男人天堂| 国产成人91sexporn| 日韩欧美精品免费久久| 亚洲国产精品国产精品| 99久久精品一区二区三区| 国产成人一区二区在线| 久久av网站| 亚洲最大成人中文| 欧美激情极品国产一区二区三区 | 国产精品人妻久久久久久| av在线观看视频网站免费| 爱豆传媒免费全集在线观看| 国产高清国产精品国产三级 | 久久久久久久亚洲中文字幕| 大又大粗又爽又黄少妇毛片口| 97在线人人人人妻| 国产成人a区在线观看| 国产成人精品一,二区| 亚洲中文av在线| 天天躁夜夜躁狠狠久久av| 青春草亚洲视频在线观看| 美女视频免费永久观看网站| 夫妻性生交免费视频一级片| 中文欧美无线码| 国产精品av视频在线免费观看| 亚洲精品aⅴ在线观看| 这个男人来自地球电影免费观看 | 国产黄色免费在线视频| 国产黄片视频在线免费观看| 亚洲国产精品国产精品| 女人久久www免费人成看片| 大话2 男鬼变身卡| 天天躁夜夜躁狠狠久久av| 久久热精品热| 日本色播在线视频| 一级毛片 在线播放| 夫妻午夜视频| 国产黄片视频在线免费观看| 一本一本综合久久| 天堂俺去俺来也www色官网| a级毛色黄片| 精品亚洲乱码少妇综合久久| 一级av片app| 老女人水多毛片| 久久亚洲国产成人精品v| 国产探花极品一区二区| 国产一区二区三区综合在线观看 | 日韩伦理黄色片| 国产精品爽爽va在线观看网站| av不卡在线播放| 秋霞在线观看毛片| 涩涩av久久男人的天堂| 免费黄网站久久成人精品| 国产爽快片一区二区三区| av在线老鸭窝| 人人妻人人看人人澡| 国产av码专区亚洲av| 街头女战士在线观看网站| 十分钟在线观看高清视频www | 久久亚洲国产成人精品v| 色5月婷婷丁香| 亚洲精品亚洲一区二区| 草草在线视频免费看| 国产精品一及| 一级二级三级毛片免费看| 亚洲婷婷狠狠爱综合网| 久久久色成人| 超碰97精品在线观看| 国产午夜精品一二区理论片| 国产免费一级a男人的天堂| 97在线视频观看| 好男人视频免费观看在线| 久久人人爽av亚洲精品天堂 | 建设人人有责人人尽责人人享有的 | 成人黄色视频免费在线看| 身体一侧抽搐| 观看av在线不卡| 最黄视频免费看| av又黄又爽大尺度在线免费看| 舔av片在线| 亚洲欧美清纯卡通| 欧美激情国产日韩精品一区| 少妇人妻 视频| 大片免费播放器 马上看| 国产淫片久久久久久久久| 久久久久精品性色| 18禁动态无遮挡网站| 亚洲av中文字字幕乱码综合| 久久精品国产亚洲av涩爱| 成年美女黄网站色视频大全免费 | 亚洲久久久国产精品| 精华霜和精华液先用哪个| 噜噜噜噜噜久久久久久91| 最近的中文字幕免费完整| 欧美3d第一页| av黄色大香蕉| 久久婷婷青草| 国产成人精品福利久久| 男女边摸边吃奶| 久久久久久伊人网av| 久久6这里有精品| 亚洲精品成人av观看孕妇| 精品人妻熟女av久视频| 欧美3d第一页| 又爽又黄a免费视频| 人人妻人人看人人澡| 亚洲,欧美,日韩| 男女边摸边吃奶| 哪个播放器可以免费观看大片| 国产成人91sexporn| 日韩强制内射视频| 99久久综合免费| 国产乱来视频区| 老师上课跳d突然被开到最大视频| 国产av一区二区精品久久 | 国产成人精品久久久久久| 亚洲美女黄色视频免费看| 深爱激情五月婷婷| 国产精品99久久99久久久不卡 | 国产亚洲一区二区精品| 又黄又爽又刺激的免费视频.| 99精国产麻豆久久婷婷| 久久久久久伊人网av| 好男人视频免费观看在线| 99热这里只有精品一区| 中文资源天堂在线| 亚洲第一区二区三区不卡| 国产成人91sexporn| 噜噜噜噜噜久久久久久91| 免费观看的影片在线观看| 亚洲精品中文字幕在线视频 | 99久久综合免费| 在线亚洲精品国产二区图片欧美 | 高清av免费在线| 欧美性感艳星| 国产69精品久久久久777片| 国产大屁股一区二区在线视频| 国产在线男女| 肉色欧美久久久久久久蜜桃| 视频区图区小说| 国产精品一区www在线观看| 网址你懂的国产日韩在线| 国产精品三级大全| 国产精品99久久99久久久不卡 | 精品国产一区二区三区久久久樱花 | 99热这里只有精品一区| 成人亚洲欧美一区二区av| 亚洲国产欧美在线一区| 久久精品国产亚洲网站| 一级毛片aaaaaa免费看小| 久久精品人妻少妇| 有码 亚洲区| 国产精品久久久久久久久免| 成人免费观看视频高清| 少妇人妻 视频| 精品一区二区三区视频在线| 大又大粗又爽又黄少妇毛片口| 91午夜精品亚洲一区二区三区| 免费久久久久久久精品成人欧美视频 | 人人妻人人看人人澡| 国产成人freesex在线| 亚洲国产高清在线一区二区三| 国产男人的电影天堂91| 女人十人毛片免费观看3o分钟| 国产免费又黄又爽又色| 亚洲精品aⅴ在线观看| 97在线视频观看| a 毛片基地| 成人一区二区视频在线观看| 大香蕉久久网| 国产精品无大码| 精品久久久久久久久av| 亚洲国产最新在线播放| 国产在线一区二区三区精| 高清午夜精品一区二区三区| 精品久久久久久久末码| 欧美人与善性xxx| 亚洲一级一片aⅴ在线观看| 久久久久网色| 久久久久久久亚洲中文字幕| 大陆偷拍与自拍| 哪个播放器可以免费观看大片| 久久久久人妻精品一区果冻| 最近2019中文字幕mv第一页| 免费人成在线观看视频色| 国产高清不卡午夜福利| 高清日韩中文字幕在线| 国产亚洲精品久久久com| 久久国产精品男人的天堂亚洲 | 欧美日本视频| 久久热精品热| 精品一区在线观看国产| 精品人妻熟女av久视频| 亚洲av中文字字幕乱码综合| 夜夜看夜夜爽夜夜摸| 亚洲国产精品专区欧美| av卡一久久| 丝瓜视频免费看黄片| 中国美白少妇内射xxxbb| 亚洲精品视频女| 激情 狠狠 欧美| 久久久色成人| 精品亚洲乱码少妇综合久久| 性高湖久久久久久久久免费观看| 色吧在线观看| 干丝袜人妻中文字幕| 亚洲成色77777| 日韩中文字幕视频在线看片 | 午夜免费男女啪啪视频观看| 插逼视频在线观看| 欧美bdsm另类| 精品少妇久久久久久888优播| 亚洲精品中文字幕在线视频 | 男人狂女人下面高潮的视频| 五月开心婷婷网| 97超视频在线观看视频| 六月丁香七月| 精品人妻一区二区三区麻豆| 欧美三级亚洲精品| 色吧在线观看| av不卡在线播放| 美女中出高潮动态图| 精品人妻偷拍中文字幕| 我要看日韩黄色一级片| 99久国产av精品国产电影| 亚洲人成网站高清观看| 日韩欧美 国产精品| 国产精品三级大全| 久久精品熟女亚洲av麻豆精品| 日韩中文字幕视频在线看片 | 国产精品国产三级国产av玫瑰| 特大巨黑吊av在线直播| 国产精品久久久久久久电影| 爱豆传媒免费全集在线观看| 午夜激情福利司机影院| 久久婷婷青草| 老司机影院毛片| 日韩制服骚丝袜av| 99久久综合免费| 精品亚洲成国产av| 高清黄色对白视频在线免费看 | 精品少妇久久久久久888优播| 国产毛片在线视频| 久久精品久久精品一区二区三区| 亚洲精华国产精华液的使用体验| 少妇人妻久久综合中文| 最近的中文字幕免费完整| av专区在线播放| 免费观看在线日韩| 国产精品一及| 国产免费一区二区三区四区乱码| 高清在线视频一区二区三区| 国产av精品麻豆| 又粗又硬又长又爽又黄的视频| 天堂俺去俺来也www色官网| 精品熟女少妇av免费看| 国产一区亚洲一区在线观看| 国内精品宾馆在线| 精品久久久噜噜| 99精国产麻豆久久婷婷| 亚洲欧洲国产日韩| 国产精品熟女久久久久浪| 91精品一卡2卡3卡4卡| 亚洲国产av新网站| 免费看av在线观看网站| 欧美bdsm另类| 极品少妇高潮喷水抽搐| 国产伦精品一区二区三区视频9| 一级av片app| 欧美日韩综合久久久久久| 黑人猛操日本美女一级片| 国产 一区 欧美 日韩| 男女下面进入的视频免费午夜| 亚洲欧美清纯卡通| 国产成人freesex在线| 久久人人爽人人爽人人片va| av国产久精品久网站免费入址| 亚洲国产精品999| 黑人高潮一二区| 成人二区视频| 日韩av不卡免费在线播放| 一二三四中文在线观看免费高清| 国产成人a区在线观看| 国产一区二区三区av在线| 插阴视频在线观看视频| 在线观看免费视频网站a站| 汤姆久久久久久久影院中文字幕| 亚洲欧美日韩卡通动漫| 亚洲精品一区蜜桃| 久久99精品国语久久久| 2021少妇久久久久久久久久久| av视频免费观看在线观看| 少妇人妻 视频| 国产精品久久久久久av不卡| 久久亚洲国产成人精品v| 久久久久久久久大av| 色视频www国产| av福利片在线观看| 成人漫画全彩无遮挡| 纯流量卡能插随身wifi吗| 在线观看免费日韩欧美大片 | 人妻 亚洲 视频| 男女下面进入的视频免费午夜| 国产在视频线精品| 成人午夜精彩视频在线观看| 成人黄色视频免费在线看| 黄色视频在线播放观看不卡| 国产精品一区二区在线观看99| 亚洲精品久久午夜乱码| 日韩欧美一区视频在线观看 | 日本爱情动作片www.在线观看| 国产成人a区在线观看| 又大又黄又爽视频免费| 亚洲欧美中文字幕日韩二区| 精品少妇黑人巨大在线播放| 又大又黄又爽视频免费| 熟女人妻精品中文字幕| 久久亚洲国产成人精品v| 最近最新中文字幕免费大全7| 免费在线观看成人毛片| 成人亚洲精品一区在线观看 | 一级爰片在线观看| 激情五月婷婷亚洲| 亚洲人成网站在线观看播放| 97超视频在线观看视频| 国产69精品久久久久777片| 国产av精品麻豆| 日本猛色少妇xxxxx猛交久久| 国产 精品1| 黑丝袜美女国产一区| 国产精品精品国产色婷婷| 亚洲国产精品国产精品| 午夜福利在线在线| 午夜福利影视在线免费观看| 欧美激情极品国产一区二区三区 | 激情五月婷婷亚洲| 亚洲图色成人| 亚洲精品久久久久久婷婷小说| 国模一区二区三区四区视频| 欧美老熟妇乱子伦牲交| 亚洲国产欧美在线一区| 亚洲人成网站在线播| 欧美老熟妇乱子伦牲交| 国产午夜精品一二区理论片| 最新中文字幕久久久久| 亚洲欧美日韩另类电影网站 | 免费播放大片免费观看视频在线观看| 久久久a久久爽久久v久久| 中文精品一卡2卡3卡4更新| 国产高清三级在线| 草草在线视频免费看| 国产精品久久久久久久久免| 久久久午夜欧美精品| 国产在视频线精品| 亚洲经典国产精华液单| 男女国产视频网站| 涩涩av久久男人的天堂| 亚洲美女黄色视频免费看| 中文乱码字字幕精品一区二区三区| 十八禁网站网址无遮挡 | 日韩大片免费观看网站| 蜜桃亚洲精品一区二区三区| 国产亚洲5aaaaa淫片| 国产无遮挡羞羞视频在线观看| 日日啪夜夜撸| 精品一区二区三区视频在线| 亚洲经典国产精华液单| 国产精品精品国产色婷婷| 男女无遮挡免费网站观看| 黄片wwwwww| 毛片女人毛片| 99re6热这里在线精品视频| 久久影院123| 夜夜爽夜夜爽视频| 亚洲真实伦在线观看| av女优亚洲男人天堂| 久久久久久久久久人人人人人人| 日韩强制内射视频| 久久久久久久亚洲中文字幕| 亚洲欧洲国产日韩| 少妇人妻精品综合一区二区| 亚洲av二区三区四区| 18禁在线播放成人免费| 日韩欧美一区视频在线观看 | 亚洲国产精品专区欧美| 亚洲不卡免费看| 91aial.com中文字幕在线观看| 水蜜桃什么品种好| 成人亚洲欧美一区二区av| 久久97久久精品| 尾随美女入室| 亚洲美女视频黄频| 超碰97精品在线观看| 婷婷色麻豆天堂久久| 在线免费观看不下载黄p国产| 久久热精品热| 99热6这里只有精品| 成人毛片a级毛片在线播放| 性色avwww在线观看| 一本色道久久久久久精品综合| 国产精品久久久久久精品电影小说 | 亚洲电影在线观看av| 嫩草影院新地址| h视频一区二区三区| 丰满人妻一区二区三区视频av| 观看免费一级毛片| 七月丁香在线播放| 亚洲av成人精品一区久久| 亚洲精品乱久久久久久| 免费黄频网站在线观看国产| 18禁在线播放成人免费| 你懂的网址亚洲精品在线观看| 国产视频内射| 一本—道久久a久久精品蜜桃钙片| 中文字幕久久专区| 激情 狠狠 欧美| 日韩,欧美,国产一区二区三区| av卡一久久| 51国产日韩欧美| 国产欧美日韩精品一区二区| 麻豆国产97在线/欧美| 国产成人午夜福利电影在线观看| 噜噜噜噜噜久久久久久91| 成人高潮视频无遮挡免费网站| 六月丁香七月| 国产精品久久久久久久电影| 亚洲综合精品二区| 交换朋友夫妻互换小说| 干丝袜人妻中文字幕| www.色视频.com| av专区在线播放| 一级二级三级毛片免费看| h日本视频在线播放| 老司机影院成人| 欧美精品一区二区大全| a 毛片基地| av播播在线观看一区| 成年av动漫网址| 十分钟在线观看高清视频www | 国产午夜精品久久久久久一区二区三区| 少妇 在线观看| 成人黄色视频免费在线看| 美女高潮的动态| 免费高清在线观看视频在线观看| 97精品久久久久久久久久精品| 欧美精品一区二区大全| 中文乱码字字幕精品一区二区三区| 欧美区成人在线视频| 一区二区三区免费毛片| 久久女婷五月综合色啪小说| 久久韩国三级中文字幕| 国产欧美亚洲国产| 免费观看无遮挡的男女| 免费大片黄手机在线观看| 久久久久久久精品精品| 精品国产露脸久久av麻豆| 亚洲精品第二区| 在线观看一区二区三区| 超碰97精品在线观看| 美女国产视频在线观看| 日韩一区二区视频免费看| 校园人妻丝袜中文字幕| 熟女电影av网| 毛片女人毛片| 99视频精品全部免费 在线| 国产淫片久久久久久久久| 国产精品福利在线免费观看| 欧美3d第一页| 亚洲综合色惰| 免费人成在线观看视频色| 久久ye,这里只有精品| www.av在线官网国产| 亚洲四区av| 日产精品乱码卡一卡2卡三| 中文字幕制服av| 亚洲自偷自拍三级| 日日摸夜夜添夜夜添av毛片| 91久久精品国产一区二区成人| 天堂中文最新版在线下载| 韩国av在线不卡| 亚洲av成人精品一二三区| 欧美一级a爱片免费观看看| 中国国产av一级| 国产男女超爽视频在线观看| 亚洲欧美清纯卡通| 亚洲精品乱码久久久久久按摩| 91精品国产九色| 黑人高潮一二区| 国产伦理片在线播放av一区| 久久午夜福利片| 日韩av在线免费看完整版不卡| 在线亚洲精品国产二区图片欧美 | 久久女婷五月综合色啪小说| 在线免费十八禁| 美女主播在线视频| 99热这里只有是精品50|