一种一维距离像的特征提取方法与流程

1.本公开涉及信号处理技术领域，尤其涉及一种一维距离像的特征提取方法。


背景技术：

2.高超声速飞行器在高空飞行时，在飞行器表面会产生包覆的等离子鞘，会影响其一维距离像在多普勒域的扩展、信号能量的分散、幅度轮廓不规则形变，导致特征提取困难；同时，目标的高超声速飞行会限制样本数量和运算时间，进一步增加了特征提取的难度。
3.现有的距离像特征提取要么不适用于高超声速环境的等离子鞘多扰动的情况，要么数据维数偏大，对数据存储量和计算量要求高，不能适应高超声速目标对计算时间的需求。因此，有必要改善上述相关技术方案中存在的一个或者多个问题。
4.需要说明的是，在上述背景技术部分公开的信息仅用于加强对本公开的背景的理解，因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。


技术实现要素：

5.本公开实施例的目的在于提供一维距离像的特征提取方法，以克服高超声速目标数据量少、计算时间短的应用限制，同时满足对高超声速目标识别的需求。
6.根据本公开实施例提供一种一维距离像的特征提取方法，包括：
7.利用宽带雷达的多个连续脉冲的脉压数据，获得目标脉冲距离像；
8.对所述目标脉冲距离像进行阈值处理，得到阈值处理后的脉冲距离像；
9.从所述阈值处理后的脉冲距离像中提取实一维距离像，并计算多个所述连续脉冲的所述实一维距离像的最大距离单元长度，得到实脉冲距离像；
10.利用所述实脉冲距离像，得到功率脉冲距离像；
11.利用所述功率脉冲距离像，构造元素熵；
12.利用构造出的所述元素熵提取一维距离像的特征参数。
13.本公开的一示例性实施例中，利用宽带雷达的多个连续脉冲的脉压数据，获得目标脉冲距离像的步骤包括：
14.将宽带雷达的多个连续脉冲观测窗内的脉压数据进行拼接，获得所述目标脉冲距离像a
m，n
；所述目标脉冲距离像a
m，n
包括m
×
n的矩阵，其中，m＝1，2，
…
，m，m表示脉冲序号，m表示脉冲总数；n＝1，2，
…
，n，n表示目标脉冲距离像a
m，n
中的距离单元序号，n表示目标脉冲距离像a
m，n
中的距离单元总数。
15.本公开的一示例性实施例中，在对所述目标脉冲距离像进行阈值处理，得到阈值处理后的脉冲距离像的步骤中：所述阈值处理包括硬阈值处理或软阈值处理。
16.本公开的一示例性实施例中，当对所述目标距离像a
m，n
进行硬阈值处理时，得到硬阈值处理后的脉冲距离像h
m，n
包括：
[0017][0018]
其中，t表示硬阈值门限。
[0019]
本公开的一示例性实施例中，从所述阈值处理后的脉冲距离像中提取实一维距离像，并计算多个所述连续脉冲的所述实一维距离像的最大距离单元长度，得到实脉冲距离像的步骤包括：
[0020]
提取第m个脉冲的第一个非零元素所处的距离单元序号b
m，1
：
[0021]bm，1
＝max[find(h
m，n
≠0)]，n＝1，2，
…
，n
ꢀꢀꢀ
(2)
[0022]
其中，find(h
m，n
≠0)表示为读取h
m，n
≠0时的距离单元序号n；
[0023]
提取第m个脉冲的最后一个非零元素所处的距离单元序号b
m，2
：
[0024]bm，2
＝min|find(h
m，n
≠0)]，n＝1，2，
…
，n
ꢀꢀꢀ
(3)
[0025]
其中，h
m，n
表示硬阈值处理后的脉冲距离像；m表示脉冲序号，m表示脉冲总数；n表示观测窗内的距离单元序号，n表示观测窗内的距离单元总数。
[0026]
计算第m个脉冲中实一维距离像的长度km：
[0027]km
＝b
m，2-b
m，1
+1
ꢀꢀ
(4)
[0028]
提取m个脉冲中实一维距离像的最大距离单元长度k
max
：
[0029]kmax
＝max(km)
ꢀꢀ
(5)
[0030]
获得m个脉冲的实一维距离像的所述实脉冲距离像r
m，k
：
[0031]rm，k
＝|a
m，n
|，b
m，1
≤n≤b
m，2
+k
max-1
ꢀꢀꢀ
(6)
[0032]
其中，所述实脉冲距离像r
m，k
包括包括m
×kmax
的矩阵，k表示实脉冲距离像r
m，k
中的距离单元序号，k＝1，2，
…
，k
max
。
[0033]
本公开的一示例性实施例中，利用所述实脉冲距离像，得到功率脉冲距离像的计算公式包括：
[0034][0035]
其中，k表示距离单元序号，k＝1，2，
…
，k
max
；k
max
表示一维距离像的最大距离单元长度；r
m，k
表示实脉冲距离像；p
m，k
表示功率脉冲距离像；m表示脉冲序号。
[0036]
本公开的一示例性实施例中，利用所述功率脉冲距离像，构造元素熵的步骤包括：
[0037]
利用所述功率脉冲距离像p
m，k
定义向量y
g1，g2
；
[0038]
构造所述向量y
g1，g2
的元素熵；
[0039]
其中，1≤g1≤g2≤m；m表示脉冲总数。
[0040]
本公开的一示例性实施例中，利用所述功率脉冲距离像p
m，k
定义向量y
g1，g2
的步骤中，所述向量y
g1，g2
的计算公式包括：
[0041]yg1，g2
(k)＝p
g1，k
p
g2，k
ꢀꢀꢀ
(8)。
[0042]
本公开的一示例性实施例中，在构造所述向量y
g1，g2
的元素熵的步骤中，构造所述元素熵的计算公式包括：
[0043]
[0044]
其中，||
·
||1表示求解1-范数，||
·
||2表示求解2-范数。
[0045]
本公开的一示例性实施例中，利用构造出的所述元素熵提取一维距离像的特征参数的计算公式包括：
[0046][0047]
其中，ln表示求解以10为底的对数。
[0048]
本公开提供的技术方案可以包括以下有益效果：
[0049]
(1)本公开对等离子鞘等强扰动、非线性变化的对象具有很强的适应性，能够在缩短计算时间的同时，保留回波信号的信号强度分布特征，能够满足雷达对高超声速目标的识别要求；
[0050]
(2)本公开对雷达回波数据的脉冲数要求低，在数据量少、存储空间小或计算时间受限的情况下，仍能获得稳定、良好的计算结果。
[0051]
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本公开。
附图说明
[0052]
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本公开的实施例，并与说明书一起用于解释本公开的原理。显而易见的，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0053]
图1示出本公开示例性实施例中一维距离像的提取方法的步骤示意图；
[0054]
图2示出本公开示例性实施例中一维距离像的另一种提取方法的步骤示意图；
[0055]
图3示出本公开示例性实施例中一维距离像的提取方法中经硬阈值处理后的第一脉冲距离像h
m，n
的示意图；
[0056]
图4示出本公开示例性实施例中一维距离像的提取方法中第二脉冲距离像r
m，k
的示意图；
[0057]
图5示出本公开示例性实施例中一维距离像的提取方法中功率脉冲距离像p
m，k
的示意图；
[0058]
图6示出本公开示例性实施例仿真实验中包含等离子鞘的高超声速飞行器元素熵的对数图；
[0059]
图7示出本公开示例性实施例仿真实验中不包含等离子鞘的高超声速飞行器元素熵的对数图；
[0060]
图8示出本公开示例性实施例仿真实验中包含等离子鞘的高超声速飞行器特征参数χe的曲线图；
[0061]
图9示出本公开示例性实施例仿真实验中不包含等离子鞘的高超声速飞行器特征参数χe的曲线图。
具体实施方式
[0062]
现在将参考附图更全面地描述示例实施方式。然而，示例实施方式能够以多种形式实施，且不应被理解为限于在此阐述的范例；相反，提供这些实施方式使得本公开将更加全面和完整，并将示例实施方式的构思全面地传达给本领域的技术人员。所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中。
[0063]
此外，附图仅为本公开的示意性图解，并非一定是按比例绘制。图中相同的附图标记表示相同或类似的部分，因而将省略对它们的重复描述。附图中所示的一些方框图是功能实体，不一定必须与物理或逻辑上独立的实体相对应。可以采用软件形式来实现这些功能实体，或在一个或多个硬件模块或集成电路中实现这些功能实体，或在不同网络和/或处理器装置和/或微控制器装置中实现这些功能实体。
[0064]
本公开实施方式的第一方面提供了一维距离像的提取方法，参照图1和图2，该提取方法包括以下步骤：
[0065]
步骤s101：利用宽带雷达的多个连续脉冲的脉压数据，获得目标脉冲距离像；
[0066]
步骤s102：对该目标脉冲距离像进行阈值处理，得到阈值处理后的脉冲距离像；
[0067]
步骤s103：从阈值处理后的脉冲距离像中提取实一维距离像，并计算多个连续脉冲的实一维距离像的最大距离单元长度，得到实脉冲距离像；
[0068]
步骤s104：利用该实脉冲距离像，得到功率脉冲距离像；
[0069]
步骤s105：利用该功率脉冲距离像，构造元素熵；
[0070]
步骤s106：利用构造出的所述元素熵提取一维距离像的特征参数。
[0071]
下面将对该提取方法的具体步骤进行说明。
[0072]
在步骤s101中，该步骤包括：将宽带雷达的多个连续脉冲观测窗内的脉压数据进行拼接，获得目标脉冲距离像a
m，n
；该目标脉冲距离像a
m，n
为m
×
n的矩阵，其中，m＝1，2，
…
，m，m表示脉冲序号，m表示脉冲总数；n＝1，2，
…
，n，n表示目标脉冲距离像a
m，n
的观测窗内的距离单元序号，n表示目标脉冲距离像a
m，n
的观测窗内的距离单元总数。
[0073]
在步骤s102中，这里对目标脉冲距离像进行阈值处理可以是硬阈值处理，也可以根据实际需要采用软阈值处理，且不会影响后续处理结果。
[0074]
在该步骤中，采用了对目标脉冲距离像a
m，n
进行硬阈值处理，得到硬阈值处理后的脉冲距离像h
m，n
包括：
[0075][0076]
其中，t表示硬阈值门限。
[0077]
这里通过采用最优阈值选择类方法可以得到硬阈值门限t，参照图3所示，对于硬阈值处理后的脉冲距离像h
m，n
，若低于硬阈值门限t，则距离单元取0；若大于等于硬阈值门限t，则距离单元保持原有数值。采用最优阈值选择类方法得到的硬阈值门限t，可直接用于本公开的硬阈值的处理。
[0078]
在步骤s103中，具体包括了以下子步骤：
[0079]
步骤s1031：提取第m个脉冲的第一个非零元素所处的距离单元序号b
m，1
：
[0080]bm，1
＝max|find(h
m，n
≠0)]，n＝1，2，
…
，n
ꢀꢀꢀ
(2)
[0081]
其中，find(h
m，n
≠0)表示为读取h
m，n
≠0时的距离单元序号n；
[0082]
步骤s1032：提取第m个脉冲的最后一个非零元素所处的距离单元序号b
m，2
：
[0083]bm，2
＝min[find(h
m，n
≠0)]，n＝1，2，
…
，n
ꢀꢀꢀ
(3)
[0084]
其中，h
m，n
表示硬阈值处理后的脉冲距离像；m表示脉冲序号，m表示脉冲总数；n表示观测窗内的距离单元序号，n表示观测窗内的距离单元总数
[0085]
步骤s1033：计算第m个脉冲中实一维距离像的长度km：
[0086]km
＝b
m，2-b
m，1
+1
ꢀꢀꢀ
(4)
[0087]
步骤s1034：提取m个脉冲中实一维距离像的最大距离单元长度k
max
：
[0088]kmax
＝max(km)
ꢀꢀꢀ
(5)
[0089]
步骤s1035：获得m个脉冲的实一维距离像的所述实脉冲距离像r
m，k
：
[0090]rm，k
＝|a
m，n
|，b
m，1
≤n≤b
m，2
+k
max-1
ꢀꢀꢀ
(6)
[0091]
其中，实脉冲距离像r
m，k
包括包括m
×kmax
的矩阵，k表示实脉冲距离像r
m，k
中的距离单元序号，k＝1，2，
…
，k
max
。
[0092]
参照图4所示，这里包括了等离子鞘的一维距离像，它的信号能量强点与脉冲维基本对齐，但这并不能形成清晰稳定的直线状航迹，这显示的是一个等离子鞘影响下一维距离像的特征。需要注意的是，这里从硬阈值处理后的脉冲距离像h
m，n
中提取，并计算得到的实脉冲距离像r
m，k
包含背景噪声。
[0093]
在步骤s104中，利用实脉冲距离像r
m，k
，计算得到功率脉冲距离像p
m，k
的计算公式为：
[0094][0095]
其中，k表示距离单元序号，k＝1，2，
…
，k
max
；k
max
表示一维距离像的最大距离单元长度；r
m，k
表示实脉冲距离像；p
m，k
表示功率脉冲距离像；m表示脉冲序号。
[0096]
这里需要说明的是，由于熵是用来衡量信号的稀疏性和相关性。如果距离单元的数值为零，则相关性远比目标信号的分量高，就会影响到特征提取。
[0097]
参照图5所示，图5中有两条稳定的航迹，右侧航迹略微倾斜，信号能量比图4更为集中，从而形成了更为清晰稳定的航迹。
[0098]
在步骤s105中，包含以下子步骤：
[0099]
步骤s1051：利用功率脉冲距离像p
m，k
定义向量y
g1，g2
；这里向量y
g1，g2
的计算公式包括：
[0100]yg1，g2
(k)＝p
g1，k
p
g2，k
ꢀꢀꢀ
(8)
[0101]
其中，1≤g1≤g2≤m；m表示脉冲总数。
[0102]
步骤s1052：构造向量y
g1，g2
的元素熵。元素熵的计算公式包括：
[0103][0104]
其中，||
·
||1表示求解1-范数，||
·
||2表示求解2-范数。
[0105]
在步骤s106中，利用构造出的元素熵提取一维距离像的特征参数的计算公式包
括：
[0106][0107]
其中，ln表示求解以10为底的对数。
[0108]
为了更好的说明本公开的技术效果，进行了如下的仿真实验：
[0109]
仿真条件：
[0110]
采用飞行器在飞行高度为30km，飞行速度为25ma的等离子鞘数据，数据来源于典型钝锥的三维高超声速cfd流场仿真数据。设定的雷达类型为x波段宽带雷达，信号类型为线性调频信号，载波载频f＝10ghz，带宽b＝1ghz，脉宽τ＝200μs，采样频率fs＝5ghz。
[0111]
仿真内容：
[0112]
仿真内容一、分别模拟仿真包含等离子鞘和不包含等离子鞘的高超声速飞行器典型钝锥的元素熵的对数图。
[0113]
这里选取起始俯仰角为60
°
的60个连续脉冲的数据作为实验数据。图6为包含等离子鞘的高超声速飞行器元素熵对数图，图7为不包含等离子鞘的高超声速飞行器元素熵对数图。在图6和图7中，x轴表示一维距离像的相对距离范围，y轴表示元素熵序号，z轴表示熵值。从图6和图7中可以看出，包含等离子鞘的一维距离像占据更宽的相对距离范围，信号能量更为分散，元素熵曲线的极小值点更多。
[0114]
仿真内容二、分别模拟仿真包含等离子鞘和不包含等离子鞘的高超声速飞行器典型钝锥的特征参数χe曲线图。
[0115]
这里选取起始俯仰角为60
°
的60个连续脉冲的数据作为实验数据。图8示出了，处理脉冲数分别为m＝5、m＝10、m＝15、m＝20，仿真实验中包含等离子鞘的高超声速飞行器特征参数χe曲线图。图9示出了，处理脉冲数分别为m＝5、m＝10、m＝15、m＝20，仿真实验中不包含等离子鞘的高超声速飞行器特征参数χe曲线。从图8和图9中可以看出：
[0116]
首先，脉冲数越多曲线越平滑，但脉冲数的变化对特征参数χe曲线的而影响细微，在计算时间受限的情况下计算结果更稳定；
[0117]
其次，包含等离子鞘的一维距离像占据更宽的距离单元范围，信号能量更为分散、元素熵曲线的极小值点更多。
[0118]
需要说明的是，尽管在附图中以特定顺序描述了本公开中方法的各个步骤，但是，这并非要求或者暗示必须按照该特定顺序来执行这些步骤，或是必须执行全部所示的步骤才能实现期望的结果。附加的或备选的，可以省略某些步骤执行等。另外，也易于理解的是，这些步骤可以是例如在多个模块/进程/线程中同步或异步执行。
[0119]
关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
[0120]
应当注意，尽管在上文详细描述中提及了用于动作执行的系统的若干单元，但是这种划分并非强制性的。实际上，根据本公开的实施方式，上文描述的两个或更多单元的特征和功能可以在一个单元中具体化。反之，上文描述的一个单元的特征和功能可以进一步划分为由多个单元来具体化。可以根据实际的需要选择其中的部分或者全部单元来实现本公开方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
[0121]
通过以上的实施方式的描述，本领域的技术人员易于理解，这里描述的示例实施
方式可以通过软件实现，也可以通过软件结合必要的硬件的方式来实现。因此，根据本公开实施方式的技术方案可以以软件产品的形式体现出来，该软件产品可以存储在一个非易失性存储介质(可以是cd-rom，u盘，移动硬盘等)中或网络上，包括若干指令以使得一台计算设备(可以是个人计算机、服务器、或者网络设备等)执行根据本公开实施方式的上述沙盘联动控制方法。
[0122]
本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本公开的其它实施方案。本技术旨在涵盖本公开的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本公开的一般性原理并包括本公开未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本公开的真正范围和精神由所附的权利要求指出。