CN104573303A - 未知辐射源信号信噪比的估计方法 - Google Patents

未知辐射源信号信噪比的估计方法 Download PDF

Info

Publication number
CN104573303A
CN104573303A CN201410221363.8A CN201410221363A CN104573303A CN 104573303 A CN104573303 A CN 104573303A CN 201410221363 A CN201410221363 A CN 201410221363A CN 104573303 A CN104573303 A CN 104573303A
Authority
CN
China
Prior art keywords
signal
formula
ratio
noise ratio
noise
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.)
Pending
Application number
CN201410221363.8A
Other languages
English (en)
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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201410221363.8A priority Critical patent/CN104573303A/zh
Publication of CN104573303A publication Critical patent/CN104573303A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Noise Elimination (AREA)

Abstract

本发明公开了一种未知辐射源信号信噪比的估计方法，该方法包括以下主要步骤：(1)对接收到的未知辐射源信号进行WVD变换处理，并取零频处的信息；(2)对零频处信息的幅度进行归一化处理；(3)通过相关数学变换计算去噪后信号的功率；(4)依据信噪比估计公式得到接收信号的信噪比。

Description

未知辐射源信号信噪比的估计方法
技术领域
本发明设计是一种未知辐射源信号信噪比的估计方法，能够在没有先验知识的情况下，估算出与信号互不相关、不同时刻相互独立的各种加性噪声的大小。
背景技术
在现代雷达和通信电子战中，对辐射源信号进行个体特征的提取与识别具有重要意义，利用接收信号的信噪比作为个体特征对辐射源进行唯一识别已成为研究热点。当前用于信噪比估计的方法可以分为两类：一类是基于先验知识的估计算法，例如M2M4算法，需预先知道信号的包络信息；另一类则是无需先验知识的估计算法，例如基于特征值分解的算法，此类算法无需任何先验知识，但噪声类型仅限于高斯白噪声。对当前的辐射源信号进行信噪比的估计需满足两个条件：一是无需任何先验知识；二是适用于高斯、非高斯等各种加性噪声。因此，亟需一种能够在没有先验知识的前提下，准确估算各类高斯、非高斯等加性噪声大小的方法。
本发明针对上述问题，对接收到的信号首先进行威格纳-维尔分布变换(WVD)，取WVD变换零频处信息，此处信息有效抑制了与信号互不相关的各种加性噪声，但中心频率变为原始信号的两倍，幅度也为原来的常数倍。然后对幅度进行归一化处理，根据数学变换可以得到去噪后信号的功率，最后依据信噪比估计公式得到接收信号的信噪比。该方法无需任何先验知识，且适用于与信号互不相关、不同时刻相互独立的各种加性噪声，并通过仿真验证了该方法的准确性和可行性。
发明内容
本发明的方法包括以下技术措施：
(1)基于计算机仿真平台，建立辐射源信号信噪比估计的处理模型；
(2)对接收到的未知辐射源信号进行WVD变换处理，并取零频处的信息；
(3)对零频处信息的幅度进行归一化处理；
(4)通过相关数学变换计算去噪后信号的功率；
(5)依据信噪比估计公式得到接收信号的信噪比。
本发明相比背景技术具有如下优点：
(1) 本发明工程实现容易，成本低；
(2) 本发明具有普适性。
附图说明
附图1是本发明的流程图。参照附图1，本发明的流程由建立模型、提取WVD变换零频处信息、归一化处理、计算去噪信号功率、信噪比估计以及获得最终的信噪比值6个部分组成。其中1用于建立辐射源信号信噪比估计的处理模型；2用于对接收到的未知辐射源信号进行WVD变换处理，并取零频处的信息；3用于对零频处信息的幅度进行归一化处理；4用于通过相关数学变换计算去噪后信号的功率；5用于依据信噪比估计公式得到接收信号的信噪比；6用于获得最终的信噪比值。
附图2、3是指以线性调频信号为例，在每种信噪比下运算100次，估计出信噪比的均值和标准偏差大小的结果图。其中，图2横坐标是指接收到未知信号的实际信噪比大小，纵坐标表示估计出的信噪比均值；图3横坐标是指接收到未知信号的实际信噪比大小，纵坐标表示估计出的信噪比偏差。
具体实施方式
实施本发明的原理如下：对接收信号首先提取其WVD变换零频处信息，进行归一化处理后通过数学变换近似得到接收信号中原始信号的功率，最后根据信噪比计算公式得到信噪比的估计值。本发明旨在没有任何先验知识的前提出，准确估算出与信号互不相关、不同时刻相互独立的各种加性噪声的大小值，为辐射源进行唯一识别打下基础。
(1) 提取WVD变换零频处信息。设信号                                                ,为调频系数，为角频率，WVD变换定义为
                    （1）
