CN105866760A - 一种动目标时空二维特性分析方法 - Google Patents

一种动目标时空二维特性分析方法 Download PDF

Info

Publication number
CN105866760A
CN105866760A CN201610333814.6A CN201610333814A CN105866760A CN 105866760 A CN105866760 A CN 105866760A CN 201610333814 A CN201610333814 A CN 201610333814A CN 105866760 A CN105866760 A CN 105866760A
Authority
CN
China
Prior art keywords
pseudo
ambiguity function
signal
noise
peak value
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
CN201610333814.6A
Other languages
English (en)
Other versions
CN105866760B (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.)
Nanjing Changyuan Intelligent Technology Co Ltd
Original Assignee
Nanjing Changyuan Intelligent Technology Co Ltd
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 Nanjing Changyuan Intelligent Technology Co Ltd filed Critical Nanjing Changyuan Intelligent Technology Co Ltd
Priority to CN201610333814.6A priority Critical patent/CN105866760B/zh
Publication of CN105866760A publication Critical patent/CN105866760A/zh
Application granted granted Critical
Publication of CN105866760B publication Critical patent/CN105866760B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details 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

Abstract

本发明公开了一种动目标时空二维特性分析方法,步骤如下:第一步,对截获得到的中频采样信号x(t)作模糊函数变换,得到模糊函数分布矩阵第二步,搜索的最大值χmax和次最大值χsub max,利用模糊函数与噪声的信噪比关系SNR=10lg[χsub max/(χmax‑χsub max)],在的中心位置去除噪声,得到连续波伪码调相的模糊函数第三步,取切面在此平面内搜索各个峰值所对应的延时τi,计算它们之间的差值,然后求平均,即得伪码调制信号的周期第四步,取切面其中n=1,2,…2N+1,此时在各个切面上求取峰值之间的差值,通过求这些差值的平均值即可得到码元速率1/Tc。该动目标时空二维特性分析方法为基于模糊函数特征参数提取法,在信噪比为dB以上时都有较高的估计精度,整个识别过程不含复杂变换,易于通过软件来实现。

Description

