CN103217670B - 一种基于pca的外辐射源微弱信号检测方法 - Google Patents
一种基于pca的外辐射源微弱信号检测方法 Download PDFInfo
- Publication number
- CN103217670B CN103217670B CN201310107865.3A CN201310107865A CN103217670B CN 103217670 B CN103217670 B CN 103217670B CN 201310107865 A CN201310107865 A CN 201310107865A CN 103217670 B CN103217670 B CN 103217670B
- Authority
- CN
- China
- Prior art keywords
- signal
- target
- time delay
- road signal
- path
- 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.)
- Expired - Fee Related
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种基于主量分析PCA的外辐射源微弱信号检测方法,包括以下步骤:求采样序列x(t)的自相关矩阵R=E{x(t)xH(t)};对自相关矩阵R进行奇异值分解SVD得到R=QΛQH,利用分解得到的对角阵Λ以及酉矩阵Q对采样序列x(t)进行线性变换得到接收信号从强到弱的排列y(t);将接收信号y(t)中第3路至第M路信号分别与第1路信号y1(t)做时频分布二维相关得到对应相关值|Ψi(τ,f)|,当有某一路信噪比大于等于门限值时,则将其峰值点对应的时延与多普勒频确定为一个目标的时延与多普勒频移,否则判定该路信号中不存在目标。本发明能够实现在强同频干扰以及低信噪比条件下提取微弱目标信号。
Description
技术领域
本发明属于信号处理技术,具体涉及微弱信号检测技术。
背景技术
基于外辐射源的双多基地雷达***能够利用地理位置已知的民用电视信号、广播信号、通信信号、卫星信号,以及手机基站等作为照射源,通过对空中运动目标反射的回波信号的检测,实现对目标的检测与定位跟踪。由于外辐射源雷达自身不发射信号,整个雷达***具有隐蔽性强、抗干扰能力强和生存能力强等优点,因而受到各国的广泛关注。
***的原理示意如图1所示,以广播信号为例,雷达到接收信号的包括有来自于广播发射塔本身的直达波、运动目标的回波、其它同频广播发射塔发送的直达波以及其它物体反射的回波,比如各种建筑反射的静物回波。因此,雷达接收信号中的目标回波信号是极其微弱的,并且淹没于各种干扰中。
在实际应用上它面临着许多的挑战,其中,最直接的一个问题就是基于民用信号的外辐射源雷达***如何更有效地实现对微弱目标的检测与跟踪,主要体现为在外辐射源雷达的接收信号中,目标回波信号淹没在直达波和多径信号等杂波以及接收机热噪声中,增加了目标检测的难度。
传统对外辐射源目标信号的检测方法是用参考通道接收直达波信号,主天线接收多路信号,其中包括直达波信号、目标回波信号以及噪声和干扰。处理过程如图2所示,首先用主通道信号和参考通道信号采用NLMS算法对直达波进行自适应对消,再通过空域滤波对同频干扰进行抑制,最后采用时频二维相关提取目标时延、多普勒频率,从而估计目标DOA(目标到达角度)参数。然而,此种方法在目标信号十分微弱的情况下将失效,同时,由于同频干扰也可能存在于同一个波束以内,这将导致无法提取目标时延、多普勒频率参数,从而无法估计目标的DOA。
发明内容
本发明所要解决的技术问题是,提供一种存在同频干扰情况下有效实现微弱目标回波信号检测的方法。
本发明为解决上述技术问题所采用的技术方案是,一种基于主量分析PCA的外辐射源微弱信号检测方法,包括以下步骤:
步骤1:接收方对接收到的待处理的主通道阵列数字调制信号进行带通采样得到采样序列x(t),t表示采样时刻;
步骤2:求采样序列x(t)的自相关矩阵R=E{x(t)xH(t)},E{·}表示求均值;
步骤3:对自相关矩阵R进行奇异值分解SVD得到R=QΛQH,利用分解得到的对角阵Λ以及酉矩阵Q对采样序列x(t)进行线性变换得到接收信号从强到弱的排列y(t),其中y(t)=[y1(t) y2(t) … yM(t)]T,M为接收天线阵元个数;
步骤4:将接收信号y(t)中第3路至第M路信号分别与第1路信号y1(t)做时频分布二维相关得到对应相关值|Ψi(τ,f)|,i=3,…,M,L为带通采样长度,τ为时延、f多普勒频移;
步骤5:分别在3至M各路信号对应的相关值|Ψi(τ,f)|中搜素峰值点,将各路信号峰值点对应的时延与多普勒频移作为一个有效信号,计算各路信噪比,当有某一路信噪比大于等于门限值时,则将其峰值点对应的时延与多普勒频确定为一个目标的时延与多普勒频移,否则判定该路信号中不存在目标。
本发明的有益效果是,本发明能够实现在强同频干扰以及低信噪比条件下提取微弱目标信号。另外本发明在在时频二维相干匹配提取时延、多普勒频移时,不需要再采用参考天线接收信号作为参考,从而起到了节约硬件成本的作用,另外本发明不用采用自适应对消等步骤,整体运算量也较传统方法有较大的减小。
附图说明
图1:基于广播外辐射源定位的***原理示意图。
图2:传统方法信号处理流程图。
图3:本发明的信号处理流程图。
图4:目标衰减-60dB,干扰衰减-10dB的情况下,传统方法目标信号的时频二维相关图。
图5:目标衰减-60dB,干扰衰减-10dB的情况下,本发明目标信号的时频二维相关图。
图6:目标衰减-60dB,干扰衰减-10dB的情况下,传统方法目标信号的时频二维相关图时延维。
图7:目标衰减-60dB,干扰衰减-10dB的情况下,本发明目标信号的时频二维相关图时延维。
图8:目标衰减-60dB,干扰衰减-10dB的情况下,传统方法目标信号的时频二维相关图频率维。
图9:目标衰减-60dB,干扰衰减-10dB的情况下,本发明目标信号的时频二维相关图频率维。
具体实施方式
本发明处理信号模型如下:
对阵列天线接收到的数字调制信号y(t)进行建模,主要针对带通信号进行处理,设s(t)表示某一个广播电视塔发射的调频广播连续的经过调制后的信号,则为目标的回波信号,其中M为接收天线阵元个数,a1,…,aM表示目标回波的衰减,τ1,…τM表示目标的时延,f1,…,fM表示目标的多普勒频移,则由阵列天线接收到的信号为 z(t)表示同频干扰信号,n(t)表示接收机的热噪声,则本发明所要做的工作是从接收到的信号中提取目标的时延、多普勒频率参数,这是为了在后续处理中能够估计目标DOA,从而最终实现对目标的定位与跟踪做准备。
实施步骤如图3所示:
步骤1:接收方对接收到的待处理的主通道阵列数字调制信号进行带通采样,得到采样序列x(t)=[x1(t) x2(t) … xM(t)]T,其中M表示阵元数,t表示采样时刻,采样长度为L,T表示矩阵转置。x(t)包含有目标时延、多普勒频率以及DOA的信息。
步骤2:求采样序列x(t)的自相关矩阵R=E{x(t)xH(t)},E{·}表示求均值。可假设采样序列在采样长度L内是平稳的随机过程,这样求统计平均可以采用求时间平均代替,则 H表示共轭转置。
步骤3:对自相关矩阵R进行奇异值分解SVD得到R=QΛQH,利用分解得到的对角阵Λ以及酉矩阵Q,对x(t)进行线性变换得到其中y(t)=[y1(t) y2(t) … yM(t)]T,y(t)是将接收信号从强到弱的排列,y1(t)为最强的直达波信号,即参考信号。
步骤4:采用时频二维相干匹配滤波提取目标的时延和多普勒频移:由于直达波信号比较强,故在进行主分量分析后所得到的第一路信号y1(t)可作为直达波信号,y1(t)=Aexp{j[2πfc(t-D1)+θ(t-D1)]}+n(t),其中n(t)为干扰项,fc为广播信号的载波频率,D1为广播信号从发射到接收所经过的时间,θ(t-D1)为音频信号调制项。而由于干扰功率仅次于直达波功率,故第二路信号y2(t)可看作是同频干扰信号。yi(t),i≠1,2为包含目标回波参数的信号,由于经过主分量分析,其中所含干扰和直达波分量已经比较小,故其表达式可写作yi(t)=Bexp{j[2π(fc+fdm)(t-Dm)+θ(t-Dm)]}+m(t),其中fdm为目标多普勒频率,Dm为广播信号从发射、经过目标的散射到达雷达接收天线所经过的时间,m(t)可看作与估计所需(τ,fd)参数无关的干扰项,对y1(t)做时延可得:
y1(t-n△t)=Aexp{j[2πfc(t-D1-n△t)+θ(t-D1-n△t)]}+n((t-n△t))
延时长度△t用采样周期,或者采样周期的整数倍。用y1(t-n△t)的共轭式与yi(t),i≠1,2相乘可得:
*表示共轭矩阵,可以看到,通过调整直达波信号的时延,当D2=D1+n△t时,上式中指数项中的两项相位将相消,因而上式将变为多普勒频移项和噪声项,即:
对上式做FFT运算,即可得到目标的时延和多普勒频移。
将步骤3得到的M路信号的第3路到第M路信号与第1路信号y1(t)做时频二维相关,i=3,…,M,(τ,f)为第i路信号对应的时延和多普勒频率。当存在目标时,第i路上的|Ψi(τ,f)|会产生一个峰值点|Ψi(τd,fd)|,搜索到峰值点后,将(τd,fd)作为一个有效信号,计算信噪比, 表示第i路信号上除|Ψi(τd,fd)|外,所有|Ψi(τ,f)|的均值,若其值大于预先设置的门限值13dB,则表示目标存在,对应的(τd,fd)即为一个目标的时延和多普勒频率;若其值小于这一门限值,则表示第i路信号上不存在目标。
为了验证本发明效果,分别采用两段独立的随机音频信号作为调制信号用以产生调频广播信号,以此作为仿真场景中的直达波信号和同频干扰,载波频率为80kHz,最大频偏为±75kHz,采样频率为200kHz。仿真平台操作***是WIN7,采用MATLAB进行编译。
表1给出混合信号的各仿真参数,并比较在不同条件下传统方法和本发明提出方法的对目标信号的检测性能,最后还比较了传统方法和本发明所提方法的运算量。
表1仿真信号参数
信号 | 时延(μs) | 功率衰减(dB) | 到达角(°) | 多普勒频率(Hz) |
主台直达波s1 | 0 | 0 | -50 | 0 |
主台多径s2 | 75 | -10 | -20 | 0 |
主台多径s3 | 40 | -20 | 30 | 0 |
临台干扰s4 | 0 | -10~-30 | 8 | 0 |
目标回波s5 | 439.12 | -60~-85 | 5 | 119.37 |
接收阵列阵元数为12。附加的噪声为-50dB,高于目标回波信号10dB~35dB。时频二维相干积累时间为1s。数字波束形成DBF采用切比雪夫加权,波束范围为从-45°到45°,波束间隔为10°。恒虚警检测器采用单元选大恒虚警检测器,设置恒虚警门限系数为4.5。对传统方法性能如表2,本发明提出方法性能如表3:
表2传统方法检测性能表
表3本发明所提方法检测性能表
设目标时延、多普勒频率处时频二维相关最大峰值为A,除通过恒虚警率处理CFAR可检测到的目标峰值以外,时频二维相关图上底噪的平均为B,则上表中最大峰值区分度定义为10log(A/B)。
通过表2和表4所提供的仿真数据可以说明本发明所提出的方法在有强同频干扰,目标信号和干扰在同一个波束以内,以及目标信号十分微弱的情况下都能有效的检测到目标信号,并提取出目标的时延、多普勒频率参数。
目标衰减-60dB,干扰衰减-10dB的条件下,传统方法和本发明的目标信号时频二维相关图分别为图4、图5;传统方法和本发明的目标信号时频二维相关图时延维分别为图6、图7;传统方法和本发明的目标信号时频二维相关图频率维分别为图8、图9。可以分别从,传统方法和本发明的时频二维相关图时延维和多普勒维观察对比目标,通过对比可以发现,传统方法时频二维相关图的底噪很高,目标的最大峰值区分度较低,而本发明时频二维相关图的底噪较小,且目标最大峰值区分度较高,检测效果好。
由于传统方法和本发明都要用到时频二维相关,故比较计算量只需比较时频二维相关前的运算量,设接收阵列阵元个数为M,快拍数为N,波束个数为L,则比较结果如表4、表5所示:
表4传统方法运算量表
表5本发明所提方法运算量表
通过以上计算对比可以看出,本发明比传统方法的运算量有一定量的减小。且本方法不需要使用参考天线对准直达波信号,在一定程度上节省了硬件的成本。
Claims (5)
1.一种基于PCA的外辐射源微弱信号检测方法,包括以下步骤:
步骤1:接收方对接收到的待处理的主通道阵列数字调制信号进行带通采样得到采样序列x(t),t表示采样时刻;
步骤2:求采样序列x(t)的自相关矩阵R;
步骤3:对自相关矩阵R进行奇异值分解SVD得到R=QΛQH,利用分解得到的对角阵Λ以及酉矩阵Q对采样序列x(t)进行线性变换得到接收信号从强到弱的排列y(t),其中y(t)=[y1(t) y2(t)…yM(t)]T,M为接收天线阵元个数;
步骤4:将接收信号y(t)中第3路至第M路信号分别与第1路信号y1(t)做时频分布二维相关得到对应相关值|Ψi(τ,f)|, L为带通采样长度,τ为时延、f多普勒频移;
步骤5:分别在3至M各路信号对应的二维相关值|Ψi(τ,f)|中搜素峰值点,将各路信号峰值点对应的时延与多普勒频移作为一个有效信号,计算各路信噪比,当有某一路信噪比大于等于门限值时,则将其峰值点对应的时延与多普勒频移确定为一个目标的时延与多普勒频移,否则判定该路信号中不存在目标。
2.如权利要求1所述一种基于PCA的外辐射源微弱信号检测方法,其特征在于,自相关矩阵R=E{x(t)xH(t)},E{·}表示求均值。
3.如权利要求1所述一种基于PCA的外辐射源微弱信号检测方法,其特征在于, H表示共轭转置。
4.如权利要求1、2或3所述一种基于PCA的外辐射源微弱信号检测方法,其特征在于,步骤5中计算各路信噪比的方法为:|Ψi(τd,fd)|为第i路上的二维相关值|Ψi(τ,f)|的峰值点,表示第i路信号上除|Ψi(τd,fd)|外,所有|Ψi(τ,f)|的均值。
5.如权利要求4所述一种基于PCA的外辐射源微弱信号检测方法,其特征在于,所述门限值为13dB。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310107865.3A CN103217670B (zh) | 2013-03-29 | 2013-03-29 | 一种基于pca的外辐射源微弱信号检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310107865.3A CN103217670B (zh) | 2013-03-29 | 2013-03-29 | 一种基于pca的外辐射源微弱信号检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103217670A CN103217670A (zh) | 2013-07-24 |
CN103217670B true CN103217670B (zh) | 2015-01-21 |
Family
ID=48815626
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310107865.3A Expired - Fee Related CN103217670B (zh) | 2013-03-29 | 2013-03-29 | 一种基于pca的外辐射源微弱信号检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103217670B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412314B (zh) * | 2013-08-01 | 2015-07-29 | 清华大学 | 一种基于天体辐射信号相干性的成像方法 |
CN106559371B (zh) | 2015-09-24 | 2021-01-12 | 索尼公司 | 用于无线通信的电子设备以及无线通信方法 |
CN105866750B (zh) * | 2016-05-23 | 2018-10-12 | 西安电子科技大学 | 一种多个gps卫星微弱回波信号检测的方法 |
CN106597409B (zh) * | 2016-12-19 | 2019-02-12 | 电子科技大学 | 一种基于外辐射源雷达的时延-多普勒频移提取方法 |
CN108287341A (zh) * | 2017-12-15 | 2018-07-17 | 空间防务技术(深圳)有限公司 | 一种低空目标探测方法及其***、电子设备 |
CN109768942B (zh) * | 2018-12-20 | 2021-05-14 | 合肥工业大学 | 一种瑞利信道最大多普勒频移估计方法 |
CN110609263B (zh) * | 2019-10-29 | 2022-11-04 | 电子科技大学 | 一种同时计算脉冲激光雷达目标回波时延和频偏的方法 |
CN113341393B (zh) * | 2021-06-10 | 2023-11-10 | 西安空间无线电技术研究所 | 一种无直达波通道外辐射源雷达微弱目标检测方法 |
CN116840775B (zh) * | 2023-06-30 | 2024-02-27 | 中国人民解放军军事科学院***工程研究院 | 一种已知强信号场景的弱信号提取方法 |
CN117741586B (zh) * | 2024-02-19 | 2024-05-07 | 中国科学院空天信息创新研究院 | 一种超宽带自适应阵列接收装置及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1653354A (zh) * | 2002-03-13 | 2005-08-10 | 雷神加拿大有限公司 | 雷达中用于频谱生成的***和方法 |
CN101741348A (zh) * | 2009-12-09 | 2010-06-16 | 北京天碁科技有限公司 | 多相滤波器、数字信号处理***和滤波方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6271786B1 (en) * | 1999-06-07 | 2001-08-07 | Raytheon Company | Random noise radar target detection device |
-
2013
- 2013-03-29 CN CN201310107865.3A patent/CN103217670B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1653354A (zh) * | 2002-03-13 | 2005-08-10 | 雷神加拿大有限公司 | 雷达中用于频谱生成的***和方法 |
CN101741348A (zh) * | 2009-12-09 | 2010-06-16 | 北京天碁科技有限公司 | 多相滤波器、数字信号处理***和滤波方法 |
Non-Patent Citations (2)
Title |
---|
基于FPGA的MLS信号处理;王宇伟等;《制导与引信》;20070630;第28卷(第2期);36-41 * |
无源定位中的弱信号检测与提取技术研究;龙宁;《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》;20051115;1-69 * |
Also Published As
Publication number | Publication date |
---|---|
CN103217670A (zh) | 2013-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103217670B (zh) | 一种基于pca的外辐射源微弱信号检测方法 | |
CN103293518B (zh) | 广播信号外辐射源定位检测方法 | |
Liu et al. | On the performance of the cross-correlation detector for passive radar applications | |
CN106125053B (zh) | 脉冲多普勒雷达极化抗干扰方法 | |
CN106909779B (zh) | 基于分布式处理的mimo雷达克拉美罗界计算方法 | |
CN103235294A (zh) | 一种基于外辐射源定位的微弱信号分离估计方法 | |
EP1972962A2 (en) | Transmitter independent techniques to extend the performance of passive coherent location | |
CN104007421B (zh) | 基于全变差和压缩感知的罗兰c被动雷达toa估计方法 | |
CN107561507B (zh) | 一种外辐射源雷达的杂波对消方法 | |
CN106872968A (zh) | 一种基于ofdm信号的外辐射源雷达微弱目标检测方法 | |
CN103018727A (zh) | 一种基于样本训练的机载雷达非平稳杂波抑制方法 | |
CN104765020A (zh) | 有源假目标干扰的极化鉴别方法 | |
CN102087354A (zh) | 无源雷达分组ls-clean微弱目标检测方法 | |
CN111948618B (zh) | 一种基于卫星外辐射源的前向散射目标探测方法及*** | |
CN106054165A (zh) | 一种实现分布式无源雷达目标检测的方法 | |
CN102353947A (zh) | 一种基于csa-mwf的无源雷达目标回波信号子空间的估计方法 | |
CN105891817A (zh) | 一种无直达波条件下分布式无源雷达目标检测的方法 | |
Bongioanni et al. | Performance analysis of a multi-frequency FM based passive bistatic radar | |
CN104142496A (zh) | 基于连通域划分的统计mimo雷达多目标定位方法 | |
CN107272030B (zh) | 一种基于多个gps卫星信息融合的目标探测方法 | |
Zhang et al. | Interrupted sampling repeater jamming recognition and suppression based on phase-coded signal processing | |
CN106597381A (zh) | 全相参全极化mimo雷达四通道融合目标检测方法 | |
Hamdollahzadeh et al. | Moving target localization in bistatic forward scatter radars: Performance study and efficient estimators | |
CN104931946A (zh) | 外辐射源雷达多普勒-时延两维双门限目标检测方法 | |
CN111007487A (zh) | 一种基于时间反演的多基地雷达目标检测方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150121 Termination date: 20190329 |