将信号写成其复数形式
                          （2）
则的WVD变换如下所示
令公式(3)中的，可得
令(4)式中的， ，可得
          （5）
上式中。同理可得，当时，有
                             (6)
由以上公式可知，信号在WVD变换零频处的幅度为原始信号的常数倍，中心频率为2倍。显然，对上述WVD变换的零频处可定义为
                         (7)
由公式(7)可知，对能量有限的信号而言，其WVD变换的零频处实际就是其对称相关函数。假设是周期为的实信号，则的为
                         (8)
令        
                                    (9)
其中为零均值的有用实信号，为与信号互不相关、不同时刻相互独立的加性随机噪声，则可得到
      （10）
由公式(10)可知，对求取其对称相关函数后，信号中的噪声基本被有效抑制，并且越大，抑制的效果越好。
(2) 归一化处理。由公式(5)、(6) 可知，求取信号的后，幅度变为原信号的常数倍，中心频率变为原信号的2倍，如何利用变换后得到的对信号进行信噪比的估计呢?假设原始信号的表达式如下
                                                                                                                        (11)
其中,为与信号互不相关、不同时刻相互独立的加性噪声，对按公式(3)、(4)进行处理，可得
                                                                                            (12)
再对进行归一化处理得到
                   (13)
（3）计算去噪信号功率。信噪比的估计方法如下式所示
                                                                                                     (14)
式中是受噪声污染后的信号的功率，是未受噪声污染后的信号的功率，令和分别为信号和经采样后的信号序列，采样点数为N，则和的表达式分别如下
                                                                                                                        (15)
                                                                                                                         (16)
对接收到的未知信号进行信噪比估计时，已知，但未知。
令，和分别为和经采样后的信号序列，采样点数为N。由数学变换可得到如下公式
                                                                                                                      (17)
由公式(17)可得
                                                                                                       (18)
依据公式(16)可得
                                                                                                                (19)
（4）信噪比估计。最终由公式(14)可得到信噪比的估计大小如下
                (20)
下面结合实例说明一下整个发明的优势。
假设原始信号为线性调频信号，中心频率为，带宽，脉宽。仿真时加上由高斯白噪声和服从韦布尔分布的非高斯噪声所产生的混合噪声。每次仿真所加的高斯白噪声和服从韦布尔分布的非高斯噪声的能量相等，与原始信号相比的信噪比分别为10dB、15dB、20dB、25dB、30dB、35dB，最终产生的混合噪声大小如表1所示。服从韦布尔分布的噪声的形状参数为4、强度参数为2。混合噪声的大小按公式(21)计算，式中和为分别加高斯白噪声和韦布尔噪声时的信噪比大小。
                              (21)
表1 原始信号所加噪声的大小
高斯白噪声/dB 10 15 20 25 30 35
韦布尔噪声/dB 10 15 20 25 30 35
混合噪声/dB 6.99 11.99 16.99 21.99 26.99 31.99
利用本发明对接收到的未知辐射源信号进行信噪比的估计，每种信噪比下计算100次，估计信噪比的均值如图2所示，标准偏差如图3所示。由图2和图3可知，利用本发明估计出的信噪比具有较高的准确度和稳定度。随着信噪比的增高，准确度和稳定度也不断提升，可以为未知辐射源信号的唯一识别提供支撑。

Claims (4)

1.对接收到的信号                                                首先进行威格纳-维尔分布 (WVD)变换，取WVD变换零频处信息，记为。
2.对零频处信息的幅度进行归一化处理，记为。
3.根据相关数学变换得到去噪后信号的功率：
信噪比的估计方法如下式所示
      (1)
式中是受噪声污染后的信号的功率，是未受噪声污染后的信号的功率，令和分别为信号和经采样后的信号序列，采样点数为N，则和的表达式分别如下
    (2)
    (3)
对接收到的未知信号进行信噪比估计时，已知，但未知，令，和分别为和经采样后的信号序列，采样点数为N，由数学变换可得到如下公式
       (4)
由公式(4)可得
       (5)
