CN1143188A - 在光谱测定法中稳定光谱的方法 - Google Patents

Abstract

一种对被测目标物中的一种特定组分作定量分析的方法，包括下列步骤：测量穿过目标物的透射光能谱或来自目标物的反射光能谱，把能谱分成多个波长区间，由此得到多个局部的能谱，在每个波长区间内的多个局部能谱用包含在每个波长区间内的一个预定波长的实测能值归一化，用已归一化的多个局部能谱通过多变量分析进行定量分析。

Description

在光谱测定法中稳定光谱的方法

本发明是关于在光谱测定法中稳定光谱的方法，特别是关于使测量稳定的方法，通过在数字计算中减少由于分光镜内温度、电源电压等变化所引起的光谱基线的波动和由于加到被测物上的光的入射角变化所导致的透射率的变化。

通常，可用一种分析曲线对一种特定成分的浓度进行定量分析，该分析曲线是预先通过测量透过已知浓度的样品所透过的光的能谱或被已知其浓度的样品反射的光的能谱而得到的.普遍公知的获得分析曲线的方法有采用单波长的方法和采用双波长或多波长的方法。

然而，在光谱中光谱基线的波动是由于分光镜内部和外部的温度变化、电源电压的变化等原因所引起的，这些基线的波动使定量分析的准确性降低，因此必须消除光谱中的误差成分。

作为用数字计算减少基线波动的方法，通常公知的有：双波长基线补偿、数字微分、低频滤波富里叶变换。下面对这些方法中的每一种都作简单的描述。

双波长基线补偿：

双波长基线补偿法法利用两个波长λ₁和λ₂使基线部分消失；溶液的任何组分对该波长的吸收率是常数。准确地说，这两个波长λ1和λ₂分别满足下面的方程(1)和(2)。

其中，α_i，λ1：对波长λ₁的吸收率，

c_i：i组分浓度，n：溶液中的组分数目，I：容器长。

其中，α_i，λ2：对波长λ₂的吸收率，

c_i：i组分浓度，

n：溶液中的组分数目，

I：容器长。

因此，如果在波长λ₁和λ₂处所测得的吸收谱f(λ)的吸收率分别是const3和const4，那么，通过从f(λ)减去(const4-const3)/(λ₂－λ₁)×(λ－λ₁)＋const3，得到补偿光谱。

数字微分

数字微分法通过关于波长对吸收率的数字微分来减少或消除低频处的偏移。微分的次数通常是一次或二次。通过一次微分可以消除与波长无关的常数项。通过二次微分可以消除关于该波长的一次偏移。数字微分时，吸收率关于波长微分，从而，高于二次的微分会使信号大大失真。因此，极少使用高于三次的数字微分。

低频滤波富里叶变换

低频滤波富里叶变换法把光谱经富里叶变换、也就是把能谱从时域转换到频域，滤去引起基线偏移的低频部分。

双波长补偿法要求用于补偿的两个波长λ₂和λ₁分别满足上面的方程(1)和方程(2)。然而，在组分多于两种的多组分溶液中，常常并不存在能满足上面方程(1)(2)的两个波长。在这种情况下，不可能用双波长补偿法。

数字微分法由于数字微分使诸如随机杂波之类的高频杂波被放大，因此，基线波动的带宽必定宽于被物体吸收的光的带宽。此外，由于数字微分而使衍生谱畸变。

在低频滤波富里叶变换法中，存在着一个问题，即由于与基线波动部分一起过滤而使被物体吸收之光的信息在低频区域内丢失。

此外，数字微分和低频滤波富里叶变换存在着一个共同的问题，这个问题是，如果能量E是时间t和波长λ这二个变量的函数，那么E(λ，t)的全微分由下面的方程(3)表示：dE＝(E/λ)_tdλ＋(E/t)_λdt.

                                     (3)

已有的数字微分和低频滤波富里叶变换法两者都是对波长λ作用，因此，不可能使作为t的函数的方程(3)第二项(E/t)_λ被减弱或消去。

因此，本发明的目的是提供在光谱测定法中稳定光谱的方法，以减少与时间有关的光谱误差部分，同时也减少由于加到被测物上的光的入射角变化所导致的透射率的变化。

