拉曼光谱的数据初步处理

1、摘要本文主要目的是熟悉拉曼光谱仪原理，并掌握拉曼光谱仪的实验测量技术以及拉曼光谱的数据初步处理。文章首先论述了拉曼光谱仪开发设计、安装调试中所应用的基本理论、设计原理与关键技术，介绍了激光拉曼光谱仪的发展动态、研究方向和国内外总体概况。其次阐述了拉曼散射的经典理论及其量子解释。并说明了分析拉曼光谱数据的各种可行的方法，包括平滑，滤波等。再次根据光谱仪器设计原理详细论述了分光光学系统的结构设计和激光拉曼光谱仪的总体设计，并且对各个部件的选择作用及原理做了详细的描述。最后，测量了几种样品的拉曼光谱，并利用文中阐述的光谱处理方法进行初步处理，并且进行了合理的分析对比。总之，本文主要从两个方面来分析拉

2、曼光谱仪的实验测量和光谱数据处理研究：一、拉曼光谱仪的结构，详细了解拉曼光谱仪的工作原理。二、拉曼光谱数据处理分析，用合理的方法处理拉曼光谱可以有效便捷的得到较为理想的实验结果。通过对四氯化碳、乙醇、正丁醇的光谱测量以及光谱数据分析，得到了较为理想实验效果，证明本文所论述方法的可行性和正确性。关键词：拉曼光谱仪光栅光谱分析AbstractPurposeofthispaperisfamiliarwithRamanSpectrometer,andmasteryofexperimentalmeasurementsofRamanspectroscopyandRamanspectroscopytechn

3、iquespreliminarydataprocessing.TliearticlefirstdiscussestheRamanspectrometerdevelopment,designjnstallationandcommissioningintheapplicationofthebasictheory,designprinciplesandkeytechnologiesJaserRamanspectrometerdevelopmentsjesearchdirectionandoverallprofileathomeandabroad.Tliesecondsectiondescribest

4、heclassicaltheoryofRamanscatteringandquantumexplanation.AndshowstheRamanspectraofthevariouspossibleways,includingsmoothingandfiltering.AgainaccordingtospectrometerdesignprinciplesdiscussedindetailthespectroscopicopticalsystemdesignandlaserRamanspectrometeroveralldesign,andthechoicefortheroleofthevar

5、iouscomponentsandtheprincipleofadetaileddescription.Finally,themeasuredRamanspectraofseveralsamples,andusepaperdescribesmethodsforspectralprocessing!nitialtreatment,andforareasonableanalysisandcomparison.Insummary,thispapermainlyfromtwoaspectstoanalyzeexperimentalmeasurementsofRamanspectroscopyandsp

6、ectraldataprocessingresearch:First,thestructureofRamanspectroscopy,Ramanspectroscopydetailedunderstandingoftheworkingprinciple.Second,Ramanspectroscopydataprocessingandanalysis,areasonableapproachtoeffectiveandconvenientRamanspectroscopycanbemoreidealresults.Throughcarbontetrachloride,ethanol,n-buta

7、nolandspectraldataanalysisspectralmeasurementsobtainedmoresatisfactoryexpeiinientalresultsdiscussedinthisarticledemonstratethefeasibilityandcorrectness.Keywords:Ramanspectrometergratingspectralanalys目录第1章引言111拉曼光谱分析技术112现代拉曼光谱技术与特点213研究拉曼光谱仪的意义214本文的主要内容3第2章基本理论421拉曼散射经典解释网422拉曼散射的量子解释61.11 散射过程的量子跃

8、迁6222量子力学结果71.12 3Placzek近似121.0.0 3拉曼光谱数据分析方法153.3 .3.1数据平滑处理154.4 32基线校正175.5 33数据求导处理176.6 .34数据增强算法187.7 .35傅里叶变换188.8 36小波变换199.9 .37数字滤波19第3章常规拉曼检测系统212.0.0 1光源213.2滤光片2333拉曼光谱仪及计算机软件241.1.1 .3.1光栅252 32光电倍增管27第4章拉曼光谱测量及数据处理和结论293.0.0 1物质的拉曼光谱测量294.0.0 2拉曼光谱数据处理与分析323 平滑处理323 低通滤波处理3543结论37第5章

9、论文总结与展望38致谢参考文献：40第1章引言拉曼光谱分析技术1928年印度实验物理学家拉曼发现了光的一种类似于康普顿效应的光散射效应，称为拉曼效应。简单地说就是光通过介质时由于入射光与分子运动之间相互作用而引起的光频率改变。拉蛀因此获得1930年的诺贝尔物理学奖，成为第一个获得这一奖项并且没有接受过西方教育的亚洲人拉曼散射遵守如下规律：散射光中在每条原始入射谱线（频率为9）两侧对称地伴有频率为。o土0（k=l,2,3,）的谱线，长波一侧的谱线称红伴线或斯托克斯线，短波一侧的谱线称紫伴线或反斯托克斯线：频率差必与入射光频率血无关，由散射物质的性质决定，每种散射物质都有H己特定的频率差，其中有些

