基于紅外光譜概率神經網絡的訶子產地鑒別

安生梅+吳啟勛+吉守祥

摘要：利用傅里葉變換紅外光譜，測定了不同產地訶子（Terminalia chebula Retz.）樣品的紅外光譜圖。采用常規預處理方法和小波變換對紅外光譜原始數據進行了預處理，并采用主成分分析進一步壓縮光譜數據，前3個主成分的累積貢獻率為98.054%。以前3個主成分作為網絡輸入，訶子產地類別作為網絡輸出，建立了概率神經網絡，同時對建立該網絡模型的擴展常數進行了分析。模型分析表明，建立的網絡模型能夠對40個訶子樣品進行產地鑒別，紅外光譜法結合神經網絡可作為中藥材產地分類鑒別的一種新的現代化方法。

關鍵詞：訶子（Terminalia chebula Retz.）;傅里葉變換紅外光譜;小波變換;概率神經網絡

中圖分類號：O657.33 ? ? ? ?文獻標識碼：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]，其為《中國藥典》收載品種，訶子原產于印度、緬甸等地，我國云南、廣西、廣東、西藏等地均有分布。訶子苦、酸、澀、平，歸肺、大腸經，具有澀腸斂肺、降火利咽之功效，用于治療久瀉久痢、便血脫肛、肺虛喘咳、久咳不止、咽痛音啞等?，F代藥理試驗證明，訶子具有抑菌、殺滅或抑制惡性腫瘤細胞等作用[2，3]，臨床應用廣泛，尤其在藏藥和蒙藥中最為常用[4，5]。本草考證表明，現今訶子品種和產地均已增多，云南省為現今訶子主產地，所產的絨毛訶子為全國主流品種[6]。但其質量鑒別至今仍處于被動狀態，僅有一些薄層、色譜及紫外光譜方法鑒別訶子質量的報道[7，8]。截至目前，針對訶子質量鑒別比較系統的現代紅外光譜方法未見報道。藥材的紅外光譜是其各個成分的光譜疊加而成的，是藥材本質特征的反映。不同產地訶子因是同一品種組成相似，紅外光譜整體上有較大的相似性;但產地不同，一些組成又存在差異，必然反映在紅外光譜的指紋區存在差異。這種細微差異很難單純通過譜圖觀察區分其特征，通過一些化學計量學方法則能夠實現對不同產地藥材的自動識別。本研究根據不同產地訶子紅外光譜數據，結合小波變換對紅外光譜數據進行消噪和壓縮，采用概率神經網絡對小波變換處理后的紅外光譜數據進行分析，并實現了對訶子的產地鑒別。endprint

1 ?訶子紅外光譜分析

1.1 ?材料

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

1.2 ?儀器及藥品

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

藥品：光譜純溴化鉀購自天津天光光學儀器有限公司。

1.3 ?不同產地訶子紅外光譜測定

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

2 ?人工神經網絡建模

2.1 ?常規紅外光譜的預處理

采用標準正態變換、多元散射校正等常規紅外光譜的預處理方法[9]，對訶子紅外光譜進行了預處理。雖然能夠校正因散射而引起的光譜誤差，但總體上看，變換后消噪并不理想，也不能用于光譜數據壓縮，效果有限。

2.2 ?紅外光譜的小波變換預處理

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

RMS=■■（A■■-A■■）2■

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

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

2.3 ?人工神經網絡建模

2.3.1 ?訶子紅外光譜的系統聚類分析 ?從訶子小波降噪和壓縮后的樣品紅外光譜數據矩陣標準化后的數據出發，選定標準化歐氏距離（seuclidean）作為樣品相似性的測度，采用離差平方和（ward）法作為類與類之間的距離，應用MATLAB編程，作系統聚類分析。結果表明，30個樣品明顯聚為4類（國外的4個樣品歸入到國內的產地類別）。其中，8個云南主產地訶子樣品和2個緬甸樣品、8個廣西訶子樣品分別聚為一類，廣東、西藏和印度訶子樣品也分別聚為一類，但樣品之間分類有交叉，其中有兩個廣東樣品，落入西藏類別。因此擬采用分類準確度高的人工神經網絡建立產地鑒別模型。

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