本发明是基于下面的考虑，由于诸如分光镜温度的变化、电源电压的变化等扰动所导致的能谱的误差通常表现为光谱线的漂移，也就是光谱的瞬时变化。因此，在波长λ和时间t下的透射光强度I(λ，t)和波长λ_r和时间t下的吸收光强度可分别由下面的方程(4)和(5)表示： $I(\mathrm{\λ},t)=k(\mathrm{\λ},t)\·I_0(\mathrm{\λ},t_0)\·p(\mathrm{\λ},\mathrm{\θ})\exp (-\underset{i}{\mathrm{\Σ}}{a}_i\left(\mathrm{\λ}\right)c_{i}l),\left(4\right)$ $I({\mathrm{\λ}}_r,t)=k({\mathrm{\λ}}_r,t)\·I_0({\mathrm{\λ}}_r,t_0)\underset{}{\·p({\mathrm{\λ}}_r,\mathrm{\θ})\·\exp \[-\underset{i}{\mathrm{\Σ}}{a}_i\left({\mathrm{\λ}}_r\right)c_{i}l\],\left(5\right)}$ 其中，I(λ，t)：波长λ和时间t下的透射光强度，

I(λ_r，t)：波长λ_r和时间t下的透射光强度，

I₀(λ，t₀)：波长λ和时间t₀下的入射光强度，

I₀(λ_r，t₀)：波长λ_r和时间t₀下的入射光强度，

P(λ，θ)：波长λ、入射角为θ的光的表面透射率，

P(λ_r，θ)：波长λ_r、入射角为θ的光的表面透射率，

α_i(λ)：组分i对波长λ的吸收系数，

α_i(λ_r)：组分i对波长λ_r的吸收系数，

θ：光与被测物体表面的入射角，

C_i：组分i的浓度，和

1：有效光程长。

此外，方程(4)中k(λ，t)是波长λ和时间t下的入射光强度与波长λ和时间t₀下的入射光强度I₀(λ，t₀)的比值。因此，

      k(λ，t₀)＝1    (6)

同样，方程(5)中k(λ_r，t)是波长λr和时间t下的入射光强度与波长λ_r和时间t₀下的入射光强度I₀(λ_r，t₀)的比值。因此，

     k(λ_r，t₀)＝1    (7)

方程(4)与方程(5)的比使我们得出下述方程(8)所表示的能量强度比I_N： $I_N=\frac{k(\mathrm{\λ},t)}{k({\mathrm{\λ}}_r,t)}\·\frac{I_0(\mathrm{\λ},t_0)}{I_0({\mathrm{\λ}}_r,t_0)}\·\frac{p(\mathrm{\λ},\mathrm{\θ})}{p({\mathrm{\λ}}_r,\mathrm{\θ})}\·\exp \[-\underset{i}{\mathrm{\Σ}}(a_i\left(\mathrm{\λ}\right)-a_i\left({\mathrm{\λ}}_r\right)c_{i}l)\]\left(8\right)$

这里，我们假定

α_i(λ)≠α_i(λ_r)    (9)

k(λ，t)的全微分由下面的方程(10)表示：

  dk(λ，t)＝(k/λ)_tdλ＋(k/t)_λdt.

                                                    (10)

这里，(k/λ)_t表示入射光强度比值k(λ，t)波长相关性，(k/t)_λ表示入射光的强度比值k(λ，t)关于时间的变化率。

如果在波长λ和时间t下入射光强度k(λ，t)×I₀(λ，t)对波长λ_r和时间t下入射光强度k(λ_r，t)×I₀(λ_r，t₀)的比率用β(λ，λ_r，t)表示，那么： $\mathrm{\β}(\mathrm{\λ},{\mathrm{\λ}}_r,t)=\frac{k(\mathrm{\λ},\mathrm{\τ})}{k({\mathrm{\λ}}_r,t)}\·\frac{I_0(\mathrm{\λ},t_0)}{I_0({\mathrm{\λ}}_r,t_0)}$

现在把波长范围分成几个子区间，把波长λ和λ_r限定在第j个子区间内(λ_j、最小，λ_j、最大)，那么可用下式(12)近似表示比率k(λ，t)波长相关性的项：

  (k/λ)_t＝o，  (12)

因此，

   k(λ，t)≈k(t)    (13)

