專利名稱：一種色譜峰匹配的分析方法
技術(shù)領(lǐng)域：
本發(fā)明屬于色譜峰的分析方法，特別中藥色譜指紋圖譜峰匹配的分析方法， 尤其屬于用計(jì)算機(jī)程序通過(guò)以目標(biāo)色譜峰對(duì)另一指紋圖譜的漸進(jìn)窗口正交投影 分析來(lái)實(shí)現(xiàn)不同指紋圖譜間色譜峰的匹配的分析方法。
技術(shù)背景色譜圖譜特別是中藥指紋圖譜就其本質(zhì)而言，應(yīng)可視為中藥的一種依賴于 不同提取方法所得的活性化學(xué)組分的相對(duì)濃度譜。由于中藥色譜指紋圖譜所含 化學(xué)組分較多，且很多化學(xué)組分還缺乏相應(yīng)的標(biāo)準(zhǔn)對(duì)照品，難于采用傳統(tǒng)的分 析方法來(lái)精確描述。如果色譜實(shí)驗(yàn)條件控制得好，同一中藥樣本的色譜指紋圖 譜的重復(fù)性應(yīng)是相對(duì)穩(wěn)定的。但是由于實(shí)驗(yàn)的時(shí)間、地點(diǎn)以及一些不可控實(shí)驗(yàn) 因素的差異，往往也會(huì)在一定程度上導(dǎo)致色譜保留時(shí)間的偏移，給中藥色譜指 紋圖譜的識(shí)別與鑒定帶來(lái)不便，如不作必要的儀器系統(tǒng)誤差的校準(zhǔn)，還可能導(dǎo) 致錯(cuò)誤的結(jié)論。盡管在色譜指紋圖譜建立過(guò)程中，是以最大信息量為目標(biāo)盡量 保證峰與峰的分離，但每一色譜圖譜中所含的化學(xué)成分較復(fù)雜，色譜峰間有時(shí) 不可避免的發(fā)生重疊、缺失或漂移，依據(jù)保留時(shí)間或相對(duì)保留時(shí)間無(wú)法進(jìn)行不 同圖譜間的色譜峰匹配。另外，中藥色譜指紋圖譜成為我國(guó)對(duì)中藥注射液進(jìn)行 質(zhì)量控制的手段，各地藥品檢驗(yàn)所對(duì)在中藥注射液進(jìn)行復(fù)核時(shí)，盡管可完全按 標(biāo)準(zhǔn)規(guī)定的分析實(shí)驗(yàn)條件進(jìn)行檢測(cè)，但是，不可能使用同一臺(tái)儀器、同一色譜 柱，從而將導(dǎo)致所得的中藥色譜圖譜發(fā)生一定的波動(dòng)和色譜保留時(shí)間的偏移， 這將給中藥色譜指紋圖譜的識(shí)別與鑒定帶來(lái)困難。目前對(duì)中藥色譜指紋圖譜組分保留時(shí)間漂移進(jìn)行校準(zhǔn)主要采用局部最小二乘校正法(李博巖，梁逸曾，胡蕓,等.分析化學(xué),2004,32(3): 313-316.)等?，F(xiàn)有方 法沒(méi)能充分利用聯(lián)用色譜儀器在色譜分離過(guò)程中所記錄的各流出組分的波譜信 息，不易進(jìn)行程序化處理并使之與HPLC儀相匹配。 發(fā)明內(nèi)容本發(fā)明所要解決的技術(shù)問(wèn)題在于提供一種聯(lián)用色譜中的組分光譜信息，通 過(guò)設(shè)有程序的數(shù)字計(jì)算機(jī)來(lái)準(zhǔn)確地進(jìn)行不同指紋圖譜間色譜峰的匹配分析。本發(fā)明解決一技術(shù)問(wèn)題的技術(shù)方案為 一種色譜峰匹配的分析方法，包括含DAD檢測(cè)器的高效液相色譜儀對(duì)物質(zhì)進(jìn)行檢測(cè)，存儲(chǔ)程序的數(shù)字計(jì)算機(jī)，該存 儲(chǔ)程序的數(shù)字計(jì)算機(jī)執(zhí)行如下步驟a)在圖譜2中提取待匹配的色譜峰，扣除色 譜峰中的噪音信號(hào)，將扣除噪音信號(hào)的待匹配的色譜峰進(jìn)行奇異值分解形成抽 象光譜2，提取抽象光譜2中的第一列，以此形成矩陣A;b) 在待鑒圖譜l中尋找與圖譜2中待匹配的色譜峰相匹配的目標(biāo)色譜峰， 設(shè)置移動(dòng)窗口 w ，在待鑒圖譜l中自第1點(diǎn)起移動(dòng)窗口，同時(shí)對(duì)待鑒圖譜l 進(jìn)行奇異值分解，并將抽象光譜中的第1至窗口尺寸w的列數(shù)提取出來(lái)，以其構(gòu)造正交投影陣M;c) 將矩陣A對(duì)M進(jìn)行投影運(yùn)算，得投影后的殘余光譜；d) 計(jì)算待識(shí)別色譜峰A與殘余光譜的夾角，判斷夾角余弦是否小于噪音信號(hào)， 若夾角余弦小于噪音信號(hào)，則說(shuō)明在待鑒圖譜1中在夾角余弦小于噪音信號(hào)的 時(shí)間段內(nèi)的色譜峰的主成分與圖譜2中的待匹配的色譜峰為相同的物質(zhì)，此峰 相匹配，否則，待鑒圖譜1中沒(méi)有與圖譜2中待匹配的色譜峰相匹配的色譜峰。所述的W為大于等于1的整數(shù)。 所述的噪音信號(hào)的值為0.05。 本發(fā)明用于決明子色譜指紋圖譜峰的匹配分析。本發(fā)明用于莪術(shù)油色譜指紋圖譜峰的匹配的分析。本發(fā)明的原理為設(shè)X為采用聯(lián)用色譜測(cè)得的樣本2的中藥色譜指紋圖譜的N XM階二維數(shù)據(jù)矩陣，它們的行方向表征了色譜的流出信息(M個(gè)流出時(shí)間點(diǎn))，列方向表征了光譜信息(N個(gè)吸收波長(zhǎng))；設(shè)Y為采用聯(lián)用色譜測(cè)得的 樣本S l的色譜指紋圖譜的iV'xM'階二維數(shù)據(jù)矩陣，它們的行方向表征了色譜的流出信息(M'個(gè)流出時(shí)間點(diǎn))，列方向表征了光譜信息(iv'個(gè)吸收波長(zhǎng))。首先確定樣本2中待匹配的基準(zhǔn)色譜峰a扃=X ( [ t！ : t2 ]， ) (1) t,、 t2為基準(zhǔn)色譜峰的起止時(shí)間。再將&1進(jìn)行奇異值分解，并將抽象光譜中的第一列即主成分提取出來(lái)[Si， vb di ] = svd (ai ) (2) A=山(，1 ) (3) 再對(duì)樣本1進(jìn)行處理，尋找與樣本2中基準(zhǔn)色譜峰相匹配的目標(biāo)色譜峰。 設(shè)置移動(dòng)窗口 w ,在Y中自第1點(diǎn)起移動(dòng)至n2 -w +1止(n2為Y的列數(shù))， 同時(shí)對(duì)Y進(jìn)行奇異值分解，并將抽象光譜中的第1至窗口尺寸w的列數(shù)提取 出來(lái)，以其構(gòu)造正交投影陣。[s2， v2， d2 ] = svd ( Y ( i : i + w - 1， ) ) (4) M = I — d2 d2+ (5) 其中，I為單位矩陣，d2+為d2的廣義逆陣。用矩陣M對(duì)A進(jìn)行投影運(yùn)算 a* = M A (6) a*表示A投影后的殘余光譜。為消除不等性噪聲的影響，令投影后所得 的殘余光譜a^與投影前的原始光譜A的夾角余弦為rj ，即并以rj作為鑒別色譜峰峰匹配的判據(jù)。由于M扣除了 A與d2相關(guān)的光譜信息，a*保留了 A與d2不相關(guān)的光 譜信息；所以，r」值越小，即經(jīng)投影所得的殘余光譜a*與投影前的原始光譜A 的夾角越大，說(shuō)明3*保留A的光譜信息越少，A與d2相關(guān)性越大；反之巧值 越大，說(shuō)明a*保留A的光譜信息越多，A與d2的光譜相關(guān)性越小。當(dāng)巧值 等于0時(shí)，a*與A完全正交，M完全扣除了A與d2相關(guān)的光譜信息，a*完 全只剩下隨機(jī)噪聲。以樣本2中待匹配的色譜峰的抽象光譜矩陣A對(duì)投影矩陣M進(jìn)行投影，可 得光譜殘差?yuàn)A角余弦曲線。隨著窗口的移動(dòng)，若該曲線的r」值為零或接近于零， 則說(shuō)明樣本1中此段時(shí)間內(nèi)的色譜峰的主成分與樣本2中的目標(biāo)色譜峰為相 同的物質(zhì)，此峰相匹配，否則，樣本1中沒(méi)有與樣本2中目標(biāo)色譜峰相匹配 的色譜峰。指紋圖譜通常情況下都是進(jìn)行主成分間的峰匹配，當(dāng)采用本發(fā)明的分析法 對(duì)不同指紋圖譜間的色譜峰進(jìn)行匹配時(shí)，通過(guò)窗口的設(shè)置可將抽象光譜中的第1 至窗口尺寸w的列數(shù)提取出來(lái)，依次進(jìn)行峰匹配。如某一色譜峰為由三種物質(zhì) 組成的混合峰，如設(shè)窗口為l，則進(jìn)行的是主成分間的峰匹配；窗口設(shè)為2，則 進(jìn)行的是次成分間的峰匹配；窗口設(shè)為3，則進(jìn)行的是第三種成分間的峰匹配。本發(fā)明與現(xiàn)有技術(shù)相比，克服了僅單憑單波長(zhǎng)數(shù)據(jù)進(jìn)行色譜峰匹配的盲目 性；無(wú)需驗(yàn)證峰純度，光譜識(shí)別率較好，可有效區(qū)分同系物等具有近似光譜及 末端吸收的成分色譜峰；可處理不復(fù)雜的重疊色譜峰問(wèn)題，既能夠快速處理數(shù) 據(jù)，還能防止數(shù)據(jù)的漏選。


