CN102798855A - 一种基于数字电视信号的直升机目标识别方法 - Google Patents

一种基于数字电视信号的直升机目标识别方法 Download PDF

Info

Publication number
CN102798855A
CN102798855A CN2012102832268A CN201210283226A CN102798855A CN 102798855 A CN102798855 A CN 102798855A CN 2012102832268 A CN2012102832268 A CN 2012102832268A CN 201210283226 A CN201210283226 A CN 201210283226A CN 102798855 A CN102798855 A CN 102798855A
Authority
CN
China
Prior art keywords
target
helicopter
signal
doppler
radar
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.)
Granted
Application number
CN2012102832268A
Other languages
English (en)
Other versions
CN102798855B (zh
Inventor
单涛
陶然
刘升恒
卓智海
张建
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201210283226.8A priority Critical patent/CN102798855B/zh
Publication of CN102798855A publication Critical patent/CN102798855A/zh
Application granted granted Critical
Publication of CN102798855B publication Critical patent/CN102798855B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及利用基于数字电视信号的外辐射源雷达检测识别直升机目标的方法,属于雷达目标探测处理领域。首先将目标回波信号与直达波信号进行自适应对消;然后对直达波信号及对消后所得的目标回波信号做互模糊函数,在互模糊函数结果中寻找直升机目标的峰值位置,然后沿其多普勒轴检测是否有周期出现的谐波,根据谐波的周期性间隔可得到直升机的螺旋桨转速;根据主峰的多普勒,计算直升机的径向飞行速度;根据主峰所在时延单元计算目标的时延,从而可实现对直升机的检测与参数测量。与传统外辐射源雷达只进行运动目标检测相比,本发明利用螺旋桨微多普勒调制引起的多次谐波副峰特征进行直升机目标识别,可以估计直升机螺旋桨转速等附加信息。

Description

一种基于数字电视信号的直升机目标识别方法
技术领域
本发明涉及一种基于数字电视信号的直升机目标识别方法,特别涉及一种外辐射源雷达中直升机检测、识别的方法,属于雷达目标探测技术领域。
背景技术
外辐射源雷达是一种双多基地雷达系统,利用非合作照射源进行目标探测,由于本身不发射信号,因而具有较低廉的成本。常用的照射源包括调频广播信号、GSM信号、模拟电视、数字电视信号等。其中数字电视信号是外辐射源雷达的理想照射源,数字电视采用数字调制技术,发射功率稳定;同时数字电视信号的带宽较宽(7.56MHz),使外辐射源雷达具有比较高的距离分辨率,有利于探测低空飞行的目标,因此本发明选取数字电视信号作为外辐射源。
如图1所示,本外辐射源雷达利用数字电视广播发射塔作为非合作发射源,通过参考信号接收天线接收直达波信号,通过回波信号接收天线接收目标回波信号,但同时接收到直达波和杂波干扰,因此外辐射源雷达信号处理机首先采用自适应对消算法抑制目标回波接收天线接收到的直达波和杂波干扰,然后利用互模糊函数实现相干积累,并进行恒虚警检测,根据互模糊函数中峰值位置测量目标回波信号与直达波信号的时延差和多普勒频移,从而实现对运动目标的检测和定位。
随着低空空域的逐步开放,航模直升机、小型民用直升机、武装直升机等各类直升机目标对城市的安全造成了巨大的潜在威胁,因此亟需开发一种适用于城市环境的检测识别直升机的雷达技术。
由于电视塔高,电视信号低空覆盖良好,同时外辐射源雷达工作时本身不发射电磁信号,在城市环境具有良好的电磁兼容性。另一方面,外辐射源雷达与脉冲体制雷达在检测直升机方面也具有一定的优势。直升机与普通固定翼飞机的一个主要区别在于其具有螺旋桨,对脉冲体制雷达,由于雷达发射脉冲的速度和直升机螺旋桨转速没有固定的关系,并且仅在雷达发射的电磁波垂直于螺旋桨时回波才强,所以螺旋桨的回波是随机的。而对数字电视信号外辐射源雷达,由于电视塔发射的是连续波信号,所以每次螺旋桨转向雷达时都会有回波,因此在各种检测识别直升机的雷达中,基于数字电视信号的外辐射源雷达有其独特的优势。
直升机螺旋桨的旋转在雷达回波中产生了额外的频率调制,表现在频谱上的特征就是在机身回波主峰周围产生谐波谱线,这就是直升机的微多普勒特征。假设电视塔发射波长为λ的连续波信号,则接收到的理想直升机目标回波信号可以表示为:
s r ( t ) ρexp [ j ( 4 π λ v h t + R + h sin β ) ] × Σ k = 1 K Lh k ( t ) · exp [ j 4 π λ · L 2 cos β sin ( ωt + k K 2 π ) ] - - - ( 1 )
其中,ρ为目标的电磁反射率,λ为外辐射源发射信号波长,vh为直升机相对于雷达的径向速度,R和h分别是直升机螺旋桨中心距雷达的初始距离和螺旋桨初始高度,β为直升机仰角,L为直升机旋翼长度,K为直升机螺旋桨叶片数,ω为旋翼转速,hk(t)为第k个旋翼叶片的周期雷达散射截面(RCS)波动且
h k ( t ) = sin c [ 4 π λ · L 2 cos β sin ( ωt + k K 2 π ) ] - - - ( 2 )
由Victor C.Chen的“Radar Signatures of Rotor Blades”一文,对式(1)做傅里叶变换,其频谱可表示为
S r ( f ) = Σ k = 1 K C k [ δ ( f - hKω / 2 π ) + δ ( f + kKω / 2 π ) ] - - - ( 3 )
其中系数Ck由ω、β、λ、K及L等参数共同确定,这是典型的第一类Bessel函数,由此可以看出,直升机旋翼的微多普勒效应使得在中心频率左右将会对称周期性地出现以ω/2π为周期的谐波谱线。
由于外辐射源雷达接收到的目标回波信号十分微弱,因此在航模直升机这类低慢小目标的预警场合,采用较长时间的相参积累以提高目标的信噪比,并提取目标的微多普勒特征是十分重要的。对此,可采用外辐射源雷达中常用的互模糊函数进行相参积累,实现直升机微多普勒特征的提取。
发明内容
本发明的目的是在实现传统外辐射源雷达的目标检测功能以外,充分利用直升机的微多普勒信息,在频域进行目标识别,并估计直升机螺旋桨转速信息,为此提出一种基于数字电视信号的直升机目标识别方法,通过对互模糊函数相参积累后的目标时延、多普勒二维相关处理结果进行特征提取分析,实现直升机目标识别。
本发明的目的是通过以下技术方案实现的。
本发明提出的一种基于数字电视信号的直升机目标识别方法,其所依托的外辐射源雷达硬件平台包括参考信号接收天线、目标回波接收天线、信号处理机,待测目标为直升机,所述外辐射源雷达通过参考信号接收天线接收直达波信号,通过目标回波接收天线接收目标反射的电视信号;直升机目标识别方法的具体步骤为:
1)外辐射源雷达将通过目标回波接收天线接收到的目标回波信号y(n)与通过参考信号接收天线接收到的直达波信号x(n)在信号处理机中进行自适应对消得到信号e(n),且
e(n)=y(n)-wT(n)x(n)                                    (4)
其中w(n)为滤波器系数,n表示时域样本数;
2)信号处理机利用互模糊函数将直达波信号x(n)与步骤1)得到的对消输出信号e(n)进行二维相关处理,得到互模糊函数结果|x(d,k)|2,且
|x(d,k)|2=|FFT{e(n)x*(n-d)}|2                          (5)
其中,d表示时延单元数,k表示多普勒频移单元数;
3)信号处理机对互模糊函数结果|x(d,k)|2进行恒虚警检测;
4)信号处理机分析步骤3)得到的恒虚警检测结果,观察是否检测到目标,根据检测结果分为以下两种情况:
4.1如果目标存在,则进一步观察是否存在沿多普勒轴的与某一中心频率对称的呈周期出现的多次谐波副峰,此时有以下两种情况:
4.1.1若存在沿多普勒轴的与某一中心频率对称的呈周期出现的多次谐波副峰,则存在直升机目标,然后提取多次谐波副峰间重复周期,计算直升机螺旋桨转速,提取主峰的多普勒频率,计算直升机的径向飞行速度;
4.1.2若不存在沿多普勒轴的与某一中心频率对称的呈周期出现的多次谐波副峰,则按一般非直升机目标处理;
4.2若不存在目标,则返回步骤1)继续接收数据,直到本次任务完成。
有益效果
与传统外辐射源雷达只进行运动目标检测相比,本发明利用螺旋桨微多普勒调制引起的多次谐波副峰特征进行直升机目标识别,可以估计直升机螺旋桨转速等附加信息。
附图说明
图1为典型外辐射源雷达系统配置示意图;
图2为基于外辐射源雷达的直升机检测算法流程图;
图3为有直升机目标存在时,通过互模糊函数得到的相参积累结果;
图4为图3中相参积累结果的多普勒维截面。
具体实施方式
下面结合附图和实施例对本发明做进一步说明。
实施例
本实施例利用仿真的手段来验证本发明的有效性。
仿真参数:某数字电视广播塔发射信号带宽为7.56MHz,信号的中心频率为794MHz,雷达信号处理机基带采样率为9MHz。测试目标为一小型航模直升机,主螺旋桨叶片共两片,翼展0.6m,转速为每分钟2100~2200转。
一种基于数字电视信号的直升机目标识别方法,方法的步骤如图2所示,在本实施例中具体为:
1)外辐射源雷达将通过目标回波接收天线接收到的目标直升机信号y(n)与通过参考信号接收天线接收到的直达波信号x(n)在信号处理机中进行自适应对消,对消滤波器阶数为1200,得到信号e(n);
2)信号处理机利用互模糊函数将直达波信号x(n)与步骤1)得到的对消输出信号e(n)进行二维相关处理,积累时间为0.4s,待积累数据长度为3.6M点,得到|x(d,k)|2,相参积累结果如图3所示;
3)信号处理机对互模糊函数结果|x(d,k)|2进行恒虚警检测;
4)信号处理机分析步骤3)得到的恒虚警检测结果,观测起始时延为1,最大时延为40,观测多普勒范围为200,观察到沿多普勒维存在与某一中心频率对称且呈周期出现的多次谐波副峰,因此判定存在直升机目标。
为了更清楚地观测具体测量结果,取图3的多普勒维横截面,如图4所示,其中第101号多普勒单元对应0多普勒速度。根据主峰的位置,得到目标的延时单元数为16,对应直升机的时延为16/(9×106)=1.78(微秒);直升机机身在第62-101=-39多普勒单元内,谐波副峰重复周期为174-146=28个多普勒单元,多普勒速度为28/0.4=70Hz。直升机机身的径向多普勒为-39/0.4=-97.5Hz,方向为远离雷达运动。由于直升机主螺旋桨叶片为2片,则得到直升机螺旋桨的转速为每秒35转,即每分钟2100转,而由直升机的说明书,其旋翼转速为每分钟2100~2200转,完全符合。
通过数据分析可以看出,基于数字电视信号的外辐射源雷达可以精确地检测到直升机目标,并且测出直升机的螺旋桨转速及径向飞行速度。
以上所述的具体描述,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种基于数字电视信号的直升机目标识别方法,其所依托的外辐射源雷达硬件平台包括参考信号接收天线、目标回波接收天线、信号处理机,待测目标为直升机,所述外辐射源雷达通过参考信号接收天线接收直达波信号,通过目标回波接收天线接收目标反射的电视信号,其特征在于,直升机目标识别方法的具体步骤为:
1)外辐射源雷达将通过目标回波接收天线接收到的目标回波信号y(n)与通过参考信号接收天线接收到的的直达波信号x(n)在信号处理机中进行自适应对消得到信号e(n),且
e(n)=y(n)-wT(n)x(n)                                    (1)
其中w(n)为滤波器系数,n表示时域样本数;
2)信号处理机利用互模糊函数将直达波信号x(n)与步骤1)得到的对消输出信号e(n)进行二维相关处理,得到互模糊函数结果|x(d,k)|2,且
|x(d,k)|2=|FFT{e(n)x*(n-d)}|2                          (2)
其中,d表示时延单元数,k表示多普勒频移单元数;
3)信号处理机对互模糊函数结果|x(d,k)|2进行恒虚警检测;
4)信号处理机分析步骤3)得到的恒虚警检测结果,观察是否检测到目标,根据检测结果分为以下两种情况:
4.1如果目标存在,则进一步观察是否存在沿多普勒轴的与某一中心频率对称的呈周期出现的多次谐波副峰,此时有以下两种情况:
4.1.1若存在沿多普勒轴的与某一中心频率对称的呈周期出现的多次谐波副峰,则存在直升机目标,然后提取多次谐波副峰间重复周期,计算直升机螺旋桨转速,提取主峰的多普勒频率,计算直升机的径向飞行速度;
4.1.2若不存在沿多普勒轴的与某一中心频率对称的呈周期出现的多次谐波副峰,则按一般非直升机目标处理;
4.2若不存在目标,则返回步骤1)继续接收数据,直到本次任务完成。
CN201210283226.8A 2012-08-09 2012-08-09 一种基于数字电视信号的直升机目标识别方法 Expired - Fee Related CN102798855B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210283226.8A CN102798855B (zh) 2012-08-09 2012-08-09 一种基于数字电视信号的直升机目标识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210283226.8A CN102798855B (zh) 2012-08-09 2012-08-09 一种基于数字电视信号的直升机目标识别方法

