CN107634768B - FPGA-based (field programmable Gate array-based) MWC (Multi-media wall) compression sampling broadband digital receiver PDW (Polymer dispersed W) forming method - Google Patents

FPGA-based (field programmable Gate array-based) MWC (Multi-media wall) compression sampling broadband digital receiver PDW (Polymer dispersed W) forming method Download PDF

Info

Publication number
CN107634768B
CN107634768B CN201710810433.7A CN201710810433A CN107634768B CN 107634768 B CN107634768 B CN 107634768B CN 201710810433 A CN201710810433 A CN 201710810433A CN 107634768 B CN107634768 B CN 107634768B
Authority
CN
China
Prior art keywords
module
signal
fpga
frequency
low
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.)
Active
Application number
CN201710810433.7A
Other languages
Chinese (zh)
Other versions
CN107634768A (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering 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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201710810433.7A priority Critical patent/CN107634768B/en
Publication of CN107634768A publication Critical patent/CN107634768A/en
Application granted granted Critical
Publication of CN107634768B publication Critical patent/CN107634768B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention provides a FPGA-based (field programmable gate array-based) MWC (multi-media wire) compression sampling broadband digital receiver PDW (pulse width modulation) forming method, belonging to the field of electronic countermeasure in information and communication engineering. Firstly, signals are moved to a baseband through a frequency mixing module, baseband compression sampling signals are obtained through a low-pass filter, and then speed reduction extraction is carried out, so that the data volume is greatly reduced, and hardware implementation is facilitated. And inputting the extracted signals into a CORDIC module to extract amplitude and phase, measuring the frequency of a sub-band by using a phase difference frequency measurement method, and calculating the sub-band of the signals by using a multi-channel signal parallel channelization mode through a frequency calculation module to obtain the absolute carrier frequency of the signals. And (4) extracting the pulse by a threshold decision method, and calculating the arrival time and the disappearance time of the signal extracted by the pulse to obtain the pulse width. The invention verifies the feasibility of forming physical realization of the PDW data of the MWC compression sampling structure-based digital receiver, and lays a theoretical and hardware realization foundation for the subsequent FPGA realization of the whole system based on the novel receiver.

Description

FPGA-based (field programmable Gate array-based) MWC (Multi-media wall) compression sampling broadband digital receiver PDW (Polymer dispersed W) forming method
Technical Field
The invention belongs to the field of electronic countermeasure in information and communication engineering, and particularly relates to a method for realizing a novel digital receiver based on MWC compression sampling by utilizing an FPGA (field programmable gate array), carrying out under-Nyquist sampling on a broadband intermediate-frequency signal and generating a corresponding PDW (product data word), wherein the generated PDW is used for sorting and identifying subsequent electronic reconnaissance.
Background
For the currently commonly used low interception probability (LPI) signals (including chirp modulation (LFM) signals, phase-coded (PSK) signals, frequency-coded (FSK) signals, etc.) and frequency-agile signals, the receiving bandwidth of the wideband digital receiver usually contains a plurality of sub-signals with different carrier frequency positions and different bandwidths, even the carrier frequency positions and bandwidths of the sub-signals are time-varying. In the traditional uniform channelization structure, the sub-band division is blind, and the sub-band bandwidth cannot be changed after the sub-band division is determined, so that the problem of channel crossing is caused when a signal is received, the digital signal processing is complex, the sub-band bandwidth is increased, the sensitivity of a receiver is greatly reduced, and the signal-to-noise ratio improvement effect is reduced.
From the literature on the current research of new channelization structures, most of the research is on how to solve the problems of uniform channelization structures by using non-uniform channelization structures, but the channel division is always fixed and has no adaptive capability. The theory of compressive sampling has been studied intensively in this year. The MWC structure proposed by Mishali m. et al in 2010 is the most successful compressive sampling structure at present, and was successfully expanded to the discrete digital domain by the cheng subject group of harabin engineering university in 2016, and used for the design of wideband digital receivers. The receiver structure utilizes the characteristic of frequency domain sparse signals to realize Sub-Nyquist sampling, saves storage space, reduces the number of branch paths, reduces the complexity of a system and successfully solves the problem of channel crossing. But the research is still in the stages of theoretical research and software simulation.
The MWC structure-based compressed sampling broadband digital receiver is successfully built on an FPGA of Xilinx company, carrier frequency, pulse width and arrival time obtained by an input signal are simulated and verified by utilizing Vivado software self-contained simulation capability, and the correctness and feasibility of realizing the MWC structure-based broadband digital receiver by using the FPGA and forming PDW data are verified through verification that the signal information of the input signal is basically consistent with that of the input signal.
Disclosure of Invention
The invention provides a method for forming PDW (pulse width modulation) of an MWC (multi-media wire) compression sampling broadband digital receiver based on an FPGA (field programmable gate array), which is applied to an electronic reconnaissance system and aims to provide a method for generating PDW data on the FPGA by the MWC compression sampling structure based broadband digital receiver and verify feasibility.
The purpose of the invention is realized as follows: (1) the signal is input into a frequency mixing module and multiplied by a Bernoulli pseudo-random sequence and then output to a low-pass filtering module; (2) extracting a baseband signal from the low-pass filtering module and outputting the baseband signal to a speed reduction extraction module; (3) performing corresponding multiple extraction on the signals in the speed reduction extraction module to reduce the data volume and outputting the signals to the CORDIC module; (4) calculating the amplitude and the phase of the signal in a CORDIC module; (5) outputting the amplitude to a pulse extraction module for pulse extraction and outputting to an arrival time measurement module for measuring the arrival time of the pulse; (6) then outputting the amplitude to a pulse width measurement module, and calculating the pulse width in the pulse width measurement module; (7) outputting the phase to a phase difference frequency measurement module, and measuring the relative frequency of a sub-band; (8) and inputting the DTFT of each sub-channel signal into an absolute frequency measurement module to calculate the sub-channel number of the signal, so as to calculate the true carrier frequency of the signal.
The invention is applicable to the following conditions:
(1) inputting signal pulse compression radar signals, wherein the signal pulse compression radar signals comprise conventional radar signals, two-phase encoding signals, four-phase encoding signals, linear frequency modulation signals, frequency agility signals, non-linear frequency modulation signals and frequency encoding signals;
(2) the sampling sequence of the analog signal x (n) satisfies the nyquist sampling theorem and can be AD-collected as a digital signal.
Compared with the prior art, the method provided by the invention has the beneficial effects that: the broadband digital receiver based on the MWC compression sampling structure is realized on the FPGA platform, corresponding PDW data is successfully obtained, feasibility of physical realization of the novel digital receiver structure is verified, and a certain foundation is laid for subsequent FPGA realization of the whole electronic reconnaissance system based on the novel receiver.
Drawings
FIG. 1 is a block diagram of the overall system flow;
FIG. 2 is a simulation diagram of the output of the FPGA mixing module;
FIG. 3 is a diagram of a low-pass FIR filter polyphase structure;
FIG. 4 is a simulation diagram of the output of the FPGA low pass filter module;
FIG. 5 is a simulation diagram of the output of the FPGA slowdown extraction module;
FIG. 6 is a simulation diagram of the output of the FPGA Cordic module;
FIG. 7 is a simulation diagram of FPGA phase difference frequency measurement output;
FIG. 8 is a simulation diagram of FPGA pulse extraction and pulse width measurement output.
Detailed Description
The method provided by the invention is further described with reference to the accompanying drawings:
the invention is realized based on FPGA, and the basic flow comprises the following nine modules: the device comprises a frequency mixing module, a low-pass filtering module, a speed reduction extraction module, a CORDIC module, a phase difference frequency measurement module, a pulse extraction module, an arrival time measurement module, a pulse width measurement module and an absolute frequency measurement module, which are shown in figure 1. The first 8 modules are described with a one-way channelization structure and the 9 th module is described with a multi-way channelization structure.
1) Frequency mixing module
The frequency mixing module has the functions of: multiplying the input signal by a random Bernoulli sequence distributed at random of +/-1, performing convolution operation on a frequency domain in reality, performing frequency spectrum shifting, and shifting the sparse input signal on the frequency domain to a baseband.
The signal to be tested is a complex exponential signal with the sampling frequency of 200MHz, the signal-to-noise ratio of 15dB, the carrier frequency of 41MHz, the sampling point of 4000 points and the pulse width of 18 mu s, and the mixing sequence is 4 paths of 4000 points +1, -1 Bernoulli pseudorandom signals PN1, PN2, PN3 and PN4 which are randomly distributed. The FPGA implementation process is as follows.
1. Firstly, inputting a complex signal to be detected into an FPGA;
2. then inputting the real part, the imaginary part and the pseudorandom series into the FPGA respectively and simultaneously;
3. multiplying the signal and the bernoulli pseudorandom number on the rising edge of each clock;
4. and obtaining a signal which is mixed with 24-bit width and then is transferred to a baseband, and outputting the signal to a low-pass filtering module.
The FPGA simulation result graph is shown in FIG. 2.
2) Low pass filter module
The low-pass filtering module has the functions of: the signal mixed to baseband is extracted.
The input of the module is mixed data with 24 bit width. The low-pass filter adopts an FIR filter, and the parameters of the filter need to be designed in consideration of sampling frequency, passband ripple, stopband attenuation, passband length and stopband length. The number of sub-bands is set to 20 to ensure the system implementation, so that various parameters of the FIR filter can be designed. In the invention, the base band bandwidth is designed to be 10MHz, so the design is carried outThe initial frequency of the transition band is 5MHz, and the cut-off frequency is 10 MHz. The order of the prototype low-pass filter is 43 order, rp=3dB,rsAnd (4) carrying out quantization, truncation and rounding on the obtained filter coefficient, converting the obtained filter coefficient into a 16-system number, and writing the 16-system number into an FPGA program to construct the low-pass filter required by the design.
And writing the polyphase form of the low-pass filter in a low-pass filtering module by adopting a Verilog language and realizing a low-pass FIR filter. The polyphase form of the low-pass filter is shown below.
Let the low-pass FIR filter be
Figure GDA0001460512570000031
Where h (n) is the unit impulse response of the low-pass filter, the polyphase representation is shown in formula (1).
Wherein the content of the first and second substances,m is the number of the extracted multiphase structure, and the multiphase structure diagram is shown in FIG. 3.
It can be seen that the polyphase structure of the low-pass FIR filter is essentially the multiplication of the signal delays by the filter coefficients, respectively, and then summing them together. The FPGA implementation of the low pass filter module is as follows.
1. Writing the low-pass filter coefficients into a low-pass filter module;
2. performing 43-pipeline delay processing by sequentially performing pipeline assignment on 43 24-bit registers, wherein each register is delay of one clock cycle;
3. inputting the mixed signals into a first register to form input signal delay in a multiphase structure;
4. multiplying each register by the corresponding quantized 13-bit filter coefficient and adding the multiplied results in an assign statement;
5. and outputting the result of the assign statement to the speed reduction extraction module.
The simulation results obtained after the mixed signal passes through the low-pass filter are shown in fig. 4.
3) Speed reduction extraction module
The speed reduction extraction module has the following functions: and performing speed reduction processing of extracting one data per 20 clock periods on the signal data output by the low-pass filtering module to obtain baseband compressed sampling data, wherein the sampling rate is 10MHz, the number of the output compressed sampling signal points of each path is reduced to 200 points, and the data volume is reduced. The method provides convenience for subsequent FPGA processing and reduces the occupation of hardware resources. The FPGA implementation method is as follows.
1. Inputting the signal output by the low-pass filtering module;
2. writing a counter in the always block, and counting once every clock period;
3. when the count reaches 20, storing the input signal into a new register, and setting the counter to be 1;
4. intercepting the 11 th bit to the 25 th bit after the speed reduction extraction to form an output with the bit width of 15 bits
5. And outputting the intercepted signal to a CORDIC module.
The FPGA simulation results are shown in fig. 5.
4) CORDIC module
The CORDIC module has the functions of: and (4) approximating and calculating to obtain the amplitude and phase output of the signal by using the real part and the imaginary part of the input signal.
The Y-axis value obtained by the 9-level 2-division pipelining operation is already sufficiently close to 0, and therefore 9 coordinate system rotation operations are determined to be performed. The FPGA implementation method is as follows.
1. Newly building 9 register groups containing X-axis, Y-axis and phase;
2. respectively inputting the real part signal data and the imaginary part signal data into an X-axis register and a Y-axis register to judge whether the highest bit of the Y-axis register is a positive number;
3. rotating counterclockwise if the highest bit is a positive number, and rotating clockwise if the highest bit is a negative number;
4. rotating for 9 times by utilizing the idea of the 2-division method;
5. calculating and updating data in a register corresponding to the pipelining group by using the rotated X-axis and Y-axis coordinates by using the rotated angle;
6. judging the output of the Y axis, and if the low-order output of the Y axis is equal to 0, proving that the CORDIC operation is completed;
7. the amplitude of the signal in the X-axis register and the phase of the signal in the phase register are output.
The output simulation waveform of the FPGA is shown in fig. 6.
5) Phase difference frequency measuring module
The phase difference frequency measurement module has the functions of: and 4-point phase difference averaging is carried out by utilizing the 10-bit wide phase output by the CORDIC module, and the signal frequency in the sub-band is measured by utilizing a phase difference frequency measurement principle. In the digital domain the frequency and phase are first order differential relationships. The required frequency can be obtained by calculating the phase obtained by the CORDIC algorithm according to equation (2). The FPGA implementation process is as follows.
Figure GDA0001460512570000051
In the formula fsFor the sampling frequency, Δ φ (n) is the phase difference, and f is the subband relative frequency.
1. Inputting the phase output in the CORDIC module;
2. performing a pipeline delay process on the phase by using 2 registers, and solving the phase difference;
3. judging that 2 pi is subtracted from the phase difference when the phase difference is larger than + pi, and 2 pi is added to the phase difference when the phase difference is smaller than-pi;
4. 4-point accumulation average is carried out on the phase difference;
5. and multiplying the obtained average value by the sampling frequency and dividing by 2 pi to obtain and output the frequency in the signal sub-band.
The FPGA simulation diagram is shown in FIG. 7. As can be seen from the figure, the measured sub-band frequency is 0.9 MHz.
6) Pulse extraction module and pulse width measurement module
The pulse extraction module has the functions of: and comparing the amplitude output by the CORDIC module with a corresponding threshold value, regarding the signal exceeding the threshold value as an effective pulse, extracting the pulse, and calculating according to the arrival time and the disappearance time to obtain the pulse width. The system realizes the calculation and selection of the threshold value through the formula (3).
Figure GDA0001460512570000052
In the formula, PfFor false alarm probability, Q (x) is a standard normal distribution function,
Figure GDA0001460512570000053
is the noise variance and N is the data length. Its FPGA implementation is as follows.
1. Inputting the amplitude output by the CORDIC module into a pulse extraction module;
2. calculating to obtain a threshold value, wherein the threshold value in the test of the invention is decimal 37;
3. comparing the amplitude value with a threshold value, setting the amplitude value to be 1 when the amplitude value is larger than the threshold value, and simultaneously using accumulation operation to make a counter;
4. and calculating and outputting the pulse width through the counter value.
The FPGA simulation result graph is shown in FIG. 8. It can be seen from fig. 8 that 182 pulse amplitudes are above the threshold, 5ns per clock, due to the sampling frequency of 200 MHz. After the down sampling the clock becomes 100 ns. Therefore, the pulse width is 182 × 0.1 μ s — 18.2 μ s.
7) Time of arrival measurement module
The function of the time measurement module is achieved: the measurement obtains the time taken from the start of the system to the time the rising edge of the pulse is extracted. The FPGA implementation method is as follows.
1. At the beginning of the system, a counter is constructed, and the counter counts once every clock period;
2. in a pulse extraction module, two registers are newly built and carry out 1-level pipeline processing;
3. judging that the pulse extraction value inverted at the next moment is one, and proving that the rising edge of the pulse is captured;
4. the counter is stopped and the time of arrival is calculated.
8) Absolute frequency measuring module
The absolute frequency measurement module has the functions of: and calculating the actual carrier frequency of the signal by using the ratio of the DTFT of the two adjacent channels to obtain the number of the sub-band where the signal is located and matching with the relative frequency of the sub-band.
We design the subband bandwidth to be larger than the bandwidth of the input signal. It is therefore assumed that the carrier frequency of the input signal is only present in the l' th sub-band of each branch of the MWC wideband digital receiver. To find the absolute frequency, multi-channel signal detection is required and mixing is performed using a cyclic shifted periodic pseudo-random sequence. The relation between the sequence of the mth path and the sequence of the first path is shown in formula (4), and the fourier series form is shown in formula (5).
Figure GDA0001460512570000061
In the formula
Figure GDA0001460512570000062
Representing the sequence of principal values p1[n]With MpPeriodic pseudorandom sequence obtained by periodic prolongation of period
Figure GDA0001460512570000063
Rectangular function
Figure GDA0001460512570000065
When n is more than or equal to 0 and less than or equal to MpAnd 1 is 1.
Figure GDA0001460512570000071
The DTFT transform of the mth mixed signal can be written as equation (6).
The ratio of the DTFT transforms of the m +1 th path and the m-th path can be derived as equation (7).
Figure GDA0001460512570000073
Simplifying equation (7) results in an estimate for subband l' as in equation (8).
Figure GDA0001460512570000074
The final absolute frequency estimate can be derived from equation (9), where Δ f is the relative frequency in the subband.
f′c=l′·fp+Δf (9)
The experimental results are as follows:
from the conclusion obtained in fig. 7, the frequency of the subband obtained by frequency measurement is 0.9MHz, and if the signal obtained by calculation of the absolute frequency module is on the 5 th subband, the absolute frequency of the signal obtained by calculation is 40.9MHz, and the input signal is 41MHz, so that the result of frequency measurement of the subband is correct.
From the conclusion of fig. 8, the measured pulse width is 18.2 mus, while the input signal pulse width is 18 mus. Thus, the pulse width measurement is correct.
Therefore, the conclusion can be drawn that the MWC compression sampling receiver built on the FPGA platform in the invention has a correct structure and can correctly output PDW data such as carrier frequency, pulse width, arrival time and the like.
The parameters and results set in the embodiments are only used for verifying the feasibility of the method provided by the invention, the method provided by the invention is not limited to the parameters and results set in the embodiments, and non-essential modifications of the invention still belong to the protection scope of the invention.

Claims (1)

1. A PDW forming method of an MWC compression sampling broadband digital receiver based on an FPGA is characterized in that: (1) the signal is input into a frequency mixing module and multiplied by a Bernoulli pseudo-random sequence and then output to a low-pass filtering module; (2) extracting a baseband signal from the low-pass filtering module and outputting the baseband signal to a speed reduction extraction module; (3) performing corresponding multiple extraction on the signals in the speed reduction extraction module to reduce the data volume and outputting the signals to the CORDIC module; (4) calculating the amplitude and the phase of the signal in a CORDIC module; (5) outputting the amplitude to a pulse extraction module for pulse extraction and outputting to an arrival time measurement module for measuring the arrival time of the pulse; (6) then outputting the amplitude to a pulse width measurement module, and calculating the pulse width in the pulse width measurement module; (7) outputting the phase to a phase difference frequency measurement module, and measuring the relative frequency of a sub-band; (8) inputting the DTFT of each sub-channel signal into an absolute frequency measurement module to calculate the sub-channel number of the signal, thereby calculating the true carrier frequency of the signal;
the frequency mixing module FPGA is realized by (1) firstly inputting a complex signal to be detected into the FPGA; (2) then inputting the real part, the imaginary part and the pseudorandom series into the FPGA respectively and simultaneously; (3) multiplying the signal and the bernoulli pseudorandom number on the rising edge of each clock; (4) obtaining a signal which is transferred to a baseband after the frequency mixing of the 24-bit wide, and outputting the signal to a low-pass filtering module;
the FPGA of the low-pass filtering module is realized by writing a multiphase form of the FPGA in a Verilog language in the low-pass filtering module and realizing a low-pass FIR filter, wherein the multiphase form of the low-pass filter is as follows, and the low-pass FIR filter is designed as
Figure FDA0002235035570000011
Where h (n) is the unit impulse response of the low-pass filter, the polyphase representation is shown in formula (1)
Figure FDA0002235035570000012
Wherein the content of the first and second substances,
Figure FDA0002235035570000013
m is a multiphase structure extraction number, and the FPGA of the low-pass filter module is realized as shown in the following (1) writing a low-pass filter coefficient into the low-pass filter module; (2) performing 43-pipeline delay processing by sequentially performing pipeline assignment on 43 24-bit registers, wherein each register is delay of one clock cycle; (3) inputting the mixed signal to a first registerTo form an input signal delay in a polyphase structure; (4) multiplying each register by the corresponding quantized 13-bit filter coefficient and adding the multiplied results in an assign statement; (5) outputting the result of the assign statement to a speed reduction extraction module;
the FPGA is implemented in the process that (1) a signal output by the low-pass filtering module is input; (2) writing a counter in the always block, and counting once every clock period; (3) when the count reaches 20, storing the input signal into a new register, and setting the counter to be 1; (4) intercepting the 11 th bit to the 25 th bit after the speed reduction extraction to form an output with a 15-bit width; (5) outputting the intercepted signal to a CORDIC module;
the CORDIC module FPGA is realized by (1) newly building 9 register groups comprising an X axis, a Y axis and a phase; (2) respectively inputting the real part signal data and the imaginary part signal data into an X-axis register and a Y-axis register to judge whether the highest bit of the Y-axis register is a positive number; (3) rotating counterclockwise if the highest bit is a positive number, and rotating clockwise if the highest bit is a negative number; (4) rotating for 9 times by utilizing the idea of the 2-division method; (5) calculating and updating data in a register corresponding to the pipelining group by using the rotated X-axis and Y-axis coordinates by using the rotated angle; (6) judging the output of the Y axis, and if the low-order output of the Y axis is equal to 0, proving that the CORDIC operation is completed; (7) outputting the amplitude of the signal in the X-axis register and the phase of the signal in the phase register;
the phase difference frequency measurement module FPGA is realized by the following steps,
Figure FDA0002235035570000021
in the formula fsIs sampling frequency, delta phi (n) is phase difference, f is subband relative frequency, (1) phase input in CORDIC module; (2) performing a pipeline delay process on the phase by using 2 registers, and solving the phase difference; (3) judging that 2 pi is subtracted from the phase difference when the phase difference is larger than + pi, and 2 pi is added to the phase difference when the phase difference is smaller than-pi; (4) 4-point accumulation average is carried out on the phase difference; (5) multiplying the obtained average value by the sampling frequency and dividing by 2 pi to obtain and outputFrequency within a signal sub-band;
the pulse extraction module and the pulse width measurement module FPGA are realized in the process that the calculation and the selection of the threshold value are realized by the formula (3)
In the formula, PfFor false alarm probability, Q (x) is a standard normal distribution function,
Figure FDA0002235035570000023
the variance of the noise is taken as N is the data length, (1) the amplitude output in the CORDIC module is input into a pulse extraction module; (2) calculating to obtain a threshold value; (3) comparing the amplitude value with a threshold value, setting the amplitude value to be 1 when the amplitude value is larger than the threshold value, and simultaneously using accumulation operation to make a counter; (4) calculating and outputting pulse width through the counter value;
the FPGA measurement module is implemented in the process that (1) when a system starts, a counter is constructed, and the counting is carried out once in each clock period; (2) in a pulse extraction module, two registers are newly built and carry out 1-level pipeline processing; (3) judging that the pulse extraction value inverted at the next moment is one, and proving that the rising edge of the pulse is captured; (4) stopping the counter, and calculating the reaching time;
the absolute frequency measurement module FPGA is realized by assuming that the carrier frequency of an input signal only exists in the l' sub-band of each branch of the MWC broadband digital receiver, performing multi-channel signal detection, and performing frequency mixing by adopting a cyclic shift periodic pseudorandom sequence, so that the relation between the sequence of the m-th path and the sequence of the first path is shown in a formula (4), and the Fourier series form is shown in a formula (5),
Figure FDA0002235035570000031
in the formulaRepresenting the sequence of principal values p1[n]With MpPeriodic pseudorandom sequence obtained by periodic prolongation of period
Figure FDA0002235035570000033
Rectangular function
Figure FDA0002235035570000034
When n is more than or equal to 0 and less than or equal to MpWhen the value of-1 is 1,
the DTFT transform of the mth mixing signal can be written as equation (6),
Figure FDA0002235035570000036
the ratio of DTFT transforms for the m +1 th path and the m-th path can be derived as equation (7),
Figure FDA0002235035570000037
simplifying equation (7) yields an estimate of subband l' as in equation (8),
Figure FDA0002235035570000038
the final absolute frequency estimate can be obtained from equation (9),
fc′=l′·fp+Δf (9)
where Δ f is the relative frequency in the subband.
CN201710810433.7A 2017-09-11 2017-09-11 FPGA-based (field programmable Gate array-based) MWC (Multi-media wall) compression sampling broadband digital receiver PDW (Polymer dispersed W) forming method Active CN107634768B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710810433.7A CN107634768B (en) 2017-09-11 2017-09-11 FPGA-based (field programmable Gate array-based) MWC (Multi-media wall) compression sampling broadband digital receiver PDW (Polymer dispersed W) forming method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710810433.7A CN107634768B (en) 2017-09-11 2017-09-11 FPGA-based (field programmable Gate array-based) MWC (Multi-media wall) compression sampling broadband digital receiver PDW (Polymer dispersed W) forming method