主要網絡參數如下：網絡結構3-30-4，即輸入神經元3個，對應于訶子訓練樣品的前3個主成分，隱層30個神經元，對應訶子訓練樣品個數，輸出層神經元4個，對應訶子樣品的類別，以1.0代表云南訶子樣品，2.0代表廣西訶子樣品，以3.0代表廣東訶子樣品，4.0代表西藏訶子樣品。2個緬甸樣品劃分為云南1.0類別，2個印度樣品劃分為西藏4.0類別。隱層的傳遞函數是徑向基函數（radbas），輸出層為競爭層，競爭傳遞函數compet的功能是找出隱層徑向基傳遞函數中各元素的最大值，并且使與最大值對應類別的神經元輸出為1，其他類別的神經元輸出為0。

3 ?小結

以測定的不同產地訶子紅外光譜為研究對象，采用小波變換技術，均方根誤差RMS作為降噪與壓縮性能評價指標，對紅外光譜進行了消噪和壓縮，小波變換重構后的壓縮數據，既保持了原紅外光譜的特征，又有效地減少了光譜變量。據此結合概率神經網絡建立了不同訶子產地鑒別模型。該網絡模型對產地鑒別預測效果很好，能鑒別訶子主產地或偽劣訶子，正確識別率達到100%。紅外光譜結合適當的化學計量學如人工神經網絡方法，具有制樣簡單方便、測定快捷、譜圖準確、重現性好等優點，適合大量樣品的快速鑒別研究，有一定的可靠性和實用性，從而為訶子分類鑒別、質量控制提供了新的適用方法，也為鑒別其他中草藥提供了思路。

參考文獻：

[1] 中華人民共和國衛生部藥典委員會.中華人民共和國藥典（一部）[M].北京：化學工業出版社，2005.129.

[2] 周同惠.紙色譜和薄層色譜[M].北京：科學出版社，1989.

[3] 王金華，孫芳云，杜冠華.訶子的藥理作用研究進展[J].中國民族醫藥志，2006，12（6）：45-47.

[4] 顏玉貞，謝培山，宋力飛，等.西青果藥材及提取物液相指紋圖譜的應用研究[J].中成藥，2004，26（8）：6066-6071.

[5] 周 寧，彭富全.中藏藥中訶子、毛訶子和余甘子的本草考證[J]. 中草藥，2001，32（4）：355-357.

[6] 陳建南，徐鴻華.訶子的產地和品種考證[J].中藥材，1996，19（10）：533-535.

[7] 丁 ?崗，陸蘊如，冀春茹，等.訶子的薄層色譜鑒別[J].北京中醫藥大學學報，2001，24（1）：45-46.

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

[9] 尼 ?珍，胡昌勤，馮 ?芳.近紅外光譜分析中光譜預處理方法的作用及其發展[J].藥物分析雜志，2008，28（5）：824-829.

[10] 李丹婷，張長江，汪 ?勁.基于小波變換的紅外光譜分析[J].光譜學與光譜分析，2006，26（11）：2024-2026.

[11] 李夢龍，康 彬，戚華溢.紅外光譜小波壓縮性能的評價[J].高等學校化學學報，2002，23（7）：1281-1284.

[12] 曹紅翠，吳啟勛.中藥甘草中微量元素的主成分分析和系統聚類分析[J].甘肅聯合大學學報（自然科學版），2010，24（2）：61-63.endprint

1 ?訶子紅外光譜分析

1.1 ?材料

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

1.2 ?儀器及藥品

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

藥品：光譜純溴化鉀購自天津天光光學儀器有限公司。

1.3 ?不同產地訶子紅外光譜測定

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

2 ?人工神經網絡建模

2.1 ?常規紅外光譜的預處理

采用標準正態變換、多元散射校正等常規紅外光譜的預處理方法[9]，對訶子紅外光譜進行了預處理。雖然能夠校正因散射而引起的光譜誤差，但總體上看，變換后消噪并不理想，也不能用于光譜數據壓縮，效果有限。

2.2 ?紅外光譜的小波變換預處理

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

RMS=■■（A■■-A■■）2■

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

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

2.3 ?人工神經網絡建模

2.3.1 ?訶子紅外光譜的系統聚類分析 ?從訶子小波降噪和壓縮后的樣品紅外光譜數據矩陣標準化后的數據出發，選定標準化歐氏距離（seuclidean）作為樣品相似性的測度，采用離差平方和（ward）法作為類與類之間的距離，應用MATLAB編程，作系統聚類分析。結果表明，30個樣品明顯聚為4類（國外的4個樣品歸入到國內的產地類別）。其中，8個云南主產地訶子樣品和2個緬甸樣品、8個廣西訶子樣品分別聚為一類，廣東、西藏和印度訶子樣品也分別聚為一類，但樣品之間分類有交叉，其中有兩個廣東樣品，落入西藏類別。因此擬采用分類準確度高的人工神經網絡建立產地鑒別模型。

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

