激光相位法测距的原理激光相位测距中，把连续的激光进行幅度调制，调制光的光强随时间做周期性变化，测定调制光往返过程中所经过的相位变化即可求出时间和距离。

图.1 相位式激光测距原理示意图如图1所示，设发射处与反射处（提升容器）的距离为x ，激光的速度为c ，激光往返它们之间的时间为t ，则有：cxt 2设调制波频率为f ，从发射到接收间的相位差为 ，则有：N cfxft 242 (2) 其中，N 为完整周期波的个数， 为不足周期波的余相位。

因此可解出：)(2)22(24N N fcN f c f c x(3) 其中，f c L s 2 称为测尺或刻度，N 即是整尺数， 2 N 为余尺。

根据测得的相位移的大小，可知道N 余尺的大小。

而整尺数N 必须通过选择多个合适的测尺频率才能确定，测尺频率的选择是提升容器精确定位的关键因素之一。

多尺测量方法测量正弦信号相移的方法都无法确定相位的整周期数，即不能确定出相位变化中 2的整倍数N ，而只能测量不足 2的相位尾数 ，因此公式（2.3）中的N 值无法确定，使该式产生多个解，距离D 就不能确定。

解决此缺陷的办法是选用一个较低的测尺频率s f ，使其测尺长度s L 稍大于该被测距离，这种状况下不会出现距离的多值解。

但是由于测相系统的测相误差，会导致测距误差，并且选用的s L 越大则测距误差越大。

因此为了得到较高的测距精度而使用较短的测尺长度，即较大的测尺频率s f ，系统的单值测定距离就相应变小。

为了解决长测程和高精度之间的矛盾，一般使用的解决办法是：当待测距离D 大于基本测尺sb L （精测测尺）时，可再使用一个或几个辅助测尺sl L （又叫粗测测尺），然后将各个测尺测得的距离值组合起来得到单一的和精确的距离信息。

由此可见，用一组测尺共同对距离D 进行测量就可以解决距离的多值解，即用短尺保证精度，用长尺保证量程。

这样就解决高精度和长测程的矛盾[4]。

本系统选用10米作为精尺，1000米作为粗尺，带入公式即可求得精尺频率和粗尺频率：精尺频率 MHz L cf 152510（4） 粗尺频率 kHz L cf 150210001000 （5） 其中，光速s m c /1038 。

上面公式计算出的只是个大概的数值，实际上光速要小于s m /1038 ，而且c 还和实际的大气条件（比如矿井温湿度、气体成分、风速等）有关，因此，这些测尺频率需要进一步调整，具体的做法是在现场标定。

混频原理及其在系统中的应用模拟相乘混频器混频是将信号频率由一个量值变换为另一个量值的过程。

如图2.2所示，信号输入和输出的关系分析如下：图2.2 模拟相乘混频器Fig. 2 Frequency mixer设输入信号分别为)cos()(s s sm s t U t U 和)cos()(l l lm l t U t U ，经过模拟乘法器相乘以后得：l s l s lm sm M l s l s lm sm M z t U U K t U U K U)(cos 21)(cos 21（2.6） 由上式可以看出，经过模拟乘法器将两个信号相乘，就实现了两个信号的差频与和频，其中M K 为增益系数。

通过带通滤波器或者低通滤波器后，即可得到差分输出：l s l s lm sm M I t U U K U)(cos 21在相位法测距中使用混频精尺频率15MHz 的正弦信号是中高频信号，对其进行测量是很困难的，这样就要求对信号波形做一定的变化，在保证相位不变的情况下降低信号频率，使后级的模数转换器采样更容易。

在本相位式测距系统中，设由DDS 发出的调制信号和APD 接收到的回波信号分别为1U 、2U ：)cos(11 t U （2.7） )cos(22 t U （2.8）其中f 2 ，f 是精尺频率，其值为15MHz ，此时两路信号的相位差是21 。

另外一个DDS 发出的本振信号)cos(313 t U ，其中112f ，1f 为本振频率，其值为MHz 985.14。

将调制信号1U 与本振信号3U 混频：31131131131)(cos )(cos 21)cos()cos( t t t t U U （2.9） 使用低通滤波器保留其低频 kHz 15的正弦信号，得到：311)(cos t U s （2.10）同理可得回波信号2U 与本振3U 混频后的信号：321)(cos t U l （2.11）此时我们可以得到s U 与l U 的相位差：213231)()( （2.12）由此可见，混频前后相位差不变，信号频率降低到了kHz 15。

同理，对于粗尺频率150kHz ，引入的本振频率为135kHz ，经过上述方法，同样可以在相位差不变的情况下将信号频率降低到kHz 15。

基于快速傅里叶变换的相位测量方法相位法激光测距系统的测量精度主要取决于测相的精度，而传统的测相方式通常大量采用模拟电路，无法解决模拟元件固有的缺点（如温漂、零漂严重，抗干扰能力差等），尤其在煤矿开采现场，不仅环境条件十分恶劣（淋水、粉尘等），而且现场有各种大功率机电设备，有很强的电磁串扰，因而这种采用模拟元器件搭建电路的测相方法在稳定性和可靠性方面都很不理想[4]。

用基于信号频谱分析的鉴相方法，需要对采样的信号进行数字信号处理，这就要求将回波信号这样的模拟量转换为数字量，系统中采用模数转换器（ADC ）实现，采样过程需要遵循一定的条件，采样后的数据进行快速傅里叶变换（FFT ）算法。

采样定理A/D 转换是相位法测距的重要组成部分，是整个数字化处理的基础。

从模拟的连续时域信号得到离散的数字信号应该遵循一定的原则，这就是在数字信号处理领域著名的采样定理。

2 Nyquist 采样定理Nyquist 采样定理是针对基带信号而言的，又称低通采样定理[8]。

设有一个频率带限信号)(t x ，其频带限制在),0(H f 内，如果以不小于H s f f 2 的采样速率对)(t x 进行等间隔采样，得到离散的采样信号)()(s nT x n x 其中s s f T 1 称为采样间隔，则原有信号)(t x 将被得到的采样值)(n x 完全地确定。

采样之后，信号频谱周期化，变为原信号频谱移频后的多个谱叠加，如果原信号)(t x 的频谱如图2.3(a)所示，那么采样后的信号频谱就如图2.3(b)所示[8]。

图3(a) 带限信号)(t x 频谱 3(b) )(t x 经采样后信号频谱Fig. 2.3(a) Spectrum of band-limited signal )(t x Fig. 2.3(b) Spectrum of sampled )(t x由图可见，)(w X s 中包含有)(w X 的频谱成分，如图2.3(b)中虚框部分所示，只要满足H s f f 2 或H s w w 2 （2.13）则虚框部分不会与其它频率部分相混叠。

这时只需要一个带宽大于等于H w 2的低通滤波器，就能滤出原来的信号)(t x 。

Nyquist 采样定理告诉我们，如果以不低于信号最高频率的两倍的采样速率对带限信号进行采样，那么所得到的离散采样值就能准确地确定原信号。

该定理的用意在于，时间上连续的模拟信号可以用时间上离散的采样值来取代，这样就为模拟信号的数字化处理提供了理论依据[8]。

带通信号采样定理带通信号采样定理又称欠采样定理、带通采样定理或中频采样定理[9]。

Nyquist 采样定理只讨论了频谱分布在),0(H f 上的基带信号的采样问题，如果信号的频率分布在某一有限的频带),(H L f f 上时，也需要遵循一定的原则。

当然，根据NyquistHH S HH S S采样定理来进行采样。

但是，当L H H f f B f 时，也就是当信号的最高频率H f 远远大于其信号带宽B 时，如果仍按Nyquist 采样率来采样的话，则其采样率会很高，甚至很难实现，或者后级处理的速度也满足不了要求。

这样的情况下，可以按照带通信号采样定理来采样。

带通信号采样定理：设一个带限信号)(t x ，其频率限制在),(H L f f 内，如果其采样率满足：12)(2n f f f H L s (2.14)式中，n 取满足)(2L H s f f f 的最大正整数（0，1，2，…），则用s f 进行等间隔采样所得到的信号采样值)(s nT x 能准确地确定原信号)(t x [9]。

式（2.14）用带通信号的中心频率0f 和频率带宽B 也可以表示为：1240n f f s (2.15) 式中 20HL f f f，n 取满足B f s 2 （B 为频带宽度）的最大整数。

显然，当20H f f 、H f B 时，取0 n ，式（2.15）就是Nyquist 采样定理，即满足H s f f 2 。

由式（2.15）可见，当确定了频带宽度B ，为了能用最低采样速率即两倍频带宽度速率)2(B f s 对带通信号进行采样，带通信号的中心频率必须满足：B n f 2120（2.16） 或 B n f f H L )12( （2.17） 也即信号的最高和最低频率相加是带宽的整数倍。

带通信号采样前后的频谱示意图如图2.4所示。

(a) 采样前 (b)采样后4 带通信号采样前后的频谱Fig.2.4 Spectrum of bandpass signal before and after samplingL 0H L 0H S S上述带通采样定理适用的先决条件是：只允许在其一个频带上存在信号，但是不允许在许多不同的频带上同时存在信号，否则就将会引起混叠。

但实际的情况是在多个频带上都有信号，为解决这一问题，一般要在采样之前先将信号通过一个带通滤波器，也称抗混叠滤波器[9]。

以上的结论为我们对正弦信号的采样提供了一个总的准则：采样的频率应为信号频率的整数倍，采样的点数应包括整数倍的周期，由于本论文采样过后的信号进行FFT 的处理，基2的FFT 算法要求输入离散数据的点数是2的整数次幂，所以我们这里的采样频率应是正弦信号频率的M 2倍（M 在工程上一般取大于等于2的正整数）[9]。

快速傅里叶变换（FFT ） 离散傅里叶变换（DFT ）离散傅里叶变换（DFT ）的定义为：设x (n )为N 点有限长序列，其DFT 为1)()(N n nk N W n x k X k =0, 1, …, N -1 （2.18）其中：NjN eW 2 ，称为蝶形因子。

一般说来，x (n )和nk N W 都是复数，)(k X 也是复数，所以每计算一个)(k X 值，需要N 次复数乘法和N-1次复数加法。

而)(k X 一共有N 个点（k 从0取到N-1），所以完成整个DFT 运算总共需要2N 次复数乘法及 1 N N 次复数加法[10]。