一种基于分式逼近的方向角航迹预测方法与流程

本发明涉及航迹预测技术领域，尤其涉及一种基于分式逼近的方向角航迹预测方法。

背景技术：

随着电磁环境不断复杂化，无源定位的要求不断提高，国内外学者对测向交叉定位的研究也在不断深入。

在无源定位系统中，测向交叉定位是无源定位中运用最广泛、研究时间最长的一种定位方法。这种定位方法需要多个基站协同工作来完成，由此，当信息源观测数据出现野值或者测量不准确时，多个基站进行交叉定位时容易出现虚假目标点或者导致定位不准确。如何剔除测量野值，并且选择合适的信息源方向角航迹用于交叉定位，直接影响定位解算的精度。

因此，基于航迹测量值与航迹预测值进行合理的方向角航迹筛选是交叉定位方法的一个重要环节，目前关于如何进行方向角航迹预测方法的研究很少。

技术实现要素：

为解决上述问题，本申请实施例提供了一种基于分式逼近的方向角航迹预测方法，所述方法包括：

对纯方向角航迹进行初始化处理，得到处理后航迹；

基于预设时间段内的处理后航迹，确定分式拟合函数；

根据所述分式拟合函数预测指定时刻的方向角，其中指定时刻非所述预设时间段内的时刻，且指定时刻晚于所述预设时间段。

可选地，所述对纯方向角航迹进行初始化处理，包括：将所述方向角由角度制转化为弧度制。

可选地，所述将所述方向角由角度制转化为弧度制，包括：

通过公式将所述方向角由角度制d°转化为弧度制d。

可选地，所述预设时间段为时间为tj-n+1,tj-n+2,…,tj-1的时间段；

其中，j为指定时刻，n为预设的时刻总数量，tj-n+1为处理后航迹中帧号为j-n+1的帧对应的时刻，tj-n+2为处理后航迹中帧号为j-n+2的帧对应的时刻，tj-1为处理后航迹中帧号为j-1的帧对应的时刻，所述处理后航迹由多张视频帧按时间顺序组合而成。

可选地，所述基于预设时间段内的处理后航迹，确定分式拟合函数，包括：

基于奇异值分解的近似方法确定分式拟合函数参数a，b，c，d的值；

基于所述a，b，c，d的值确定分式拟合函数

其中，t为所述预设时间段内的任一时刻。

可选地，所述基于奇异值分解的近似方法确定分式拟合函数参数a，b，c，d的值，包括：

计算矩阵a＝[a1a2a3a4]的奇异值分解svd(a)＝[usv]；

根据分解结果确定拟合函数参数a，b，c，d的值；

其中，a1＝[11…1]^t，a2＝[tj-n+1tj-n+2…tj-1]^t，a4＝[-dj-n+1tj-n+1-dj-n+2tj-n+2…-dj-1tj-1]^t；

dj-n+1为处理后航迹中帧号为j-n+1的帧对应的方向角弧度值，dj-n+2为处理后航迹中帧号为j-n+2的帧对应的方向角弧度值，dj-1为处理后航迹中帧号为j-1的帧对应的方向角弧度值，u为奇异值分解得到的酉矩阵，s为奇异值分解得到的对角矩阵，v为奇异值分解得到的酉矩阵的共轭转置。

可选地，所述根据分解结果确定拟合函数参数a，b，c，d的值，包括：

[a,b,c,d]＝v*x；

其中，x＝b(i)/s(i,i)；b(i)为矩阵b中的元素，b＝u^t*y，y＝[dj-n+1dj-n+2…dj-1]^t，i为元素标识。

可选地，所述根据所述分式拟合函数预测指定时刻方向角，包括：

指定时刻j的方向角

可选地，得到指定时刻方向角后，通过将其由弧制转化为角度制

可选地，所述方向角包括：方位角和俯仰角。

有益效果如下：

本申请提供的方法对于纯方向角航迹进行初始化处理后，通过基于预设时间段内处理后航迹计算分式拟合函数，然后采用拟合函数计算指定时刻的状态预测值，可用于纯方向角航迹的机动检测以及航迹质量评估，也可用于剔除航迹野值。

附图说明

下面将参照附图描述本申请的具体实施例，其中：

图1示出了本申请一实施例提供的一种基于分式逼近的方向角航迹预测方法的流程示意图；

图2示出了本申请一实施例提供的另一种基于分式逼近的方向角航迹预测方法的流程示意图。

具体实施方式

在无源定位系统中，测向交叉定位是无源定位中运用最广泛、研究时间最长的一种定位方法。这种定位方法需要多个基站协同工作来完成，由此，当信息源观测数据出现野值或者测量不准确时，多个基站进行交叉定位时容易出现虚假目标点或者导致定位不准确。如何剔除测量野值，并且选择合适的信息源方向角航迹用于交叉定位，直接影响定位解算的精度。