10、与介质的红外吸收频率相一致。拉曼光谱即拉曼散射的光谱。靠近瑞利散射线两侧的谱线称为小拉曼光谱:远离瑞利散射线的两侧出现的谱线称为大拉曼光谱。拉曼散射的强度比瑞利散射要弱得多。瑞利散射线的强度只有入射光强度的千分之一，拉曼光谱强度大约只有瑞利线的千分之一。与入射光频率uo相同的成分称为瑞利散射，频率对称分布在为两侧的谱线或谱带为称为拉曼散射u拉曼光谱的理论解释是：入射光子与分子发生非弹性散射，分子吸收频率为小的光子，发射-%的光子，同时分子从低能态跃迁到高能态（斯托克斯线）：分子吸收频率为4的光子，发射i）o+阪的光子，同时分子从高能态跃迁到低能态（反斯托克斯线）与分子红外光谱不同，极性分子和非

11、极性分子都能产生拉曼光谱。拉曼光谱为研究晶体或分子的结构提供了重要手段，在光谱学中形成了拉曼光谱学的一个分支。用拉曼散射的方法可迅速定出分子振动的固有频率，并可决定分子的对称性、分子内部的作用力等。但拉曼光谱本身有一定的局限性，比如拉曼散射的强度较弱，对样品进行拉曼放射研究时有强大的荧光及瑞利散射干扰等等。因此它在相当长一段时间里未真正成为一种有实际应用价值的工具，直到激光器的问世，提供了优质高强度单色光，有力推动了拉曼散射的研究及其应用。激光使拉曼光谱获得了新生，因为激光的高强度极大地提高了包含双光子过程的拉曼光谱分辨率和实用性。此外强激光引起的非线性效应导致了新的拉曼散射现象。为了进一步提

12、高拉曼散射的强度，人们先后发展了傅立叶变换拉曼光谱、表面增强拉曼光谱、超位拉曼光谱、共振拉曼光谱、时间分辨拉曼光谱等新技术，使光谱仪的效率和灵敏度得到更大的提高。目前拉曼光谱的应用范围遍及化学、物理学、生物学和医学等各个领域，对于定性分析、高度定量分析和测定分子结构都有很大价值。随着拉曼光谱学研究的深入，拉曼光谱的应用必将愈来愈广泛。现代拉曼光谱技术与特点30年代拉曼光谱曾是研究分子结构的主要手段，此时的拉曼光谱仪是以汞弧灯为光源，物质产生的拉曼散射谱线极其微弱，因此应用受到限制，尤其是红外光谱的出现，使得拉曼光谱在分子结构分析中的地位一落千丈。直至60年代激光光源的问世，以及光电信号转换器件

13、的发展才给拉曼光谱带来新的转机。世界上各大仪器厂家相继推出了激光拉曼光谱仪，此时拉曼光谱的应用领域不断拓宽。70年代中期，激光拉曼探针的出现，给微区分析注入活力。80年代以来，随着科学技术的飞速发展，激光拉曼光谱仪在性能方面H臻完善,如：美国Spex公司和英国Reinshow公司相继推出了拉曼探针共焦激光拉曼光谱仪，低功率的激光光源的使用使激光器的使用寿命大大延长，共焦显微拉曼的引入可以进行类似生物切片的激光拉曼扫描，从而得出样品在不同深度时的拉曼光谱。EG&GDil”公司推出多测点在线工业用拉曼系统，采用的光纤可达200m,从而使拉曼光谱的应用范围更加广阔。90年代初，由于社会生产活

14、动的需要，人们乂探索出多项技术并应用于拉曼光谱仪中，使小型便携式拉亚光谱仪出现并不断发展起来成为可能。这些技术包括：引进光纤对远距离或危险的样品进行测量：用声光调制器代替光栅作为分光元件测量拉曼光谱；利用全息带阻滤光片滤除瑞利散射的干扰；研制开发出便携激光器等L3研究拉曼光谱仪的意义由于拉曼光谱具有制样简单、水的干扰少、拉曼光谱分辨率较高等特点，故其可以广泛应用于有机物、无机物以及生物样品的应用分析中。拉曼光谱技术己广泛应用于医药、文物、宝石鉴定和法庭科学等领域。对文物样品的无损分析研究。使文物的鉴定、年代的测定及文物的恢复和保存的方法更安全可靠:对爆炸物、毒品、墨迹等的痕迹无损检测为法庭提供

15、科学证据的有力手段：对宝石的光谱分析研究对认识各地宝石中的包含物差异性。并使宝石的鉴别与评价有了科学依据。近年来该技术在细胞和组织的癌变方面的检测也取得了很大的进展，随着分析方法完善和研究病例的增多以及对于病变组织差异性的规律性认识深化。拉曼光谱发展成诊断肿瘤方法的可行性将得到确认.总之，随着激光技术的发展和检测装置的改进。拉曼光谱技术在当代工业生产和科学研究中必将得到越来越广泛的应用。1.4本文的主要内容本文主要论述了拉曼光谱仪开发设计、安装调试中所应用的基本理论、设计原理与关键技术，介绍了激光拉曼光谱仪的发展动态、研究方向和国内外总体概况。阐述了拉曼散射原理及其量子解释。以具体说明了分析拉

16、曼光谱数据的各种可行的方法，包括平滑，滤波等方法的使用。根据光谱仪器设计原理详细论述了分光计光学系统的结构设计、激光拉曼光谱仪的总体设计。并且对各个部件的选择作用及原理分析，做了详细的描述。最后，测量了几种样品的拉曼光谱，并对光谱利用文中阐述的光谱分析方法进行分析对比，并且进行了合理的分析。拉曼光谱仪的实验测量和光谱数据处理研究主要从两个方面来分析：一、拉曼光谱仪的结构，详细了解拉曼光谱仪的工作原理。二、拉曼光谱数据处理分析，用合理的方法分析拉曼光谱可以有效便捷的得到较为理想的实验结果。通过对四氯化碳、乙醇、正丁醇的测量光谱以及光谱数据分析，得到较为理想实验效果，证明本文所论述方法的可行性和正