Publications (2)

Publication Number Publication Date
CN102798855A true CN102798855A (zh) 2012-11-28
CN102798855B CN102798855B (zh) 2014-01-29

Family

ID=47198011

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210283226.8A Expired - Fee Related CN102798855B (zh) 2012-08-09 2012-08-09 一种基于数字电视信号的直升机目标识别方法

Country Status (1)

Country Link
CN (1) CN102798855B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103197300A (zh) * 2013-03-26 2013-07-10 中国科学院电子学研究所 一种基于gpu的外辐射源雷达直达波杂波对消实时处理方法
CN104077787A (zh) * 2014-07-08 2014-10-01 西安电子科技大学 一种基于时域和多普勒域的飞机目标分类方法
CN104133211A (zh) * 2014-07-07 2014-11-05 中国船舶重工集团公司第七二四研究所 一种多普勒频率变换雷达目标分类识别方法
CN105022038A (zh) * 2015-08-07 2015-11-04 武汉大学 一种基于形态分量分析的外辐射源雷达风场杂波抑制方法
CN105093215A (zh) * 2015-08-31 2015-11-25 西安电子科技大学 基于多普勒信息的雷达对低空慢速小目标的跟踪方法
CN106291484A (zh) * 2016-10-13 2017-01-04 武汉大学 基于副峰识别的外辐射源雷达虚假目标抑制方法及装置
CN107049235A (zh) * 2013-01-05 2017-08-18 高通股份有限公司 处理皮肤电导信号以减轻噪声和检测信号特征
CN107728119A (zh) * 2017-10-18 2018-02-23 潘卫军 机场终端区无人机预警与反制系统及方法
CN108287341A (zh) * 2017-12-15 2018-07-17 空间防务技术(深圳)有限公司 一种低空目标探测方法及其系统、电子设备
CN113238288A (zh) * 2021-05-20 2021-08-10 桂林电子科技大学 一种基于差值谱线的旋翼目标特征提取方法
CN115436905A (zh) * 2022-09-30 2022-12-06 中国科学院空天信息创新研究院 基于无源接收系统的基线估计方法、系统及电子设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6633254B1 (en) * 2002-08-15 2003-10-14 Bae Systems Information And Electronics Systems Integration Inc. Self-modulating remote moving target detector
CN101208983A (zh) * 2005-06-28 2008-06-25 索尼爱立信移动通讯股份有限公司 用于根据联播通信信号来确定移动装置的位置的方法、系统和装置
CN101414866A (zh) * 2007-10-17 2009-04-22 中卫视讯卫星科技(北京)有限公司 移动目标卫星跟踪定位、识别、语音视频、数据传输终端机

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6633254B1 (en) * 2002-08-15 2003-10-14 Bae Systems Information And Electronics Systems Integration Inc. Self-modulating remote moving target detector
CN101208983A (zh) * 2005-06-28 2008-06-25 索尼爱立信移动通讯股份有限公司 用于根据联播通信信号来确定移动装置的位置的方法、系统和装置
CN101414866A (zh) * 2007-10-17 2009-04-22 中卫视讯卫星科技(北京)有限公司 移动目标卫星跟踪定位、识别、语音视频、数据传输终端机

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
万显荣等: "基于低频段数字广播电视信号的外辐射源雷达发展现状与趋势", 《雷达学报》 *
鞠晓燕等: "基于数字电视的无源雷达信号建模与分析", 《无线电通讯技木》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107049235A (zh) * 2013-01-05 2017-08-18 高通股份有限公司 处理皮肤电导信号以减轻噪声和检测信号特征
CN103197300B (zh) * 2013-03-26 2015-07-01 中国科学院电子学研究所 一种基于gpu的外辐射源雷达直达波杂波对消实时处理方法
CN103197300A (zh) * 2013-03-26 2013-07-10 中国科学院电子学研究所 一种基于gpu的外辐射源雷达直达波杂波对消实时处理方法
CN104133211A (zh) * 2014-07-07 2014-11-05 中国船舶重工集团公司第七二四研究所 一种多普勒频率变换雷达目标分类识别方法
CN104077787A (zh) * 2014-07-08 2014-10-01 西安电子科技大学 一种基于时域和多普勒域的飞机目标分类方法
CN104077787B (zh) * 2014-07-08 2017-03-08 西安电子科技大学 一种基于时域和多普勒域的飞机目标分类方法
CN105022038A (zh) * 2015-08-07 2015-11-04 武汉大学 一种基于形态分量分析的外辐射源雷达风场杂波抑制方法
CN105093215A (zh) * 2015-08-31 2015-11-25 西安电子科技大学 基于多普勒信息的雷达对低空慢速小目标的跟踪方法
CN106291484A (zh) * 2016-10-13 2017-01-04 武汉大学 基于副峰识别的外辐射源雷达虚假目标抑制方法及装置
CN106291484B (zh) * 2016-10-13 2018-11-23 武汉大学 基于副峰识别的外辐射源雷达虚假目标抑制方法及装置
CN107728119A (zh) * 2017-10-18 2018-02-23 潘卫军 机场终端区无人机预警与反制系统及方法
CN108287341A (zh) * 2017-12-15 2018-07-17 空间防务技术(深圳)有限公司 一种低空目标探测方法及其系统、电子设备
CN113238288A (zh) * 2021-05-20 2021-08-10 桂林电子科技大学 一种基于差值谱线的旋翼目标特征提取方法
CN113238288B (zh) * 2021-05-20 2022-07-01 桂林电子科技大学 一种基于差值谱线的旋翼目标特征提取方法
CN115436905A (zh) * 2022-09-30 2022-12-06 中国科学院空天信息创新研究院 基于无源接收系统的基线估计方法、系统及电子设备

