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

    基于紅外光譜概率神經(jīng)網(wǎng)絡(luò)的訶子產(chǎn)地鑒別

    2014-12-12 09:57安生梅吳啟勛吉守祥
    湖北農(nóng)業(yè)科學(xué) 2014年20期
    關(guān)鍵詞:小波變換訶子

    安生梅+吳啟勛+吉守祥

    摘要:利用傅里葉變換紅外光譜,測定了不同產(chǎn)地訶子(Terminalia chebula Retz.)樣品的紅外光譜圖。采用常規(guī)預(yù)處理方法和小波變換對紅外光譜原始數(shù)據(jù)進(jìn)行了預(yù)處理,并采用主成分分析進(jìn)一步壓縮光譜數(shù)據(jù),前3個主成分的累積貢獻(xiàn)率為98.054%。以前3個主成分作為網(wǎng)絡(luò)輸入,訶子產(chǎn)地類別作為網(wǎng)絡(luò)輸出,建立了概率神經(jīng)網(wǎng)絡(luò),同時對建立該網(wǎng)絡(luò)模型的擴展常數(shù)進(jìn)行了分析。模型分析表明,建立的網(wǎng)絡(luò)模型能夠?qū)?0個訶子樣品進(jìn)行產(chǎn)地鑒別,紅外光譜法結(jié)合神經(jīng)網(wǎng)絡(luò)可作為中藥材產(chǎn)地分類鑒別的一種新的現(xiàn)代化方法。

    關(guān)鍵詞:訶子(Terminalia chebula Retz.);傅里葉變換紅外光譜;小波變換;概率神經(jīng)網(wǎng)絡(luò)

    中圖分類號:O657.33 ? ? ? ?文獻(xiàn)標(biāo)識碼:A ? ? ? ?文章編號:0439-8114(2014)20-4977-03

    DOI:10.14088/j.cnki.issn0439-8114.2014.20.056

    Identifying the Origin of Terminalia chebula Retz. Based on Infrared Spectrum

    and Probability Neural Networks

    AN Sheng-mei,WU Qi-xun,JI Shou-xiang

    (College of Chemistry and Life Sciences,Qinghai University for Nationalities,Xining ?810007,China)

    Abstract: Infrared spectrum of Terminalia chebula Retz. from different fields were determined by Fourier transform infrared spectroscopy (FTIR). The original data matrix of FTIR were pretreated with common preprocessing and wavelet transform. The spectra variables were compressed through the wavelet transformation. The principal component analysis (PCA) method was used to compress the spectral data. The PCA results showed that the first three principal components had the cumulative reliability of 98.054%. The first 3 principal components were used for input nods of probability neural network model and the area sorts of Terminalia chebula Retz.were used for parameters of export. The spread of the probability neural networks model were studied in detail. The model distinguished the producing area of the 40 samples of Terminalia chebula Retz. correctly. The infrared spectral technology combined with the artificial neural networks was proved to be a reliable and practical method for identifying of geographical origin of Lycium barbarum L. It will provide a new reference for the identification of traditional chinese medicine.

    Key words:Terminalia chebula Retz.; fourier transform infrared spectroscopy; wavelet transform; probability neural networks

    訶子為使君子科植物訶子(Terminalia chebula Retz.)或絨毛訶子[T.bellerica (Gaertn.) Roxb.]的干燥成熟果實[1],其為《中國藥典》收載品種,訶子原產(chǎn)于印度、緬甸等地,我國云南、廣西、廣東、西藏等地均有分布。訶子苦、酸、澀、平,歸肺、大腸經(jīng),具有澀腸斂肺、降火利咽之功效,用于治療久瀉久痢、便血脫肛、肺虛喘咳、久咳不止、咽痛音啞等。現(xiàn)代藥理試驗證明,訶子具有抑菌、殺滅或抑制惡性腫瘤細(xì)胞等作用[2,3],臨床應(yīng)用廣泛,尤其在藏藥和蒙藥中最為常用[4,5]。本草考證表明,現(xiàn)今訶子品種和產(chǎn)地均已增多,云南省為現(xiàn)今訶子主產(chǎn)地,所產(chǎn)的絨毛訶子為全國主流品種[6]。但其質(zhì)量鑒別至今仍處于被動狀態(tài),僅有一些薄層、色譜及紫外光譜方法鑒別訶子質(zhì)量的報道[7,8]。截至目前,針對訶子質(zhì)量鑒別比較系統(tǒng)的現(xiàn)代紅外光譜方法未見報道。藥材的紅外光譜是其各個成分的光譜疊加而成的,是藥材本質(zhì)特征的反映。不同產(chǎn)地訶子因是同一品種組成相似,紅外光譜整體上有較大的相似性;但產(chǎn)地不同,一些組成又存在差異,必然反映在紅外光譜的指紋區(qū)存在差異。這種細(xì)微差異很難單純通過譜圖觀察區(qū)分其特征,通過一些化學(xué)計量學(xué)方法則能夠?qū)崿F(xiàn)對不同產(chǎn)地藥材的自動識別。本研究根據(jù)不同產(chǎn)地訶子紅外光譜數(shù)據(jù),結(jié)合小波變換對紅外光譜數(shù)據(jù)進(jìn)行消噪和壓縮,采用概率神經(jīng)網(wǎng)絡(luò)對小波變換處理后的紅外光譜數(shù)據(jù)進(jìn)行分析,并實現(xiàn)了對訶子的產(chǎn)地鑒別。endprint

    1 ?訶子紅外光譜分析

    1.1 ?材料

    本研究選用40個不同產(chǎn)地的訶子樣品,采自國內(nèi)云南、廣西、廣東、西藏等地區(qū)及國外緬甸、印度等地。經(jīng)青海藏醫(yī)藥研究所庫進(jìn)良研究員鑒定為使君子科植物訶子或絨毛訶子的干燥成熟果實即訶子。樣品來源見表1,將采集的干果置于普通干燥箱內(nèi)60 ℃干燥至恒重,取出用微型植物粉碎機粉碎裝袋密封備用。

    1.2 ?儀器及藥品

    儀器:IRPrestge-21型傅里葉變換紅外光譜儀,購自日本島津公司,分辨率為0.5~1.0 cm,信噪比為40 000∶1。檢測器采用新型高靈敏度DLATGS檢測器,可實現(xiàn)FTIR分析的高靈敏度和良好的穩(wěn)定性。

    藥品:光譜純溴化鉀購自天津天光光學(xué)儀器有限公司。

    1.3 ?不同產(chǎn)地訶子紅外光譜測定

    取已烘干至恒重粉碎的訶子樣品,過100目篩,并取1 mg過篩后的訶子粉與100 mg光譜純的溴化鉀在研缽中混合均勻后壓片,用IRPrestge-21型傅里葉變換紅外光譜儀測定紅外光譜。測定范圍4 000~400 cm-1,掃描信號累加16次,分辨率4 cm-1。為了降低誤差,譜圖基線采用儀器自帶軟件自動進(jìn)行多點基線校正和平滑處理。每個樣品平行測定3次,取其平均值作為該樣品的原始光譜,得到不同產(chǎn)地訶子的紅外光譜數(shù)據(jù)及圖像(圖1)。

    2 ?人工神經(jīng)網(wǎng)絡(luò)建模

    2.1 ?常規(guī)紅外光譜的預(yù)處理

    采用標(biāo)準(zhǔn)正態(tài)變換、多元散射校正等常規(guī)紅外光譜的預(yù)處理方法[9],對訶子紅外光譜進(jìn)行了預(yù)處理。雖然能夠校正因散射而引起的光譜誤差,但總體上看,變換后消噪并不理想,也不能用于光譜數(shù)據(jù)壓縮,效果有限。

    2.2 ?紅外光譜的小波變換預(yù)處理

    本研究將小波變換用于訶子紅外光譜的降噪和壓縮[10],均在在MATLAB小波工具箱下通過編程實現(xiàn)。采用均方根誤差RMS評價小波信號降噪和壓縮性能[11]:

    RMS=■■(A■■-A■■)2■

    其中,A■■為降噪或壓縮后數(shù)據(jù)重構(gòu)光譜第k波長的吸光度,A■■為原始光譜第k波長的吸光度。原始信號與消噪或壓縮信號的均方根誤差RMS越小,則消噪或壓縮信號就越接近于原始信號,消噪或壓縮效果越好。

    通過對常用小波基函數(shù)的篩選,選取高階正交緊支集Daubechies(db5)小波作為小波變換的母函數(shù),因為支集長度太短不消失矩太低,不利于信號能量的集中;而支集太長,通常會產(chǎn)生邊界問題,不利于消噪,選階數(shù)N=5是合適的。wmaxlev函數(shù)確定小波分解時的最高分解層數(shù)(本研究選定的層數(shù)N=3),采用thselect函數(shù)對比分析4種不同的小波閾值降噪模型對光譜數(shù)據(jù)進(jìn)行降噪的效果(RMS評價),確定采用啟發(fā)式SURE(heursure)方法獲取閾值,然后對光譜信號從第一層到第N層的高頻系數(shù)進(jìn)行閾值量化處理。實際小波變換后,原始信號消噪與壓縮信號的均方根誤差RMS很小,平均達(dá)到10-6~10-5數(shù)量級,與光譜分析誤差相當(dāng),表明降噪、壓縮效果良好。光譜經(jīng)壓縮,數(shù)據(jù)長度由1 868壓縮至241,為原始數(shù)據(jù)的1/8。小波變換重構(gòu)后的壓縮數(shù)據(jù),消去了噪聲,但仍保留了原始光譜數(shù)據(jù)的信息(圖2),可以用于人工神經(jīng)網(wǎng)絡(luò)方法鑒別訶子產(chǎn)地。

    2.3 ?人工神經(jīng)網(wǎng)絡(luò)建模

    2.3.1 ?訶子紅外光譜的系統(tǒng)聚類分析 ?從訶子小波降噪和壓縮后的樣品紅外光譜數(shù)據(jù)矩陣標(biāo)準(zhǔn)化后的數(shù)據(jù)出發(fā),選定標(biāo)準(zhǔn)化歐氏距離(seuclidean)作為樣品相似性的測度,采用離差平方和(ward)法作為類與類之間的距離,應(yīng)用MATLAB編程,作系統(tǒng)聚類分析。結(jié)果表明,30個樣品明顯聚為4類(國外的4個樣品歸入到國內(nèi)的產(chǎn)地類別)。其中,8個云南主產(chǎn)地訶子樣品和2個緬甸樣品、8個廣西訶子樣品分別聚為一類,廣東、西藏和印度訶子樣品也分別聚為一類,但樣品之間分類有交叉,其中有兩個廣東樣品,落入西藏類別。因此擬采用分類準(zhǔn)確度高的人工神經(jīng)網(wǎng)絡(luò)建立產(chǎn)地鑒別模型。

    2.3.2 ?人工神經(jīng)網(wǎng)絡(luò)建模 ?盡管對樣品紅外光譜進(jìn)行了壓縮,但每個樣品仍然有241個數(shù)據(jù),不宜直接作為人工神經(jīng)網(wǎng)絡(luò)的輸入數(shù)據(jù),可采用主成分分析壓縮數(shù)據(jù)[12]。以訶子小波降噪和壓縮后的紅外光譜數(shù)據(jù)矩陣標(biāo)準(zhǔn)化后的數(shù)據(jù)作主成分分析。前3個主成分方差累積貢獻(xiàn)達(dá)98.054%(λ1=161.52,λ2= 67.277,λ3=7.517 2),即取前3個主成分就能保留原紅外光譜數(shù)據(jù)矩陣的98.054%信息。因涉及多產(chǎn)地鑒別,選用以分類見長的概率神經(jīng)網(wǎng)絡(luò)進(jìn)行訶子產(chǎn)地的鑒別。MATLAB神經(jīng)網(wǎng)絡(luò)工具箱中newpnn函數(shù)可以方便地創(chuàng)建概率神經(jīng)網(wǎng)絡(luò)。

    主要網(wǎng)絡(luò)參數(shù)如下:網(wǎng)絡(luò)結(jié)構(gòu)3-30-4,即輸入神經(jīng)元3個,對應(yīng)于訶子訓(xùn)練樣品的前3個主成分,隱層30個神經(jīng)元,對應(yīng)訶子訓(xùn)練樣品個數(shù),輸出層神經(jīng)元4個,對應(yīng)訶子樣品的類別,以1.0代表云南訶子樣品,2.0代表廣西訶子樣品,以3.0代表廣東訶子樣品,4.0代表西藏訶子樣品。2個緬甸樣品劃分為云南1.0類別,2個印度樣品劃分為西藏4.0類別。隱層的傳遞函數(shù)是徑向基函數(shù)(radbas),輸出層為競爭層,競爭傳遞函數(shù)compet的功能是找出隱層徑向基傳遞函數(shù)中各元素的最大值,并且使與最大值對應(yīng)類別的神經(jīng)元輸出為1,其他類別的神經(jīng)元輸出為0。

    3 ?小結(jié)

    以測定的不同產(chǎn)地訶子紅外光譜為研究對象,采用小波變換技術(shù),均方根誤差RMS作為降噪與壓縮性能評價指標(biāo),對紅外光譜進(jìn)行了消噪和壓縮,小波變換重構(gòu)后的壓縮數(shù)據(jù),既保持了原紅外光譜的特征,又有效地減少了光譜變量。據(jù)此結(jié)合概率神經(jīng)網(wǎng)絡(luò)建立了不同訶子產(chǎn)地鑒別模型。該網(wǎng)絡(luò)模型對產(chǎn)地鑒別預(yù)測效果很好,能鑒別訶子主產(chǎn)地或偽劣訶子,正確識別率達(dá)到100%。紅外光譜結(jié)合適當(dāng)?shù)幕瘜W(xué)計量學(xué)如人工神經(jīng)網(wǎng)絡(luò)方法,具有制樣簡單方便、測定快捷、譜圖準(zhǔn)確、重現(xiàn)性好等優(yōu)點,適合大量樣品的快速鑒別研究,有一定的可靠性和實用性,從而為訶子分類鑒別、質(zhì)量控制提供了新的適用方法,也為鑒別其他中草藥提供了思路。

    參考文獻(xiàn):

    [1] 中華人民共和國衛(wèi)生部藥典委員會.中華人民共和國藥典(一部)[M].北京:化學(xué)工業(yè)出版社,2005.129.

    [2] 周同惠.紙色譜和薄層色譜[M].北京:科學(xué)出版社,1989.

    [3] 王金華,孫芳云,杜冠華.訶子的藥理作用研究進(jìn)展[J].中國民族醫(yī)藥志,2006,12(6):45-47.

    [4] 顏玉貞,謝培山,宋力飛,等.西青果藥材及提取物液相指紋圖譜的應(yīng)用研究[J].中成藥,2004,26(8):6066-6071.

    [5] 周 寧,彭富全.中藏藥中訶子、毛訶子和余甘子的本草考證[J]. 中草藥,2001,32(4):355-357.

    [6] 陳建南,徐鴻華.訶子的產(chǎn)地和品種考證[J].中藥材,1996,19(10):533-535.

    [7] 丁 ?崗,陸蘊如,冀春茹,等.訶子的薄層色譜鑒別[J].北京中醫(yī)藥大學(xué)學(xué)報,2001,24(1):45-46.

    [8] 尹春梅,韓世柳,田成旺,等.訶子HPLC指紋圖譜及模式識別研究[J].中草藥,2010,41(9):1467-1469.

    [9] 尼 ?珍,胡昌勤,馮 ?芳.近紅外光譜分析中光譜預(yù)處理方法的作用及其發(fā)展[J].藥物分析雜志,2008,28(5):824-829.

    [10] 李丹婷,張長江,汪 ?勁.基于小波變換的紅外光譜分析[J].光譜學(xué)與光譜分析,2006,26(11):2024-2026.

    [11] 李夢龍,康 彬,戚華溢.紅外光譜小波壓縮性能的評價[J].高等學(xué)?;瘜W(xué)學(xué)報,2002,23(7):1281-1284.

    [12] 曹紅翠,吳啟勛.中藥甘草中微量元素的主成分分析和系統(tǒng)聚類分析[J].甘肅聯(lián)合大學(xué)學(xué)報(自然科學(xué)版),2010,24(2):61-63.endprint

    1 ?訶子紅外光譜分析

    1.1 ?材料

    本研究選用40個不同產(chǎn)地的訶子樣品,采自國內(nèi)云南、廣西、廣東、西藏等地區(qū)及國外緬甸、印度等地。經(jīng)青海藏醫(yī)藥研究所庫進(jìn)良研究員鑒定為使君子科植物訶子或絨毛訶子的干燥成熟果實即訶子。樣品來源見表1,將采集的干果置于普通干燥箱內(nèi)60 ℃干燥至恒重,取出用微型植物粉碎機粉碎裝袋密封備用。

    1.2 ?儀器及藥品

    儀器:IRPrestge-21型傅里葉變換紅外光譜儀,購自日本島津公司,分辨率為0.5~1.0 cm,信噪比為40 000∶1。檢測器采用新型高靈敏度DLATGS檢測器,可實現(xiàn)FTIR分析的高靈敏度和良好的穩(wěn)定性。

    藥品:光譜純溴化鉀購自天津天光光學(xué)儀器有限公司。

    1.3 ?不同產(chǎn)地訶子紅外光譜測定

    取已烘干至恒重粉碎的訶子樣品,過100目篩,并取1 mg過篩后的訶子粉與100 mg光譜純的溴化鉀在研缽中混合均勻后壓片,用IRPrestge-21型傅里葉變換紅外光譜儀測定紅外光譜。測定范圍4 000~400 cm-1,掃描信號累加16次,分辨率4 cm-1。為了降低誤差,譜圖基線采用儀器自帶軟件自動進(jìn)行多點基線校正和平滑處理。每個樣品平行測定3次,取其平均值作為該樣品的原始光譜,得到不同產(chǎn)地訶子的紅外光譜數(shù)據(jù)及圖像(圖1)。

    2 ?人工神經(jīng)網(wǎng)絡(luò)建模

    2.1 ?常規(guī)紅外光譜的預(yù)處理

    采用標(biāo)準(zhǔn)正態(tài)變換、多元散射校正等常規(guī)紅外光譜的預(yù)處理方法[9],對訶子紅外光譜進(jìn)行了預(yù)處理。雖然能夠校正因散射而引起的光譜誤差,但總體上看,變換后消噪并不理想,也不能用于光譜數(shù)據(jù)壓縮,效果有限。

    2.2 ?紅外光譜的小波變換預(yù)處理

    本研究將小波變換用于訶子紅外光譜的降噪和壓縮[10],均在在MATLAB小波工具箱下通過編程實現(xiàn)。采用均方根誤差RMS評價小波信號降噪和壓縮性能[11]:

    RMS=■■(A■■-A■■)2■

    其中,A■■為降噪或壓縮后數(shù)據(jù)重構(gòu)光譜第k波長的吸光度,A■■為原始光譜第k波長的吸光度。原始信號與消噪或壓縮信號的均方根誤差RMS越小,則消噪或壓縮信號就越接近于原始信號,消噪或壓縮效果越好。

    通過對常用小波基函數(shù)的篩選,選取高階正交緊支集Daubechies(db5)小波作為小波變換的母函數(shù),因為支集長度太短不消失矩太低,不利于信號能量的集中;而支集太長,通常會產(chǎn)生邊界問題,不利于消噪,選階數(shù)N=5是合適的。wmaxlev函數(shù)確定小波分解時的最高分解層數(shù)(本研究選定的層數(shù)N=3),采用thselect函數(shù)對比分析4種不同的小波閾值降噪模型對光譜數(shù)據(jù)進(jìn)行降噪的效果(RMS評價),確定采用啟發(fā)式SURE(heursure)方法獲取閾值,然后對光譜信號從第一層到第N層的高頻系數(shù)進(jìn)行閾值量化處理。實際小波變換后,原始信號消噪與壓縮信號的均方根誤差RMS很小,平均達(dá)到10-6~10-5數(shù)量級,與光譜分析誤差相當(dāng),表明降噪、壓縮效果良好。光譜經(jīng)壓縮,數(shù)據(jù)長度由1 868壓縮至241,為原始數(shù)據(jù)的1/8。小波變換重構(gòu)后的壓縮數(shù)據(jù),消去了噪聲,但仍保留了原始光譜數(shù)據(jù)的信息(圖2),可以用于人工神經(jīng)網(wǎng)絡(luò)方法鑒別訶子產(chǎn)地。

    2.3 ?人工神經(jīng)網(wǎng)絡(luò)建模

    2.3.1 ?訶子紅外光譜的系統(tǒng)聚類分析 ?從訶子小波降噪和壓縮后的樣品紅外光譜數(shù)據(jù)矩陣標(biāo)準(zhǔn)化后的數(shù)據(jù)出發(fā),選定標(biāo)準(zhǔn)化歐氏距離(seuclidean)作為樣品相似性的測度,采用離差平方和(ward)法作為類與類之間的距離,應(yīng)用MATLAB編程,作系統(tǒng)聚類分析。結(jié)果表明,30個樣品明顯聚為4類(國外的4個樣品歸入到國內(nèi)的產(chǎn)地類別)。其中,8個云南主產(chǎn)地訶子樣品和2個緬甸樣品、8個廣西訶子樣品分別聚為一類,廣東、西藏和印度訶子樣品也分別聚為一類,但樣品之間分類有交叉,其中有兩個廣東樣品,落入西藏類別。因此擬采用分類準(zhǔn)確度高的人工神經(jīng)網(wǎng)絡(luò)建立產(chǎn)地鑒別模型。

    2.3.2 ?人工神經(jīng)網(wǎng)絡(luò)建模 ?盡管對樣品紅外光譜進(jìn)行了壓縮,但每個樣品仍然有241個數(shù)據(jù),不宜直接作為人工神經(jīng)網(wǎng)絡(luò)的輸入數(shù)據(jù),可采用主成分分析壓縮數(shù)據(jù)[12]。以訶子小波降噪和壓縮后的紅外光譜數(shù)據(jù)矩陣標(biāo)準(zhǔn)化后的數(shù)據(jù)作主成分分析。前3個主成分方差累積貢獻(xiàn)達(dá)98.054%(λ1=161.52,λ2= 67.277,λ3=7.517 2),即取前3個主成分就能保留原紅外光譜數(shù)據(jù)矩陣的98.054%信息。因涉及多產(chǎn)地鑒別,選用以分類見長的概率神經(jīng)網(wǎng)絡(luò)進(jìn)行訶子產(chǎn)地的鑒別。MATLAB神經(jīng)網(wǎng)絡(luò)工具箱中newpnn函數(shù)可以方便地創(chuàng)建概率神經(jīng)網(wǎng)絡(luò)。

    主要網(wǎng)絡(luò)參數(shù)如下:網(wǎng)絡(luò)結(jié)構(gòu)3-30-4,即輸入神經(jīng)元3個,對應(yīng)于訶子訓(xùn)練樣品的前3個主成分,隱層30個神經(jīng)元,對應(yīng)訶子訓(xùn)練樣品個數(shù),輸出層神經(jīng)元4個,對應(yīng)訶子樣品的類別,以1.0代表云南訶子樣品,2.0代表廣西訶子樣品,以3.0代表廣東訶子樣品,4.0代表西藏訶子樣品。2個緬甸樣品劃分為云南1.0類別,2個印度樣品劃分為西藏4.0類別。隱層的傳遞函數(shù)是徑向基函數(shù)(radbas),輸出層為競爭層,競爭傳遞函數(shù)compet的功能是找出隱層徑向基傳遞函數(shù)中各元素的最大值,并且使與最大值對應(yīng)類別的神經(jīng)元輸出為1,其他類別的神經(jīng)元輸出為0。

    3 ?小結(jié)

    以測定的不同產(chǎn)地訶子紅外光譜為研究對象,采用小波變換技術(shù),均方根誤差RMS作為降噪與壓縮性能評價指標(biāo),對紅外光譜進(jìn)行了消噪和壓縮,小波變換重構(gòu)后的壓縮數(shù)據(jù),既保持了原紅外光譜的特征,又有效地減少了光譜變量。據(jù)此結(jié)合概率神經(jīng)網(wǎng)絡(luò)建立了不同訶子產(chǎn)地鑒別模型。該網(wǎng)絡(luò)模型對產(chǎn)地鑒別預(yù)測效果很好,能鑒別訶子主產(chǎn)地或偽劣訶子,正確識別率達(dá)到100%。紅外光譜結(jié)合適當(dāng)?shù)幕瘜W(xué)計量學(xué)如人工神經(jīng)網(wǎng)絡(luò)方法,具有制樣簡單方便、測定快捷、譜圖準(zhǔn)確、重現(xiàn)性好等優(yōu)點,適合大量樣品的快速鑒別研究,有一定的可靠性和實用性,從而為訶子分類鑒別、質(zhì)量控制提供了新的適用方法,也為鑒別其他中草藥提供了思路。

    參考文獻(xiàn):

    [1] 中華人民共和國衛(wèi)生部藥典委員會.中華人民共和國藥典(一部)[M].北京:化學(xué)工業(yè)出版社,2005.129.

    [2] 周同惠.紙色譜和薄層色譜[M].北京:科學(xué)出版社,1989.

    [3] 王金華,孫芳云,杜冠華.訶子的藥理作用研究進(jìn)展[J].中國民族醫(yī)藥志,2006,12(6):45-47.

    [4] 顏玉貞,謝培山,宋力飛,等.西青果藥材及提取物液相指紋圖譜的應(yīng)用研究[J].中成藥,2004,26(8):6066-6071.

    [5] 周 寧,彭富全.中藏藥中訶子、毛訶子和余甘子的本草考證[J]. 中草藥,2001,32(4):355-357.

    [6] 陳建南,徐鴻華.訶子的產(chǎn)地和品種考證[J].中藥材,1996,19(10):533-535.

    [7] 丁 ?崗,陸蘊如,冀春茹,等.訶子的薄層色譜鑒別[J].北京中醫(yī)藥大學(xué)學(xué)報,2001,24(1):45-46.

    [8] 尹春梅,韓世柳,田成旺,等.訶子HPLC指紋圖譜及模式識別研究[J].中草藥,2010,41(9):1467-1469.

    [9] 尼 ?珍,胡昌勤,馮 ?芳.近紅外光譜分析中光譜預(yù)處理方法的作用及其發(fā)展[J].藥物分析雜志,2008,28(5):824-829.

    [10] 李丹婷,張長江,汪 ?勁.基于小波變換的紅外光譜分析[J].光譜學(xué)與光譜分析,2006,26(11):2024-2026.

    [11] 李夢龍,康 彬,戚華溢.紅外光譜小波壓縮性能的評價[J].高等學(xué)?;瘜W(xué)學(xué)報,2002,23(7):1281-1284.

    [12] 曹紅翠,吳啟勛.中藥甘草中微量元素的主成分分析和系統(tǒng)聚類分析[J].甘肅聯(lián)合大學(xué)學(xué)報(自然科學(xué)版),2010,24(2):61-63.endprint

    1 ?訶子紅外光譜分析

    1.1 ?材料

    本研究選用40個不同產(chǎn)地的訶子樣品,采自國內(nèi)云南、廣西、廣東、西藏等地區(qū)及國外緬甸、印度等地。經(jīng)青海藏醫(yī)藥研究所庫進(jìn)良研究員鑒定為使君子科植物訶子或絨毛訶子的干燥成熟果實即訶子。樣品來源見表1,將采集的干果置于普通干燥箱內(nèi)60 ℃干燥至恒重,取出用微型植物粉碎機粉碎裝袋密封備用。

    1.2 ?儀器及藥品

    儀器:IRPrestge-21型傅里葉變換紅外光譜儀,購自日本島津公司,分辨率為0.5~1.0 cm,信噪比為40 000∶1。檢測器采用新型高靈敏度DLATGS檢測器,可實現(xiàn)FTIR分析的高靈敏度和良好的穩(wěn)定性。

    藥品:光譜純溴化鉀購自天津天光光學(xué)儀器有限公司。

    1.3 ?不同產(chǎn)地訶子紅外光譜測定

    取已烘干至恒重粉碎的訶子樣品,過100目篩,并取1 mg過篩后的訶子粉與100 mg光譜純的溴化鉀在研缽中混合均勻后壓片,用IRPrestge-21型傅里葉變換紅外光譜儀測定紅外光譜。測定范圍4 000~400 cm-1,掃描信號累加16次,分辨率4 cm-1。為了降低誤差,譜圖基線采用儀器自帶軟件自動進(jìn)行多點基線校正和平滑處理。每個樣品平行測定3次,取其平均值作為該樣品的原始光譜,得到不同產(chǎn)地訶子的紅外光譜數(shù)據(jù)及圖像(圖1)。

    2 ?人工神經(jīng)網(wǎng)絡(luò)建模

    2.1 ?常規(guī)紅外光譜的預(yù)處理

    采用標(biāo)準(zhǔn)正態(tài)變換、多元散射校正等常規(guī)紅外光譜的預(yù)處理方法[9],對訶子紅外光譜進(jìn)行了預(yù)處理。雖然能夠校正因散射而引起的光譜誤差,但總體上看,變換后消噪并不理想,也不能用于光譜數(shù)據(jù)壓縮,效果有限。

    2.2 ?紅外光譜的小波變換預(yù)處理

    本研究將小波變換用于訶子紅外光譜的降噪和壓縮[10],均在在MATLAB小波工具箱下通過編程實現(xiàn)。采用均方根誤差RMS評價小波信號降噪和壓縮性能[11]:

    RMS=■■(A■■-A■■)2■

    其中,A■■為降噪或壓縮后數(shù)據(jù)重構(gòu)光譜第k波長的吸光度,A■■為原始光譜第k波長的吸光度。原始信號與消噪或壓縮信號的均方根誤差RMS越小,則消噪或壓縮信號就越接近于原始信號,消噪或壓縮效果越好。

    通過對常用小波基函數(shù)的篩選,選取高階正交緊支集Daubechies(db5)小波作為小波變換的母函數(shù),因為支集長度太短不消失矩太低,不利于信號能量的集中;而支集太長,通常會產(chǎn)生邊界問題,不利于消噪,選階數(shù)N=5是合適的。wmaxlev函數(shù)確定小波分解時的最高分解層數(shù)(本研究選定的層數(shù)N=3),采用thselect函數(shù)對比分析4種不同的小波閾值降噪模型對光譜數(shù)據(jù)進(jìn)行降噪的效果(RMS評價),確定采用啟發(fā)式SURE(heursure)方法獲取閾值,然后對光譜信號從第一層到第N層的高頻系數(shù)進(jìn)行閾值量化處理。實際小波變換后,原始信號消噪與壓縮信號的均方根誤差RMS很小,平均達(dá)到10-6~10-5數(shù)量級,與光譜分析誤差相當(dāng),表明降噪、壓縮效果良好。光譜經(jīng)壓縮,數(shù)據(jù)長度由1 868壓縮至241,為原始數(shù)據(jù)的1/8。小波變換重構(gòu)后的壓縮數(shù)據(jù),消去了噪聲,但仍保留了原始光譜數(shù)據(jù)的信息(圖2),可以用于人工神經(jīng)網(wǎng)絡(luò)方法鑒別訶子產(chǎn)地。

    2.3 ?人工神經(jīng)網(wǎng)絡(luò)建模

    2.3.1 ?訶子紅外光譜的系統(tǒng)聚類分析 ?從訶子小波降噪和壓縮后的樣品紅外光譜數(shù)據(jù)矩陣標(biāo)準(zhǔn)化后的數(shù)據(jù)出發(fā),選定標(biāo)準(zhǔn)化歐氏距離(seuclidean)作為樣品相似性的測度,采用離差平方和(ward)法作為類與類之間的距離,應(yīng)用MATLAB編程,作系統(tǒng)聚類分析。結(jié)果表明,30個樣品明顯聚為4類(國外的4個樣品歸入到國內(nèi)的產(chǎn)地類別)。其中,8個云南主產(chǎn)地訶子樣品和2個緬甸樣品、8個廣西訶子樣品分別聚為一類,廣東、西藏和印度訶子樣品也分別聚為一類,但樣品之間分類有交叉,其中有兩個廣東樣品,落入西藏類別。因此擬采用分類準(zhǔn)確度高的人工神經(jīng)網(wǎng)絡(luò)建立產(chǎn)地鑒別模型。

    2.3.2 ?人工神經(jīng)網(wǎng)絡(luò)建模 ?盡管對樣品紅外光譜進(jìn)行了壓縮,但每個樣品仍然有241個數(shù)據(jù),不宜直接作為人工神經(jīng)網(wǎng)絡(luò)的輸入數(shù)據(jù),可采用主成分分析壓縮數(shù)據(jù)[12]。以訶子小波降噪和壓縮后的紅外光譜數(shù)據(jù)矩陣標(biāo)準(zhǔn)化后的數(shù)據(jù)作主成分分析。前3個主成分方差累積貢獻(xiàn)達(dá)98.054%(λ1=161.52,λ2= 67.277,λ3=7.517 2),即取前3個主成分就能保留原紅外光譜數(shù)據(jù)矩陣的98.054%信息。因涉及多產(chǎn)地鑒別,選用以分類見長的概率神經(jīng)網(wǎng)絡(luò)進(jìn)行訶子產(chǎn)地的鑒別。MATLAB神經(jīng)網(wǎng)絡(luò)工具箱中newpnn函數(shù)可以方便地創(chuàng)建概率神經(jīng)網(wǎng)絡(luò)。

    主要網(wǎng)絡(luò)參數(shù)如下:網(wǎng)絡(luò)結(jié)構(gòu)3-30-4,即輸入神經(jīng)元3個,對應(yīng)于訶子訓(xùn)練樣品的前3個主成分,隱層30個神經(jīng)元,對應(yīng)訶子訓(xùn)練樣品個數(shù),輸出層神經(jīng)元4個,對應(yīng)訶子樣品的類別,以1.0代表云南訶子樣品,2.0代表廣西訶子樣品,以3.0代表廣東訶子樣品,4.0代表西藏訶子樣品。2個緬甸樣品劃分為云南1.0類別,2個印度樣品劃分為西藏4.0類別。隱層的傳遞函數(shù)是徑向基函數(shù)(radbas),輸出層為競爭層,競爭傳遞函數(shù)compet的功能是找出隱層徑向基傳遞函數(shù)中各元素的最大值,并且使與最大值對應(yīng)類別的神經(jīng)元輸出為1,其他類別的神經(jīng)元輸出為0。

    3 ?小結(jié)

    以測定的不同產(chǎn)地訶子紅外光譜為研究對象,采用小波變換技術(shù),均方根誤差RMS作為降噪與壓縮性能評價指標(biāo),對紅外光譜進(jìn)行了消噪和壓縮,小波變換重構(gòu)后的壓縮數(shù)據(jù),既保持了原紅外光譜的特征,又有效地減少了光譜變量。據(jù)此結(jié)合概率神經(jīng)網(wǎng)絡(luò)建立了不同訶子產(chǎn)地鑒別模型。該網(wǎng)絡(luò)模型對產(chǎn)地鑒別預(yù)測效果很好,能鑒別訶子主產(chǎn)地或偽劣訶子,正確識別率達(dá)到100%。紅外光譜結(jié)合適當(dāng)?shù)幕瘜W(xué)計量學(xué)如人工神經(jīng)網(wǎng)絡(luò)方法,具有制樣簡單方便、測定快捷、譜圖準(zhǔn)確、重現(xiàn)性好等優(yōu)點,適合大量樣品的快速鑒別研究,有一定的可靠性和實用性,從而為訶子分類鑒別、質(zhì)量控制提供了新的適用方法,也為鑒別其他中草藥提供了思路。

    參考文獻(xiàn):

    [1] 中華人民共和國衛(wèi)生部藥典委員會.中華人民共和國藥典(一部)[M].北京:化學(xué)工業(yè)出版社,2005.129.

    [2] 周同惠.紙色譜和薄層色譜[M].北京:科學(xué)出版社,1989.

    [3] 王金華,孫芳云,杜冠華.訶子的藥理作用研究進(jìn)展[J].中國民族醫(yī)藥志,2006,12(6):45-47.

    [4] 顏玉貞,謝培山,宋力飛,等.西青果藥材及提取物液相指紋圖譜的應(yīng)用研究[J].中成藥,2004,26(8):6066-6071.

    [5] 周 寧,彭富全.中藏藥中訶子、毛訶子和余甘子的本草考證[J]. 中草藥,2001,32(4):355-357.

    [6] 陳建南,徐鴻華.訶子的產(chǎn)地和品種考證[J].中藥材,1996,19(10):533-535.

    [7] 丁 ?崗,陸蘊如,冀春茹,等.訶子的薄層色譜鑒別[J].北京中醫(yī)藥大學(xué)學(xué)報,2001,24(1):45-46.

    [8] 尹春梅,韓世柳,田成旺,等.訶子HPLC指紋圖譜及模式識別研究[J].中草藥,2010,41(9):1467-1469.

    [9] 尼 ?珍,胡昌勤,馮 ?芳.近紅外光譜分析中光譜預(yù)處理方法的作用及其發(fā)展[J].藥物分析雜志,2008,28(5):824-829.

    [10] 李丹婷,張長江,汪 ?勁.基于小波變換的紅外光譜分析[J].光譜學(xué)與光譜分析,2006,26(11):2024-2026.

    [11] 李夢龍,康 彬,戚華溢.紅外光譜小波壓縮性能的評價[J].高等學(xué)?;瘜W(xué)學(xué)報,2002,23(7):1281-1284.

    [12] 曹紅翠,吳啟勛.中藥甘草中微量元素的主成分分析和系統(tǒng)聚類分析[J].甘肅聯(lián)合大學(xué)學(xué)報(自然科學(xué)版),2010,24(2):61-63.endprint

    猜你喜歡
    小波變換訶子
    永德縣訶子產(chǎn)業(yè)發(fā)展對策研究
    訶子的臨床應(yīng)用及其用量探究
    永德縣訶子產(chǎn)業(yè)化開發(fā)現(xiàn)狀及發(fā)展對策
    基于雙樹四元數(shù)小波變換的圖像降噪增強
    訶子化學(xué)成分、生物活性及分析方法研究進(jìn)展
    金訶子低溫提取物對神經(jīng)膠質(zhì)瘤C6細(xì)胞的抑制作用
    訶子與毛訶子的鑒別
    亚洲aⅴ乱码一区二区在线播放| 亚洲av国产av综合av卡| 国内少妇人妻偷人精品xxx网站| 26uuu在线亚洲综合色| 精品久久久久久久久av| 在线免费十八禁| 中文字幕免费在线视频6| 国产伦理片在线播放av一区| 国产精品福利在线免费观看| 午夜免费鲁丝| 国产精品国产三级国产专区5o| 久久精品熟女亚洲av麻豆精品| 欧美激情极品国产一区二区三区 | 777米奇影视久久| 综合色丁香网| 一本久久精品| 亚洲第一av免费看| 1000部很黄的大片| 777米奇影视久久| 夜夜骑夜夜射夜夜干| 熟妇人妻不卡中文字幕| 国产高清不卡午夜福利| 一本—道久久a久久精品蜜桃钙片| 国产免费视频播放在线视频| 一级黄片播放器| 99热这里只有是精品50| 国产高潮美女av| 日韩大片免费观看网站| 国产老妇伦熟女老妇高清| 免费黄频网站在线观看国产| 亚洲欧美精品自产自拍| 国产成人免费观看mmmm| 色婷婷久久久亚洲欧美| 日韩一区二区视频免费看| 身体一侧抽搐| 国产一区亚洲一区在线观看| 熟女人妻精品中文字幕| 日韩在线高清观看一区二区三区| 欧美日韩视频精品一区| 亚洲成人手机| 国产 一区 欧美 日韩| 黄色日韩在线| 亚洲精品乱久久久久久| 国产精品久久久久成人av| 又爽又黄a免费视频| 免费av不卡在线播放| 一区二区av电影网| 夜夜骑夜夜射夜夜干| 国产色爽女视频免费观看| 最近中文字幕2019免费版| 日韩制服骚丝袜av| 亚洲国产最新在线播放| 最近最新中文字幕大全电影3| 久久女婷五月综合色啪小说| 亚洲精品色激情综合| 春色校园在线视频观看| 秋霞伦理黄片| 汤姆久久久久久久影院中文字幕| 观看美女的网站| 久久精品久久久久久久性| 欧美日本视频| 欧美bdsm另类| 国模一区二区三区四区视频| 欧美一级a爱片免费观看看| av.在线天堂| 亚洲自偷自拍三级| 少妇人妻精品综合一区二区| 内射极品少妇av片p| 亚洲欧美日韩无卡精品| 卡戴珊不雅视频在线播放| 午夜精品国产一区二区电影| 大陆偷拍与自拍| 美女国产视频在线观看| 丝袜脚勾引网站| 国产伦精品一区二区三区四那| 男的添女的下面高潮视频| 九色成人免费人妻av| 国产成人aa在线观看| 80岁老熟妇乱子伦牲交| 久久99蜜桃精品久久| 午夜免费鲁丝| 久久久久久人妻| 亚洲经典国产精华液单| 少妇人妻久久综合中文| 亚洲国产色片| 麻豆成人av视频| 男女啪啪激烈高潮av片| 精品国产三级普通话版| 大片电影免费在线观看免费| 丰满人妻一区二区三区视频av| 午夜激情久久久久久久| 丰满乱子伦码专区| 精品一区在线观看国产| 91久久精品国产一区二区三区| 日韩欧美精品免费久久| 日本av手机在线免费观看| 最近手机中文字幕大全| 免费黄色在线免费观看| 亚洲精品456在线播放app| 国产久久久一区二区三区| 亚洲国产毛片av蜜桃av| 一级爰片在线观看| 亚洲欧美成人精品一区二区| 老司机影院毛片| 国产女主播在线喷水免费视频网站| 国产男女内射视频| 婷婷色av中文字幕| 免费大片黄手机在线观看| 一区二区av电影网| 国产精品99久久99久久久不卡 | 九九爱精品视频在线观看| 国产黄片美女视频| 国产亚洲欧美精品永久| 亚洲欧美一区二区三区国产| 色综合色国产| 欧美精品国产亚洲| 妹子高潮喷水视频| 欧美三级亚洲精品| 天堂8中文在线网| 18禁动态无遮挡网站| 久久久久久久久久久免费av| 亚洲欧洲日产国产| 久久久久久久国产电影| 我的老师免费观看完整版| 丰满乱子伦码专区| 亚洲av在线观看美女高潮| 高清午夜精品一区二区三区| 中国国产av一级| 国产一级毛片在线| 在线看a的网站| a级毛片免费高清观看在线播放| 人妻系列 视频| 免费看光身美女| 搡女人真爽免费视频火全软件| 亚洲va在线va天堂va国产| 精品一区二区三卡| 欧美日韩精品成人综合77777| 毛片一级片免费看久久久久| 精品亚洲成国产av| 免费高清在线观看视频在线观看| 最近中文字幕高清免费大全6| av线在线观看网站| 亚洲av中文av极速乱| 一级毛片久久久久久久久女| 亚洲av在线观看美女高潮| 高清黄色对白视频在线免费看 | 99热6这里只有精品| 人妻一区二区av| 亚洲欧美成人综合另类久久久| 夜夜骑夜夜射夜夜干| 久久久精品免费免费高清| 国产精品成人在线| 高清欧美精品videossex| 日韩欧美一区视频在线观看 | 观看美女的网站| 亚洲精品一区蜜桃| 亚洲av国产av综合av卡| 午夜福利影视在线免费观看| 欧美+日韩+精品| 亚洲国产av新网站| 99热6这里只有精品| 色综合色国产| 久久ye,这里只有精品| 中文字幕免费在线视频6| 欧美xxⅹ黑人| 国产免费视频播放在线视频| 久久婷婷青草| 久久99热这里只有精品18| 久久久久精品性色| 亚洲色图av天堂| 精品午夜福利在线看| 国产精品国产三级国产专区5o| 看十八女毛片水多多多| 97在线视频观看| 国产无遮挡羞羞视频在线观看| 波野结衣二区三区在线| 搡老乐熟女国产| 久久久久久人妻| 国产女主播在线喷水免费视频网站| 一级二级三级毛片免费看| 亚洲国产精品国产精品| 国产精品偷伦视频观看了| videossex国产| 午夜免费观看性视频| 日本午夜av视频| 我要看黄色一级片免费的| 欧美高清成人免费视频www| 男女下面进入的视频免费午夜| 久久久色成人| 亚洲性久久影院| 一级毛片 在线播放| 久久精品夜色国产| 国产有黄有色有爽视频| 天美传媒精品一区二区| 国产伦精品一区二区三区四那| av国产免费在线观看| av卡一久久| 日韩电影二区| 亚洲怡红院男人天堂| 久久久久国产精品人妻一区二区| 精品国产露脸久久av麻豆| 多毛熟女@视频| 亚洲人成网站在线播| 国产高潮美女av| 女人十人毛片免费观看3o分钟| 欧美97在线视频| 亚洲熟女精品中文字幕| 亚洲精品国产av成人精品| 免费看不卡的av| 菩萨蛮人人尽说江南好唐韦庄| 少妇 在线观看| 久久女婷五月综合色啪小说| 天堂8中文在线网| 国产高清三级在线| 欧美bdsm另类| 国产亚洲一区二区精品| 午夜福利网站1000一区二区三区| 久久午夜福利片| 蜜桃在线观看..| 美女国产视频在线观看| 18禁动态无遮挡网站| 边亲边吃奶的免费视频| 精品久久久精品久久久| 亚洲,欧美,日韩| 日韩三级伦理在线观看| 高清毛片免费看| 男的添女的下面高潮视频| 欧美日韩视频高清一区二区三区二| 免费看av在线观看网站| 尾随美女入室| 久久人妻熟女aⅴ| 最黄视频免费看| 在线 av 中文字幕| 九九在线视频观看精品| 国产黄片视频在线免费观看| av视频免费观看在线观看| 亚洲伊人久久精品综合| 国产精品久久久久久精品古装| 成年美女黄网站色视频大全免费 | 又大又黄又爽视频免费| 久久精品久久久久久久性| 欧美xxⅹ黑人| 亚洲婷婷狠狠爱综合网| 亚洲在久久综合| 自拍欧美九色日韩亚洲蝌蚪91 | 国产成人精品福利久久| 深爱激情五月婷婷| 日韩一区二区视频免费看| 久久久成人免费电影| 韩国av在线不卡| av.在线天堂| av在线播放精品| h视频一区二区三区| 日本wwww免费看| 亚洲精品第二区| 大陆偷拍与自拍| 青青草视频在线视频观看| 国产成人a区在线观看| 观看av在线不卡| 午夜激情福利司机影院| 亚洲av日韩在线播放| 美女福利国产在线 | 国产欧美另类精品又又久久亚洲欧美| 久久精品国产亚洲av天美| 久久久午夜欧美精品| 毛片女人毛片| 老司机影院毛片| 少妇裸体淫交视频免费看高清| 青春草视频在线免费观看| 日本wwww免费看| 欧美丝袜亚洲另类| 日本-黄色视频高清免费观看| 欧美zozozo另类| 日韩在线高清观看一区二区三区| 国产精品欧美亚洲77777| 国产在线一区二区三区精| 国产一区二区三区av在线| 人妻少妇偷人精品九色| 我的老师免费观看完整版| 美女内射精品一级片tv| 成人免费观看视频高清| 在线观看美女被高潮喷水网站| 久久精品熟女亚洲av麻豆精品| 精品一区二区免费观看| 蜜桃久久精品国产亚洲av| 日本与韩国留学比较| 亚洲国产色片| 熟女电影av网| 性色avwww在线观看| 一个人看视频在线观看www免费| 欧美国产精品一级二级三级 | 免费观看在线日韩| 91在线精品国自产拍蜜月| 亚洲av中文字字幕乱码综合| 日本猛色少妇xxxxx猛交久久| 精品亚洲成a人片在线观看 | 亚洲av中文字字幕乱码综合| 最近最新中文字幕免费大全7| 亚洲国产成人一精品久久久| 2022亚洲国产成人精品| 一级黄片播放器| 人妻夜夜爽99麻豆av| 欧美成人一区二区免费高清观看| .国产精品久久| 久久99热这里只频精品6学生| 国产精品人妻久久久久久| 高清av免费在线| av黄色大香蕉| 国产成人91sexporn| 欧美极品一区二区三区四区| 精品国产乱码久久久久久小说| 成年免费大片在线观看| 蜜臀久久99精品久久宅男| 国产人妻一区二区三区在| 国产有黄有色有爽视频| 亚洲综合精品二区| .国产精品久久| 国精品久久久久久国模美| 久久av网站| 一级毛片电影观看| 少妇的逼水好多| 大片电影免费在线观看免费| 中文字幕久久专区| 久久人人爽人人爽人人片va| 美女中出高潮动态图| 一本—道久久a久久精品蜜桃钙片| 亚洲美女搞黄在线观看| 日韩精品有码人妻一区| 99热这里只有精品一区| 视频区图区小说| 晚上一个人看的免费电影| 国产亚洲5aaaaa淫片| 最近2019中文字幕mv第一页| 我的老师免费观看完整版| 中文字幕精品免费在线观看视频 | 熟女人妻精品中文字幕| 亚洲av日韩在线播放| 我要看黄色一级片免费的| 最近中文字幕2019免费版| 欧美成人a在线观看| 国产一区二区三区综合在线观看 | 欧美最新免费一区二区三区| videossex国产| 久久久久久久久久久免费av| 亚洲高清免费不卡视频| 国产男女超爽视频在线观看| 少妇猛男粗大的猛烈进出视频| 久久鲁丝午夜福利片| 午夜免费鲁丝| 黄色配什么色好看| 国产精品久久久久久av不卡| 黑丝袜美女国产一区| 亚洲成人手机| 大陆偷拍与自拍| videossex国产| 欧美最新免费一区二区三区| 久久午夜福利片| 亚洲经典国产精华液单| 能在线免费看毛片的网站| 久久国产精品大桥未久av | 伊人久久国产一区二区| 99久久中文字幕三级久久日本| 日日啪夜夜爽| 一本—道久久a久久精品蜜桃钙片| 国产乱人偷精品视频| 精品一区二区免费观看| 国产色爽女视频免费观看| 亚洲一级一片aⅴ在线观看| 国产大屁股一区二区在线视频| 99热网站在线观看| 亚洲av在线观看美女高潮| 九九久久精品国产亚洲av麻豆| 国产91av在线免费观看| 少妇人妻久久综合中文| 亚洲人成网站高清观看| 2018国产大陆天天弄谢| av专区在线播放| 有码 亚洲区| 日韩一区二区视频免费看| 麻豆国产97在线/欧美| 日韩一区二区视频免费看| 夫妻午夜视频| 国产伦理片在线播放av一区| 久久精品久久精品一区二区三区| 三级经典国产精品| 日韩免费高清中文字幕av| 欧美三级亚洲精品| 丝袜脚勾引网站| 国内少妇人妻偷人精品xxx网站| 水蜜桃什么品种好| 汤姆久久久久久久影院中文字幕| 婷婷色av中文字幕| 成人国产麻豆网| 亚洲av在线观看美女高潮| 久久久久久伊人网av| 国产亚洲av片在线观看秒播厂| 免费不卡的大黄色大毛片视频在线观看| 国产大屁股一区二区在线视频| 18禁裸乳无遮挡动漫免费视频| 在线观看一区二区三区激情| 国国产精品蜜臀av免费| 色5月婷婷丁香| 校园人妻丝袜中文字幕| 男女无遮挡免费网站观看| 在线观看一区二区三区| 熟女人妻精品中文字幕| www.av在线官网国产| 一本色道久久久久久精品综合| 97精品久久久久久久久久精品| 九九在线视频观看精品| 久久精品久久精品一区二区三区| 久久婷婷青草| 九九久久精品国产亚洲av麻豆| 三级国产精品片| 永久网站在线| 亚洲aⅴ乱码一区二区在线播放| 好男人视频免费观看在线| 人妻一区二区av| 日韩av免费高清视频| 久久久久久久久久久丰满| 欧美日韩在线观看h| 水蜜桃什么品种好| 大陆偷拍与自拍| 免费观看性生交大片5| 久热这里只有精品99| 亚洲一级一片aⅴ在线观看| 欧美3d第一页| 22中文网久久字幕| 亚洲精品第二区| 黄片无遮挡物在线观看| 在线亚洲精品国产二区图片欧美 | 国产亚洲5aaaaa淫片| 男人狂女人下面高潮的视频| 亚洲国产av新网站| 成人综合一区亚洲| 国产高清国产精品国产三级 | 性色avwww在线观看| 亚洲欧美精品自产自拍| 少妇熟女欧美另类| 26uuu在线亚洲综合色| 精品人妻一区二区三区麻豆| 高清av免费在线| 国产成人午夜福利电影在线观看| 啦啦啦中文免费视频观看日本| 夫妻性生交免费视频一级片| 男人和女人高潮做爰伦理| 男的添女的下面高潮视频| 国产乱来视频区| 久久国产乱子免费精品| 高清在线视频一区二区三区| 日韩av免费高清视频| 在线播放无遮挡| 日韩成人av中文字幕在线观看| 一个人看的www免费观看视频| 插阴视频在线观看视频| 日韩av不卡免费在线播放| 热re99久久精品国产66热6| 最新中文字幕久久久久| 99热这里只有是精品50| 狂野欧美激情性xxxx在线观看| 免费黄色在线免费观看| 成人毛片60女人毛片免费| 在线 av 中文字幕| 国产精品成人在线| 高清在线视频一区二区三区| 国产亚洲午夜精品一区二区久久| 亚洲性久久影院| 91精品伊人久久大香线蕉| 插逼视频在线观看| 熟女人妻精品中文字幕| 日本午夜av视频| 97在线人人人人妻| 国产亚洲91精品色在线| 熟女av电影| 亚洲精品久久午夜乱码| 亚洲中文av在线| 久久青草综合色| 欧美人与善性xxx| 国产免费一区二区三区四区乱码| 插逼视频在线观看| 亚洲欧美清纯卡通| 肉色欧美久久久久久久蜜桃| 亚洲精品第二区| 久久精品国产自在天天线| av在线app专区| 欧美国产精品一级二级三级 | 亚洲美女黄色视频免费看| 久久精品国产亚洲av天美| 午夜激情久久久久久久| 99re6热这里在线精品视频| 欧美成人一区二区免费高清观看| 天天躁夜夜躁狠狠久久av| 国产精品久久久久久久久免| 国产精品国产三级国产av玫瑰| 小蜜桃在线观看免费完整版高清| 91久久精品国产一区二区成人| 日本免费在线观看一区| 人妻系列 视频| 亚洲美女视频黄频| 国产伦精品一区二区三区四那| 免费观看性生交大片5| 亚洲国产高清在线一区二区三| 秋霞在线观看毛片| 99久久精品热视频| 在线观看免费视频网站a站| 天美传媒精品一区二区| 精品熟女少妇av免费看| 性色av一级| 日韩制服骚丝袜av| 赤兔流量卡办理| 国产国拍精品亚洲av在线观看| 国产伦理片在线播放av一区| 最后的刺客免费高清国语| 亚洲激情五月婷婷啪啪| 欧美成人午夜免费资源| 国产精品欧美亚洲77777| 岛国毛片在线播放| 久久久久网色| 美女主播在线视频| 青春草视频在线免费观看| 色视频在线一区二区三区| 大陆偷拍与自拍| 国产中年淑女户外野战色| 日韩伦理黄色片| 亚洲第一av免费看| 天美传媒精品一区二区| 制服丝袜香蕉在线| 久热久热在线精品观看| 国产亚洲91精品色在线| 亚洲欧美日韩东京热| 国产精品嫩草影院av在线观看| 日本av免费视频播放| 涩涩av久久男人的天堂| 国产精品精品国产色婷婷| 久久女婷五月综合色啪小说| 少妇人妻精品综合一区二区| 国产伦精品一区二区三区视频9| 日韩欧美一区视频在线观看 | 九九爱精品视频在线观看| 精品久久久久久久久av| 伦理电影大哥的女人| 熟女人妻精品中文字幕| 国产精品av视频在线免费观看| 80岁老熟妇乱子伦牲交| 王馨瑶露胸无遮挡在线观看| 午夜日本视频在线| 亚洲av在线观看美女高潮| 日韩强制内射视频| 国产男女超爽视频在线观看| 日韩一区二区三区影片| 最近最新中文字幕免费大全7| 尾随美女入室| 高清欧美精品videossex| 欧美日韩在线观看h| 国产成人精品一,二区| 人人妻人人添人人爽欧美一区卜 | 国产免费视频播放在线视频| 日韩av不卡免费在线播放| 久久久久久久久久久丰满| 国产精品一区二区在线不卡| 你懂的网址亚洲精品在线观看| 亚洲成人一二三区av| 国产亚洲91精品色在线| 五月开心婷婷网| 午夜激情福利司机影院| 国产精品一及| 久久精品国产亚洲av涩爱| 久久久久网色| 日韩精品有码人妻一区| 午夜免费男女啪啪视频观看| 性色avwww在线观看| 91aial.com中文字幕在线观看| 一个人免费看片子| 国产高清有码在线观看视频| 国产高潮美女av| 丝瓜视频免费看黄片| 舔av片在线| 欧美日韩一区二区视频在线观看视频在线| 妹子高潮喷水视频| 大片电影免费在线观看免费| 三级经典国产精品| 日韩一区二区三区影片| 亚州av有码| www.av在线官网国产| 亚洲人与动物交配视频| 国产精品一及| 十八禁网站网址无遮挡 | 人妻 亚洲 视频| 丰满少妇做爰视频| 国产真实伦视频高清在线观看| 大又大粗又爽又黄少妇毛片口| 国产成人a∨麻豆精品| 国产亚洲5aaaaa淫片| 欧美一级a爱片免费观看看| 卡戴珊不雅视频在线播放| 91午夜精品亚洲一区二区三区| videos熟女内射| 国产黄片视频在线免费观看| 中文字幕av成人在线电影| 老师上课跳d突然被开到最大视频| 少妇裸体淫交视频免费看高清| 欧美精品一区二区大全| 97在线人人人人妻| 天美传媒精品一区二区| 日本-黄色视频高清免费观看| 久久精品夜色国产| 一区在线观看完整版| 日日啪夜夜爽| 丰满人妻一区二区三区视频av| 美女中出高潮动态图| 一级二级三级毛片免费看| 亚洲色图av天堂| 免费观看的影片在线观看| 亚洲精品乱码久久久久久按摩| 如何舔出高潮|