17、确性。第2章基本理论当一束频率为U0的单色光照射到样品上后，分子可以使入射光发生散射。大部分光只是改变方向发生散射，而光的频率仍与激发光的频率相同，这种散射称为瑞利散射；约占总散射光强度的10610一叫不仅改变了光的传播方向，而且散射光的频率也改变了，不同于激发光的频率，称为拉曼散射。拉曼散射中频率减少的称为斯托克斯散射，频率增加的散射称为反斯托克斯散射，斯托克斯散射通常要比反斯托克斯散射强得多，拉曼光谱仪通常测定的大多是斯托克斯散射，也统称为拉曼散射。散射光与入射光之间的频率差V称为拉曼位移，拉曼位移与入射光频率无关，它只与散射分子本身的结构有关。拉曼散射是由于分子极化率的改变而产生的。拉曼

18、位移取决于分子振动能级的变化，不同化学键或基团有特征的分子振动，AE反映了指定能级的变化，因此与之对应的拉曼位移也是特征的。这是拉曼光谱可以作为分子结构定性分析的依据。拉曼散射经典解释光照射到物质上发生弹性散射和非弹性散射。弹性散射的散射光是与激光光波波长相同的成分，非弹性散射的散射光有比激发光波长长的和短的成分,统称为拉曼效应。角频率为uo的光入射到一个分子上，可以感应产生电偶极矩。一级近似下，所产生的感应电偶极矩P与入射光波电场E的关系可表达为下式：P二A-E式中，A是一个二阶张量，通常称A为极化率张量。如果角频率为四的入射光波只感生振荡角频率为叫4的感应电偶极矩，该感生电偶极矩会辐射出与

19、入射光角频率相同的散射光，也就是瑞利散射。但若考虑到分子内部本身有振动和转动，各有其特征频率，导致激发光每个周期所遇的分子振动和转动相位不同，相应的极化率也不同，分子的感生偶极发射受自身振动和转动频率调制，会辐射出异于入射光频率的散射光，其中波长比瑞利光长的拉曼光叫斯托克斯线，比瑞利光短的叫反斯托克斯线。考虑分子中的原子由于热运动而在平衡位置附近振动，那么，P二A-E可以写作：PxfEx+aoEy+aJz同理：Pyna.Ex+a/y+a.EN（2-1）PLa.Ex+aqEy+aNEz可知Pe=Ae|E|,其中Ae是极化率A的E分量。（21）式中的x,y,z是固定在分子上的坐标系的三个坐标轴，由

20、于假设没有转动，这个坐标系也是固定在空间上的。询是和P与E的方向无关的常数，也就是分子极化率张量A的分量。可以知道：知=知。一般情况下，当各个原子核从其平衡位置有一对位移时，极化率的六个分量中的每个分量都会发生改变。对于小位移的情形，可以把展开并保留到一级项式中（）0表示分子处于平衡状态时物理量的值，Q，0是引入的振动简正坐标，求和遍及全部简正坐标。由于考虑的是分子内部振动小位移的情况，振动可近似为简谐，于是得:(2-3)1=或cos(6t+%)=或cos(2ft+外)其中4表示振动的幅度，W，9k表示振动的频率和初相位。又,Ex=E°2t,Ey=Ey2t,Ez=E°2n)

21、t(2-4)将(22)(23)(24)代入(21),得:(2-5)E=®)(X+(%)国+(aJE口C0S27U+Z向。理+(就中+(短)。E*cos2(qj+14)t+cos2(qj-14)12同理，对于P,Pz也能得到类似的式子。综上所述，感生偶极矩的振动情况如下：（1）以入射辐射的频率刖振动，结果也就是瑞利（Rayleigh）散射；（2）以频率D。士咏振动，结果也就是拉曼散射，频率为1）0-%的散射光是斯托克斯线，频率为小+%的散射光是反斯托克斯线。从Q5）式还可以断言，不同分子间瑞利散射光彼此之间是相干的。而因为公式中含为土”的项只是纯粹的叠加而没有交叉项，所以对于多分子体系

22、，其拉曼散射总强度是各个分子拉曼散射强度的代数和，拉曼散射光不相干。拉曼散射的量子解释散射过程的量子跃迁拉曼散射（弹性）瑞利反射斯托克斯散射斯托克斯放射电=0kL=ksCOL=%+，°=0一（0Lsq施二与+4kL=ksq图21瑞利、拉曼散射过程中的量子跃迁拉曼散射的完善解释需用量子力学理论，不仅可解释散射光的频率差，还可解决强度和偏振等问题。图21给出散射过程量子跃迁的三能级图，其中3L、此分别表示激光入射光子的频率和波矢，o）s、ks分别表示散射光子的频率和波矢，sq和q分别表示散射过程中伴随产生或湮没的元激发的频率和波矢。在入射光（量）子被吸收后，使电子和晶格振动从初态（He，

23、nq）跃迁到一个虚中间态；随即辐射出散射光子（3s，&）由中间虚态回到终态，与此同时，产生（或淹没）了一个频率为3q而波矢为q元激发。量子力学结果核与电子组成的系统遵从的薛定丹方程为二访与P°)(r,t)(2-6)式中r代表各粒子的所有坐标，它的通解为甲。(二名")©"(2-7)X对不含时薛定娉方程How=耳。()的本征值和本征值函数分别是Er和Ho,对k态，B|Jk=(e,n),e和n分别是电子量子数及核量子数集合。对于，/k、0r=0和二1,其通解为K°)(r，t)=Z4(r)e*t(2-8)因系统受到的微扰来自于光波电虚场，而光波

24、波长远大于原子间距。显然，这些理论对可见光、紫外、红外光都是正确的。对X射线的结论则不适用。为了简单起见，先不考虑共振现象，则光波电磁场可以写成如下的形式：EnAef+Ae%(2-9)式中A是复振幅，则光波场与系统的微扰互作用能为H=-EM而M=£>rj，即系统中的电子偶极矩。此时微扰系统的薛定谓方程为J(Ho-E-M)T(r,t)=i/-T(r,t)(2-10)di若k态中未受微扰的系统由平°)(1)=4中(雄-卬r所描述，则(2-10)方程的微扰解为T(r,t)=鹫)(r,t)+K)(r,t)=甲式，0(2-11)将(2T1)式代入(2-10)方程中，略去(EM)

