[发明专利]一种基于Fabry-Perot标准具的透过率曲线参数估计方法

说明书

技术领域

本发明属于光电探测领域，具体说是一种基于Fabry-Perot标准具的透过率曲线参数估计方法。

背景技术

测风激光雷达是以激光为光源向大气发射激光脉冲，接收大气气溶胶粒子和大气分子的后向散射信号，通过分析发射激光的径向多普勒频移来反演风速的。激光器产生的相干激光与大气粒子(米散射)和分子(瑞利散射)的相互作用是用激光雷达测量大气参数的基础。

目前，直接探测测风激光雷达系统中的鉴频器件大多数采用的是Fabry-Perot标准具。Fabry-Perot标准具主要由两块平行放置的平面玻璃板或石英板组成，两板的内表面镀反射膜。为了得到尖锐的条纹，两镀膜面应保持平行，其平行度一般要求达到(1/20～1/100)λ。如果两板之间的光程固定，称为固定式Fabry-Perot标准具；如果两板之间的光程可以调节，则称为可调谐式Fabry-Perot标准具。

激光雷达测风误差的控制一直把重心放在减小激光器产生的激光抖动和跳模上，少有考虑从提高标准具频率选择特性估计上提高风速估计精度。一般情况下，标准具的参数是符合其设计值，并有一定误差。然而在实际使用中，我们发现由于标准具是高精度装置，实验的温度、湿度对其都有一定影响，在使用一定时间时候都需要进行重新估计和校对，同时对于米散射的激光雷达系统，由于激光谱宽因为没有瑞丽效应展宽，透过率函数相对分子散射的情况会变的非常陡峭，所以在使用边缘技术时，对标准具的精度要求非常高。

另一方面，近年来无论是相干还是非相干激光雷达，对于激光器抖动都是采用锁定追踪的方法进行误差消减，但是这种方案是依赖最初的透过率函数参数估计拟合，并以此为基准对激光器的频率进行跟踪。因此，为了提高激光雷达的测风精度必须对标准具的各项参数进行精确的估计，以形成较为准确的透过率曲线，以此形成的频率响应函数才能反演出正确的风速数据。这不仅在理论研究方案有重要的意义，对于推进激光雷达在军事、民用领域的应用都具有重要意义。

发明内容

本发明的发明目的在于提供一种估计精度高，抗激光器频率抖动能力强，实现简单的基于Fabry-Perot标准具的透过率曲线参数估计方法。

本发明的技术目的是采用如下技术方案实现的：

一种基于Fabry-Perot标准具的透过率曲线参数估计方法，其特征在于：其包括按顺序进行的下列步骤：

(1)首先按照较大间隔步长移动标准具腔长，获得整个标准具腔长量程下激光的透过率；

(2)利用透过率曲线的线谱形式，减小噪声影响，然后谱估计的方法，分析其频谱

(3)利用线性最小二乘算法得到每步对应频率；

(4)按照较小步长移动标准具腔长，获得一个峰体的透过率函数；

(5)将步骤(3)得到的每步对应频率带入，利用非线性最小二乘方法对峰值透过率和有效透过率进行估计，绘制标准具透过率函数。

所述步骤(2)将标准具的透过率曲线视为正弦信号，使用的谱估计方法是MUSIC算法，则其基频为每步对应频率和自由谱间距的比值。

所述的步骤(3)从线谱多峰值中求取基频的方法为线性最小二乘法。

所述的步骤(5)使用的透过率曲线公式是带入了步骤(3)中估计得到每步对应频率，参数估计方法为非线性最小二乘法，数据为步骤(4)中小步长，长时间扫描的结果，最后使用估计得到的三个参数，完成标准具在当时工作条件下的透过率曲线。

本发明提供的基于Fabry-Perot标准具的透过率曲线参数估计方法，是利用透过率函数的线谱形式，使用信号处理的方法对其中参数进行估计。其中，使用MUSIC算法和线性最小二乘计算每步对应频率，得到频率参数以后，利用第二次扫描获进行幅度估计，本发明估计精度高，抗激光器频率抖动能力强，实现简单。

附图说明

图1是本发明的流程示意图；

图2是本发明方法中大步长扫描标准具全量程透过率曲线；

图3是使用MUSIC算法得到的透过率曲线的伪谱；

图4是得到的谱峰的谱峰位置；

图5-1是本发明中较小步长扫描脉冲位置的透过率曲线(每步对应频率估计)；

图5-2是本发明中较小步长扫描脉冲位置的透过率曲线(峰值透过率估计)；

图5-3是本发明中较小步长扫描脉冲位置的透过率曲线(有效反射率估计)；

图5-4是本发明中较小步长扫描脉冲位置的透过率曲线(峰值位置估计)；

图6是本发明对四个参数的估计结果在不同信噪比下的均方根误差。