高速运动目标深海落点测量精度仿真分析 *
张志伟1,2,孙 翱1,张 旭1,辛 健1
(1.中国人民解放军91550部队43分队,辽宁 大连 116023;2.中国人民解放军92493部队博士后科研工作站,辽宁 大连 116023)
摘 要: 为深入研究基于TDOA体制下高速运动目标深海落点测量基站布设方式、目标落点位置、参考基站选取等变化对落点测量精度的影响规律,基于泰勒级数展开法建立了TDOA体制下无源定位测量模型,以站址误差、时延误差、声速误差等作为主要误差源变量,采用了蒙特卡洛数值模拟方法通过大子样抽样将呈正态分布的随机数代入定位测量模型,复现多基站条件下时延差定位测量精度特性。分析结果表明,在综合误差源影响下,测量模型能够对基站布设、落点位置、参考基站选取等变化产生合理响应,采用蒙特卡洛数值模拟方法能够有效估计深海高速运动目标落点位置测量精度,可为多基站无源定位测量布站设计与精度评定问题提供参考。
关键词: TDOA体制;蒙特卡洛方法;基站布设;落点位置;参考基站;精度分析
随着我国海洋研究和海洋探测从近海浅水区域向远洋深水区域拓展,对高速运动目标深海落点高精度定位测量需求日益凸显。由于海水对光波和无线电波的强烈吸收特性,无法满足水下目标探测、导航和定位需求。相比之下,声波在水中的传播性能要好得多,因此声波成为水下探测、识别与信息传递的主要手段[1]。利用目标入水声对落点进行定位属于水下无源定位问题,若已知水声信道和水听器阵的接收信息,就可推算声源信息[2]。常用方法包括多基站交会定位、目标运动分析方法(TMA)和匹配场处理技术(MFP)等[3-4]。利用多个基站进行交会测量,通常采用纯方位(bearings-only,BO)交会定位、时间定位(time of arrival,TOA)或时延差(time difference of arrival,TDOA)定位[5-6]。随着深海测量范围的增加,方位交会因误差随距离增大较快变得不适用;TOA方式对时间同步要求较高,需声源与阵元之间严格的时间同步,在无源定位中是不可行的[7],而TDOA主要是通过测量待测声源与各阵元之间的到达时间差而实现定位,仅需要各阵元之间的时间同步,在海面声目标有显著识别特征条件下可通过广义互相关等方法获得较高的时延估计精度,特别对于实时声信号有较快、较强的处理能力[8-10]。因此,基于TDOA体制多站时差定位方法更加适用于深海大区域入水目标落点测量。
目前，关于基阵构型与定位精度关系等问题研究较多[11-17],但对于TDOA体制下深海目标落点位置变化、参考基站选取等对落点定位精度影响规律缺乏一定认识。高速运动目标深海落点测量中涉及到的误差源主要包括站址误差、时延误差、声速误差、水深误差等,虽然通过某些技术手段可以将上述误差降到最低,但是测量误差始终存在。此外,基站布设方式、目标落点位置及参考基站选取等变化对落点测量精度也存在一定影响。
针对上述问题,本文拟在各种测量误差源综合影响下,通过分析基站布设方式、目标落点位置、参考基站选取等变化对落点位置测量精度影响规律,以期为高速运动目标深海落点测量系统设计及相关领域的工程应用提供理论支撑。
1 无源定位测量模型建立
基于TDOA体制的多站时差定位方法本质上是利用一组双曲线交会的方法确定声源位置,即利用每一组时延差测量值确定一条双曲线,最后由多组双曲线相交确定待测声源位置。图1为基于TDOA体制的无源定位示意图,其中黑色圆点为多个接收基阵,红色三角形符号表示双曲交会得到的声源位置。
互联网企业只有建立科学有效的薪酬体系，才能充分发挥薪酬的激励作用，从而提高员工的工作积极性，增强员工履行职责的责任心。公司对员工付出劳动的回报主要体现在其薪酬体系上，通过发放工资、奖金、津贴等来满足员工的物质需求，从而增强员工对其工作的认可程度，降低员工的跳槽率。
因此，本文中利用磁控溅射法制备的Ce3+:YAG荧光薄膜，产生的上转换荧光的发光机理为双光子同时吸收。即Ce3+的外层4f基态电子在同时吸收两个光子的能量后，直接跃迁到5d激发态能级上, 然后通过无辐射弛豫到5d轨道的低能级，最终通过辐射弛豫回基态能级，并产生Ce3+离子的上转换荧光。
图1 基于TDOA体制的无源定位示意图
2)当基站布局完全对称时,如图1中的4基站和8基站,各坐标点X 方向和Y 方向定位精度相当;
(3) 对于溶洞被淤泥或其他软弱土体充填的情况，利用泡沫混凝土作桩，穿越溶洞，与桩间土形成复合地基共同承担荷载。桩长增加，桩侧摩阻力增大，单桩承载力提高，地基承载力得到提高。但桩长继续增加，地基承载力提高有限，存在“有效桩长”。
(1)
设在T 0时刻,待测目标发出声信号,平均声速为C m ,到达各基站时间测量值可由以下公式求得
血常规检测是非常重要的检测指标,给患者进行血常规检查的时候,血细胞检测仪是非常重要的设备,能够让患者短时间内获得较准确的结果[1]。血涂片分析使用和研究当前并没得到足够的重视,临床中误检和漏检的问题仍然存在。血常规检验是可以对早期病变进行发现的一种有效的检测方式,对骨髓造血、血循环等等问题可以尽早的发现。部分仪器会将小型粒细胞判断为中间细胞,把大淋巴细胞当成粒细胞,因此而引发了检验的结果偏差。因此临床中对于血常规检测的重视度必须足够高,尤其是血涂片分析,能够对血细胞指标进行观察,此次我院就血常规检验中使用血涂片检验的临床价值进行探讨分析,现根据研究开展以下报道。
牛血清白蛋白(BSA)，电泳纯（进口分装）；人血清白蛋白(HSA)，电泳纯（进口分装）；氯霉素(CAP)，琥珀酸酐(HS，化学纯），三正丁胺（化学纯），氯甲酸乙酯（化学纯），其他试剂均为AR级。
(2)
以目标落点测量区域中心为原点，建立站心直角坐标系,Y 轴指向真北,X 轴指向东,Z 轴垂直于OXY 平面向上为正,与X 轴、Y 轴构成右手坐标系,测量区域范围为20 km×20 km。将误差源参数设置如下:基站站址误差为20 m,水深测量误差为1.0 m,声速测量误差为2 m/s,时延误差设为200 ms,基站守时误差通常小于1 μs,相对于时延估计误差可忽略。
水资源在日常生活中和人们的利益密切关联，且其是生命之源。水文工作者需牢牢把握水质检测关卡，保证用水安全，从而实现对水中微生物滋生导致的传染病爆发、疾病流行，微生物实验室检验及管理，且相对于化学检测来说，而水质检测的质控、处理，以及采集需注重多环节。随着新出台及实施的《生活饮用水标准检验方法》[1]，不断提升了微生物检测的要求，同时获得人们更多的思考及关注度。结果的分析，以及样品采集、检测、管理等多环节，均包含于水质监测微生物实验室样品管理及检测中，彼此之间互为依靠，紧密联系，但是当有一些不足或者是缺陷出现的时候，则会对最终的检验结果产生极大的影响。
(3)
将式(3)代入式(2),可得
（2018理数20，12分）某工厂的某种产品成箱包装，每箱200件，每一箱产品在交付用户之前要对产品作检验，如检验出不合格品，则更换为合格品．检验时，先从这箱产品中任取20件作检验，再根据检验结果决定是否对余下的所有产品作检验．设每件产品为不合格品的概率都为p(0
(4)
设
为初始位置至各基站到达时延,由于待测声源发出信号时刻T 0无法确定,可利用多基站TDOA定位体制,即通过测量声信号到达各个基站的时延差确定待测目标落点位置,由式(4)可得到达时延差为
(5)
上式可进一步改写成如下形式
ΔT =B ΔX
(6)
其中,
由最小二乘原理可得[18]
(7)
B =
(8)
紫趾综合征是华法林的罕见不良反应，属于华法林引起的皮肤坏死的一种病症，其发生率在0.01%～0.1%。起病初期，足趾患部皮肤感觉异常，并可伴有灼烧样疼痛，随着病程进展，疼痛可能加剧，足趾患部皮肤出现网状红（紫）斑，逐渐破溃，最终足趾患部皮肤出现坏疽[1]。该病起病隐匿，且预后较差，接近一半的患者需要外科手术治疗[2]。今有1例如前述状况，临床药师参与紫趾综合征的治疗，获得明显好转。
ΔX =(B TB )-1B TΔT
(9)
(10)
为保证计算精度,需进行多次迭代,在每一次递归中通过求解TDOA测量误差的局部最小二乘解来改进估计位置。当
(设门限ε 为一小正数)时,迭代停止,此时
即为待测点最终估值。
2 蒙特卡洛数值仿真过程
1)假设高速入水目标落点的真实坐标值为X =[x ,y ,z ]T,在理论值附近选取迭代初值,即分别以x ,y ,z 为均值,方差为
的正态分布
生成迭代坐标初值X 0=[x 0,y 0,z 0]T;
2)设n 个接收基站的真实坐标值为X i =[x i ,y i ,z i ]T,待测目标落点至各基站的真实距离为R i =‖X -X i ‖,若给出平均声速
可得到目标至各基站的真实传播时延为T i ,按
正态分布生成到达时刻测量值
进一步得到待测点到各基站与到参考基站的测量时延差为
3)分别按
正态分布生成随机误差,并叠加至X i 上,得到基站坐标测量值为
按
正态分布生成平均声速测量值
（45）疏叶假护蒴苔Metacalypogeia alternifolia（Nees.） Grolle 余夏君等（2018）
4)根据
得到初始位置到各基站的距离为
结合
可求得到达时刻初始值
进一步得到初始时延差
5)由声源初始位置X 0、基站坐标
平均声速
初始位置到各基站距离
按式(8)计算得到系数矩阵B ,由测量时延差和初始时延差得到观测向量ΔT ;
6)由公式(9)、(10)计算得出第一次迭代后的目标位置
用替换X 0。重复上述迭代过程,当时,停止迭代,此时即为目标位置最终解;
通过分析表1中数据可得如下结论:
3 落点测量精度仿真分析
将式(1)进行线性化处理,即在待测点附近选取某一初始位置X 0=(x 0,y 0,z 0)T进行泰勒级数展开,忽略二阶及以上分量，只保留前两项,利用该值模拟高斯-马尔柯夫定理，有
3.1 基站布设方式对定位精度影响
参考基站选择[10 000 m, 3 000 m,-10 000 m],基站个数选择3～8时,观察落点定位精度随基站数变化规律。图2为基站布设示意图,表1为选择不同基站数时部分坐标点定位精度变化情况。
本次调查样本男女比例基本均衡，平均年龄约为23.7岁，月消费水平从800到10 000元不等，集中于2 500元以下；多为来自东南亚、南亚的留学生群体，有极小部分华裔；来华时长从1个月到7年不等，以小于1年和2.5～4年居多，调查样本的基本情况如表1所示。
7)根据蒙特卡洛数值方法,将上述步骤重复抽样试验10 000次,并对所有测量结果按照正态分布进行统计分析,即可得到X 方向定位误差σ X 、Y 方向定位误差σ Y 及总位置误差σ 总
2.面向大数据精准营销，重构企业信息化系统。汽车经销商信息化系统在初期往往缺乏全景规划，一个企业建立多套系统，包括ERP 、CRM、OA、财务、人力资源、决策平台等，数据在每个系统里面孤立流动。通过构建以客户为核心的全生命周期的大数据精准营销平台，实现技术融合、数据融合、业务融合，为数据的共享、流通、挖掘和应用提供基础，才能不同系统、不同技术、不同部门间的数据壁垒，提高整体的运行效率。在系统建设中前，要首先对自身数字化能力进行诊断，先易后难，先基础后全面。特别是要根据企业规模和自身信息化团队的能力水平，以及对需求和方案的理解，选择自建、第三方服务商定制开发或成型整体解决方案的快速应用。
1)随着基站个数的增加,总体定位精度得到了提高。3基站时定位精度较差,且X 方向定位精度要明显低于Y 方向定位精度,4基站以上定位精度提高幅度不大;
假设在高速运动目标落点海域一定范围内布设n (n ≥3)个基站，且位置坐标X i 已知,X i =(x i ,y i ,z i )T (i =1,2,…,n ),待测目标位置设为X =(x ,y ,z )T ,T i 为待测目标到达各阵元的时间测量值,待测目标至各基站距离为
3)从5基站到6基站,Y 方向定位精度基本保持不变,X 方向双曲交会性能提高,故X 方向精度提高明显;从7基站到8基站,X 方向定位精度基本保持不变,Y 方向双曲交会性能提高,故Y 方向精度提高明显。5基站是在4基站的基础上,通过在X 方向增加了一个基站,由于在X 方向上双曲线得到拉伸,交会性能下降,故Y 方向精度要优于X 方向;同理7基站中,由于Y 方向上双曲线得到拉伸,故X 方向精度要优于Y 方向。
图2 基站布设示意图
表1 选择不同基站个数时目标点的定位精度
3.2 目标落点位置对定位精度影响
当四个基站A、B、C、D位置确定后,在定位解算过程中,选择A点为参考基站,各误差源参数设置保持不变。当待测目标位于测量区域不同位置时,如图3所示,按照沿X 轴方向、沿Y 轴方向及沿两个对角线方向分别考察定位精度的变化情况。表2为待测目标沿坐标轴及对角线方向分布时定位精度变化情况。
图3 待测目标位置分布
通过分析表2中数据可得到如下结论:
1)沿X 轴方向分布的坐标点定位误差与沿Y 轴方向分布的坐标点定位误差关于对角线Ⅱ对称;沿坐标轴及对角线方向分布的坐标点定位总误差先减小后增大,在坐标(0,0)处精度最高;
2)沿直线Y =0,X 方向误差先减小后增大,Y 方向误差在靠近基线AD 方向误差较小,如图2中的4、5两点。由于沿X 方向及Y 方向各坐标点定位误差关于对角线Ⅱ对称,故沿直线X =0,Y 方向误差先减小后增大,X 方向误差在靠近基线AC 方向误差较小;
3)沿对角线Ⅱ方向上各坐标点X 和Y 方向定位精度相当,对角线Ⅰ的点定位误差关于对角线Ⅱ对称。
具体方法为，采用比正常安装支架多上两架支架的方法，让过渡液压支架最外侧的两架安装到运料巷位置，然后用增高装置改造这两架支架，提高支护高度，其余支架作用不变。为保证支架增高装置的稳定性，设计高度不能过高，改造后的过渡支架的支护高度为4.3m，92623工作面煤厚为5.5m，因此还需要在支架增高装置上方摆设木垛进行接顶。在拉架过程中随采高变化增减木垛高度，并及时补充拉架过程中掉落的部分木料，视情况使用长把工具回收掉落在支架后方的木料。为确保顶板的支护强度，必要时可以挑长钢梁加强增高装置至外段超前支护的顶板，确保摆木垛的人员的顶板安全。
上述分析仅仅针对沿测量区域四个方向坐标点的误差分布规律,并不能代表整个测量区域的误差分布规律,因而必须从整个测量区域来分析不同落点误差分布规律。
3.3 参考基站选取对定位精度影响
在建立TDOA体制定位模型过程中,需要获取各基站至参考基站的到达时延差,因此建模过程涉及到参考基站选取问题。当选取的参考基站变化时,为便于观察区域内落点处于不同位置定位精度变化规律,可将测量区域进行网格化处理,具体方法如下:
选取横坐标范围[-8 000 m,8 000 m]、纵坐标范围[-8 000 m,8 000 m]的区域,按1 000 m×1 000 m间距将测量区域进行网格化处理,将该区域划分为289个坐标点。当落点处于不同位置时,区域内各坐标点X 方向、Y 方向及总位置误差分布情况如图4～图7所示,表3为部分坐标点定位误差精度统计情况。
党的十八届四中全会决定设立“国家宪法日”，建立宪法宣誓制度，充分显示党中央尊重宪法、遵守宪法的决心，使得全国人民认识到宪法的重要和至高无上的权威。依法治国、宪法先行，全国各地修建宪法公园，宪法进学校、进课堂，进行宪法知识讲座，全民学宪法、领导干部带头学宪法，使得宪法信仰的旗帜再次高高飘扬在共和国的上空。
由图4～7结合表3可知,当选择不同参考基站且目标落点处于不同位置时,落点的X 方向、Y 方向及总位置误差分布存在如下规律:
1)由于基站布设存在对称性,定位误差分布也存在明显对称性。当选择某一参考基站并得到该区域落点测量误差分布规律,根据误差分布的对称性,很容易得到选择其他参考基站时误差分布规律。上述分析说明了定位模型能够对参考基站选取及目标落点位置变化产生合理响应,进一步验证了模型的准确性;
表2 沿坐标轴及对角线方向分布时目标点的定位精度
图4 参考基站选择A点时定位误差变化情况
图5 参考基站选择B点时定位误差变化情况
图6 参考基站选择C点时定位误差变化情况
图7 参考基站选择D点时定位误差变化情况
表3 部分坐标点定位精度统计
2)当目标落点位于测量区域中心附近时,参考基站选取对定位精度造成的影响很小,例如坐标(0 m,0 m),不同参考基站造成的定位误差趋于一致,可忽略其影响。区域内各坐标点总位置误差分布精度较高区域为图中呈近似等腰三角形部分,三角形的顶角指向参考基站方向,三角形的两底角指向对角线方向。若落点位于测量区域边缘对角线附近时,参考基站选取对定位精度影响明显,例如坐标(8 km,8 km),不同参考基站造成X 方向最大误差值相差82.2 m,Y 方向最大误差值相差80.4 m,总位置误差最大相差80.4 m,总位置误差最大相差达到56.2 m,不可忽略其影响。
4 结束语
基于TDOA的多站时差定位体制作为一种有效方法被广泛用于水下无源定位测量,而其定位测量精度则受多种复杂因素影响。为了深入分析高速运动目标深海落点测量精度,本文首先建立了基于TDOA体制定位测量模型,给出了蒙特卡洛数值仿真过程,然后在测量误差源综合影响下,分析了基站布设方式、目标落点位置、参考基站选取对定位测量的影响。通过研究落点定位误差影响规律,对高速运动目标深海落点位置测量基站布设、误差控制、精度评定及相关工程应用具有一定指导意义。
在深海落点测量模型建立过程中,将声速视为常数,但深海中待测目标至基站距离较远时,声波在深海声场沿会聚区声信道传播过程中声线路径弯曲效应明显,若仍将声速视为常量,则对定位精度会产生一定影响,下一步将改进和完善定位模型,建立变声速条件下测量定位模型,同时深入分析声速环境起伏对定位测量精度的影响。
参考文献:
[1]孙大军, 郑翠娥.水声导航、定位技术发展趋势探讨[J].海洋技术学报, 2015, 34(3): 64-67.
[2]刘伯胜, 雷家煜.水声学原理[M].2版.哈尔滨: 哈尔滨工程大学出版社, 2009.
[3]Hinich M J, Sullivan E L.Maximum Likelihood Signal Processing for a Vertical Array[J].Journal of the Acoustical Society of America, 1973, 54(1): 499-503.
[4]Bucker H P.Use of Calculated Sound Fields and Matched Field Detection to Locate Sound Sources in Shallow Water[J].Journal of the Acoustical Society of America, 1976, 59: 368-373.
[5]刘聪锋.无源定位与跟踪[M].西安: 西安电子科技大学出版社, 2011.
[6]Y.Wang and D.Ho.Unified near-field and Far-field Localization For AOA and Hybrid AOA-TDOA Positionings[J].Wireless Commun, 2017, 17(2): 1242-1254.
[7]焦小涛.水下目标被动定位系统关键技术研究[D].西安：西安工业大学, 2013.
[8]张学森, 冯海泓, 黄敏燕，等.两步法提高时延差估计精度的分析和实验验证[J].声学学报, 2013, 38(2): 153-159.
[9]C.Hue, J.-P.Le Cadre, and P.Perez, Sequential Monte Carlo Methods for Multiple Target Tracking and Data Fusion[J].Signal Process, 2002, 50(2): 309-325.
[10]A.O.Hero, S.C.Schwartz.A New Generalized Cross Correlator[J].Signal Processing, 1985, 33(4): 38-45.
[11]Williams S M, Frampton K D, Amundson I, et al.Decentralized Acoustic Source Localization in A distributed Sensor Network[J].Applied Acoustics, 2006, 67(5): 996-1008.
[12]Kozick R J, Sadler B M.Source Localization with Distributed Sensor Arrays and Partial Spatial Coherence[J].IEEE Transactions on Signal Processing, 2004, 52(3): 601-616.
[13]金磊磊, 马艳.任意四元阵的定位盲区讨论及误差影响[J].探测与控制学报, 2015, 37(2): 90-94.
[14]张振中, 夏小华, 马康.舰炮射击海上弹着点声学测量方法及精度仿真[J].指挥控制与仿真, 2014, 36(4): 99-102.
[15]刘德耀, 吴军波, 武翰文.弹丸海上落点声学测量方法研究[J].测试技术学报, 2016, 30(6): 512-516.
[16]陈连, 于国栋, 宋永志，等.声源定位算法的精度分析[J].测绘科学, 2016, 41(5): 21-25.
[17]张旭, 韩旭, 辛健.垂向航行体水下多元阵测量精度仿真算法[J].探测与控制学报, 2017, 39(3): 86-90.
[18]武汉大学测绘学院测量平差学科组.误差理论与测量平差基础[M].武汉: 武汉大学出版社, 2014.
Simulation Analysis of Deep-sea Falling Point Measurement Precision for the High Speed Moving Target
ZHANG Zhi-wei1,2, SUN Ao1, ZHANG Xu1, XIN Jian1
(1.43 Teams in 91550 Troops of PLA, Dalian 116023; 2.Postdoctoral Scientific Research Station in 92493 Troops of PLA, Dalian 116023, China)
Abstract :In order to further research influence rules of the deep-sea falling point measurements accuracy of high-speed moving target which changes with base station layout, falling point position, reference base station, etc., the paper establishes a passive location model based on the Taylor-series expansion method under TDOA system, then takes the site error, delay error and acoustic velocity error as the main error sources variable, plugs the normally distributed random numbers into the measurement model by architectural kind of sampling with the Monte Carlo numerical simulation method, and finally reproduces the accuracy characteristics of time delay differential positioning under the condition of the multi-base stations.The analysis results show that the measurement model can reasonably respond to the changes of base station layout, falling point positon, reference base station selection under the influence of composition error sources, the precision of the deep-sea falling point measurement of the high-speed moving target can be estimated with Monte Carlo numerical simulation method effectively, which can provide references for the design of the multi-base stations in passive location and accuracy evaluation.
Key words :TDOA system; Monte Carlo method; base station layout; falling point position; reference base station; precision analysis
中图分类号: P229.5；E992
文献标志码: A
DOI :10.3969/j.issn.1673-3819.2019.06.017
文章编号： 1673-3819(2019)06-0091-07
收稿日期： 2019-03-11
修回日期： 2019-04-28
*基金项目： 国家自然科学基金(61701504)
作者简介： 张志伟(1987—),男,江苏连云港人,博士研究生,工程师,研究方向为水下测量技术。
孙 翱(1965—),男,高级工程师。
(责任编辑:胡志强)
标签：TDOA体制论文; 蒙特卡洛方法论文; 基站布设论文; 落点位置论文; 参考基站论文; 精度分析论文; 中国人民解放军91550部队43分队论文; 中国人民解放军92493部队博士后科研工作站论文;