CN113673312A - 基于深度学习的雷达信号脉内调制识别方法 - Google Patents
基于深度学习的雷达信号脉内调制识别方法 Download PDFInfo
- Publication number
- CN113673312A CN113673312A CN202110762981.3A CN202110762981A CN113673312A CN 113673312 A CN113673312 A CN 113673312A CN 202110762981 A CN202110762981 A CN 202110762981A CN 113673312 A CN113673312 A CN 113673312A
- Authority
- CN
- China
- Prior art keywords
- signal
- time
- radar
- pulse modulation
- intra
- 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 43
- 238000013135 deep learning Methods 0.000 title claims abstract description 19
- 230000006870 function Effects 0.000 claims abstract description 58
- 238000012545 processing Methods 0.000 claims abstract description 31
- 230000005855 radiation Effects 0.000 claims abstract description 11
- 238000005070 sampling Methods 0.000 claims abstract description 11
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 238000012549 training Methods 0.000 claims description 26
- 238000012360 testing method Methods 0.000 claims description 25
- 238000010606 normalization Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000005315 distribution function Methods 0.000 claims description 6
- 238000011176 pooling Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 15
- 238000007781 pre-processing Methods 0.000 abstract description 6
- 230000002401 inhibitory effect Effects 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 11
- 230000008569 process Effects 0.000 description 10
- 230000004913 activation Effects 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000013527 convolutional neural network Methods 0.000 description 4
- 230000009977 dual effect Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 2
- 230000001965 increasing effect Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001149 cognitive effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000006467 substitution reaction 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/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
-
- 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/414—Discriminating targets with respect to background clutter
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Biomedical Technology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于雷达信号处理技术领域,公开了一种基于深度学习的雷达信号脉内调制识别方法,包括以下步骤:S1、对雷达采样信号进行滤波和去噪预处理;S2、将预处理后的采样信号进行Cohen类时频分布处理,得到时频图像;S3、对时频图像进行处理后,输入训练好的DCNN‑C网络模型中,通过网络模型自动判断输入雷达脉内调制信号类型,完成识别;所述DCNN‑C网络模型包括DCNN网络和与之拼接的分类网络。本发明利用双球型核函数进行Choi‑Williams时频分布处理,对雷达信号抑制交叉项的效果更好,信号鲁棒性特征更加明显;本发明利用DCNN网络对时频图像实现背景去噪处理,可以有效免去时频预处理造成的信号能量的信息损失,从而提高辐射源识别的准确度,而且简单实用且有效。
Description
技术领域
本发明属于雷达信号处理技术领域,具体涉及一种基于深度学习的雷达信号脉内调制识别方法。
背景技术
雷达信号脉内调制识别在信息化战役中发挥重要作用,针对敌人的雷达通信系统进行电子对抗,侦查地方情报,对增强军事武器装备的抗敌能力具有重要意义,同时在民用网络安全接入、认知无线电以及公共安全等领域具有较高的应用价值。
对于雷达信号脉内调制识别,传统做法是通过提取的到达时间(Time OfArrival,TOA)、到达角(Angle Of Arrival,AOA)、载波频率(Carrier Frequency,CF)、脉冲幅度(Pulse Amplitude,PA)和脉冲宽度(Pulse Width,PW)序列参数,构成脉冲描述字(Pulse Description Word,PDW),然后利用分类方法,如模板匹配法、支持向量等完成识别。随着电磁环境日益复杂,这五类参数,作为传统的识别方法已经不能满足当下所需要的识别精度。
为了有效改善信号特征提取依赖人工提取的缺点,随着人工智能的发展,深度学习逐渐被应用于雷达信号脉内调制识别方法中,并取得良好效果。有学者提出基于时频图像和卷积神经网络的信号自动识别系统,该方法是将信号通过时频分析转化为图像特征提取问题,通过对图像去噪等方式提高对图像的识别精度。时频图像与卷积神经网络相结合实现雷达信号脉内调制识别的技术考量了图像统计特征,但是提取的特征并不能完全表征信号的本质特征,且真实环境中的雷达回波信号背景杂波较为复杂,通过仿真实验得到的分析结果并不能完全表征信号的本质信息,尤其是在低信噪比的环境中受噪声影响大。同时,为去除时频图像中的噪声而采取过多的图像预处理去噪会导致部分有效信息丢失,降低识别的精确度,而随着卷积神经网络深度的增加同时也会造成参数大幅增加,提高计算复杂度。
发明内容
本发明提供了一种基于深度学习的雷达信号脉内调制识别方法,旨在解决现有技术中的信号处理效率低,原始信号信息丢失和信号识别准确度低的技术问题。
为了解决上述技术问题,本发明采用的技术方案为:一种基于深度学习的雷达信号脉内调制识别方法,包括以下步骤:
S1、对雷达采样信号进行滤波和去噪预处理;
S2、将预处理后的采样信号进行Cohen类时频分布处理,得到时频图像;所述进行Cohen类时频分布处理具体包括以下步骤:
S201、计算采样信号的模糊函数,计算公式为:
其中,s1为预处理后的雷达脉内调制信号表达式,其中,S1 *表示s1的共轭函数,τ为时延,v为频移,AF(v,τ)表示信号s1的模糊函数,t表示自变量时间;
S202、计算Choi-Williams时频分布函数,计算公式为:
其中,Φ(v,τ)表示核函数,其表达式为:
其中,α,β表示指数系数,用来调节核函数的大小,d表示距离中心的距离;
S3、对时频图像进行处理后,输入训练好的DCNN-C网络模型中,通过网络模型自动判断输入雷达脉内调制信号类型,完成识别;所述DCNN-C网络模型包括DCNN网络和与之拼接的分类网络。
所述的基于深度学习的雷达信号脉内调制识别方法,还包括对所述DCNN-C网络模型进行训练的步骤。
对所述DCNN-C网络模型进行训练的步骤具体为:
建立雷达脉内调制信号的样本数据库,所述雷达脉内调制信号的样本数据库中包含12种类别雷达脉内调制信号的灰度时频图像;将样本数据分为测试集和训练集;
建立DCNN-C网络;
根据所述训练集样本数据以及目标雷达脉内调制信号源信息对上述网络进行训练,
通过测试集数据以及目标雷达脉内调制信号源信息对所述雷达辐射源识别模型进行测试,直至测试集的精度达到90%以上。
建立雷达脉内调制信号的样本数据库的具体方法为:
在-10dB~10dB范围内,间隔2dB共11个信噪比,每个信噪比下分别对12类雷达脉内调制信号采集100个调制信号作为原始信号,得到共13200个随机雷达脉内调制信号;
对所述随机雷达脉内调制信号进行基于双球型核函数的Choi-Williams时频分布处理得到时频图像,对时频图像进行归一化和灰度处理得到对应的时频域灰度图像,将得到的13200张灰度时频图像作为DCNN-C网络模型的训练样本数据;
在-10dB~10dB范围内,间隔2dB共11个信噪比,每个信噪比下分别对12类雷达脉内调制信号采集20个调制信号作为原始信号,得到共2640个随机雷达脉内调制信号;
同样对随机雷达脉内调制信号进行基于双球型核函数的Choi-Williams时频分布处理得到时频图像,对时频图像进行归一化和灰度处理得到对应的时频域灰度图像,得到共2640张时频域灰度图像作为DCNN-C网络模型的测试集数据。
所述DCNN网络包括7个卷积层,每个卷积层后端都有一个ReLU函数进行非线性化,所述分类网络层包含256个尺寸为4x4的卷积核、池化层、全连接层和softmax层。
所述雷达信号包括NS信号、LFM信号、EQFM信号、2FSK信号、4FSK信号、DLFM信号、BPSK信号、Frank信号、LFM-BPSK信号和2FSK-BPSK信号。
所述步骤S3中,对时频图像进行处理具体包括:对其进行尺寸和幅度归一化调整,以及灰度化处理。
本发明与现有技术相比具有以下有益效果:
本发明提供了一种基于深度学习的雷达信号脉内调制识别方法,通过将待识别的雷达脉内调制信号进行预处理,然后将预处理后的信号进行Cohen类时频分析方法得到时频图像,对其进行归一化和灰度化,将得到的时频灰度图像输送到基于改进DCNN的分类网络训练,并保留效果良好时的网络参数;将测试集输入到已经训练好的网络模型,完成雷达脉内调制信号的识别分类。其中在时频分析部分,采用改进核函数的Choi-Williams分布时频分析方法,相对传统Cohen类时频分析中采用的指数型核函数对信号在模糊域进行能量提取,本发明采用双球型核函数,在时频分析能力上有显著的提升,即在低信噪比得到更加清晰的时频图像,这意味着对后续使用改进DCNN的分类网络时可以有效提高网络模型的识别率。
附图说明
图1为本发明实施例提供基于深度学习的雷达信号脉内调制信号识别方法的流程图;
图2为本发明实施例中采用的DCNN-C网络模型结构示意图;
图3为本发明实施例提供的时频域灰度图像(SNR=-6dB)的示意图;
图4为本发明实施例中的双球型核函数的示意图;图中,横轴代表时延,纵轴代表频移;
图5为正弦调制信号在模糊域的能量分布示意图;
图6为现有技术与本发明的时频分析方法的对比示意图;a为通过传统Choi-Williams时频分布的指数型核函数对正弦调频雷达脉内调制信号进行时频分析得到的图像,图中b为通过本发明Choi-Williams时频分布引入双球型核函数后得到的时频图像;
图7为采用本发明实施例的DCNN网络对SNR=-6dB下的4FSK信号的图像背景的去噪效果对比图,其中(a)为去噪前的图像,(b)为对应的去噪后的图像;
图8为采用本发明实施例的DCNN网络对SNR=-6dB下的SFM信号的图像背景的去噪效果对比图,其中(a)为去噪前的图像,(b)为对应的去噪后的图像;
图9为采用本发明实施例的DCNN网络对SNR=-6dB下的EQFM信号的图像背景的去噪效果对比图,其中(a)为去噪前的图像,(b)为对应的去噪后的图像。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例;基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为解决现有技术存在的雷达信号脉内调制识别处理过程的局限性,识别效果不佳和针对性不强的技术问题,本发明实施例提供了一种基于深度学习的雷达信号脉内调制识别方法,通过将待识别的雷达脉内调制信号进行滤波和去噪预处理,然后将预处理后的信号进行Cohen类时频分析方法得到时频图像,对其进行归一化和灰度化处理得到64*64的灰度时频图像,将其输送到DCNN-C(deep convolutional neural network-classification)网络模型进行训练,并保留效果良好时的网络参数;将测试集输入到已经训练好的网络模型,完成雷达脉内调制信号的识别分类。其中在时频分析部分,采用改进核函数的Cohen类时频分布方法,相对传统Choi-Williams时频分布中采用的指数型核函数对信号在模糊域进行能量提取,本发明采用双球型核函数,在时频分析能力上有显著的提升,即在低信噪比得到更加清晰的时频图像,这意味着对后续使用DCNN-C网络时可以有效提高网络模型的识别率。
具体地,如图1所示,本发明实施例提供的一种基于深度学习的雷达信号脉内调制识别方法,包括以下步骤:
S1、将截获的雷达信号s(t)进行预处理,对采样信号进行滤波和去噪预处理得到s1(t);
S2、将预处理后s1(t)信号进行Cohen类时频分布处理,得到时频图像;
本实施例中,使用的时频分析方法采用Choi-Williams时频分布,所述时频分布方法是引入一种指数型核函数,该核函数可以有效抑制那些由不同频率中心和不同时间中心的自项所产生的交叉项。所述Cohen类时频分布进行时频分析的过程如下:
AFs1(v,τ)为模糊函数,是对时变自相关函数作傅里叶变换,得到一个联合时频表达式,表达式如下:
上述式中,s1(t)为预处理后的雷达辐射源信号表达式,其中,S1 *表示s1的共轭函数,τ为时延,v为频移,AF(v,τ)表示信号s1的模糊函数,t表示时间。
Choi-Williams时频分布函数的表达式为:
其中,Φ(v,τ)表示核函数,C(t,w)表示雷达脉内调制信号的Choen类时频分布函数,通过上述时频分布函数,可以输出雷达脉内调制信号的时频图像。本实施例中,通过式(1)计算雷达脉内调制信号的模糊函数,可以将雷达脉内调制信号转化到模糊域,然后,通过式(2)的计算,即利用模糊函数与核函数相乘,实现在模糊域内提取能量,再将该部分进行傅里叶逆变换,得到时频图像。
具体地,本实施例采用的Choi-Williams时频分布处理方法中,采用双球型核函数,其表达式为:
其中,α,β用来调节核函数的大小,d表示距离中心的距离,核函数的宽度用四倍高斯核函数标准差来估计。图5所示,本发明实施例中的双球型核函数的示意图;图中,横轴代表时延,纵轴代表频移。
具体地,本实施例中,通过mesh()函数,并且取绝对值,可以将类时频分布函数C(t,w)转化为时频图像。
S3、将上述得到的时频图像进行归一化和灰度化处理得到灰度时频图像,将测试集灰度时频图像输送至已经训练好的DCNN-C网络模型;所述训练好的DCNN-C网络保留了训练效果良好时的权重参数,并且使用该参数对测试网络进行初始化。
具体地,将时频图像进行归一化和灰度化处理后,输送至DCNN-C网络模型,如图3所示,为归一化和灰度化处理后得到的时频域灰度图像(SNR=-6dB)的示意图。由于imagesc函数处理生成的时频域图像为彩色图像,可以通过灰度化处理将所述时频域图像处理生成单通道的时频域灰度图像,相比于处理三通道的彩色图像,卷积神经网络算法在处理时频域灰度图像时所需处理的数据量大大降低,但是同时不会影响最终网络模型对图像的识别准确率。
具体地,DCNN-C网络即为改进DCNN的分类网络,包括DCNN网络和与之拼接的分类网络,其训练方法如下:
S301、建立雷达脉内调制信号的样本数据库,所述雷达脉内调制信号的样本数据库中包含12种类别雷达脉内调制信号的灰度时频图像。也就是说,雷达脉内调制信号的样本数据库中包含时频灰度图像以及其相对应的目标雷达脉内调制信号源信息。
作为本发明的一种可行实施例,在-10dB到10dB信噪比下每隔2dB产生数量相等的不同调制方式的雷达样本信号,所述不同脉内调制方式的雷达样本信号包括常规脉冲信号,线性调频信号,正弦调频信号,频率编码信号、相位编码信号以及混合调制信号,上述雷达辐射源样本信号的参数如表1所示。
表1 仿真雷达信号参数表:
S302、将样本数据分为测试集和训练集。
本发明实例中,为使仿真的雷达脉内调制信号更接近与真实环境接收到的信号,本发明在仿真的雷达信号中添加高斯白噪声来模拟实际环境中的噪声影响,在-10dB~10dB范围内,间隔2dB共11个信噪比,每个信噪比下分别对12类雷达脉内调制信号采集100个调制信号作为原始信号,共13200个随机雷达脉内调制信号,通过双球型核函数的Choi-Williams时频分布处理得到时频图像,对其时频图像进行归一化和灰度处理后,得到了对应的时频域灰度图像,将得到的13200张灰度时频图像作为训练DCNN-C网络模型的样本数据。根据上述的样本选取规则,在-10dB~10dB范围内,间隔2dB共11个信噪比,每个信噪比下分别对12类雷达脉内调制信号采集20个调制信号作为原始信号,得到共2640个随机雷达脉内调制信号,再将通过基于双球型核函数的Choi-Williams时频分布处理得到的时频图像进行归一化和灰度化处理,得到时频域灰度图像共2640张作为测试集数据。
S303、建立基于改进DCNN的分类网络模型,即DCNN-C网络模型。
具体的,本实施例中,基于改进DCNN的分类网络中,DCNN网络包括7个卷积层,每个卷积层后端都有一个激活函数:ReLU函数进行非线性化,在DCNN后拼接分类网络以最终达到分类的效果。DCNN通过从污染图像中分离出噪声图像可以实现潜在的清晰图像。
给定网络层的深度d的DCNN,每一层都使用相同的卷积算子。卷积运算实际上是一种线性滤波,而线性特性并不能代表非平凡问题,因而在DCNN中,将每个卷积层的输出后都添加一个非线性激活函数,将其输出结果作为下一层的输入。
在上述DCNN网络中,如图2所示,使用整流线性单元ReLU作为每一层的激活函数,其是一种非线性映射,与其他激活函数相比,ReLU可以在大型复杂数据集上快速有效地训练深度神经网络结构。
DCNN-C网络模型中,损失函数的数学表达式为:
其中,θ=[w,b]是DCNN的网络参数,w是权重矩阵,b为偏置,在模型训练过程中,最小化损失函数以找到DCNN的最优参数θ。
本实施例中,DCNN-C网络模型的结构图如图2所示,其输入层大小是64*64*1;每个卷积层都采用64个5*5的卷积核;分类网络层包含256个尺寸为4x4的卷积核、池化层、全连接层和softmax进行最终的分类。
其中分类网络的池化层是非线性下采样的一种形式,主要作用是通过减少网络的参数来减小计算量,并且在一定程度上可以控制过拟合,该部分使用的平均池化;而全连接层用来计算向量和权重向量之间的点积,再加上一个偏置,结果通过sigmoid函数输出,根据12种雷达脉内调制信号,该部分包括12个节点。
S304、根据所述训练集样本数据以及目标雷达脉内调制信号源信息对上述网络进行训练,本申请中设置的迭代次数为400次,每次迭代读取33张数据图像。迭代过程中使用交叉熵函数计算损失,反向传播过程中梯度下降采用Adam优化算法,在计算过程中会对初始化的参数进行调整,同时计算精度,训练完成后,保存权重参数。
S305、通过测试集数据以及目标雷达脉内调制信号源信息对所述雷达辐射源识别模型进行测试,在测试过程中,设置的迭代次数为400次,初始化参数采用训练网络保留的权重参数,将2640张时频域灰度图像作为测试集输入到DCNN-C网络中,通过交叉熵函数计算损失,反向传播过程中梯度下降采用Adam优化算法,同时计算精度,整体正确率达到预期值时训练完成,即已达到本申请所需精度的要求。
重复上述所有步骤并完善参数,直到测试的整体正确率达到预期值时,完成训练。
具体地,上述预期值指测试集的精度达到90%以上。
上述模型训练过程中,当损失值较大时,表明DCNN-C网络的识别准确率不高,识别出的雷达脉内调制信息与实际匹配的信息存在较大的差异,需要使用反向传播算法对所述雷达辐射源模型的参数进行调整;当损失值较小时,表明DCNN-C网络模型的识别准确率足够,对输入到模型中的雷达脉内调制信息识别与实际信息基本相符;
本发明实例通过上述13200个样本数据训练的雷达辐射源模型可以达到对测试集90%以上的精度,且易错率低。
训练完成后的DCNN-C网络模型,可以对时频域灰度图像的进行调制类型判断,也就是说,将时频域灰度图像直接输入到DCNN-C网络模型中,可以输出对应的响应雷达脉内调制信号类型。
所述步骤S3中,所述DCNN-C网络中进行类别判断的雷达脉内调制信号包括表1中的12种类型,具体包括NS信号、LFM信号、EQFM信号、2FSK信号、4FSK信号、DLFM信号、BPSK信号、Frank信号、LFM-BPSK信号和2FSK-BPSK信号。
本实施例中,时频分析函数的核函数在设计时根据信号在模糊域的能量分布特性,把双球型核函数引入到Choi-Williams时频分布中,用于将预处理后的采样信号进行时频处理,以便得到时频图像。
如图5所示,为正弦调制信号(SFM)在模糊域的能量分布示意图,SFM信号的模糊域能量分布集中在中间部分,该信号的模糊域能量分布说明,本发明实施例请中选取的双球型核函数可以更好的提取模糊域能量。
图6为本发明的识别方法与现有技术的对比示意图,图中采用的数据源为根据表1中的参数随机产生的正弦调频雷达脉内调制信号,图中a为通过传统Choi-Williams时频分布的指数型核函数进行时频分析得到的图像,图中b为通过本发明Choi-Williams时频分布引入双球型核函数后得到的时频图像。图6说明,本发明的时频处理方法采用双球型核函数,可以得到背景更清晰的时频图像,通过后续利用DCNN-C网络进行进一步去噪和分类,具有更好的识别效果。
图7~9所示,为本发明实施例采用的DCNN网络对SNR=-6dB下的4FSK,SFM,EQFM三类图像背景去噪效果示意图,从图中可以看出,本发明采用的DCNN可以较好地对雷达信号时频图像背景去噪。
现有技术方法采用Choi-Williams时频分布、快速傅里叶变换以及Wigner-Ville分布得到各类调制信号的时频图像,然后将时频图像输入到卷积神经网络中进行识别。例如,基于扩张残差网络的雷达辐射源信号识别(参见秦鑫,黄洁,查雄,等.基于扩张残差网络的雷达辐射源信号识别[J].电子学报,2020,48(3):456-462.DOI:10.3969/j.issn.0372-2112.2020.03.006),采用Choi-Williams时频分布得到时频图像,该方法在信噪比为-6dB时可实现对雷达辐射源信号整体识别率达到98.2%。综合来看,现有技术的最大识别率约为97%,平均识别率约为90%,而本发明为有效针对低信噪比(SNR=-10dB、SNR=-8dB)的雷达信号,采用的双球型核函数在不影响其他脉内调制类型信号时频图像效果的情况下,可以提高对部分信号(例如正弦调频信号SFM)的分析效果,采用DCNN-C网络模型对包括低信噪比(SNR=-10dB、SNR=-8dB)下信号的整体最大识别率可以达到98.5%,多次测试的平均识别率能够达到95.8%,实验数据说明本发明可以有效应对低信噪比下的雷达脉内调制信号识别。上述识别率η的计算方法为:
其中N0指全部数据集,本发明测试过程中N0=2640,N指本发明DCNN-C网络模型能够识别准确的数据。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (7)
1.一种基于深度学习的雷达信号脉内调制识别方法,其特征在于,包括以下步骤:
S1、对雷达采样信号进行滤波和去噪预处理;
S2、将预处理后的采样信号进行Cohen类时频分布处理,得到时频图像;所述行Cohen类时频分布处理具体包括以下步骤:
S201、计算采样信号的模糊函数,计算公式为:
其中,s1为预处理后的雷达脉内调制信号表达式,其中, s1 *表示s1的共轭函数,τ为时延, v为频移,AF(v,τ)表示信号s 1的模糊函数,t表示自变量时间;
S202、计算Choi-Williams时频分布函数,计算公式为:
其中,α,β表示指数系数,用来调节核函数的大小,d表示距离中心的距离;
S3、对时频图像进行处理后,输入训练好的DCNN-C网络模型中,通过网络模型自动判断输入雷达脉内调制信号类型,完成识别;所述DCNN-C网络模型包括DCNN网络和与之拼接的分类网络。
2.根据权利要求1所述的基于深度学习的雷达信号脉内调制识别方法,其特征在于,还包括对所述DCNN-C网络模型进行训练的步骤。
3.根据权利要求2所述的基于深度学习的雷达信号脉内调制识别方法,其特征在于,对所述DCNN-C网络模型进行训练的步骤具体为:
建立雷达脉内调制信号的样本数据库,所述雷达脉内调制信号的样本数据库中包含12种类别雷达脉内调制信号的灰度时频图像;将样本数据分为测试集和训练集;
建立DCNN-C网络;
根据所述训练集样本数据以及目标雷达脉内调制信号源信息对上述网络进行训练,
通过测试集数据以及目标雷达脉内调制信号源信息对所述雷达辐射源识别模型进行测试,直至测试集的精度达到90%以上。
4.根据权利要求3所述的基于深度学习的雷达信号脉内调制识别方法,其特征在于,建立雷达脉内调制信号的样本数据库的具体方法为:
在-10dB~10dB范围内,间隔2dB共11个信噪比,每个信噪比下分别对12类雷达脉内调制信号采集100个调制信号作为原始信号,得到共13200个随机雷达脉内调制信号;
对所述随机雷达脉内调制信号进行基于双球型核函数的Choi-Williams时频分布处理得到时频图像,对时频图像进行归一化和灰度处理得到对应的时频域灰度图像,将得到的13200张灰度时频图像作为DCNN-C网络模型的训练样本数据;
在-10dB~10dB范围内,间隔2dB共11个信噪比,每个信噪比下分别对12类雷达脉内调制信号采集20个调制信号作为原始信号,得到共2640个随机雷达脉内调制信号;
同样对随机雷达脉内调制信号进行基于双球型核函数的Choi-Williams时频分布处理得到时频图像,对时频图像进行归一化和灰度处理得到对应的时频域灰度图像,得到共2640张时频域灰度图像作为DCNN-C网络模型的测试集数据。
5.根据权利要求1所述的基于深度学习的雷达信号脉内调制识别方法,其特征在于,所述DCNN网络包括7个卷积层,每个卷积层后端都有一个ReLU函数进行非线性化,所述分类网络层包含256个尺寸为4x4的卷积核、池化层、全连接层和softmax层。
6.根据权利要求1所述的基于深度学习的雷达信号脉内调制识别方法,其特征在于,所述雷达信号包括NS信号、LFM信号、EQFM信号、2FSK信号、4FSK信号、DLFM信号、BPSK信号、Frank信号、LFM-BPSK信号和2FSK-BPSK信号。
7.根据权利要求1所述的基于深度学习的雷达信号脉内调制识别方法,其特征在于,所述步骤S3中,对时频图像进行处理具体包括:对其进行尺寸和幅度归一化调整,以及灰度化处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110762981.3A CN113673312B (zh) | 2021-07-06 | 2021-07-06 | 基于深度学习的雷达信号脉内调制识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110762981.3A CN113673312B (zh) | 2021-07-06 | 2021-07-06 | 基于深度学习的雷达信号脉内调制识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113673312A true CN113673312A (zh) | 2021-11-19 |
CN113673312B CN113673312B (zh) | 2023-08-22 |
Family
ID=78538616
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110762981.3A Active CN113673312B (zh) | 2021-07-06 | 2021-07-06 | 基于深度学习的雷达信号脉内调制识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113673312B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115392325A (zh) * | 2022-10-26 | 2022-11-25 | 中国人民解放军国防科技大学 | 一种基于CycleGan的多特征降噪调制识别方法 |
CN115659162A (zh) * | 2022-09-15 | 2023-01-31 | 云南财经大学 | 雷达辐射源信号脉内特征提取方法、系统及设备 |
CN115856987A (zh) * | 2023-02-28 | 2023-03-28 | 西南科技大学 | 一种复杂环境下的核脉冲信号与噪声信号甄别方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107301432A (zh) * | 2017-07-11 | 2017-10-27 | 哈尔滨工程大学 | 基于时频分析的自适应辐射源调制识别方法 |
US20180357542A1 (en) * | 2018-06-08 | 2018-12-13 | University Of Electronic Science And Technology Of China | 1D-CNN-Based Distributed Optical Fiber Sensing Signal Feature Learning and Classification Method |
CN109871757A (zh) * | 2019-01-11 | 2019-06-11 | 中国船舶重工集团公司第七二四研究所 | 一种基于联合时频特征的雷达信号脉内调制类型识别方法 |
CN110175560A (zh) * | 2019-05-24 | 2019-08-27 | 哈尔滨工程大学 | 一种雷达信号脉内调制识别方法 |
CN112487933A (zh) * | 2020-11-25 | 2021-03-12 | 中国人民解放军海军工程大学 | 一种基于自动化深度学习的雷达波形识别方法及系统 |
-
2021
- 2021-07-06 CN CN202110762981.3A patent/CN113673312B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107301432A (zh) * | 2017-07-11 | 2017-10-27 | 哈尔滨工程大学 | 基于时频分析的自适应辐射源调制识别方法 |
US20180357542A1 (en) * | 2018-06-08 | 2018-12-13 | University Of Electronic Science And Technology Of China | 1D-CNN-Based Distributed Optical Fiber Sensing Signal Feature Learning and Classification Method |
CN109871757A (zh) * | 2019-01-11 | 2019-06-11 | 中国船舶重工集团公司第七二四研究所 | 一种基于联合时频特征的雷达信号脉内调制类型识别方法 |
CN110175560A (zh) * | 2019-05-24 | 2019-08-27 | 哈尔滨工程大学 | 一种雷达信号脉内调制识别方法 |
CN112487933A (zh) * | 2020-11-25 | 2021-03-12 | 中国人民解放军海军工程大学 | 一种基于自动化深度学习的雷达波形识别方法及系统 |
Non-Patent Citations (1)
Title |
---|
秦鑫;黄洁;查雄;骆丽萍;胡德秀;: "基于扩张残差网络的雷达辐射源信号识别", 电子学报, no. 03 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115659162A (zh) * | 2022-09-15 | 2023-01-31 | 云南财经大学 | 雷达辐射源信号脉内特征提取方法、系统及设备 |
CN115659162B (zh) * | 2022-09-15 | 2023-10-03 | 云南财经大学 | 雷达辐射源信号脉内特征提取方法、系统及设备 |
CN115392325A (zh) * | 2022-10-26 | 2022-11-25 | 中国人民解放军国防科技大学 | 一种基于CycleGan的多特征降噪调制识别方法 |
CN115392325B (zh) * | 2022-10-26 | 2023-08-18 | 中国人民解放军国防科技大学 | 一种基于CycleGan的多特征降噪调制识别方法 |
CN115856987A (zh) * | 2023-02-28 | 2023-03-28 | 西南科技大学 | 一种复杂环境下的核脉冲信号与噪声信号甄别方法 |
CN115856987B (zh) * | 2023-02-28 | 2023-05-02 | 西南科技大学 | 一种复杂环境下的核脉冲信号与噪声信号甄别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113673312B (zh) | 2023-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107220606B (zh) | 基于一维卷积神经网络的雷达辐射源信号的识别方法 | |
CN113673312B (zh) | 基于深度学习的雷达信号脉内调制识别方法 | |
CN109597043B (zh) | 基于量子粒子群卷积神经网络的雷达信号识别方法 | |
CN114429156B (zh) | 雷达干扰多域特征对抗学习与检测识别方法 | |
CN109993280A (zh) | 一种基于深度学习的水下声源定位方法 | |
CN108169708B (zh) | 模块化神经网络的直接定位方法 | |
CN110222748B (zh) | 基于1d-cnn多域特征融合的ofdm雷达信号识别方法 | |
CN109902715B (zh) | 一种基于上下文聚合网络的红外弱小目标检测方法 | |
CN109188414A (zh) | 一种基于毫米波雷达的手势运动探测方法 | |
CN114595732B (zh) | 基于深度聚类的雷达辐射源分选方法 | |
CN111722199A (zh) | 一种基于卷积神经网络的雷达信号检测方法 | |
CN109711314B (zh) | 基于特征融合与sae的雷达辐射源信号分类方法 | |
CN116047427B (zh) | 一种小样本雷达有源干扰识别方法 | |
CN114726692B (zh) | 一种基于SEResNet-LSTM的辐射源调制方式识别方法 | |
CN112859012A (zh) | 一种基于级联卷积神经网络的雷达欺骗干扰识别方法 | |
CN112749633A (zh) | 分离与重构的个体辐射源识别方法 | |
CN115630276A (zh) | 一种基于YOLOv5的雷达脉内调制类型快速识别方法 | |
CN112307926B (zh) | 一种基于生成对抗网络的声学被动式舰船目标分类方法 | |
CN116482618B (zh) | 基于多损失特征自校准网络的雷达有源干扰识别方法 | |
CN112213697A (zh) | 一种基于贝叶斯决策理论用于雷达欺骗干扰识别的特征融合方法 | |
CN113359091B (zh) | 基于深度学习的多核函数混叠雷达辐射源识别方法和系统 | |
CN115223033A (zh) | 一种合成孔径声呐图像目标分类方法及系统 | |
Fan et al. | Weak target detection based on deep neural network under sea clutter background | |
CN112666528B (zh) | 一种基于卷积神经网络的多站雷达系统干扰鉴别方法 | |
CN116112932B (zh) | 一种数据知识双驱动的射频指纹识别方法与系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |