本发明涉及微波雷达技术领域，尤其涉及一种改进的基于de模型的频带拼接方法、装置、设备及介质。
背景技术：
雷达实测获取的回波数据带宽通常较窄，不能满足对目标进行高分辨特征提取的需求。为此，需要对多段不相邻的频域数据进行频带拼接。目前常用的频带拼接方法主要包括基于gtd模型的频带拼接方法和基于de模型的频带拼接方法。其中，gtd模型参数提取较为复杂，且对回波数据的拟合能力不及de模型。而de模型因为忽视了散射中心随频率的变化规律，因此对于大带宽的数据难以达到良好的拼接效果。
技术实现要素：
本发明要解决的技术问题在于，针对现有基于de模型的频带拼接方法对于大带宽的数据的拼接效果不佳的缺陷，提供一种改进的基于de模型的频带拼接方法、装置、设备及介质。
为了解决上述技术问题，本发明第一方面，提供了一种改进的基于de模型的频带拼接方法，该方法包括：
采用旋转不变技术估计信号参数方法联合各频段数据估计散射中心位置；
使用第一段频带的数据估计散射中心幅度系数；
计算第二段频带至第p段频带的幅度校正因子，进行p-1次多项式拟合，将缺失频段中心频率代入获得缺失频带的幅度校正因子；
根据估计的散射中心位置、缺失频带的幅度校正因子和散射中心幅度系数生成缺失频带数据，完成频带拼接。
在根据本发明的改进的基于de模型的频带拼接方法中，优选地，所述根据估计的散射中心位置、缺失频带的幅度校正因子和散射中心幅度系数生成缺失频带数据，包括：
通过以下公式生成缺失频带数据：
其中，yq(n)为待估计的第q段缺失频带中第n个频点对应的雷达回波数据，fqn对应着待估计的第q段缺失频带的第n个频点，1≤q≤p-1；
在根据本发明的改进的基于de模型的频带拼接方法中，优选地，所述采用旋转不变技术估计信号参数方法联合各频段数据估计散射中心位置，包括：
联合各频段数据生成信号子空间es：
对ψ进行特征值分解，求得散射中心位置rm的估计：
其中，e1表示es去掉最后一行后形成的矩阵，e2表示es去掉第一行后形成的矩阵，λm表示ψ的第m个特征值，angle(·)表示求相角运算，c代表光速，δf为步进频率间隔。
在根据本发明的改进的基于de模型的频带拼接方法中，优选地，所述使用第一段频带的数据估计散射中心幅度系数，包括：
由最小二乘法得到对第一段频带的幅度系数向量amp1的估计为：
y1＝[y1(1),y1(2),…,y1(n)]，y1(n)为第一段频带的第n个频点对应的回波数据，1≤n≤n，n为第一段频带的回波总数，a＝[a(r1)t,a(r2)t,…a(rm)t]，a(rm)为第一段频带的第m个散射中心对应的导向矢量；
根据amp1＝[a1,a2,…,am]估计第1段频带的数据中第m个散射中心的散射复幅度am，1≤m≤m。
在根据本发明的改进的基于de模型的频带拼接方法中，优选地，所述计算第二段频带至第p段频带的幅度校正因子，包括：
用最小二乘法求第二段频带至第p段频带的幅度系数向量ampi的估计，i＝2,…,p；
其中，a＝[a(r1)t,a(r2)t,…a(rm)t]，a(rm)为第一段频带的第m个散射中心对应的导向矢量，yi＝[yi(1),yi(2),…,yi(n)]，yi(n)为第i段频带的第n个频点对应的回波数据；
由此上式获得对幅度校正因子的估计：
kim为第i段频带的第m个散射中心的幅度校正因子。
在根据本发明的改进的基于de模型的频带拼接方法中，优选地，所述进行p-1次多项式拟合，将缺失频段中心频率代入获得缺失频带的幅度校正因子，包括：
对第m个散射中心，估计该散射中心在缺失频带的幅度校正因子，p段数据的幅度校正因子组成向量k＝[k1m,k2m,…kpm],对k进行p-1次多项式拟合：
其中，fci为第i段频带的中心频率；
求解上述p元一次方程组，共有p个方程，求解bi的值，将第q段缺失频带的中心频率代入上式，获得第q段频段的幅度校正因子
本发明第二方面，提供了一种改进的基于de模型的频带拼接装置，包括：
散射中心位置估计单元，用于采用旋转不变技术估计信号参数方法联合各频段数据估计散射中心位置；
散射幅度系数估计单元，用于使用第一段频带的数据估计散射中心幅度系数；
幅度校正因子计算单元，用于计算第二段频带至第p段频带的幅度校正因子，进行p-1次多项式拟合，将缺失频段中心频率代入获得缺失频带的幅度校正因子；
频带拼接单元，用于根据估计的散射中心位置、缺失频带的幅度校正因子和散射中心幅度系数生成缺失频带数据，完成频带拼接。
在根据本发明的改进的基于de模型的频带拼接装置中，优选地，所述频带拼接单元通过以下公式生成缺失频带数据：
其中，yq(n)为待估计的第q段缺失频带中第n个频点对应的雷达回波数据，fqn对应着待估计的第q段缺失频带的第n个频点，1≤q≤p-1；
本发明第三方面，提供了一种执行改进的基于de模型的频带拼接方法的设备，包括：至少一个处理器、至少一个存储器以及存储在存储器中的计算机程序指令，当计算机程序指令被处理器执行时实现如前所述的方法。
本发明第四方面，提供了一种计算机可读存储介质，其上存储有计算机程序指令，其特征在于，当所述计算机程序指令被处理器执行时实现如前所述的方法。
实施本发明的改进的基于de模型的频带拼接方法、装置、设备及计算机可读存储介质，具有以下有益效果：本发明对de模型进行改进，引入幅度校正因子来估计de模型的参数，重构缺失频带数据，完成频带拼接，相对于基于gtd模型的频带拼接方法而言对数据的拟合更为准确，相对于传统的de模型频带拼接方法而言对大带宽的数据的拼接效果更好。
附图说明
图1为根据本发明第一实施例的改进的基于de模型的频带拼接方法的流程图；
图2为根据本发明第二实施例的改进的基于de模型的频带拼接方法的流程图；
图3为根据本发明优选实施例的改进的基于de模型的频带拼接装置的模块框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明对de模型进行改进，引入幅度校正因子来估计de模型的参数。首先将各段频带数据联合起来估计散射中心的位置参数，再分别估计各个频带数据散射中心的幅度以及幅度校正因子。接下来估计出缺失频带的幅度校正因子，重构缺失频带数据，完成频带拼接。
请参阅图1，为根据本发明第一实施例的改进的基于de模型的频带拼接方法的流程图。如图1所示，该实施例提供改进的基于de模型的频带拼接方法包括以下步骤：
首先，在步骤s101中，执行散射中心位置估计步骤，采用旋转不变技术估计信号参数方法(esprit)联合各频段数据估计散射中心位置。
随后，在步骤s102中，执行散射幅度系数估计步骤，使用第一段频带的数据估计散射中心幅度系数。
随后，在步骤s103中，执行幅度校正因子计算步骤，计算第二段频带至第p段频带的幅度校正因子，进行p-1次多项式拟合，将缺失频段中心频率代入获得缺失频带的幅度校正因子。
最后，在步骤s104中，执行频带拼接步骤，根据估计的散射中心位置、缺失频带的幅度校正因子和散射中心幅度系数生成缺失频带数据，完成频带拼接。
本发明采用改进的de模型，引入幅度校正因子来估计de模型的参数，从而完成频带拼接，可应用于雷达目标特征提取方面。
基于现有的de模型，认为雷达后向散射回波观测序列可以表示为：
其中，y(n)代表雷达的第n个频点对应的回波观测数据；m为散射中心个数；am是第m个散射中心的复散射强度系数；rm表示第m个散射中心相对于零相位参考面的距离；fn是频率，c代表光速；en为测量噪声，这里假设为独立的高斯白噪声。
为了表示散射中心幅度随频率的变化关系，本发明引入了幅度校正因子k，由于散射中心幅度系数的模值随频率变化明显，而幅度系数的相位与频率相关性较小，因此，设定k为实数。则p段带宽、频率间隔均相等的同一目标不同频带的回波数据可表示为：
其中，y1(n)，y2(n)，…，yp(n)数据的中心频率依次增加且频带互不相邻，k1m＝1。
将上式改写为矩阵形式：
y＝as+e(3)
其中，y＝[y1,y2,…,yp]，yi＝[yi(1),yi(2),…,yi(n)]t，y1(n)为第i段频带的第n个频点所对应的回波数据；
a＝[a(r1)t,a(r2)t,…a(rm)t]，a(rm)为第m个散射中心对应的导向矢量：a(rm)＝[1,exp(wm),exp(2wm),…,exp(nwm)]，wm＝-j4πnδfrm/c；δf为步进频率间隔，n为第1段频带的回波数据总数，c代表光速；
s＝[s1,s2,…sp]，
e＝[e1(n),e2(n),…,ep(n)]；ei(n)＝[ei1,ei2,…ein]t,ein为第i段频带的第n个频点所对应的测量噪声。
因此，本发明基于上述改进的de模型，即公式(2)完成频带拼接。
请结合参阅图2，为根据本发明第二实施例的改进的基于de模型的频带拼接方法的流程图。如图2所示，该第二实施例提供的改进的基于de模型的频带拼接方法具体包括以下步骤s201～s206：
其中，步骤s201～s203可以任意顺序先后执行或者同步执行：
在步骤s201中，执行散射中心位置估计步骤，采用旋转不变技术估计信号参数方法(esprit)联合各频段数据估计散射中心位置。
优选地，该步骤中联合各频段数据生成信号子空间es：
对ψ进行特征值分解，求得散射中心位置rm的估计：
其中e1表示es去掉最后一行后形成的矩阵，e2表示es去掉第一行后形成的矩阵，λm表示ψ的第m个特征值，angle(·)表示求相角运算，c代表光速，δf为步进频率间隔。
在步骤s202中，执行散射幅度系数估计步骤，使用第一段频带的数据估计散射中心幅度系数。
优选地，该步骤中由最小二乘法得到对第一段频带的幅度系数向量amp1的估计为：
y1＝[y1(1),y1(2),…,y1(n)]，y1(n)为第一段频带的回波数据，1≤n≤n，n为第一段频带的回波总数，a＝[a(r1)t,a(r2)t,…a(rm)t]，a(rm)为第一段频带的第m个散射中心对应的导向矢量；
根据amp1＝[a1,a2,…,am]估计散射中心幅度系数即第1段频带的数据中第m个散射中心的散射复幅度am，1≤m≤m。
在步骤s203～s204中，执行，计算第二段频带至第p段频带的幅度校正因子，进行p-1次多项式拟合，将缺失频段中心频率代入获得缺失频带的幅度校正因子。
其中，在步骤s203中，计算第二段频带至第p段频带的幅度校正因子；
优选地，该步骤s203用最小二乘法求第二段频带至第p段频带的幅度系数向量ampi的估计，i＝2,…,p；
其中，a＝[a(r1)t,a(r2)t,…a(rm)t]，a(rm)为第一段频带的第m个散射中心对应的导向矢量，yi＝[yi(1),yi(2),…,yi(n)]，yi(n)为第i段频带的第n个频点的回波数据；
由此上式获得对幅度校正因子的估计：
kim为第i段频带的第m个散射中心的幅度校正因子。
其中，在步骤s204中，进行p-1次多项式拟合，将缺失频段中心频率代入获得缺失频带的幅度校正因子
优选地，该步骤s204对第m个散射中心，估计该散射中心在缺失频带的幅度校正因子，p段数据的幅度校正因子组成向量k＝[k1m,k2m,…kpm],对k进行p-1次多项式拟合：
其中，fci为第i段频带的中心频率；
求解上述p元一次方程组，共有p个方程，求解bi的值，将第q段缺失频带的中心频率代入上式，获得第q段频段的幅度校正因子
随后，在步骤s205和步骤s206中，执行频带拼接步骤；
在步骤s205中，根据估计的散射中心位置、缺失频带的幅度校正因子和散射中心幅度系数生成缺失频带数据；
优选地，前述步骤已经得到已缺失频段的散射中心位置rm，幅度校正因子
其中，yq(n)为待估计的第q段缺失频带中第n个频点对应的雷达回波数据，fqn对应着待估计的第q段缺失频带的第n个频点，1≤q≤p-1，这里假设观测雷达为频率步进雷达，则fqn＝fqc+n·δf(1≤n≤n),fqc为待估计的第q段缺失频带的中心频率，δf为步进频率间隔，n为频点总数。
基于相同的发明构思，本发明还提供了一种改进的基于de模型的频带拼接装置。请参阅图3，为根据本发明优选实施例的改进的基于de模型的频带拼接装置的模块框图。如图3所示，该实施例提供的装置300包括：散射中心位置估计单元301、散射幅度系数估计单元302、幅度校正因子计算单元303和频带拼接单元304。
其中，散射中心位置估计单元301用于采用旋转不变技术估计信号参数方法联合各频段数据估计散射中心位置。该散射中心位置估计单元301的处理过程与前述方法中散射中心位置估计步骤一致，在此不再赘述。
散射幅度系数估计单元302用于使用第一段频带的数据估计散射中心幅度系数。该散射幅度系数估计单元302的处理过程与前述方法中散射幅度系数估计步骤一致，在此不再赘述。
幅度校正因子计算单元303用于计算第二段频带至第p段频带的幅度校正因子，进行p-1次多项式拟合，将缺失频段中心频率代入获得缺失频带的幅度校正因子。该幅度校正因子计算单元303的处理过程与前述方法中幅度校正因子计算步骤一致，在此不再赘述。
频带拼接单元304用于根据估计的散射中心位置、缺失频带的幅度校正因子和散射中心幅度系数生成缺失频带数据，完成频带拼接。该频带拼接单元304的处理过程与前述方法中频带拼接步骤一致，在此不再赘述。
本发明实施例还提供了一种执行改进的基于de模型的频带拼接方法的设备，包括：至少一个处理器、至少一个存储器以及存储在存储器中的计算机程序指令，当计算机程序指令被处理器执行时实现如上述实施方式中改进的基于de模型的频带拼接方法。
本发明实施例提供了一种计算机可读存储介质，其上存储有计算机程序指令，当计算机程序指令被处理器执行时实现如上述实施方式中改进的基于de模型的频带拼接方法。
综上所述，本发明相对于基于gtd模型的频带拼接方法而言对数据的拟合更为准确，相对于传统的de模型频带拼接方法而言，能够避免传统de模型由于忽视了散射中心随频率的变化规律而难以对大带宽数据进行准确拼接的缺点，使得对于大带宽的数据的拼接效果更好。
应该理解地是，本发明中改进的基于de模型的频带拼接方法及装置的原理相同，因此对改进的基于de模型的频带拼接方法的实施例的详细阐述也适用于改进的基于de模型的频带拼接装置。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。