主要網絡參數如下：網絡結構3-30-4，即輸入神經元3個，對應于訶子訓練樣品的前3個主成分，隱層30個神經元，對應訶子訓練樣品個數，輸出層神經元4個，對應訶子樣品的類別，以1.0代表云南訶子樣品，2.0代表廣西訶子樣品，以3.0代表廣東訶子樣品，4.0代表西藏訶子樣品。2個緬甸樣品劃分為云南1.0類別，2個印度樣品劃分為西藏4.0類別。隱層的傳遞函數是徑向基函數（radbas），輸出層為競爭層，競爭傳遞函數compet的功能是找出隱層徑向基傳遞函數中各元素的最大值，并且使與最大值對應類別的神經元輸出為1，其他類別的神經元輸出為0。

3 ?小結

以測定的不同產地訶子紅外光譜為研究對象，采用小波變換技術，均方根誤差RMS作為降噪與壓縮性能評價指標，對紅外光譜進行了消噪和壓縮，小波變換重構后的壓縮數據，既保持了原紅外光譜的特征，又有效地減少了光譜變量。據此結合概率神經網絡建立了不同訶子產地鑒別模型。該網絡模型對產地鑒別預測效果很好，能鑒別訶子主產地或偽劣訶子，正確識別率達到100%。紅外光譜結合適當的化學計量學如人工神經網絡方法，具有制樣簡單方便、測定快捷、譜圖準確、重現性好等優點，適合大量樣品的快速鑒別研究，有一定的可靠性和實用性，從而為訶子分類鑒別、質量控制提供了新的適用方法，也為鑒別其他中草藥提供了思路。

參考文獻：

[1] 中華人民共和國衛生部藥典委員會.中華人民共和國藥典（一部）[M].北京：化學工業出版社，2005.129.

[2] 周同惠.紙色譜和薄層色譜[M].北京：科學出版社，1989.

[3] 王金華，孫芳云，杜冠華.訶子的藥理作用研究進展[J].中國民族醫藥志，2006，12（6）：45-47.

[4] 顏玉貞，謝培山，宋力飛，等.西青果藥材及提取物液相指紋圖譜的應用研究[J].中成藥，2004，26（8）：6066-6071.

[5] 周 寧，彭富全.中藏藥中訶子、毛訶子和余甘子的本草考證[J]. 中草藥，2001，32（4）：355-357.

[6] 陳建南，徐鴻華.訶子的產地和品種考證[J].中藥材，1996，19（10）：533-535.

[7] 丁 ?崗，陸蘊如，冀春茹，等.訶子的薄層色譜鑒別[J].北京中醫藥大學學報，2001，24（1）：45-46.

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

[9] 尼 ?珍，胡昌勤，馮 ?芳.近紅外光譜分析中光譜預處理方法的作用及其發展[J].藥物分析雜志，2008，28（5）：824-829.

[10] 李丹婷，張長江，汪 ?勁.基于小波變換的紅外光譜分析[J].光譜學與光譜分析，2006，26（11）：2024-2026.

[11] 李夢龍，康 彬，戚華溢.紅外光譜小波壓縮性能的評價[J].高等學?；瘜W學報，2002，23（7）：1281-1284.

[12] 曹紅翠，吳啟勛.中藥甘草中微量元素的主成分分析和系統聚類分析[J].甘肅聯合大學學報（自然科學版），2010，24（2）：61-63.endprint

1 ?訶子紅外光譜分析

1.1 ?材料

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

1.2 ?儀器及藥品

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

藥品：光譜純溴化鉀購自天津天光光學儀器有限公司。

1.3 ?不同產地訶子紅外光譜測定

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

2 ?人工神經網絡建模

2.1 ?常規紅外光譜的預處理

采用標準正態變換、多元散射校正等常規紅外光譜的預處理方法[9]，對訶子紅外光譜進行了預處理。雖然能夠校正因散射而引起的光譜誤差，但總體上看，變換后消噪并不理想，也不能用于光譜數據壓縮，效果有限。