Also Published As

Publication number Publication date
CN102798855B (zh) 2014-01-29

Similar Documents

Publication Publication Date Title
CN102798855B (zh) 一种基于数字电视信号的直升机目标识别方法
CN108957443B (zh) 一种基于双发双收相干雷达的无人机旋翼长度和转速的估计方法
CN103176178B (zh) 雷达动目标Radon-分数阶傅里叶变换长时间相参积累检测方法
CN103323829B (zh) 基于Radon-分数阶模糊函数的雷达动目标长时间相参积累检测方法
CN102707269B (zh) 一种机载雷达距离走动校正方法
CN104898103B (zh) 基于多通道杂波图的低速目标检测方法
CN103207387B (zh) 一种机载相控阵pd雷达杂波的快速模拟方法
CN102749621B (zh) 一种双基地合成孔径雷达频域成像方法
US20120105272A1 (en) Method for filtering the radar echoes produced by wind turbines
Kong et al. Wind turbine radar interference studies by polarimetric measurements of a scaled model
CN104076362B (zh) 一种铁路列车运行安全防控雷达
CN104730518A (zh) 一种基于高斯拟合的雷达多普勒谱估计海面流场的方法
CN102650689A (zh) 一种步进频率脉冲雷达位移测量方法
JP2014182010A (ja) レーダ装置
Zhang et al. Extraction of vibrating features with dual-channel fixed-receiver bistatic SAR
CN110082729A (zh) 天气雷达补盲平滑处理系统及方法
CN103064084A (zh) 基于距离频域的解模糊方法
CN102937709B (zh) 基于直达波的双基地sar无同步回波聚焦成像方法
CN102830394A (zh) 基于多谱线积累的弱目标探测方法
Shi et al. Experimental demonstration for ionospheric sensing and aircraft detection with a HF skywave multistatic radar
CN104391288B (zh) 拟合干涉相位的虚警剔除方法
Wu et al. Parameter estimation for SAR moving target detection using fractional Fourier transform
CN105487071A (zh) 一种基于步进频雷达的多运动目标信号处理方法
Tikkinen et al. Utilization of long coherent integration time in helicopter recognition by passive coherent location (PCL) radar
Xiaoliang et al. Analysis on the range of wind farm's impact area for secondary surveillance radar

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140129

Termination date: 20140809

EXPY Termination of patent right or utility model