从而，β(λ，λ_r，t)可以用下面的方程(14)表示，也可以表示为方程(15)： $\mathrm{\β}(\mathrm{\λ},{\mathrm{\λ}}_r,t)=\frac{k(\mathrm{\λ},t)}{k({\mathrm{\λ}}_r,t)}\·\frac{I_0(\mathrm{\λ},t_0)}{I_0({\mathrm{\λ}}_r,t_0)}$ $=\frac{k\left(t\right)}{k\left(t\right)}\·\frac{I_0(\mathrm{\λ},t_0)}{I_0({\mathrm{\λ}}_r,t_0)}$ $=\frac{I_0(\mathrm{\λ},t_0)}{I_0({\mathrm{\λ}}_r,t_0)}---\left(14\right)$

     β(λ，λ_r，t)－β(λ，λ_r).

                                          (15)

从上面的方程(15)可以看到，可使由于时间与光谱中基线偏移的依赖关系所引起的误差部分被减少或消除。

其次，它表示，可通过得出在子区间内的能谱比率减小因光的入射角变化所引起的透射率变化。

根据菲涅耳定律，当入射到被测物体的波长λ的光入射角为θi时的表面透射率P(λ，θ)，以及当入射到被测物体的波长λ_r的光入射角为θ_i时的表面透射率P(λ_r，θ)可以分别由下面的方程(16)、(17)、(18)、(19)表示： $p{(\mathrm{\λ},\mathrm{\θ})}_{\left|\right|}$ $\frac{n_2\·\mathrm{\λ}}{n_1\·\mathrm{\λ}}$ $\frac{{\cos \mathrm{\θ}}_t}{{\cos \mathrm{\θ}}_i}$ ${\left(\frac{{2n}_{1,\mathrm{\λ}}{\cos \mathrm{\θ}}_i}{n_{2,\mathrm{\λ}}{\cos \mathrm{\θ}}_i+n_{1,\mathrm{\λ}}\cos {\mathrm{\θ}}_t}\right)}^2$ $\left(16\right)$ ${p(\mathrm{\λ},\mathrm{\θ})}_{\⊥}=\frac{n_{2,\mathrm{\λ}}}{n_{1,\mathrm{\λ}}}\·\frac{\cos {\mathrm{\θ}}_t}{\cos {\mathrm{\θ}}_i}\·{\left(\frac{{2n}_{1,\mathrm{\λ}}\cos {\mathrm{\θ}}_1}{n_{1,\mathrm{\λ}}\cos {\mathrm{\θ}}_i+n_{2,\mathrm{\λ}}\cos {\mathrm{\θ}}_t}\right)}^2\left(17\right)$ ${p({\mathrm{\λ}}_r,\mathrm{\θ})}_{\left|\right|}=\frac{n_{2,{\mathrm{\λ}}_r}}{n_{1,{\mathrm{\λ}}_r}}\·\frac{\cos {\mathrm{\θ}}_t}{\cos {\mathrm{\θ}}_i}\·{\left(\frac{{2n}_{1,{\mathrm{\λ}}_r}\cos {\mathrm{\θ}}_i}{n_{2,{\mathrm{\λ}}_r}\cos {\mathrm{\θ}}_i+n_{1,{\mathrm{\λ}}_r}\cos {\mathrm{\θ}}_t}\right)}^2\left(18\right)$ ${p({\mathrm{\λ}}_r,\mathrm{\θ})}_{\⊥}=\frac{n_{2,{\mathrm{\λ}}_r}}{n_{1,{\mathrm{\λ}}_r}}\·\frac{\cos {\mathrm{\θ}}_t}{\cos {\mathrm{\θ}}_i}\·{\left(\frac{{2n}_{1,{\mathrm{\λ}}_r}\cos {\mathrm{\θ}}_i}{n_{1,{\mathrm{\λ}}_r}\cos {\mathrm{\θ}}_i+n_{2,{\mathrm{\λ}}_r}\cos {\mathrm{\θ}}_t}\right)}^2\left(19\right)$ 其中：P(λ，θ_i)‖：当水平偏振光以θ_i的入射角入射在被测物体上时，波长为λ的光表面透射率。

P(λ，θ)：当竖直偏振光以θ_i的入射角入射在被测物体上时，波长为λ的光表面透射率。

P(λ_r，θ)‖：当水平偏振光以θ_i的入射角入射在被测物体上时，波长为λ_r的光表面透射率。

