CN108761441A - 一种雷达信号相位编码规律识别方法及装置 - Google Patents
一种雷达信号相位编码规律识别方法及装置 Download PDFInfo
- Publication number
- CN108761441A CN108761441A CN201810615792.1A CN201810615792A CN108761441A CN 108761441 A CN108761441 A CN 108761441A CN 201810615792 A CN201810615792 A CN 201810615792A CN 108761441 A CN108761441 A CN 108761441A
- Authority
- CN
- China
- Prior art keywords
- radar
- interval
- pulse signal
- phase comparing
- journal
- 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
Links
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
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
本发明涉及一种雷达信号相位编码规律识别方法及装置,该方法采集雷达脉冲信号;对采集的雷达脉冲信号进行小波滤波处理,得到滤波后的雷达脉冲信号;确定滤波后的雷达脉冲信号的突变点以及各个突变点对应的位置;计算相邻两个突变点之间的间隔,顺序记录间隔;根据顺序记录的间隔,识别雷达信号的编码规律。本发明简单可靠,使得雷达号识别的准确率较高。
Description
技术领域
本发明属于电子侦察信号分析技术领域,具体涉及一种雷达信号相位编码规律识别方法及装置。
背景技术
对接收到的雷达辐射源信号进行精确处理是电子侦察的重要内容,也是信息保障的重要手段。其中,对雷达信号辐射源的编码规律进行精确分析与处理是其中的重要组成部分,能够获得雷达辐射源信号的精确参数,是进行雷达信号识别和分选处理的重要依据,具有重要的意义。
传统的基于脉冲描述字的五种特征参数识别方法的可靠性急剧下降,已很难适应需求。随着技术不断进步,新体制雷达层出不穷,雷达信号形式日益复杂,复杂体制雷达辐射源迅速增加并逐渐占据主导地位。其中,相位编码雷达信号是一类重要的调制信号,对其进行编码规律的识别和判决是重要的信号精确分析技术。现有的识别算法只能在高信噪比条件下进行处理,难以适应低信噪比雷达信号的识别需求,存在识别错误率高的问题。
发明内容
本发明的目的在于提供一种雷达信号相位编码规律识别方法及装置,用以解决现有技术中雷达信号识别错误率高的问题。
为解决上述技术问题,本发明的技术方案为:
本发明提供了一种雷达信号相位编码规律识别方法,包括如下步骤:
采集雷达脉冲信号;
对采集的雷达脉冲信号进行小波滤波处理,得到滤波后的雷达脉冲信号;
确定滤波后的雷达脉冲信号的突变点以及各个突变点对应的位置;
计算相邻两个突变点之间的间隔,顺序记录间隔;
根据顺序记录的间隔,识别雷达信号的编码规律。
本发明还提供了一种雷达信号相位编码规律识别装置,包括处理器,所述处理器用于执行指令实现如下方法:
采集雷达脉冲信号;
对采集的雷达脉冲信号进行小波滤波处理,得到滤波后的雷达脉冲信号;
确定滤波后的雷达脉冲信号的突变点以及各个突变点对应的位置;
计算相邻两个突变点之间的间隔,顺序记录间隔;
根据顺序记录的间隔,识别雷达信号的编码规律。
本发明的有益效果:
本发明的雷达信号相位编码规律识别方法及装置,对采集的雷达脉冲信号进行小波处理,从而可准确获得突变点以及突变点对应的位置,从而可得到相应的雷达编码规律。本发明简单可靠,使得雷达号识别的准确率较高。
作为方法及装置的进一步改进,根据顺序记录的间隔,将顺序记录的间隔统一进行归一化的量化处理,得到对应的量化值,所述量化值为相应数量的码元;根据得到的相应数量的码元,识别雷达信号的编码规律。
作为方法及装置的进一步改进,所述码元为0或1。
作为方法及装置的进一步改进,为了方便识别编码,相邻间隔对应的码元不同。
作为方法及装置的进一步改进,所述将顺序记录的间隔统一进行归一化的量化处理,得到对应的量化值包括:选取所有间隔中的最小值,顺序求取每个间隔与间隔中的最小值的比值,并将得到的比值四舍五入,得到对应的量化值。
作为方法及装置的进一步改进,通过观察滤波后的雷达脉冲信号的形状,并根据选择的检测门限,来确定滤波后的雷达脉冲信号的突变点以及各个突变点对应的位置。
附图说明
图1是一种雷达信号相位编码规律识别方法流程图;
图2是接收采集后的数字信号示意图;
图3是信号频谱示意图;
图4是信号滤波输出以及门限检测的示意图;
图5是一种编码规律识别结果图;
图6是另一种编码规律识别结果图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚,下面结合附图及实施例,对本发明作进一步的详细说明。
本发明提供了一种雷达信号相位编码规律识别装置,该装置包括处理器,该处理器用于执行指令实现本发明的雷达信号相位编码规律识别方法。下面结合图1,对该方法做详细的说明。
首先,采集雷达脉冲信号。假设接收到的雷达脉冲信号如图2所示,为s(t):
s(t)=a(t)exp{j2πf0t}+w(t),0≤t≤T
上式中,a(t)表示雷达辐射源信号的基带信号,f0表示信号的载频,w(t)表示接收机噪声信号,T表示脉冲信号的时长。
假设采样的时间间隔为Ts,采样的总点数为N,则得到的数字信号为:
s(n)=a(nTs)exp{j2πf0nTs}+w(nTs),0≤n≤N
然后,生成小波滤波器,小波滤波器的表达式为:
上式中,Δ为带宽参数,δt为取样间隔,a为尺度参数,M为滤波器的长度,f为频率参数。
在实际的场景中,可根据实际情况选择合适的δt、Δ、a、f和M,以达到最优的处理结果。一种典型的δt、Δ、a、f和M的取值如下:
接着,在生成上述介绍的小波滤波器后,对采集到的雷达脉冲信号进行滤波处理。滤波系数由小波滤波器系数x(n)决定。该滤波过程,可以在时域实现,也可以在频域实现。
在时域实现的表达式为:
上式中,y(n)为小波滤波器的输出。
在频域实现的表达式为:
y(n)=DFTL -1[DFTL(s(n))⊙DFTL(x(n))]
上式中,DFTL(s(n))表示对信号s(n)进行L点的离散傅里叶变换(DFT),DFT-1 L(s(n))表示对信号s(n)进行L点的离散傅里叶逆变换,⊙表示对应位置相乘,L=max(M,N)表示M、N中的较大值。
其次,对y(n)的突变点进行检测。通过观察|y(n)|形状,选择合适的检测门限,记录突变点的位置N1、N2、…、Np,其中,p为突变点的总个数。
最后,根据上述突变点,计算相邻两个突变点之间的间隔,分别为:L1=N2-N1,L2=N3-N2,…,Lp-1=Np-Np-1;按照顺序记录间隔,根据间隔来识别雷达脉冲信号的编码规律。具体步骤如下:
1)选取计算的各个间隔中的最小值Lmin,Lmin=min{L1,L2,…,Lp};
2)求取每个位置差Li(i=1,…,p-1)与最小值Lmin的比值,分别为:k1=L1/Lmin,k2=L2/Lmin,…,kp-1=Lp-1/Lmin;
3)将上述比值进行四舍五入,得到对应的量化值,m1=floor(k1),m2=floor(k2),…,mp-1=floor(kp-1),这里的量化值为相应数量的码元,码元为0或1;
4)根据得到的相应数量的码元,识别雷达信号的编码规律,且相邻两个间隔对应的码元不同。
其代码可编写如下:
m=floor(L1/Lmin),其中floor()表示四舍五入,则在G中连续写入m个1
For i=2:p-1
m=floor(Li/Lmin),其中floor()表示四舍五入。
若上一次写入的是0,则在G的结尾连续写入m个1;
若上一次写入的是1,则在G的结尾连续写入m个0;
END
例如,得到的量化值分别为5,4,1,2,那么对应的雷达信号的编码规律为111110000100,其对应的波形图可如图5所示,图5中上边的波形图为最终的编码规律图,下边的图为采样波形图,采用适当的采样周期便可将其识别出来;当然,也可将编码规律识别为000001111011。
在上述步骤1)~步骤4)中,其主体思想为:根据顺序记录的间隔,将顺序记录的间隔统一进行归一化的量化处理,来得到对应的量化值,这里的量化值即为相应数量的码元。上述步骤中,通过求取所有间隔与间隔的最小值的比值来进行归一化处理,当然,还可采用其他的归一化的量化处理方法来得到对应的量化值。
而且,在进行编码时,为了方便识别,设置的相邻两个间隔对应的码元不同。当然,也可设置的相同,此时需要人为的在各个间隔对应的码元中间添加与之不同的码元。例如,得到的量化值分别为5,4,1,2,那么对应的雷达信号的编码规律可为111110011110010011,在这里,人为地添加了两个0来表示各个间隔之间的空隙,添加0的个数可根据需求设定,对应的波形图可如图6所示;当然,也可识别为00000110000110110,在该识别规律里,间隙设置为两个1。
例如,若进行滤波处理后的波形如图3所示。
在图3中检测门限为y=2800,故可确定突变点分别为:N1=0,N2=900,N3=1460,N4=2000四个突变点。则计算的各个突变点之间的位置差分别为:L1=900,L2=560,L3=540,则Lmin=540。计算的各个比值分别为k1=900/540=1.667,k2=560/540=1.037,k3=540/540=1,将这三个比值四舍五入,得到的编码个数数分别为:m1=2,m1=1,m1=1。则编码规律为:0010或者1101。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (10)
1.一种雷达信号相位编码规律识别方法,其特征在于,包括如下步骤:
采集雷达脉冲信号;
对采集的雷达脉冲信号进行小波滤波处理,得到滤波后的雷达脉冲信号;
确定滤波后的雷达脉冲信号的突变点以及各个突变点对应的位置;
计算相邻两个突变点之间的间隔,顺序记录间隔;
根据顺序记录的间隔,识别雷达信号的编码规律。
2.根据权利要求1所述的雷达信号相位编码规律识别方法,其特征在于,根据顺序记录的间隔,将顺序记录的间隔统一进行归一化的量化处理,得到对应的量化值,所述量化值为相应数量的码元;根据得到的相应数量的码元,识别雷达信号的编码规律。
3.根据权利要求2所述的雷达信号相位编码规律识别方法,其特征在于,所述码元为0或1。
4.根据权利要求3所述的雷达信号相位编码规律识别方法,其特征在于,相邻间隔对应的码元不同。
5.根据权利要求2所述的雷达信号相位编码规律识别方法,其特征在于,所述将顺序记录的间隔统一进行归一化的量化处理,得到对应的量化值包括:选取所有间隔中的最小值,顺序求取每个间隔与间隔中的最小值的比值,并将得到的比值四舍五入,得到对应的量化值。
6.根据权利要求1所述的雷达信号相位编码规律识别方法,其特征在于,通过观察滤波后的雷达脉冲信号的形状,并根据选择的检测门限,来确定滤波后的雷达脉冲信号的突变点以及各个突变点对应的位置。
7.一种雷达信号相位编码规律识别装置,其特征在于,包括处理器,所述处理器用于执行指令实现如下方法:
采集雷达脉冲信号;
对采集的雷达脉冲信号进行小波滤波处理,得到滤波后的雷达脉冲信号;
确定滤波后的雷达脉冲信号的突变点以及各个突变点对应的位置;
计算相邻两个突变点之间的间隔,顺序记录间隔;
根据顺序记录的间隔,识别雷达信号的编码规律。
8.根据权利要求7所述的雷达信号相位编码规律识别装置,其特征在于,根据顺序记录的间隔,将顺序记录的间隔统一进行归一化的量化处理,得到对应的量化值,所述量化值为相应数量的码元;根据得到的相应数量的码元,识别雷达信号的编码规律。
9.根据权利要求8所述的雷达信号相位编码规律识别装置,其特征在于,所述码元为0或1。
10.根据权利要求9所述的雷达信号相位编码规律识别装置,其特征在于,相邻间隔对应的码元不同。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810615792.1A CN108761441A (zh) | 2018-06-14 | 2018-06-14 | 一种雷达信号相位编码规律识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810615792.1A CN108761441A (zh) | 2018-06-14 | 2018-06-14 | 一种雷达信号相位编码规律识别方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108761441A true CN108761441A (zh) | 2018-11-06 |
Family
ID=64021146
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810615792.1A Pending CN108761441A (zh) | 2018-06-14 | 2018-06-14 | 一种雷达信号相位编码规律识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108761441A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112213693A (zh) * | 2020-09-21 | 2021-01-12 | 扬州健行电子科技有限公司 | 相位编码雷达信号码元宽度的解算方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105277923A (zh) * | 2015-11-20 | 2016-01-27 | 清华大学 | 一种单通道雷达信号分选方法 |
CN106908773A (zh) * | 2017-03-31 | 2017-06-30 | 中国科学院电子学研究所 | 连续脉冲雷达的信号时域循环编码和恢复方法 |
-
2018
- 2018-06-14 CN CN201810615792.1A patent/CN108761441A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105277923A (zh) * | 2015-11-20 | 2016-01-27 | 清华大学 | 一种单通道雷达信号分选方法 |
CN106908773A (zh) * | 2017-03-31 | 2017-06-30 | 中国科学院电子学研究所 | 连续脉冲雷达的信号时域循环编码和恢复方法 |
Non-Patent Citations (2)
Title |
---|
张鑫 等: "基于小波变换模值的二相编码信号调制特征分析", 《通信技术》 * |
陈兵 等: "一种快速的多二相编码信号参数估计方法", 《电子信息对抗技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112213693A (zh) * | 2020-09-21 | 2021-01-12 | 扬州健行电子科技有限公司 | 相位编码雷达信号码元宽度的解算方法及系统 |
CN112213693B (zh) * | 2020-09-21 | 2024-02-06 | 扬州健行电子科技有限公司 | 相位编码雷达信号码元宽度的解算方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109034127B (zh) | 一种频谱异常检测方法、装置和电子设备 | |
US4559602A (en) | Signal processing and synthesizing method and apparatus | |
US6832172B2 (en) | Apparatus and method for spectrum analysis-based serial data jitter measurement | |
CN108345033A (zh) | 一种微地震信号时频域初至检测方法 | |
CN110503060B (zh) | 一种光谱信号去噪方法及其系统 | |
WO2016060620A1 (en) | Frame based spike detection module | |
CN105429719A (zh) | 基于功率谱和多尺度小波变换分析强干扰信号检测方法 | |
CN115031875B (zh) | 基于频谱反馈的振弦传感器自适应扫频激励方法及系统 | |
CN110647656A (zh) | 一种利用变换域稀疏化和压缩降维的音频检索方法 | |
CN112364296A (zh) | 一种基于深度学习的p波到时自动拾取方法 | |
CN114760176B (zh) | 电力线通信自适应脉冲噪声抑制方法和装置、存储介质 | |
CN108761441A (zh) | 一种雷达信号相位编码规律识别方法及装置 | |
CN102988041A (zh) | 心磁信号噪声抑制中的信号选择性平均方法 | |
CN109450829B (zh) | 数字调制信号码速率估计方法和装置 | |
CN112509601B (zh) | 一种音符起始点检测方法及系统 | |
CN110632563B (zh) | 一种基于短时傅里叶变换的脉内频率编码信号参数测量方法 | |
CN103995799B (zh) | 基于fft谱校正的频率相位脑机接口解码方法及其装置 | |
CN110287853B (zh) | 一种基于小波分解的暂态信号去噪方法 | |
CN101849823B (zh) | 一种基于排列组合熵的神经元动作电位特征提取方法 | |
CN110944336A (zh) | 一种基于有限新息率的时频谱感知方法 | |
CN106802293A (zh) | 波形峰值检测方法及装置 | |
CN108596383B (zh) | 储层分类的方法及装置 | |
CN112488092B (zh) | 基于深度神经网络的导航频段信号类型识别方法及系统 | |
CN114239757A (zh) | 一种电磁时间序列数据的去噪方法及系统 | |
CN111008356B (zh) | 一种基于WTSVD算法扣除背景的γ能谱集分析方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181106 |