CN110907912B - 谐振区目标窄带模式下的极点特征求解方法 - Google Patents
谐振区目标窄带模式下的极点特征求解方法 Download PDFInfo
- Publication number
- CN110907912B CN110907912B CN201911208345.5A CN201911208345A CN110907912B CN 110907912 B CN110907912 B CN 110907912B CN 201911208345 A CN201911208345 A CN 201911208345A CN 110907912 B CN110907912 B CN 110907912B
- Authority
- CN
- China
- Prior art keywords
- frequency domain
- target
- pole
- polynomial
- domain response
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/418—Theoretical aspects
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种谐振区目标窄带模式下的极点特征求解方法,所述方法包括如下步骤:步骤一、建立电磁响应传输函数并利用稀疏信号表示多项式主极点数学模型;步骤二、采用二阶锥寻优框架构造系数向量,通过调整窄带频域响应采样点数,计算频域响应偏差合理选定滤波器阶数;步骤三、利用系数向量计算多项式主极点,选定第四象限极点作为目标特征的主极点,并拟合目标频域响应。本发明利用很少的RCS频域数据,构造系数向量拟合多点RCS回波特征,在此基础上,提取表征目标结构的主极点信息,可用于提高目标分类识别概率。同时,引入凸优化寻优策略增强对频域响应鲁棒自适应性,由于频域采样数据量减少,有效提高了系统的实时性和环境适应性。
Description
技术领域
本发明属于雷达目标特性和识别技术领域,涉及一种谐振区目标窄带模式下的极点特征求解方法。
背景技术
极点是谐振区雷达目标识别的重要特征,且对方位不敏感。然而极点的提取需要宽带信号的时域晚期响应或者频域响应,这在实际系统中很难满足。在目标识别的研究过程中,人们多年来都困惑于复杂目标的散射特性随着极化、方位等因素剧烈变化。而将极点用于目标识别则可以克服这些缺点。当目标的尺寸L和雷达波长λ相当时(λ≤L<10λ),目标处在谐振区。极点是谐振区目标电磁散射的重要特征,它受目标本身的特性如形状、尺寸、材料等决定,与目标姿态以及雷达极化方式无关,将极点应用于目标识别可以克服极化、方位改变导致的剧烈变化。
目前极点提取的方式主要分为时域法和频域法,时域法主要利用E脉冲法、Prony法,频域法主要包含MPM和Cauchy方法,其中MPM利用频域数据通过逆傅立叶变换获得时域晚期响应进行求解,而Cauchy直接利用频域数据通过矩阵分解获得传输函数系数,进而获得极点信息。因为时域法对于雷达回波信号提取难度较大,目前,主要利用频域法进行极点提取。然而,频域法通常需要在很宽的频域范围内对目标回波进行采样,因此需要大量不同频点回波数据,这对于窄带雷达系统来说难以满足。
从理论上来说,极点数目有无穷多个,但一般只需考虑影响目标谐振的有限个主极点。一般来说,起主要作用的具有较大留数的极点称为主极点。其中距虚轴较近的极点又是决定中、长时间响应的主要因素,由于越靠近虚轴,衰减越缓慢。而目标的主体结构又是由距离实轴较近的极点决定,越靠近实轴,振荡频率越低。所以靠近原点而且留数较大的极点最能反映目标的结构特征。这对于窄带雷达系统如何通过更少的频域数据提取目标的特征提供了解决方向。
发明内容
为了解决背景技术中存在的问题,特别是针对窄带雷达系统利用更少的频域响应提取目标极点特征的技术问题,本发明提供了一种谐振区目标窄带模式下的极点特征求解方法。该方法能够利用很少的RCS频域数据,构造系数向量拟合多点RCS回波特征,在此基础上,提取表征目标结构的主极点信息,可用于提高目标分类识别概率。同时,引入凸优化寻优策略增强对频域响应鲁棒自适应性,由于频域采样数据量减少,有效提高了系统的实时性和环境适应性。
本发明的目的是通过以下技术方案实现的:
一种谐振区目标窄带模式下的极点特征求解方法,包括如下步骤:
步骤一、建立电磁响应传输函数并利用稀疏信号表示多项式主极点数学模型,分析电磁响应传输函数系数满足稀疏条件以及采样矩阵满足重构条件。
本步骤中,晚期响应可以表征目标的结构特性,并且能够利用一系列衰减正弦信号和的形式表示其电磁响应传输函数h(t),电磁响应传输函数h(t)可以转化为频域响应函数。电磁响应传输函数h(t)表示为:
本步骤中,多项式求解可表示为矩阵形式,并可转化为稀疏信号的求解模型进行计算,则多项式主极点数学模型表示为:
式中,ai是多项式分子系数,bj是多项式分母系数,s是频率点(即:采样频率),H(s)是频域响应函数,i=0,1,…,N,j=1,2…,M,N是H(s)的分子多项式阶数,M是H(s)分母多项式阶数(即:滤波器阶数),通常M=N+1,b0是分母多项式常数项。
本步骤中,采样矩阵表示为:
Ax=y;
x=[a0…aN b1…bM]T;
y=b0[H(s1)…H(sK)]T;
K是频域数据采样点数。
本步骤中,从采样矩阵中抽取的每K个列向量构成的矩阵是非奇异的。采样矩阵列是非线性相关的,满足重构条件。
步骤二、采用二阶锥寻优框架构造系数向量
xj=[a0…aN b1…bM]T,计算在滤波器阶数j=1,2,…,M和系数xj下的频域响应Hj(s),通过调整窄带频域响应采样点数,计算频域响应偏差Bias(j)合理选定滤波器阶数,如果Bias(j)≤0.001,则选定阶数j;否则重复步骤一。
本步骤中,求解系数向量是范数求解问题,既是非线性的也是非二次的,将其转换为二阶锥优化问题求解,可以通过调整Newton迭代次数控制内点法计算复杂度,并且迭代次数独立于求解问题本身,保持相对较低的计算量。
本步骤中,二阶锥寻优框架的公式如下:
式中,zk是残差,λ是正则化参数,K是频域响应采样点数。
式中,aj是多项式分子系数,bj是多项式分母系数,s表示采样频率,j=1,…,M选定的滤波器阶数。
本步骤中,频域响应偏差Bias(j)的计算公式如下:
本步骤中,系数的数量与选择滤波器的阶数成正比。阶数太小,传输函数不能拟合窄带数据;阶数太大,矩阵A是病态矩阵,无法得到正确的解。选择合理的阈值,可以将系数向量等效为稀疏信号。
步骤三、利用系数向量xj=[a0…aN b1…bM]T计算多项式主极点,因为主极点σl+jωl的实部σl是衰减因子,表征电磁波沿目标表面的能量损失,是负值;ωl是谐振角频率,是正值;因此,选定第四象限极点作为目标特征的主极点,并拟合目标频域响应。
本步骤中,通过系数向量,求解多项式主极点及谐振频率,拟合目标的频域响应,进行目标的分类识别。
相比于现有技术,本发明具有如下优点:
1、本发明采用的是目标频域数据,对实际雷达系统而言更容易获得,可替代实际雷达测试或微波暗室缩比试验,能够节约目标识别成本。
2、本发明基于窄带雷达频域数据实现,需要的频域响应数据量大大减少,符合实际雷达系统工作条件,实时性更好。
3、本发明基于稀疏信号重构理论,利用二阶锥寻优手段提取主极点特征信息,相比于同类方法鲁棒性更好。
附图说明
图1是本发明谐振区目标窄带模式下的极点特征求解方法(SRM)的处理流程框图;
图2是本发明SRM方法与现有MPM、Cauchy方法的处理流程对比框图;
图3是本发明细导线目标的RCS频域响应;
图4是本发明细导线目标在第一个谐振点附近窄带等间隔采样点;
图5是本发明利用SRM方法针对图3采样点构造的系数向量;
图6是本发明利用SRM方法针对图4系数向量恢复的窄带频域响应与原始频域响应的拟合结果对比;
图7是本发明细导线目标针对第二个谐振点采样构造的系数向量;
图8是本发明利用SRM方法针对图6系数向量恢复的窄带频域响应与原始频域响应的拟合结果对比;
图9是本发明细导线目标极点特征提取结果和现有技术提取结果的比较。
具体实施方式
下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。
本发明提供了一种谐振区目标窄带模式下的极点特征提取方法(SRM),如图1所示,所述方法包括以下步骤:
步骤一、选定滤波器阶数j=j0。
步骤二、构造采样矩阵A和频域响应数据y。
在谐振区,目标被电磁波照射后,目标谐振现象产生的反射电磁波叫做晚期响应。晚期响应可以表征目标的结构特性,并且能够利用一系列衰减正弦信号和的形式表示其电磁响应传输函数h(t)。
对电磁响应传输函数h(t)做Laplace变换,即:
式中,l=1,2,…,L表示极点个数,上角标‘*’表示共轭;sl是第l个极点,Rl是第l个极点对应的留数。
时域信号构成是复指数函数,因此转换后的频域系统是线性时不变系统,且传输函数是以极点和零点为特征表示的。
式中,ai表示多项式分子系数,bj表示多项式分母系数,s表示采样频率。
令M=N+1表示滤波器阶数。
转换成矩阵的形式为:
Ax=y (5);
x=[a0…aN b1…bM]T;
y=b0[H(s1)…H(sK)]T;
K表示频域数据采样点数,b0为常数项。
计算系数向量x=[a0…aN b1…bM]T后,通过公式(3)可以获得多项式极点分布,并且由于系数向量是稀疏的,那么频域响应y=b0[H(s1)…H(sK)]T就只需要很少的几个频点就可以实现极点的重构。
对窄带信号,公式(5)是一个欠定方程,即K<M+N+1。此时,如果待恢复信号满足稀疏表示条件,则可以较大概率通过求解欠定方程来获得。
步骤三、按照公式(6)二阶锥寻优框架,计算出系数向量xj=[a0…aN b1…bM]T;
步骤五、滤波器阶数选择,计算传输函数与频域响应的偏差来判断阶数选择的合理性。利用公式(7)计算频域响应偏差,如果Bias(j)≤0.001,则选定阶数j;否则重复步骤一;
步骤六:根据公式(3)计算主极点以及拟合的窄带频域响应。
为了验证构造的系数向量是否满足稀疏条件,我们利用更多频域采样数据获得满秩矩阵进行精确求解,首先对公式(5)中的采样矩阵进行分块,
可知:[R22][b]=0,[a]=-[R11]-1[R12][b]。对[R22]做SVD分解[R22]=[U][Σ][V]H,b的解向量可以表示为[b]=[V],进一步可确定a。从a,b的解可以确定x的分布特点。计算结果如表1所示,分子多项式很小,分母多项式系数中只有低阶的几个值较大。因此,可以将分子和分母组成的系数向量当作稀疏信号进行表示,在这个前提下,目标的回波特性可以利用稀疏表示理论进行求解。
表1细导线计算系数表
为了验证本发明可以利用更少的频域响应采样点获得同现有技术可比拟极点提取精度,对细导线目标进行了极点提取仿真实验。
1、细导线目标极点提取与验证
细导线长度L=0.15m(L/D=150,D为细导线的直径),频率范围取50MHz~7GHz,共计696个频点。矩量法计算得到细导线的回波频域数据如图3所示,由于细导线在设定频率范围内具有明显的谐振特点,从频域数据上可以清晰的找到4个谐振频率点。取窄带频率范围为840MHz~1.34GHz(对应角频率0.53e10~0.84e10),对于SRM方法,按间隔100MHz抽取频点数据,即得到11个频点RCS数据,如图4所示。
因为MPM和Cauchy方法要求更多的频域采样数据点,在这里最优阶数为10,按照前述理论分析,MPM和Cauchy方法要求采样频点数量要大于等于20个。SRM利用抽取的11个频点计算系数值如图5所示,从系数值上可以看出分子多项式系数很小,分母多项式系数相对较大,这有利于反映极点的特征信息。将计算的系数值代入公式(3)得到反演的RCS曲线,如图6所示。从结果上可以看出,在窄带限定范围内(两条绿色虚线之间),计算得到的RCS曲线能够非常准确的拟合原始的RCS曲线;在窄带边界区域,也能够符合RCS的变化趋势,实际上这与抽选的频点数量是相关的,当抽选数量较多时,边界拟合延伸性更好,而抽选数量不足时,边界拟合越差;在离窄带边界更远的区域,完全失真。
同样,获取窄带范围为2.74GHz~3.24GHz(对应角频率1.72e10~2.03e10)可获得第二个谐振点的拟合曲线,如图8所示。
2.与MPM、Cauchy方法比较
利用MPM、Cauchy和SRM三种方法对细导线提取极点,得到结果如图9所示。SRM与其他两种方法极点提取结果是一致的,但是SRM提取每个极点用了11个采样点,MPM和Cauchy提取每个极点用了51个采样点数据。因此,本发明提出的SRM方法在满足相同重构精度的前提下,需要的频域响应数量远低于现有技术方法。
通过对谐振区目标窄带模式下的极点提取结果、现有文献结果和理论结果对比验证,说明了本发明方法在满足极点提取精度的前提下,相对需要更少的频域响应采样点。
Claims (3)
1.一种谐振区目标窄带模式下的极点特征求解方法,其特征在于所述方法包括如下步骤:
步骤一、建立电磁响应传输函数并利用稀疏信号表示多项式主极点数学模型,分析电磁响应传输函数系数满足稀疏条件以及采样矩阵满足重构条件,其中:
所述电磁响应传输函数h(t)表示为:
所述多项式主极点数学模型表示为:
式中,ai是多项式分子系数,bj是多项式分母系数,s是频率点,H(s)是频域响应函数,i=0,1,…,N,j=1,2…,M,N是H(s)的分子多项式阶数,M是H(s)分母多项式阶数,即:滤波器阶数,b0是分母多项式常数项;
所述采样矩阵表示为:
Ax=y;
x=[a0…aN b1…bM]T;
y=b0[H(s1)…H(sK)]T;
K是频域数据采样点数;
步骤二、采用二阶锥寻优框架构造系数向量xj=[a0…aN b1…bM]T,计算在滤波器阶数j=1,2,…,M和系数xj下的频域响应Hj(s),通过调整窄带频域响应采样点数,计算频域响应偏差Bias(j)合理选定滤波器阶数,如果Bias(j)≤0.001,则选定阶数j;否则重复步骤一,其中:
所述二阶锥寻优框架的公式如下:
式中,zk是残差,λ是正则化参数;
步骤三、利用系数向量xj=[a0…aN b1…bM]T计算多项式主极点,选定第四象限极点作为目标特征的主极点,并拟合目标频域响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911208345.5A CN110907912B (zh) | 2019-11-30 | 2019-11-30 | 谐振区目标窄带模式下的极点特征求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911208345.5A CN110907912B (zh) | 2019-11-30 | 2019-11-30 | 谐振区目标窄带模式下的极点特征求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110907912A CN110907912A (zh) | 2020-03-24 |
CN110907912B true CN110907912B (zh) | 2022-09-30 |
Family
ID=69821386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911208345.5A Active CN110907912B (zh) | 2019-11-30 | 2019-11-30 | 谐振区目标窄带模式下的极点特征求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110907912B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112526476A (zh) * | 2020-12-21 | 2021-03-19 | 上海无线电设备研究所 | 基于谐振散射回波信号的雷达目标识别方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5973998A (en) * | 1997-08-01 | 1999-10-26 | Trilon Technology, Llc. | Automatic real-time gunshot locator and display system |
JP2005167820A (ja) * | 2003-12-04 | 2005-06-23 | Osaka Industrial Promotion Organization | 超音波センサ素子及び超音波アレイセンサ装置の共振周波数の調整方法及び調整装置 |
CN101776748A (zh) * | 2009-01-12 | 2010-07-14 | 哈尔滨工业大学 | 复杂目标极点特征提取方法和装置 |
CN102279394A (zh) * | 2011-05-17 | 2011-12-14 | 西安电子科技大学 | 低频窄带高分辨超声探测成像方法 |
CN104280722A (zh) * | 2014-10-17 | 2015-01-14 | 南京理工大学 | 基于稀疏表示理论的雷达目标极点提取方法 |
EP2872913A1 (en) * | 2012-09-18 | 2015-05-20 | Koninklijke Philips N.V. | Magnetic resonance guided linac |
CN104950936A (zh) * | 2015-07-13 | 2015-09-30 | 浙江工业大学 | 基于稳幅的谐振频率跟踪和振动载荷振幅联合控制系统 |
CN107976659A (zh) * | 2017-10-17 | 2018-05-01 | 上海无线电设备研究所 | 基于S2method时频分析的目标极点特征提取与检测方法 |
-
2019
- 2019-11-30 CN CN201911208345.5A patent/CN110907912B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5973998A (en) * | 1997-08-01 | 1999-10-26 | Trilon Technology, Llc. | Automatic real-time gunshot locator and display system |
JP2005167820A (ja) * | 2003-12-04 | 2005-06-23 | Osaka Industrial Promotion Organization | 超音波センサ素子及び超音波アレイセンサ装置の共振周波数の調整方法及び調整装置 |
CN101776748A (zh) * | 2009-01-12 | 2010-07-14 | 哈尔滨工业大学 | 复杂目标极点特征提取方法和装置 |
CN102279394A (zh) * | 2011-05-17 | 2011-12-14 | 西安电子科技大学 | 低频窄带高分辨超声探测成像方法 |
EP2872913A1 (en) * | 2012-09-18 | 2015-05-20 | Koninklijke Philips N.V. | Magnetic resonance guided linac |
CN104280722A (zh) * | 2014-10-17 | 2015-01-14 | 南京理工大学 | 基于稀疏表示理论的雷达目标极点提取方法 |
CN104950936A (zh) * | 2015-07-13 | 2015-09-30 | 浙江工业大学 | 基于稳幅的谐振频率跟踪和振动载荷振幅联合控制系统 |
CN107976659A (zh) * | 2017-10-17 | 2018-05-01 | 上海无线电设备研究所 | 基于S2method时频分析的目标极点特征提取与检测方法 |
Non-Patent Citations (4)
Title |
---|
"A narrowband high frequency distributed power transformer model for partial discharge location";S.D. Mitchell 等;《2007 Australasian Universities Power Engineering Conference 》;20071231;第1-7页 * |
"Accurate Interpolation of Amplitude-Only Frequency Domain Response Based on an Adaptive Cauchy Method";J. Yang;《IEEE Transactions on Antennas and Propagation》;20161231;第1005-1013页 * |
"基于谐振区频域数据的雷达目标极点提取方法";伍光新 等;《系统工程与电子技术》;20110415;第774-778页 * |
"超宽带雷达目标电磁特征抽取与识别方法研究";王党卫;《中国博士论文全文数据库》;20070515;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110907912A (zh) | 2020-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106646344B (zh) | 一种利用互质阵的波达方向估计方法 | |
CN109683161B (zh) | 一种基于深度admm网络的逆合成孔径雷达成像的方法 | |
CN110045323B (zh) | 一种基于矩阵填充的互质阵稳健自适应波束形成算法 | |
CN108710102B (zh) | 基于互质阵列二阶等价虚拟信号离散傅里叶逆变换的波达方向估计方法 | |
CN106443621B (zh) | 基于正交匹配稀疏重构的相干信源动态doa追踪方法 | |
CN107192878A (zh) | 一种基于压缩感知的电力系统谐波检测方法及装置 | |
CN107843875A (zh) | 基于奇异值分解降噪的贝叶斯压缩感知雷达数据融合方法 | |
Georgiou | Solution of the general moment problem via a one-parameter imbedding | |
CN105137373B (zh) | 一种指数信号的去噪方法 | |
CN104965188A (zh) | 一种阵列误差下的波达方向估计方法 | |
CN109444842A (zh) | 一种目标电磁散射特性数据重构方法和装置 | |
CN103605121B (zh) | 基于快速稀疏贝叶斯学习算法的宽带雷达数据融合方法 | |
CN104780127B (zh) | 一种时延‑多普勒r‑l解卷多路径信道估计方法 | |
CN102879824A (zh) | 一种基于迭代收缩的快速稀疏Radon变换方法 | |
Li et al. | A deep learning method for denoising based on a fast and flexible convolutional neural network | |
CN112560342A (zh) | 一种基于dnn的大气波导参数估计方法 | |
CN110907912B (zh) | 谐振区目标窄带模式下的极点特征求解方法 | |
CN110413939B (zh) | 一种基于原子范数的到达角估计方法 | |
CN104459635B (zh) | 基于迭代收缩加权融合的自适应空滤滤波方法 | |
CN105022025A (zh) | 基于稀疏处理的信号波达方向估计方法 | |
CN113484842B (zh) | 基于目标姿态分帧与散射字典截断的rcs重构方法 | |
CN102509268B (zh) | 基于免疫克隆选择的非下采样轮廓波域图像去噪方法 | |
CN110412535B (zh) | 一种序贯的空时自适应处理参数估计方法 | |
CN107315713B (zh) | 一种基于非局部相似性的一维信号去噪增强方法 | |
CN105426916A (zh) | 图像相似度计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |