CN111257838A - Multichannel signal preprocessing method based on broadband receiver - Google Patents
Multichannel signal preprocessing method based on broadband receiver Download PDFInfo
- 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
Links
- 238000007781 pre-processing Methods 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000005070 sampling Methods 0.000 claims abstract description 13
- 238000006243 chemical reaction Methods 0.000 claims description 37
- 238000001914 filtration Methods 0.000 claims description 20
- 239000013307 optical fiber Substances 0.000 claims description 13
- 238000000605 extraction Methods 0.000 claims description 11
- 238000010183 spectrum analysis Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000035485 pulse pressure Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 238000009432 framing Methods 0.000 claims 1
- 238000012545 processing Methods 0.000 abstract description 22
- 238000010586 diagram Methods 0.000 description 9
- 238000013461 design Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000002699 waste material Substances 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/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting 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
本发明属于雷达信号处理领域,具体涉及一种基于宽带接收机多通道信号预处理方法,满足对大带宽信号的采样处理亦能完成雷达信号的预处理,进而提高系统的抗干扰性能。
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.
Description
技术领域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,需满足:其中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: 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,需满足:其中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: 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,即保证内的信号不混叠。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 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由中频采样定理得:其中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: 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:
fc为中频,fs为采样频率,mod为求余函数,round为四舍五入函数,需要说明,由于数控振荡器NCO的位宽仅有12bit,所以得到的频率调谐字并不精确,也说明了在现场可编程门阵列FPGA中完成第二级数字下变频的必要性。ADC内置的滤波器为半带滤波器,最多为四级半带滤波器级联,可完成最多16倍的抽取,原理上可以在ADC内部完成16倍的抽取,减轻优化现场可编程门阵列FPGA运算所需资源,但由于半带滤波器频响特性较差,过渡带较宽,容易发生频谱混叠,所以不在ADC内部完成大倍数的抽取,该部分数字下变频的主要目的为降低数据率,减轻现场可编程门阵列FPGA的运算负担。 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核产生混频系数,与输入数据完成复乘。计算公式: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: 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
在完成上述流程之后,多路数据会被打包由光纤发送模块将数据发送至信号处理分机,整套数据接收及预处理流程至此就算完成。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)
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)
| 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)
| 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 |
-
2020
- 2020-01-15 CN CN202010042806.2A patent/CN111257838B/en active Active
Patent Citations (4)
| 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)
| Title |
|---|
| 王利华: "宽带分数抽取数字下变频设计", 《雷达科学与技术》 * |
Cited By (1)
| 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 |