CN103760531A - 一种窄带雷达的空中目标回波的噪声抑制方法 - Google Patents
一种窄带雷达的空中目标回波的噪声抑制方法 Download PDFInfo
- Publication number
- CN103760531A CN103760531A CN201410032365.2A CN201410032365A CN103760531A CN 103760531 A CN103760531 A CN 103760531A CN 201410032365 A CN201410032365 A CN 201410032365A CN 103760531 A CN103760531 A CN 103760531A
- Authority
- CN
- China
- Prior art keywords
- matrix
- noise
- signal
- aerial target
- narrow
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000002592 echocardiography Methods 0.000 title abstract description 5
- 230000000452 restraining effect Effects 0.000 title abstract 4
- 239000011159 matrix material Substances 0.000 claims abstract description 54
- 238000012545 processing Methods 0.000 claims abstract description 24
- 239000013598 vector Substances 0.000 claims abstract description 17
- 230000001629 suppression Effects 0.000 claims description 20
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 11
- 238000004088 simulation Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012549 training Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
Images
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/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
- G01S7/2928—Random or non-synchronous interference pulse cancellers
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
本发明属于窄带雷达空中目标回波的噪声抑制技术领域,公开了一种窄带雷达的空中目标回波的噪声抑制方法。该窄带雷达的空中目标回波的噪声抑制方法,包括以下步骤:窄带雷达接收空中目标的回波信号,对空中目标的回波信号进行滑窗处理,得出滑窗处理后矩阵;估计噪声信号在时域中的平均功率M;得出滑窗处理后矩阵E的自相关矩阵R;计算R的特征向量;用R的特征向量组成特征向量矩阵U;设置噪声门限thresh_α,用λr减去thresh_α;如果r取P时,λr-1-thresh_α>0,λr-thresh_α<0,则用U的前P列构成矩阵Us;根据矩阵Es中第一列元素和最后一行的元素,得出经噪声抑制处理后的空中目标的回波信号x。
Description
技术领域
本发明属于窄带雷达空中目标回波的噪声抑制技术领域,特别涉及一种窄带雷达的空中目标回波的噪声抑制方法。
背景技术
在现代战争中,由于作战任务和环境的需要,直升机、螺旋桨飞机和喷气式飞机承担着不同的任务,具有的威胁程度不同,因此实现对这三类飞机目标的探测具有非常重要的意义。
目前为止,现有的空中目标回波的噪声抑制方法中通常使用训练数据和测试数据,而训练数据和测试数据的信噪比条件相差不大。但在实际情况中,训练数据通常都是通过一些合作性的实验或者通过仿真来获取的高信噪比数据;然而,测试数据多数情况下都是在一些非合作的条件下获取的,例如在战争过程中,战场的复杂环境、气候条件以及目标距离雷达的远近都会造成测试数据中噪声强度的不同。这样,就会给空中目标回波的噪声抑制带来一定的困难。
发明内容
本发明的目的在于提出一种窄带雷达的空中目标回波的噪声抑制方法。该窄带雷达的空中目标回波的噪声抑制方法能够对空中目标回波进行有效地噪声抑制,可用于对窄带雷达空中目标回波进行预处理。
为实现上述技术目的,本发明采用如下技术方案予以实现。
一种窄带雷达的空中目标回波的噪声抑制方法,包括以下步骤:
S1:所述窄带雷达接收空中目标的回波信号,所述空中目标的回波信号的脉冲积累数为N,对所述空中目标的回波信号进行滑窗处理,窗长为得出滑窗处理后矩阵E,表示N/2向下取整;估计噪声信号在时域中的平均功率M;
S2:得出滑窗处理后矩阵E的自相关矩阵R,其中,EH表示E的复共轭转置矩阵;然后计算R的特征值和特征向量,R的特征值从大到小依次为λ1至用R的特征向量组成特征向量矩阵U,U的第r列为R的第r特征向量ur,r取1至ur为与λr对应的特征向量;
S3:设置噪声门限thresh_α,用λr减去thresh_α;如果r取P时,λr-1-thresh_α>0,λr-thresh_α<0,则用U的前P列构成矩阵Us;
本发明的特点和进一步改进在于:
在步骤S1中,所述窄带雷达接收到的空中目标的回波信号表示为矩阵S,S=[S1,...,SN]T,其中,Sk空中目标的回波信号在第k时刻的信号,k取1至N;滑窗处理后矩阵其中,为矩阵S进行第L次滑窗处理后得到的信号,L取1至
在步骤S1中,将所述空中目标的回波信号进行离散傅里叶变换,使所述空中目标的回波信号为频域信号,根据所述频域信号,估计噪声信号在频域中的平均功率;然后将所述噪声信号在频域中的平均功率除以离散傅里叶变换点数,得出噪声信号在时域中的平均功率M。
本发明的有益效果为:本发明的窄带雷达的空中目标回波的噪声抑制方法能够对空中目标回波进行有效地噪声抑制,可用于对窄带雷达空中目标回波进行预处理。
附图说明
图1为本发明实施例的一种窄带雷达的空中目标回波的噪声抑制方法的流程示意图;
图2a为喷气式飞机在不含噪声时的多普勒频谱示意图;
图2b为喷气式飞机在含噪声时的多普勒频谱示意图;
图2c为在使用本发明对噪声信号进行抑制后喷气式飞机的多普勒频谱示意图;
图3a为螺旋桨飞机在不含噪声时的多普勒频谱示意图;
图3b为螺旋桨飞机在含噪声时的多普勒频谱示意图;
图3c为在使用本发明对噪声信号进行抑制后螺旋桨飞机的多普勒频谱示意图;
图4a为直升机在不含噪声时的多普勒频谱示意图;
图4b为直升机在含噪声时的多普勒频谱示意图;
图4c为在使用本发明对噪声信号进行抑制后直升机的多普勒频谱示意图;
图5为不同门限条件下平均重构误差随信噪比的稳健性曲线示意图。
具体实施方式
下面结合附图对本发明作进一步说明:
参照图1,为本发明实施例的一种窄带雷达的空中目标回波的噪声抑制方法的流程示意图。该窄带雷达的空中目标回波的噪声抑制方法包括以下步骤:
S1:通过滑窗法处理窄带雷达接收的空中目标的回波信号,并估计噪声信号在时域中的平均功率。具体过程如下:
窄带雷达接收空中目标的回波信号,例如窄带雷达接收到的空中目标的回波信号为低信噪比信号。将窄带雷达接收到的空中目标的回波信号表示为矩阵S,S=[S1,...,SN]T,其中,Sk空中目标的回波信号在第k时刻的信号,k取1至N,N为空中目标的回波信号的脉冲积累数。
对上述空中目标的回波信号进行滑窗处理,得出滑窗处理后矩阵E;窗长为得出滑窗处理后矩阵E,表示N/2向下取整;滑窗处理后矩阵其中,为矩阵S进行第L次滑窗处理后得到的信号,L取1至[.]T表示矩阵的转置。
将上述空中目标的回波信号进行离散傅里叶变换,使上述空中目标的回波信号为频域信号,根据上述频域信号,估计对应的噪声信号在频域中的平均功率,上述噪声信号位于上述空中目标的回波信号中。将上述噪声信号在频域中的平均功率除以离散傅里叶变换点数,得出噪声信号在时域中的平均功率M,M=σ2。
S2:得出滑窗处理后矩阵的自相关矩阵,并得出该自相关矩阵的特征值和特征向量。其具体过程如下:
首先根据滑窗处理后矩阵E,得出滑窗处理后矩阵E的自相关矩阵R,其中,EH为E的复共轭转置矩阵。然后对自相关矩阵R作特征值分解,得出该自相关矩阵R的特征值和特征向量。自相关矩阵R的特征值从大到小依次为λ1至此时,也可以表示为c和d均取1至且c>d。
S3:设置噪声门限,并根据该噪声门限确定主分量的个数。其具体过程如下:
首先设置噪声门限thresh_α,噪声门限thresh_α可以采取多种方式,例如表示自由度为2的卡方分布在水平α上的上侧分位数(即χ2(2)在水平α上的上侧分位数)。下面说明将thresh_α设置为的理由:假设上述噪声信号的实部和虚部均服从零均值的高斯分布,w=w1+jw2,其中,w1为上述噪声信号的实部,w2为上述噪声信号的虚部。则 经推导后可得出:
由此可知,上述噪声信号服从自由度为2的χ2(2)分布,根据χ2(2)分布来设置噪声门限thresh_α较为合理。
然后,依次用λ1至减去thresh_α,直到得出一个小于thresh_α的特征值λp,即λp-1-thresh_α>0,λp-thresh_α<0,这样,P为主分量的个数(或信号子空间维数)。
然后,用U的前P列构成矩阵Us,Us=[u1,...,uP]。
S4:得出经噪声抑制处理后的空中目标的回波信号(也就是将信号投影到信号子空间),其具体过程如下:
下面通过仿真实验对本发明作进一步说明:
仿真实验中所用数据包含常规窄带雷达接收的直升机、螺旋桨飞机和喷气式飞机的含噪实测信号。在仿真实验中,使用本发明对上述三类飞机的回波信号进行噪声抑制,其中,在设置噪声门限时,α=0.10或α=0.01。在使用本发明对上述三类飞机的回波信号进行噪声抑制之后,计算所有样本的平均重构误差。参照图2a,为喷气式飞机在不含噪声时的多普勒频谱示意图,参照图2b,为喷气式飞机在含噪声时的多普勒频谱示意图,参照图2c,为在使用本发明对噪声信号进行抑制后喷气式飞机的多普勒频谱示意图;参照图3a,为螺旋桨飞机在不含噪声时的多普勒频谱示意图,参照图3b,为螺旋桨飞机在含噪声时的多普勒频谱示意图,参照图3c,为在使用本发明对噪声信号进行抑制后螺旋桨飞机的多普勒频谱示意图;参照图4a,为直升机在不含噪声时的多普勒频谱示意图,参照图4b,为直升机在含噪声时的多普勒频谱示意图,参照图4c,为在使用本发明对噪声信号进行抑制后直升机的多普勒频谱示意图。从上述图中,可以看出:通过本发明的窄带雷达的空中目标回波的噪声抑制方法,可以对上述三类飞机的回波信号中含有的噪声进行较好的抑制。
参照图5,为不同门限条件下平均重构误差随信噪比的稳健性曲线示意图。图5中,门限1对应的α为0.10,门限2对应的α为0.01。从图5中可以看出,利用本发明得到的经过噪声抑制处理后的信号更加接近于不含噪声信号,即平均重构误差更小。同时,从图5还可以看出,门限的选择对结果影响较大,因此,在使用本发明介绍的方法时,要选择合适的门限才能得到更加理想的结果。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (4)
1.一种窄带雷达的空中目标回波的噪声抑制方法,其特征在于,包括以下步骤:
S1:所述窄带雷达接收空中目标的回波信号,所述空中目标的回波信号的脉冲积累数为N,对所述空中目标的回波信号进行滑窗处理,窗长为得出滑窗处理后矩阵E,表示N/2向下取整;估计噪声信号在时域中的平均功率M;
S2:得出滑窗处理后矩阵E的自相关矩阵R,其中,EH表示E的复共轭转置矩阵;然后计算R的特征值和特征向量,R的特征值从大到小依次为λ1至用R的特征向量组成特征向量矩阵U,U的第r列为R的第r特征向量ur,r取1至ur为与λr对应的特征向量;
S3:设置噪声门限thresh_α,用λr减去thresh_α;如果r取P时,λr-1-thresh_α>0,λr-thresh_α<0,则用U的前P列构成矩阵Us;
3.如权利要求1所述的一种窄带雷达的空中目标回波的噪声抑制方法,其特征在于,在步骤S1中,将所述空中目标的回波信号进行离散傅里叶变换,使所述空中目标的回波信号为频域信号,根据所述频域信号,估计噪声信号在频域中的平均功率;然后将所述噪声信号在频域中的平均功率除以离散傅里叶变换点数,得出噪声信号在时域中的平均功率M。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410032365.2A CN103760531B (zh) | 2014-01-23 | 2014-01-23 | 一种窄带雷达的空中目标回波的噪声抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410032365.2A CN103760531B (zh) | 2014-01-23 | 2014-01-23 | 一种窄带雷达的空中目标回波的噪声抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103760531A true CN103760531A (zh) | 2014-04-30 |
CN103760531B CN103760531B (zh) | 2016-01-20 |
Family
ID=50527798
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410032365.2A Active CN103760531B (zh) | 2014-01-23 | 2014-01-23 | 一种窄带雷达的空中目标回波的噪声抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103760531B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109655794A (zh) * | 2018-11-14 | 2019-04-19 | 上海无线电设备研究所 | 一种窄带自卫噪声压制干扰的检测识别方法 |
CN115877342A (zh) * | 2023-01-09 | 2023-03-31 | 北京海兰信数据科技股份有限公司 | 一种目标对象的确定方法、装置及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003079045A2 (en) * | 2002-03-13 | 2003-09-25 | Raytheon Canada Limited | System and method for spectral generation in radar |
CN102608598A (zh) * | 2012-03-19 | 2012-07-25 | 西安电子科技大学 | 基于子空间投影的实孔径前视成像方法 |
CN103412287A (zh) * | 2013-09-01 | 2013-11-27 | 西安电子科技大学 | 基于lvd的线性调频信号参数估计方法 |
-
2014
- 2014-01-23 CN CN201410032365.2A patent/CN103760531B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003079045A2 (en) * | 2002-03-13 | 2003-09-25 | Raytheon Canada Limited | System and method for spectral generation in radar |
CN102608598A (zh) * | 2012-03-19 | 2012-07-25 | 西安电子科技大学 | 基于子空间投影的实孔径前视成像方法 |
CN103412287A (zh) * | 2013-09-01 | 2013-11-27 | 西安电子科技大学 | 基于lvd的线性调频信号参数估计方法 |
Non-Patent Citations (2)
Title |
---|
F. COLONE等: "Cancellation of clutter and multipath in passive radar using a sequential approach", 《RADAR, 2006 IEEE CONFERENCE ON》 * |
田光明等: "基于Wigner分布时频遮隔的信号分解算法", 《电子学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109655794A (zh) * | 2018-11-14 | 2019-04-19 | 上海无线电设备研究所 | 一种窄带自卫噪声压制干扰的检测识别方法 |
CN109655794B (zh) * | 2018-11-14 | 2020-09-08 | 上海无线电设备研究所 | 一种窄带自卫噪声压制干扰的检测识别方法 |
CN115877342A (zh) * | 2023-01-09 | 2023-03-31 | 北京海兰信数据科技股份有限公司 | 一种目标对象的确定方法、装置及设备 |
CN115877342B (zh) * | 2023-01-09 | 2023-08-18 | 江苏途索海洋技术服务有限公司 | 一种目标对象的确定方法、装置及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN103760531B (zh) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104569948B (zh) | 海杂波背景下子带自适应glrt‑ltd检测方法 | |
WO2019061006A1 (en) | METHOD AND DEVICE FOR DIAGNOSING BEARING FAULT, READABLE STORAGE MEDIUM, AND ELECTRONIC DEVICE | |
CN103926572B (zh) | 一种机载非正侧阵雷达用自适应子空间的杂波抑制方法 | |
CN106019256B (zh) | 基于自回归模型的雷达信号自适应检测方法 | |
CN103454622B (zh) | 基于稀疏约束的宽带雷达目标复回波去噪方法 | |
CN106646406B (zh) | 基于改进小波阈值去噪的外弹道测速雷达功率谱检测方法 | |
CN104076360B (zh) | 基于压缩感知的二维sar稀疏目标成像方法 | |
CN108931766B (zh) | 一种基于稀疏重构的非均匀stap干扰目标滤除方法 | |
CN103926599A (zh) | 基于emd迭代阈值滤波的gnss多径效应抑制方法 | |
CN104111449B (zh) | 一种改进的基于广义内积的空时二维自适应处理方法 | |
CN105699950B (zh) | 基于自适应迭代前后向平滑共轭梯度的雷达杂波抑制方法 | |
CN110221256A (zh) | 基于深度残差网络的sar干扰抑制方法 | |
CN112612006B (zh) | 基于深度学习的机载雷达非均匀杂波抑制方法 | |
CN109001687A (zh) | 基于广义旁瓣相消结构的机载雷达空时自适应滤波方法 | |
CN107664759A (zh) | 一种基于稀疏非均匀样本选优的机载雷达动目标检测方法 | |
CN105116388A (zh) | 基于鲁棒主成分分析的天波超视距雷达瞬态干扰抑制方法 | |
EP4083654A1 (en) | Radar communication with interference suppression | |
CN103777189A (zh) | 一种基于信息几何复自回归模型的雷达弱目标检测方法 | |
CN106772302A (zh) | 一种复合高斯背景下的知识辅助stap检测方法 | |
CN106772275A (zh) | 一种基于Alpha稳定分布的低秩矩阵恢复检测方法 | |
CN105158749A (zh) | 高频雷达海杂波幅度统计分布检验方法 | |
CN104793210B (zh) | 基于压缩感知的机载相控阵雷达低空风切变风速估计方法 | |
CN103760531B (zh) | 一种窄带雷达的空中目标回波的噪声抑制方法 | |
CN106950546B (zh) | 基于马氏距离重加权的非均匀杂波抑制方法 | |
CN107132518A (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 |