啁啾脉冲技术(chirped-pulse amplification, CPA)的出现^[1-2]，极大地提高了激光脉冲的峰值功率。目前激光聚焦功率密度已达到10²⁰W/cm²以上^[3]，在激光等离子体物理、激光加速、强场下团簇光物理、超短X射线和γ射线脉冲等研究领域具有重要的应用^[4-9]。对于飞秒激光脉冲，脉冲宽度是非常重要的激光参量，它决定了激光脉冲的时间特性，由此指导激光脉冲可具体应用在超快领域。除此之外，通过激光能量和脉宽信息可得知激光峰值功率，结合激光峰值功率和激光聚焦焦斑可以获取激光的光强，由此指导激光可应用的强场领域。因此实时、准确、快捷地测量激光脉冲宽度是非常必要的。

对于强场飞秒激光，脉冲的重复频率通常为赫兹量级甚至是单发^[10-11]，对于脉宽低于100fs的脉冲测量常用单次自相关方法^[12]。该方法的基本思路是通过利用线阵或面阵CCD测量激光脉冲的自相关信号，通过自相关信号的宽度信息除以反卷积因子D从而获取激光脉冲宽度信息。传统的单次自相关方法都是通过CCD采集到信号以后，再离线分析数据^[13-14]，这无法实时获取激光的脉冲宽度信息。对于大型强场飞秒激光系统，考虑到激光系统需要在超净恒温环境中运行，且激光系统抽运源的有限寿命等因素，这导致激光运行成本极高。实时获取激光脉冲信息能够让激光系统在优化运行时实时反映激光脉宽状态，从而及时地提供反馈调整激光脉宽，这不仅能最大程度地节约激光的优化调整时间，而且极大地提高激光系统的运行机时，降低运行成本。针对北京大学正在建设的超小型激光粒子加速器(compact laser plasma accelerator，CLAPA)^[15-16]，其基础是200TW的强场钛宝石飞秒激光系统，快速精确地诊断激光脉宽信息可以极大的提高CLAPA中激光系统机时，提高激光的使用效率、降低激光使用成本，从而提高CLAPA系统的运行效率。

本文中通过单次自相关方法，基于LabVIEW实现对CCD采集到的信号进行实时处理，在图像处理时通过减小计算区域以及对图像数据积分的方法降低噪声对数据的影响，结合在线自标定，从而实现了对强场飞秒激光脉冲的实时、准确测量。

图 1是单次自相关仪测量原理图。两束具有一定宽度的激光脉冲在非线性晶体(例如β相偏硼酸钡(BaB₂O₄，BBO))中产生和频光(second-harmonic generation，SHG)信号，激光脉冲的宽度与产生的和频光信号的空间尺寸相关^[17]：

Figure 1.  Principle of single-shot autocorrelator

式中，Δx为和频光信号的空间宽度，τ为激光脉冲宽度，v_g为激光脉冲的群速度，φ为两束激光脉冲夹角的半角。由(1)式可见，激光脉冲的时间宽度通过和频光信号转换到光束的空间尺寸，通过测量产生的和频光信号空间宽度以及参与和频的两脉冲夹角可以获得激光脉冲宽度信息。实际操作中，产生和频的两个激光脉冲的夹角不易准确测量，而且随着光路调节，夹角会有略微变化，因此通常不是直接测量夹角来得到激光脉冲宽度，而是采用自标定的方法。自标定原理如图 2所示。在T₀时刻，两路激光延时相同，给两路激光的其中一路激光脉冲1提供延时t后，两路激光在晶体中重叠的位置将发生变化，因此产生的和频光信号会有位置的偏移，利用偏移量(x₂-x₁)除以延时t，获得和频光信号单位偏移量对应的时间尺度Δt。

Figure 2.  Principle of self-calibration

在干涉仪的基础上设计单次相关仪光路，如图 3所示。入射飞秒激光经过分束片(beam splitter, BS)分成2束激光，其中一束激光经过M₁和M₂反射后入射到M₆镜片上，另一束光经过M₃, M₄, M₅反射后也入射到M₆上，M₁和M₂放置在平移台上，该平移台行程为10cm，精度为1.25μm，可为光路提供精确延时。两束激光经过M₆反射后在晶体BBO中产生和频光信号，和频光信号由面阵CCD(型号Basler acA1600-60gm)接收。通过调节M₅可以方便地使两路激光在BBO晶体中空间重合，调节M₁和M₂放置的平移台可以使两路激光在BBO晶体中时间重合。考虑到色散对飞秒激光的展宽效应，装置中尽可能使用低色散元件：分束片BS厚度为2mm，镜片M₁~M₆均为银镜，BBO厚度为50μm。测量的激光源是CLAPA系统中的掺钛蓝宝石飞秒激光，中心波长为800nm，重复频率为5Hz。

