• <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)
    在现免费观看毛片| av在线天堂中文字幕| 噜噜噜噜噜久久久久久91| 国产在视频线精品| 午夜精品国产一区二区电影 | 国产在线男女| 午夜福利网站1000一区二区三区| 亚洲av成人av| 国内精品宾馆在线| 免费人成在线观看视频色| 亚洲成人av在线免费| 久久精品国产自在天天线| or卡值多少钱| 熟女人妻精品中文字幕| 免费观看精品视频网站| 一区二区三区高清视频在线| 亚洲欧洲日产国产| 亚洲av免费高清在线观看| 极品少妇高潮喷水抽搐| 三级男女做爰猛烈吃奶摸视频| 久久精品久久久久久噜噜老黄| 偷拍熟女少妇极品色| 国产成人aa在线观看| 国产精品麻豆人妻色哟哟久久 | 男人爽女人下面视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产毛片a区久久久久| 亚洲欧美精品专区久久| 啦啦啦中文免费视频观看日本| 99热这里只有是精品在线观看| 日韩电影二区| 午夜日本视频在线| 日韩强制内射视频| 麻豆av噜噜一区二区三区| 国产成人午夜福利电影在线观看| 亚洲精品亚洲一区二区| 国产一区二区亚洲精品在线观看| 麻豆乱淫一区二区| 少妇熟女aⅴ在线视频| 黄片无遮挡物在线观看| 久久99热这里只频精品6学生| 国产精品一二三区在线看| 午夜免费男女啪啪视频观看| 日本与韩国留学比较| a级毛片免费高清观看在线播放| 亚洲精品日本国产第一区| 国产亚洲av片在线观看秒播厂 | 亚洲av一区综合| 午夜福利视频精品| 午夜爱爱视频在线播放| 在现免费观看毛片| 一级毛片久久久久久久久女| 狂野欧美激情性xxxx在线观看| 天美传媒精品一区二区| 国产视频首页在线观看| 国产高清有码在线观看视频| 男插女下体视频免费在线播放| 国产成人aa在线观看| 熟妇人妻不卡中文字幕| 欧美一级a爱片免费观看看| 亚洲性久久影院| 亚洲成人av在线免费| 亚洲精品日韩在线中文字幕| 26uuu在线亚洲综合色| 国产亚洲91精品色在线| 久久97久久精品| 成年免费大片在线观看| 久久鲁丝午夜福利片| 久久国产乱子免费精品| 在线观看免费高清a一片| 一级av片app| 特级一级黄色大片| 日韩制服骚丝袜av| 亚洲精品日本国产第一区| 大又大粗又爽又黄少妇毛片口| 国内精品美女久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲av嫩草精品影院| 久久久a久久爽久久v久久| 亚洲国产高清在线一区二区三| 免费不卡的大黄色大毛片视频在线观看 | 精品久久久噜噜| 国产毛片a区久久久久| 欧美激情国产日韩精品一区| av在线亚洲专区| 国产一区有黄有色的免费视频 | 最近中文字幕2019免费版| 日韩欧美 国产精品| 亚洲av中文av极速乱| 久久精品夜夜夜夜夜久久蜜豆| 国产 亚洲一区二区三区 | 国产一区二区亚洲精品在线观看| 在线观看免费高清a一片| 免费看日本二区| 99久国产av精品| 国产午夜精品一二区理论片| 在线播放无遮挡| 99久久中文字幕三级久久日本| 天美传媒精品一区二区| 亚洲人与动物交配视频| 久99久视频精品免费| 国产精品蜜桃在线观看| 国产午夜精品论理片| 寂寞人妻少妇视频99o| 免费av毛片视频| 听说在线观看完整版免费高清| 欧美成人a在线观看| 韩国高清视频一区二区三区| 美女被艹到高潮喷水动态| 肉色欧美久久久久久久蜜桃 | 91aial.com中文字幕在线观看| 麻豆国产97在线/欧美| 狂野欧美激情性xxxx在线观看| 男人和女人高潮做爰伦理| 免费在线观看成人毛片| 久久久久久久久久久丰满| 国产精品不卡视频一区二区| 在线天堂最新版资源| 亚洲成人精品中文字幕电影| 日韩欧美三级三区| 欧美一区二区亚洲| 狂野欧美白嫩少妇大欣赏| 精品99又大又爽又粗少妇毛片| 少妇熟女aⅴ在线视频| 自拍偷自拍亚洲精品老妇| 亚洲国产精品专区欧美| 九色成人免费人妻av| 亚洲欧美日韩卡通动漫| 午夜免费激情av| 又爽又黄无遮挡网站| 亚洲av成人精品一区久久| 国产精品99久久久久久久久| 国产乱人视频| 麻豆国产97在线/欧美| 毛片一级片免费看久久久久| 99久久精品国产国产毛片| 最近最新中文字幕大全电影3| 亚洲av二区三区四区| 2021天堂中文幕一二区在线观| 97在线视频观看| 三级毛片av免费| 日韩电影二区| av在线亚洲专区| 一区二区三区高清视频在线| 久久人人爽人人爽人人片va| 日本一二三区视频观看| 午夜福利成人在线免费观看| 波多野结衣巨乳人妻| 欧美精品一区二区大全| 国产精品久久久久久精品电影| 青春草亚洲视频在线观看| 国产高清国产精品国产三级 | 亚洲国产日韩欧美精品在线观看| 一个人看视频在线观看www免费| 中文字幕av成人在线电影| 色综合亚洲欧美另类图片| 18禁在线播放成人免费| 非洲黑人性xxxx精品又粗又长| 日韩电影二区| av在线亚洲专区| 久久热精品热| 日韩大片免费观看网站| 波多野结衣巨乳人妻| 三级男女做爰猛烈吃奶摸视频| 日日啪夜夜撸| 久久精品夜夜夜夜夜久久蜜豆| 日韩av免费高清视频| 久久精品久久精品一区二区三区| 超碰97精品在线观看| 可以在线观看毛片的网站| 在线a可以看的网站| 中文欧美无线码| 亚洲av福利一区| 午夜免费激情av| 国产日韩欧美在线精品| 网址你懂的国产日韩在线| 国产av在哪里看| 精品一区二区免费观看| 国精品久久久久久国模美| 国产精品久久久久久久久免| 久久精品熟女亚洲av麻豆精品 | 国产 一区精品| 青春草亚洲视频在线观看| 婷婷色av中文字幕| 久久久欧美国产精品| 免费观看av网站的网址| 亚洲欧美精品专区久久| 校园人妻丝袜中文字幕| 午夜精品国产一区二区电影 | 午夜爱爱视频在线播放| 大又大粗又爽又黄少妇毛片口| 99久久精品一区二区三区| 尾随美女入室| av网站免费在线观看视频 | 日日摸夜夜添夜夜爱| 国产精品美女特级片免费视频播放器| 成人欧美大片| 国产高潮美女av| 中文精品一卡2卡3卡4更新| 欧美日韩国产mv在线观看视频 | 久久久久精品久久久久真实原创| 激情 狠狠 欧美| 国产精品伦人一区二区| 久久久久免费精品人妻一区二区| 蜜臀久久99精品久久宅男| 男女国产视频网站| 精品一区二区三区视频在线| 国产成人a区在线观看| 日本熟妇午夜| 婷婷色综合大香蕉| 亚洲av电影在线观看一区二区三区 | 亚洲精品成人久久久久久| 69av精品久久久久久| 青春草国产在线视频| 十八禁国产超污无遮挡网站| 日韩,欧美,国产一区二区三区| 亚洲乱码一区二区免费版| 亚洲18禁久久av| 最近视频中文字幕2019在线8| 亚洲四区av| 九九爱精品视频在线观看| 波野结衣二区三区在线| 黄片wwwwww| 国产精品日韩av在线免费观看| 国产人妻一区二区三区在| 亚洲三级黄色毛片| 欧美日韩综合久久久久久| 亚洲国产精品专区欧美| av在线天堂中文字幕| 99热这里只有是精品在线观看| 亚洲人成网站在线观看播放| 男女边摸边吃奶| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲欧美日韩无卡精品| 18禁在线无遮挡免费观看视频| 精品99又大又爽又粗少妇毛片| 精品酒店卫生间| 久久久久久伊人网av| 最近中文字幕高清免费大全6| 久久人人爽人人片av| a级一级毛片免费在线观看| 最近中文字幕2019免费版| 亚洲美女视频黄频| 97热精品久久久久久| 中文字幕人妻熟人妻熟丝袜美| 国产毛片a区久久久久| 我的老师免费观看完整版| 丰满乱子伦码专区| 麻豆乱淫一区二区| 久久精品国产亚洲av天美| 日日啪夜夜爽| 三级经典国产精品| 久久久久久久久大av| 亚洲国产精品成人久久小说| 啦啦啦韩国在线观看视频| 亚洲精品自拍成人| 69人妻影院| 99热6这里只有精品| 91狼人影院| 日产精品乱码卡一卡2卡三| av在线播放精品| 又黄又爽又刺激的免费视频.| 中文天堂在线官网| 亚洲综合色惰| 国产在视频线精品| 99热这里只有是精品50| 熟女电影av网| 亚洲精品乱久久久久久| 人妻少妇偷人精品九色| 免费观看精品视频网站| 午夜福利高清视频| 精品久久国产蜜桃| 97精品久久久久久久久久精品| 国产在线男女| 免费看a级黄色片| 永久网站在线| 黑人高潮一二区| 能在线免费观看的黄片| 欧美激情国产日韩精品一区| av国产免费在线观看| 国产精品久久久久久av不卡| 熟妇人妻不卡中文字幕| 美女脱内裤让男人舔精品视频| 在线免费观看的www视频| 亚洲av免费高清在线观看| 波多野结衣巨乳人妻| 少妇人妻一区二区三区视频| 丝袜美腿在线中文| 老师上课跳d突然被开到最大视频| 夜夜看夜夜爽夜夜摸| 少妇的逼好多水| 日韩一区二区三区影片| 久久鲁丝午夜福利片| 久久久成人免费电影| 老司机影院成人| 中文字幕久久专区| 只有这里有精品99| 大片免费播放器 马上看| 精品人妻一区二区三区麻豆| 蜜臀久久99精品久久宅男| 成年人午夜在线观看视频 | 欧美另类一区| 免费人成在线观看视频色| 久久久久久伊人网av| 中文字幕亚洲精品专区| 婷婷色综合大香蕉| 欧美丝袜亚洲另类| 国产黄色免费在线视频| 午夜精品在线福利| 午夜福利在线观看免费完整高清在| 一二三四中文在线观看免费高清| 免费看光身美女| 亚洲av电影在线观看一区二区三区 | 国产探花在线观看一区二区| 五月玫瑰六月丁香| 一级爰片在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久久久久免| 哪个播放器可以免费观看大片| 99热6这里只有精品| 欧美精品一区二区大全| 成人性生交大片免费视频hd| av网站免费在线观看视频 | 国产精品人妻久久久久久| 丰满人妻一区二区三区视频av| 国产v大片淫在线免费观看| 国内少妇人妻偷人精品xxx网站| 国语对白做爰xxxⅹ性视频网站| 精品少妇黑人巨大在线播放| 舔av片在线| 国产精品麻豆人妻色哟哟久久 | 日日啪夜夜爽| 久久久久精品久久久久真实原创| 成人二区视频| 只有这里有精品99| 97精品久久久久久久久久精品| av.在线天堂| 日日摸夜夜添夜夜爱| 欧美高清成人免费视频www| 亚洲精品乱久久久久久| av在线老鸭窝| 国产成人91sexporn| 成人鲁丝片一二三区免费| 午夜福利高清视频| 国产老妇伦熟女老妇高清| 免费av不卡在线播放| 天美传媒精品一区二区| 亚洲精品一二三| 国产男人的电影天堂91| 精品人妻一区二区三区麻豆| 婷婷六月久久综合丁香| 欧美日韩亚洲高清精品| 日本色播在线视频| 成人鲁丝片一二三区免费| 国产永久视频网站| 午夜福利视频1000在线观看| 尾随美女入室| 日韩伦理黄色片| 亚洲高清免费不卡视频| 内地一区二区视频在线| 国产成人freesex在线| 国产真实伦视频高清在线观看| 啦啦啦中文免费视频观看日本| 色视频www国产| 亚洲最大成人中文| a级毛片免费高清观看在线播放| 亚洲人成网站高清观看| 美女国产视频在线观看| 内地一区二区视频在线| 男人舔女人下体高潮全视频| 国产精品伦人一区二区| 国产午夜精品久久久久久一区二区三区| 国产精品久久久久久精品电影小说 | 777米奇影视久久| 天堂影院成人在线观看| 亚洲国产精品专区欧美| 久久久久久久久久久丰满| 精品久久久久久久久亚洲| 韩国高清视频一区二区三区| 国产亚洲5aaaaa淫片| 亚洲欧美中文字幕日韩二区| 纵有疾风起免费观看全集完整版 | 中文字幕久久专区| 国产精品一区二区三区四区久久| 最近最新中文字幕大全电影3| 久久这里只有精品中国| 丝瓜视频免费看黄片| 日本色播在线视频| 久热久热在线精品观看| 美女大奶头视频| 国产单亲对白刺激| 精品国产三级普通话版| 天堂影院成人在线观看| 亚洲av福利一区| 日韩电影二区| 亚洲最大成人av| 久久久欧美国产精品| 最近视频中文字幕2019在线8| 日韩成人伦理影院| 搡老妇女老女人老熟妇| 岛国毛片在线播放| 日韩欧美精品免费久久| 天堂中文最新版在线下载 | 少妇人妻精品综合一区二区| av.在线天堂| 丰满乱子伦码专区| 亚洲综合色惰| 亚洲国产色片| 中国美白少妇内射xxxbb| 自拍偷自拍亚洲精品老妇| 六月丁香七月| 国产成人一区二区在线| 久久久久久久午夜电影| 免费少妇av软件| 亚洲一级一片aⅴ在线观看| 精品久久久久久成人av| 亚洲熟女精品中文字幕| 99热全是精品| 国产永久视频网站| 久久久久九九精品影院| 免费观看精品视频网站| 国产精品一及| 欧美成人午夜免费资源| 日本午夜av视频| 超碰97精品在线观看| 蜜臀久久99精品久久宅男| 久久人人爽人人片av| 成人毛片60女人毛片免费| 亚洲欧美清纯卡通| 小蜜桃在线观看免费完整版高清| xxx大片免费视频| 国产 一区 欧美 日韩| 亚洲av.av天堂| a级毛色黄片| 亚洲激情五月婷婷啪啪| 简卡轻食公司| 欧美三级亚洲精品| 亚洲国产欧美在线一区| 久久精品国产亚洲网站| 中文天堂在线官网| 五月天丁香电影| 精品熟女少妇av免费看| 亚洲婷婷狠狠爱综合网| 黄色日韩在线| 免费观看无遮挡的男女| www.av在线官网国产| 久久精品久久精品一区二区三区| 高清视频免费观看一区二区 | 免费看av在线观看网站| 天天躁夜夜躁狠狠久久av| 国产爱豆传媒在线观看| 国产成人精品婷婷| 蜜桃亚洲精品一区二区三区| av播播在线观看一区| 亚洲av免费在线观看| a级毛色黄片| 最近中文字幕2019免费版| 国产视频首页在线观看| 日韩电影二区| 久久久成人免费电影| 特大巨黑吊av在线直播| 只有这里有精品99| 好男人在线观看高清免费视频| 成年版毛片免费区| 国产av国产精品国产| 伊人久久精品亚洲午夜| 91在线精品国自产拍蜜月| 18+在线观看网站| 欧美性猛交╳xxx乱大交人| 亚洲三级黄色毛片| 国产 一区精品| 亚洲国产色片| 一区二区三区高清视频在线| 一二三四中文在线观看免费高清| 精品久久久噜噜| 又粗又硬又长又爽又黄的视频| 天堂网av新在线| 国产精品人妻久久久影院| 国产大屁股一区二区在线视频| 久久久久精品久久久久真实原创| 一级av片app| 国产黄频视频在线观看| 国产精品一区www在线观看| 亚洲av二区三区四区| 一区二区三区四区激情视频| kizo精华| 久久久久久久午夜电影| 联通29元200g的流量卡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 夫妻午夜视频| 亚洲在线自拍视频| 七月丁香在线播放| 99久久人妻综合| 日韩伦理黄色片| 晚上一个人看的免费电影| 日韩欧美精品v在线| 麻豆av噜噜一区二区三区| 久久99精品国语久久久| 亚洲精品第二区| 中文字幕av成人在线电影| 成人av在线播放网站| 国产一级毛片七仙女欲春2| 毛片女人毛片| 国产免费福利视频在线观看| 毛片女人毛片| 国产午夜福利久久久久久| videossex国产| 免费观看av网站的网址| 建设人人有责人人尽责人人享有的 | 成人二区视频| 在线免费观看的www视频| 亚洲精华国产精华液的使用体验| 六月丁香七月| 男女啪啪激烈高潮av片| 九九爱精品视频在线观看| 床上黄色一级片| 国产一区二区三区综合在线观看 | 18禁在线播放成人免费| 精品国产一区二区三区久久久樱花 | 神马国产精品三级电影在线观看| 内地一区二区视频在线| 色哟哟·www| 国产精品一区二区三区四区免费观看| 国产老妇伦熟女老妇高清| 亚洲国产精品sss在线观看| 中文字幕人妻熟人妻熟丝袜美| 日韩强制内射视频| 亚洲成人av在线免费| 非洲黑人性xxxx精品又粗又长| 亚洲av成人精品一二三区| av专区在线播放| 女人十人毛片免费观看3o分钟| 日本黄大片高清| 国产一级毛片在线| 亚州av有码| 尾随美女入室| 日本一二三区视频观看| av播播在线观看一区| 久久久a久久爽久久v久久| 国产精品女同一区二区软件| 亚洲电影在线观看av| 干丝袜人妻中文字幕| 亚洲精品色激情综合| 亚洲性久久影院| 99久久精品国产国产毛片| 日韩欧美精品v在线| 一级爰片在线观看| 久久99热这里只频精品6学生| 国产免费福利视频在线观看| 欧美 日韩 精品 国产| 天堂√8在线中文| 成人高潮视频无遮挡免费网站| 麻豆成人午夜福利视频| 成人国产麻豆网| 波多野结衣巨乳人妻| 久久午夜福利片| 免费高清在线观看视频在线观看| 老女人水多毛片| 少妇裸体淫交视频免费看高清| 国产一级毛片七仙女欲春2| 国产v大片淫在线免费观看| 国产国拍精品亚洲av在线观看| 黄色一级大片看看| 国产麻豆成人av免费视频| 精品一区二区免费观看| 中国美白少妇内射xxxbb| 日本一二三区视频观看| 哪个播放器可以免费观看大片| 亚洲三级黄色毛片| 午夜福利视频精品| 国产黄色免费在线视频| 观看免费一级毛片| 免费看美女性在线毛片视频| 亚洲电影在线观看av| 极品教师在线视频| 夫妻午夜视频| 大陆偷拍与自拍| 少妇熟女欧美另类| 美女黄网站色视频| 日韩,欧美,国产一区二区三区| 久久久久精品性色| 精品久久久久久电影网| 国产伦理片在线播放av一区| 一本一本综合久久| 成人国产麻豆网| 国产精品不卡视频一区二区| 国产亚洲av嫩草精品影院| 婷婷色综合大香蕉| 欧美变态另类bdsm刘玥| 丰满少妇做爰视频| 亚洲精品日韩av片在线观看| 成人欧美大片| 欧美日韩国产mv在线观看视频 | av在线蜜桃| 国产久久久一区二区三区| 舔av片在线| 亚洲av日韩在线播放| 高清毛片免费看| 国产精品精品国产色婷婷| 午夜福利在线观看免费完整高清在| 中文字幕亚洲精品专区| 91久久精品国产一区二区三区| 日本免费在线观看一区| 特大巨黑吊av在线直播| 综合色av麻豆| 国产成人精品久久久久久| 免费黄网站久久成人精品| 中国国产av一级| 国产一级毛片七仙女欲春2| 内地一区二区视频在线| 美女xxoo啪啪120秒动态图| 亚洲精品aⅴ在线观看| 高清在线视频一区二区三区| 日本av手机在线免费观看| 男女啪啪激烈高潮av片|