CN116381613A - Signal environment statistical method for broadband array multi-beam receiver - Google Patents
Signal environment statistical method for broadband array multi-beam receiver Download PDFInfo
- Publication number
- CN116381613A CN116381613A CN202310030623.2A CN202310030623A CN116381613A CN 116381613 A CN116381613 A CN 116381613A CN 202310030623 A CN202310030623 A CN 202310030623A CN 116381613 A CN116381613 A CN 116381613A
- Authority
- CN
- China
- Prior art keywords
- frequency
- amplitude
- long pulse
- pulse frequency
- beat
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000007619 statistical method Methods 0.000 title abstract description 5
- 239000011159 matrix material Substances 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 20
- 238000001228 spectrum Methods 0.000 claims description 9
- 238000007499 fusion processing Methods 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 4
- 230000004927 fusion Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 230000011218 segmentation Effects 0.000 claims description 2
- 230000005855 radiation Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
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)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a signal environment statistical method for a broadband array multi-beam receiver, which comprises the following steps: generating 'frequency and wave beam' two-dimensional matrix amplitude data based on the array element data in the broadband array multi-wave beam receiver; based on the amplitude data of the 'frequency and wave beam' two-dimensional matrix, obtaining a peak searching result formed by a single beat: estimating azimuth resolution units corresponding to the peak searching results; obtaining a long pulse frequency description word based on a peak searching result formed by adjacent beats and a corresponding azimuth resolution unit thereof; reporting the long pulse frequency description word based on a preset minimum beat number threshold and a preset maximum beat number threshold of the long pulse frequency description word; and graphically displaying the reported long pulse frequency description word so that an operator can know the signal environment in the broadband array multi-beam receiver. The invention can output the space-time frequency and energy domain statistical result with moderate data quantity, realizable engineering and higher resolution, and is used for operators to know the signal environment.
Description
Technical Field
The invention relates to the technical field of passive radar reception, in particular to a signal environment statistical method for a broadband array multi-beam receiver.
Background
The passive radar does not emit electromagnetic signals, but detects and discovers the target by intercepting radar signals of a target radiation source. With the development and popularization of various electronic devices, the electromagnetic environment becomes more and more complex, radar signals are often submerged in various communication, interference and clutter signals, and great difficulty is brought to operator judgment and closed loop treatment. Therefore, the current signal environment (energy distribution of various radiation source signals in the dimensions of time, space, frequency and the like) is effectively known, an operator can be helped to focus on a current work task, the work strategy can be timely adjusted, and the signal interception quality of a receiver is improved.
The broadband array multi-beam receiver is more and more commonly applied to passive radars due to the advantages of large bandwidth, wide airspace, high gain and the like, but the real-time processing data volume is particularly huge due to the high sampling rate, the large number of beams and the like. Taking 1GHz sampling, 1024-point FFT and 32 beams as examples, after frequency domain beam forming, output 16 x 10 per second 9 And (5) point frequency domain amplitude-phase data. These data most closely characterize the energy distribution of the current signal environment in the time-frequency dimension (frequency spacing of 1GHz/1024≡1MHz, time spacing of 1024×1/1 ghz=1.024 us, 32 beams). However, such a large amount of data is not acceptable nor practical if the data transmission and graphic display processes are performed directly.
The reduction processing is carried out on the amplitude-phase data after the frequency domain wave beam forming in engineering, and the data quantity is reduced, so that the method is a feasible method. The maximum maintained spectrum is a typical method, and the processing mode is as follows: firstly, forming amplitude data of a frequency domain beam formed by each beat, and taking a maximum value along a beam dimension; then, by accumulating for a long time (a typical value of 200ms including many beats), the maximum hold processing is performed on the amplitude value of each frequency point at a different beat, and the data amount can be effectively compressed. In the literature "autocorrelation method for extracting an antenna scanning period based on a maximum hold spectrum" (electronic information countermeasure technique, 2018, 01, fan Shengzhao, etc.), an antenna scanning period of a stronger signal can be extracted based on maximum hold spectrum data, but the maximum hold processing also causes a large amount of weak signal loss. In addition, the maximum maintenance spectrum also has low time resolution (accumulation duration, typically 200 ms) for the radiation source signal description, and the azimuth/beam dimension distribution information of the signal is lost, so that an operator can know the signal environment in a very rough way.
Disclosure of Invention
In view of the above, the present invention provides a signal environment statistics method for a broadband array multi-beam receiver to solve the above-mentioned technical problems.
The invention discloses a signal environment statistical method for a broadband array multi-beam receiver, which comprises the following steps:
step 1: generating 'frequency and wave beam' two-dimensional matrix amplitude data based on the array element data in the broadband array multi-wave beam receiver;
step 2: based on the frequency and beam two-dimensional matrix amplitude data, obtaining a peak searching result formed by a single beat:
step 3: estimating azimuth resolution units corresponding to the peak searching results;
step 4: obtaining a long pulse frequency description word based on a peak searching result formed by adjacent beats and a corresponding azimuth resolution unit thereof;
step 5: reporting the long pulse frequency description word based on a preset minimum beat number threshold and a preset maximum beat number threshold of the long pulse frequency description word;
step 6: and graphically displaying the reported long pulse frequency description word so that an operator can know the signal environment in the broadband array multi-beam receiver.
Further, the step 1 includes:
and acquiring the array element data in the broadband array multi-beam receiver at a high sampling rate, and generating 'frequency and beam' two-dimensional matrix amplitude data through broadband signal narrowband segmentation and beam forming processing.
The "frequency, beam" two-dimensional matrix amplitude data is expressed as:
k is the number of effective frequency points of the FFT output frequency spectrum; l is the number of multi-beams, y i,j Representing the signal amplitude received by the ith frequency bin, the jth beam.
Further, the step 2 includes:
step 21: performing amplitude peak searching processing based on the amplitude data of the 'frequency and wave beam' two-dimensional matrix;
step 22: the signal amplitude data with known frequency and azimuth are excluded from peak searching processing through beam shadow hiding and frequency shadow hiding processing;
step 23: and sorting the amplitude peak searching results according to the amplitude, and outputting a preset number of peak searching results with larger amplitude.
Further, the step 3 includes:
dividing the azimuth coverage range of the broadband array multi-beam receiver into a plurality of azimuth resolution units, and estimating the corresponding azimuth resolution units according to the peak searching result output in the step 2.
Further, in dividing the azimuth coverage of the broadband array multi-beam receiver into a plurality of azimuth resolution units:
the broadband array multi-beam receiver adopts a contrast method for angle measurement, and can be used as an azimuth resolution unit or further divided into smaller azimuth resolution units in the azimuth range pointed by adjacent beams.
Further, the step 4 includes:
the peak amplitude generated by the adjacent beats is subjected to association fusion processing based on a frequency and azimuth resolution unit, and a multi-beat fusion result is called a long pulse frequency descriptor, wherein the elements of the peak amplitude comprise { beginning beat number, ending beat number, frequency, maximum frequency point, minimum frequency point, azimuth resolution unit number, amplitude mean value, amplitude standard deviation, amplitude maximum value and amplitude minimum value }.
Further, the step 5 includes:
setting a minimum beat number threshold and a maximum beat number threshold reported by the long pulse frequency descriptive word;
in the fusion process, if the long pulse frequency descriptive word exceeds the maximum beat number threshold, immediately reporting the processing;
when the long pulse frequency description word is not associated with the peak value in the next beat, comparing the existing beat number of the long pulse frequency description word with a minimum beat number threshold, and if the existing beat number exceeds the minimum beat number threshold, reporting and outputting.
Further, the minimum beat number threshold is a key parameter for controlling the data amount of the long pulse frequency descriptor, and the larger the minimum beat number threshold of the long pulse frequency descriptor is, the smaller the probability of false-reporting the long pulse frequency descriptor due to random noise is, and the smaller the whole data amount of the long pulse frequency descriptor is in probability.
Further, the step 6 includes:
LFDW graphical display: the display console software receives the reported long pulse frequency description word, and displays the time-frequency and amplitude information contained in the long pulse frequency description word in a graphical mode for an operator to know the signal environment.
Due to the adoption of the technical scheme, the invention has the following advantages:
1. based on the minimum beat number threshold of the reported LFDW, the data quantity of the generated LFDW can be flexibly controlled: when only white noise exists, the probability of searching for a peak value for multiple times is greatly reduced due to the randomness of the amplitude on different frequencies and different beams, and the higher the minimum beat number threshold is, the smaller the corresponding probability is; when a signal with a single frame length/time length larger exists, the multi-section beat spectrum data contained in the signal is described by one piece of LFDW data; compared with the original frequency domain beam forming amplitude-phase data, the total data volume is greatly compressed, and corresponding data transmission and graphic display processing become completely feasible in engineering;
2. compared with the rough description of the maximum maintenance frequency spectrum on the signal environment (the typical value of the time resolution is 200ms, and no azimuth/beam distribution information exists), the LFDW extracted by the invention comprises more complete time-space frequency and energy distribution description information of the signal environment: the time resolution is consistent with the FFT beat length; by dividing the azimuth resolution unit, the azimuth distribution of the signal is described; describing the frequency point distribution information of the LFDW, wherein the frequency resolution is consistent with the FFT frequency resolution; the amplitude information comprises a mean value, a standard deviation, a maximum value, a minimum value and the like, and the categories are more abundant.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments described in the embodiments of the present invention, and other drawings may be obtained according to these drawings for those skilled in the art.
Fig. 1 is a flow chart of a signal environment statistics method for a broadband array multi-beam receiver according to an embodiment of the present invention;
fig. 2 is a schematic diagram of a graphical display sample of LFDW on a display console according to an embodiment of the invention.
Detailed Description
The present invention will be further described with reference to the accompanying drawings and examples, wherein it is apparent that the examples described are only some, but not all, of the examples of the present invention. All other embodiments obtained by those skilled in the art are intended to fall within the scope of the embodiments of the present invention.
The invention is described in further detail below with reference to fig. 1, taking a wideband DBF receiver employing frequency domain beamforming and FFT frequency domain detection as an example:
step 1: the receiver transforms the acquired multi-array element data through FFT, channel amplitude-phase correction, frequency domain beam forming and the like to form 'frequency and beam' two-dimensional amplitude matrix data;
k is the number of effective frequency points of the FFT output frequency spectrum; l is the number of multi-beams, y i,j Representing the signal amplitude received by the ith frequency point and the jth wave beam;
step 2: carrying out peak searching processing on two-dimensional amplitude matrix data frequency by frequency points, wherein each frequency point only selects beam data with the largest amplitude as peak amplitude data, records the amplitude data of left and right adjacent beams, sorts the amplitude peaks of all frequency points according to the amplitude, and selects at most 2-4 peak amplitudes with the largest amplitude to enter subsequent processing, and the rest is discarded;
meanwhile, the operator can set beam shadow and frequency shadow, and strong signals with known frequency and azimuth are not included in peak searching processing, so that other unknown signals can be found through focusing;
step 3, the azimuth range between the adjacent beam directions is used as an azimuth resolution unit, the amplitude values of the peak beam and the adjacent beam are compared, and the azimuth resolution unit where the peak amplitude signal is located can be determined, so that four-element statistical information of each peak amplitude is determined, wherein the four-element statistical information comprises { current section mark number, azimuth resolution unit number, frequency and amplitude };
step 5, judging whether the current LFDW finishes the fusion processing and accords with the output condition according to the LFDW reporting minimum beat number threshold and the maximum beat number threshold set by the user: comprises { beginning section mark, ending section mark, maximum frequency point, minimum frequency point, azimuth resolution unit number, amplitude mean value, amplitude standard deviation, amplitude maximum value, amplitude minimum value }; when the next beat is not associated with a peak value, comparing the existing beat number of the LFDW with a minimum beat number threshold, reporting the output only when the existing beat number of the LFDW exceeds the minimum beat number threshold, discarding the LFDW without exceeding the minimum beat number threshold, and controlling the data quantity output by the LFDW by manually changing the minimum beat number threshold so as to achieve the aim of reducing the data quantity;
Finally, it should be noted that: the above embodiments are only for illustrating the technical aspects of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the above embodiments, it should be understood by those of ordinary skill in the art that: modifications and equivalents may be made to the specific embodiments of the invention without departing from the spirit and scope of the invention, which is intended to be covered by the claims.
Claims (10)
1. A signal environment statistics method for a wideband array multi-beam receiver, comprising the steps of:
step 1: generating 'frequency and wave beam' two-dimensional matrix amplitude data based on the array element data in the broadband array multi-wave beam receiver;
step 2: based on the frequency and beam two-dimensional matrix amplitude data, obtaining a peak searching result formed by a single beat:
step 3: estimating azimuth resolution units corresponding to the peak searching results;
step 4: obtaining a long pulse frequency description word based on a peak searching result formed by adjacent beats and a corresponding azimuth resolution unit thereof;
step 5: reporting the long pulse frequency description word based on a preset minimum beat number threshold and a preset maximum beat number threshold of the long pulse frequency description word;
step 6: and graphically displaying the reported long pulse frequency description word so that an operator can know the signal environment in the broadband array multi-beam receiver.
2. The method according to claim 1, wherein the step 1 comprises:
and acquiring the array element data in the broadband array multi-beam receiver at a high sampling rate, and generating 'frequency and beam' two-dimensional matrix amplitude data through broadband signal narrowband segmentation and beam forming processing.
3. The method of claim 2, wherein the "frequency, beam" two-dimensional matrix amplitude data is represented as:
k is the number of effective frequency points of the FFT output frequency spectrum; l is the number of multi-beams, y i,j Representing the signal amplitude received by the ith frequency bin, the jth beam.
4. The method according to claim 1, wherein the step 2 comprises:
step 21: performing amplitude peak searching processing based on the amplitude data of the 'frequency and wave beam' two-dimensional matrix;
step 22: the signal amplitude data with known frequency and azimuth are excluded from peak searching processing through beam shadow hiding and frequency shadow hiding processing;
step 23: and sorting the amplitude peak searching results according to the amplitude, and outputting a preset number of peak searching results with larger amplitude.
5. The method according to claim 1, wherein the step 3 comprises:
dividing the azimuth coverage range of the broadband array multi-beam receiver into a plurality of azimuth resolution units, and estimating the corresponding azimuth resolution units according to the peak searching result output in the step 2.
6. The method of claim 5, wherein in dividing the azimuth coverage of the wideband array multi-beam receiver into azimuth resolution cells:
the broadband array multi-beam receiver adopts a contrast method for angle measurement, and can be used as an azimuth resolution unit or further divided into smaller azimuth resolution units in the azimuth range pointed by adjacent beams.
7. The method according to claim 1, wherein the step 4 comprises:
the peak amplitude generated by the adjacent beats is subjected to association fusion processing based on a frequency and azimuth resolution unit, and a multi-beat fusion result is called a long pulse frequency descriptor, wherein the elements of the peak amplitude comprise { beginning beat number, ending beat number, frequency, maximum frequency point, minimum frequency point, azimuth resolution unit number, amplitude mean value, amplitude standard deviation, amplitude maximum value and amplitude minimum value }.
8. The method according to claim 1, wherein the step 5 comprises:
setting a minimum beat number threshold and a maximum beat number threshold reported by the long pulse frequency descriptive word;
in the fusion process, if the long pulse frequency descriptive word exceeds the maximum beat number threshold, immediately reporting the processing;
when the long pulse frequency description word is not associated with the peak value in the next beat, comparing the existing beat number of the long pulse frequency description word with a minimum beat number threshold, and if the existing beat number exceeds the minimum beat number threshold, reporting and outputting.
9. The method of claim 8, wherein the minimum beat count threshold is a key parameter for controlling the amount of data of the long pulse frequency descriptors, and wherein the larger the minimum beat count threshold of the long pulse frequency descriptors, the smaller the probability of false reporting the long pulse frequency descriptors due to random noise, and the smaller the probability of the overall amount of data of the long pulse frequency descriptors.
10. The method according to claim 1, wherein the step 6 comprises:
LFDW graphical display: the display console software receives the reported long pulse frequency description word, and displays the time-frequency and amplitude information contained in the long pulse frequency description word in a graphical mode for an operator to know the signal environment.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310030623.2A CN116381613A (en) | 2023-01-09 | 2023-01-09 | Signal environment statistical method for broadband array multi-beam receiver |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310030623.2A CN116381613A (en) | 2023-01-09 | 2023-01-09 | Signal environment statistical method for broadband array multi-beam receiver |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116381613A true CN116381613A (en) | 2023-07-04 |
Family
ID=86973865
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310030623.2A Pending CN116381613A (en) | 2023-01-09 | 2023-01-09 | Signal environment statistical method for broadband array multi-beam receiver |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116381613A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116599606A (en) * | 2023-07-19 | 2023-08-15 | 中国电子科技集团公司第二十九研究所 | Spread spectrum signal receiving method and system based on channelized weighted cross-correlation processing |
-
2023
- 2023-01-09 CN CN202310030623.2A patent/CN116381613A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116599606A (en) * | 2023-07-19 | 2023-08-15 | 中国电子科技集团公司第二十九研究所 | Spread spectrum signal receiving method and system based on channelized weighted cross-correlation processing |
CN116599606B (en) * | 2023-07-19 | 2023-09-19 | 中国电子科技集团公司第二十九研究所 | Spread spectrum signal receiving method and system based on channelized weighted cross-correlation processing |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109444820B (en) | Method for detecting target after interference suppression of multi-channel radar when clutter and interference coexist | |
US6271787B1 (en) | System and method for detecting signals across radar and communications bands | |
CN109655795B (en) | Meter-wave radar co-frequency narrow pulse interference suppression method and system based on waveform entropy | |
CN110988831B (en) | Parameter adjustable detector for signal mismatch in clutter and interference coexistence environment | |
CN110031729A (en) | Detection method, system and the data fusion analytical unit in local discharge signal source | |
CN116381613A (en) | Signal environment statistical method for broadband array multi-beam receiver | |
CN107561509B (en) | Airborne millimeter wave radar power line detection method | |
KR101765239B1 (en) | Method and Apparatus for eliminating false detected Pulse Description Word due to short pulse caused by inter-channel interference in broadband multi-channel receiver | |
CN112904287B (en) | Tracking radiation source parameter estimation method based on passive radar finder and application thereof | |
CN112180331A (en) | Adaptive radio frequency shielding pulse frequency point strategy scheduling method | |
JP2934974B2 (en) | Digital statistical processing apparatus and method for determining signal parameters | |
US20010040525A1 (en) | System and methid for detecting signals across radar and communications bands | |
CN113835077B (en) | Search radar target detection method and system based on variable pulse repetition frequency | |
EP1756983B1 (en) | Method for detecting targets | |
CN113406578B (en) | Distributed unmanned aerial vehicle radar target detection method, device and storage medium | |
US6859161B1 (en) | System for time thresholding | |
CN111123216B (en) | Frequency guiding method based on radar signal histogram | |
CN113589274A (en) | Kick signal identification method and device and terminal | |
CN112098952A (en) | Radar reconnaissance clutter suppression method based on time domain statistical processing | |
CN116299204A (en) | High duty cycle signal discrimination and suppression method for wide bandwidth airspace receiver | |
CN114844514B (en) | Broadband and wideband signal monitoring and evaluating system | |
CN118191756B (en) | Radar signal detection method and system | |
CN115951311A (en) | High duty cycle signal discrimination and suppression method for broadband narrow beam receiver | |
CN118311502A (en) | Radar scanning period estimation method, device, equipment and medium | |
CN117493917A (en) | Parameter overlapping multi-target sorting method based on azimuth quadratic fitting |
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 |