圖1為本發(fā)明的計(jì)算機(jī)流程圖。圖2為決明子空白溶液的Agilent 1100高效液相色譜儀色譜柱1的HPLC指 紋圖譜。圖3為標(biāo)準(zhǔn)品溶液的Agilent 1100高效液相色譜儀色譜柱1的HPLC指紋圖譜。圖4為大決明子生品的Agilent 1100高效液相色譜儀色譜柱1的HPLC指紋圖譜。圖5為小決明子生品的HP1100高效液相色譜儀色譜柱2的HPLC指紋圖譜。 圖6為大黃酚標(biāo)準(zhǔn)品與大決明子指紋圖譜全部二維數(shù)據(jù)正交投影前后光譜 矢量夾角余弦圖。圖7為大黃素標(biāo)準(zhǔn)品與大決明子指紋圖譜全部二維數(shù)據(jù)正交投影前后光譜 矢量夾角余弦圖。圖8為大決明子的405號(hào)峰與小決明子指紋圖譜全部二維數(shù)據(jù)正交投影前 后光譜矢量夾角余弦圖。圖9為大決明子的415號(hào)峰與小決明子指紋圖譜全部二維數(shù)據(jù)正交投影前 后光譜矢量夾角余弦圖。圖IO為吉馬酮標(biāo)準(zhǔn)品的HPLC圖譜。圖11為蓬莪術(shù)揮發(fā)油的HPLC指紋圖譜。圖12為廣西莪術(shù)揮發(fā)油的HPLC指紋圖譜。圖13為溫郁金揮發(fā)油的HPLC指紋圖譜。圖14為吉馬酮標(biāo)準(zhǔn)品與溫郁金揮發(fā)油指紋圖譜全部二維數(shù)據(jù)正交投影前后光譜矢量夾角余弦圖。圖15為廣西莪術(shù)的10號(hào)峰與蓬莪術(shù)指紋圖譜全部二維數(shù)據(jù)正交投影前后光譜矢量夾角余弦圖。圖16為溫郁金的C6號(hào)峰與廣西莪術(shù)指紋圖譜全部二維數(shù)據(jù)正交投影前后 光譜矢量夾角余弦圖。圖17為蓬莪術(shù)的A6號(hào)峰與溫郁金指紋圖譜全部二維數(shù)據(jù)正交投影前后光譜 矢量夾角余弦圖在圖3中，31為蘆薈大黃素，32為大黃酸，33為大黃素，34為大黃酚，35為 大黃素甲醚；在圖4中，401—418為大決明子生品中待匹配的色譜峰；在圖5中， l一26為小決明子生品中待匹配的色譜峰；在圖11中，A1—A12為蓬莪術(shù)揮發(fā)油 的色譜吸收峰；在圖12中，Bl-B13為廣西莪術(shù)揮發(fā)油的色譜吸收峰；在圖13中， Cl-C9為溫郁金揮發(fā)油的色譜吸收峰；具體實(shí)施方式
下面結(jié)合附圖詳細(xì)說(shuō)明本發(fā)明的優(yōu)選實(shí)施例。由圖1所示，本發(fā)明所述的色譜峰匹配的分析方法，包括含DAD檢測(cè)器的高 效液相色譜儀對(duì)物質(zhì)進(jìn)行檢測(cè)，形成了圖譜的全部信息，利用MATLAB語(yǔ)言進(jìn) 行處理，得到圖譜的全部信息，所述的信息包括相應(yīng)峰的峰面積，峰高、保留 時(shí)間，將圖譜的全部信息輸入到存儲(chǔ)程序的數(shù)字計(jì)算機(jī)，該存儲(chǔ)程序的數(shù)字計(jì)算機(jī)執(zhí)行如下步驟步驟201—207為在圖譜2中提取待匹配的色譜峰，扣除待匹配色譜峰中的 噪音信號(hào)，將扣除噪音信號(hào)的待匹配色譜峰進(jìn)行奇異值分解形成抽象光譜2，提 取抽象光譜2中的第一列，以此形成矩陣A;步驟101 —102為選取欲與圖譜2中提取的待匹配色譜峰進(jìn)行峰匹配的圖 譜，作為圖譜l，扣除圖譜l中的噪音信號(hào)，避免不同儀器或者同一儀器在不同 時(shí)間所產(chǎn)生的漂移對(duì)檢測(cè)結(jié)果的影響，常用的噪音信號(hào)為0.05。步驟103 — 107為在扣除噪音信號(hào)的圖譜l中設(shè)置標(biāo)準(zhǔn)窗口，標(biāo)準(zhǔn)窗口的大小為W，在Y中自第1點(diǎn)起移動(dòng)至n2 —w+l止(n2為Y的列數(shù))，同時(shí)對(duì)Y進(jìn) 行奇異值分解，形成抽象光譜l，并將抽象光譜l中的第1至窗口尺寸w的列 數(shù)提取出來(lái)，以其構(gòu)造正交投影陣。步驟108 — 109為將矩陣A對(duì)標(biāo)準(zhǔn)圖譜中的正交投影矩陣進(jìn)行投影，扣除所 取窗口處的光譜數(shù)據(jù)，得殘余光譜，計(jì)算待識(shí)別色譜峰與殘余光譜的夾角。步驟110為判斷夾角余弦是否小于噪音信號(hào)，若夾角余弦小于噪音信號(hào)， 則為步驟lll，說(shuō)明在待測(cè)圖譜1中在夾角余弦小于噪音信號(hào)的時(shí)間段內(nèi)的色 譜峰的主成分與圖譜2中的待匹配色譜峰為相同的物質(zhì)，此峰相匹配；若夾角 余弦大于噪音信號(hào)，則為步驟U2，待測(cè)圖譜1中沒(méi)有與圖譜2中待匹配色譜 峰相匹配的色譜峰。將本發(fā)明用于決明子色譜指紋圖譜峰的匹配分析。1.儀器與試藥1.1儀器Agilent 1100高效液相色譜儀(包括AgilentllOO四元泵，Agilent 1100 二極管 陣列檢測(cè)器，Agilent 1100色譜化學(xué)工作站)；HP1100高效液相色譜儀(包括 AgilentllOO四元泵，HP1100 二極管陣列檢測(cè)器，HP色譜化學(xué)工作站)。1.2樣品與試劑大決明子藥材，產(chǎn)地為河北，小決明子藥材，產(chǎn)地為安徽，經(jīng)佳木斯大學(xué) 生藥學(xué)教研室劉娟教授鑒定分別為豆科植物決明(Ca肌'"L.)和小 決明(C似W"tomL.)的成熟種子；蘆薈大黃素對(duì)照品(批號(hào)766-200210)、大 黃素對(duì)照品(批號(hào)756-200110)、大黃酸對(duì)照品(批號(hào)774-200311)、大黃酚對(duì) 照品(批號(hào)796-200204)、大黃素甲醚對(duì)照品(批號(hào)758-2004080)均購(gòu)自中國(guó)藥 品生物制品檢定所。乙腈、甲醇為色譜純，鹽酸、磷酸及氯仿為分析純。2.方法與結(jié)果 2.1供試品制備2丄1對(duì)照品溶液稱取蘆薈大黃素、大黃酸、大黃素、大黃酚和大黃素 甲醚對(duì)照品適量，以甲醇為溶劑，配制成200 ng'mL"的蘆薈大黃素、大黃酸、 大黃素、大黃酚和大黃素甲醚五種成分的混合對(duì)照品溶液。2丄2供試品溶液取決明子樣品，粉碎，過(guò)80目篩，稱取2.0g，力n 1.5 moH/1的鹽酸溶液30 mL，加熱回流1.5h，放冷，加氯仿25 mL，回流20min, 取出，分取氯仿層，酸水層再加氯仿25 mL，回流20min，分取氯仿層，同上 法操作，加氯仿再回流提取二次，合并四次氯仿提取液，加水20mL振搖洗脫， 棄去水層，將氯仿液蒸干，殘?jiān)蛹状?5mL使溶解，過(guò)濾，定量轉(zhuǎn)移至25mL 容量瓶中，加甲醇至刻度，搖勻，經(jīng)0.45pm微孔濾膜過(guò)濾，即得。2.2色譜條件色譜柱l: Zorbax Eclipse XDB - C18(4.6 mm x 250 mm ， 5 pm )，保護(hù)柱 Zorbax SB - <:18柱(4,6 mm x 12.5 mm ， 5 );色譜柱2: Hypersil ODS (4,0 mm x 250 mm ， 5 ,)，保護(hù)柱Scienhome (28柱(3.0 mm x 20 mm ， 5 jam );流動(dòng) 相乙腈(A ) -0.1%磷酸水溶液(B)作梯度洗脫，Omin， 10%A-90% B; 5 min， 15 % A - 85 % B; 10 min， 40 % A - 60 % B; 40 min， 70 % A - 30 % B; 41min后100。/。A。流速l.OmL'min";柱溫室溫；檢測(cè)波長(zhǎng)278 nm， 參比波長(zhǎng)550 nm;光譜數(shù)據(jù)記錄范圍220 500 nm，波長(zhǎng)間隔2nm;色譜記錄 時(shí)間0 50min;進(jìn)樣量IO jaL。2.3色譜圖的確定2.3.1特征峰選定按3.1.2供試品溶液的制備方法做空白試驗(yàn)，結(jié)果見(jiàn)圖 2。檢出峰中峰面積最大值及峰高最大值分別在200A及5毫吸收度以下，因此10選取AreaRejec＝200， Height Reject=5對(duì)所有供試品進(jìn)行積分，選取有積分值的 峰作為決明子的特征峰。2.3.2取蘆薈大黃素、大黃酸、大黃素、大黃酚和大黃素甲醚對(duì)照品混合溶 液于Agilent 1100高效液相色譜儀、色譜柱1上按選定色譜條件同法測(cè)定，結(jié) 果見(jiàn)圖3。2.3.3另取大決明子、小決明子供試品溶液分別于Agilent 1100高效液相色 譜儀、色譜柱1及HP 1100高效液相色譜儀、色譜柱2上按選定色譜條件測(cè)定， 記錄其色譜圖，按規(guī)定方法積分，結(jié)果見(jiàn)圖4 圖5。2.4色譜峰的匹配對(duì)不同條件下的標(biāo)準(zhǔn)品溶液、大、小決明子供試品溶液的色譜指紋圖譜間 的峰進(jìn)行匹配分析。如大黃酚標(biāo)準(zhǔn)品與大決明子指紋圖譜間的峰匹配，見(jiàn)圖6; 大黃素標(biāo)準(zhǔn)品與大決明子指紋圖譜間的峰匹配，見(jiàn)圖7;大決明子的5號(hào)峰與小決明子指紋圖譜間的峰匹配，見(jiàn)圖8;大決明子的15號(hào)峰與小決明子指紋圖譜間的峰匹配，見(jiàn)圖9。由圖6可知，大黃酚標(biāo)準(zhǔn)品的色譜峰與大決明子指紋圖譜中時(shí)間點(diǎn)在 40.3955 40.0488的色譜峰即第415號(hào)色譜峰相匹配；由圖7可知，在大決明子 的色譜指紋圖譜中沒(méi)有與大黃素標(biāo)準(zhǔn)品色譜峰相匹配的色譜峰；由圖8可知， 大決明子的第405號(hào)色譜峰與小決明子指紋圖譜中時(shí)間點(diǎn)在14.9471 14.9956 的色譜峰即第7號(hào)色譜峰相匹配；由圖9可知，大決明子的第415號(hào)色譜峰與 小決明子指紋圖譜中時(shí)間點(diǎn)在40.9990 41.3033的色譜峰即第22號(hào)色譜峰相匹 配。同理可知，大決明子的401、 402、 403、 404、 405、 406、 407、 408、 409、 410、 411、 412、 413、 414、 415、 416、 417、 418號(hào)峰分別與小決明子的1、 2、3、 4、 7、 9、 10、 13、 14、 15、 17、 20、 21、 22、 23、 24、 25、 26號(hào)峰相匹配; 大決明子的412、 415、 416號(hào)峰，小決明子的19、 22、 23號(hào)峰分別為大黃素、 大黃酚和大黃素甲醚。本發(fā)明用于莪術(shù)油色譜指紋圖譜峰的匹配的分析。3.實(shí)驗(yàn)部分3.1儀器與試劑Agilent 1100高效液相色譜儀(包括Agilent 1100四元泵，Agilent 1100 二極 管陣列檢測(cè)器，Agilent 1100色譜化學(xué)工作站)。蓬莪術(shù)、廣西莪術(shù)和溫郁金購(gòu)于 北京興中中藥飲片廠，產(chǎn)地浙江，經(jīng)佳木斯大學(xué)生藥教研室史權(quán)教授鑒定，符 合中國(guó)藥典(2005版) 一部項(xiàng)下的規(guī)定。吉馬酮對(duì)照品(批號(hào)111665-200401)購(gòu) 自中國(guó)藥品生物制品檢定所。甲醇為色譜純，無(wú)水硫酸鈉、磷酸為分析純。3.2供試品制備3.2.1揮發(fā)油制備取經(jīng)粉碎過(guò)30目篩的蓬莪術(shù)、廣西莪術(shù)、溫郁金粗粉 各100g，分別置于2000mL圓底燒瓶中，加8倍量的水浸泡過(guò)夜，按中國(guó)藥典 (2005版) 一部附錄X D中揮發(fā)油測(cè)定的甲法，水蒸氣蒸餾提取8h。分別收 集揮發(fā)油部分，加入無(wú)水硫酸鈉脫水后，密封、暗處冷藏。3.2.2供試液制備取蓬莪術(shù)、廣西莪術(shù)、溫郁金揮發(fā)油各20iiL，分別置 10mL量瓶中，加入甲醇至刻度，搖勻。經(jīng)0.45pm微孔濾膜過(guò)濾，制成蓬莪術(shù)、 廣西莪術(shù)、溫郁金揮發(fā)油的供試液。取吉馬酮對(duì)照品約10 mg，精密稱定，置 100mL容量瓶中，加入甲醇至刻度，搖勻，經(jīng)0.45 pm微孔濾膜過(guò)濾，吸取濾 液lmL置10mL容量瓶中，加入甲醇至刻度，配制吉馬酮對(duì)照品溶液。3.3色譜條件色譜柱Zorbax Eclipse XDB - C18(4.6mmx250mm， 5jim)，保護(hù)柱ZorbaxSB-ds柱(4.6mmxl2.5mm， 5pm);流動(dòng)相甲醇(A) -0,2%磷酸水溶液(B) 作梯度洗脫，0min， 55%A-45%B; 5 min， 55%A-45%B; 25 min， 75% A-25%B; 40min， 75%A-25%B; 50min, 100 %A。流速l.Oml.min—、 柱溫室溫；檢測(cè)波長(zhǎng)210nm;光譜數(shù)據(jù)記錄范圍190 350 nm，波長(zhǎng)間隔2 nm;色譜記錄時(shí)間0 50min;進(jìn)樣體積20 pL。取吉馬酮對(duì)照品按選定色譜條件測(cè)定，結(jié)果見(jiàn)圖IO。另取蓬莪術(shù)、廣西莪 術(shù)、溫郁金供試液分別按選定色譜條件同法測(cè)定，記錄其色譜圖，結(jié)果見(jiàn)圖11 圖13。4色譜峰的匹配對(duì)吉馬酮標(biāo)準(zhǔn)品、蓬莪術(shù)、廣西莪術(shù)、溫郁金揮發(fā)油色譜指紋圖譜間的峰 進(jìn)行匹配分析。如吉馬酮標(biāo)準(zhǔn)品與溫郁金揮發(fā)油指紋圖譜間的峰匹配，見(jiàn)圖14; 廣西莪術(shù)的IO號(hào)峰與蓬莪術(shù)指紋圖譜間的峰匹配，見(jiàn)圖15;溫郁金的6號(hào)峰與 廣西莪術(shù)指紋圖譜間的峰匹配，見(jiàn)圖16;蓬莪術(shù)的6號(hào)峰與溫郁金指紋圖譜間 的峰匹配。由圖14可知，吉馬酮標(biāo)準(zhǔn)品的色譜峰與溫郁金揮發(fā)油指紋圖譜中時(shí)間點(diǎn)在 28.1160 28.1360的色譜峰即第C6號(hào)色譜峰相匹配；由圖15可知，在蓬莪術(shù)的 色譜指紋圖譜中沒(méi)有與廣西莪術(shù)的B10號(hào)色譜峰相匹配的色譜峰；由圖16可知， 溫郁金的第C6號(hào)色譜峰與廣西莪術(shù)指紋圖譜中時(shí)間點(diǎn)在24.5048 24.5248的色 譜峰即第B9號(hào)色譜峰相匹配；由圖17可知，在溫郁金指紋圖譜中沒(méi)有與蓬莪 術(shù)的第A6號(hào)色譜峰相匹配的色譜峰。
權(quán)利要求
1、一種色譜峰匹配的分析方法，包括含DAD檢測(cè)器的高效液相色譜儀對(duì)物質(zhì)進(jìn)行檢測(cè)，將檢測(cè)結(jié)果輸入到存儲(chǔ)程序的數(shù)字計(jì)算機(jī)，其特征在于該存儲(chǔ)程序的數(shù)字計(jì)算機(jī)執(zhí)行如下步驟a)在圖譜2中提取待匹配的色譜峰，扣除色譜峰中的噪音信號(hào)，將扣除噪音信號(hào)的待匹配的色譜峰進(jìn)行奇異值分解形成抽象光譜2，提取抽象光譜2中的第一列，以此形成矩陣A；b)在待鑒圖譜1中尋找與圖譜2中待匹配的色譜峰相匹配的目標(biāo)色譜峰，設(shè)置移動(dòng)窗口w，在待鑒圖譜1中自第1點(diǎn)起移動(dòng)窗口，同時(shí)對(duì)待鑒圖譜1進(jìn)行奇異值分解，并將抽象光譜中的第1至窗口尺寸w的列數(shù)提取出來(lái)，以其構(gòu)造正交投影陣M；c)將矩陣A對(duì)M進(jìn)行投影運(yùn)算，得投影后的殘余光譜；d)計(jì)算待識(shí)別色譜峰A與殘余光譜的夾角，判斷夾角余弦是否小于噪音信號(hào)，若夾角余弦小于噪音信號(hào)，則說(shuō)明在待鑒圖譜1中在夾角余弦小于噪音信號(hào)的時(shí)間段內(nèi)的色譜峰的主成分與圖譜2中的待匹配的色譜峰為相同的物質(zhì)，此峰相匹配，否則，待鑒圖譜1中沒(méi)有與圖譜2中待匹配的色譜峰相匹配的色譜峰。
2、 根據(jù)權(quán)利要求l所述的一種色譜峰匹配的分析方法，其特征在于所述的W為大于等于1的整數(shù)。
3、 根據(jù)權(quán)利要求l所述的一種色譜峰匹配的分析方法，其特征在于所述的噪音信號(hào)的值為O. 05。
4、 權(quán)利要求l所述的一種色譜峰匹配的分析方法用于決明子色譜指紋圖譜 峰的匹配分析。
5、 權(quán)利要求l所述的一種色譜峰匹配的分析方法用于莪術(shù)油色譜指紋圖譜 峰的匹配的分析。
全文摘要
本發(fā)明公開(kāi)一種色譜峰匹配的分析方法，包括含DAD檢測(cè)器的高效液相色譜儀對(duì)物質(zhì)進(jìn)行檢測(cè)，將檢測(cè)結(jié)果輸入存儲(chǔ)程序的數(shù)字計(jì)算機(jī)，計(jì)算機(jī)執(zhí)行如下步驟選取待鑒別的色譜峰，將待鑒別色譜峰形成抽象光譜2，提取抽象光譜2中的第一列形成矩陣A；選取欲峰匹配的圖譜1，設(shè)置窗口大小為W，移動(dòng)窗口形成抽象光譜1，構(gòu)造成正交投影矩陣；將矩陣A對(duì)正交投影矩陣進(jìn)行投影，得殘余光譜；判斷殘余光譜與原始光譜A的夾角余弦是否小于噪音信號(hào)，若小于噪音信號(hào)，此峰相匹配；否則不匹配。本發(fā)明與現(xiàn)有技術(shù)相比，無(wú)需驗(yàn)證峰純度，光譜識(shí)別率較好，可有效區(qū)分同系物等具有近似光譜及末端吸收的成分色譜峰，既能夠快速處理數(shù)據(jù)，還能防止數(shù)據(jù)的漏選。
文檔編號(hào)G01N30/86GK101256176SQ200710190088
公開(kāi)日2008年9月3日 申請(qǐng)日期2007年11月21日 優(yōu)先權(quán)日2007年11月21日
發(fā)明者鄒純才, 鄢海燕 申請(qǐng)人:鄒純才