拉曼光谱的数据初步处理Word格式.docx
《拉曼光谱的数据初步处理Word格式.docx》由会员分享，可在线阅读，更多相关《拉曼光谱的数据初步处理Word格式.docx（34页珍藏版）》请在冰豆网上搜索。
但拉曼光谱本身有一定的局限性，比如拉曼散射的强度较弱，对样品进行拉曼散射研究时有强大的荧光及瑞利散射干扰等等。
因此它在相当长一段时间里未真正成为一种有实际应用价值的工具，直到激光器的问世，提供了优质高强度单色光，有力推动了拉曼散射的研究及其应用。
激光使拉曼光谱获得了新生，因为激光的高强度极提高了包含双光子过程的拉曼光谱分辨率和实用性。
此外强激光引起的非线性效应导致了新的拉曼散射现象。
为了进一步提高拉曼散射的强度，人们先后发展了傅立叶变换拉曼光谱、表面增强拉曼光谱、超位拉曼光谱、共振拉曼光谱、时间分辨拉曼光谱等新技术，使光谱仪的效率和灵敏度得到更大的提高。
目前拉曼光谱的应用围遍及化学、物理学、生物学和医学等各个领域，对于定性分析、高度定量分析和测定分子结构都有很大价值。
随着拉曼光谱学研究的深入，拉曼光谱的应用必将愈来愈广泛[4]。
1.2现代拉曼光谱技术与特点
30年代拉曼光谱曾是研究分子结构的主要手段，此时的拉曼光谱仪是以汞弧灯为光源，物质产生的拉曼散射谱线极其微弱，因此应用受到限制，尤其是红外光谱的出现，使得拉曼光谱在分子结构分析中的地位一落千丈。
直至60年代激光光源的问世，以及光电信号转换器件的发展才给拉曼光谱带来新的转机。
世界上各大仪器厂家相继推出了激光拉曼光谱仪，此时拉曼光谱的应用领域不断拓宽。
70年代中期，激光拉曼探针的出现，给微区分析注入活力。
80年代以来，随着科学技术的飞速发展，激光拉曼光谱仪在性能方面日臻完善，如：
美国Spex公司和英国Reinshow公司相继推出了拉曼探针共焦激光拉曼光谱仪，低功率的激光光源的使用使激光器的使用寿命大大延长，共焦显微拉曼的引入可以进行类似生物切片的激光拉曼扫描，从而得出样品在不同深度时的拉曼光谱。
EG＆GDilor公司推出多测点在线工业用拉曼系统，采用的光纤可达200m，从而使拉曼光谱的应用围更加广阔。
90年代初，由于社会生产活动的需要，人们又探索出多项技术并应用于拉曼光谱仪中，使小型便携式拉曼光谱仪出现并不断发展起来成为可能。
这些技术包括：
引进光纤对远距离或危险的样品进行测量；
用声光调制器代替光栅作为分光元件测量拉曼光谱；
利用全息带阻滤光片滤除瑞利散射的干扰；
研制开发出便携激光器等[5][6]。
1.3研究拉曼光谱仪的意义
由于拉曼光谱具有制样简单、水的干扰少、拉曼光谱分辨率较高等特点，故其可以广泛应用于有机物、无机物以及生物样品的应用分析中。
拉曼光谱技术己广泛应用于医药、文物、宝石鉴定和法庭科学等领域。
对文物样品的无损分析研究。
使文物的鉴定、年代的测定及文物的恢复和保存的方法更安全可靠；
对爆炸物、毒品、墨迹等的痕迹无损检测为法庭提供科学证据的有力手段：
对宝石的光谱分析研究对认识各地宝石中的包含物差异性。
并使宝石的鉴别与评价有了科学依据。
近年来该技术在细胞和组织的癌变方面的检测也取得了很大的进展，随着分析方法完善和研究病例的增多以及对于病变组织差异性的规律性认识深化。
拉曼光谱发展成诊断肿瘤方法的可行性将得到确认．总之，随着激光技术的发展和检测装置的改进。
拉曼光谱技术在当代工业生产和科学研究中必将得到越来越广泛的应用[7]。
1.4本文的主要容
本文主要论述了拉曼光谱仪开发设计、安装调试中所应用的基本理论、设计原理与关键技术，介绍了激光拉曼光谱仪的发展动态、研究方向和国外总体概况。
阐述了拉曼散射原理及其量子解释。
以具体说明了分析拉曼光谱数据的各种可行的方法，包括平滑，滤波等方法的使用。
根据光谱仪器设计原理详细论述了分光计光学系统的结构设计、激光拉曼光谱仪的总体设计。
并且对各个部件的选择作用及原理分析，做了详细的描述。
最后，测量了几种样品的拉曼光谱，并对光谱利用文中阐述的光谱分析方法进行分析对比，并且进行了合理的分析。
拉曼光谱仪的实验测量和光谱数据处理研究主要从两个方面来分析：
二、拉曼光谱数据处理分析，用合理的方法分析拉曼光谱可以有效便捷的得到较为理想的实验结果。
通过对四氯化碳、乙醇、正丁醇的测量光谱以及光谱数据分析，得到较为理想实验效果，证明本文所论述方法的可行性和正确性。
第2章基本理论
当一束频率为的单色光照射到样品上后，分子可以使入射光发生散射。
大部分光只是改变方向发生散射，而光的频率仍与激发光的频率相同，这种散射称为瑞利散射；
约占总散射光强度的，不仅改变了光的传播方向，而且散射光的频率也改变了，不同于激发光的频率，称为拉曼散射。
拉曼散射中频率减少的称为斯托克斯散射，频率增加的散射称为反斯托克斯散射，斯托克斯散射通常要比反斯托克斯散射强得多，拉曼光谱仪通常测定的大多是斯托克斯散射，也统称为拉曼散射。
散射光与入射光之间的频率差v称为拉曼位移，拉曼位移与入射光频率无关，它只与散射分子本身的结构有关。
拉曼散射是由于分子极化率的改变而产生的。
拉曼位移取决于分子振动能级的变化，不同化学键或基团有特征的分子振动，ΔE反映了指定能级的变化，因此与之对应的拉曼位移也是特征的。
这是拉曼光谱可以作为分子结构定性分析的依据。
2.1拉曼散射经典解释[8]
光照射到物质上发生弹性散射和非弹性散射。
弹性散射的散射光是与激光光波波长相同的成分，非弹性散射的散射光有比激发光波长长的和短的成分，统称为拉曼效应。
角频率为的光入射到一个分子上，可以感应产生电偶极矩。
一级近似下，所产生的感应电偶极矩P与入射光波电场E的关系可表达为下式：
P=A·
E式中，A是一个二阶量，通常称A为极化率量。
如果角频率为的入射光波只感生振荡角频率为叫的感应电偶极矩，该感生电偶极矩会辐射出与入射光角频率相同的散射光，也就是瑞利散射。
但若考虑到分子部本身有振动和转动，各有其特征频率，导致激发光每个周期所遇的分子振动和转动相位不同，相应的极化率也不同，分子的感生偶极发射受自身振动和转动频率调制，会辐射出异于入射光频率的散射光，其中波长比瑞利光长的拉曼光叫斯托克斯线，比瑞利光短的叫反斯托克斯线。
考虑分子中的原子由于热运动而在平衡位置附近振动，那么，P=A·
E可以写作：
同理：
（2-1）
可知，其中是极化率A的E分量。
（2-1）式中的x，y，z是固定在分子上的坐标系的三个坐标轴，由于假设没有转动，这个坐标系也是固定在空间上的。
是和P与E的方向无关的常数，也就是分子极化率量A的分量。
可以知道：
。
一般情况下，当各个原子核从其平衡位置有一位移时，极化率的六个分量中的每个分量都会发生改变。
对于小位移的情形，可以把展开并保留到一级项
（2-2）
式中（）0表示分子处于平衡状态时物理量的值，，是引入的振动简正坐标，求和遍及全部简正坐标。
由于考虑的是分子部振动小位移的情况，振动可近似为简谐，于是得：
（2-3）
其中表示振动的幅度，，表示振动的频率和初相位。
又，
（2-4）
将（2-2）（2-3）（2-4）代入（2-1），得：
（2-5）
同理，对于Py，PZ也能得到类似的式子。
综上所述，感生偶极矩的振动情况如下：
（1）以入射辐射的频率振动，结果也就是瑞利（Rayleigh）散射；
（2）以频率振动，结果也就是拉曼散射，频率为的散射光是斯托克斯线，频率为的散射光是反斯托克斯线。
从（2-5）式还可以断言，不同分子间瑞利散射光彼此之间是相干的。
而因为公式中含的项只是纯粹的叠加而没有交叉项，所以对于多分子体系，其拉曼散射总强度是各个分子拉曼散射强度的代数和，拉曼散射光不相干。
2.2拉曼散射的量子解释
2.2.1散射过程的量子跃迁
图2-1瑞利、拉曼散射过程中的量子跃迁
拉曼散射的完善解释需用量子力学理论，不仅可解释散射光的频率差，还可解决强度和偏振等问题。
图2-1给出散射过程量子跃迁的三能级图，其中、分别表示激光入射光子的频率和波矢，、分别表示散射光子的频率和波矢，ωq和q分别表示散射过程中伴随产生或湮没的元激发的频率和波矢。
在入射光（量）子被吸收后，使电子和晶格振动从初态（,nq）跃迁到一个虚中间态；
随即辐射出散射光子（,）由中间虚态回到终态，与此同时，产生（或淹没）了一个频率为ωq而波矢为q元激发。
2.2.2量子力学结果
核与电子组成的系统遵从的薛定谔方程为
（2-6）
式中r代表各粒子的所有坐标，它的通解为
（2-7）
对不含时薛定谔方程的本征值和本征值函数分别是和,对k态，即k=（e，n），e和n分别是电子量子数及核量子数集合。
对于，r≠k、=0和=1，其通解为
（2-8）
因系统受到的微扰来自于光波电磁场，而光波波长远大于原子间距。
显然，这些理论对可见光、紫外、红外光都是正确的。
对X射线的结论则不适用。
为了简单起见，先不考虑共振现象，则光波电磁场可以写成如下的形式：
（2-9）
式中A是复振幅，则光波场与系统的微扰互作用能为
而
，即系统中的电子偶极矩。
此时微扰系统的薛定谔方程为
（2-10）
若k态中未受微扰的系统由
所描述，则（2-10）方程的微扰解为
（2-11）
将（2-11）式代入（2-10）方程中，略去
的二次项，再由（2-6）方程得以下方程：
（2-12）
对（2-12）做求解处理，取
（2-13）
式中
，
（2-14）
，
（2-15）
受到微扰系统的矩阵元为
（2-16）
由（2-11）到（2-16）式可得
（2-17）
式中的和分别为
（2-18）
（2-19）
又因为=，对k=m的条件可得
（2-20）
是实的，它是k态中偶极子动量的期待值，且与入射辐射有相同的时间关系。
因此，偶极子辐射的强度仍有以下的经典表示式：
（2-21）
由（2-20）式可得
（2-22）（2-21）式给出了（偶）极矩为的偶极子的瑞利散射光强。
需特别注意的是：
与相反，是复的。
要找出与（2-17）式中个别真实偶极子经典辐射相关联的情况，必须用到克莱因（Klein）的结果：
若,即,初态能量小于末态能量，k、m别为初、末态，则
分量的辐射发射为零，就等价于真实偶极子的经典辐射，则有
若，即,初态能量小于末态能量，k,m分别为初、末态。
为了能应用（2-17）式，必先考虑构成真实偶极子的情况，即
（2-23）
与（2-11）式相同，散射光光强由下式决定：
（2-24）
在对时间取平均时消去交叉项后有
（2-25）
由辐射发射的原理知：
仅对,和的条件下才能产生辐射，下面对（2-25）中各项的意义作以讨论。
表示式（2-25）中的的第一项初态能量（），大于末态能量（）。
它描述了与外来激光频率无关的伴随k→m跃迁的自发辐射，见图2-1（a）。
（2-25）式中的第二项是正常拉曼散射，即。
末态（m）的能量比始态（k）的能量大，也可以比它小：
（1）>
末态能量大于始态能量；
散射辐射能量小于激光（单）光子能量,
即这对应于斯托克斯过程。
（2）<
态能量小于始态能量；
散射辐射能量大于激光（单）光子能量，即这对应于反斯托克斯过程。
需要指出的是：
光谱仪所接收到的信号是，见图（2-2）（b）。
（2-25）式中的第三项表示伴有两个量子感应发射，即k→m的跃迁。
这类发射只有在受激粒子数剧增时才能被观测到。
图（2-2）（c）给出了这一过程的量子跃迁，特别注意激光器的能量与初末能态相应能级的关系。
图2-2斯托克斯跃迁（<
），，和是跃迁过程中的中间能级
与正常拉曼效应相关的第二项系数由（2-18）式确定，其中跃迁矩的求和是从初态k到所有未受微扰系统的r态及跃迁矩从r态回到末态m的求和。
并不是说真实的散射过程中存在如上的跃迁过程，完全是因为考虑到数学上处理微扰问题的需要，即受到微扰系统的波函数完全用未受到微扰波函数来表示，也就是在的表示中包括了矩跃迁和的积，而不是单独的某一个跃迁矩。
跃迁矩既可以是正的，也可以是负的；
既属于不同的态r，也可以附加在另一个态上。
当它们所处的态被湮没为另一个使变为零的态时，就不能产生拉曼散射，常称之为拉曼散射的禁戒条件。
在对（2-18）式的求和中，不仅包括了初态k之上的中间态r，也包括了低于k态的任何r态。
因为中间态是在吸收了入射光子后产生的高于初态的激发态，所以（2-18）式求和过程中包含的低于初态k的概念显然是不合适的。
还应强调的是：
确定自发初末态间跃迁概率的跃迁矩并没有全部写入的表达式。
因此，自发辐射（或吸收）与拉曼散射强度间无直接的联系，它们的选择规律也全然不同。
由上述的讨论可知：
整个空间4兀立体角拉曼散射强度为
（2-26）
对于k=m，就是瑞利散射强度。
的分量可以写成如下的形式：
（2-27）
式中代表Portor表示中入、散射的偏振方向，式（2-27）中
（2-28）
（2-28）式就是散射量，通常它是复数，而且是非对称的量。
对于k=m，则有
（2-29）
因
（2-30）
因此，若
是实的，则它是对称的，该结论不仅对=0的静电场微扰是正确的，而且对哈密顿量为实的系统也是适用的。
若不再考虑（2-20）式中的永久偶极矩，而利用在（2-9）式中的系数关系
，再从（2-20）式可得矩阵元表示式，为
（2-31）式中是实的，就是k态的电极化率，将（2-27）式代入（2—26）式中有
（2-32）
再用，是入射光偏振方向的单位失，而入射光光强为，则有
（2-33）
被定义为km跃迁的拉曼跃迁截面。
的量纲是，和的单位分别是erg*/s和erg。
若入射光沿方向偏振，在沿方向用分析检测器观察散射光，则单位立体角
中散射强度为
（2-34）
用
直接从（2-31）到（2-33）式计算散射强度只有在简单系统（如谐振子、自由电子、和某些简单原子）中是可行的。
对于分子、晶体及复杂能级系统中，（2-28）式中出现的受激中间态r都无法进行直接的计算。
Plaezek’s近似为直接计算一般结果提供了新的途径[9][10]。
2.2.3Placzek近似
首先考虑电子态不发生变化的散射过程，即初、末态是相同的基态。
该情况中，仅是振动态发生了变化，而且满足能量守恒条件，即
因为光散射是由于系统中的电子引起的，光能转移到各个核上，反之亦然，通过核和电子运动间的耦合可以产生拉曼散射效应。
假设电子的基态是非简并的，而且原子核被固定在仅产生瑞利（辐）散射的位置，则散射光的强度由电极化率量
确定，（2-31）式下标中的k=0表示电子的基态。
电子的极化率是实的，而且具有对称性：
在（2-29）式中，对具有
固定核系的本征频率
和本征函数影
不仅取决于核的位置，而且也取决于电子极化率，该分量是核组态R的函数，即）。
在以下的假设中，认为振动着的核系统中散射强度是一样的。
Plaezek假设有以下三点容：
1．电予的基态必须是非简并的；
2．绝热近似必须是有效的；
3．激发光源的频率必须小于任何一电子的跃迁频率，但远大于振动的频率，即
若具有振动核系统从
态跃迁到
态，其中O表示电子的基态而
、
是振动态。
根据近似条件，
可认为矩阵元是由电子极化率分量矩阵元
所确定，也就是由下式表示：
（2-35）
是电子基态O和振动态
的振动波函数。
将（2-35）式代入（2-34）式，而且略去下指标O，就可以得到伴随
振动跃迁光散射强度的表示式为
（2-36）
式中的
分别是两个振动态确定的布局数即
和
，而
由下式确定：
（2-37）
式中的q式简正坐标的脚码，表示可能出现的声子振动频率的个数。
在Placzck近似条件下由n态跃迁到态的拉曼散射的光强还需做进一步的讨论。
若入射光沿着方向偏振，沿方向观测散射光，则按（2-34）式可得到单位立体角散射光强，为
（2-38）
因为极化率取决于核的配置R，所以可以将极化率按简正坐标展成以下的级数：
（2-39）
这是所有量子数假设不变伴有跃迁的极化率矩阵元。
若略去式中的正比项，则有
（2-40）
若、
，则有，这是瑞利散射。
由（2-38）式知可以得到它的散射光强，为
（2-41）
因为在式（2-40）中略去了（2-39）式中的二次项，所以瑞利散射光强与温度T无关。
以下讨论拉曼散射光强：
对第一级拉曼散射有，而斯托克斯散射有。
由此可得
（2-42）
对于反斯托克斯散射有，由此可得
（2-43）
由（2-38）、（2-42）和（2-43）式可得拉曼散射强度表示式。
由（2-37）式可知：
对斯托克斯和反斯托克斯散射分布有和相应于这两个散射的光强分别为
（2-44）
（2-45）
是某一元激发q的平均布居数，若元激发是（热）声子，则它满足玻色一爱因斯坦分布即
（2-46）
图2-3给出了声子平均布居数随温度增加而变化的情况：
由于温度升高，处于较高能级的布居数也随之增加，分布状况发生了变化，使可能参与跃迁的声子“种类”有所增加，也就是拉曼谱峰中峰的数目增加，由原来的两种，增加到四种。
在小于100K的低温区，仅实现了（）的声子，随着温度增加又激发了（），（）和（），三种声子且。
图2-3玻色．爱因斯坦分布得出的平均布居数随温度的变化
图2-3是低、高温能级布居数变化情况的示意。
拉曼散射光强与受激后跃迁的元激发数成正比，因此可以得出以下结论：
这就是拉曼谱中振动“热带”产生的原因。
由（2-46）式可得出参与热激发的声子数，为
（2-47）
较高温度下，“热”声子数几乎与温度成正比（图2-3）。
由（2-44）和（2-45）式有
（2-48）
解决了经典电磁理论在解释斯托克斯与反斯托克斯散射光强比时的困难[11][12]。
2.3拉曼光谱数据分析方法
光谱分析技术的数据处理主要涉及两个方面的容：
一是光谱预处理方法的研究，目的是针对特定的样品体系，通过对光谱的适当处理，减弱和消除各种非目标因素对光谱的影响，净化谱图信息，为校正模型的建立和未知样品组成或性质的预测奠定基础；
二光谱定性和定量方法的研究，目的在于建立稳定、可靠的定性或定量分析模型，并最终确定未知样品和对其定量[13][14]。
MATLAB是Mathworks公司开发的一种主要用于数值计算及可视化图形处理的高级计算语言。
它将数值分析、矩阵计算、图形、图象处理和仿真等诸多功能集成在一个极易使用的交互式环境之中，为科学研究、工程设计及数据处理和数值计算提供了一种高效率的编程工具”。
在这种编程环境下，任何复杂的计算问题及其解的描述均符合人们的科学思维方式和数学表达习惯，而不像Fortran、Basic、C等高级程序设计语言那样难以学习和掌握。
MATLAB允许用户根据数值计算的复杂程度，对问题进行分段甚至逐句编程处理，显然，这是与C、Fortran等传统高级语言完全不同的。
此外，用MATLAB求解问题一般不需要用户考虑采用何种算法以及怎样具体实现等低层问题，更不必深入了解相应算法的具体细节，因而对用户算法语言方面的要求比较低。
仪器采集的原始光谱中除包含与样品组成有关的信息外，同时也包含来自各方面因素所产生的噪音信号。
这些噪音信号会对谱图信息产生干扰，有些情况下还非常严重，从而影响校正模型的建立和对未知样品组成或性质的预测。
因此，光谱数据预处理主要解决光谱噪音的滤除、数据的筛选、光谱围的优化及消除其他因素对数据信息的影响，为下步校正模型的建立和未知样品的准确预测打下基础。
常用的数据预处理方法有光谱数据的平滑、基线校正、求导、归一化处理等。
2.3.1数据平滑处理
数据平滑处理：
信号平滑是消除噪声最常用的一种方法，其基本假设是光谱含有的噪声为零均随机白噪声，若多次测量取平均值可降低噪声提高信噪比。
平滑处理常用方法有邻近点比较法、移动平均法、指数平均法等[13][14][15]。
（1）邻近点比较法
对于许多干扰性的脉冲信号，将每一个数据点和它旁边邻近的数据点的值进行比较可以测得其存在。
如果与邻近点的数值相差太大，超过给定的阈值，便可认为该数据是一个脉冲干扰，并通过邻近数据点的平均值来取代这一数据点值，就可以把这一干扰脉冲去掉，这样不影响信号的其它部分。
在这一数据点处理过程中，需注意选择调节参数，也就是考虑邻近数据点值，以及判断一个数据点和邻近数据点之间不同的阈值。
这个阈值一般定义为噪音测量偏差的倍数，以免把必要的有用信号去掉。
这一方法有时也称为邻近点平滑法，也叫做单点平滑法。
（2）移动平均法
由于平滑是通过对信号进行平均而减小噪音，因而多点平滑效果更好。
移动平均法是多点平滑中最简单的一种。
先选择在数据序列中相邻的奇数个数据点，这奇数个数据点即构成一个窗口。
计算在窗口奇数个数据点的平均值，然后
升级会员 