一种高精度船只目标光电定位系统设计与实现

李英杰, 李绍军, 蒋鹏, 卢嘉成, 陈力, 梁冬冬, 周磊

李英杰, 李绍军, 蒋鹏, 卢嘉成, 陈力, 梁冬冬, 周磊. 一种高精度船只目标光电定位系统设计与实现[J]. 红外技术, 2023, 45(10): 1090-1095.
引用本文: 李英杰, 李绍军, 蒋鹏, 卢嘉成, 陈力, 梁冬冬, 周磊. 一种高精度船只目标光电定位系统设计与实现[J]. 红外技术, 2023, 45(10): 1090-1095.
LI Yingjie, LI Shaojun, JIANG Peng, LU Jiacheng, CHEN Li, LIANG Dongdong, ZHOU Lei. Design and Implementation of a High-precision Ship Target Photoelectric Positioning System[J]. Infrared Technology , 2023, 45(10): 1090-1095.
Citation: LI Yingjie, LI Shaojun, JIANG Peng, LU Jiacheng, CHEN Li, LIANG Dongdong, ZHOU Lei. Design and Implementation of a High-precision Ship Target Photoelectric Positioning System[J]. Infrared Technology , 2023, 45(10): 1090-1095.

一种高精度船只目标光电定位系统设计与实现

详细信息
    作者简介:

    李英杰(1987-),男,博士,从事光电技术研究。E-mail:lsj_gdsf@163.com

    通讯作者:

    李绍军(1989-),男,硕士,从事云台伺服控制技术研究。E-mail:lsj_gdsf@163.com

  • 中图分类号: TN06