2.2 ?紅外光譜的小波變換預處理

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

RMS=■■（A■■-A■■）2■

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

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

2.3 ?人工神經網絡建模

2.3.1 ?訶子紅外光譜的系統聚類分析 ?從訶子小波降噪和壓縮后的樣品紅外光譜數據矩陣標準化后的數據出發，選定標準化歐氏距離（seuclidean）作為樣品相似性的測度，采用離差平方和（ward）法作為類與類之間的距離，應用MATLAB編程，作系統聚類分析。結果表明，30個樣品明顯聚為4類（國外的4個樣品歸入到國內的產地類別）。其中，8個云南主產地訶子樣品和2個緬甸樣品、8個廣西訶子樣品分別聚為一類，廣東、西藏和印度訶子樣品也分別聚為一類，但樣品之間分類有交叉，其中有兩個廣東樣品，落入西藏類別。因此擬采用分類準確度高的人工神經網絡建立產地鑒別模型。

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

主要網絡參數如下：網絡結構3-30-4，即輸入神經元3個，對應于訶子訓練樣品的前3個主成分，隱層30個神經元，對應訶子訓練樣品個數，輸出層神經元4個，對應訶子樣品的類別，以1.0代表云南訶子樣品，2.0代表廣西訶子樣品，以3.0代表廣東訶子樣品，4.0代表西藏訶子樣品。2個緬甸樣品劃分為云南1.0類別，2個印度樣品劃分為西藏4.0類別。隱層的傳遞函數是徑向基函數（radbas），輸出層為競爭層，競爭傳遞函數compet的功能是找出隱層徑向基傳遞函數中各元素的最大值，并且使與最大值對應類別的神經元輸出為1，其他類別的神經元輸出為0。

3 ?小結

以測定的不同產地訶子紅外光譜為研究對象，采用小波變換技術，均方根誤差RMS作為降噪與壓縮性能評價指標，對紅外光譜進行了消噪和壓縮，小波變換重構后的壓縮數據，既保持了原紅外光譜的特征，又有效地減少了光譜變量。據此結合概率神經網絡建立了不同訶子產地鑒別模型。該網絡模型對產地鑒別預測效果很好，能鑒別訶子主產地或偽劣訶子，正確識別率達到100%。紅外光譜結合適當的化學計量學如人工神經網絡方法，具有制樣簡單方便、測定快捷、譜圖準確、重現性好等優點，適合大量樣品的快速鑒別研究，有一定的可靠性和實用性，從而為訶子分類鑒別、質量控制提供了新的適用方法，也為鑒別其他中草藥提供了思路。

參考文獻：

[1] 中華人民共和國衛生部藥典委員會.中華人民共和國藥典（一部）[M].北京：化學工業出版社，2005.129.

[2] 周同惠.紙色譜和薄層色譜[M].北京：科學出版社，1989.

[3] 王金華，孫芳云，杜冠華.訶子的藥理作用研究進展[J].中國民族醫藥志，2006，12（6）：45-47.

[4] 顏玉貞，謝培山，宋力飛，等.西青果藥材及提取物液相指紋圖譜的應用研究[J].中成藥，2004，26（8）：6066-6071.

[5] 周 寧，彭富全.中藏藥中訶子、毛訶子和余甘子的本草考證[J]. 中草藥，2001，32（4）：355-357.

[6] 陳建南，徐鴻華.訶子的產地和品種考證[J].中藥材，1996，19（10）：533-535.

[7] 丁 ?崗，陸蘊如，冀春茹，等.訶子的薄層色譜鑒別[J].北京中醫藥大學學報，2001，24（1）：45-46.

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

[9] 尼 ?珍，胡昌勤，馮 ?芳.近紅外光譜分析中光譜預處理方法的作用及其發展[J].藥物分析雜志，2008，28（5）：824-829.

[10] 李丹婷，張長江，汪 ?勁.基于小波變換的紅外光譜分析[J].光譜學與光譜分析，2006，26（11）：2024-2026.

[11] 李夢龍，康 彬，戚華溢.紅外光譜小波壓縮性能的評價[J].高等學?；瘜W學報，2002，23（7）：1281-1284.

[12] 曹紅翠，吳啟勛.中藥甘草中微量元素的主成分分析和系統聚類分析[J].甘肅聯合大學學報（自然科學版），2010，24（2）：61-63.endprint