CN105445716B - 一种基于旋翼回波信号时域自相关的悬停直升机检测方法 - Google Patents
一种基于旋翼回波信号时域自相关的悬停直升机检测方法 Download PDFInfo
- Publication number
- CN105445716B CN105445716B CN201510831078.2A CN201510831078A CN105445716B CN 105445716 B CN105445716 B CN 105445716B CN 201510831078 A CN201510831078 A CN 201510831078A CN 105445716 B CN105445716 B CN 105445716B
- Authority
- CN
- China
- Prior art keywords
- signal
- time domain
- data
- sequence
- rotary blade
- 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
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/411—Identification of targets based on measurements of radar reflectivity
- G01S7/412—Identification of targets based on measurements of radar reflectivity based on a comparison between measured values and known or stored values
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于旋翼回波信号时域自相关的悬停直升机检测方法,先对回波信号数据进行中频采样,将采样后的数据进行模数转换,并对转换后的数字信号进行数字正交化处理,转换为基带信号,计算基带信号的模值,进行恒虚警检测,根据恒虚警检测结果进行时域窗位置搜索,检测出闪烁脉冲信号数据,提取闪烁脉冲信号,并计算闪烁周期,将提取后的闪烁脉冲信号组合成新的序列,对闪烁脉冲信号新序列做自相关处理,得到自相关序列,对自相关序列进行恒虚警检测,根据恒虚警检测结果,完成慢速或悬停直升机检测识别。本发明设计简单,计算量小,通用性好,易于工程实现。
Description
技术领域
本发明涉及一种基于旋翼回波信号时域自相关的悬停直升机检测方法,特别适用于单一频域检测体制的脉冲多普勒雷达系统对悬停直升机的检测识别。
背景技术
传统脉冲多普勒雷达系统的直升机回波中频采样主要由机身回波、旋翼回波、地杂波和接收机热噪声四部分构成。对于运动直升机,其机身回波存在多普勒,雷达系统利用多普勒信息,提取制导信息。对于慢速或悬停直升机探测的主要难点在于直升机飞行高度低,其机身回波的多普勒频率接近于零频与地杂波叠加在一起,受地海杂波影响较大,传统脉冲多普勒雷达系统无法利用机身回波多普勒信息进行检测识别,导致对于慢速或悬停直升机脉冲多普勒雷达系统探测性能严重下降。
发明内容
本发明提供一种基于旋翼回波信号时域自相关的悬停直升机检测方法,设计简单,计算量小,通用性好,易于工程实现。
为了达到上述目的,本发明提供一种基于旋翼回波信号时域自相关的悬停直升机检测方法,包含以下步骤:
步骤S1、对回波信号数据进行中频采样,将采样后的数据进行模数转换,并对转换后的数字信号进行数字正交化处理,转换为基带信号;
步骤S2、计算基带信号的模值,对模值进行恒虚警检测;
步骤S3、根据恒虚警检测结果进行时域窗位置搜索,检测出闪烁脉冲信号数据;
步骤S4、提取闪烁脉冲信号,并计算闪烁周期,将提取后的闪烁脉冲信号按照出现的先后顺序组合成新的序列;
步骤S5、对闪烁脉冲信号新序列做自相关处理,得到自相关序列;
步骤S6、对自相关序列进行恒虚警检测;
步骤S7、根据恒虚警检测结果,完成慢速或悬停直升机检测识别。
所述的步骤S3中,在进行时域窗位置搜索时,每个时域窗包含3个子窗,子窗长度为3帧数据长度,子窗间隔为11帧数据长度,数据窗滑动5次,每次滑动2帧数据长度,每次统计窗内闪烁信号出现的次数,闪烁信号出现最多的窗所在的位置即为闪烁信号出现的位置。
所述的步骤S5中,在自相关处理中,对于周期性信号序列,自相关函数定义为:
公式中:x(m)为自相关序列,N为数据长度,k为延迟点数;
自相关函数具有周期性,即,R(k)=R(k+Tp),其中,Tp为信号周期。
所述的步骤S5中,在频域计算自相关序列:先将原序列做快速傅里叶变换,同时将原序列逆序后做快速傅里叶变换,接着将两个快速傅里叶变换的结果点乘,然后将结果做快速傅里叶逆变换,即可得到自相关序列。
所述的步骤S6中,对自相关序列进行恒虚警检测,当自相关序列中的数据信号超过恒虚警门限时,记录该数据信号。
所述的步骤S7中,将恒虚警检测获得的数据信号与恒虚警检测中的参考噪声门限比较,超过参考噪声门限的数据信号则认为是直升机目标。
本发明设计简单,计算量小,通用性好,易于工程实现。
附图说明
图1是本发明的流程图。
图2是某型号直升机旋翼回波经过数字正交化处理后的基带信号。
图3是旋翼回波基带信号频谱。
图4是进行二进制检测即时域窗位置搜索的示意图。
具体实施方式
以下根据图1~图4,具体说明本发明的较佳实施例。
普通雷达目标回波能量在时域上是持续的,在多普勒频域上是积累集中的,与此相反,直升机旋翼回波在时域上表现为Sinc函数形状的周期闪烁脉冲,是周期性、小时宽闪烁脉冲,脉峰过后迅速衰减,在下一个闪烁脉冲出现之前几乎没有旋翼回波能量,而在多普勒频域是分散的,存在展宽,作相参积累可能性不大。就旋翼回波检测而言,慢速或悬停直升机机身的雷达回波可被视为目标杂波。旋翼回波时域检测时,信号具有最大瞬时功率,具有一定的优势,若能有效滤除地杂波和目标杂波,可以提高旋翼回波检测时的信杂比,提高旋翼回波检测概率。
如图1所示,本发明提供一种基于旋翼回波信号时域自相关的悬停直升机检测方法,包含以下步骤:
步骤S1、对回波信号数据进行中频采样,将采样后的数据进行模数转换,并对转换后的数字信号进行数字正交化处理,转换为基带信号;
图2为某型号直升机旋翼回波经过数字正交化处理后的基带信号,由图2可以看出旋翼回波为周期性小时宽闪烁脉冲,脉峰过后迅速衰减,在下一个闪烁脉冲出现之前几乎没有旋翼回波能量,对慢速或悬停直升机旋翼回波信号的检测应有一定的观测时间,以保证观测期间内至少2~3次有效闪烁脉冲;
图3为旋翼回波基带信号频谱,由于旋翼旋转的调制作用,直升机旋翼回波脉冲在时间上呈现sin c函数调制形状,其瞬时频谱将被展宽,当叶片朝向雷达方向运动时,形成正的多普勒谱,而当叶片背向雷达方向运动时,形成负的多普勒谱,自然,当叶片数目为偶数时,处在脉峰期间的回波脉冲将同时具有正、负两个多普勒谱;
本实施例中,采用中频接收机进行中频采样,采用模拟数字转换器ADC进行模数转换,中频接收机由微波-前中组件和中频接收信道共同构成,微波-前中组件的作用为将接收到的高频回波信号经过频率搬移,转换为中频信号,中频接收信道的作用为将接收的中频信号放大、滤波,使中频信号获得适当的增益并抑制带外干扰,中频接收信道本身具有一定的增益,同时在中频接收机的输出端采用数控衰减器,通过控制衰减量,使中频输出信号始终保持在模拟数字转换器ADC的动态范围内;
步骤S2、计算基带信号的模值,然后进行恒虚警检测;
将回波基带信号采样数据排列成恒虚警检测器所需的形式,接着将回波基带信号的I/Q数据取模值,然后对模值进行恒虚警检测,当被检测量超过恒虚警门限(门限设置为10dB)时,认为闪烁信号存在,记为“1”,同时记录闪烁信号的幅值和位置;
步骤S3、根据恒虚警检测结果进行时域窗位置搜索,检测出闪烁脉冲信号数据;
如图4所示,为进行二进制检测即时域窗位置搜索的示意图,每个时域窗包含3个子窗,子窗长度为3帧数据长度,子窗间隔为11帧数据长度,数据窗滑动5次,每次滑动2帧数据长度,每次统计窗内闪烁信号出现的次数,闪烁信号出现最多的窗所在的位置即为闪烁信号出现的位置;
步骤S4、提取闪烁脉冲信号,并计算闪烁周期,将提取后的闪烁脉冲信号按照出现的先后顺序组合成新的序列;
最多提取3个闪烁脉冲数据;
步骤S5、对闪烁脉冲信号新序列做自相关处理,得到自相关序列;
对于周期性信号序列,自相关函数定义为:
公式中:
x(m)为自相关序列,N为数据长度,k为延迟点数。
自相关反映了信号和其自身作了一段延迟后的相似程度,一个周期信号的自相关函数也是周期的,并且和原信号的周期相同,即R(k)=R(k+Tp),Tp为信号周期,由于调制回波基带信号具有闪烁周期,故其自相关函数具有周期性;
为了减少计算量,将在频域计算自相关序列,即,先将原序列做FFT(快速傅里叶变换),同时将原序列逆序后做FFT,两个FFT结果点乘,然后将结果做IFFT(快速傅里叶逆变换),即可得到自相关序列;
根据基带信号自相关结果可以看出,闪烁脉冲之间存在一定的相关性,而闪烁脉冲与噪声并没有相关性,故而提高了检测信噪比,有利于信号检测,从自相关波形中还可以看出,自相关函数具有明显的周期性,且周期等于基带信号闪烁周期;
步骤S6、对自相关序列进行恒虚警检测(门限设置为10dB);
对自相关序列进行恒虚警检测,当自相关序列中的数据信号超过恒虚警门限(门限设置为10dB)时,记录该数据信号;
步骤S7、根据恒虚警检测结果,完成慢速或悬停直升机检测识别;
将步骤S6中恒虚警检测获得的数据信号与恒虚警检测中的参考噪声门限比较,超过参考噪声门限的数据信号则认为是直升机目标。本发明采用外环式速度跟踪回路,通过控制VCO来控制速度门的位置,可将目标和地杂波滤除,中频模拟信号数字化接收技术既避免了传统模拟接收机固有的温度漂移、增益变化和直流电平漂移等问题,还具备硬件可编程、参数可配置、调试方便等优点,本发明设计简单,计算量小,通用性好,易于工程实现,是一种适合于脉冲多普勒雷达系统使用的慢速或悬停直升机检测技术。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (6)
1.一种基于旋翼回波信号时域自相关的悬停直升机检测方法,其特征在于,包含以下步骤:
步骤S1、对回波信号数据进行中频采样,将采样后的数据进行模数转换,并对转换后的数字信号进行数字正交化处理,转换为基带信号;
步骤S2、计算基带信号的模值,对模值进行恒虚警检测;
步骤S3、根据恒虚警检测结果进行时域窗位置搜索,检测出闪烁脉冲信号数据;
步骤S4、提取闪烁脉冲信号,并计算闪烁周期,将提取后的闪烁脉冲信号按照出现的先后顺序组合成新的序列;
步骤S5、对闪烁脉冲信号新序列做自相关处理,得到自相关序列;
步骤S6、对自相关序列进行恒虚警检测;
步骤S7、根据恒虚警检测结果,完成慢速或悬停直升机检测识别。
2.如权利要求1所述的基于旋翼回波信号时域自相关的悬停直升机检测方法,其特征在于,所述的步骤S3中,在进行时域窗位置搜索时,每个时域窗包含3个子窗,子窗长度为3帧数据长度,子窗间隔为11帧数据长度,数据窗滑动5次,每次滑动2帧数据长度,每次统计窗内闪烁信号出现的次数,闪烁信号出现最多的窗所在的位置即为闪烁信号出现的位置。
3.如权利要求1所述的基于旋翼回波信号时域自相关的悬停直升机检测方法,其特征在于,所述的步骤S5中,在自相关处理中,对于周期性信号序列,自相关函数定义为:
公式中:x(m)为自相关序列,N为数据长度,k为延迟点数;
自相关函数具有周期性,即,R(k)=R(k+Tp),其中,Tp为信号周期。
4.如权利要求3所述的基于旋翼回波信号时域自相关的悬停直升机检测方法,其特征在于,所述的步骤S5中,在频域计算自相关序列:先将原序列做快速傅里叶变换,同时将原序列逆序后做快速傅里叶变换,接着将两个快速傅里叶变换的结果点乘,然后将结果做快速傅里叶逆变换,即可得到自相关序列。
5.如权利要求1所述的基于旋翼回波信号时域自相关的悬停直升机检测方法,其特征在于,所述的步骤S6中,对自相关序列进行恒虚警检测,当自相关序列中的数据信号超过恒虚警门限时,记录该数据信号。
6.如权利要求1所述的基于旋翼回波信号时域自相关的悬停直升机检测方法,其特征在于,所述的步骤S7中,将恒虚警检测获得的数据信号与恒虚警检测中的参考噪声门限比较,超过参考噪声门限的数据信号则认为是直升机目标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510831078.2A CN105445716B (zh) | 2015-11-25 | 2015-11-25 | 一种基于旋翼回波信号时域自相关的悬停直升机检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510831078.2A CN105445716B (zh) | 2015-11-25 | 2015-11-25 | 一种基于旋翼回波信号时域自相关的悬停直升机检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105445716A CN105445716A (zh) | 2016-03-30 |
CN105445716B true CN105445716B (zh) | 2018-04-17 |
Family
ID=55556135
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510831078.2A Active CN105445716B (zh) | 2015-11-25 | 2015-11-25 | 一种基于旋翼回波信号时域自相关的悬停直升机检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105445716B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106597411B (zh) * | 2016-12-30 | 2019-06-04 | 无锡市雷华科技有限公司 | 雷达信号处理方法 |
CN111953624A (zh) * | 2019-05-17 | 2020-11-17 | 中兴通讯股份有限公司 | 一种回波时延估计方法及装置 |
CN111198358B (zh) * | 2020-01-11 | 2022-06-03 | 桂林电子科技大学 | 一种基于相关检测的多旋翼无人机回波积累方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104239894A (zh) * | 2014-09-03 | 2014-12-24 | 西安电子科技大学 | 基于时域相关性特征的飞机目标分类方法 |
CN104803008A (zh) * | 2015-04-27 | 2015-07-29 | 中国直升机设计研究所 | 一种直升机模型旋翼旋转信号采集装置 |
CN104965215A (zh) * | 2014-02-18 | 2015-10-07 | 波音公司 | 监测辐射的系统和方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6870502B1 (en) * | 2003-08-29 | 2005-03-22 | Raytheon Company | Advanced asynchronous pulse detector |
-
2015
- 2015-11-25 CN CN201510831078.2A patent/CN105445716B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104965215A (zh) * | 2014-02-18 | 2015-10-07 | 波音公司 | 监测辐射的系统和方法 |
CN104239894A (zh) * | 2014-09-03 | 2014-12-24 | 西安电子科技大学 | 基于时域相关性特征的飞机目标分类方法 |
CN104803008A (zh) * | 2015-04-27 | 2015-07-29 | 中国直升机设计研究所 | 一种直升机模型旋翼旋转信号采集装置 |
Non-Patent Citations (1)
Title |
---|
飞机目标雷达回波微多普勒特性分析及应用研究;廖欢;《中国优秀硕士学位论文全文数据库 信息科技辑》;20130715(第7期);正文第8,18,21,41,44,45,52,62页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105445716A (zh) | 2016-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106597429A (zh) | 一种基于雷达信号处理系统的速度跟踪子系统 | |
CN102680948B (zh) | 一种线性调频信号调频率和起始频率估计方法 | |
CN106814353A (zh) | 一种雷达信号处理系统 | |
CN105445716B (zh) | 一种基于旋翼回波信号时域自相关的悬停直升机检测方法 | |
CN105429719B (zh) | 基于功率谱和多尺度小波变换分析强干扰信号检测方法 | |
CN106814358A (zh) | 一种用于超宽带穿墙雷达中多人体目标检测的系统及方法 | |
CN102890272A (zh) | 一种毫米波云雷达信号处理方法 | |
CN106443178A (zh) | 一种基于IQuinn‑Rife综合的正弦信号频率估计方法 | |
CN109617631B (zh) | 基于数字信道化瞬时参数测量的侦察系统自适应接收方法 | |
CN111693944A (zh) | 雷达有源干扰信号参数提取与干扰样式识别方法及装置 | |
CN109387813A (zh) | 一种基于宽带雷达信号接收的高精度频率测量方法 | |
CN110133632B (zh) | 一种基于cwd时频分析的复合调制信号识别方法 | |
CN109587089B (zh) | 一种提升无人机信号识别准确度的方法 | |
CN104991239B (zh) | 一种基于脉冲压缩雷达的距离旁瓣抑制方法 | |
CN106249209B (zh) | 一种抗速度欺骗干扰的自适应迭代估计方法 | |
CN103888389B (zh) | 一种时频重叠信号的幅度估计方法 | |
CN106526566A (zh) | 基于fpga高速预处理的信号脉内特征实时分析处理方法 | |
CN108333568A (zh) | 冲击噪声环境下基于Sigmoid变换的宽带回波Doppler和时延估计方法 | |
CN114578296B (zh) | 基于相位编码信号的间歇采样干扰抑制方法 | |
CN111983567A (zh) | 一种基于数字滤波的雷达信号微特征分析系统 | |
CN113259027B (zh) | 一种基于Haar变换的敌我识别信号码速率计算方法 | |
Tian et al. | Waveform design and processing in OFDM based radar-communication integrated systems | |
CN105137404A (zh) | 基于预脉冲间处理的雷达压缩采样方法及系统 | |
CN115390020A (zh) | Mcpc信号间歇采样干扰抑制方法、装置及雷达设备 | |
Liu et al. | Countermeasure for interrupted-sampling repeater jamming based on fractional Fourier transformation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |