[go: up one dir, main page]

CN111257838A - Multichannel signal preprocessing method based on broadband receiver - Google Patents

Multichannel signal preprocessing method based on broadband receiver Download PDF

Info

Publication number
CN111257838A
CN111257838A CN202010042806.2A CN202010042806A CN111257838A CN 111257838 A CN111257838 A CN 111257838A CN 202010042806 A CN202010042806 A CN 202010042806A CN 111257838 A CN111257838 A CN 111257838A
Authority
CN
China
Prior art keywords
frequency
data
conversion
programmable gate
gate array
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
Application number
CN202010042806.2A
Other languages
Chinese (zh)
Other versions
CN111257838B (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN202010042806.2A priority Critical patent/CN111257838B/en
Publication of CN111257838A publication Critical patent/CN111257838A/en
Application granted granted Critical
Publication of CN111257838B publication Critical patent/CN111257838B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers
    • G01S7/292Extracting wanted echo-signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明属于雷达信号处理领域,具体涉及一种基于宽带接收机多通道信号预处理方法,满足对大带宽信号的采样处理亦能完成雷达信号的预处理,进而提高系统的抗干扰性能。

Figure 202010042806

The invention belongs to the field of radar signal processing, and in particular relates to a multi-channel signal preprocessing method based on a wideband receiver, which satisfies the sampling processing of large bandwidth signals and can also complete radar signal preprocessing, thereby improving the anti-jamming performance of the system.

Figure 202010042806

Description

一种基于宽带接收机多通道信号预处理方法A Multi-Channel Signal Preprocessing Method Based on Wideband Receiver

技术领域technical field

本发明属于雷达信号处理领域,具体涉及一种基于宽带接收机多通道信号预处理方法。The invention belongs to the field of radar signal processing, in particular to a multi-channel signal preprocessing method based on a wideband receiver.

背景技术Background technique

随着雷达技术及现代宽带通信技术的发展,系统对模拟输入带宽的要求越来越高,进而需要更高采样率的模数转换器(ADC)来满足这一需求,但高速的模数转换器(ADC)得到的大带宽信号需要后续处理器有极高的处理能力,对资源消耗极大,特别是多通道并行处理更是难以实现,针对这一现状可利用模数转换器(ADC)中内置的数字下变频模块(DDC)对原始信号进行数字下变频进而将数据率降低,得到的数据便可易于现场可编程门阵列(FPGA)做进一步的信号处理。With the development of radar technology and modern broadband communication technology, the system has higher and higher requirements for the analog input bandwidth, and then requires a higher sampling rate analog-to-digital converter (ADC) to meet this demand, but the high-speed analog-to-digital conversion The large bandwidth signal obtained by the ADC (ADC) requires the subsequent processor to have extremely high processing capabilities, which consumes a lot of resources, especially multi-channel parallel processing is difficult to achieve. The built-in digital down-conversion module (DDC) digitally down-converts the original signal to reduce the data rate, and the obtained data can be easily processed by a field programmable gate array (FPGA).

另外,由于雷达工作于越来越恶劣的电磁干扰环境中,所以该接收机会对数字信号进行预处理提高系统的抗干扰性能,进一步减轻信号处理机的负担,因此,研制高效的宽带数字接收机对于完成宽带通信接收系统的数字化改造,提高雷达,遥测等通信接收系统的性能,实现最终的软件无线电接收系统具有重要意义。In addition, since the radar works in an increasingly harsh electromagnetic interference environment, the receiver will preprocess the digital signal to improve the anti-interference performance of the system and further reduce the burden of the signal processor. Therefore, an efficient broadband digital receiver is developed. It is of great significance to complete the digital transformation of the broadband communication receiving system, improve the performance of the radar, telemetry and other communication receiving systems, and realize the final software radio receiving system.

发明内容SUMMARY OF THE INVENTION

为了解决现有技术中存在的上述问题,本发明提供了一种基于宽带接收机多通道信号预处理方法。本发明要解决的技术问题通过以下技术方案实现:In order to solve the above problems existing in the prior art, the present invention provides a multi-channel signal preprocessing method based on a wideband receiver. The technical problem to be solved by the present invention is realized by the following technical solutions:

一种基于宽带接收机多通道信号预处理方法,包括:A multi-channel signal preprocessing method based on a wideband receiver, comprising:

先确定回波信号的中频f0及调制带宽B,进而可得ADC的采样频率fs,需满足:

Figure BDA0002368341480000021
其中fH=f0+B/2,fL=f0-B/2,fh称为上限频率,fl称为下限频率;由此我们可以确定采样频率fs;由于需要利用模数转换器ADC的内置数字下变频模块,所以还需要确定抽取倍数及滤波器带宽,参数确定后,便可以将模数转换器ADC配置在相应的模式下进行工作;First determine the intermediate frequency f 0 and modulation bandwidth B of the echo signal, and then obtain the sampling frequency f s of the ADC, which must satisfy:
Figure BDA0002368341480000021
Where f H = f 0 +B/2, f L =f 0 -B/2, f h is called the upper limit frequency, and f l is called the lower limit frequency; thus we can determine the sampling frequency f s ; The built-in digital down-conversion module of the converter ADC, so it is necessary to determine the decimation multiple and filter bandwidth. After the parameters are determined, the analog-to-digital converter ADC can be configured to work in the corresponding mode;

多片的模数转换器ADC工作在同一参数下,将多片多通道模数转换器ADC的数据经204b协议编码后的数据帧送至现场可编程门阵列FPGA中进行解帧处理;经解帧后就可以得到第一级数字下变频后的数据;The multi-chip analog-to-digital converter ADC works under the same parameters, and the data frame encoded by the multi-chip multi-channel analog-to-digital converter ADC is sent to the field programmable gate array FPGA for de-frame processing; After the frame, the data after the first-level digital down-conversion can be obtained;

上述处理后数据已被混频到零频附近,且经过两级半带滤波后数据率大大降低,已适合现场可编程门阵列FPGA做信号的预处理,这一步会对其中某一通道的数据进行频谱分析,由现场可编程门阵列FPGA完成快速傅里叶变换FFT运算,得到的FFT计算结果会由光纤发送给信号处理机,由信号处理机进一步找出干扰最小的频点,该频点可用于下次混频的中频值,从而达到频率捷变,可引导雷达工作频率到干扰频谱的空隙或弱区;After the above processing, the data has been mixed to near zero frequency, and the data rate is greatly reduced after two-stage half-band filtering, which is suitable for field programmable gate array FPGA for signal preprocessing. This step will affect the data of one of the channels. Perform spectrum analysis, complete the fast Fourier transform FFT operation by the field programmable gate array FPGA, and the obtained FFT calculation result will be sent to the signal processor by the optical fiber, and the signal processor will further find the frequency point with the least interference. The intermediate frequency value that can be used for the next mixing, so as to achieve frequency agility, and can guide the radar operating frequency to the gap or weak area of the interference spectrum;

雷达工作在变频的模式中,需要完成第二级数字下变频DDC,将信号混频到基带,并进一步完成更大倍数的抽取;混频系数的中心频率由光纤接收的模式字中解析而来;根据频点值,现场可编程门阵列FPGA内利用DDS_IP核产生混频系数,与输入数据完成复乘;将与上一步对应产生的m个混频系数Δθ存储在ROM中,通过解析模式字完成频点选择;使用现场可编程门阵列FPGA中FIR_IP核完成滤波抽取,滤波器的系数可由MATLAB的FDA工具生成,同样将生成的系数存储于rom中,系统启动时该系数自动会被读到FIR_IP中,完成运算;The radar works in the frequency conversion mode, and needs to complete the second-stage digital down-conversion DDC, mix the signal to the baseband, and further complete the extraction of larger multiples; the center frequency of the mixing coefficient is analyzed from the mode word received by the fiber ; According to the frequency point value, the DDS_IP core is used in the field programmable gate array FPGA to generate the mixing coefficients, and the multiplication is completed with the input data; the m mixing coefficients Δθ corresponding to the previous step are stored in the ROM, and the mode word is analyzed by analyzing the mode word. Complete the frequency point selection; use the FIR_IP core in the field programmable gate array FPGA to complete the filter extraction, the coefficients of the filter can be generated by the FDA tool of MATLAB, and the generated coefficients are also stored in the rom, and the coefficients will be automatically read when the system starts. In FIR_IP, complete the operation;

外界强窄脉冲表现为在数字下变频DDC后有少数几个较大的值随着窄脉冲宽度与强度的不同,数字下变频DDC后扩展点的个数不同,会导致在脉压后形成一个台阶,出现虚警;因此需要在完成数字下变频DDC后进行窄脉冲剔除;The external strong and narrow pulses appear as a few larger values after the digital down-conversion DDC. With the difference in the width and intensity of the narrow pulses, the number of expansion points after the digital down-conversion DDC is different, which will lead to the formation of a pulse pressure after the DDC. Steps, false alarms occur; therefore, narrow pulse rejection needs to be performed after the completion of the digital down-conversion DDC;

对各种预处理结果通过光纤发送给信号处理机,同时光纤也会接受由信号处理机发送来的模式字,解析其中关于接收机的各部分工作参数。The various preprocessing results are sent to the signal processor through the optical fiber, and the optical fiber also accepts the mode word sent by the signal processor, and analyzes the working parameters of each part of the receiver.

附图说明Description of drawings

图1是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法流程示意图;1 is a schematic flowchart of a multi-channel signal preprocessing method based on a wideband receiver provided by an embodiment of the present invention;

图2是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法高速模数转换器内置DDC结构框图;2 is a structural block diagram of a built-in DDC of a high-speed analog-to-digital converter based on a multi-channel signal preprocessing method of a wideband receiver provided by an embodiment of the present invention;

图3是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法现场可编程门阵列内部DDC实现结构;3 is an internal DDC implementation structure of a field programmable gate array based on a multi-channel signal preprocessing method for a wideband receiver provided by an embodiment of the present invention;

图4是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法窄带信号滤波器示意图;4 is a schematic diagram of a narrowband signal filter based on a multi-channel signal preprocessing method for a wideband receiver provided by an embodiment of the present invention;

图5是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法窄脉冲剔除实现结构;5 is an implementation structure of narrow pulse elimination based on a multi-channel signal preprocessing method for a wideband receiver provided by an embodiment of the present invention;

图6是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法抗异步干扰实现结构图。FIG. 6 is a structural diagram of an implementation of anti-asynchronous interference based on a multi-channel signal preprocessing method based on a wideband receiver provided by an embodiment of the present invention.

具体实施方式Detailed ways

下面结合具体实施例对本发明做进一步详细的描述,但本发明的实施方式不限于此。The present invention will be described in further detail below with reference to specific embodiments, but the embodiments of the present invention are not limited thereto.

请参见图1,图1是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法流程示意图,包括:Please refer to FIG. 1. FIG. 1 is a schematic flowchart of a multi-channel signal preprocessing method based on a wideband receiver provided by an embodiment of the present invention, including:

先确定回波信号的中频f0及调制带宽B,进而可得ADC的采样频率fs,需满足:

Figure BDA0002368341480000041
其中fH=f0+B/2,fL=f0-B/2,fh称为上限频率,fl称为下限频率;由此我们可以确定采样频率fs;由于需要利用模数转换器ADC的内置数字下变频模块,所以还需要确定抽取倍数及滤波器带宽,参数确定后,便可以将模数转换器ADC配置在相应的模式下进行工作;First determine the intermediate frequency f 0 and modulation bandwidth B of the echo signal, and then obtain the sampling frequency f s of the ADC, which must satisfy:
Figure BDA0002368341480000041
Where f H = f 0 +B/2, f L = f 0 -B/2, f h is called the upper limit frequency, and f l is called the lower limit frequency; thus we can determine the sampling frequency f s ; The built-in digital down-conversion module of the converter ADC, so it is necessary to determine the decimation multiple and filter bandwidth. After the parameters are determined, the analog-to-digital converter ADC can be configured to work in the corresponding mode;

多片的模数转换器ADC工作在同一参数下,将多片多通道模数转换器ADC的数据经204b协议编码后的数据帧送至现场可编程门阵列FPGA中进行解帧处理;经解帧后就可以得到第一级数字下变频后的数据;The multi-chip analog-to-digital converter ADC works under the same parameters, and the data frame encoded by the multi-chip multi-channel analog-to-digital converter ADC is sent to the field programmable gate array FPGA for de-frame processing; After the frame, the data after the first-level digital down-conversion can be obtained;

上述处理后数据已被混频到零频附近,且经过两级半带滤波后数据率大大降低,已适合现场可编程门阵列FPGA做信号的预处理,这一步会对其中某一通道的数据进行频谱分析,由现场可编程门阵列FPGA完成快速傅里叶变换FFT运算,得到的FFT计算结果会由光纤发送给信号处理机,由信号处理机进一步找出干扰最小的频点,该频点可用于下次混频的中频值,从而达到频率捷变,可引导雷达工作频率到干扰频谱的空隙或弱区;After the above processing, the data has been mixed to near zero frequency, and the data rate is greatly reduced after two-stage half-band filtering, which is suitable for field programmable gate array FPGA for signal preprocessing. This step will affect the data of one of the channels. Perform spectrum analysis, complete the fast Fourier transform FFT operation by the field programmable gate array FPGA, and the obtained FFT calculation result will be sent to the signal processor by the optical fiber, and the signal processor will further find the frequency point with the least interference. The intermediate frequency value that can be used for the next mixing, so as to achieve frequency agility, and can guide the radar operating frequency to the gap or weak area of the interference spectrum;

雷达工作在变频的模式中,需要完成第二级数字下变频DDC,将信号混频到基带,并进一步完成更大倍数的抽取;混频系数的中心频率由光纤接收的模式字中解析而来;根据频点值,现场可编程门阵列FPGA内利用DDS_IP核产生混频系数,与输入数据完成复乘;将与上一步对应产生的m个混频系数Δθ存储在ROM中,通过解析模式字完成频点选择;使用现场可编程门阵列FPGA中FIR_IP核完成滤波抽取,滤波器的系数可由MATLAB的FDA工具生成,同样将生成的系数存储于rom中,系统启动时该系数自动会被读到FIR_IP中,完成运算;The radar works in the frequency conversion mode, and needs to complete the second-stage digital down-conversion DDC, mix the signal to the baseband, and further complete the extraction of larger multiples; the center frequency of the mixing coefficient is analyzed from the mode word received by the fiber ; According to the frequency point value, the DDS_IP core is used in the field programmable gate array FPGA to generate the mixing coefficients, and the multiplication is completed with the input data; the m mixing coefficients Δθ corresponding to the previous step are stored in the ROM, and the mode word is analyzed by analyzing the mode word. Complete the frequency point selection; use the FIR_IP core in the field programmable gate array FPGA to complete the filter extraction, the coefficients of the filter can be generated by the FDA tool of MATLAB, and the generated coefficients are also stored in the rom, and the coefficients will be automatically read when the system starts. In FIR_IP, complete the operation;

外界强窄脉冲表现为在数字下变频DDC后有少数几个较大的值随着窄脉冲宽度与强度的不同,数字下变频DDC后扩展点的个数不同,会导致在脉压后形成一个台阶,出现虚警;因此需要在完成数字下变频DDC后进行窄脉冲剔除;The external strong and narrow pulses appear as a few larger values after the digital down-conversion DDC. With the difference in the width and intensity of the narrow pulses, the number of expansion points after the digital down-conversion DDC is different, which will lead to the formation of a pulse pressure after the DDC. Steps, false alarms occur; therefore, narrow pulse rejection needs to be performed after the completion of the digital down-conversion DDC;

实现算法:在判断点左右各取PT个保护单元,然后各取若干个点求平均,用相对门限和绝对门限进行判断。若同时满足,则通过置零将该点剔除,否则保留。其中,相对门限的作用是判断该点是否为窄脉冲点,其高位控制该模块的开启与禁用;绝对门限的作用是避免剔除噪声和信号。Realization algorithm: Take PT protection units on the left and right of the judgment point, and then take several points to average, and use relative threshold and absolute threshold to judge. If both are satisfied, remove the point by setting zero, otherwise keep it. Among them, the function of the relative threshold is to judge whether the point is a narrow pulse point, and its high position controls the opening and disabling of the module; the function of the absolute threshold is to avoid eliminating noise and signals.

对各种预处理结果通过光纤发送给信号处理机,同时光纤也会接受由信号处理机发送来的模式字,解析其中关于接收机的各部分工作参数。The various preprocessing results are sent to the signal processor through the optical fiber, and the optical fiber also accepts the mode word sent by the signal processor, and analyzes the working parameters of each part of the receiver.

本发明一种基于宽带接收机多通道信号预处理系统及方法整体实现框图。该系统的处理流程:The present invention is an overall implementation block diagram of a multi-channel signal preprocessing system and method based on a wideband receiver. The processing flow of the system:

由于使用高采样率fs 1GHz左右进行数据采集,高数据率难以直接在现场可编程门阵列FPGA内处理。因此需要在数据进入现场可编程门阵列FPGA前进行预处理以降低数据率。即整个设计方案分两步完成信号的数字下变频。Due to the use of high sampling rate f s around 1GHz for data acquisition, high data rates are difficult to handle directly within a field programmable gate array FPGA. Therefore, it is necessary to preprocess the data before entering the field programmable gate array FPGA to reduce the data rate. That is, the entire design is divided into two steps to complete the digital down-conversion of the signal.

首先在ADC中,利用内置的数字下变频DDC模块,基于全频带完成数据滤波和抽取,其滤波器系数全部固定。First, in the ADC, the built-in digital down-conversion DDC module is used to complete data filtering and decimation based on the full frequency band, and its filter coefficients are all fixed.

如果只使用1级半带滤波器,完成2倍抽取,经仿真验证不会发生频谱混叠。此时数据率降为fs/2为500MHz,虽然可以在现场可编程门阵列FPGA中处理该数据率,但仍然较高,而且会使现场可编程门阵列FPGA后续处理中资源和功耗增加。If only 1-stage half-band filter is used to complete decimation by 2, it is verified by simulation that spectral aliasing will not occur. At this time, the data rate is reduced to f s /2, which is 500MHz. Although this data rate can be processed in the field programmable gate array FPGA, it is still relatively high, and it will increase the resources and power consumption in the subsequent processing of the field programmable gate array FPGA. .

如果使用2级半带滤波器,完成4倍抽取,经仿真验证不会发生频谱混叠。此时数据率降约为250MHz,适合FPGA处理。If a 2-stage half-band filter is used, decimation by a factor of 4 is completed, and it is verified that no spectral aliasing occurs. At this time, the data rate is reduced to about 250MHz, which is suitable for FPGA processing.

在现场可编程门阵列FPGA内进一步完成数字下变频,实现更为精确的混频,且需将数据抽取到更小倍数。在ADC中利用其内置的DDC完成第一次下变频称为DDC_1,理想混频频率为fc,但实际由于NCO位宽太小仅有12bit,导致频率略大于fc,记为fc_actual;信号带宽为B,进行4倍抽取,将数据率降到fs/4,即保证

Figure BDA0002368341480000061
内的信号不混叠。The digital down-conversion is further completed in the field programmable gate array FPGA to achieve more accurate mixing, and the data needs to be extracted to smaller multiples. In the ADC, the built-in DDC is used to complete the first down-conversion, which is called DDC_1. The ideal mixing frequency is f c , but in reality, the NCO bit width is too small and only 12 bits, resulting in a frequency slightly larger than f c , denoted as fc_actual; signal The bandwidth is B, and 4 times decimation is performed to reduce the data rate to f s /4, which is guaranteed
Figure BDA0002368341480000061
The signal inside is not aliased.

将数据输入到现场可编程门阵列FPGA后进行第二次下变频称为DDC_2,全频带不做此处理,根据信号所处频率进行向下移频,或者向上移频再进一步抽取降低数据率。After the data is input to the field programmable gate array FPGA, the second down-conversion is called DDC_2. This processing is not performed in the whole frequency band, and the frequency is shifted down according to the frequency of the signal, or up-shifted and further extracted to reduce the data rate.

请参见图2,图2是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法高速模数转换器内置DDC结构框图,为高速模/数转换器内置数字下变频DDC结构框图;经前端模拟处理后的数据在这一步完成模数转换,并完成第一级数字下变频。采样频率fs由中频采样定理得:

Figure BDA0002368341480000062
其中fH=f0+B/2,fL=f0-B/2确定,模拟数据经过内置ADC模块后完成模数转换,得到数字信号,该信号需混频到基带上,混频所需的本振信号由其内部数控振荡器NCO产生,NCO的频率调谐字可通过下式计算:Please refer to FIG. 2. FIG. 2 is a structural block diagram of a built-in DDC of a high-speed analog-to-digital converter based on a multi-channel signal preprocessing method of a wideband receiver provided by an embodiment of the present invention, which is a structural block diagram of a built-in digital down-conversion DDC of a high-speed analog-to-digital converter. ; The data processed by the front-end analog completes the analog-to-digital conversion in this step, and completes the first-level digital down-conversion. The sampling frequency f s is obtained by the intermediate frequency sampling theorem:
Figure BDA0002368341480000062
Where f H =f 0 +B/2, f L =f 0 -B/2 is determined, the analog data is converted through the built-in ADC module to complete the analog-to-digital conversion, and a digital signal is obtained. The signal needs to be mixed to the baseband. The required local oscillator signal is generated by its internal numerically controlled oscillator NCO, and the frequency tuning word of the NCO can be calculated by the following formula:

Figure BDA0002368341480000063
fc为中频,fs为采样频率,mod为求余函数,round为四舍五入函数,需要说明,由于数控振荡器NCO的位宽仅有12bit,所以得到的频率调谐字并不精确,也说明了在现场可编程门阵列FPGA中完成第二级数字下变频的必要性。ADC内置的滤波器为半带滤波器,最多为四级半带滤波器级联,可完成最多16倍的抽取,原理上可以在ADC内部完成16倍的抽取,减轻优化现场可编程门阵列FPGA运算所需资源,但由于半带滤波器频响特性较差,过渡带较宽,容易发生频谱混叠,所以不在ADC内部完成大倍数的抽取,该部分数字下变频的主要目的为降低数据率,减轻现场可编程门阵列FPGA的运算负担。
Figure BDA0002368341480000063
fc is the intermediate frequency, fs is the sampling frequency, mod is the remainder function, and round is the rounding function. It needs to be explained that since the bit width of the numerical control oscillator NCO is only 12 bits, the frequency tuning word obtained is not accurate, which also shows that in the field The necessity of completing the second stage of digital downconversion in programmable gate array FPGAs. The built-in filter of the ADC is a half-band filter, which can be cascaded up to four stages of half-band filters, which can complete up to 16 times of decimation. However, due to the poor frequency response characteristics of the half-band filter and the wide transition band, spectrum aliasing is prone to occur, so the decimation of large multiples is not completed inside the ADC. The main purpose of this part of the digital down-conversion is to reduce the data rate. , reduce the computational burden of field programmable gate array FPGA.

得到ADC所需的工作参数后,由现场可编程门阵列FPGA通过SPI协议完成配置,配置成功后ADC便可以成功与现场可编程门阵列FPGA完成通信,数据会在ADC内部完成204b协议编码,再经吉比特收发器GTX链路层将数据送到现场可编程门阵列FPGA后中,再由现场可编程门阵列FPGA完成数据解析处理和后续预处理流程。After obtaining the working parameters required by the ADC, the field programmable gate array FPGA completes the configuration through the SPI protocol. After the configuration is successful, the ADC can successfully communicate with the field programmable gate array FPGA, and the data will be encoded in the 204b protocol inside the ADC. The data is sent to the field programmable gate array FPGA through the GTX link layer of the gigabit transceiver, and then the field programmable gate array FPGA completes the data analysis processing and subsequent preprocessing process.

本次设计中,在ADC完成第一级数字下变频DDC后,将数据传给现场可编程门阵列FPGA,在现场可编程门阵列FPGA中对其中某一通道数据进行干扰谱分析。实现算法:选择在每个脉冲重复周期的中后段截取数据以脉冲上升沿为基准,每256个数据为一个单位,截取数据具体位置由模式字确定,先对该段数据进行加窗处理,目的在于增大对副瓣的抑制,然后连续做8组256点快速傅里叶变换FFT,将每组快速傅里叶变换FFT结果累加,并做FFTSHIFT操作,得到256个点。将得到的256个点经光纤发送给信号处理分机,用于信号处理分机进行抗干扰处理,建立查找表,完成256个快速傅里叶变换FFT值与m个频点的对应关系,再由信号处理分机找到干扰最小的频段,添加到模式字中,由查找表的对应关系找到对应的频点,完成频率捷变。In this design, after the ADC completes the first-stage digital down-conversion DDC, the data is transmitted to the field programmable gate array FPGA, and the interference spectrum analysis is performed on the data of one of the channels in the field programmable gate array FPGA. Implementation algorithm: choose to intercept the data in the middle and back sections of each pulse repetition period based on the rising edge of the pulse, and every 256 data is a unit. The purpose is to increase the suppression of side lobes, and then continuously perform 8 groups of 256-point fast Fourier transform FFT, accumulate the results of each group of fast Fourier transform FFT, and perform FFTSHIFT operation to obtain 256 points. The obtained 256 points are sent to the signal processing extension through the optical fiber, which is used for the signal processing extension to perform anti-interference processing, establish a look-up table, and complete the corresponding relationship between the 256 fast Fourier transform FFT values and m frequency points, and then use the signal The processing extension finds the frequency band with the least interference, adds it to the mode word, and finds the corresponding frequency point from the correspondence of the look-up table to complete the frequency agility.

请参见图3,图3是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法现场可编程门阵列内部DDC实现结构,为现场可编程门阵列内部数字下变频DDC实现结构;ADC通过串行高速线将数据传给现场可编程门阵列FPGA,现场可编程门阵列FPGA内经过JESD204B解码后,经过数据重组,将数据恢复为fs/4,32bit。Please refer to FIG. 3 , FIG. 3 is an internal DDC implementation structure of a field programmable gate array based on a multi-channel signal preprocessing method of a wideband receiver provided by an embodiment of the present invention, which is an internal digital down-conversion DDC implementation structure of the field programmable gate array; The ADC transmits the data to the field programmable gate array FPGA through the serial high-speed line. After the JESD204B decoding in the field programmable gate array FPGA, the data is restored to f s /4, 32bit after data recombination.

设计中窄带和宽带切换做,全频带一直处理。图3所示一个通道的处理流程。在FPGA设计中DDS产生的混频系数只与信号的实际中频有关,与信号带宽没有关系,因此混频模块可以宽/窄带复用。但宽窄带滤波结构不同,所以该模块不能复用。对于全频带,取低波束信号在现场可编程门阵列FPGA内暂时只做快速傅里叶变换FFT谱分析处理。现场可编程门阵列FPGA中多路DDC_2的处理方式相同。In the design, the narrowband and wideband switching is done, and the full frequency band is always processed. Figure 3 shows the processing flow of one channel. In the FPGA design, the mixing coefficient generated by DDS is only related to the actual intermediate frequency of the signal, and has nothing to do with the signal bandwidth, so the mixing module can be used for wide/narrowband multiplexing. However, the structure of wide and narrowband filtering is different, so this module cannot be reused. For the full frequency band, the low beam signal is temporarily processed in the field programmable gate array FPGA for fast Fourier transform FFT spectral analysis. The processing method of multiple DDC_2 in the field programmable gate array FPGA is the same.

混频系数的中心频率取决于:由光纤传给现场可编程门阵列FPGA的模式字中的频点值。根据频点值,现场可编程门阵列FPGA内利用DDS_IP核产生混频系数,与输入数据完成复乘。计算公式:

Figure BDA0002368341480000081
rem为求余函数,round四舍五入函数,fs采样频率,f0为中频。经仿真,当N=48时,DDS核产生的混频系数是准确的。计算f0时将ADC产生的频偏考虑在内,将对应产生的m个Δθ存储在ROM中。通过改变Δθ即可产生不同的混频系数。The center frequency of the mixing coefficient depends on the frequency value in the mode word transmitted by the optical fiber to the field programmable gate array FPGA. According to the value of the frequency point, the DDS_IP core is used in the field programmable gate array FPGA to generate the mixing coefficient and complete the complex multiplication with the input data. Calculation formula:
Figure BDA0002368341480000081
rem is the remainder function, round is the rounding function, fs sampling frequency, and f 0 is the intermediate frequency. Through simulation, when N=48, the mixing coefficient generated by the DDS core is accurate. When calculating f 0 , the frequency offset generated by the ADC is taken into account, and the corresponding generated m Δθ are stored in the ROM. Different mixing coefficients can be generated by changing Δθ.

在完成多路数据混频后,进入滤波模块完成多路数据的低通滤波,主要目的是滤除混频引入的“镜频”分量。实现低通滤波时,使用FPGA的FIR_IP核。窄带滤波器系数这里设定为1426阶或2015阶,由于抽取倍数太高,采用改进型“多相滤波结构”,滤波与抽取同时完成。After completing the multi-channel data mixing, enter the filtering module to complete the low-pass filtering of the multi-channel data, the main purpose is to filter out the "image frequency" component introduced by the mixing. When implementing low-pass filtering, use the FIR_IP core of the FPGA. The narrowband filter coefficient is set to 1426th order or 2015th order here. Because the extraction multiple is too high, an improved "polyphase filter structure" is adopted, and the filtering and extraction are completed at the same time.

窄带滤波的低通滤波器LPF设计为1426阶系数为1427个,通带1M,过渡带宽0.5M,带外抑制85dB,窄带滤波器系数使用MATLAB中Fdatool工具设计。窄带滤波如果直接使用多相滤波的结构,需要在1427个系数后补零到1440个,假设需要完成90倍的抽取需分为90路每路16个系数,使用量分析如下:每路I或者Q各用1个DSP48,配置16个系数,使用90+90个DSP48,5路就需要900个DSP48,会浪费很多DSP48。The low-pass filter LPF of narrow-band filtering is designed with 1426 order coefficients, 1427 coefficients, a passband of 1M, a transition bandwidth of 0.5M, and an out-of-band rejection of 85dB. The narrowband filter coefficients are designed using the Fdatool tool in MATLAB. If the narrowband filter directly uses the structure of polyphase filtering, it needs to add zeros to 1440 after 1427 coefficients. Assuming that 90 times of extraction needs to be completed, it needs to be divided into 90 channels with 16 coefficients per channel. The usage analysis is as follows: each channel I or Each Q uses 1 DSP48, configures 16 coefficients, uses 90+90 DSP48, and 5 channels requires 900 DSP48, which will waste a lot of DSP48.

窄带滤波采用新的结构。滤波等效于数据与滤波器系数进行相关,如果每一个数据系数滑动一次完成相关运算,则需要1427个DSP48,然后在进行90倍抽取,得到最终数字下变频DDC结果。本次设计中每90个数据系数滑动一次完成相关运算,滤波、抽取同时完成,在1427个系数后补零到1440个,需要1440/90=16个DSP48,工作速度为180MHz,每个通道I+Q需要32个DSP48。请参见图4,图4是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法窄带信号滤波器示意图,窄带信号滤波器设计示意图如图4所示实虚部处理相同:Narrowband filtering uses a new structure. Filtering is equivalent to correlating data with filter coefficients. If each data coefficient is slid once to complete the correlation operation, 1427 DSP48s are required, and then 90-fold decimation is performed to obtain the final digital down-conversion DDC result. In this design, every 90 data coefficients slide once to complete the correlation operation, and the filtering and decimation are completed at the same time. After 1427 coefficients are filled with zeros to 1440, 1440/90=16 DSP48s are required, the working speed is 180MHz, and each channel I +Q requires 32 DSP48s. Please refer to FIG. 4. FIG. 4 is a schematic diagram of a narrowband signal filter based on a multi-channel signal preprocessing method of a wideband receiver provided by an embodiment of the present invention. The design schematic diagram of the narrowband signal filter is shown in FIG. 4. The real and imaginary part processing is the same:

在宽带模式下,混频后的数据进行开窗处理,窗长及开窗位置由模式字解析而来,将数据率从fs/4稀释至fs/8,以便后续宽带滤波处理。在完成多路数据开窗后,进入宽带滤波模块完成多路数据的低通滤波,主要目的是滤除混频引入的“镜频”分量。实现低通滤波时,使用现场可编程门阵列FPGA的FIR_IP核。宽带滤波器系数由FDA工具产生,由于滤波器系数阶数高,采用“多相滤波结构”,滤波与抽取同时完成。In the broadband mode, the mixed data is windowed, and the window length and window position are parsed from the mode word, and the data rate is diluted from f s /4 to f s /8 for subsequent broadband filtering. After completing the multi-channel data windowing, enter the broadband filtering module to complete the low-pass filtering of the multi-channel data, the main purpose is to filter out the "image frequency" component introduced by the mixing. When implementing low-pass filtering, use the FIR_IP core of the field programmable gate array FPGA. The broadband filter coefficients are generated by the FDA tool. Due to the high order of the filter coefficients, a "polyphase filter structure" is adopted, and the filtering and extraction are completed at the same time.

外界强窄脉冲表现为在数字下变频DDC后有少数几个较大的值随着窄脉冲宽度与强度的不同,DDC后扩展点的个数不同,会导致在脉压后形成一个台阶,出现“虚警”。因此需要在完成DDC后进行窄脉冲剔除。The external strong and narrow pulses appear as a few large values after the digital down-conversion DDC. With the difference in the width and intensity of the narrow pulses, the number of expansion points after DDC is different, which will lead to the formation of a step after the pulse pressure. "False alarm". Therefore, narrow pulse rejection needs to be performed after the DDC is completed.

实现算法:在判断点左右各取PT个保护单元,然后各取n个点求平均,用相对门限和绝对门限进行判断。若同时满足,则通过置零将该点剔除,否则保留。其中,相对门限的作用是判断该点是否为窄脉冲点,其高位控制该模块的开启与禁用;绝对门限的作用是避免剔除噪声和信号。Implementation algorithm: Take PT protection units on the left and right of the judgment point, and then take n points to average, and use relative threshold and absolute threshold to judge. If both are satisfied, remove the point by setting zero, otherwise keep it. Among them, the function of the relative threshold is to judge whether the point is a narrow pulse point, and its high position controls the opening and disabling of the module; the function of the absolute threshold is to avoid eliminating noise and signals.

请参见图5,图5是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法窄脉冲剔除实现结构,其中T表示进行判断的点,PT表示保护单元个数。如果保护单元设置为1,可以剔除强度小于-31dB、宽度为1.0us窄脉冲,但对于近区1.5us的强窄脉冲回波,会剔除中间点,形成凹口。为避免剔除整机发射的近区1.5us的窄脉冲回波,本次设计中只剔除“宽脉冲回波区”的窄脉冲干扰,对应距离为15Km之外。Please refer to FIG. 5 . FIG. 5 is an implementation structure of narrow pulse elimination based on a multi-channel signal preprocessing method of a wideband receiver provided by an embodiment of the present invention, where T represents the point of judgment, and PT represents the number of protection units. If the protection unit is set to 1, narrow pulses with an intensity of less than -31dB and a width of 1.0us can be rejected, but for the strong and narrow pulse echoes of 1.5us in the near area, the middle point will be rejected to form a notch. In order to avoid eliminating the 1.5us narrow pulse echo in the near area emitted by the whole machine, only the narrow pulse interference in the "wide pulse echo area" is excluded in this design, and the corresponding distance is beyond 15Km.

请参见图6,图6是本发明实施例提供的一种基于宽带接收机多通道信号预处理方法抗异步干扰实现结构图,为抗异步干扰实现结构;为防止其它雷达发射信号对信号处理造成影响,进行抗异步处理。由于其它雷达与本机雷达重频不同,因此在相邻脉冲重复周期(PRT)上干扰出现的距离门位置不同;而且重频相差越大,距离门相差越多。Please refer to FIG. 6. FIG. 6 is a structure diagram of an implementation of anti-asynchronous interference based on a multi-channel signal preprocessing method of a wideband receiver provided by an embodiment of the present invention, which is an implementation structure of anti-asynchronous interference; Influence, anti-asynchronous processing. Since the repetition frequency of other radars is different from that of the local radar, the position of the range gate is different in the adjacent pulse repetition period (PRT); and the greater the repetition frequency difference, the greater the range gate difference.

实现算法:用RAM存储一个PRT的数据,每次将1/8*PRT+7/8*RAMPRT表示当前PRT的值,RAM表示上个PRT存到RAM中的值更新到RAM中。用相对门限和绝对门限判断当前PRT和RAM中相同方位的值。若两个门限同时满足,则通过置零将该点剔除,否则保留。相对门限判断该点是异步干扰点,其高位控制该模块的开启与禁用;绝对门限的作用是为避免剔除噪声和信号。Implementation algorithm: use RAM to store the data of a PRT, each time 1/8*PRT+7/8*RAMPRT represents the value of the current PRT, and RAM represents the value of the previous PRT stored in the RAM to be updated to the RAM. Use the relative threshold and the absolute threshold to judge the value of the same orientation in the current PRT and RAM. If the two thresholds are satisfied at the same time, the point is eliminated by setting zero, otherwise it is retained. The relative threshold determines that the point is an asynchronous interference point, and its high bit controls the opening and disabling of the module; the function of the absolute threshold is to avoid eliminating noise and signals.

在完成上述流程之后,多路数据会被打包由光纤发送模块将数据发送至信号处理分机,整套数据接收及预处理流程至此就算完成。After the above process is completed, the multi-channel data will be packaged and sent to the signal processing extension by the optical fiber transmission module, and the entire data reception and preprocessing process is completed at this point.

以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。The above content is a further detailed description of the present invention in combination with specific preferred embodiments, and it cannot be considered that the specific implementation of the present invention is limited to these descriptions. For those of ordinary skill in the technical field of the present invention, without departing from the concept of the present invention, some simple deductions or substitutions can be made, which should be regarded as belonging to the protection scope of the present invention.

Claims (1)

1. A multichannel signal preprocessing method based on a broadband receiver is characterized by comprising the following steps:
determining the intermediate frequency f of the echo signal0And modulating the bandwidth B to obtain the sampling frequency f of the ADCsThe following requirements are met:
Figure FDA0002368341470000011
wherein f isH=f0+B/2,fL=f0-B/2,fhReferred to as the upper limit frequency, flReferred to as the lower limit frequency; from which we can determine the sampling frequency fs(ii) a Because a built-in digital down-conversion module of an analog-to-digital converter (ADC) is required to be utilized, an extraction multiple and a filter bandwidth are required to be determined, and after parameters are determined, the ADC can be configured to work in a corresponding mode;
the multi-chip analog-to-digital converter (ADC) works under the same parameter, and a data frame obtained by encoding data of the multi-chip multi-channel analog-to-digital converter (ADC) through a 204b protocol is sent to a Field Programmable Gate Array (FPGA) for de-framing; after the frame is decoded, data after the first-stage digital down-conversion can be obtained;
the processed data are mixed to be near zero frequency, the data rate is greatly reduced after two-stage half-band filtering, the data are suitable for being preprocessed by a Field Programmable Gate Array (FPGA), the data of a certain channel are subjected to frequency spectrum analysis in the step, the Field Programmable Gate Array (FPGA) completes fast Fourier transform ((FFT) operation, the obtained FFT calculation result is sent to a signal processor by an optical fiber, the signal processor further finds out a frequency point with minimum interference, and the frequency point can be used for an intermediate frequency value of next mixing, so that frequency agility is achieved, and the working frequency of the radar can be guided to a gap or a weak area of an interference frequency spectrum;
when the radar works in a frequency conversion mode, a second-stage Digital Down Conversion (DDC) is required to be completed, signals are mixed to a baseband, and further, extraction of a larger multiple is completed; the center frequency of the mixing coefficient is resolved from the mode word received by the optical fiber; according to the frequency point value, generating a mixing coefficient by using a DDS _ IP core in a Field Programmable Gate Array (FPGA), and completing complex multiplication with input data; storing m mixing coefficients delta theta generated corresponding to the previous step in a ROM, and completing frequency point selection through analyzing mode words; the FIR _ IP core in a Field Programmable Gate Array (FPGA) is used for finishing filtering extraction, the coefficient of the filter can be generated by an FDA tool of MATLAB, the generated coefficient is stored in rom, and the coefficient can be automatically read into the FIR _ IP to finish operation when the system is started;
the external strong and narrow pulse is represented as a few larger values after digital down-conversion (DDC) (along with the difference of the width and the strength of the narrow pulse, the number of extension points after the digital down-conversion (DDC) is different), which can cause a step to be formed after pulse pressure and a false alarm to appear; therefore, narrow pulse rejection is required after Digital Down Conversion (DDC) is completed;
the various preprocessing results are sent to the signal processor through the optical fiber, and the optical fiber also receives the mode words sent by the signal processor and analyzes the working parameters of all parts of the receiver.
CN202010042806.2A 2020-01-15 2020-01-15 Multichannel signal preprocessing method based on broadband receiver Active CN111257838B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010042806.2A CN111257838B (en) 2020-01-15 2020-01-15 Multichannel signal preprocessing method based on broadband receiver

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010042806.2A CN111257838B (en) 2020-01-15 2020-01-15 Multichannel signal preprocessing method based on broadband receiver

Publications (2)

Publication Number Publication Date
CN111257838A true CN111257838A (en) 2020-06-09
CN111257838B CN111257838B (en) 2023-02-10

Family

ID=70948916

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010042806.2A Active CN111257838B (en) 2020-01-15 2020-01-15 Multichannel signal preprocessing method based on broadband receiver

Country Status (1)

Country Link
CN (1) CN111257838B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115436884A (en) * 2022-09-07 2022-12-06 四川九洲电器集团有限责任公司 Multichannel data transmission defect detection and repair method and system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2007352488A1 (en) * 2007-04-25 2008-11-06 Telekom Malaysia Berhad Transceiver front end for software radio systems
US20170010345A1 (en) * 2015-07-10 2017-01-12 Mark Resources, Inc. Direct Sampling of Received Signals in Radar
CN108008359A (en) * 2017-11-08 2018-05-08 武汉滨湖电子有限责任公司 A kind of cascade digital based on pattern-band radio frequency sampling filters anti-Communication Jamming method
CN108627808A (en) * 2017-03-15 2018-10-09 武汉玉航科技有限公司 Radar jammer ultra-wideband digital signal processing method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2007352488A1 (en) * 2007-04-25 2008-11-06 Telekom Malaysia Berhad Transceiver front end for software radio systems
US20170010345A1 (en) * 2015-07-10 2017-01-12 Mark Resources, Inc. Direct Sampling of Received Signals in Radar
CN108627808A (en) * 2017-03-15 2018-10-09 武汉玉航科技有限公司 Radar jammer ultra-wideband digital signal processing method
CN108008359A (en) * 2017-11-08 2018-05-08 武汉滨湖电子有限责任公司 A kind of cascade digital based on pattern-band radio frequency sampling filters anti-Communication Jamming method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王利华: "宽带分数抽取数字下变频设计", 《雷达科学与技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115436884A (en) * 2022-09-07 2022-12-06 四川九洲电器集团有限责任公司 Multichannel data transmission defect detection and repair method and system

Also Published As

Publication number Publication date
CN111257838B (en) 2023-02-10

Similar Documents

Publication Publication Date Title
CN102355281B (en) Combinable frequency hopping wireless transmitting-receiving system and operation method thereof
CN101378263B (en) Multi-carrier digital receiver based on digital intermediate frequency and multi-carrier digital receive method
CN107976572A (en) A kind of channelizing high-precision real Time-frequency Spectrum Analysis method and its system
CN111510255A (en) Frequency hopping signal blind detection and parameter estimation method based on broadband frequency spectrum data
CN108008359A (en) A kind of cascade digital based on pattern-band radio frequency sampling filters anti-Communication Jamming method
CN116466305A (en) A Method for Integrated Detection of Radar Signal Interference
CN111257838B (en) Multichannel signal preprocessing method based on broadband receiver
CN111521981A (en) Multichannel intermediate frequency signal generation method for radar signal source
CN110365351B (en) Method and device for realizing channelized reception based on filter multiplexing architecture
CN109981136A (en) A kind of catching method of Frequency Hopping Communication Signal, device and computer storage medium
CN106872948A (en) Radar working channel acquisition system and its method based on two-stage digital down coversion
CN117955523A (en) Anti-interference method based on cluster unmanned aerial vehicle platform in multi-mode ad hoc network scene
CN202190265U (en) Hybrid Frequency Hopping Wireless Transceiver System
CN107294628A (en) A kind of frequency spectrum intelligent perception system and its method based on USRP
CN106712805B (en) A shortwave pre-selector
CN112379337A (en) Short-distance false alarm eliminating method for long-short pulse combined pulse compression radar
CN105959019B (en) Big bandwidth communication signal receiving analyzer based on binary channels Phase Stacking and method
CN105915193A (en) Improved generation method of multiphase filter
CN111585606B (en) Coherent FFH/DS mixed spread spectrum system
CN110677216B (en) Digital radio frequency front end facing electronic countermeasure and radio frequency signal frequency detection method
CN201114162Y (en) A Multicarrier Digital Receiver System Based on Digital IF Technology
CN113114166A (en) High-speed parallel DDC (direct digital control) and FIR (finite impulse response) filtering processing method based on FPGA (field programmable Gate array)
CN115250153B (en) Digital channelization method and system for adaptive tracking filtering
CN102685049A (en) Fractional order channelized separation method for reaching two linear frequency modulation (LFM) signals at the same time
CN206341215U (en) Device is selected in a kind of shortwave prognosis

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