CN107632291A - 一种基于极化雷达抗角反射器干扰的方法 - Google Patents
一种基于极化雷达抗角反射器干扰的方法 Download PDFInfo
- Publication number
- CN107632291A CN107632291A CN201710712673.3A CN201710712673A CN107632291A CN 107632291 A CN107632291 A CN 107632291A CN 201710712673 A CN201710712673 A CN 201710712673A CN 107632291 A CN107632291 A CN 107632291A
- Authority
- CN
- China
- Prior art keywords
- radar
- polarization
- polarized
- target
- corner reflector
- 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.)
- Pending
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于极化雷达抗角反射器干扰的方法,包含以下过程:外场对船模型和角反射器进行雷达回波采样,对回波数据进行处理,提取其全极化散射矩阵数据;根据所述全极化散射矩阵数据计算每个频率点下的极化不变量;将所述计算得到的极化不变量数据放在数据库中,作为SVM分类器的训练样本;雷达发射单频信号,极化雷达实时采集回波数据,计算极化特征,形成测试样本,在所述数据库中找到与所述雷达实时采集回波数据频率相对应的训练样本,运用SVM分类器进行鉴别舰船和角反射器。本发明具有受环境因素小,相对极化分解计算量较少,且具有实时性,简单易行,且可以推广到识别其他假目标干扰。
Description
技术领域
本发明涉及全极化体制雷达领域,特别涉及一种极化雷达抗角反射器干扰的方法。
背景技术
雷达目标的极化信息与其结构、材料、形状、姿态取向等有着直接的联系。通过极化信息的提取,可以获取目标表面粗糙度、对称性和取向等其他参数难以表征的特征,是完整刻画目标特性所不可或缺的。
直接利用目标散射矩阵或者从功率矩阵中提取目标的极化特征。由散射矩阵获得的极化参数和极化不变量确定了在给定姿态下目标的后向散射特征,是雷达从视线方向所能获得的目标最大信息量。利用极化参数的变换特征和极化不变量作为特征来识别目标,是一条有效的技术途径。
角反射器干扰在当前的电子对抗系统中应用非常普遍。但是关于如何鉴别舰船和角的文献非常少。一方面一直把角反射器替代舰船当作打靶的工具,并没有把角反射器当作干扰来对待;另一方面,没有找到较合适的方法对抗角反射器。
现有技术中,通过采用散射中心模型,推导了海面目标横摇运动引入的雷达微多普勒信号模型,利用时频分析方法和一阶矩方法,提取了目标回波的瞬时多普勒和微多普勒。通过定义目标雷达回波的微多普勒主频率特征实现了角反射器干扰的有效鉴别。但是此方法受海面环境的影响因素较大。
在现有技术中,还可利用目标宽带条件下极化分解比重邻频差异的统计特性,将Krogager极化分解算法用于舰船和角反射器阵列的鉴别。该方法适用于单个角反射器和角反射器阵列,可以达到较好的鉴别效果。目标分解方法主要适用于目标散射特性差别较大的情况,多用于极化SAR图像分类中;类型相同的目标可能会具有相同的散射分量,很难用简单的方法对其区分;需要配合其他的特征进行识别,且此算法很难做到实时性。
目前,专利方面,电子科技大学提出“一种基于波形设计的主被动雷达协同抗干扰方法”(申请号:CN 201410468299),其是利用雷达波形设计技术,采用主被动雷达协同工作的模式,利用被动雷达实时侦测的回波信息,对干扰进行检测和识别,估计干扰的相关参数,主动雷达利用被动雷达提供的干扰信息,设计最优的发射波形,并对回波信息进行处理。该方法采用主被动雷达同时工作,增加了系统的复杂性;同时,要设计最优发射波形,不能做到实时性。
深圳光启高等理工研究院提出“一种雷达抗干扰的方法”(申请号:CN201110060728),其是在接收雷达的反射面和馈源之间放置一个超材料极化转换器,用以转换发射雷达发出的干扰电磁波的极化方式。然后,接收雷达改用经极化转换器转换后的极化方式通信。该方法适用于两种极化方式不同的雷达之间工作,无法对抗角反射器无源干扰。
合肥佳瑞林电子技术有限公司提出“一种抗干扰雷达”(申请号: CN201610617067),该方法适用于杂波抑制和对经过调制的有源干扰起以一定的抗干扰效果,但是无法对抗角反射器这种无源干扰。
发明内容
本发明目的是提供一种基于极化雷达抗角反射器干扰的方法,通过在极化识别的基础上提出基于目标回波极化特征和宽带雷达目标极化特性相结合的识别方法,此方法首先提取目标回波的极化散射矩阵,通过散射矩阵计算其极化不变量,通过极化不变量初步鉴别船和角反射器,为提高识别率,通过宽带雷达目标的极化特性进行进一步分类识别,从而实现达到快速有效识别角反射器和船的目的。
为了实现以上目的,本发明通过以下技术方案实现:
一种基于极化雷达抗角反射器干扰的方法,包含以下过程:外场对船模型和角反射器进行雷达回波采样,对回波数据进行处理,提取其全极化散射矩阵数据。根据所述全极化散射矩阵数据计算每个频率点下的极化不变量。将所述计算得到的极化不变量数据放在数据库中,作为SVM分类器的训练样本。雷达发射单频信号,极化雷达实时采集回波数据,计算极化特征,形成测试样本,在所述数据库中找到与所述雷达实时采集回波数据频率相对应的训练样本,运用SVM分类器进行鉴别舰船和角反射器。
优选地,所述通过回波数据计算极化散射矩阵进一步包含以下过程:主天线接收电压:
vm1(t)=gm 2(SHH+2εSHV+ε2SVV)χ(t)+nm1(t)
式中,gm为主天线的电压增益,SHH、SVV为极化散射矩阵共极化分量,ε为天线的交叉极化分量,SHV为为雷达主天线发射极化方式为垂直极化且接收的极化方式为水平极化时极化散射矩阵的极化分量,nm1(t)为主天线接收通道的接收机噪声,服从正态分布。
若为垂直极化方式,那么主天线的接收电压为:
式中,Tr为雷达脉冲重复周期,nm2(t)为主天线接收通道的接收机噪声,服从正态分布。
当雷达发射水平极化信号时,辅天线的接收电压:
vc1(t)=gmgc[SHV(1+ε2)+ε(SHV+SVV)]χ(t)+nc1(t)
式中,gc为辅天线的电压方向图,nc1(t)为当雷达发射水平极化信号时,辅助天线接收通道接收机噪声。
在雷达发射垂直极化信号时,辅天线的接收电压:
vc2(t)=gmgc[SVV+2εSHV+ε2SHH]χ(t)+nc2(t)
式中,nc2(t)为当雷达发射垂直极化信号时,辅助天线接收通道接收机噪声。
令mA=SHH-SVV
则
同理,由(1)和(3)得
SHH(1-2ε2)-ε2SVV=mB
极化散射矩阵共极化分量的估计值分别为:
上述极化散射分量SHV的估计值为:
式中, 和为极化散射矩阵共极化分量,为极化散射矩阵交叉极化分量。
根据雷达目标的互易性:SHV=SVH,得到回波全极化散射矩阵。
优选地,所述极化不变量包含:
行列式模值|△|:
散射矩阵的迹S1:
S1=SVV+SHH;
功率矩阵迹P1:
P1=|SHH|2+|SVV|2+2|SHV|2;
功率矩阵行列式值|△|P:
组合不变量系数η:
目标极化形状因子γ:
散射总能量σS:
目标观测长度Ls:
Ls=(max{l|H(l)>q}-min{l|H(l)>q})△r
式中,H(l)为距离像;q为门限阈,其与噪声电平有关;△r为距离分辨单元大小;
频域响应特性为:雷达发射电磁波频率为X波段,船和角反射器在各个极化通道的频域特性;
频率敏感性P与Q:
P=max{|FFT-1[η1(f)]|}
Q=max{|FFT-1[λ2(f)]|}
式中:
优选地,所述特征参量P、Q用于表示雷达目标极化散射随频率的变化特征,度量雷达目标电磁散射的极化色散程度和其极化散射结构的稳定程度。
优选地,对于所述目标分类识别进一步包含以下过程:基于目标的行列式模值、散射矩阵的迹、功率矩阵迹、功率矩阵行列式值、组合不变量系数、目标极化形状因子通过SVM分类器进行分类识别;
当船识别率低于第一预设值,角反识别率大于第二预设值时,则认为雷达距离门内没有目标,无需进行二级分类识别;
当角反射器识别率低于第一预设值,舰船识别率大于第二预设值,则认为雷达距离门内船目标,则无进行二级分类识别;
当角反射器和船识别率大于第一预设值并且小于第二预设值;则进行二级识别。
优选地,所述二级识别是对目标的目标散射总能量、径向长度、频域响应特性、频率敏感性特征参量P、Q进行分类识别。
本发明与现有技术相比具有以下优点:
本发明通过在极化识别的基础上提出基于目标回波极化特征和宽带雷达目标极化特性相结合的识别方法,此方法首先提取目标回波的极化散射矩阵,通过散射矩阵计算其极化不变量,通过极化不变量初步鉴别船和角反射器,为提高识别率,通过宽带雷达目标的极化特性进行进一步分类识别,从而具有能够达到快速有效识别角反射器和船的优点。
本发明比微多普勒鉴别角反射器受环境因素小,相对极化分解计算量较少;且本发明分两级分类易做到实时性,简单易行,且可以推广到识别其他假目标干扰。
附图说明
图1为本发明一种基于极化雷达抗角反射器干扰的方法的流程图。
具体实施方式
以下结合附图,通过详细说明一个较佳的具体实施例,对本发明做进一步阐述。
如图1所示,本发明一种基于极化雷达抗角反射器干扰的方法,包含以下过程:
步骤S1、外场对船模型和角反射器进行雷达回波采样,对回波数据进行处理,提取其全极化散射矩阵数据。
步骤S2、根据上述全极化散射矩阵数据计算每个频率点下的极化不变量;上述极化不变量包含:行列式模值、功率矩阵迹、组合不变量系数、目标极化形状因子、散射总能量、径向长度、频域响应特性、频率敏感性两个参量P 和Q。
步骤S3、将上述计算得到的极化不变量数据放在数据库中,作为SVM分类器的训练样本。
步骤S4、雷达发射单频信号,极化雷达实时采集到在此频率下的回波数据,计算极化特征,形成测试样本,在特征库中找到与上述频率的频点对应的训练样本,运用SVM分类器进行分类识别,从而达到鉴别舰船和角反射器的目的。
通过SVM分类器对测试样本进行分类识别时,进一步包含以下过程:首先基于目标的行列式模值、散射矩阵的迹、功率矩阵迹、功率矩阵行列式值、组合不变量系数、目标极化形状因子运用SVM进行分类识别。如果船识别率低于20%,角反识别率大于80%,则认为雷达距离门内没有舰船目标,不需要进行二级分类识别;如果角反射器识别率低于20%,舰船识别率大于 80%,则认为雷达距离门内有船目标,则不需要进行二级分类识别;如果角反射器和船识别率大于20%并且小于80%;则进行二级识别。二级识别对目标的目标散射总能量、径向长度、频域响应特性、频率敏感性特征参量P、 Q进行分类识别。
在本实施例中,在步骤S1中,外场获得舰船和角反射器的频率处于X 波段,实验发射信号采用步进频体制实现,步进频率为2MHz,俯仰角为 40~60°,方位角为0~180°,角度步进1°的回波数据。
在步骤S2中,所述每个频率点的俯仰角为40~60°,方位为0~180°角度步进1°。
在本实施例中,假设角反射器边长为1m,舰船模型长为10米,宽为3m, 高为1米,角反射器距离船模型5m。且船和角反射器的RCS相近。接收天线放在距离船和角反射器20m处。雷达步进跳频测量系统的工作频率范围为 X波段,跳频间隔为2MHz,步进跳频测量系统与舰船和角反射器连线垂直放置,雷达测量的俯仰角设为40~60°,采样间隔1°,方位角设置从0°~180°,雷达分时发射水平、垂直极化波,接收时两个正交极化通道同时接收信号,同时分别测量HH、HV、VV、VH4种线极化波回波。
在步骤S1中,通过回波数据计算极化散射矩阵,具体计算过程如下:
雷达主天线接收电压:
vm1(t)=gm 2(SHH+2εSHV+ε2SVV)χ(t)+nm1(t) (1)
式中,gm为主天线的电压增益,SHH、SVV为极化散射矩阵共极化分量,ε为天线的交叉极化分量,SHV为为雷达主天线发射极化方式为垂直极化且接收的极化方式为水平极化时极化散射矩阵的极化分量,nm1(t)为主天线接收通道的接收机噪声,服从正态分布。
发射的电磁波若为垂直极化,那么主天线的接收电压为:
式中,Tr为雷达脉冲重复周期,nm2(t)为主天线接收通道的接收机噪声,服从正态分布。
当雷达发射水平极化信号时,辅天线的接收电压:
vc1(t)=gmgc[SHV(1+ε2)+ε(SHV+SVV)]χ(t)+nc1(t) (3)
式中,gc为辅天线的电压方向图,nc1(t)为当雷达发射水平极化信号时,辅助天线接收通道接收机噪声。
在雷达发射垂直极化信号时,辅天线的接收电压:
vc2(t)=gmgc[SVV+2εSHV+ε2SHH]χ(t)+nc2(t) (4)
式中,nc2(t)为当雷达发射垂直极化信号时,辅助天线接收通道接收机噪声。式(1)减式(2)得
令mA=SHH-SVV (6)
则由式5和式6得到:
同理得:
由式6和式8得极化散射矩阵共极化分量的估计值分别为:
极化散射分量SHV的估计值为:
式中, 和为极化散射矩阵共极化分量,为极化散射矩阵交叉极化分量;根据雷达目标的互易性:SHV=SVH,得到回波全极化散射矩阵。
本实施例中,在步骤S2,由极化散射矩阵得到每个频率点和每个角度的特征值即计划不变量进一步为:
行列式模值|△|
散射矩阵的迹S1:
S1=SVV+SHH (13)
功率矩阵迹P1:
P1=|SHH|2+|SVV|2+2|SHV|2 (14)
功率矩阵行列式值|△|P:
组合不变量系数η:
目标极化形状因子γ:
散射总能量σS:
目标观测长度Ls:
Ls=(max{l|H(l)>q}-min{l|H(l)>q})△r (19)
式中,H(l)为距离像;q为门限阈,其与噪声电平有关;△r为距离分辨单元大小。
频域响应特性为:雷达发射电磁波频率为X波段,船和角反射器在各个极化通道的频域特性。
频率敏感性P与Q:
P=max{|FFT-1[η1(f)]|} (20)
Q=max{|FFT-1[λ2(f)]|} (21)
式中:
上述特征参量P、Q侧重反映了雷达目标极化散射随频率的变化特征,度量了雷达目标电磁散射的极化色散程度和其极化散射结构的稳定程度。
在本实施例中,将计算的特征参量放入数据库中。将船的11个特征量作为SVM分类器的正训练样本,将角反射器的11个特征量作为SVM分类器的负训练样本。雷达采集的回波数据的11个特征量作为测试样本。SVM分类器的核函数选择高斯核函数,并选择对应的核参数和惩罚因子。
对于目标分类识别,分两级操作。第一步:先基于目标的行列式模值、散射矩阵的迹、功率矩阵迹、功率矩阵行列式值、组合不变量系数、目标极化形状因子进行分类识别。在实际应用时,雷达发射单频信号,获得在此频率下的回波数据,计算极化特征,形成测试样本,在特征库中找到本频点的训练样本,运用SVM进行分类识别。
如果船识别率低于20%,角反识别率大于80%,则认为雷达距离门内没有目标,不需要进行二级分类识别;如果角反射器识别率低于20%,舰船识别率大于80%,则认为雷达距离门内船目标,则不需要进行二级分类识别;如果角反射器和船识别率大于20%并且小于80%;则进行二级识别。二级识别对目标的目标散射总能量、径向长度、频域响应特性、频率敏感性特征参量P、Q进行分类识别。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (6)
1.一种基于极化雷达抗角反射器干扰的方法,其特征在于,包含以下过程:外场对船模型和角反射器进行雷达回波采样,对回波数据进行处理,提取其全极化散射矩阵数据;
根据所述全极化散射矩阵数据计算每个频率点下的极化不变量;
将所述计算得到的极化不变量数据放在数据库中,作为SVM分类器的训练样本;
雷达发射单频信号,极化雷达实时采集回波数据,计算极化特征,形成测试样本,在所述数据库中找到与所述雷达实时采集回波数据频率相对应的训练样本,运用SVM分类器进行鉴别舰船和角反射器。
2.如权利要求1所述一种基于极化雷达抗角反射器干扰的方法,其特征在于,所述通过回波数据计算极化散射矩阵进一步包含以下过程:
主天线接收电压:
vm1(t)=gm 2(SHH+2εSHV+ε2SVV)χ(t)+nm1(t)
式中,gm为主天线的电压增益,SHH、SVV为极化散射矩阵共极化分量,ε为天线的交叉极化分量,SHV为为雷达主天线发射极化方式为垂直极化且接收的极化方式为水平极化时极化散射矩阵的极化分量,nm1(t)为主天线接收通道的接收机噪声,服从正态分布。
若为垂直极化方式,那么主天线的接收电压为:
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
<mi>A</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>g</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&epsiv;S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>+</mo>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>&chi;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
式中,Tr为雷达脉冲重复周期,nm2(t)为主天线接收通道的接收机噪声,服从正态分布;
当雷达发射水平极化信号时,辅天线的接收电压:
vc1(t)=gmgc[SHV(1+ε2)+ε(SHV+SVV)]χ(t)+nc1(t)
式中,gc为辅天线的电压方向图,nc1(t)为当雷达发射水平极化信号时,辅助天线接收通道接收机噪声;
在雷达发射垂直极化信号时,辅天线的接收电压:
vc2(t)=gmgc[SVV+2εSHV+ε2SHH]χ(t)+nc2(t)
式中,nc2(t)为当雷达发射垂直极化信号时,辅助天线接收通道接收机噪声;
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>m</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
<mi>A</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>g</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>&chi;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,令mA=SHH-SVV,则
令
SHH(1-2ε2)-ε2SVV=mB,则
极化散射矩阵共极化分量的估计值分别为:
<mrow>
<mover>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>^</mo>
</mover>
<mo>&ap;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msubsup>
<mi>g</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mi>&chi;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>m</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&alpha;&epsiv;v</mi>
<mrow>
<mi>c</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>v</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
<mi>A</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mover>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>^</mo>
</mover>
<mo>&ap;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msubsup>
<mi>g</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mi>&chi;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
<mi>A</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&alpha;&epsiv;v</mi>
<mrow>
<mi>c</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
极化散射分量SHV的估计值为:
<mrow>
<mover>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>^</mo>
</mover>
<mo>&ap;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msubsup>
<mi>g</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mi>&chi;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>&alpha;v</mi>
<mrow>
<mi>c</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&epsiv;v</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
<mi>A</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
式中, 和为极化散射矩阵共极化分量,为极化散射矩阵交叉极化分量;
根据雷达目标的互易性:SHV=SVH,得到回波全极化散射矩阵。
3.如权利要求2所述一种基于极化雷达抗角反射器干扰的方法,其特征在于,所述极化不变量包含:
行列式模值|△|:
<mrow>
<mo>|</mo>
<mi>&Delta;</mi>
<mo>|</mo>
<mo>=</mo>
<mo>|</mo>
<mi>d</mi>
<mi>e</mi>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>S</mi>
<mo>&rsqb;</mo>
<mo>|</mo>
<mo>=</mo>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>|</mo>
<mo>;</mo>
</mrow>
散射矩阵的迹S1:
S1=SVV+SHH;
功率矩阵迹P1:
P1=|SHH|2+|SVV|2+2|SHV|2;
功率矩阵行列式值|△|P:
<mrow>
<mo>|</mo>
<mi>&Delta;</mi>
<msub>
<mo>|</mo>
<mi>P</mi>
</msub>
<mo>=</mo>
<mo>|</mo>
<mi>d</mi>
<mi>e</mi>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>G</mi>
<mo>&rsqb;</mo>
<mo>|</mo>
<mo>=</mo>
<mo>|</mo>
<mi>d</mi>
<mi>e</mi>
<mi>s</mi>
<mo>&lsqb;</mo>
<msup>
<mi>S</mi>
<mi>H</mi>
</msup>
<mi>S</mi>
<mo>&rsqb;</mo>
<mo>|</mo>
<mo>=</mo>
<mo>|</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>H</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>|</mo>
<mo>;</mo>
</mrow>
组合不变量系数η:
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>|</mo>
<mi>&Delta;</mi>
<mo>|</mo>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>2</mn>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>|</mo>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
目标极化形状因子γ:
<mrow>
<mi>&gamma;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<msubsup>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<msubsup>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
散射总能量σS:
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>S</mi>
</msub>
<mo>=</mo>
<munder>
<mo>&Integral;</mo>
<mi>&Omega;</mi>
</munder>
<mo>&lsqb;</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mi>d</mi>
<mi>f</mi>
<mo>;</mo>
</mrow>
目标观测长度Ls:
Ls=(max{l|H(l)>q}-min{l|H(l)>q})△r
式中,H(l)为距离像;q为门限阈,其与噪声电平有关;△r为距离分辨单元大小;
频域响应特性为:雷达发射电磁波频率为X波段,船和角反射器在各个极化通道的频域特性;
频率敏感性P与Q:
P=max{|FFT-1[η1(f)]|}
Q=max{|FFT-1[λ2(f)]|}
式中:
4.如权利要求3所述一种基于极化雷达抗角反射器干扰的方法,其特征在于,所述特征参量P、Q用于表示雷达目标极化散射随频率的变化特征,度量雷达目标电磁散射的极化色散程度和其极化散射结构的稳定程度。
5.如权利要求4所述一种基于极化雷达抗角反射器干扰的方法,其特征在于,
对于所述目标分类识别进一步包含以下过程:基于目标的行列式模值、散射矩阵的迹、功率矩阵迹、功率矩阵行列式值、组合不变量系数、目标极化形状因子通过SVM分类器进行分类识别;
当船识别率低于第一预设值,角反识别率大于第二预设值时,则认为雷达距离门内没有目标,无需进行二级分类识别;
当角反射器识别率低于第一预设值,舰船识别率大于第二预设值,则认为雷达距离门内船目标,则无进行二级分类识别;
当角反射器和船识别率大于第一预设值并且小于第二预设值;则进行二级识别。
6.如权利要求5所述一种基于极化雷达抗角反射器干扰的方法,其特征在于,所述二级识别是对目标的目标散射总能量、径向长度、频域响应特性、频率敏感性特征参量P、Q进行分类识别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710712673.3A CN107632291A (zh) | 2017-08-18 | 2017-08-18 | 一种基于极化雷达抗角反射器干扰的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710712673.3A CN107632291A (zh) | 2017-08-18 | 2017-08-18 | 一种基于极化雷达抗角反射器干扰的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107632291A true CN107632291A (zh) | 2018-01-26 |
Family
ID=61100552
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710712673.3A Pending CN107632291A (zh) | 2017-08-18 | 2017-08-18 | 一种基于极化雷达抗角反射器干扰的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107632291A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110109068A (zh) * | 2019-04-24 | 2019-08-09 | 西安电子科技大学 | 基于极化特征向量的箔条干扰识别方法 |
CN110857979A (zh) * | 2018-08-26 | 2020-03-03 | 吴成彬 | 一种基于相参信号的极化频选特性材料抗干扰辨识系统及方法 |
CN110865340A (zh) * | 2019-12-04 | 2020-03-06 | 上海无线电设备研究所 | 一种基于极化特性辅助的海面角反射器干扰对抗方法 |
CN111427035A (zh) * | 2020-03-25 | 2020-07-17 | 湖南国科锐承电子科技有限公司 | 一种基于子带微多普勒差异的雷达实测数据集扩充方法 |
CN113608184A (zh) * | 2021-08-04 | 2021-11-05 | 上海无线电设备研究所 | 一种相控阵天线发射自检方法 |
CN113759325A (zh) * | 2021-07-29 | 2021-12-07 | 西安电子科技大学 | 基于极化不变量和极化分解的海面角反射器干扰鉴别方法 |
CN114076924A (zh) * | 2020-08-20 | 2022-02-22 | 西安电子科技大学 | 一种基于多频回波数据实现大型复杂目标分类识别的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101996405A (zh) * | 2010-08-30 | 2011-03-30 | 中国科学院计算技术研究所 | 一种快速的玻璃图像缺陷检测及分类方法及其装置 |
WO2012038781A1 (en) * | 2010-09-23 | 2012-03-29 | Nokia Corporation | State change sensing based on individual location patterns |
EP2527791A1 (en) * | 2011-05-25 | 2012-11-28 | CSR Technology Holdings Inc. | Hierarchical context detection method to determine location of a mobile device on a person's body |
CN102831442A (zh) * | 2011-06-13 | 2012-12-19 | 索尼公司 | 异常行为检测设备和方法及生成该检测设备的设备和方法 |
CN105160353A (zh) * | 2015-08-18 | 2015-12-16 | 西安电子科技大学 | 基于多特征集的极化sar数据地物分类方法 |
-
2017
- 2017-08-18 CN CN201710712673.3A patent/CN107632291A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101996405A (zh) * | 2010-08-30 | 2011-03-30 | 中国科学院计算技术研究所 | 一种快速的玻璃图像缺陷检测及分类方法及其装置 |
WO2012038781A1 (en) * | 2010-09-23 | 2012-03-29 | Nokia Corporation | State change sensing based on individual location patterns |
EP2527791A1 (en) * | 2011-05-25 | 2012-11-28 | CSR Technology Holdings Inc. | Hierarchical context detection method to determine location of a mobile device on a person's body |
CN102831442A (zh) * | 2011-06-13 | 2012-12-19 | 索尼公司 | 异常行为检测设备和方法及生成该检测设备的设备和方法 |
CN105160353A (zh) * | 2015-08-18 | 2015-12-16 | 西安电子科技大学 | 基于多特征集的极化sar数据地物分类方法 |
Non-Patent Citations (3)
Title |
---|
朱珍珍 等: "基于极化分解的舰船和角反射器鉴别方法", 《舰船电子对抗》 * |
李永祯 等: "基于辅助天线的目标极化散射矩阵估算方法", 《现代雷达》 * |
王福友 等: "基于极化不变量特征的雷达目标识别技术", 《雷达科学与技术》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110857979A (zh) * | 2018-08-26 | 2020-03-03 | 吴成彬 | 一种基于相参信号的极化频选特性材料抗干扰辨识系统及方法 |
CN110857979B (zh) * | 2018-08-26 | 2024-05-14 | 吴成彬 | 一种基于相参信号的极化频选特性材料抗干扰辨识系统及方法 |
CN110109068A (zh) * | 2019-04-24 | 2019-08-09 | 西安电子科技大学 | 基于极化特征向量的箔条干扰识别方法 |
CN110109068B (zh) * | 2019-04-24 | 2022-12-02 | 西安电子科技大学 | 基于极化特征向量的箔条干扰识别方法 |
CN110865340A (zh) * | 2019-12-04 | 2020-03-06 | 上海无线电设备研究所 | 一种基于极化特性辅助的海面角反射器干扰对抗方法 |
CN110865340B (zh) * | 2019-12-04 | 2022-03-29 | 上海无线电设备研究所 | 一种基于极化特性辅助的海面角反射器干扰对抗方法 |
CN111427035A (zh) * | 2020-03-25 | 2020-07-17 | 湖南国科锐承电子科技有限公司 | 一种基于子带微多普勒差异的雷达实测数据集扩充方法 |
CN114076924A (zh) * | 2020-08-20 | 2022-02-22 | 西安电子科技大学 | 一种基于多频回波数据实现大型复杂目标分类识别的方法 |
CN113759325A (zh) * | 2021-07-29 | 2021-12-07 | 西安电子科技大学 | 基于极化不变量和极化分解的海面角反射器干扰鉴别方法 |
CN113759325B (zh) * | 2021-07-29 | 2023-12-26 | 西安电子科技大学 | 基于极化不变量和极化分解的海面角反射器干扰鉴别方法 |
CN113608184A (zh) * | 2021-08-04 | 2021-11-05 | 上海无线电设备研究所 | 一种相控阵天线发射自检方法 |
CN113608184B (zh) * | 2021-08-04 | 2023-09-22 | 上海无线电设备研究所 | 一种相控阵天线发射自检方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107632291A (zh) | 一种基于极化雷达抗角反射器干扰的方法 | |
Shui et al. | Tri-feature-based detection of floating small targets in sea clutter | |
Hu et al. | Optimal signal processing in ground-based forward scatter micro radars | |
CN106707247B (zh) | 一种基于紧凑天线阵的高频海洋雷达目标检测方法 | |
CN111220955B (zh) | 基于垂直阵列协方差矩阵特征分解的机载气象雷达地杂波抑制方法 | |
CN112068119A (zh) | 一种被动雷达导引头对雷达和诱饵干扰的识别方法 | |
Zhan et al. | Recognition and mitigation of micro-Doppler clutter in radar systems via support vector machine | |
Chen et al. | Radar signal processing for low-observable marine target-challenges and solutions | |
CN108872947B (zh) | 一种基于子空间技术的海杂波抑制方法 | |
CN102269807A (zh) | 基于ais信息的雷达抗干扰技术的实现方法 | |
Hu et al. | Statistic characteristic analysis of forward scattering surface clutter in bistatic radar | |
Chen et al. | A method against DRFM dense false target jamming based on jamming recognization | |
Fang et al. | E 2 DTF: An End-to-End Detection and Tracking Framework for Multiple Micro-UAVs With FMCW-MIMO Radar | |
CN109061598A (zh) | 一种stap杂波协方差矩阵估计方法 | |
CN109581366B (zh) | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 | |
Pasculli et al. | DVB-T passive radar tracking on real data using Extended Kalman Filter with DOA estimation | |
Shen et al. | Efficient adaptive approach for airborne radar short-range clutter suppression | |
CN114325599B (zh) | 一种针对不同环境的自动门限检测方法 | |
CN110488239A (zh) | 基于调频连续波雷达的目标检测方法 | |
CN112684425A (zh) | 一种恒虚警检测后的目标二次甄别方法 | |
Latthe et al. | A review on ISAR imaging techniques for low RCS targets | |
CN113759325B (zh) | 基于极化不变量和极化分解的海面角反射器干扰鉴别方法 | |
Gao et al. | Estimating the Location of HFSWR Extended Targets Based on Clustering Algorithm | |
Zhang et al. | Modified Space-Time Adaptive Processing with first-order bragg lines kept in HFSWR | |
Björklund et al. | Clutter properties for STAP with smooth and faceted cylindrical conformal antennas |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180126 |