搜索
，姚放吾1，张正宏21南京邮电大学计算机学院，南京（210003）2上海宏皓科技有限公司，上海（200136）E-mail：luoweicawy@yahoo.com.cn摘要：船用导航雷达属于边扫描边跟踪雷达系统。它要求对多目标录取的同时，进行目标跟踪，并为各个目标建立相应的航迹。α-β-γ滤波理论是一种基于匀加速运动目标的线性滤波方法。本文提出了利用ARM处理器实现该理论算法的船用导航雷达跟踪设计方法。并对跟踪设计中遇到的目标录取、机动检测和航迹相关进行研究分析，提出相应的解决方法。最后，对跟踪性能进行了实验仿真。关键词：航海雷达；α-β-γ滤波；自动跟踪；机动检测1.引言船用导航雷达系统可以同时对多目标（100以上）进行实时录取、跟踪、自动计算出目标的动态信息，以避免碰撞。这要求跟踪系统能对目标进行预测和滤波。由于船舶机动比较缓慢，雷达天线扫描周期(即：目标一次采样时间，一般在2.5秒以内)比较短，所以可以近似看作是线性运动。目前常应用于雷达的线性滤波的理论有：维纳滤波、卡尔曼滤波，α-β滤波、α-β-γ滤波。其中维纳滤波是所有线性滤波中最佳的，但需要丰富的验前信息。卡尔曼滤波是采用递推式滤波方法，只需要前一次历史数据信息，适合于计算机处理，但它要求不断的更新计算新息。α-β滤波、α-β-γ滤波都是卡尔曼滤波的简化模型，避免了对新息的大量计算，更适用于多目标的边扫描边跟踪系统。其中前者应用于匀速运动模型，后者应用于匀加速运动模型，并能够对加速度进行预测和滤波，有利于机动检测。本文着重于将α-β-γ滤波理论应用于ARM的具体设计，并对跟踪技术中的目标录取、机动检测和航迹相关具体设计进行分析，提出适合于ARM处理器的解决方法。2.跟踪工作模型雷达接收机接收到回波信号后，经杂波处理及数字量化送录取模块进行目标录取和航迹相关处理，再把录取目标的距离、方位信息送录取跟踪模块，对目标进行跟踪和机动检测，最后计算输出目标的参数信息。如下图1所示图1跟踪工作原理其中：1.数字化回波信号模块根据跟踪精度，把原始回波信号数据经杂波抑制处理后，进行数字量化处理。2.目标录取模块负责录取波门中的点集，进而找出波门中的所有可能目标。3.航迹相关模块对波门中的所有可疑目标进行航迹相关，从中选一个最合适的作为真实跟踪目标。4.目标跟踪模块依据目标的测量数据，对跟踪波门进行α-β-γ滤波和机动检测，预测波门的位置,确定波门的大小，以保证对目标的正确跟踪，减少目标跟踪丢失率。5.目标参数输出模块计算目标属性参数，如：速度、航向、TCPA、CPA等跟踪数据。3.基于ARM硬件设计设计时考虑设计成本和性能，整体设计硬件接口（如图2）。图2硬件接口为了充分利用硬件资源，跟踪器和雷达人机界面由单块ARM实现。设计采用无操作系统方案，有利于人机界面和跟踪模块的有效工作。并且FPGA和ARM之间的ARPA回波数据和图形数据均采用DMA（行中断触发）方式进行数据交互。试验表明，可以满足60M的图形界面刷新率和ARPA数据的实时跟踪。由于ARM处理速度有限，这就要求跟踪算法简单有效，使用α-β-γ滤波作为跟踪器，并在跟踪过程中做机动检测和波门位置大小调整，有利于提高跟踪速度和精确度。跟踪设计过程为：1．ARM负责把手动，自动录取波门及实时变化的目标跟踪波门送FPGA。2．FPGA负责对雷达回波信号中波门内的目标进行录取，并把目标的方位和距离送ARM进行α-β-γ滤波和目标跟踪。4.基于ARM软件设计4.1ARPA目标的录取常用的跟踪方式[1]分为：边缘跟踪，形心跟踪，重心跟踪，峰值跟踪等。其中，边缘跟踪是以目标边缘作为跟踪参考点的跟踪方法，在导航避碰雷达中，一般取目标的最前缘作为边缘跟踪的参考点，并且多用于近目标跟踪。形心跟踪是重心跟踪方法的简化，它把目标看作是密度均匀的物体，并以此方式下求出的重心作为跟踪参考点的跟踪方法，当目标姿态变化时,形心点浮动范围较小,所以用形心跟踪比较平稳，而且抗杂波干扰的能力强。形心跟踪常应用于多目标跟踪，由于它把目标看作是密度均匀的物体，所以对目标点集的描述上只用一位来表示，即：目标点集的有无。这样可以大大提高目标的录取速度。峰值跟踪是以目标图像上最亮点或最暗点作为跟踪参考点的一种跟踪方法。也可以把目标图像划分为均匀的N个区域，然后对N个区域分别进行亮度统计，把统计值最大或最小的那个区域作为峰值所在的区域，最后再在这个区域中取一个点作为峰值。峰值跟踪法能跟踪任意大小的目标图像，但它更适合跟踪小目标，尤其是跟踪远处目标或运动较快的目标。在ARM设计中，动态地选择不同的跟踪点有利于提高跟踪速度，优化跟踪性能。近距离（如3海里以内）慢速目标，采用前缘点作为跟踪点。跟踪较远距离（如12海里以内）目标，采用形心作为跟踪点。跟踪远距离（如12海里以外）目标，采用峰值作为跟踪点。4.2航迹相关当几个录取波门交叠或相隔较近时，新点集会同时与多个波门相关，容易造成目标航迹混淆，导致目标跟踪丢失。解决方法可以先对交叠或间隔较近的多个录取波门采取联合录取，再把录取后的目标分配给相应录取波门。分配原则可以充分依赖历史航迹信息，但在船用导航雷达系统中，由于同时跟踪目标多，实时性较强，必须降低数据处理量。一种有效的分配原则是：最近优先（距波门中心），质量匹配其次。4.3目标跟踪4.3.1α-β-γ滤波理论α-β-γ滤波[2]是一种常增益卡尔曼滤波，目标的状态分量(x，v，a)为：位置，速度和加速度。利用α-β-γ滤波可以预测目标下次可能出现的位置。状态预测方程为：X(n+1|n)=F(n)X(n|n)(1)状态更新方程：X(n+1|n+1)=X(n+1|n)+G(n+1)[z(n+1)–H(n+1)X(n+1|n)](2)常增益：G(n+1)=(α,β/T,γ/2T)′(3)状态转移矩阵:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡= １ ／ ２100T02T1F(n)T(4)量测矩阵：H(n+1)=(1，0，0)(5)各分量常增益(α从n=1开始取值，β从n=2开始取值，γ从n=3开始取值)：α=3(32n–3n+2)/(n(n+1)(n+2))(6)β=8(2n-1)/(n(n+1)(n+2))(7)γ=60/(n(n+1)(n+2))(8)其中：n为采样序数1，2，3…z(n+1)为第n+1次测量数据。T为测量数据的采样周期，在船用雷达系统中表现为：天线旋转周期。α,β,γ分别为状态位置、速度和加速度分量的常增益。4.3.2录取波门位置预测、大小调节为了实现同时对多目标进行实时跟踪，必须尽可能地减少目标的搜索范围，提高目标的录取速度。常用的方法是：确定目标录取波门后，在波门内进行目标录取。它要求目标以相对较大的概率落在波门内，否则，容易造成目标录取失败。采用α-β-γ滤波技术可以预测目标录取波门的中心位置，通过合理地动态设置波门大小，可以大大提高目标落在波门内的概率，减少目标搜索失败率。4.3.2.1波门位置预测和波门外推在对回波扫描线点集录取过程中，由于测量噪声、机内噪声和杂波等数据的存在，导致获取的测量数据并不一定是目标的真实数据，和实际数据会有一定的误差。所以在对目标录取波门预测时，要对目标位置，速度及加速度进行平滑滤波。根据3.3.1中的α-β-γ滤波方程可以得出下列计算公式[2]（式中参数含义同3.3.1中的α-β-γ滤波方程中的参数,α，β，γ计算方法同(6),(7),(8)）：x(n|n)=x(n|n-1)+α[z(n)–x(n|n-1)](9)v(n|n)=v(n|n-1)+β[z(n)–x(n|n-1)]/T(10)a(n|n)=a(n|n-1)+γ[z(n)–x(n|n-1)]/2T(11)x(n+1|n)=x(n|n)+v(n|n)T+a(n|n)2T/2(12)v(n+1|n)=v(n|n)+a(n|n)*T(13)a(n+1|n)=a(n|n)(14)运用公式(6)、(7)、(8)可以分别计算距离、速度和加速度的平滑值，公式(9)根据平滑值可以估算下一次波门中心的位置。利用以上公式时，首先录取获得前三次的量测数据z(1)、z(2)和z(3)，再利用z(1)、z(2)和z(3)初始化x(3|2)=z(3)，v(3|2)=(z(3)-z(2))/T，a(3|2)=((z(3)-z(2))/T-(z(2)-z(1))/T)/T，以后逐个递归算出波门中心位置的估值。在船用导航雷达系统中，雷达扫描的范围为海平面，从雷达回波扫描线录取得到的目标测量数据（即：波门中心位置的测量数据）为（γ,θ）平面极坐标（即：于本船的距离和方位），而目标在距离和方位上不满足匀加速运动模型，也就不能在极坐标体系中运用α-β-γ滤波。解决方法是：先把测量的极坐标（γ,θ）数据映射成平面直角坐标（x，y）数据，然后在直角坐标体系中，对x和y方向分别采用α-β-γ滤波估算目标在两个方向上相对本船的位置(即：波门中心位置估值)，速度等信息。由于扫描线是实时处理的，ARM单位时间内处理的数据量也就非常大。故，在进行坐标转换时，应该采用查表方式，以提高雷达跟踪性能。录取测量数据时，船舶打舵、杂波等因素，会导致目标测量数据录取失败，从而无法进行α-β-γ滤波。这时需要以上次的x(n-1|n-1)、v(n-1|n-1)和a(n-1|n-1)值对录取波门进行外推。即：x(n|n)=x(n-1|n-1)+v(n-1|n-1)T+a(n-1|n-1)2T/2v(n|n)=v(n-1|n-1)a(n|n)=a(n-1|n-1)4.3.2.2波门大小调节通过α-β-γ滤波，可以估算波门中心在下一次采样周期可能出现的位置，但估值和目标实际出现的位置存在一定的误差，尤其是在目标做机动时，误差更为明显。所以，每次滤波结束后，应该合理地设置波门大小，波门越大，目标点集录取成功率越高，否则，越低。可以通过下式计算波门的最大尺寸：W=MVT(15)其中：W,V，T分别为：波门的最大尺寸，目标最大速度和目标采样周期(天线旋转周期)M：波门放大系数，与噪声（机动作为噪声处理），杂波有关，取值大小根据数据统计得出。录取第二次量测数据时，取第一次量测数据作为第二次录取波门的中心，波门大小则设置为W。船舶的航速一般不大于50节，目标采样周期设为3秒，m取10，由式(15)计算出W大概是1500米。故：船用导航雷达录取波门最大尺寸可设置为1500米*1500米。随着采样序列的增加，α，β，γ取值越小，目标跟踪逐渐趋于稳定，波门大小逐步缩小。当采样序列达到一定次数K时，它们几乎接近于0，波门中心位置估值基本上与实际目标位置吻合，此时波门大小设置为目标的L(L1)倍。其中剩下的(1–1/L)波门大小空间保证目标机动时，能够对目标录取成功，以便对目标做机动检测。L可以通过实验数据得出，由于船舶机动比较缓慢，L一般可以取值为1倍到2倍之间。根据式（6）计算得出，K取值一般要大于25。所以，从开始目标跟踪到跟踪次数达到K次过程中，可以分阶段地调节波门大小。4.4机动检测α-β-γ滤波只实用于均加速模型，而船舶运动不可避免的会出现机动，所以，在采用α-β-γ滤波过程中，需要对目标运动做机动检测，当目标被检测做机动时，为了减少目标录取失败率，需要通过增加录取波门尺寸，扩大目标点集的搜索区域，通过增大α，β，γ增益系数，使波门中心预测值靠近测量值。一种机动检测的方法是计算目标新息变化[3]。计算公式如下：v(n)=|z(n)–x(n|n-1)|(16)其中：v(n)为：第n次目标新息。z(n)为：第n次测量数据。x(n|n-1)为：第n-1次对第n次估值。利