一种动目标时空二维特性分析方法
技术领域
本发明涉及动目标时空二维特性分析方法技术领域,具体为一种动目标时空二维特性分析方法。
背景技术
前已有不少文献提出各种对伪码调相信号特征参数提取方法,取得较好的结果,然而这些识别方法大都是采用谱相关技术或者Cohen类时频分析技术,算法复杂度高、运算量大,不适合应用于近程雷达的实时侦察识别。目标的所有特性都体现在回波信号的特征参数上,通过测量回波信号与发射信号的频率和时间差值等参数,能够确定目标的速度、目标与信源间的距离等。因此对信号进行参数提取是目标识别的前提,具有重要意义。周期模糊函数是对低截获概率信号进行分析研究的有效工具,它描述了信号的全部特性,不同形式的信号具有不同的周期模糊函数。伪码调相连续波近程雷达是利用伪随机m序列的自相关函数类似于白噪声相关函数的随机特性而工作的。该体制雷达抗干扰能力强,具有良好的距离分辨率和距离截止特性,现已成为近程探测雷达的重要发展方向之一,因此研究伪码调相连续波信号的识别和特征参数提取具有重要的意义
发明内容
本发明的目的在于提供一种动目标时空二维特性分析方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:一种动目标时空二维特性分析方法,在时间延时τ轴上,M个周期的连续波伪码调相信号的周期模糊函数在延时轴为伪码周期的整数倍上产生峰值,且在ξ=0的时延轴上,产生的峰值个数为(2M+1),通过这些峰值之间的差值来估计伪码调相信号的周期Tp; 在多普勒频移轴ξ上,周期模糊函数在码元速率1/Tc的整数倍上产生峰值,在对应的τ为伪码调相信号周期整数倍时都会有这样一组峰值出现,这些峰值的多普勒频移值可以求出伪码的码元宽度Tc;通过在截获信号的模糊函数中心位置去除噪声后,即得到连续波伪码调相的模糊函数;对基于续波伪码调相信号的动目标特征参数提取步骤如下:
第一步,对截获得到的中频采样信号x(t)作模糊函数变换,得到模糊函数分布矩阵
第二步,搜索的最大值χmax和次最大值χsub max,利用模糊函数与噪声的信噪比关系SNR=10lg[χsub max/(χmaxsub max)],在的中心位置去除噪声,得到连续波伪码调相的模糊函数
第三步,取切面在此平面内搜索各个峰值所对应的延时τl,计算它们之间的差值,然后求平均,即得伪码调制信号的周期
第四步,取切面其中n=1,2,…2N+1,此时在各个切面上求取峰值之间的差值,通过求这些差值的平均值即可得到码元速率1/Tc
与现有技术相比,本发明的有益效果是:该动目标时空二维特性分析方法为基于模糊函数特征参数提取法,在信噪比为dB以上时都有较高的估计精度。相比较于传统的参数提取算法,这些特征参数提取算法复杂度降低了很多,大大降低了计算量,体现了近程雷达侦察干扰实时性的要求;且整个识别过程不含复杂变换,易于通过软件来实现。
附图说明
图1为本发明对6个周期连续波伪码调相信号参数估计归一化均方误差分析示意图。
图2为本发明对30个周期连续波伪码调相信号参数估计归一化均方误差分析示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供一种技术方案:
实施例1:
一种动目标时空二维特性分析方法,在时间延时τ轴上,M个周期的连续波伪码调相信号的周期模糊函数在延时轴为伪码周期的整数倍上产生峰值,且在ξ=0的时延轴上,产生的峰值个数为(2M+1),通过这些峰值之间的差值来估计伪码调相信号的周期Tp;在多普勒频移轴ξ上,周期模糊函数在码元速率1/Tc的整数倍上产生峰值,在对应的τ为伪码调相信号周期整数倍时都会有这样一组峰值出现,这些峰值的多普勒频移值可以求出伪码的码元宽度Tc;通过在截获信号的模糊函数中心位置去除噪声后,即得到连续波伪码调相的模糊函数;对基于续波伪码调相信号的动目标特征参数提取步骤如下:
第一步,对截获得到的中频采样信号x(t)作模糊函数变换,得到模糊函数分布矩阵
第二步,搜索的最大值χmax和次最大值χsub max,利用模糊函数与噪声的信噪比关系SNR=10lg[χsub max/(χmaxsub max)],在的中心位置去除噪声,得到连续波伪码调相的模糊函数
第三步,取切面在此平面内搜索各个峰值所对应的延时τl,计算它们之间的差值,然后求平均,即得伪码调制信号的周期
第四步,取切面其中n=1,2,…2N+1,此时在各个切面上求取峰值之间的差值,通过求这些差值的平均值即可得到码元速率1/Tc
实施例2:
根据实施例1所述的一种动目标时空二维特性分析方法,从图中可以看出,在信噪比为0dB以上都有较高的估计精度。显然,码元宽度和伪码周期的估计误差随着信噪比的增加而减小,而且截取6个伪码周期的信号进行参数估计时,码元宽度的估计方差在信噪比为5dB时才趋于稳定,而截取30个伪码周期信号进行参数估计时,码元宽度的估计方差在0dB时基本趋于稳定;同样的,伪码周期的估计方差在6个伪码周期时在0dB趋于稳定,而在30个伪码周期时在-5dB趋于稳定。显然,随着采样数据的增加,噪声对参数估计精度的影响越来越小。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。

Claims (1)

1.一种动目标时空二维特性分析方法,其特征在于:在时间延时τ轴上,M个周期的连续波伪码调相信号的周期模糊函数在延时轴为伪码周期的整数倍上产生峰值,且在ξ=0的时延轴上,产生的峰值个数为(2M+1),通过这些峰值之间的差值来估计伪码调相信号的周期Tp;在多普勒频移轴ξ上,周期模糊函数在码元速率1/Tc的整数倍上产生峰值,在对应的τ为伪码调相信号周期整数倍时都会有这样一组峰值出现,这些峰值的多普勒频移值可以求出伪码的码元宽度Tc;通过在截获信号的模糊函数中心位置去除噪声后,即得到连续波伪码调相的模糊函数;对基于续波伪码调相信号的动目标特征参数提取步骤如下:
第一步,对截获得到的中频采样信号x(t)作模糊函数变换,得到模糊函数分布矩阵
第二步,搜索的最大值χmax和次最大值χsub max,利用模糊函数与噪声的信噪比关系SNR=10lg[χsub max/(χmaxsub max)],在的中心位置去除噪声,得到连续波伪码调相的模糊函数
第三步,取切面在此平面内搜索各个峰值所对应的延时τi,计算它们之间的差值,然后求平均,即得伪码调制信号的周期
第四步,取切面其中n=1,2,…2N+1,此时在各个切面上求取峰值之间的差值,通过求这些差值的平均值即可得到码元速率1/Tc
CN201610333814.6A 2016-05-20 2016-05-20 一种动目标时空二维特性分析方法 Active CN105866760B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610333814.6A CN105866760B (zh) 2016-05-20 2016-05-20 一种动目标时空二维特性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610333814.6A CN105866760B (zh) 2016-05-20 2016-05-20 一种动目标时空二维特性分析方法

