室内空间测量定位系统共模误差分析与消除 下载: 1019次
1 引言
大尺寸空间精密测量技术是制造工艺环节中关键的质量控制手段,随着高端、高性能产品的不断出现,测量对象尺度大幅扩展,需同时兼顾工程测量的大尺度与传统精密测量的高精度,同时,制造数量和产能需求不断增加,测量效率亟待提高[1-2],这对大空间精密测量技术提出了新的要求与挑战。
近年来,分布式测量系统发展迅速,以室内空间测量定位系统(wMPS)为典型代表,wMPS将旋转激光扫描角度传感原理与大尺寸空间整体网络式测量理念相结合,摆脱了以单个测站为核心的测量模式,通过多个测量单元构建出测量尺度可拓展的大型测量网,测量网内可实现并行测量,是一种高精度、高测量效率的三维坐标测量系统[3]。
wMPS在发展与完善中逐渐衍生出一些关键技术,包括系统硬件性能的提升[4-6]、算法的优化与改进[7-9]及辅助测量配件的开发[10-11]等,这些技术可以有效提高系统的整体性能。但是,wMPS信号产生于发射基站内部的激光器,经由物理空间中传播并通过接收器光敏单元接收,最后在某一时间、空间等客观条件下以时间序列形式表现出来,系统通过时间脉冲之间相位关系计算发射基站扫描旋转角度,从而实现时间-空间角度变换。激光信号在各个环节中可能受到包括传播路径中的温湿度梯度变化、软硬件固有缺陷等的多种因素影响,这些误差混叠在一起并一同呈现于时间序列中,最终影响测量的可靠性与准确性。
本文以wMPS时间信号测量原理为基础,分析时间序列产生过程中可能存在的共模误差,利用时间序列统计理论对wMPS的测量信号进行预处理以及时、频域分析,针对不同时间序列中存在的共模噪声进行分析与消除。
2 wMPS测量误差分析
wMPS主要包括发射节点与接收节点两个基本单元,发射站由旋转平台与固定基座组成,两台线型激光器成一定角度安装于旋转平台,发射站向外发射同步脉冲信号以及扫描激光信号。发射站的同步脉冲信号与光平面激光信号均被接收器光敏单元接收,通过硬件计时方式对各个信号到达接收器的时间进行记录,继而可以得到发射站坐标系原点到接收器的水平角
激光信号计时原理如
式中:
当电机调速不稳定时,转速的提高或降低将在时间信息中额外增加或减小时间误差Δ
由(5)式和(6)式可知,两个扫描角具有相同的测量误差2π·Δ
3 共模噪声分析
理论上,wMPS时间序列中存在由多种因素混叠而成的公共噪声,称之为共模误差,通过对时间序列进行分析可有效验证共模误差的存在,其中具体措施包括时间序列预处理、时域分析以及频域分析等。
3.1 预处理
预处理分析可以消除不符合统计规律的异常数据,修正不平稳的wMPS时间序列,保证时频域分析的准确性。
3.1.1 消除粗差
时间序列中不可避免地存在异常粗差点,可通过时间序列的峰度系数判断整体数据的离散程度,并通过3
式中:
数学上已证明,符合正态分布的数据峰度系数近似等于3,若数据中存在粗差,峰度系数将大于3,此时分布曲线中存在明显的厚尾现象,表明数据中有较大的离群点,即粗差[13]。时间序列的峰度系数可以作为粗差判断的指标,当序列中存在粗差,可通过3
3
则可认为数据含有粗大误差,应予以消除,在实际判断过程中,考虑到中位数的抗差作用,可选择使用序列的中位数替代全部数据的均值。
3.1.2 剔除直流量
时间序列可以认为是真值与噪声的混叠序列,真值一般为序列中的直流分量,占据大部分能量,而噪声为交流分量。为了后续更好地对噪声进行处理与分析,需要对时间序列内直流分量进行剔除;同时,接收器在空间中可能会因为蠕变等因素出现一定程度的运动趋势,继而导致时间序列不平稳。因此,在消除数据粗差后,需要在原始时间序列中消除直流趋势项以获得残差时间序列,一般情况下,直流趋势项为一维直线,可以通过最小二乘法拟合获得。
3.2 时域相关性分析
在测量系统中,单个发射站与单个接收器之间存在3个时间序列,扫描角1对应相对时间序列
相关系数指标可以定量分析时间序列之间的相关程度,本文选择皮尔森(Pearson)相关系数[14]进行相关性分析,皮尔森相关性系数定义为
式中:
3.3 频域分析
信号的频域特征分析是通过频率角度分析并认识信号,此方法可以得到信号频率方面的特性,较时域分析方法能够提供更为丰富的信息,一般频域特征分析中,采用傅里叶变换将信号从时域映射至频域。
傅里叶变换[15-16]是以周期信号傅里叶级数为基础,当周期信号满足狄里赫利条件时,可以展开为傅里叶级数,即周期时间信号可以通过一组完备的正交函数集表示为
式中:系数
现实中无法获得长度无限的连续信号,因此需要将连续函数改换为离散数据,且把有限长序列作为周期性离散信号的一个周期进行处理,即离散时间傅里叶变换,通过离散傅里叶变换可进行各个残差时间序列的时频域转换。离散傅里叶变换定义为
式中:
4 共模误差消除
若时间序列中存在共模噪声,则需要选择合适的方法对其进行消除。主成分分析是一种统计分析方法,主要通过正交变换将一组多维度相关随机向量转化为一组对应的不相关的随机向量。在原多维向量集合中,由于冗余信息的存在,向量之间具有相关关系,主成分分析方法通过一定规则的正交转换将原向量集合重新组合成新的无相关关系的随机向量集合,并选取较少的主要的维度信息来尽可能多地完整地反映原随机向量集合,实现维度的降低。这种统计方法已经广泛应用于地球物理信号处理、噪声消除等方面[17-20]。
主成分分析思想可以应用于共模误差消除。假设空间中布置有若干接收器,存在
式中:
主成分分析通过协方差判断向量之间的相关性,而在时间序列预处理过程中,已经通过多项式拟合剔除时间序列的趋势项,因此,每个时间序列的平均值均为0,实现了中心化,则时间序列矩阵
协方差矩阵
式中:
假设存在变换矩阵
若矩阵
因此,在主成分分析中,通过变换矩阵
由(10)式可知,协方差矩阵
计算协方差矩阵
根据特征值的贡献率选择合适的正交基进行矩阵变换,假设选择前
原始数据
若变换矩阵选择全部正交基,则还原后
5 实验验证与分析
5.1 实验平台搭建
搭建实验平台对wMPS时间序列进行分析,实验平台搭建在8 m×6 m×3.5 m的室内场地上,在空间中布置发射站以及接收器,利用3个接收器采集发射站信号,实验平台场景如
5.2 时间序列分析
首先对信号进行预处理,消除较大粗差并去除直流分量,得到最终的残差时间序列,通过皮尔森相关系数对相关性进行验证,计算其中一组发射站与接收器对应的残差时间序列的相关系数,得到时间序列相关系数如
表 1. 时间序列相关系数
Table 1. Correlation parameters of time series
|
由
由此可见,接收器处于不同的空间位置也会影响残差时间序列之间的相关系数。考虑到实验所在的空间具有相同的温度与湿度,故不同位置处相关系数的不一致与传播链路无关。为了进一步判断相关性存在差异的原因,使用新的接收器替换原接收器,且保证新接收器放置位置与原位置保持一致。结果表明,接收器的替换不影响相关系数,相关性存在差异的主要原因为对于接收器来说发射站旋转过程不一致。
表 2. 部分时间序列相关系数
Table 2. Correlation parameters of partial time series
|
将每个信号生成残差时间序列转变至频域,其中接收器1的残差时间序列频谱图如
图 3. 残差时间序列频谱图。 (a)残差时间序列1的幅度谱; (b)残差时间序列2的幅度谱; (c)扫描周期的幅度谱
Fig. 3. Spectra of residual time series. (a) Magnitude spectrum of residual time series 1; (b) magnitude spectrum of residual time series 2; (c) magnitude spectrum of scanning period
图 4. 主成分分析滤波(a)前(b)后结果。(1)~(3)对应接收器1~3
Fig. 4. Results (a) before and (b) after principal component analysis filtering. (1)-(3) corresponding to receivers 1-3
综上所述,在残差时间序列中普遍存在一种同时影响各个时间序列的误差,即共模噪声,其主要误差源来自发射站。
5.3 共模误差消除
针对3个接收器中残差时间序列进行主成分分析滤波以消除共模误差,首先构造原始时间序列矩阵
6 结论
针对时间测量信号中的共模噪声提出一种噪声滤波方法,基于时间序列统计理论对系统时间序列信号时、频域特征进行分析,引入主成分分析滤波方法对测量信号中共模噪声进行消除。实验结果表明,主成分分析滤波方法可在一定程度上降低时间序列噪声。此外,共模误差存在的本质是发射站调速不稳等因素共同导致,可在误差产生的根源进行控制,例如选择高精度的电机伺服系统维持转速稳定;传播路径上由于反光等情况引起的多路径问题,无法通过滤波手段进行去除,需要加入完好性分析算法进行判断。这些问题有待在今后的工作中进一步研究。
[1] 叶声华, 邾继贵, 张滋黎, 等. 大空间坐标尺寸测量研究的现状与发展[J]. 计量学报, 2008, 29(4): 1-6.
[2] FranceschiniF, GalettoM, MaisanoD, et al. Distributed large-scale dimensional metrology[M]. London: Springer London, 2011.
[3] 杨凌辉. 基于光电扫描的大尺度空间坐标测量定位技术研究[D]. 天津: 天津大学, 2010.
Yang LH. Research on large-scale space coordinate measurement location technology based on optical scanning[D]. Tianjin: Tianjin University, 2010.
[4] 端木琼. wMPS系统的硬件平台优化及动态坐标测量关键技术研究[D]. 天津: 天津大学, 2011.
Duan MQ. Research on hardware platform optimization and key technologies of dynamic metrology for wMPS[D]. Tianjin: Tianjin University, 2011.
[5] 刘哲旭. 面向大型复杂现场的测量定位方法研究[D]. 天津: 天津大学, 2014.
Liu ZX. Research on measurement and positioning method in large complex field[D]. Tianjin: Tianjin University, 2014.
[6] 劳达宝. 室内测量定位系统发射站优化及校准技术研究[D]. 天津: 天津大学, 2011.
Lao DB. Research on transmitter optimization and calibration technology of indoor measurement positioning system[D]. Tianjin: Tianjin University, 2011.
[7] 黄喆, 杨凌辉, 赵子越, 等. 基于扩展卡尔曼滤波的光电扫描动态坐标测量算法研究[J]. 激光与光电子学进展, 2016, 53(5): 051201.
[8] 王姣, 杨凌辉, 黄喆, 等. 基于光电扫描/捷联惯导系统的室内组合导航算法研究[J]. 激光与光电子学进展, 2016, 53(10): 101201.
[11] 孙卿, 杨凌辉. 基于测距传感器的WMPS非接触测量方法[J]. 激光技术, 2016, 40(5): 670-675.
[12] 费业泰. 误差理论与数据处理[M]. 7版. 北京: 机械工业出版社, 2015.
Fei YT. Error theory and data processing[M]. 7th ed. Beijing: China Machine Press, 2015.
[13] Decarlo L T. On the meaning and use of kurtosis[J]. Psychological Methods, 1997, 2(3): 292-307.
[14] 王建军, 宋香荣. 现代应用统计学[M]. 北京: 机械工业出版社, 2016.
Wang JJ, Song XY. Modern applied statistics[M]. Beijing: China Machine Press, 2016.
[15] 吴大正. 信号与线性系统分析[M]. 北京: 高等教育出版社, 2005.
Wu DZ. Analysis of signals and linear systems[M]. Beijing: Higher Education Press, 2005.
[16] 周子云. 离散傅里叶级数在五棱镜扫描检测中的应用[J]. 激光与光电子学进展, 2015, 52(7): 071209.
[17] 符养. 中国大陆现今地壳形变与 GPS 坐标时间序列分析[D]. 上海: 中国科学院上海天文台, 2002.
FuY. Present-day crustal deformation in China and GPS-derived coordinate time series analysis[D]. Shanghai: Shanghai Observatory,Chinese Academy of Sciences, 2002.
[18] 袁治灵, 陈俊波, 黄伟源, 等. 基于稳健性主成分分析算法的光学相干层析成像去除散斑噪声的研究[J]. 光学学报, 2018, 38(5): 0511002.
[19] 张欢, 池越, 周亚同, 等. 结合主成分分析的四维块匹配协同滤波三维地震信号去噪[J]. 激光与光电子学进展, 2018, 55(4): 041007.
[20] 于海超, 刘刚, 欧全宏, 等. 月季花叶片的二维相关红外光谱研究[J]. 激光与光电子学进展, 2016, 53(5): 053003.
Article Outline
郭思阳, 林嘉睿, 杨凌辉, 邾继贵. 室内空间测量定位系统共模误差分析与消除[J]. 激光与光电子学进展, 2019, 56(4): 041201. Siyang Guo, Jiarui Lin, Linghui Yang, Jigui Zhu. Analysis and Elimination of Common-Mode Errors in Workspace Measurement and Positioning System[J]. Laser & Optoelectronics Progress, 2019, 56(4): 041201.