专利名称：基于任意极值的低相干干涉解调方法
技术领域：
本发明涉及一种基于低相干干涉条纹的任意一个波峰或波谷的低相干干涉解调方法，适用于压力，温度，应变和位移传感等领域。
背景技术：
低相干干涉是一种能够高精度测量绝对位移的有效方法，应用于表面三维轮廓检测、透明材料的厚度测量和光学相干层析技术等领域中。通常的做法是通过定位低相干干涉条纹的零光程位置获取测量信息，因此零光程位置定位的准确度和重复性会直接影响解调装置的精度和稳定度。在低相干干涉发展的起步阶段，利用干涉条纹的对比度或调制度来确定零光程差位置是应用最为广泛的一种方法，理论上干涉条纹对比度最大处就是零光程差位置，这种方法简单但是抗噪声能力差，在信噪比比较差的条件下，很难达到很高的精度。随着低相干干涉技术的发展，提取干涉条纹包络并利用最大值定位零光程位置的方法得到重视，提取低相干干涉条纹包络的方法包括最小二乘曲线拟合法、傅里叶变换滤波法、 小波算法和希尔波特算法等，最大值搜索法和重心法被用于直接计算零光程位置。
但是，当低相干干涉条纹畸变较大，信噪比较低时，上述方法零光程位置定位误差较大，导致测量精度受限，难于满足测量精度要求高于0. 5% F. S.的领域。发明内容
本发明的目的在于克服现有技术存在的上述不足，提供一种简单有效的基于任意极值的低相干干涉解调方法。
本发明提供的基于任意极值的低相干干涉解调方法，首先对接收的低相干干涉信号进行自适应i^ourier滤波，然后从滤波后的干涉信号中提取出有效的波峰或波谷以获取干涉条纹的整体平移信息，进而得到F-P传感器的腔长变化实现解调，具体实现步骤为
第1、对干涉条纹信号进行离散Rmrier变换；
第2、根据离散Rmrier变换后的幅频信息确定干涉信号所处的频率段，并以该频率段内幅值最大处对应的频率为中心，通过设置阈值幅度来自动地确定干涉信号的起始频率和截止频率，保持起始频率和截止频率段内对应的幅值和相位不变，其余频率对应的幅值和相位都设置为0 ；
第3、进行离散Rmrier反变换，反变换后得到的复数的实部即为自适应Rmrier 滤波后的滤波信号。
第4、搜寻出滤波后干涉信号所有的极值点，然后通过设置一个上限阈值和下限阈值，以提取出有效的波峰和波谷，其中任意一个有效的波峰或波谷的极值位置都随外界气体压强做单调平移。
第5、对提取出的有效波峰或波谷按照从左到右的顺序进行编号，每一个有效的波峰或波谷对应一个特定的编号，解调时始终追踪其中一个特定编号的波峰或波谷的极值位置，从而实现对干涉条纹的整体平移信息的获取，进而获取F-P传感器的腔长变化实现解调。
第2步所述的自适应Rmrier滤波的阈值幅度的确定方法是阈值幅度的设定依干涉信号的信噪比而不同，但阈值赋值的设定只要能够满足区分出干涉信号的频率和噪声的频率即可。
第4步所述的上限阈值和下线阈值的确定方法是上限阈值的设定需要满足保证高于上限阈值的有效波峰数目不变，下限阈值的设定需要满足保证低于下限阈值的有效波谷数目不变。
当需要更高解调精度时，同时追踪多个编号的波峰和波谷，并进行信号平均处理获取干涉条纹的整体平移信息。
本发明中Fourier滤波的理论依据
在工程实际中，任何一种实际的宽带光源，如常见的发光二极管(LED)、超辐射发光二极管(SLD)，卤素灯等，它们的光谱一般都呈近似高斯分布，为了便于分析，我们假设光源光谱成理想高斯分布，则低相干干涉条纹的光强可以表示为
I(z) = rela(z-z")]2 cos[y (z-z0)](1)
其中，α、β、γ是与光路系统有关的常数，、是低相干干涉信号的零光程差对应在光楔中的坐标。不妨设标准高斯函数f(z) = Yexp[-[az]2]的连续傅立叶变换为 F(co)，F(co)同样是高斯函数，那么根据连续傅立叶变换的时移和频移特性，可以得出⑴ 式的连续傅立叶变换G(co)为( (《) = 1/2八《 + / )^。+ 1/2八《-。(2)
从上式可以看出，干涉信号的傅立叶变换的幅频特性曲线由两个对称的高斯函数组成，函数的对称轴分别位于ω = β和ω = -β处，系统确定后β为一常量，所以理论上两个高斯函数的位置固定不变，只截取正频率部分的高斯函数，即对l/2F( -/ )e—_。进行Rmrier反变换可以得到
g (Z) = f (Z-Z0) exp [j β (ζ_ζ0) ] = Yexp [-[a (ζ-ζ0) ]2] exp [j β (ζ_ζ0) ] (3)
然后对式(3)取实部即得原始干涉信号
I (z) = Re[^(z)] = γcos[y (z - z0)](4)
对式(3)进行取模运算可以得到包络函数ye--z。)]2。
本发明的有益效果及优点在于
1、本发明不需要寻找低相干干涉条纹零光程差位置，色散、光学误差，入射角度误差和噪声等因素会对对称的低相干干涉条纹产生极大畸变，形成的非对称低相干干涉条纹，且非对称性随机变化，零光程差位置判定误差很大，本发明通过跟踪任意一个特定的波峰或波谷，消除了低相干干涉条纹对称性畸变的影响，而且计算更简单，能够实现高速、高精度解调。
2、当需要更高解调精度时，能够通过同时追踪多个编号的波峰和波谷，并进行信号平均处理达到。
图1为空间扫描型低相干干涉光纤传感大气压力解调装置示意图2为实际解调装置中采集的干涉信号的非对称性以及包络法的缺陷分析图，(a)中⑴为实际采集的干涉信号，(2)为经过自适应Rnirier滤波后的滤波信号，(3)为滤波信号的包络；(b)为零光程差位置和低相干干涉包络峰值点之间的偏移示意图3为干涉信号经过离散Rmrier后的部分幅频特性曲线；
图4为本发明中提出的阈值法提取有效波峰和有效波谷示意图5为空间扫描型低相干干涉系统使用的光源光谱图6为压强和探测的峰值点之间的关系曲线，其中(a)采用包络法，(b)采用本发明中提出的基于任意一个有效波峰或波谷的方法；
图7为解调误差曲线，其中(a)为包络法，(b)为本发明中提出的基于任意一个有效波峰或波谷的方法；
图中1-宽带光源，2-光耦合器，3-(F-P)传感器，4_自聚焦准直透镜，5_光楔， 6-线阵CCD。
具体实施方式
实施例
附图1是本发明方法涉及的一个基于任意极值的低相干干涉解调装置图，结合外界大气压力的测量，对本方法进行说明本发明的解调方法通过实验进行了验证，参见附图 2-附图7。
宽带光源1 (其光源光谱见附图5)发出的光导入法布里-珀罗传感器3，从3的两个端面反射的光发生近似双光束干涉，光程差为腔长的两倍，腔长随外界大气压强呈线性变化。被调制过的光信号从耦合器2的出口导出，通过自聚焦准直透镜4并透过光楔5，最终到达线阵CCD 6，光楔5可以看作为空间分布的法布里-珀罗传感器，当5引起的光程差和3引起的光程差相匹配时，会在6相应的局部区域产生明显的低相干干涉条纹。
实验中需要的大气压强通过高精度、高稳定压力源产生，该压力源可以达到10 的控制精度，实验中控制压强以IUa为间隔从60kPa单调递增到200kPa，压强间隔时间为5分钟。线阵CXD的有效像元数为3000点，利用CXD的像元同步脉冲作为NI采集卡 (NI-USB6361)的外部采样时钟，进行模数转换后得到的数字信号由3000个离散数据点构成，每两个相邻数据之间代表一个C⑶像素的宽度7 μ m。附图2的(a)中(1)为下 C⑶输出的一帧未经任何处理的干涉信号，(2)为经过自适应Rmrier滤波后的信号，(3)滤波后的信号的包络，附图2的(b)为实际的零光程差位置和通过包络提取法得到的峰值点位置之间的偏移，这里需要说明，由于光从法布里-珀罗传感器的第二个反射面反射时(空气到硅片，光疏介质到光密介质)会引入半波损失，所以零光程应该是中心波谷位置，从附图2的(b)中可以看出，由于光楔的色散、制造误差，入射角度误差偏差以及噪声等因素，干涉条纹和提取的干涉包络已经失去了明显的对称性，包络峰值位置和零光程位置已发生了较大的偏移，在这种情况下，如果采用包络峰值探测法就会带来很大的解调误差。然而，这种不对称性对单个波峰或波谷的影响很小，而且影响主要体现在波峰或波谷的幅度上，基本上不会改变其极值点的位置，这也是利用单个有效波峰或波谷能够实现对畸变造成的非对称干涉条纹获得高精度解调的原因。
基于任意极值的低相干干涉解调方法具体步骤
首先对经过采集卡得到的数字干涉信号进行Rmrier变换，将时域函数χ (η)转换到频域X(k)，其中Z㈨=，这里N = 3000,离散Rmrier变换后得
权利要求
1.一种基于任意极值的低相干干涉解调方法，其特征在于该方法首先对接收的低相干干涉信号进行自适应i^ourier滤波，然后从滤波后的干涉信号中提取出有效的波峰或波谷以获取干涉条纹的整体平移信息，进而得到F-P传感器的腔长变化实现解调，具体实现步骤为第1、对干涉条纹信号进行离散i^urier变换；第2、根据离散Rmrier变换后正频率的幅频信息确定干涉信号所处的频率段，并以该频率段内幅值最大处对应的频率为中心，通过设置阈值幅度来自动地确定干涉信号的起始频率和截止频率，保持起始频率和截止频率段内对应的幅值和相位不变，其余频率对应的幅值和相位都设置为0 ；第3、进行离散Rmrier反变换，反变换后得到的复数的实部即为自适应Rmrier滤波后的滤波信号；第4、搜寻出滤波后干涉信号所有的极值点，然后通过设置一个上限阈值和下限阈值， 以提取出有效的波峰和波谷，其中任意一个有效的波峰或波谷的极值位置都随外界气体压强做单调平移；第5、对提取出的有效波峰或波谷按照从左到右的顺序进行编号，每一个有效的波峰或波谷对应一个特定的编号，解调时始终追踪其中一个特定编号的波峰或波谷的极值位置， 从而实现对干涉条纹的整体平移信息的获取，进而获取F-P传感器的腔长变化实现解调。
2.根据权利要求书1所述的方法，其特征在于第2步所述的自适应Rmrier滤波中阈值幅度的确定方法是阈值幅度的设定依干涉信号的信噪比而不同，但阈值赋值的设定只要能够满足区分出干涉信号的频率和噪声的频率即可。
3.根据权利要求书1所述的方法，其特征在于第4步所述的上限阈值和下线阈值的确定方法是上限阈值的设定需要满足保证高于上限阈值的有效波峰数目不变，下限阈值的设定需要满足保证低于下限阈值的有效波谷数目不变。
4.根据权利要求书1所述的方法，其特征在于当需要更高解调精度时，同时追踪多个编号的波峰和波谷，并进行信号平均处理以获取干涉条纹的整体平移信息。
全文摘要
基于任意极值的低相干干涉解调方法，属于光纤传感领域，也适用于三维形貌检测以及光学层析技术等领域。本方法选择法布里-珀罗(F-P)传感器作为感受外界气体压强变化的传感元件，光楔作为光程差空间扫描元件，在零光程差的局部区域形成干涉条纹。首先对干涉信号进行自适应Fourier滤波，然后设置合适的上限阈值和下限阈值以提取出有效的波峰和波谷，并保证波峰和波谷的数目不变，对这些波峰和波谷分别按照从左到右的顺序进行编号，解调时只追踪其中任意一个编号的波谷或波峰获取干涉条纹的位移，当需要更高精度时,同时追踪多个编号的波峰和波谷进行信号平均处理。该方法能够对畸变造成的非对称干涉条纹获得高精度解调。
文档编号G01D5/26GK102519498SQ20111040303
公开日2012年6月27日 申请日期2011年12月7日 优先权日2011年12月7日
发明者刘琨, 刘铁根, 孟祥娥, 尹金德, 江俊峰, 王少华 申请人:天津大学