Publications (2)

Publication Number Publication Date
CN105866760A true CN105866760A (zh) 2016-08-17
CN105866760B CN105866760B (zh) 2018-09-11

Family

ID=56635276

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610333814.6A Active CN105866760B (zh) 2016-05-20 2016-05-20 一种动目标时空二维特性分析方法

Country Status (1)

Country Link
CN (1) CN105866760B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3811283A1 (de) * 1988-04-02 1989-10-12 Licentia Gmbh Pulsdopplerradar
CN101173985A (zh) * 2006-11-01 2008-05-07 中国科学院国家天文台 利用卫星信号探测低空目标的被动雷达探测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3811283A1 (de) * 1988-04-02 1989-10-12 Licentia Gmbh Pulsdopplerradar
CN101173985A (zh) * 2006-11-01 2008-05-07 中国科学院国家天文台 利用卫星信号探测低空目标的被动雷达探测方法

Also Published As

Publication number Publication date
CN105866760B (zh) 2018-09-11

Similar Documents

Publication Publication Date Title
CN107561508B (zh) 一种用于匀加速运动目标的相参积累检测方法
CN104020469B (zh) 一种mimo雷达距离-角度二维超分辨率成像算法
CN106125053B (zh) 脉冲多普勒雷达极化抗干扰方法
CN102628937B (zh) 基于广义keystone变换和非相参积累的雷达检测方法
CN102230961B (zh) 基于相位补偿处理的扩展目标检测方法
CN107015205B (zh) 一种分布式mimo雷达检测的虚假目标消除方法
CN106383340A (zh) 一种随机脉冲初始相位雷达的速度假目标识别方法
CN103018727A (zh) 一种基于样本训练的机载雷达非平稳杂波抑制方法
CN110133632B (zh) 一种基于cwd时频分析的复合调制信号识别方法
CN104502898B (zh) 将修正rft和修正mdcft相结合的机动目标参数估计方法
CN103439692A (zh) 基于协方差矩阵广对称特性的stap方法
CN113221631B (zh) 基于卷积神经网络的序列脉冲抗干扰目标检测方法
CN109581318A (zh) 基于时间反转非均匀采样的雷达高机动目标相参积累检测方法
CN109799488A (zh) 一种非参数搜索的雷达机动目标长时间相参积累方法
CN109088838A (zh) 一种高动态下直扩dpsk信号的伪码—多普勒快速捕获方法
CN107255814A (zh) 一种基于lfmsk波形的雷达目标检测方法
CN109324322A (zh) 一种基于被动相控阵天线的测向与目标识别方法
CN104991239B (zh) 一种基于脉冲压缩雷达的距离旁瓣抑制方法
CN104375121A (zh) 基于先验信息的mimo雷达波形与有偏估计器的联合优化方法
CN110398730A (zh) 基于坐标旋转和非均匀傅里叶变换机动目标相参检测方法
CN114755639B (zh) 一种面向跟踪对抗的多样式复合欺骗干扰自适应生成方法
CN109782249B (zh) 一种两目标相关时延估计算法
CN103344948A (zh) 一种利用稀疏傅里叶变换计算外辐射源雷达互模糊函数的方法
CN101436877A (zh) 一种抗多径干扰的pn码自适应门限捕获方法
CN105445767A (zh) 基于平均模糊函数的boc信号参数盲估计方法

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