[發明專利]一種光譜分析方法

有效

• 申請號: 201010106805.6
• 申請日: 2010-02-03
• 公開（公告）號: CN101788507A
• 公開（公告）日: 2010-07-28
• 發明（設計）人: 李杰;徐寧;周俊武;趙建軍;繆天宇;高揚;卞寧
• 申請（專利權）人: 北京礦冶研究總院
• 主分類號: G01N23/223
• 分類號: G01N23/223;G06F17/00
• 代理公司: 北京凱特來知識產權代理有限公司 11260
• 代理人: 鄭立明;孟麗娟

1.一種光譜分析方法，其特征在于，包括：

分析標樣光譜，采用多項式擬合法和最小二乘法確定譜峰波長和峰位之 間的關系，以及譜峰波長和譜峰寬度之間的關系；

利用確定的所述譜峰波長和峰位之間的關系及譜峰波長和譜峰寬度之間 的關系，計算待測樣品擬包含的元素的激發光譜中可能存在的所有譜峰的峰 位以及譜峰寬度；

在待測樣品全譜范圍內進行尋峰，將尋峰得到的峰位與所述計算的峰位 進行比較，以及將尋峰得到的譜峰寬度與所述計算的譜峰寬度比較，判斷是 否匹配，如果均匹配，則確定不存在完全重疊峰，根據相鄰的譜峰峰位和所 述計算的譜峰寬度確定重疊區域寬度；如果存在不匹配的情況，則確定存在 完全重疊峰，根據完全重疊峰峰位和寬度確定重疊區域寬度；

在各重疊區域內進行譜峰分解，統計各重疊區域內存在的所有所述計算 的峰位和譜峰寬度，得到高斯函數均值初始值和方差初始值，將若干高斯函 數與本底函數相加得到原始譜峰，采用最小二乘法，得出高斯函數系數初始 值及本底函數的數學表達式，所述本底函數為二項式函數；

利用得到的高斯函數和本底函數進行光譜分析；

所述在待測樣品全譜范圍內進行尋峰包括：

執行第一次尋峰，在待測樣品光譜數組Spe[i]中找出滿足如下條件的 Spe[i]對應的i，組成初始峰位數組Temp_Peak[i]，i為光譜道址，取值范圍為從0 到chan，chan為正整數，包括滿足：大于所有峰位下限的加權平均值以及大于 相鄰的前后峰位的峰值；

執行第二次尋峰，在所述初始峰位數組Temp_Peak[i]中找出滿足如下條件的 Spe[i]對應的i作為鎖定的峰位，包括滿足：（Spe[Temp_Peak[i]]-(Spe[Temp_Peak[i]- Ini_PW]+Spe[Temp_Peak[i]+Ini_PW])/2)/Ini_PW>1.5，以及滿足(Temp_Peak[i+1]- Temp_Peak[i])<Ini_PW并且(Temp_Peak[i-1]-Temp_Peak[i])>-Ini_PW，其中Ini_PW為譜峰寬 度。

2.如權利要求1所述的方法，其特征在于，所述采用多項式擬合法和最 小二乘法確定譜峰波長和峰位之間的關系包括：

利用公式wavelength[j]=a1*peakpos_cal[j]²+b1*peakpos_cal[j]+c1，用最小二乘 法計算a1、b1和c1，其中wavelength[j]為標樣光譜波長，peakpos_cal[j]為標樣光 譜峰位。

3.如權利要求1所述的方法，其特征在于，所述采用多項式擬合法和最 小二乘法確定譜峰波長和譜峰寬度之間的關系包括：

利用公式wavewidth[j]=a2*peakpos_cal[j]²+b2*peakpos_cal[j]+c²，用最小二乘 法計算a2、b2和c2，其中wavewidth[j]為標樣光譜譜峰寬度，peakpos_cal[j]為標 樣光譜峰位。

4.如權利要求1所述的方法，其特征在于，所述Ini_PW初始值為chan/100。

5.如權利要求1所述的方法，其特征在于，所述均匹配包括：

尋峰得到的峰位PeakPos[i]與所述計算的峰位peakpos_lib[i]差值的最小值小 于chan/1000+1；chan為正整數；以及

尋峰得到的譜峰寬度PeakWidth[i]與所述計算的譜峰寬度peakwidth_cal[i]差 值的最小值不小于chan/100。