P(λ_r，θ)：当竖直偏振光以θ_i的入射角入射在被测物体上时，波长为λ_r的光表面透射率。

n_1,λ：介质1对波长λ的折射率，

n_2,λ：介质2对波长λ的折射率，

n_1、λr：介质1对波长λ_r的折射率，

n_2、λr：介质2对波长λ_r的折射率，

θ_i：对被测物体表面的入射角。和

θ_t：对被测物体表面的折射角。

此外，根据斯涅耳折射定律，受方程(20)的约束：

n₁sinθ_i＝n₂sinθ_t    (20)

现在，用γ(λ，λ_r，θ)来表示P(λ，θ)与P(λ_r，θ)的比，也就是：

     γ(λ，λ_r，θ)＝P(λ，θ)/P(λ_r，θ)    (21)

如果Δλ＝λ_r －λ很小，那么γ(λ，λ_r，θ)可以近似为下式：

  γ(λ，λ_r，θ)≈γ(λ，λ_r)    (22)

从方程(8)、(15)、(22)可得到方程(23)： $I_N=\mathrm{\β}(\mathrm{\λ},{\mathrm{\λ}}_r)\·\mathrm{\γ}(\mathrm{\λ},{\mathrm{\λ}}_r)\·\exp \[-\underset{i}{\mathrm{\Σ}}(a_i\left(\mathrm{\λ}\right)-a_i\left({\mathrm{\λ}}_r\right)c_i)\]\left(23\right)$

从上面方程(23)可看出，可同时减少因时间关系所引起之光谱基线波动的误差和因入射光的入射角变化所引起之透射率起伏的误差。

本发明构思基于上面所描述的思想。本发明的方法是一种光谱测定法，即用光照射被测物体，并测量透射光的能谱或反射光的能谱，从而对被测物体中一种特定的组分进行定量分析。具体地说，该方法把测得的透射光能谱或反射光能谱的波长范围分成多个区间。然后，在每个区间内，求得每个实测能量值与所测的预定波长的能值之比，以稳定光谱的测量值。用这种方法，可实现以稳定的光谱值对特定组分进行定量分析。

此外，按照本发明的一个实施例，提供一种在光谱测定法中稳定光谱的方法，该方法用光照射被测物体，并测量透射光的能谱或反射光的能谱，从而对被测物体的一种特定组分进行定量分析。特别是，该方法把测得的透射光或反射光能谱的波长范围分成多个区间。再在每个区间内求得每个实测能量值与所测的预定波长的能值之比，以稳定光谱的测量值。然后对该比值光谱进行多变量分析。用这种方法，可实现以稳定的光谱值对特定组分进行定量分析。

再有，本发明中被测物体的特定组分是葡萄糖。

本发明的方法把测得的透射光的或反射光的能谱的波长范围分成多个区间。再在每个区间.内求得每个实测能量值与所测的预定波长的能值之比，以稳定光谱的测量值。用这种方法，可同时补偿因时间关系所引起的光谱基线波动误差和因入射光的入射角变化所引起的透射率起伏的误差。

更进一步，本发明的方法是对各子区间内由能谱比所得到的光谱进行多变量分析。

具体地说，通过在各子区间内取比值补偿葡萄糖的透射光能谱或反射光能谱。

按照本发明，可同时减少因时间关系所引起之光谱基线波动的误差和因入射光的入射角变化所引起的透射率变化的误差。由此能够改进用光谱分析法进行定量分析的准确性和可信度。

此外，按照本发明，可通过多变量分析使定量分析的准确性和可信度得到改进。

更进一步，本发明是在葡萄糖作为特定组分的情况下，使用本发明的方法能够以较高的准确性测得葡萄糖的浓度，以使对人体内的血糖的测量可获得很高的准确度和可信性。

下面参考附图，结合对实施例的描述将使本发明的上述以及其他目的和特征通过变得愈为清楚。其中：

图1表示每隔5分钟测得的商品牛奶的11个透射光谱；

图2表示图1中光谱的CV谱；

图3表示取各实测的能量值与测得的波长6000cm^-1的能量值的比所得的能谱；

图4表示图3中所示光谱的CV谱；

图5表示在各子区间内取比值所得到的能谱；

图6表示图5中所示光谱的CV谱；

图7表示通过对图1所示能谱的一次微分所得的能谱；

图8表示图7中所示光谱的CV谱；