Figure 3.  System of autocorrelator

基于LabVIEW，将单次自相关脉冲宽度的测量集成在如图 4所示的实时测量显示的程序中，具体包括CCD图像显示、自标定以及脉冲宽度计算。程序分为Setup和Data两部分显示界面。其中图 4a是程序Setup界面，主要用于显示利用所搭建的装置进行测量时CCD接收到的信号，CCD采用Basler面阵CCD，每个像素的尺寸是4.5μm ×4.5μm。图 4b是程序Data界面，主要用于图像数据处理、自标定及脉宽计算显示部分。

Figure 4.  Real-time measurement program of sing-shot autocorrelator

在Setup界面中对CCD采集的图像进行显示并对图像计算区域进行选择。在进行脉宽测量时，首先在Setup界面中观察采集到的和频光信号，并调节单次自相关仪装置：通过调节装置中的延时使光信号最强，然后调节装置中的镜片M₆使光信号在CCD中间部分，再调整CCD的姿态使采集到的光信号呈竖直状态。最终在相关仪CCD中采集到的图像如图 4a所示，和频光信号为竖直的长条形状。

为了减小CCD采集信号中噪声的影响，需要将采集到的图像选择进行计算的区域。图 4a中点击Reset ROI，拖动图像中的4条游标，选定将要计算的区域，游标中间的方形区域便是选定的计算区域，该区域应包含采集到的和频光信号，选择好游标的位置后点击Set as RIO确定选择的区域。选定区域后，程序将对选定区域的每一列像素数值进行积分获取自相关曲线，该在Data界面中实时显示，如图 4b所示，其中图中横坐标是选定区域的每一列的坐标值。如该图中进行计算的图像选择了500列像素，纵坐标是每一列像素积分的数值。相对于直接选取某一行的像素作为自相关曲线计算，采取对每列像素积分这样的数据处理能够极大程度地减小光斑强度不均匀对测量精确性的影响。

在Data界面进行自相关曲线的显示、测量的自标定以及脉宽的计算显示。为了获取脉冲宽度的具体数值，需要利用光路中的延时并结合Data界面中自相关曲线的显示对相关仪进行自标定。当移动光路中的平移台为脉冲提供延时，此时图 4a中Setup界面的光信号将左右移动，转换到Data界面中，图 4b中的相关曲线会沿着x轴移动。由于此时Data界面显示的是Setup界面选定区域的图像数据，因此为了能够进行自标定，Setup界面中选定的区域应该尽可能地包含多列数据。具体标定如下：先移动平移台将相关曲线移动到Data显示界面的一端，如图 5a所示，并在软件中的MicroPos1输入此时平移台微分尺的读数5.27，再反向移动光路中的平移台将相关曲线移动到Data显示界面的另一端，如图 5b所示，并在软件中的MicroPos1输入此时平移台微分尺的读数5.05。通过平移台的2次读数可以获取时间延迟t，程序获取当提供延时t时光信号在CCD上移动的像素个数P，利用t除以像素个数P获取单个像素对应的时间尺度Y。如图 4b所示，在Calibocof中显示的是该次实验中CCD一个像素对应的时间尺度Y，为3.67fs。由于作者采用的CCD像素尺寸小，因此极大地缩小了每个像素对应的时间尺度，提高了脉冲测量的精度。

Figure 5.  Self-calibration of single-shot autocorrelator

在完成自标定后，程序计算出自相关曲线半峰全宽的像素个数B，利用B乘以每个像素对应的时间尺度Y来获得这次测量中自相关信号的宽度。利用自相关信号的宽度除以反卷积因子D获得激光脉冲的实际宽度^[18]。图 4b中Gausspulsewidth和Sech2pulse-width分别显示了此次测量中飞秒激光的高斯型和双曲正割型脉冲的脉宽。

基于LabVIEW程序实现了飞秒激光脉冲单次自相关仪的实时测量，对数据进行实时采集、处理，及时获得激光脉宽信息。利用该自相关仪对强场飞秒激光进行了实时测量，准确、快捷地获得脉冲的宽度信息。飞秒激光脉宽实时、精确、便捷的测量更有利于激光系统的应用和发展，为激光的应用提供了可靠的工具。该自相关仪将用于北京大学激光加速器CLAPA中200TW激光系统脉宽的测量研究。