Publications (2)

Publication Number Publication Date
CN107634768A CN107634768A (en) 2018-01-26
CN107634768B true CN107634768B (en) 2020-02-14

Family

ID=61100204

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710810433.7A Active CN107634768B (en) 2017-09-11 2017-09-11 FPGA-based (field programmable Gate array-based) MWC (Multi-media wall) compression sampling broadband digital receiver PDW (Polymer dispersed W) forming method

Country Status (1)

Country Link
CN (1) CN107634768B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108549078B (en) * 2018-03-30 2021-11-02 西安电子科技大学 Cross-channel combination and detection method for radar pulse signals
CN108832945A (en) * 2018-06-19 2018-11-16 哈尔滨工程大学 A kind of implementation method of the efficient MWC compression sampling digital receiver structure based on heterogeneous structure
CN108921110A (en) * 2018-07-06 2018-11-30 电子科技大学 Radar signal classification method of the novel convolutional neural networks in conjunction with Wigner-Ville distribution
CN109031289B (en) * 2018-07-13 2020-06-09 清华大学 Cognitive agile frequency conversion radar waveform design method and device
CN109975771B (en) * 2019-03-14 2022-12-09 湖南红船科技有限公司 Broadband digital channelization method based on signal third-order phase difference
CN111447019B (en) * 2020-03-05 2021-12-31 中国电子科技集团公司第二十九研究所 Device for fusing pulse signals among multiple modules
CN111565087B (en) * 2020-04-09 2023-05-02 哈尔滨工程大学 Integrated reconnaissance interference system
CN111769845B (en) * 2020-06-30 2022-03-22 电子科技大学 Weighted superposition channelization method
CN112014810B (en) * 2020-08-07 2024-04-05 西安电子科技大学 High-precision electronic reconnaissance signal parameter measurement method based on FPGA
CN114663990B (en) * 2020-12-22 2024-01-30 天津科畅慧通信息技术有限公司 Intermediate frequency combining method suitable for improving sensitivity of ETC antenna

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003049310A1 (en) * 2001-12-05 2003-06-12 Sicom Systems Ltd. Apparatus and method for a digital, wideband, intercept and analysis processor for frequency hopping signals
CN101499775A (en) * 2009-02-20 2009-08-05 武汉大学 Method and apparatus for pulse compression processing the linear frequency modulation signal by CORDIC
CN102739272A (en) * 2012-06-26 2012-10-17 哈尔滨工程大学 Channelized receiver sub-channel real-time frequency spectrum synthesis method based on field programmable gate array (FPGA)
CN104678364A (en) * 2015-03-13 2015-06-03 哈尔滨工程大学 S-band passive radar interception receiver and signal processing method thereof based on FPGA (Field Programmable Gate Array)
CN104901708A (en) * 2015-01-30 2015-09-09 哈尔滨工程大学 Compressive sampling broadband digital receiver and signal processing method thereof
CN105515695A (en) * 2015-12-04 2016-04-20 哈尔滨工程大学 Compressed sampling signal detection method based on modulated wideband converter

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003049310A1 (en) * 2001-12-05 2003-06-12 Sicom Systems Ltd. Apparatus and method for a digital, wideband, intercept and analysis processor for frequency hopping signals
CN101499775A (en) * 2009-02-20 2009-08-05 武汉大学 Method and apparatus for pulse compression processing the linear frequency modulation signal by CORDIC
CN102739272A (en) * 2012-06-26 2012-10-17 哈尔滨工程大学 Channelized receiver sub-channel real-time frequency spectrum synthesis method based on field programmable gate array (FPGA)
CN104901708A (en) * 2015-01-30 2015-09-09 哈尔滨工程大学 Compressive sampling broadband digital receiver and signal processing method thereof
CN104678364A (en) * 2015-03-13 2015-06-03 哈尔滨工程大学 S-band passive radar interception receiver and signal processing method thereof based on FPGA (Field Programmable Gate Array)
CN105515695A (en) * 2015-12-04 2016-04-20 哈尔滨工程大学 Compressed sampling signal detection method based on modulated wideband converter

