CN116819185A - 一种复杂电磁环境下基于信道化的数字鉴相方法及系统 - Google Patents
一种复杂电磁环境下基于信道化的数字鉴相方法及系统 Download PDFInfo
- Publication number
- CN116819185A CN116819185A CN202311093974.4A CN202311093974A CN116819185A CN 116819185 A CN116819185 A CN 116819185A CN 202311093974 A CN202311093974 A CN 202311093974A CN 116819185 A CN116819185 A CN 116819185A
- Authority
- CN
- China
- Prior art keywords
- pulse
- sub
- intermediate frequency
- frequency
- information
- 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
- 238000012850 discrimination method Methods 0.000 title claims description 11
- 238000001514 detection method Methods 0.000 claims abstract description 47
- 238000005070 sampling Methods 0.000 claims abstract description 43
- 238000005259 measurement Methods 0.000 claims abstract description 40
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000003775 Density Functional Theory Methods 0.000 claims abstract description 14
- 238000012544 monitoring process Methods 0.000 claims abstract description 11
- 230000004927 fusion Effects 0.000 claims description 16
- 238000001914 filtration Methods 0.000 claims description 10
- 230000003044 adaptive effect Effects 0.000 claims description 9
- 238000003708 edge detection Methods 0.000 claims description 6
- 230000003139 buffering effect Effects 0.000 claims description 5
- 230000010355 oscillation Effects 0.000 claims description 5
- 238000001228 spectrum Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000000926 separation method Methods 0.000 abstract 1
- 230000035945 sensitivity Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 240000007711 Peperomia pellucida Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0864—Measuring electromagnetic field characteristics characterised by constructional or functional features
- G01R29/0892—Details related to signal analysis or treatment; presenting results, e.g. displays; measuring specific signal features other than field strength, e.g. polarisation, field modes, phase, envelope, maximum value
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/165—Spectrum analysis; Fourier analysis using filters
- G01R23/167—Spectrum analysis; Fourier analysis using filters with digital filters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R25/00—Arrangements for measuring phase angle between a voltage and a current or between voltages or currents
-
- 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
-
- 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
Abstract
本发明公开了一种复杂电磁环境下基于信道化的数字鉴相方法及系统,属于雷达电磁环境监测技术领域。包括以下步骤:S1、对多路中频信号进行采样并缓存;对其中一路的中频信号进行数字信道化处理;S2、得到各子信道前、后沿信息;S3、编码产生前沿测量引导信息;S4、获取对应时刻开始缓存的多路中频数据;S5、进行DFT计算,得到鉴相结果;S6、对于复杂电磁环境下时域交叠信号,信道化检测和编码时同时产生多个前沿测量引导信息,产生的多个前沿引导信息重复S4和S5直至所有脉冲信号相位测完。本发明通过数字信道化将信号滤波到对应的子信道,实现将复杂电磁环境下的多输入信号从频域上分离。
Description
技术领域
本发明属于雷达电磁环境监测技术领域,具体涉及一种复杂电磁环境下基于信道化的数字鉴相方法及系统。
背景技术
如今电磁环境中包含的电磁辐射源种类和数量不断增加,如何在复杂电磁环境下准确有效地对各关注辐射源进行监测,对电磁环境监测提出了挑战。
现有的雷达电磁环境监测系统接收机普遍采用超外差结构,现有的鉴相系统大多基于时域检波包络引导,导致天线波束范围内时域交叠但空域频域分开信号会因检波包络出现时域混叠仅能测到功率较大信号的相位,丢失其他功率较小的信号。
发明内容
发明目的:为了解决上述问题,本发明提供了一种复杂电磁环境下基于信道化的数字鉴相方法及系统。
技术方案:一种复杂电磁环境下基于信道化的数字鉴相方法,包括以下步骤:
S1、对多路中频信号进行采样并缓存,得到多路的中频数据;对其中一路的中频信
号进行数字信道化处理,划分出K个子信道;
S2、对K个子信道同时进行检测,得到各子信道中的脉冲视频包络;对各子信道的
脉冲检波包络分别进行测量得到对应的前、后沿信息;其中,所述各子信道的前、后沿信息
至少包括:脉冲前、后沿到达时间、频率、幅度,其中0≤k≤K;
S3、对K个子信道的前、后沿信息分别进行编码融合,得到融合结果;根据所述融合
结果产生前沿测量引导信息和后沿脉冲描述字;其中,前沿测量引导信息中包含
以下信息:脉冲到达时间、前沿频率,其中0≤i≤500000;
S4、基于脉冲到达时间,获取对应时刻缓存的多路的中频数据;
S5、基于前沿频率对所述中频数据进行DFT计算,得到鉴相结果;其中,所
述鉴相结果包括:相位和数字幅度;
S6、对于复杂电磁环境下的时域交叠信号,在信道化检测和编码时同时产生多个
前沿测量引导信息,多个前沿测量引导信息数量为4~16个;将产生的前沿引导信息重
复S4和S5直至所有脉冲信号相位测完;
S7、基于根据脉冲到达时间、前沿频率,将每个鉴相结果与其对应的后
沿脉冲描述字进行融合。
进一步地,S1中进行数字信道化处理具体包括以下步骤:
S101、设置低通FIR原型滤波器,系数为,n=0,1,…,N-1;其中,N为滤波器阶
数;设置滤波器阶数N与信道数K满足倍数关系;
对滤波器的通带宽度进行扩展,滤波器的通带宽度计算公式如下:
;其中为采样时钟频率,为扩展带宽;
S102、将低通FIR原型滤波器累计进行K次移频,其中第k次移频,0≤k≤K,对
应得到K个滤波器;将K个滤波器组成滤波器组H,所述
滤波器组H与其对称频率的滤波器组覆盖完整的全部频谱;
S103、基于频率对称性确定第二奈奎斯特区的滤波器组中各个滤波器的
编号为:
;
S104、将多路中频信号分别与本振信号进行混频,变换到第二奈奎斯特区进行中
频采样,采样信号信号经过滤波器组H中的各个滤波器分别对应输出得到:
;
其中,为输入第n个点采样信号;为采样信号信号经过滤波器组H
中第k个滤波器输出的信号;为滤波器组H中第k个滤波器的系数;为H中
通带0滤波器系数的倒序;j为虚数单位;N为滤波器阶数。
进一步地,S2中得到各子信道的脉冲视频包络,前、后沿信息采用自适应门限检波方法,自适应门限检波方法至少包括以下步骤:
S201、获取K个子信道的滤波数据,并基于K个子信道的滤波数据对应计算各子信道的脉冲幅度,计算公式如下:
;
其中,表示第m个子信道的脉冲幅度;和分别为第m个滤波器
输出的实部和虚部;
S202、对各子信道的脉冲幅度进行检波,并测量实时噪底;基于噪底和脉冲幅度计算各子信道的浮动门限;比较每个子信道的脉冲幅度与各自浮动门限并得到脉冲视频包络VP;
S203、采用高精度时钟对各子信道在脉冲视频包络VP的前沿、后沿进行测量并得
到脉冲到达时间、前沿频率。
进一步地,S3中采用跟踪编码方式对K个子信道的前、后沿信息分别进行编码融合,包括如下步骤:
S301、根据各子信道的脉冲到达时间建立编码器,基于融合准则对各子信道的前、后沿信息进行融合;
所述融合准则为:若任意两个子信道的脉冲检波包络到达时间相差为-50~100ns、前沿频率相差为-150~150kHz、以及两个子信道顺序相邻,则表示两个子信道所对应的脉冲检波包络属于同一个信号,则对两个子信道的前、后沿信息进行融合;反之,则表示两个子信道所对应的脉冲检波包络为不同的信号,基于不同的信号各自建立对应的编码器;
S302、当其中一个编码器收到第一个前沿信息并等待预定时间后没有新的前沿信
息到达并更新编码准则时,产生前沿测量引导信息;
当一个编码器收到后沿信息并结束编码时产生后沿脉冲描述字。
进一步地,S4中采用内部RAM缓存多路的中频数据,基于脉冲到达时间,获取对
应时刻缓存的多路的中频数据包括如下步骤:
S401、在完成多路中频信号采样后,对多路的中频数据进行实时缓存,缓存的空间
以数字信道化检测产生前沿测量引导信息的最大延时为下限;
S402、多路的中频数据在缓存时,RAM使用计时器Timer的计时作为存储地址;
S403、多路的中频数据在缓存时,RAM使用前沿引导信息中的脉冲到达时间信息减去固定延时作为读取起始地址,读取预定长度的中频数据,进行鉴相。
进一步地,S5中鉴相结果的获得至少包括如下步骤:
S501、根据中的前沿频率对多路的中频数据进行对应频点的DFT计算;公式如
下:
;
其中,M表示DFT点数,;为输入第n个点采样信号;
为信号对应频点的旋转因子:;
S502、对DFT计算得到的虚部和实部作反正切得到信号的相位;
;
其中,分别表示DFT运算结果的虚部和实部。
进一步地,当多个信号在时域相互交叠频域分离时,对每个信号前沿均完成数字鉴相。
在另一个技术方案中,提供了一种复杂电磁环境下基于信道化的数字鉴相系统,用于实现如上述的一种复杂电磁环境下基于信道化的数字鉴相方法,所述系统包括:
参数控制模块,用于设置工作参数;所述工作参数至少包括本振频点、校准时间、校准源频点、工作时长;
时钟管理模块,用于产生采样时钟、系统工作时钟,进行时钟频率监控及系统复位;
中频采样模块,用于完成多路ADC进行配置和同步,对多路中频信号进行同步采样;
校准模块,用于控制校准源频点产生固定频点和脉宽脉冲信号,并将校准源脉冲的鉴相结果进行记录和标记形成校准表供测向系统使用;
信道化滤波模块,用于生成信道化滤波器组并对输入中频进行多相滤波,将信号划分到相应的信道内;
信道化检测模块,用于对各子信道滤波结果进行幅度检波;
跟踪编码和测量引导模块,用于根据各子信道前、后沿检测信息,并以到达时间、频率、幅度、以及信道编号完成后沿脉冲描述字编码,同时产生前沿测量开始引导和后沿测量结束引导信息;
中频缓存模块,用于以计时器作为写地址存入中频,并根据测量引导信息中的信号到达时间读取固定长度中频信号;
数字鉴相模块,用于根据输入引导频率对中频做定点DFT和反正切运算获得多路中的中频的相位。
有益效果:
(1)对中频采样时钟进行生成和监控并产生相应的同步复位;对多路ADC中频进行同步;对多路中频进行信号采样和缓存,进行数字信道化将信号划分成若干子信道,对信号进行数字检波和编码形成信号前沿描述字和视频检波包络;根据信号前沿的到达时间,频率等信息进行引导读取相应时刻的中频数据;对中频数据做对应频点的DFT得到相应的相位结果与信道化检测结果;
(2)本发明对复杂电磁环境适应能力强,数字信道化把复杂电磁环境下时域重叠的信号在频域上分离,对分离的信号进行前沿固定时长鉴相,实现在脉冲持续期间对多个信号完成鉴相,极大增加系统对时域重叠多信号的测向能力;而且数字信道化将信号滤波到对应的子信道并从频域上分离,能够有效应对不同空域同时到达接收机的复杂信号;
(3)本方法灵敏度高,由于子信道带宽一般在MHz量级,相比于宽开设备工作带宽GHz量级,灵敏度更高;鉴相精度高,信道化的精测频精度能达到KHz级,引导鉴相能够有效地提高鉴相精度至3°~5°。
附图说明
图1是基于本发明的数字鉴相系统的原理框图。
图2是本发明的数字信道化组成框图。
图3是原型滤波器理论幅频响应图。
图4是原型滤波器实际幅频响应图。
图5是第1个子信道滤波器的幅频响应图。
图6是滤波器组的幅频响应图。
图7是多相滤波计算结构实现框图。
具体实施方式
实施例1
如图1至7所示,本实施例提供了一种复杂电磁环境下基于信道化的数字鉴相方法,其特征在于,包括以下步骤:
S1、对多路中频信号进行采样并缓存,得到多路的中频数据;对其中一路的中频信
号进行数字信道化处理,划分出K个子信道。
以下以4路中频信号来举例说明:对中频采样时钟进行生成和监控并产生相应的同步复位,并对4路模数转换器件(ADC)进行中频采样同步;对4路中频信号进行采样,并用计时器作为内部RAM写地址进行缓存。对第1路中频信号进行数字信道化处理。
进一步地,S 1中进行数字信道化处理,具体包括以下步骤:
S101、设置低通FIR原型滤波器,系数为,n=0,1,…,N-1;其中,N为滤波器阶
数;设置滤波器阶数N与信道数K满足倍数关系;
对滤波器的通带宽度进行扩展,滤波器的通带宽度计算公式如下:
;其中,为扩展带宽;扩展的宽度与系统适应的最小脉宽有
关,通常大于;表示最小脉宽。
S102、将低通FIR原型滤波器累计进行K次移频,其中第 k次移频,0≤k≤K,
对应得到K个滤波器;将K个滤波器组成滤波器组H,所
述滤波器组H与其对称频率的滤波器组覆盖完整的全部频谱。
换言之,根据实信号频谱具有对称性,以上滤波器组及其对称频率的滤波器组通道频率刚好可以覆盖完整的0~Fs全部频谱。
S103、基于频率对称性确定第二奈奎斯特区的滤波器组中各个滤波器的
编号为:
。
S104、将多路中频信号分别与本振信号进行混频,变换到第二奈奎斯特区进行中
频采样,采样信号信号经过滤波器组中的各个滤波器分别对应输出得到:
;
其中,为输入第n个点采样信号;为采样信号信号经过滤波器组H
中第k个滤波器输出的信号;为滤波器组H中第k个滤波器的系数;为H中
通带0滤波器系数的倒序;j为虚数单位;N为滤波器阶数。
如图 1的系统结构,中频信号混频后变换到第二奈奎斯特区中频,由实信号的对称性得出第二奈奎斯特区滤波器组编号。
S2、对K个子信道同时进行检测,得到各子信道中的脉冲视频包络;对各子信道的
脉冲检波包络分别进行测量得到对应的前、后沿信息;其中,所述各子信道的前、后沿信息
至少包括:脉冲前、后沿到达时间、频率、幅度,其中0≤k≤K。
对K个子信道同时进行信号检测,得到各子信道的信号包络,并在前后沿测量参数得到前、后沿信息。
进一步地,S2中得到各子信道的脉冲视频包络,前、后沿信息采用自适应门限检波方法,自适应门限检波方法至少包括以下步骤:
S201、获取K个子信道的滤波数据,并基于K个子信道的滤波数据对应计算各子信道的脉冲幅度,计算公式如下:
;
其中,表示第m个子信道的脉冲幅度;和分别为第m个滤波器
输出的实部和虚部。
S202、对各子信道的脉冲幅度进行检波,并测量实时噪底;基于噪底和脉冲幅度计算各子信道的浮动门限;比较每个子信道的脉冲幅度与各自浮动门限并得到脉冲视频包络VP。
S203、采用高精度时钟对各子信道在脉冲视频包络VP的前沿、后沿进行测量并得
到脉冲到达时间、前沿频率。
所述自适应门限检波方法,门限生成需实时统计每个子信道一段时间内的极小值和极大值,极值统计采用滑动窗口统计法,窗口的长度需大于系统适应的最大脉宽,通常为2.5ms。自适应门限为极小值放大n倍和极大值缩小m倍两者中的大者。自适应门限可以满足系统的灵敏度和脉宽测量精度要求。
S3、对K个子信道的前、后沿信息分别进行编码融合,得到融合结果;根据所述融合
结果产生前沿测量引导信息和后沿脉冲描述字;其中,前沿测量引导信息中包含
以下信息:脉冲到达时间、前沿频率,其中0≤i≤500000。
进一步地,S3中采用跟踪编码方式对K个子信道的前、后沿信息分别进行编码融合,包括如下步骤:
S301、根据各子信道的脉冲到达时间建立编码器,基于融合准则对各子信道的前、后沿信息进行融合;
所述融合准则为:若任意两个子信道的脉冲检波包络到达时间相差为-50~100ns、前沿频率相差为-150~150kHz、以及两个子信道顺序相邻,则表示两个子信道所对应的脉冲检波包络属于同一个信号,则对两个子信道的前、后沿信息进行融合;反之,则表示两个子信道所对应的脉冲检波包络为不同的信号,基于不同的信号各自建立对应的编码器。
S302、当其中一个编码器收到第一个前沿信息并等待预定时间后没有新的前沿信
息到达并更新编码准则时,产生前沿测量引导信息;
当一个编码器收到后沿信息并结束编码时产生后沿脉冲描述字。
对于复杂电磁环境下可能存在时域交叠的情况,即和时间上交叠。
所述使用跟踪编码对各子信道检波信息融合,编码器的数量为系统适应的最大同时到达信号数,一般为4~8,在系统硬件资源足够的情况下可以适当地增加。
S4、基于脉冲到达时间,获取对应时刻缓存的多路的中频数据。
根据测量引导信息Gi中的到达时间Ti读取对应时刻开始缓存的4路中频数据。
进一步地, S4中采用内部RAM缓存多路的中频数据,基于脉冲到达时间,获取对
应时刻缓存的多路的中频数据包括如下步骤:
S401、在完成多路中频信号采样后,对多路的中频数据进行实时缓存,缓存的空间
以数字信道化检测产生前沿测量引导信息的最大延时为下限。
S402、多路的中频数据在缓存时,RAM使用计时器Timer的计时作为存储地址。
S403、多路的中频数据在缓存时,RAM使用前沿引导信息中的脉冲到达时间信息减去固定延时作为读取起始地址,读取预定长度的中频数据,进行鉴相。
S5、基于前沿频率对所述中频数据进行DFT计算,得到鉴相结果;其中,所
述鉴相结果包括:相位和数字幅度。
根据测量引导信息Gi中的前沿频率Fi对中频数据计算DFT,得到相位和幅度信息。
进一步地,S5中鉴相结果的获得至少包括如下步骤:
S501、根据中的前沿频率对多路的中频数据进行对应频点的DFT计算;公式如
下:
;
其中,M表示DFT点数,;为输入第n个点采样信号;为信号对应频
点的旋转因子:;j为虚数单位。
S502、对DFT计算得到的虚部和实部作反正切得到信号的相位;
;
其中,分别表示DFT运算结果的虚部和实部。
S6、对于复杂电磁环境下的时域交叠信号,在信道化检测和编码时同时产生多个
前沿测量引导信息,多个前沿测量引导信息数量为4~16个;将产生的前沿引导信息重
复S4和S5直至所有脉冲信号相位测完。
S7、基于根据脉冲到达时间、前沿频率,将每个鉴相结果与其对应的后
沿脉冲描述字进行融合。
实施例2
本实施例提供了一种复杂电磁环境下基于信道化的数字鉴相系统,系统包括:
参数控制模块,用于设置工作参数;所述工作参数至少包括频综频点、校准时间、校准源频点、工作时长;
时钟管理模块,用于产生采样时钟、系统工作时钟,进行时钟频率监控及系统复位;
中频采样模块,用于完成多路ADC进行配置和同步,对多路中频信号进行同步采样;
校准模块,用于控制校准源频点产生固定频点和脉宽脉冲信号,并将校准源脉冲的鉴相结果进行记录和标记形成校准表供测向系统使用;
信道化滤波模块,用于生成信道化滤波器组并对输入中频进行多相滤波,将信号划分到相应的信道内;
信道化检测模块,用于对各子信道滤波结果进行幅度检波,功能具有:浮动门限生成、脉冲检波、兔耳抑制、连续波检波、频率精测;
跟踪编码和测量引导模块,用于根据各子信道前、后沿检测信息,并以到达时间、频率、幅度、以及信道编号完成后沿脉冲描述字编码,同时产生前沿测量开始引导和后沿测量结束引导信息;该模块能够适应常规、调频、调相等信号;
中频缓存模块,用于以计时器作为写地址存入中频,并根据测量引导信息中的信号到达时间读取固定长度中频信号;
数字鉴相模块,用于根据输入引导频率对中频做定点DFT和反正切运算获得多路中的中频的相位。
该系统的工作流程如下:
参数控制模块设置工作参数,在工作周期开始前控制接收开关在自检或校准模式下选择天线信号和校准信号,控制频综频点将中频信号混频至ADC可采样的中频频域。
时钟管理模块根据ADC随路时钟利用锁相环产生适合FPGA运行的采样时钟,并对时钟频率实时监控及系统工作复位。
中频采样模块在上电时对4路ADC进行配置和同步复位,后根据采样时钟和ADC接口协议对4路中频进行同步采样并以计时器作为写地址存入中频缓存RAM;选取第1路中频输入信道化多相滤波器组进行滤波,将信号划分到各子信道中。
信道化检波模块对各子信道数据计算幅度和动态自适应门限,完成视频包络检波,对检波包络前后沿测量到达时间,频率和幅度形成前后沿描述字。
编码模块对前后沿描述字以到达时间,频率,幅度,信道编号等信息完成后沿脉冲描述字编码,同时产生前沿测量开始引导和后沿测量结束引导信息并存入FIFO;依次读取前沿测量引导信息根据测量引导信息中的信号到达时间读取中频缓存模块中固定长度中频信号。
数字鉴相模块根据输入引导频率对引导中频做定点DFT和反正切运算获得4路中频的相位直至所有引导信息处理完成;根据到达时间和频率将测得的相位与编码产生的PDW进入融合。
实施例3
在实施例2的基础之上,结合附图1~6,以一种宽带环境监测系统(详见图1)为例,对本发明实施例中的技术方案进行清楚、完整的描述。
(1):参照图1所示的一种宽带环境监测系统,工作频段为1~18GHz,采用1GHz瞬时中频带宽,通过频综可划分成步进为500MHz的33个频段,中频混频至1.4~2.4GHz。
(2):ADC采样率为2.5GHz采样第二奈奎斯特区中频信号;数字信道化使用FPGA实现,FPGA运行时钟选用312.5MHz,则并行8路处理,信道数K设计为128,信道带宽19.53125MHz,信道通带为12.5MHz,可适应最小200ns脉宽脉冲信号。
(3):设置校正时间窗Wjz(毫秒级),侦收时间窗Wr(秒级),校正周期内校准源在瞬时带宽内按50MHz步进校正,每个点校正500us,对校正脉冲的处理与常规脉冲相同,仅在结果中打上校正频点标记。
(4):对4路中频进行最长200us时长缓存。
(5):对第1路中频进行信道化滤波,检测,编码,产生前沿测量引导信息。
(6):根据测量引导信息依次读取4路中频,进行数字鉴相直至所有信号前沿鉴相完成。
(7):将编码结果与数字鉴相结果进行融合输出至测向系统进行测向。
实施例4
基于上述实施例1,步骤S1数字信道化的设计(如图 2所示),步骤如下:
(1)射频混频选用ADC的第二奈奎斯特区1.4~2.4GHz,为综合考虑系统检测灵敏度和硬件资源划分出128个子信道,实际设计滤波器通带带宽要略大于理论计算选用12.5MHz。
(2)使用matlab工具设计一个低通FIR原型低通滤波器,系数,n=0,1,…,N-1;
半带幅频响应理论和实际设计如图 3、图 4所示。
(3)将原型滤波器频移累计K次得到一个均匀滤波器组,在第 k次中滤波器频移,0≤k≤K,根据实信号谱线对称性,易得均匀滤波器组的等效幅频响应如图 5、图 6所
示。
(4)由滤波器组推导多相滤波结构;
第k个滤波器的输出描述为:
其中为输入第n个点采样信号,N为滤波器阶数,上式中N=KL,K为信道数,L表
示常数,L取值范围为2*(6~8)。设计滤波器时考虑到硬件资源通常将N设置为信道数的6~8
倍即N=2K*(6~8),p=0,1,2……K-1,l=0,1,2,…,L-1。
由(3)可知:
滤波器具有对称性,则有:
令:
矩阵
矩阵
相应的变换域表示为,则多相滤波的计算结构如图 7所示。
实施例5
基于上述实施例1,步骤S5数字鉴相方法,步骤如下:
(1)根据前沿测量引导信息中的到达时间由中频缓存模块中读取脉冲前沿512
点中频数据。
(2)根据前沿测量引导信息中的信号频率Fi对4路中频做对应频点的DFT:
;
M表示DFT点数,即读取中频点数;;为输入第n个点采样信号;
为信号对应频点的旋转因子;
实现时通过DDS模块产生相应的系数。
(3)对DFT计算得到的实部和虚部作反正切得到信号的相位。
实现时通过CORDIC算法可同时得到相位和幅度。
(4)重复上述操作直至所有的信号前沿均完成鉴相。
Claims (8)
1.一种复杂电磁环境下基于信道化的数字鉴相方法,其特征在于,包括以下步骤:
S1、对多路中频信号进行采样并缓存,得到多路的中频数据;对其中一路的中频信号进行数字信道化处理,划分出K个子信道;
S2、对K个子信道同时进行检测,得到各子信道中的脉冲视频包络;对各子信道的脉冲检波包络分别进行测量得到对应的前、后沿信息;其中,所述各子信道的前、后沿信息至少包括:脉冲前、后沿到达时间、频率/>、幅度/>,其中0≤k≤K;
S3、对K个子信道的前、后沿信息分别进行编码融合,得到融合结果;根据所述融合结果产生前沿测量引导信息和后沿脉冲描述字/>;其中,前沿测量引导信息/>中包含以下信息:脉冲到达时间/>、前沿频率/>,其中0≤i≤500000;
S4、基于脉冲到达时间,获取对应时刻缓存的多路的中频数据;
S5、基于前沿频率对所述中频数据进行DFT计算,得到鉴相结果/>;其中,所述鉴相结果/>包括:相位和数字幅度;
S6、对于复杂电磁环境下的时域交叠信号,在信道化检测和编码时同时产生多个前沿测量引导信息,多个前沿测量引导信息/>数量为4~16个;将产生的前沿引导信息/>重复S4和S5直至所有脉冲信号相位测完;
S7、基于根据脉冲到达时间、前沿频率/>,将每个鉴相结果/>与其对应的后沿脉冲描述字/>进行融合。
2.如权利要求1所述的一种复杂电磁环境下基于信道化的数字鉴相方法,其特征在于,S1中进行数字信道化处理具体包括以下步骤:
S101、设置低通FIR原型滤波器,系数为/>,n=0,1,…,N-1;其中,N为滤波器阶数;设置滤波器阶数N与信道数K满足倍数关系;
对滤波器的通带宽度进行扩展,滤波器的通带宽度计算公式如下:
;其中/>为采样时钟频率,/>为扩展带宽;
S102、将低通FIR原型滤波器累计进行K次移频,其中第k次移频/>,0≤k≤K,对应得到K个滤波器/>;将K个滤波器/>组成滤波器组H,所述滤波器组H与其对称频率的滤波器组覆盖完整的/>全部频谱;
S103、基于频率对称性确定第二奈奎斯特区的滤波器组/>中各个滤波器的编号为:
;
S104、将多路中频信号分别与本振信号进行混频,变换到第二奈奎斯特区进行中频采样,采样信号信号经过滤波器组H中的各个滤波器分别对应输出得到/>:
;
其中,为输入第n个点采样信号;/>为采样信号/>信号经过滤波器组H中第k个滤波器输出的信号;/>为滤波器组H中第k个滤波器的系数;/>为H中通带0滤波器系数的倒序;j为虚数单位;N为滤波器阶数。
3.如权利要求1所述的一种复杂电磁环境下基于信道化的数字鉴相方法,其特征在于,
S2中得到各子信道的脉冲视频包络,前、后沿信息采用自适应门限检波方法,自适应门限检波方法至少包括以下步骤:
S201、获取K个子信道的滤波数据,并基于K个子信道的滤波数据对应计算各子信道的脉冲幅度,计算公式如下:
;
其中,表示第m个子信道的脉冲幅度;/>和/>分别为第m个滤波器输出的实部和虚部;
S202、对各子信道的脉冲幅度进行检波,并测量实时噪底;基于噪底和脉冲幅度计算各子信道的浮动门限;比较每个子信道的脉冲幅度与各自浮动门限并得到脉冲视频包络VP;
S203、采用高精度时钟对各子信道在脉冲视频包络VP的前沿、后沿进行测量并得到脉冲到达时间、前沿频率/>。
4.如权利要求3所述的一种复杂电磁环境下基于信道化的数字鉴相方法,其特征在于,
S3中采用跟踪编码方式对K个子信道的前、后沿信息分别进行编码融合,包括如下步骤:
S301、根据各子信道的脉冲到达时间建立编码器,基于融合准则对各子信道的前、后沿信息进行融合;
所述融合准则为:若任意两个子信道的脉冲检波包络到达时间相差为-50~100ns、前沿频率相差为-150~150kHz、以及两个子信道顺序相邻,则表示两个子信道所对应的脉冲检波包络属于同一个信号,则对两个子信道的前、后沿信息进行融合;反之,则表示两个子信道所对应的脉冲检波包络为不同的信号,基于不同的信号各自建立对应的编码器;
S302、当其中一个编码器收到第一个前沿信息并等待预定时间后没有新的前沿信息到达并更新编码准则时,产生前沿测量引导信息;
当一个编码器收到后沿信息并结束编码时产生后沿脉冲描述字。
5.如权利要求1所述的一种复杂电磁环境下基于信道化的数字鉴相方法,其特征在于,
S4中采用内部RAM缓存多路的中频数据,基于脉冲到达时间,获取对应时刻缓存的多路的中频数据包括如下步骤:
S401、在完成多路中频信号采样后,对多路的中频数据进行实时缓存,缓存的空间以数字信道化检测产生前沿测量引导信息的最大延时为下限;
S402、多路的中频数据在缓存时,RAM使用计时器Timer的计时作为存储地址;
S403、多路的中频数据在缓存时,RAM使用前沿引导信息中的脉冲到达时间信息/>减去固定延时作为读取起始地址,读取预定长度的中频数据,进行鉴相。
6.如权利要求1所述的一种复杂电磁环境下基于信道化的数字鉴相方法,其特征在于,S5中鉴相结果的获得至少包括如下步骤:
S501、根据中的前沿频率/>对多路的中频数据进行对应频点的DFT计算;公式如下:
;
其中,M表示DFT点数,;/>为输入第n个点采样信号;
为信号对应频点/>的旋转因子:/>;
S502、对DFT计算得到的虚部和实部作反正切得到信号的相位;
;
其中,分别表示DFT运算结果的虚部和实部。
7.如权利要求1所述的一种复杂电磁环境下基于信道化的数字鉴相方法,其特征在于,
当多个信号在时域相互交叠频域分离时,对每个信号前沿均完成数字鉴相。
8.一种复杂电磁环境下基于信道化的数字鉴相系统,其特征在于,用于实现如权利要求1至7任意一项所述的一种复杂电磁环境下基于信道化的数字鉴相方法,所述系统包括:
参数控制模块,用于设置工作参数;所述工作参数至少包括本振频点、校准时间、校准源频点、工作时长;
时钟管理模块,用于产生采样时钟、系统工作时钟,进行时钟频率监控及系统复位;
中频采样模块,用于完成多路ADC进行配置和同步,对多路中频信号进行同步采样;
校准模块,用于控制校准源频点产生固定频点和脉宽脉冲信号,并将校准源脉冲的鉴相结果进行记录和标记形成校准表供测向系统使用;
信道化滤波模块,用于生成信道化滤波器组并对输入中频进行多相滤波,将信号划分到相应的信道内;
信道化检测模块,用于对各子信道滤波结果进行幅度检波;
跟踪编码和测量引导模块,用于根据各子信道前、后沿检测信息,并以到达时间、频率、幅度、以及信道编号完成后沿脉冲描述字编码,同时产生前沿测量开始引导和后沿测量结束引导信息;
中频缓存模块,用于以计时器作为写地址存入中频,并根据测量引导信息中的信号到达时间读取固定长度中频信号;
数字鉴相模块,用于根据输入引导频率对中频做定点DFT和反正切运算获得多路中的中频的相位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311093974.4A CN116819185B (zh) | 2023-08-29 | 2023-08-29 | 一种复杂电磁环境下基于信道化的数字鉴相方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311093974.4A CN116819185B (zh) | 2023-08-29 | 2023-08-29 | 一种复杂电磁环境下基于信道化的数字鉴相方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116819185A true CN116819185A (zh) | 2023-09-29 |
CN116819185B CN116819185B (zh) | 2023-11-24 |
Family
ID=88122488
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311093974.4A Active CN116819185B (zh) | 2023-08-29 | 2023-08-29 | 一种复杂电磁环境下基于信道化的数字鉴相方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116819185B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117111016A (zh) * | 2023-10-20 | 2023-11-24 | 南京航天工业科技有限公司 | 复杂电磁环境下基于信道化的实时脉内分析方法及系统 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2497958A1 (fr) * | 1981-01-09 | 1982-07-16 | Thomson Csf | Dispositif de determination de l'instant d'arrivee d'impulsions, utilisation dans un equipement de mesure de distances et equipement de mesure comportant un tel dispositif |
US5406288A (en) * | 1992-06-30 | 1995-04-11 | Thomson-Csf | Method and device to improve the probability of validity of the codes of secondary radar responses |
WO1998029759A1 (en) * | 1997-01-02 | 1998-07-09 | Raytheon Company | Uniform discrete fourier transform filter parameter encoder |
US6043771A (en) * | 1998-12-21 | 2000-03-28 | Trw Inc. | Compact, sensitive, low power device for broadband radar detection |
US20050046609A1 (en) * | 2003-08-25 | 2005-03-03 | Richard Wasiewicz | Phased null radar |
CN104333737A (zh) * | 2014-10-28 | 2015-02-04 | 北京大学软件与微电子学院无锡产学研合作教育基地 | 近距离2-5兆带宽模拟无线图传系统及信号处理方法 |
CN106771668A (zh) * | 2017-01-05 | 2017-05-31 | 西南交通大学 | 一种电磁辐射参数测试系统 |
US20170153319A1 (en) * | 2015-11-30 | 2017-06-01 | Luminar Technologies, Inc. | Lidar system with distributed laser and multiple sensor heads |
CN110376559A (zh) * | 2019-09-05 | 2019-10-25 | 中国人民解放军空军工程大学 | 单通道雷达主瓣多源干扰分离方法、装置和设备 |
CN111427016A (zh) * | 2020-04-21 | 2020-07-17 | 北京航天长征飞行器研究所 | 基于虚拟信道化的多雷达对抗方法和系统 |
CN112859019A (zh) * | 2021-01-11 | 2021-05-28 | 北京无线电计量测试研究所 | 一种脉内调制类型参数提取系统和使用方法 |
CN113447893A (zh) * | 2021-09-01 | 2021-09-28 | 湖南艾科诺维科技有限公司 | 一种雷达脉冲信号频谱自动检测方法、系统及介质 |
CN115494497A (zh) * | 2022-09-30 | 2022-12-20 | 上海无线电设备研究所 | 一种相位编码脉冲雷达信号测向方法 |
CN115993580A (zh) * | 2022-12-19 | 2023-04-21 | 成都富元辰科技有限公司 | 一种被动侦察系统信号处理方法 |
-
2023
- 2023-08-29 CN CN202311093974.4A patent/CN116819185B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2497958A1 (fr) * | 1981-01-09 | 1982-07-16 | Thomson Csf | Dispositif de determination de l'instant d'arrivee d'impulsions, utilisation dans un equipement de mesure de distances et equipement de mesure comportant un tel dispositif |
US5406288A (en) * | 1992-06-30 | 1995-04-11 | Thomson-Csf | Method and device to improve the probability of validity of the codes of secondary radar responses |
WO1998029759A1 (en) * | 1997-01-02 | 1998-07-09 | Raytheon Company | Uniform discrete fourier transform filter parameter encoder |
US6043771A (en) * | 1998-12-21 | 2000-03-28 | Trw Inc. | Compact, sensitive, low power device for broadband radar detection |
US20050046609A1 (en) * | 2003-08-25 | 2005-03-03 | Richard Wasiewicz | Phased null radar |
CN104333737A (zh) * | 2014-10-28 | 2015-02-04 | 北京大学软件与微电子学院无锡产学研合作教育基地 | 近距离2-5兆带宽模拟无线图传系统及信号处理方法 |
US20170153319A1 (en) * | 2015-11-30 | 2017-06-01 | Luminar Technologies, Inc. | Lidar system with distributed laser and multiple sensor heads |
CN106771668A (zh) * | 2017-01-05 | 2017-05-31 | 西南交通大学 | 一种电磁辐射参数测试系统 |
CN110376559A (zh) * | 2019-09-05 | 2019-10-25 | 中国人民解放军空军工程大学 | 单通道雷达主瓣多源干扰分离方法、装置和设备 |
CN111427016A (zh) * | 2020-04-21 | 2020-07-17 | 北京航天长征飞行器研究所 | 基于虚拟信道化的多雷达对抗方法和系统 |
CN112859019A (zh) * | 2021-01-11 | 2021-05-28 | 北京无线电计量测试研究所 | 一种脉内调制类型参数提取系统和使用方法 |
CN113447893A (zh) * | 2021-09-01 | 2021-09-28 | 湖南艾科诺维科技有限公司 | 一种雷达脉冲信号频谱自动检测方法、系统及介质 |
CN115494497A (zh) * | 2022-09-30 | 2022-12-20 | 上海无线电设备研究所 | 一种相位编码脉冲雷达信号测向方法 |
CN115993580A (zh) * | 2022-12-19 | 2023-04-21 | 成都富元辰科技有限公司 | 一种被动侦察系统信号处理方法 |
Non-Patent Citations (2)
Title |
---|
HEMING WANG 等: "Design of Radar Electromagnetic Environment Simulator", 《UNIFYING ELECTRICAL ENGINEERING AND ELECTRONICS ENGINEERING》, pages 1435 * |
王平: "复杂电磁环境下动态背景信号生成系统的设计与研究", 《中国优秀硕士学位论文全文数据库 (信息科技辑)》, pages 136 - 364 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117111016A (zh) * | 2023-10-20 | 2023-11-24 | 南京航天工业科技有限公司 | 复杂电磁环境下基于信道化的实时脉内分析方法及系统 |
CN117111016B (zh) * | 2023-10-20 | 2023-12-19 | 南京航天工业科技有限公司 | 复杂电磁环境下基于信道化的实时脉内分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN116819185B (zh) | 2023-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116819185B (zh) | 一种复杂电磁环境下基于信道化的数字鉴相方法及系统 | |
CN104297738B (zh) | 多通道接收机同步校准装置及同步校准与误差补偿方法 | |
CN104316913B (zh) | 多通道接收机实时校准装置及校准与误差补偿方法 | |
Pfammatter et al. | A software-defined sensor architecture for large-scale wideband spectrum monitoring | |
EP0107050B1 (en) | Signal generator for digital spectrum analyzer | |
CN105656485B (zh) | 一种多通道时间交错adc测量校准方法和装置 | |
CN107329127B (zh) | 一种用于雷达系统dbf功能检测的相位线性分析方法及系统 | |
CN113965217B (zh) | 一种双通道、单通道s/c/x三波段宽带单比特数字测频接收机 | |
JP2012078187A (ja) | パラメータ検出器、レーダ装置、誘導装置、及びパラメータ検出方法 | |
US9450598B2 (en) | Two-stage digital down-conversion of RF pulses | |
KR101840828B1 (ko) | 연속파 레이더의 도플러 주파수를 검지하는 에스디알 수신기 및 그 동작 방법 | |
JP2006525503A (ja) | デジタル電子支援装置 | |
CN109655775B (zh) | 一种任意波形发生器幅度扫频多尺度校准方法和装置 | |
CN106199187B (zh) | 一种多音信号相对相位的测试方法 | |
CN102323482A (zh) | 数字中频频谱分析仪在网络分析测量中测量相频特性的方法 | |
EP1640734B1 (en) | A method and apparatus for measuring the input frequency response of a digital receiver | |
CN116436549A (zh) | 一种自适应抑制谐杂波信号的高精度频谱分析优化方法 | |
CN115361741A (zh) | 一种高精度通道信号延迟自动校准装置及方法 | |
Song et al. | Digital estimation and compensation of analog errors in frequency-interleaved ADCs | |
Szilvási et al. | Configurable hardware-based radio interferometric node localization | |
US9350470B1 (en) | Phase slope reference adapted for use in wideband phase spectrum measurements | |
KR101404085B1 (ko) | 무선 통신기, 그를 이용한 거리 측정 장치 및 방법 | |
RU110574U1 (ru) | Устройство экспресс-анализа радиоизлучений с псевдослучайной перестройкой рабочей частоты | |
RU2763583C2 (ru) | Способ определения частоты в приемнике с субдискретизацией | |
US5418540A (en) | Method and apparatus for eliminating signal sidelobes from a received signal frame in a multichannel receiver |
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 |