依据公式(3)可得
(6)
公式（6）即为的估计值。
4.基于信噪比估计公式求得未知信号的信噪比。
CN201410221363.8A 2014-05-25 2014-05-25 未知辐射源信号信噪比的估计方法 Pending CN104573303A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410221363.8A CN104573303A (zh) 2014-05-25 2014-05-25 未知辐射源信号信噪比的估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410221363.8A CN104573303A (zh) 2014-05-25 2014-05-25 未知辐射源信号信噪比的估计方法

Publications (1)

Publication Number Publication Date
CN104573303A true CN104573303A (zh) 2015-04-29

Family

ID=53089356

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410221363.8A Pending CN104573303A (zh) 2014-05-25 2014-05-25 未知辐射源信号信噪比的估计方法

Country Status (1)

Country Link
CN (1) CN104573303A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106059722A (zh) * 2016-05-20 2016-10-26 南京长源智能科技有限公司 一种盲信号的信噪比估计方法
CN106125049A (zh) * 2016-06-08 2016-11-16 陈昌孝 雷达辐射源信号个体差异分析方法
CN108918927A (zh) * 2018-08-24 2018-11-30 广东石油化工学院 一种能量分解中功率信号滤波方法及系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102466795A (zh) * 2010-11-09 2012-05-23 何明浩 基于多指标的雷达信号脉内特征参数评估方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102466795A (zh) * 2010-11-09 2012-05-23 何明浩 基于多指标的雷达信号脉内特征参数评估方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吕明久，等.: "基于对称相关法的低信噪比信号时频分析", 《空间雷达学院学报》 *
朱振波，等.: "一种低信噪比下雷达脉内特征提取方法", 《航天电子对抗》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106059722A (zh) * 2016-05-20 2016-10-26 南京长源智能科技有限公司 一种盲信号的信噪比估计方法
CN106125049A (zh) * 2016-06-08 2016-11-16 陈昌孝 雷达辐射源信号个体差异分析方法
CN106125049B (zh) * 2016-06-08 2019-01-01 陈昌孝 雷达辐射源信号个体差异分析方法
CN108918927A (zh) * 2018-08-24 2018-11-30 广东石油化工学院 一种能量分解中功率信号滤波方法及系统

Similar Documents

Publication Publication Date Title
CN106501800B (zh) 基于代价参考粒子滤波的mimo雷达目标检测前跟踪方法
CN107124762B (zh) 一种高效消减非视距误差的无线定位方法
CN104573303A (zh) 未知辐射源信号信噪比的估计方法
CN104363194A (zh) 基于波形变换的psk调制识别方法
Wang et al. Parameter estimation of chirp signal under low SNR
CN105929380A (zh) 一种卫星激光高度计全波形激光雷达数据去噪方法
Hao et al. An efficiency-improved Tdoa-based direct position determination method for multiple sources
Zeng et al. The de-noising algorithm based on intrinsic time-scale decomposition
CN110568415B (zh) 混合高斯模型下基于Arctan函数的信号检测方法
CN107607913B (zh) 基于对数累积量的海杂波Pareto分布参数估计方法
CN109460614B (zh) 基于瞬时带宽的信号时间-频率分解方法
CN103647633A (zh) 适用尖峰噪声环境的时延估计方法
CN104793197A (zh) 基于ifft频谱相除法和梯度自适应格型滤波的直达波抑制方法
CN106226752A (zh) 一种扩展海杂波Pareto分布参数估计范围的方法
CN112462344B (zh) 通过时频域变换处理提取雷达信号脉内特征参数的方法
CN114358046A (zh) 一种多复杂度等级复杂电磁干扰环境模拟生成方法及系统
CN112014811B (zh) 一种雷达载波频率的精细估计方法
Ananieva et al. Design of a device for optimal reception of signals against the background of a two-component Markov interference
CN105403774B (zh) 一种基于排除偏差的雷电定向定位优化方法
Qiu et al. Multiple BPSK/LFM hybrid modulated signals parameter estimation and analysis intercepted by non-cooperative radar receiver
CN110198196B (zh) 通信系统中基于信号强度的振幅估计方法
Tian et al. Modified Cauchy distribution model of high-order passive intermodulation
CN102819012B (zh) 一种基于改进的非负ica的相控阵风廓线雷达信号处理方法
CN109683142A (zh) 基于差分包络检波的三角线性调频连续信号参数估计方法
CN104182971A (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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20150429