Also Published As

Publication number Publication date
CN107634768A (en) 2018-01-26

Similar Documents

Publication Publication Date Title
CN107634768B (en) FPGA-based (field programmable Gate array-based) MWC (Multi-media wall) compression sampling broadband digital receiver PDW (Polymer dispersed W) forming method
CN103941087B (en) The frequency measurement method of the high-frequency cosine signal under lack sampling speed and device thereof
JP6026531B2 (en) Radar pulse detection using a digital receiver for radar
CN109889231B (en) Pulse train signal undersampling method based on random demodulation and finite new information rate
CN103607361A (en) Time frequency overlap signal parameter estimation method under Alpha stable distribution noise
CN112180320B (en) Unmanned aerial vehicle passive positioning system and method
CN109975771B (en) Broadband digital channelization method based on signal third-order phase difference
CA2664537A1 (en) Apparatus, method and computer program product for synthesizing arbitrary waveforms using convolution processors
JP2014077791A (en) Method of indicating correlation between multiple signals, and test measurement device
Ekre Polarity coincidence correlation detection of a weak noise source
Gaarder Scattering function estimation
CN111224672A (en) Multi-harmonic signal undersampling method based on multi-channel time delay
US9331681B2 (en) System and method for gaussian random noise generation
Mishali et al. Generic sensing hardware and real-time reconstruction for structured analog signals
CN108572352A (en) A kind of method for parameter estimation of the phase-coded signal based on lack sampling
CN108696468B (en) Parameter estimation method of two-phase coding signal based on undersampling
CN104218954A (en) Method and device for compressed sampling of broadband array antenna
Orduyılmaz et al. Real-time pulse compression radar waveform generation and digital matched filtering
CN108169715B (en) Method and system for determining in-phase channel phase imbalance degree and quadrature channel phase imbalance degree
Chen et al. A frequency estimation method based on MWC discrete compressed sampling structure
Chaplyha et al. Using non-uniform sampling in real-time correlation processing of authentication signals
Chen et al. A frequency estimation method based on improved MWC discrete compressed sampling structure
Yang et al. Adaptive median threshold algorithm used in FDIS of DSSS receivers
Naeini et al. High speed under-sampling frequency measurements on FPGA
Yang et al. An optimized circulant measurement matrix construction method used in modulated wideband converter for wideband spectrum sensing

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