基于此，本申请提供一种方法，对于纯方向角航迹进行初始化处理后，通过基于预设时间段内处理后航迹计算分式拟合函数，然后采用拟合函数计算指定时刻的状态预测值，可用于纯方向角航迹的机动检测以及航迹质量评估，也可用于剔除航迹野值。

参见图1，本实施例提供的基于分式逼近的方向角航迹预测方法实现流程如下：

101，对纯方向角航迹进行初始化处理，得到处理后航迹。

本步骤对纯方向角航迹进行初始化处理，包括：将方向角由角度制转化为弧度制。

具体的，通过公式将方向角由角度制d°转化为弧度制d。

在具体实施时，方向角包括：方位角和俯仰角。

因此，本步骤分别将方位角由角度制转化为弧度制，将俯仰角由角度制转化为弧度制。

例如，通过公式将方位角由角度制a°转化为弧度制a。通过公式将俯仰角由角度制e°转化为弧度制e。

102，基于预设时间段内的处理后航迹，确定分式拟合函数。

其中，预设时间段为时间为tj-n+1,tj-n+2,…,tj-1的时间段。

其中，j为指定时刻，指定时刻非预设时间段内的时刻，且指定时刻晚于预设时间段。n为预设的预设时间段时刻总数量，tj-n+1为处理后航迹中帧号为j-n+1的帧对应的时刻，tj-n+2为处理后航迹中帧号为j-n+2的帧对应的时刻，tj-1为处理后航迹中帧号为j-1的帧对应的时刻，处理后航迹由多张视频帧按时间顺序组合而成。

本步骤的实现过程为：基于奇异值分解的近似方法确定分式拟合函数参数a，b，c，d的值。基于a，b，c，d的值确定分式拟合函数

其中，t为预设时间段内的任一时刻。

由于方向角包括：方位角和俯仰角，因此本步骤中确定的分式拟合函数包括方位角的分式拟合函数和俯仰角的分式拟合函数

对于基于奇异值分解的近似方法确定分式拟合函数参数a，b，c，d的值的过程如下：

1、计算矩阵a＝[a1a2a3a4]的奇异值分解svd(a)＝[usv]。

其中，a1＝[11…1]^t，a2＝[tj-n+1tj-n+2…tj-1]^t，a4＝[-dj-n+1tj-n+1-dj-n+2tj-n+2…-dj-1tj-1]^t。

dj-n+1为处理后航迹中帧号为j-n+1的帧对应的方向角弧度值，dj-n+2为处理后航迹中帧号为j-n+2的帧对应的方向角弧度值，dj-1为处理后航迹中帧号为j-1的帧对应的方向角弧度值，u为奇异值分解得到的酉矩阵，s为奇异值分解得到的对角矩阵，v为奇异值分解得到的酉矩阵的共轭转置。

由于方向角包括：方位角和俯仰角，因此本步骤中计算方位角的矩阵aa＝[aa1aa2aa3aa4]的奇异值分解svd(aa)＝[uasava]和俯仰角的矩阵ae＝[ae1ae2ae3ae4]的奇异值分解svd(ae)＝[ueseve]。

其中，aa1＝[11…1]^t，aa2＝[tj-n+1tj-n+2…tj-1]^t，aa4＝[-aj-n+1tj-n+1-aj-n+2tj-n+2…-aj-1tj-1]^t。

aj-n+1为处理后航迹中帧号为j-n+1的帧对应的方位角弧度值，aj-n+2为处理后航迹中帧号为j-n+2的帧对应的方位角弧度值，aj-1为处理后航迹中帧号为j-1的帧对应的方位角弧度值，ua为方位角奇异值分解得到的酉矩阵，sa为方位角奇异值分解得到的对角矩阵，va为方位角奇异值分解得到的酉矩阵的共轭转置。

ae1＝[11…1]^t，ae2＝[tj-n+1tj-n+2…tj-1]^t，ae4＝[-ej-n+1tj-n+1-ej-n+2tj-n+2…-ej-1tj-1]^t。

ej-n+1为处理后航迹中帧号为j-n+1的帧对应的俯仰角弧度值，ej-n+2为处理后航迹中帧号为j-n+2的帧对应的俯仰角弧度值，ej-1为处理后航迹中帧号为j-1的帧对应的俯仰角弧度值，ue为俯仰角奇异值分解得到的酉矩阵，se为俯仰角奇异值分解得到的对角矩阵，ve为俯仰角奇异值分解得到的酉矩阵的共轭转置。

本步骤中的奇异值分解方案采用现有方案，此处不再赘述。

2、根据分解结果确定拟合函数参数a，b，c，d的值。

具体的，[a,b,c,d]＝v*x。

其中，x＝b(i)/s(i,i)。b(i)为矩阵b中的元素，b＝u^t*y，y＝[dj-n+1dj-n+2…dj-1]^t，i为元素标识。

由于方向角包括：方位角和俯仰角，因此本步骤中确定方位角拟合函数参数aa，ba，ca，da的值[aa,ba,ca,da]＝va*xa，确定俯仰角拟合函数参数ae，be，ce，de的值[ae,be,ce,de]＝ve*xe。