图9表示了通过对图1所示能谱的二次微分所得的能谱；

图10表示图9中所示光谱的CV谱；

图11表示由低频滤波富里叶变换所得的能谱；

图12表示图11中所示光谱的CV谱；

图13表示在不经处理情况下定量分析的结果；

图14表示在对单个波长取比值的情况下定量分析的结果；

图15表示在子区间内取比值的情况下定量分析的结果；

图16表示在一次微分情况下定量分析的结果；

图17表示在二次微分情况下定量分析的结果；

图18表示在进行低频滤波富里叶变换情况下定量分析的结果。

下面将参照附图描述本发明的优选实施例。

第一实施例

把商品牛奶注入一石英器皿中，利用Parkin-Elmer公司生产的红外分光镜，在4000cm^-1至8000cm^-1波长范围内每5分钟共11次测量透射光的能谱。然后对透射光能谱的测量结果作下面6种处理：(1)不经处理，(2)取每次实测的能值与所测的单个波长的能值之比，(3)取每次实测的能值与在每个子区间内测得的预定波长的能值之比，(4)一次微分，(5)二次微分，和(6)富里叶变换。然后再对波长在4000cm^-1至8000cm^-1范围内所选定的每一个波长计算透射光能谱的CV值(＝100×标准误差/平均)。下面把通过绘制每个选定波长的CV值而获得的光谱称为CV光谱。

(1)不经处理：

附图1中显示了每隔5分钟测量的商品牛奶的11个透射光能谱。附图2中所示的是它们的CV光谱。

(2)对单个波长的比值计算：

附图3中所示的是通过计算每次实测的能值与所测的波长6000cm^-1的能值之比而获得的附图1所示的透射光的11能谱的比值光谱。图4表示它们的CV谱。

(3)子区间内的比值计算：

把每隔5分钟11次测量商品牛奶的透射光能谱的波长范围4000cm^-1至8000cm^-1分成下面8个波长子区间：

4000cm^-1至4500cm^-1(4250cm^-1)

4500cm^-1至5000cm^-1(4750cm^-1)

5000cm^-1至5500cm^-1(5250cm^-1)

5500cm^-1至6000cm^-1(5750cm^-1)

6000cm^-1至6500cm^-1(6250cm^-1)

6500cm^-1至7000cm^-1(6750cm^-1)

7000cm^-1至7500cm^-1(7250cm^-1)，

7500cm^-1至8000cm^-1(7750cm^-1)

然后，属于每个子区间内的每个透射光能谱除以括弧中所示的波长的能值，以获得如附图5所示的整个区域的比值光谱。图6表示它们的CV谱。

(4)一次微分：

图8表示从图1所示每隔5分钟11次测得的商品牛奶透射光能谱所得的一次微分光谱。图9表示它们的CV光谱。

(5)二次微分：

图9表示从图1所示每隔5分钟11次测得的商品牛奶透射光能谱所得的二次微分光谱。图10表示它们的CV光谱。

(6)富里叶变换：

图11表示从图1所示每隔5分钟11次测得的商品牛奶透射光能谱经富里叶变换所得的光谱。图12中表示它们的CV光谱。

如果我们把图6的CV光谱与图2、4、10、12的CV光谱比较，那么作为一个整体，经过子区间内比值计算所获得的CV光谱值比经不处理、对单个波长的比值计算、一次微分、二次微分、和富里叶变换所获得的CV光谱值小。因此可以看出，经子区间内的比值计算所获得的能谱比较稳定。

第二实施例

在被测者喝了由Shimizu Pharmaceutical公司生产的一种含糖的水溶液(试样名：Trelan75)之后，用近红外光照射该被测者的手指，使用由Parkin-Elmar公司生产的近红外分光镜(型号：System2000)每隔12分钟测量反射光的能谱。在测量反射光能谱的同时，我们用Kyoto Daiichi Kagaku有限公司生产的一种葡萄糖监视器(型号：GM-1320)测量被测者的血糖。我们共进行了44次测量。然后对用近红外分光镜获得的反射光能谱的测量结果作下面6种处理：(1)不经处理，(2)取每次实测的能值与所测的单个波长的能值之比，(3)取每次实测的能值与在每个子区间内测得的预定波长的能值之比，(4)一次微分，(5)二次微分，和(6)富里叶变换。这里，在上面(2)对单个波长的比值计算中，用波长6000cm^-1下的能量值除下面括弧外的每个能量值。在上面(3)子区间内的比值计算中，用括弧内所示波长下的能量值除下面实测的括弧外的能量值。