25、K)的二次项，再由(2-6)方程得以下方程：H0-i方一空7(r,t)=(EM)甲骨(2-12)I对(2-12)做求解处理，取(2-13)驾)=exp-i(+牝)t+;exp-i(-)t式中小+1V.AMh不g-AM：=-Yo,卜=-yo(2-14)Mh=M；=;MTkch,%=Q-)受到微扰系统的矩阵元为(2-15)M(t)=/T：MTk(r,t)dr由(2T1)到(2T6)式可得M*(t)=exp(i%t)+expT(%+QX(2-16)+Dkmexp-i(%-旗)t式中的Ckm和Dm1分别为(2-17)Ckm=/i?(AMQM皿MJAM-)(2-18)用一在(dMQMgM.(ATM.)、

26、Jg/rmJg、km，(2-19)乂因为Dknq；很，对k=m的条件可得式中(2-20)(AMQMijMh(AMQ知一在lim-JMkk(t)dti°ro2/k(1)IfMH(t)(t)是实的，它是k态中偶极子动量的期待值，且与入射辐射有相同的时间关系。因此，偶极子M匿(t)辐射的强度仍有以下的经典表示式：(2-21)由(220)式可得3C(2-22)(2-21)式给出了(偶)极矩为M2(t)的偶极子的瑞利散射光强。需特别注意的是：与(t)相反，M*(t)是复的。要找出与(217)式中个别真实偶极子经典辐射相关联的情况，必须用到克莱因（Klein）的结果：若3kmV。，即硒<一

27、3m<0.初态能量小于末态能量，k、m别为初、末态，则分量的辐射发射为零，就等价于真实偶极子的经典辐射，则有Mgm=M1tmexp（-int）+exp（i%t）若3km>0,即u）k3m>。，初态能量小于末态能量，k,m分别为初、末态。为了能应用（2T7）式，必先考虑构成真实偶极子的情况，即=M-exp(_g)+MLe*(i*t)，即+1川+叫代""纽）t（2-23）（2-24）+Dexp-1(6-)t+Dexpi(n-X与（211）式相同，散射光光强由下式决定:工Mkco(t)在对时间取平均时消去交叉项后有京4|MJ2+(%+5cd+(%-叱y(2-25

28、)由辐射发射的原理知：仅对3km>。，3k-3mV0和3k-3m>0的条件下才能产生辐射，下面对（2-25）中各项的意义作以讨论。ikm表示式（225）中的的第一项初态能量（Ek=ho）k）,大于末态能量（%=h3m）。它描述了与外来激光频率3l无关的伴随k-m跃迁的自发辐射，见图2-1（a）o（225）式中的第二项是正常拉曼散射，即EVEk子h3k。末态（m）的能量比始态（k）的能量大，也可以比它小：EQEk末态能量大于始态能量；散射辐射能量小于激光（单）光子能量,即h%>h3L这对应于斯托克斯过程。J/Ek态能量小于始态能量；散射辐射能量大于激光（单）光子能量，即h3sv

29、hM这对应于反斯托克斯过程。需要指出的是：光谱仪所接收到的信号是h3s，见图（2-2）（b）o（2-25）式中的第三项表示伴有两个量子感应发射，即k-m的跃迁。这类发射只有在受激粒子数剧增时才能被观测到。图（22）（c）给出了这一过程的量子跃迁，特别注意激光器的能量与初末能态相应能级的关系。Emk到m的反斯托克斯跃迁伴有两个量子感应发射图22斯托克斯跃迁（Ek），口，2和3是跃迁过程中的中间能级与正常拉曼效应相关的第二项系数C,由（2T8）式确定，其中跃迁矩的求和是从初态k到所有未受微扰系统的r态及跃迁矩从r态回到末态m的求和。并不是说真实的散射过程中存在如上的跃迁过程，完全是因为考虑到数学上

30、处理微扰问题的需要，即受到微扰系统的波函数完全用未受到微扰波函数来表示，也就是在Cg】的表示中包括了矩跃迁Mm和相训的积，而不是单独的某一个跃迁矩。跃迁矩既可以是正的，也可以是负的；既属于不同的态r,也可以附加在另一个态上。当它们所处的态被湮没为另一个使。皿变为零的态时，就不能产生拉曼散射，常称之为拉曼散射的禁戒条件。在对（218）式的求和中，不仅包括了初态k之上的中间态r,也包括了低于k态的任何r态。因为中间态是在吸收了入射光子后产生的高于初态的激发态，所以（218）式求和过程中包含的低于初态k的概念显然是不合适的。还应强调的是：确定自发初末态间跃迁概率的跃迁矩网并没有全部写入品”的表达式。

31、因此，自发辐射（或吸收）与拉曼散射强度间无直接的联系，它们的选择规律也全然不同。由上述的讨论可知：整个空间4兀立体角内拉曼散射强度为八44%=/（%+）QJ（2-26）JV对于k=m,就是瑞利散射强度。4m的分量可以写成如下的形式：(2-27)式中0、Q代表Portor表示中入、散射的偏振方向，式（227）中(C)km=S（MPJMsLi（MQKMPm（如一在）（一在）(2-28)（228）式就是散射张量，通常它是复数，而且是非对称的张量。对于k=m,则有（%）收1号）"（MQ/MJk工（%一在）（如一”）(2-29)(2-30)因此，若（C即）欣是实的，则它是对称的，该结论不仅对3

32、广0的静电场微扰是正确的，而且对哈密顿量为实的系统也是适用的。若不再考虑（220）式中的永久偶极矩Mm,而利用在（2-9）式中的系数关系（C即h二（C加鼠，再从（220）式可得矩阵元表示式，为（M况（2-31）<p<p式中邓，）是实的，就是k态的电极化率，将（227）式代入（226）式中有21km=*（%+砒）2ZG/fanA（232）再用A=|A|6,6是入射光偏振方向的单位失，而入射光光强为Io=：|A|2,则LII有2A%=TT（%+在）>£（%）>=QJ。（2-33）3cp9式中Qm=N3km+/JZZ（C0°）kF被定义为kfU跃迁的拉曼跃

33、迁截面。Q'3c。的量纲是1?,Im和I。的单位分别是erg*/s和erg-s-nr2。若入射光沿B方向偏振，在沿（P方向用分析检测器观察散射光，则单位立体角dQ中散射强度为=/（%+J），仁/w门。（2-34）用仁如卜直接从（231）到（233）式计算散射强度只有在简单系统（如谐振子、自由电子、和某些简单原子）中是可行的。对于分子、晶体及复杂能级系统中，（228）式中出现的受激中间态r都无法进行直接的计算。Plaezek*s近似为直接计算一般结果提供了新的途径“°）。2.2.3Placzek近似首先考虑电子态不发生变化的散射过程，即初、末态是相同的基态。该情况中，仅是振动态

34、发生了变化，而且满足能量守恒条件，即在=千/。因为光散射是由于系统中的电子引起的，光能转移到各个核上，反之亦然，通过核和电子运动间的耦合可以产生拉曼散射效应。假设电子的基态是非简并的，而且原子核被固定在仅产生瑞利（辐）散射的位置，则散射光的强度由电极化率张量（C的）0G确定，（231）式下标中的k=0表示电子的基态。电子的极化率是实的，而且具有对称性：在（229）式中，对具有（C2L固定核系的本征频率外和本征函数影以不仅取决于核的位置，而且也取决于电子极化率a昵，该分量是核组态R的函数，即。阿=«P（P（R）。在以下的假设中，认为振动着的核系统中散射强度是一样的。Plaezek假设有

35、以下三点内容：.电予的基态必须是非简并的：.绝热近似必须是有效的：.激发光源的频率必须小于任何一电子的跃迁频率，但远大于振动的频率，即3«3激光«3电子衣迁若具有振动核系统从k=（Ou）态跃迁到k=（Ou）态，其中0表示电子的基态而以“是振动态。根据近似条件，可认为矩阵元是由电子极化率分量矩阵元陶（R）og；所确定，也就是由下式表示：（GK。广J（R）（R）（R）cIR=（2-35）式中于ou（R）是电子基态0和振动态u的振动波函数。将（235）式代入（234）式，而且略去下指标0,就可以得到伴随Ufl；振动跃迁光散射强度的表示式为L（*）=+）'|（%兀。（2T6

36、）V式中的u、炉分别是两个振动态确定的布局数即u=|%）和u'=|q”而CO,由下式确定：UU方，=E。一E=Z（q-%）方/（2-37）q式中的q式简正坐标的脚码，表示可能出现的声子振动频率的个数。在Placzck近似条件下由n态跃迁到n态的拉曼散射的光强还需做进一步的讨论。若入射光沿着。方向偏振，沿中方向观测散射光，则按（234）式可得到单位立体角内散射光强，为I.=常数3.+J），（a2）（2-38）nnnnxxn-因为极化率0=（1的（R）取决于核的配置R,所以可以将极化率按简正坐标Qs展成以下的级数：a的=aP+Z%6,qQ+RQ+（239）q_qq这是所有量子数假设不变伴有

37、内-】、跃迁的极化率a的矩阵元。若略去式中的QqQq，正比项，则有（%）*=唱"*+"M（Q）、（240）着%=%、（%）*=a*，则有电in=0，这是瑞利散射。由（238）式知可以得到它的散射光强，为鼠=常数（叱）4|（"）门0（2-41）因为在式（240）中略去了（239）式中的二次项，所以瑞利散射光强与温度T无关。以下讨论拉曼散射光强：对第一级拉曼散射有出=为±1，而斯托克斯散射有=nq+lo由此可得（%）%中二=（£严明式4+1产（2-42）4飞对于反斯托克斯散射有1,由此可得（%入中=%&）、.,.1=（/产唯q（%产（243

38、）一q由（238）、（242）和（243）式可得拉曼散射强度表示式。由（237）式可知:对斯托克斯和反斯托克斯散射分布有3mf=-%和%/=Sq相应于这两个散射的光强分别为%*=常数（"-”二。口（1+4）（244）IAnh=常数®+63寸喇4（2-45）式中的%是某一元激发q的平均布居数，若元激发是（热）声子，则它满足玻色一爱因斯坦分布即%=exp（方,/T）-l（2-46）图23给出了声子平均布居数随温度增加而变化的情况：由于温度升高，处于较高能级的布居数也随之增加，分布状况发生了变化，使可能参与跃迁的声子“种类”有所增加，也就是拉曼谱峰中峰的数目增加，由原来的必，与两

39、种，增加到31，C）2»,四种。在小于100K的低温区，仅实现了5（6q=50K）的声子，随着温度增加乂激发J'c）2（6q=100K）,0）3（8q=200K）和0）4（6q=400K）,三种声子且%>S2>03>34。图23玻色.爱因斯坦分布得出的平均布居数随温度的变化(2-47)图23是低、高温能级布居数变化情况的示意。拉曼散射光强与受激后跃迁的元激发数成正比，因此可以得出以下结论：这就是拉曼谱中振动“热带”产生的原因。由（2-46）式可得出参与热激发的声子数，为=5_.Tq力.$toke5【Anti较高温度下，“热”声子数几乎与温度成正比（图23）。

40、由（244）和（245）式有(2-48)解决了经典电磁理论在解释斯托克斯与反斯托克斯散射光强比时的困难gg。2.3拉曼光谱数据分析方法光谱分析技术的数据处理主要涉及两个方面的内容：一是光谱预处理方法的研究，目的是针对特定的样品体系，通过对光谱的适当处理，减弱和消除各种非目标因素对光谱的影响，净化谱图信息，为校正模型的建立和未知样品组成或性质的预测奠定基础：二光谱定性和定最方法的研究，目的在于建立稳定、可靠的定性或定量分析模型，并最终确定未知样品和对其定量冏口可。MATLAB是Mathworks公司开发的一种主要用于数值计算及可视化图形处理的高级计算语言。它将数值分析、矩阵计算、图形、图象处理和

41、仿真等诸多功能集成在一个极易使用的交互式环境之中，为科学研究、工程设计及数据处理和数值计算提供了一种高效率的编程工具”。在这种编程环境下，任何复杂的计算问题及其解的描述均符合人们的科学思维方式和数学表达习惯，而不像Fortran、Basic、C等高级程序设计语言那样难以学习和掌握。MATLAB允许用户根据数值计算的复杂程度，对问题进行分段甚至逐句编程处理，显然，这是与C、Fortran等传统高级语言完全不同的。此外，用MATLAB求解问题一般不需要用户考虑采用何种算法以及怎样具体实现等低层问题，更不必深入了解相应算法的具体细节，因而对用户算法语言方面的要求比较低。仪器采集的原始光谱中除包含与样

42、品组成有关的信息外，同时也包含来自各方面因素所产生的噪音信号。这些噪音信号会对谱图信息产生干扰，有些情况下还非常严重，从而影响校正模型的建立和对未知样品组成或性质的预测。因此，光谱数据预处理主要解决光谱噪音的滤除、数据的筛选、光谱范围的优化及消除其他因素对数据信息的影响，为下步校正模型的建立和未知样品的准确预测打下基础。常用的数据预处理方法有光谱数据的平滑、基线校正、求导、归一化处理等。数据平滑处理数据平滑处理：信号平滑是消除噪声最常用的一种方法，其基本假设是光谱含有的噪声为零均随机白噪声，若多次测量取平均值可降低噪声提高信噪比。平滑处理常用方法有邻近点比较法、移动平均法、指数平均法等口习。(

43、1)邻近点比较法对于许多干扰性的脉冲信号，将每一个数据点和它旁边邻近的数据点的值进行比较可以测得其存在。如果与外近点的数值相差太大，超过给定的阈值，便可认为该数据是一个脉冲干扰，并通过邻近数据点的平均值来取代这一数据点值，就可以把这一干扰脉冲去掉，这样不影响信号的其它部分。在这一数据点处理过程中，需注意选择调节参数，也就是考虑邻近数据点值，以及判断一个数据点和邻近数据点之间不同的阈值。这个阈值一般定义为噪音测量偏差的倍数，以免把必要的有用信号去掉。这一方法有时也称为邻近点平滑法，也叫做单点平滑法。(2)移动平均法由于平滑是通过对信号进行平均而减小噪音，因而多点平滑效果更好。移动平均法是多点平滑

44、中最简单的一种。先选择在数据序列中相邻的有数个数据点，这奇数个数据点即构成一个窗口。计算在窗口内奇数个数据点的平均值，然后用求得的平均值代替奇数个数据点中的中心数据点的数据值，这样我们就得到了数据平滑后的一个新的数据点。接着去掉窗口内的第一个数据点，并添加上紧接着窗口的下一个数据点，形成移动后的一个新窗口，其中的总数据个数不变。同样地，用窗口内的奇数个数据点求平均值，并用它来代替窗口中心的一个数据点.如此移动并平均直到最后。在niatlab中可以调用的平滑函数一般为：yy=smooth(y,span,method)yy=smooth(y,sgolay,degree)其中平滑的方法有：'

45、moving、1owessloess'、*sgolay*rlowess*>Tloess'。可以设置不同的span步长来改变平滑的效果。本次实验处理数据使用的平滑法就是移动平均法。或丫=medfHtl(x,n)中值平滑方法可以通过改变n来得到不同的平滑效果。(3)指数平均法指数平均法是计算在一个具有m个数据点的移动窗口中的各数据点的加权平均.在窗口的最后一个点p1即为要平滑的点，它的权重最大，而前面的每个点分配到的权重依次递减。权重系数由平滑时间常数为T的指数函数eJ1（j标志i前面第j个点，即j=一（m2）一1,0（要平滑的点1的上=0）的形状来决定。pi后点的权重为0,

46、这一过滤函数是用点i前面的点对第i个数据点进行平滑。这一过程和用电子RC滤波器（阻容滤波器）的实时平滑类似。由于该平滑函数是不对称的，故在平滑后的数据中引入了单向失真，这一点也和实时RC滤波器一样。除了获得期望的信噪比降低外，指数平均的结果是峰的最大值下降，同时发生移动。由于用平滑常数T对峰值进行指数平滑和具有时间常数TX=T的仪器测量该峰的效果相同，因此T和峰宽比值函数的强度下降值从实验测量和理论计算都可得到口习。基线校正由于仪器背景、样品粒度和其它因素的影响，近红外分析中常常出现基线漂移和倾斜现象。采用基线校正可有效地消除这些影响。操作时可选用峰谷点扯平、偏移扣减、微分处理和基线倾斜等方法

47、，其中最常用的是一阶微分和二阶微分，但在微分处理时，要注意微分级数和微分数据点的选择。2.3.3据求导处理近红外分析中，对于样品不同组分之间的相互干扰导致吸收光谱谱线重登的现象，可采用求导的方法进行处理。其中常用的是一阶导数和二阶导数。一阶导数表示为：y/=yi+gyi-g二阶导数表示为：y1"=yi+2g2-+%一2g式中：g为光谱间隔，大小可视具体情况设定。对光谱求导一般有两种方法：直接差分法和SavitzkyGolay求导法。对于分辨率高、波长采样点多的光谱，直接差分法求取的导数光谱与实际相差不大，但对于稀疏波长采样点的光谱，该方法所求的导数则存有较大误差，这时可采用Savit

48、zkyGolay卷积求导法计算。导数光谱可有效地消除基线和其它背景的干扰，分辨重叠峰，提高分辨率和灵敏度。但它同时会引入噪声，降低信噪比。在使用时，差分宽度的选择是十分重要的：如果差分宽度太小，噪声会很大，影响所建分析模型的质量;如果差分宽度太大，平滑过度，会失去大量的细节信息。可通过差分宽度与校正标准偏差（SEP）或预测标准偏差（SEC）作图来选取最佳值，一般认为差分宽度不应超过光谱吸收峰半峰宽的L5倍。2.3.4据增强算法在使用多元校正方法建立近红外光谱分析模型时，将光谱的变动（而非光谱的绝对量）与待测性质或组成的变动进行关联。基于以上特点，在建立NIR定量或定性模型前，往往采用一些数据增

49、强算法来消除多余信息，增加样品间的差异，从而提高模型的稳健性和预测能力。常用的算法有均值中心化、标准化和归一化等，其中均值中心化和标准化是最常用的两种方法，在用这两种方法对光谱数据进行处理的同时，往往对性质或组成数据也进行同样的变换。用于消除光程变化或样品稀释等变化对光谱产生的影响。有三种光谱归一化方法：最小/最大归一化、矢量归一化、回零校正。其中常用的是矢量归一化，它是先计算出光谱的y平均值，再用光谱减去该平均值，这样光谱的中值为零，计算所有的y值的平方和，然后用光谱除以该平方和的平方根，结果光谱的矢量归一化是1。回零校正是将光谱减去最小的y值，使得最小y值变为Oo傅里叶变换傅里叶变换FT是

50、一种十分重要的信号处理技术，它能够实现频域函数与时域函数之间的转换，其实质是把原光谱分解成许多不同频率的正弦波的叠加和。根据需要可通过FT对原始光谱数据进行平滑、插值、滤波、拟合及提高分辨率等运算，或用FT频率谱即权系数（傅里叶系数）直接参与模型的建立。在光谱分析中，傅立叶变换可用来对光谱进行平滑去噪、数据压缩以及信息的提取。在matlab中，傅里叶变化的调用函数为：Y=fft(x)w=conv(u,v)y=ifft(X)通过FFT变化后，在卷积其他的滤波函数，通过逆变换得到滤波后的函数。小波变换近年来，小波变换(WT)在信号和图像处理中的应用逐渐广泛和成熟起来。与FT相比，WT具有时频局部化

51、特性。WT理论在80年代得到了迅速发展。其思想起源于工程学、物理学及纯数学领域，被认为是泛函分析、傅里叶分析、样条分析、调和分析、数值分析的完美结品。自1989年以来，由于大量基函数的出现及快速算法的建立，WT在许多领域得到了广泛的应用。90年代初，WT被引入化学领域并形成了化学小波分析。WT能够将化学信号根据频率的不同，分解成多种尺度成分，并对大小不同的尺度成分采取相应粗细的取样步长，从而能够聚焦于信号中的任何部分，分析化学信号的平滑滤噪、去除背景、数据压缩以及重叠信号解析16171819o在matlab中小波变化的主要调用函数为：wname=,sym8r;lev=4;C,LJ=wavede

52、c(X,N/wname,)thr,nkeep=wdcbm(C,L,N);synth=wdencnip('lvd',C,L,wnanieJev,thr,'s');通过选取不同的小波基和分解的层数，进行小波变化。数字滤波按功能分：低通、高通、带通、带阻滤波器：按滤波器的网络结构分：IIR滤波器和FIR滤波器；先按照一定规则将给出的数字滤波器的技术指标转换成模拟低通滤波器的技术指标，据此产生模拟滤波器原型，然后把模拟低通滤波器原型转换成数字滤波器。N阶HR滤波器：MH(z)=T1卜xakz-kk=lN-1阶FIR滤波器：N-1H(z)=Zh(n)zf用MATLAB进行

53、典型的数字滤波器的设计步骤：(1)按一定规则将给出的数字滤波器的技术指标转换成模拟低通滤波器的技术指标：(2)根据转换后的技术指标使用滤波器阶数选择函数，确定最小阶数N和固有频率Wllo(3)运用最小阶数N产生模拟滤波器原型。(4)运用固有频率Wn把模拟低通滤波器原型转换成模拟低通、高通、带通、带阻滤波器。(5)运用冲击响应不变法或双线性变换法来把模拟滤波器转换成数字滤波器。其中FIR数字滤波器的窗函数法设计步骤：用窗函数法设计FIR滤波器时，先根据wc和N求出相应的理想滤波器单位脉冲响应hd(n)；(b)选择合适的窗函数w(n)来截取hd(n)的适当长度(即阶数)，以保证实现要求的阻带衰减：

54、(c)最后得到FIR滤波器单位脉冲响应：h(n)=hd(n)*w(n).第3章常规拉曼检测系统科学技术的进步使拉曼光谱仪得到迅速的发展，人们已研制出适用于多种用途的拉曼光谱仪。从分光系统上来看，现代拉曼光谱仪器一般可分为色散型和以傅里叶变换拉曼光谱仪为代表的非色散型两种。拉曼检测系统的光路主要由光源、外光路系统、样品池、拉曼光谱仪、计算机等五部分组成。光路结构图如图31所示。Lser图3-1常规拉曼检测系统的光路图图中Laser为波长532nm的绿光激光器；L为激光器后的聚光透镜；M为平面镜，作用是使汇聚后的激光直接作用于样品上：F为窄带灌光片；V为校准透镜；Fiber为光纤：S为光谱仪；C为

55、计算机终端。在激光器之后，光谱仪之前的一套光学系统是外光路系统。它的作用是为了有效地利用光源能量、消除瑞利散射光、减少光化学反应和杂散光以及最大程度地收集拉曼散射光。光源门1928年印度物理学家拉曼发现拉曼光谱以后，拉曼光谱分析技术迅速成为了研究分了-结构的主要手段，但是以汞弧灯为光源的拉曼光谱仪由于分辨率低、照相感光纸记录拉曼光谱效果不明显、拉曼散射谱线强度低、曝光时间长等缺点，应用受到很大限制，尤其是红外光谱的出现，使得拉曼光谱的发展进入低潮。上世纪60年代以后，激光技术的问世，拉曼光谱乂进入了新的飞速发展期。世界上各大仪器厂家相继推出了激光拉曼光谱仪（刈。激光具有光强大、方向性强和单色性

56、好的特点，非常适合用于拉曼光谱仪。在拉曼光谱仪中，激光光源必须具有足够的输出功率和稳定性。由于光源辐射功率的波动和电源功率的变化成指数关系，因此往往需要用稳定电源以保证稳定，或者用参比束的方法来减少光源输出的波动对测定所产生的影响。和瑞利散射一样，拉曼散射的强度反比于波长的4次方，因此使用较短波长的激光可以获得较大的散射强度（刈、本论文的拉曼光谱仪采用平均功率为lOOmw,可连续工作，功率稳定，强度变化小，波长为532rlm的激光器。在激光器之后，为了减少杂散光以及最大程度地收集拉曼散射光，在激光器出射部位放置一聚光透镜，使激光束准确地聚焦在样品上。在图3-2中，激光束经透镜L1聚焦后，成为直

57、径很细的细光束。该光束在焦点附近形成一直径为D长度为L的焦柱。图3-2透镜前后激光束变化示意图设未聚焦的激光束直径为d,L1的焦距为f,激光波长为义。则(3-1)(3-2)n4/fD=7rdt16/f这样聚焦后的光束是未聚焦时的千分之一，而照射到样品每单位上的功率提高了约一千倍。滤光片当激发光照射到被测样品上时，除了有很弱的拉曼散射光以外，还有比拉曼散射强倍的瑞利散射光，它严重地影响着拉曼光谱的采集，甚至将拉曼光谱信号湮灭。双联单色仪是消除杂散光的有效办法之一，但这种方法使仪器复杂化，成本也随之增加。因此，反射全息窄带带阻波光片在光谱仪中的应用收到了广泛的重视。其透过率曲线如图3.3所示。这种滤光片是基于反射体全息图的波长选择性制成，具有效率高、带宽窄、单峰特性好、成本低、有效滤除瑞利杂散光，并可同时测到斯托克斯和反斯托克斯谱线等许多优点。在拉曼光谱仪中，主要是利用了全息滤光片的波长选择性、衍射效率高和带宽窄等三个主要特性图3-3滤光片透过率曲线滤光片的波长选择性：体积全息图的再现过程满足布拉格条件，因此，全息滤光片的峰值波长久