Design and Implementation of a High-precision Ship Target Photoelectric Positioning System

  • 摘要: 针对跨海桥梁或沿海机场等场所需要实时测算获取移动船只目标绝对位置和高程的现实需要,文中介绍一种高精度光电船只目标定位系统设计与实现方案。通过将高清高精度光学探测组件集成于高精度云台设备内部,可实时解算出目标点的绝对位置信息和高程信息并上报监控系统。结合雷达系统引导可实现无人值守情况下自动测算过往船只经纬度坐标及高程功能。经实际运行测试,该系统展现出良好的可靠性,具有较好的测量精度和应用前景。
    Abstract: Herein, the design and implementation of a high-precision photoelectric ship target positioning system is introduced to obtain the absolute position and elevation of a moving ship target for real-time measurements in places such as sea-crossing bridges or coastal airports. By integrating high-definition and high-precision optical detection components into the high-definition PTZ equipment, the absolute position and altitude information of the target point can be calculated in real time and reported to the monitoring system. Combined with the guidance of a radar system, it can automatically measure the longitude, latitude coordinates, and altitude of passing ships under unattended conditions. An actual operational test showed that the system had good reliability, measurement accuracy, and application prospects.
  • 随着科技的发展,激光技术不仅在测距、遥感、通信等方面得到广泛的应用,而且在军事领域得到各国的重视,各类激光武器相继推出,例如激光制导武器、激光雷达等。激光近感探测根据激光束来感知目标,通过目标的回波信号来确定目标的距离和方位,其特点是方向性强、探测精度高、抗电磁干扰能力突出。战场环境中,烟雾对激光有散射和吸收的作用,从而引起能量的衰减,出现虚警和漏警的问题[1]。因此,对于激光在烟雾环境下后向散射特性的研究十分重要。

    针对该问题,国内外科研人员进行了大量的研究。冯继青等[2]利用比尔朗伯定律和经典扩散方程建立烟雾环境下激光透过率模型,分析不同激光波长的透过率,但是该方法只考虑了单次散射,具有局限性。王红霞等[3]建立模型计算1.06 μm脉冲激光在烟雾中的传输,分析得到透过率与粒子粒径、烟雾厚度的关系,并且数值仿真脉冲激光在烟雾中的时间展宽特性。类成新等[4]研究激光在随机分布的烟尘团簇粒子的衰减特性,分析激光波长、入射角和粒子密度等参数对在烟尘中激光衰减的影响。李晓峰等[5]模拟研究在烟雾环境下不同波长激光在各个复折射率条件下的吸收、衰减和散射效应。Mori等[6]分析了非对称因子和Mie散射系数在烟雾中单次散射的变化特点。孟祥盛[7]利用偏振特性设计一种激光引信,该系统可以降低引信对烟雾后向散射信号的接收能力。陈慧敏等[8]建立烟雾后向散射模型,分析回波特性,将仿真结果与实测数值进行对比,验证模型的准确性。

    本文根据Mie散射理论,运用Monte Carlo方法建立脉冲激光近感探测模型,设置不同距离的大小目标,在无干扰和烟雾干扰条件下仿真905 nm脉冲激光,分析回波波形特征。从而为激光近感探测抗烟雾干扰提供理论基础和新的思路。

    战场上环境十分复杂,爆炸产生的烟雾粒子的主要成分是硫、碳、磷及其混合物。粒子的直径大小与爆炸强度、爆炸物成分和气候条件有关,爆炸产生的烟雾是瞬时的。烟雾也可以看作是气溶胶微粒,不仅爆炸会产生烟雾颗粒,人为释放烟雾气溶胶颗粒对制导武器系统是一种干扰[9]。本文选取发烟材料粒子的粒径大致分布在3~21 μm之间,烟雾粒子粒径分布如图 1所示。

    图  1  烟雾粒子粒径分布
    Figure  1.  Particle size distribution of smoke particles

    Mie散射理论可用于各个方向同性的球体,但是对于形状不规则的粒子同样适用。Mie散射理论是研究大气中的气溶胶微粒与辐射光发生散射的经典理论,其散射的强度与频率二次方成正比,方向性较明显。假设入射光的强度为I0,散射距离为l,则散射光强I[10]

    $$ I{\text{ = }}\frac{{{\lambda ^2}}}{{8{\pi ^2}}}\frac{{{i_1} + {i_2}}}{{{l^2}}}{I_0} $$ (1)

    式中:i1i2为强度函数,表达式为:

    $$ \left\{ \begin{array}{l} {i_1} = {s_1}(m, \theta , \alpha ) \times {s_1}^ * (m, \theta , \alpha ) \hfill \\ {i_2} = {s_2}(m, \theta , \alpha ) \times {s_2}^ * (m, \theta , \alpha ) \hfill \\ \end{array} \right. $$ (2)

    式中:m为散射体相对折射率;θ为散射角;s1s2为散射光振幅函数,s1s2分别为s1s2的共轭函数,散射体尺度参数α的表达式为[11]

    $$ \alpha {\text{ = }}\frac{{2\pi r}}{\lambda } $$ (3)

    式中:r是散射体的半径;λ为入射光波长。散射光振幅函数是无穷级数,可以取表达式的前10项来推演结果。因此,s1s2具体展开式为:

    $$ \left\{ \begin{array}{l} {s_1} = \sum\limits_{k = 1}^\infty {\frac{{2k + 1}}{{k(k + 1)}}[{a_k}{\pi _k} + {b_k}{\tau _k}]} \hfill \\ {s_2} = \sum\limits_{k = 1}^\infty {\frac{{2k + 1}}{{k(k + 1)}}[{a_k}{\tau _k} + {b_k}{\pi _k}]} \hfill \\ \end{array} \right. $$ (4)

    式中:akbk表示为Mie散射系数,该系数和散射体相对折射率m及散射体尺度参数α相关。

    烟雾粒子的散射系数Qsca和消光系数Qext的表达式分别为:

    $$ \left\{ \begin{array}{l} {Q_{{\rm{sca}}}} = \frac{2}{{{\alpha ^2}}}\sum\limits_{k = 1}^\infty {(2k + 1)({{\left| {{a_k}} \right|}^2} + {{\left| {{b_k}} \right|}^2})} \hfill \\ {Q_{{\rm{ext}}}} = \frac{2}{{{\alpha ^2}}}\sum\limits_{k = 1}^\infty {(2k + 1){{\rm{Re}}} ({a_k} + {b_k})} \hfill \\ \end{array} \right. $$ (5)

    不同相对折射率消光系数随尺度参数分布如图 2所示。

    图  2  不同相对折射率消光系数随尺度参数分布
    Figure  2.  Distribution of extinction coefficient of different relative refractive indices with scale parameter

    图 2所示,在选取的3种相对折射率下,消光系数随尺度参数的增加呈振荡衰减分布,最终趋于稳定值。相对折射率越大,震荡幅度越大。

    光子与烟雾粒子发生碰撞后各个方向的散射强度用散射相函数来表示,该函数表达式为:

    $$ P(\theta )=\frac{{\left|{S}_{1}(\theta )\right|}^{2}+{\left|{S}_{2}(\theta )\right|}^{2}}{{\displaystyle \sum _{k=1}^{\infty }(2k+1)({\left|{a}_{k}\right|}^{2}+{\left|{b}_{k}\right|}^{2})}} $$ (6)

    式中:S1(θ)、S2(θ)为散射光振幅函数。单个粒子散射相位函数与散射角关系如图 3所示。

    图  3  粒子散射相位函数与散射角关系
    Figure  3.  Relationship between scattering phase function and scattering angle of particles

    构建本模型的主要思路是将发射的脉冲激光分解成大量光子,根据Mie散射理论和Monte Carlo方法模拟光子在烟雾中的运动轨迹,统计出发生散射后的抵达光电探测器的光子。脉冲激光近感探测模型分为3部分:激光发射模型、激光在烟雾中的传输模型、激光接收模型。

    激光器发出的脉冲激光为高斯脉冲,功率表达式为:

    $$ P(t) = {P_0}\exp [ - \frac{{{{(t - \tau /2)}^2}}}{{{\tau ^2}/4\ln 2}}] $$ (7)

    式中:P0为峰值功率;τ为高斯脉冲持续的时间。光子的发射点选择在激光的束腰处,该位置的光子服从高斯分布,因此可得光子的位置为:

    $$ \left\{ \begin{array}{l} {x_t} = {\omega _0}{\xi _1} \hfill \\ {y_t} = {\omega _0}{\xi _2} \hfill \\ {z_t} = 0 \hfill \\ \end{array} \right. $$ (8)

    式中:$ {\omega _0} = {\left( {\lambda {z_0}/\pi } \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} $为束腰半径;z0为瑞利长度;ξ1ξ2为标准正态分布随机数。光子起始发射方向为:

    $$ \left\{ \begin{array}{l} {u_{xt}} = \sin {\theta _t}\cos {\varphi _t} \hfill \\ {u_{yt}} = \sin {\theta _t}\sin {\varphi _t} \hfill \\ {u_{zt}} = \cos {\theta _t} \hfill \\ \end{array} \right. $$ (9)

    式中:${\theta _t} = \left| {\left( {{\theta _0}/2} \right) \cdot {\zeta _3}} \right|$为光子发射方向的天顶角;θ0为光束发散角;ξ3为标准正态分布随机数;ϕt=2π⋅ξ4为光子发射方向的方位角;ξ4为[0, 1]区间上的均匀分布随机数。

    光子在烟雾环境中会与烟雾粒子发生碰撞,碰撞后光子的能量会发生变化,其变化为[12]

    $$ {E_1}{\text{ = }}\frac{{{Q_{{\rm{sca}}}}}}{{{Q_{{\rm{ext}}}}}}{E_0} $$ (10)

    式中:E0为散射前光子能量;E1为散射后光子能量;QscaQext分别为烟雾粒子的散射系数和消光系数,具体表达式参考1.2节。碰撞后,光子的方向也发生变化,其变化为:

    $$ \left\{\begin{array}{l} u_{x s}^{\prime}=\frac{\sin \theta_{\text {sca }}}{\sqrt{1-u_{z s}^2}}\left(u_{x s} u_{z s} \cos \varphi_{\text {sca }}-u_{y s} \sin \varphi_{\text {sca }}\right)+u_{x s} \cos \theta_{\text {sca }} \\ u_{y s}^{\prime}=\frac{\sin \theta_{\text {sca }}}{\sqrt{1-u_{z s}^2}}\left(u_{y s} u_{z s} \cos \varphi_{\text {sca }}+u_{x s} \sin \varphi_{\text {sca }}\right)+u_{y s} \cos \theta_{\text {sca }} \\ u_{z s}^{\prime}=-\sin \theta_{\text {sca }} \cos \varphi_{\text {sca }} \sqrt{1-u_{z s}^2}+u_{z s} \cos \theta_{\text {sca }} \end{array}\right. $$ (11)

    式中:(uxs, uys, uzs)为散射前的光子移动方向;(uxs′, uys′, uzs′)为散射后的光子移动方向;ϕsca为[0, 2π]均匀分布的散射方位角;θsca为散射天顶角。光子与烟雾粒子发生碰撞后,如果没有消亡(能量小于阈值),则继续朝新的方向移动,移动的距离为:

    $$ \Delta s = - \frac{{\ln \varepsilon }}{{{\mu _t}}} $$ (12)

    式中:ε为[0, 1]区间上均匀分布的随机数;μt为烟雾衰减系数。

    光子离开烟雾环境后,朝接收端光学系统移动,有一定的比例被光电探测器接收。若光子进入接收窗口,则有[13]

    $$ {({x_{\rm{f}}} - {d_{{\rm{tr}}}})^2} + y_{\rm{f}}^2 \leqslant R_{\rm{r}}^2 $$ (13)

    式中:xfyf为光子最后一次散射的位置;dtr为收发光轴间距;Rr为接收端镜头半径。同时,光子在进入接收端光学系统时,入射角需要满足接收视场角要求:

    $$ {\theta _{{\rm{in}}}} \leqslant \frac{{{\theta _{{\rm{view}}}}}}{2} $$ (14)

    式中:θin为光子入射角;θview为接收视场角。若满足上式,光子可看作是被光电探测器成功接收,成为回波光子。

    烟雾环境下脉冲激光近感探测模型仿真流程图如图 4所示。大致流程如下:输入相关参数,对脉冲激光收发系统及烟雾模型初始化,光子与粒子发生碰撞后计算出光子的能量和位置,若光子在烟雾边界内且光子存活,重复碰撞直到光子进入光电探测器或者消失。当最后一个光子完成循环流程,计算出激光回波幅值。

    图  4  激光近感探测模型仿真流程
    Figure  4.  Simulation flowchart of laser proximity detection model

    选取大小两种目标,大目标为武装直升机和小型固定翼飞机。武装直升机机体长12.5 m,宽3.4 m,高3.94 m,主旋翼直径16.35 m;小型固定翼飞机长3.3 m,机身直径0.28 m,机翼长1.56 m,高为0.7 m。激光经过该目标的回波在一个周期内距离变化量大,实验中用反射率为0.9的白板代替;小目标为小尺寸靶弹,长为2 m,直径约为12 cm,激光经过该目标的回波在一个周期内距离变化量小,实验中用反射率为0.3的灰板代替。环境选取无干扰和烟雾干扰两种环境,仿真参数如表 1所示。

    表  1  仿真参数
    Table  1.  Simulation parameters
    Simulation parameters Value
    Laser wavelength/nm 905
    Emission pulse width/ns 30
    Emission beam divergence angle/mrad 5
    Receiving field of view angle/mrad 21
    Launching system diameter/mm 10
    Receiving lens diameter/mm 30
    Transmit-receive spacing/mm 35
    Simulated photon number 106
    Smoke particle size range/μm 3-18
    Smoke complex index 1.75-0.43i
    Target surface Bloom
    Target reflectance 0.3(small target)
    0.9(big target)
    Target distance/m 3(small target)
    7(big target)
    下载: 导出CSV 
    | 显示表格

    图 5可知,取小目标和大目标的距离分别为3 m和7 m,比较小目标和大目标,作用距离增大,探测信号回波的幅值减小,即发射接收系统与目标之间的距离和探测信号回波幅值呈负相关。两者探测回波的前沿上升速率呈递增趋势。

    图  5  无干扰时探测回波信号
    Figure  5.  Detection of echo signal without interference

    图 6可知,在烟雾干扰的环境下,对小目标和大目标取相同质量浓度的烟雾,探测回波信号和图 5相比有了明显的变化。脉冲激光会先探测到烟雾,因为烟雾对激光的反射率低,所以接收信号的幅值相对较小;当脉冲激光穿过烟雾到达目标表面,探测回波幅值相对较大,但是由于烟雾环境中粒子对激光的散射和吸收作用,引起能量的衰减,相比较于无干扰条件下,大小目标回波幅值有所降低。烟雾回波和目标回波的脉冲宽度相对于发射激光波形均有一定的展宽,但是前者的展宽程度大于后者。烟雾回波波形呈现前沿陡峭,后沿平缓的非对称特征,对于大目标而言,作用距离的增加,该特征变化得更加明显。因此激光近感探测系统在探测目标时,如果不加入任何抑制后向散射信号方法,烟雾后向散射信号和目标反射信号将会混合在一起,导致探测系统信噪比降低,进而造成系统虚警、漏警等一系列问题。

    图  6  烟雾干扰时探测回波信号
    Figure  6.  Detection of echo signal when smoke interferes

    本文根据Mie散射理论,运用Monte Carlo方法建立脉冲激光近感探测模型,设置参数,仿真得到大小目标在有无烟雾干扰条件下的回波,分析回波的波形特征,得到如下结论:

    ① 无干扰情况下,发射接收系统与目标之间的距离和探测信号回波幅值呈负相关,目标回波前沿的上升速率均呈递增趋势。

    ② 烟雾干扰情况下,脉冲激光会先探测到烟雾回波后探测到目标回波且烟雾回波幅值小于目标回波幅值。烟雾回波和目标回波的脉冲宽度相对于发射激光波形均有一定的展宽,但前者的展宽程度要大于后者,烟雾回波波形呈现前沿陡峭,后沿平缓的非对称特征,对于大目标而言,作用距离的增加,该特征变化得更加明显。

  • 图  1   目标定位系统示意

    Figure  1.   Target positioning system diagram

    图  2   目标位置信息解算流程

    Figure  2.   Target location information calculation process

    图  3   目标与云台坐标系关系

    Figure  3.   The relationship between the target and PTZ coordinate system

    图  4   目标与偏移后的云台坐标系关系

    Figure  4.   The relationship between the target and offset PTZ coordinate system

    图  5   云台伺服系统方案

    Figure  5.   PTZ servo system scheme

    图  6   速度环LADRC控制器结构

    Figure  6.   Speed loop LADRC controller structure

    图  7   速度环响应曲线对比

    Figure  7.   Comparison of speed loop response curves Angle/(°)

    图  8   位置环响应曲线对比

    Figure  8.   Comparison of position loop response curves

    图  9   系统工作流程

    Figure  9.   System workflow diagram

    表  1   电机和编码器参数

    Table  1   Motor and encoder parameter

    J180LWX005 motor Value EAC90F-0M21S encoder Value
    Rated torque/(N⋅m) 8 Resolution 21 bits
    Rated current/A 4 Accuracy < 30″
    No-load speed/rpm 110 Maximum speed/rpm 6000
    Torque coefficient/(N⋅m/A) 2
    Phase resistor/Ω 1.6
    Phase inductance/mH 4.5
    下载: 导出CSV

    表  2   目标经纬度测量对照

    Table  2   Target longitude and latitude measurement comparison

    Actual value of target longitude and latitude/° Measurement value of target longitude and latitude/° Distance/m Deviation/m
    Longitude: 118.7693622; Latitude: 31.87598539 Longitude: 118.7693379; Latitude: 31.8759465 2638 4.77
    Longitude: 118.7841865; Latitude: 31.896963 Longitude: 118.7841426; Latitude: 31.8969342 3825 4.86
    Longitude: 118.7775752; Latitude: 31.9166531 Longitude: 118.777571; Latitude: 31.9166863 6025 3.71
    Longitude: 118.86845383; Latitude: 31.89987727 Longitude: 118.868425; Latitude: 31.8998692 8224 2.84
    下载: 导出CSV

    表  3   目标高程测量对照

    Table  3   Target height measurement comparison

    Actual value of target height/m Measurement value/m Distance/m Deviation/m
    87.944 88.2 2638 0.256
    113.346 112.4 3825 -0.946
    116.26 115.2 6025 -1.06
    230.2 228.7 8224 -1.5
    下载: 导出CSV
  • [1] 张艳军. 通航桥梁对船撞击动力响应研究[J]. 铁道工程学报, 2016(10): 77-82. https://www.cnki.com.cn/Article/CJFDTOTAL-TDGC201610017.htm

    ZHANG Yanjun. Research on the dynamic response of isolated bridge to barge impact[J]. Journal of Railway Engineering Society, 2016(10): 77-82. https://www.cnki.com.cn/Article/CJFDTOTAL-TDGC201610017.htm

    [2] 张攀科. 水上跑道侵入风险演化机理及预警模型研究[D]. 武汉: 武汉理工大学, 2018.

    ZHANG Panke. Study on Evolution Mechanism and the Early Warning Model of Water Runway Incursion Risks[D]. Wuhan: Wuhan University of Technology, 2018.

    [3] 王丹, 董亚力. 桥区水域船舶安全控制高度测量系统设计[J]. 船舶科学技术, 2017, 39(11A): 34-36. https://www.cnki.com.cn/Article/CJFDTOTAL-JCKX201722013.htm

    WANG Dan, DONG Yali. Design of height measurement system for ship navigation safety control in bridge area[J]. Ship Science and Technology, 2017, 39(11A): 34-36. https://www.cnki.com.cn/Article/CJFDTOTAL-JCKX201722013.htm

    [4] 王冬华, 吴壮志. 海防监控系统中雷达引导功能的设计与实现[J]. 计算机工程与科学, 2008, 30(2): 42-45. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJK200802015.htm

    WANG Donghua, WU Zhuangzhi. Design and implementation of the radar guidance in the coastal defense video surveillance systems[J]. Computer Engineering & Science, 2008, 30(2): 42-45. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJK200802015.htm

    [5] 岳尚武, 季诚胜, 孙德新. 永磁同步电机伺服系统控制中的自抗扰控制策略[J]. 红外技术, 2020, 42(2): 121-126. http://hwjs.nvir.cn/article/id/hwjs202002003

    YUE Shangwu, JI Chengsheng, SUN Dexin. Auto disturbance rejection control strategy in servo system controlling with permanent magnet synchronous motor[J]. Infrared Technology, 2020, 42(2): 121-126. http://hwjs.nvir.cn/article/id/hwjs202002003

    [6] 朱斌. 自抗扰入门[M]. 北京: 北京航空航天大学出版社, 2017.

    ZHU Bin. Introduction to Automatic Disturbance Rejection Control[M]. Beijing: Beijing University of Aeronautics and Astronautics Press, 2017.

    [7] 曾岳南, 周斌, 郑雷, 等. 永磁同步电机一阶线性自抗扰控制器的设计[J]. 控制工程, 2017, 24(9): 1818-1822. https://www.cnki.com.cn/Article/CJFDTOTAL-JZDF201709013.htm

    ZENG Yuenan, ZHOU Bin, ZHENG Lei, et al. Design of 1st-order linear active disturbance rejection controller for PMSMs[J]. Control Engineering of China, 2017, 24(9): 1818-1822. https://www.cnki.com.cn/Article/CJFDTOTAL-JZDF201709013.htm

  • 期刊类型引用(2)

    1. 贾云娟,许艳玲,杨君霞. 激光脉冲回波信号特性的检测模型. 激光杂志. 2025(02): 245-250 . 百度学术
    2. 陈印,李翔,张仁永,孙燕斌. 基于谱特征提取的光信号特性分析. 激光杂志. 2025(02): 204-209 . 百度学术

    其他类型引用(0)

图(9)  /  表(3)
计量
  • 文章访问数:  151
  • HTML全文浏览量:  45
  • PDF下载量:  44
  • 被引次数: 2
出版历程
  • 收稿日期:  2022-09-02
  • 修回日期:  2022-09-24
  • 刊出日期:  2023-10-19

目录

/

返回文章
返回
x 关闭 永久关闭

尊敬的专家、作者、读者:

端午节期间因系统维护,《红外技术》网站(hwjs.nvir.cn)将于2024年6月7日20:00-6月10日关闭。关闭期间,您将暂时无法访问《红外技术》网站和登录投审稿系统,给您带来不便敬请谅解!

预计6月11日正常恢复《红外技术》网站及投审稿系统的服务。您如有任何问题,可发送邮件至编辑部邮箱(irtek@china.com)与我们联系。

感谢您对本刊的支持!

《红外技术》编辑部

2024年6月6日