4250cm^-1(4300cm^-1)，    4350cm^-1(4400cm^-1)，

4450cm^-1(4700cm^-1)，    4500cm^-1(4700cm^-1)，

4550cm^-1(4700cm^-1)，    4600cm^-1(4700cm^-1)，

4650cm^-1(4700cm^-1)，    4750cm^-1(4700cm^-1)，

4800cm^-1(4700cm^-1)，    4850cm^-1(4700cm^-1)，

5500cm^-1(6300cm^-1)，    5600cm^-1(6300cm^-1)，

5700cm^-1(6300cm^-1)，    5800cm^-1(6300cm^-1)，

5900cm^-1(6300cm^-1)，    6000cm^-1(6300cm^-1)，

6100cm^-1(6300cm^-1)，    6200cm^-1(6300cm^-1)，

6400cm^-1(6300cm^-1)，    6500cm^-1(6300cm^-1)，

6700cm^-1(7100cm^-1)，    6800cm^-1(7100cm^-1)，

6900cm^-1(7100cm^-1)，    7000cm^-1(7100cm^-1)，

7200cm^-1(7100cm^-1)，    7300cm^-1(7100cm^-1)，

7400cm^-1(7100cm^-1)，    7500cm^-1(7100cm^-1)。

对于上面同时进行的44次测量所得的数据中，我们用在偶数次所得的22组数据作分析曲线；用在奇数次所得的22组数据计算浓度。具体地说，我们采用上面括弧外所示的28个波长下的吸收率值，并根据PLS多变量分析处理这些吸收率值，进行定量分析。

图13、14、15、16、17、18中分别表示根据(1)不经处理，(2)取每次实测的能值与所测的单个波长的能值之比，(3)取每次实测的能值与在每个子区间内测得的预定波长的能值之比，(4)一次微分，(5)二次微分，和(6)富里叶变换得到的定量分析结果。浓度的相关系数R和计算误差SEP表示在下面的表1中。这里，浓度的计算误差由下面方程(24)定义： $\mathrm{SEP}=\sqrt{\frac{\underset{n=1}{\overset{n}{\mathrm{\Σ}}}{(C_{h,i,\mathrm{ref}}-C_{h,i,\mathrm{obs}})}^2}{n}}$ $\left(24\right)$ 其中，Ch.i.ref：第h次测量时组分i浓度的实际测量值，

  Ch.i.obs：第h次测量时组分i浓度的计算值，

  n：浓度计算所用的数据值。

            表1分离的定量分析的结果

	最小二乘法回归方程	相关系数R	浓度的计算误差SEP(mg/dl)
				不处理	Y＝-13.5231＋1.08041x	0.923082	14.5359
单个波长归一化	Y＝-10.6919＋1.05423x	0.919328	14.6042
				子区间归一化	Y＝-19.4918＋1.11672x	0.947856	12.7113
一次微分	Y＝35.7489＋0.657103x	0.669839	28.3899
				二次微分	Y＝20.5183＋0.998778x	0.882189	26.3624
富里叶变换	Y＝-12.8506＋1.07155x	0.93135	13.6366

把附图15与附图13、14、16、17、和18相比较，我们看到，取子区间内能谱的比所得到的血糖计算值比较接近血糖的实际测量值，由此使浓度的定量分析得到改进。而且，从表1我们看到，血糖的计算值和观察值之间的相关系数用子区间内能谱的比值方法比通过(1)不经处理，(2)单个波长的比值计算，(4)一次微分，(5)二次微分，和(6)富里叶变换后计算的各个相关系数大。此外，子区间内比值计算法使浓度的计算误差SEP最小。

第三实施例

我们研究了光线从空气入射到水上时由于其入射角θi的变化所引起的表面透射率P(λ，θ_i)的变化，并对进行了比值计算的情况和没有进行比值计算的情况作了比较。

当没有进行比值计算时，我们得到的由于入射角θ_i的变化量Δθ_i＝θ_i、1-θ_i、2所引起的表面透射率P(λ，θ_i)的相对变化量ΔP(λ，θ_i)/P(λ，θ_i)由下面的方程(25)表示：