其中，xa＝ba(ia)/sa(ia,ia)。ba(ia)为矩阵ba中的元素，ya＝[aj-n+1aj-n+2…aj-1]^t，ia为元素标识。xe＝be(ie)/se(ie,ie)。be(ie)为矩阵be中的元素，ie为元素标识。

103，根据分式拟合函数预测指定时刻的方向角。

其中，指定时刻可以为当前时刻，也可以为下一时刻，还可以为用户指定的非预设时间段内，且于预设时间段的任一时刻。

指定时刻j的方向角

另外，得到指定时刻方向角后，通过将其由弧制转化为角度制

由于方向角包括：方位角和俯仰角，因此本步骤中根据分式拟合函数预测指定时刻的方位角和俯仰角。

具体的，指定时刻j的方位角指定时刻j的俯仰角

相应的，得到指定时刻方位角后，通过将其由弧制转化为角度制得到指定时刻俯仰角后，通过将其由弧制转化为角度制

本实施例提供的方法，对于纯方向角(包括方位角和俯仰角)航迹，采样一段时间窗的历史测量值通过分式函数拟合的方法进行方向角航迹分段分式逼近，然后基于拟合的分式函数计算预测时刻的方向角(方位角与俯仰角)。可用于纯方向角航迹的机动检测以及航迹质量评估，也可用于剔除航迹野值。

下面参照图2，对本实施例提供的方法再次进行说明。

第一步，航迹拟合的初始化处理

将测量方向角由角度制转化为弧度制，记方向角航迹的方位角、俯仰角的角度值分别为a°，e°，进行航迹的分式函数拟合前首先将相应的角度值转化为弧度制：

第二步，基于前n帧方向角航迹测量值，计算该时间段的分式拟合函数

其中，n为预设的预设时间段时刻总数量，由于处理后航迹由多张视频帧按时间顺序组合而成，每一时刻对应处理后航迹中的一帧，因此，本步骤实际是基于预设时间段(tj-n+1,tj-n+2,…,tj-1的时间段)内的处理后航迹，确定分式拟合函数。

具体的，设前n个时刻的方向角测量时刻分别为tj-n+1,tj-n+2,…,tj-1，方位角的测量值为aj-n+1,aj-n+2,…,aj-1，俯仰角的测量值为ej-n+1,ej-n+2,…,ej-1，基于这n个时间窗的测量值，分别计算该时间段方位角航迹与俯仰角航迹的拟合函数，具体步骤如下：

a)设方位角的分式拟合函数具有如下形式

b)然后采用近似方法计算方位角航迹分式拟合函数的参数，记coea＝[aa,ba,ca,da]：

计算中间参数：

aa1＝[11…1]^t

aa2＝[tj-n+1tj-n+2…tj-1]^t

aa4＝[-aj-n+1tj-n+1-aj-n+2tj-n+2…-aj-1tj-1]^t

记aa＝[aa1aa2aa3aa4]，ya＝[aj-n+1aj-n+2…aj-1]^t，下面基于奇异值分解的近似求解拟合参数。

1)计算矩阵aa的奇异值分解：

svd(aa)＝[uasava]

2)然后求解拟合系数：

xa＝ba(ia)/sa(ia,ia)

coea＝[aa,ba,ca,da]＝va*xa

c)基于前n帧俯仰角测量值，计算俯仰角航迹的分式拟合函数

对于该时间窗的俯仰角航迹进行拟合，设俯仰角航迹的分式拟合函数为

拟合函数的参数计算步骤同上，不再赘述。设解算的俯仰角航迹参数为coee＝[ae,be,ce,de]。

第三步，基于拟合函数，计算指定时刻j的预测方位角与俯仰角。

由于coea＝[aa,ba,ca,da]，因此，coee＝[ae,be,ce,de]，因此，

第四步，将预测方位角由弧度制转化为角度制

对于指定时刻(如下一时刻)，仍然基于前n个时间窗的测量值，计算方向角航迹的分式拟合函数，并基于拟合函数计算该时刻的方位角与俯仰角预测值。

例如对指定时刻的方向角航迹预测时，1.首先进行航迹拟合的初始化处理；2.然后基于前n帧方向角航迹测量值，分别计算该时间段的方位角与俯仰角分式拟合函数；3.采用拟合函数，计算下一预测时刻的预测方位角与俯仰角；4.将预测方位角由弧度制转化为角度制，并进行一下周期的计算。

本实施例提供的方法，通过基于一定时间窗的历史测量值计算分段分式拟合函数，然后采用拟合函数计算指定时刻的状态预测值。可用于纯方向角航迹的机动检测以及航迹质量评估，也可用于剔除航迹野值。

有益效果：

本申请提供的方法对于纯方向角航迹进行初始化处理后，通过基于预设时间段内处理后航迹计算分式拟合函数，然后采用拟合函数计算指定时刻的状态预测值，可用于纯方向角航迹的机动检测以及航迹质量评估，也可用于剔除航迹野值。