ΔP(λ，θ_i)/P(λ，θ_i)

＝{P(λ，θ_i、1)_‖-P(λ，θ_i、2)_‖)/P(λ，θ_i、1)    (25)

另一方面，当进行了比值计算时，我们得到的由于入射角θ_i的变化量Δθ_i＝θ_i、1-θ_i、2所引起的表面透射率P_N(λ，θ_i)的相对变化量ΔP_N(λ，_i)/P_N(λ，θ_i)，由下面的方程(26)表示：

ΔP_N(λ，θ_i)/P_N(λ，θ_i)

＝{P(λ，θ_i、1)_‖/P(λ_r，θ_i、1)_‖-P(λ，θ_i、2)_‖/P(λ_r，θ_i、2)_‖}/{P(λ，θ_i、1)_‖/P(λ_r，θ_i、1)_‖}    (26)

下面的表2和表3表示计算结果。

                        表2

波长h(μm)	n1	n2	θ；(rad)	P(λ、θ)‖
					1.256	1.00027	1.3210	0	0.980910
1.256	1.00027	1.3210	5	0.981035
					1.256	1.00027	1.3210	10	0.981416
1.256	1.00027	1.3210	15	0.982063

1.256	1.00027	1.3210	20	0.982997
					1.256	1.00027	1.3210	25	0.984242
0.6563	1.00028	1.3311	0	0.979864
					0.6563	1.00028	1.3311	5	0.979996
0.6563	1.00028	1.3311	10	0.980396
					0.6563	1.00028	1.3311	15	0.981076
0.6563	1.00028	1.3311	20	0.982056
					0.6563	1.00028	1.3311	25	0.983363

这里，n1：空气的折射率，引自Rika Nenpyo、NationalAstronomical Observatory(Maruzen公司)，

n₂：水的折射率，根据在Rika Nenpyo中描述的下面的方程(27)计算：

(n_s－1)×10⁸＝6432.8＋2949810/(146-1/λ²)＋25540/(41－1/λ²)(27)

                          表    3

θi，1(rad)	θi，2(rad)	ΔP(λ，θi)/P(λ，θi)(rad-1)	ΔPN(λ，θi)/PN(λ，θi)(rad-1)
				0	5	1.2743×10-4	7.2789×10-6
0	10	5.1584×10-4	2.7070×10-5
				10	15	6.5925×10-4	3.4321×10-5
20	25	1.2665×10-3	6.2460×10-5

如表3所见，表面透射率P_N(λ，θ_i)的相对变化量ΔP_N(λ，θ_i)/P_N(λ，θ_i)比表面透射率P(λ，θ_i)的相对变化量ΔP(λ，θ_i)/P(λ，θ_i)小。因此，通过比值计算减少了由于入射角θ_i的变化而引起的表面透射率P(λ，θ)的变化。也就是说，可以减小由于光对被测物体的入射角的变化而引起的透射率的波动。

虽然，对本发明已参考附图结合实施例作了充分的描述，但应该注意到，它的各种变化和修改对本领域的普通技术人员来说是显而易见的。这些变化和修改都被认为是包括在本发明的权利要求书所限定的范围内。

Claims (4)

1、一种在光谱测定法中稳定光谱的方法，使光照射到被测目标物，并从透射光能谱、散射光能谱和折射光能谱中选定一种光谱进行测量，从而对所说的目标物中的一种特定组分作定量分析，包括下列步骤：

测量穿过所说的目标物的透射光能谱或来自所说目标物的反射光能谱；

把所说的能谱分成多个波长区间，由此得到多个局部的能谱；

对每个波长区间内的所说的多个局部能谱中的每一个能谱用包含在所说的每个波长区间内的一个预定波长的实测能归一化；

用所说的已归一化的多个局部能谱对所说的目标物中的一种特定组分进行定量分析。

2、按照权利要求1所说的在光谱测定法中稳定光谱的方法，其中，通过对所说的多个归一化的局部能谱的多变量分析，进行所说的一种特定组分的定量分析。

3、按照权利要求1或2所说的在光谱测定法中稳定光谱的方法，其中，所说的预定的波长是每个波长区间的中间值、

4、按照权利要求1或2所述的在光谱测定法中稳定光谱的方法，其中所说的特定组分是葡萄糖。

Images