A complex environment multi-radiation source pulse angle correlation method
By employing multi-channel acquisition and angle correlation methods, the problem of sorting and tracking radiation sources in complex electromagnetic environments was solved, enabling accurate direction finding and acquisition of characteristic parameters of radiation sources.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- SHANGHAI RADIO EQUIP RES INST
- Filing Date
- 2023-11-23
- Publication Date
- 2026-07-10
AI Technical Summary
In complex electromagnetic environments, traditional time-frequency domain methods cannot accurately sort and track parameter-agile radiation sources or frequency-modulated signals, especially low-intercept signals.
The radiation source signal is acquired through multi-channel acquisition, AD sampling and frequency domain conversion are performed to obtain amplitude and phase data, threshold detection is performed to generate pulse descriptors, the elevation and azimuth angles of the pulse descriptors are used to perform angle correlation to generate correlation centers, and the characteristic parameters of the radiation source are calculated.
It enables accurate direction finding, sorting, and tracking of multiple radiation sources in complex environments, and obtains angular, time, and frequency domain information of the radiation sources.
Smart Images
Figure CN117706476B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of radiation source signal sorting and tracking, and in particular to a method for associating the angles of multiple radiation source targets in complex environments. Background Technology
[0002] The environment contains a large number of communication signals and radar signals. In a certain space, time and spectrum, electromagnetic signals crisscross, intersect, overlap densely and have uneven energy distribution, forming an extremely complex electromagnetic environment.
[0003] Most current broadband anti-radiation receivers have large instantaneous bandwidth and wide field of view, receiving pulse densities reaching 1 million pulses per second. Furthermore, radar signals (especially low-intercept signals) are highly variable. Through processing methods such as high-sensitivity receivers, high-speed acquisition circuits, and software-defined radio, information in multiple dimensions, including the time, frequency, and angular domains, can be obtained. Sorting and identifying this multi-dimensional information allows for the identification of radiation source characteristics and stable tracking.
[0004] However, in complex environments, it is no longer possible to sort and identify radiation source targets using traditional time-frequency domain information, especially for parameter-agile radiation sources or frequency-modulated signals, where time-frequency domain parameters cannot accurately sort and identify signals. Since traditional time-frequency domain features are no longer adequate for today's complex electromagnetic environment, researchers are actively seeking new features for signal sorting. Summary of the Invention
[0005] The purpose of this invention is to provide a method for pulse angle correlation of multiple radiation sources in complex environments. This invention enables the direction finding, sorting, and tracking of multiple radiation sources.
[0006] To achieve the above objectives, the present invention provides a method for correlation of pulse angles from multiple radiation sources in complex environments, comprising the following steps:
[0007] S1. Acquire signals from multiple radiation sources through multiple channels; perform AD sampling on the multi-radiation source signals acquired by each channel; perform frequency domain conversion on the AD sampling signals of each channel to obtain the amplitude and phase data of each sampling signal under each channel; sum the amplitudes of the corresponding AD sampling signals under multiple channels to obtain multiple amplitude sum signals;
[0008] S2. Threshold detection is performed on the multiple amplitudes and signals to obtain multiple pulse descriptors, wherein the pulse descriptors contain the pitch angle and azimuth angle of the corresponding radiation source signal in the inertial frame;
[0009] S3. Based on the pitch angle and azimuth angle in the pulse descriptor, perform angular association to generate multiple association centers for multiple pulse descriptors;
[0010] S4. Based on the multiple association centers, obtain the characteristics of multiple radiation sources.
[0011] Optionally, the pulse descriptor corresponds to multiple consecutive amplitudes and signals; the amplitude of any one of the multiple consecutive amplitudes and signals exceeds a set system threshold.
[0012] Optionally, let P i ={t i ,pw i ,f i ,α i ,β i} represents the i-th pulse descriptor word, i = 0, 1, 2…N, where N is the total number of pulse descriptor words; among the multiple consecutive amplitude and signal signals corresponding to the i-th pulse descriptor word, the acquisition times of the first and last amplitude and signal signals exceeding the system threshold are t and t, respectively. i t i ′;pw i pw is the width of the i-th pulse descriptor word. i =t i ′-t i ;f i Let α be the frequency of the radiation source signal corresponding to the i-th pulse descriptor. i β i Let be the elevation and azimuth angles of the radiation source signal corresponding to the i-th pulse descriptor in the inertial coordinate system.
[0013] Optionally, step S3 includes:
[0014] S31. Establish a first set A1, which initially contains all pulse descriptors; let k be a count value, with an initial value of 0;
[0015] S32, Let C k ={(α) k ′,β k ′),m k B k} represents the association center; B k Indicates with C k The set of associated pulse descriptors, the initial B k It includes arbitrarily selecting a pulse descriptor from the first set A1, and denoting this pulse descriptor as P. k,0 ;(α k ′,β k ′) is the correlation center C in the inertial frame. k The pitch angle, azimuth angle, α k ′、β k The initial value of ′ is the pulse description word P in the inertial coordinate system. k,0 The corresponding elevation and azimuth angles of the radiation source signal, m k For set Bk The number of pulse description words in m k The initial value is 1; let the count value j = 1;
[0016] S33, based on (α) k ′,β k ′), and set the second set A2 = A1 - B k The pulse descriptors in the text are respectively associated with the center C. k Perform angle correlation. If a pulse descriptor in the second set A2 satisfies the correlation condition, denote the pulse descriptor as P. k,j Update the associated center C k B k Updated to B k ∪{P k,j}; will m k Updated to m k +1; will (α) k ′,β k Updated to inertial coordinate system, pulse description word P k,j The elevation and azimuth angles of the corresponding radiation source signals; update j to j+1;
[0017] Repeat step S33 until there are no pulse descriptors that satisfy the association condition in the second set A2, then proceed to S34;
[0018] S34. Update the first set A1 to A1-B. k Determine whether the updated first set A1 is an empty set; if not, update k to k+1 and proceed to step S32; if yes, proceed to step S4.
[0019] Optional, in inertial coordinates, the pulse description word P k,j The elevation and azimuth angles of the corresponding radiation source signal are denoted as α. k,j β k,j If (α) k,j ,β k,j ) and (α k ′,β k The Euclidean distance d of ′) k,j If the distance is less than the set distance threshold, then the pulse description word P k,j With the associated center C k The association conditions are met, where:
[0020]
[0021] Optionally, the system threshold is obtained by superimposing a preset fixed threshold onto the system noise.
[0022] Optional, let the association center C be denoted as C. kThe corresponding radiation source is the k-th radiation source; the characteristics of the k-th radiation source include:
[0023] Minimum pulse width of the k-th radiation source for The minimum value in the width;
[0024] The maximum pulse width of the kth radiation source for The maximum value in the width;
[0025] Central pulse width value of k-radiation source for and The average value;
[0026] The minimum frequency of the k-th radiation source for The minimum value in the corresponding radiation source signal frequency;
[0027] The maximum frequency of the k-th radiation source for The minimum value in the corresponding radiation source signal frequency;
[0028] The frequency center value of the kth radiation source for and The average value.
[0029] Compared with the prior art, the beneficial effects of the present invention are as follows:
[0030] The present invention provides a method for pulse angle correlation of multiple radiation sources in complex environments. This method involves AD sampling of multi-channel acquired multi-radiation source signals, followed by time-frequency domain conversion to obtain corresponding amplitude and phase data. The amplitudes of the AD sampled signals from each channel are summed to obtain multiple amplitude sum signals. Threshold detection and amplitude-phase direction finding are then performed on these multiple amplitude sum signals to obtain pulse descriptors. Finally, using angle as a reference, Euclidean distance is used to correlate the pulse descriptors. The latest elevation and azimuth angles (in inertial coordinates) of the radiation source signal corresponding to the most recently correlated pulse descriptor are used as the elevation and azimuth angles of the correlation center. Radiation source characteristic parameters (including minimum pulse width, maximum pulse width, pulse width center value, minimum frequency, maximum frequency, and frequency center value) are calculated based on the pulse descriptors that meet the correlation requirements, thereby obtaining accurate information about the radiation source in the angle, time, and frequency domains. Attached Figure Description
[0031] To more clearly illustrate the technical solution of the present invention, the accompanying drawings used in the description will be briefly introduced below. Obviously, the drawings in the following description are one embodiment of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort:
[0032] Figure 1 , Figure 2 This is a flowchart of a method for correlating pulse angles from multiple radiation sources in a complex environment, as described in an embodiment of the present invention.
[0033] Figure 3 This is a schematic diagram of multiple pulse descriptors in an embodiment of the present invention; Detailed Implementation
[0034] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.
[0035] It should be understood that, when used in this specification and the appended claims, the term "comprising" indicates the presence of the described features, integrals, steps, operations, elements and / or components, but does not exclude the presence or addition of one or more other features, integrals, steps, operations, elements, components and / or collections thereof.
[0036] It should also be understood that the terminology used in this specification is for the purpose of describing particular embodiments only and is not intended to limit the scope of the application. As used in this specification and the appended claims, the singular forms “a,” “an,” and “the” are intended to include the plural forms unless the context clearly indicates otherwise.
[0037] It should also be further understood that the term “and / or” as used in this application specification and the appended claims means any combination of one or more of the associated listed items and all possible combinations, and includes such combinations.
[0038] As used in this specification and the appended claims, the term "if" may be interpreted, depending on the context, as "when," "once," "in response to determination," or "in response to detection." Similarly, the phrases "if determined" or "if [described condition or event] is detected" may be interpreted, depending on the context, as "once determined," "in response to determination," "once [described condition or event] is detected," or "in response to detection of [described condition or event]."
[0039] Furthermore, in the description of this application, the terms "first," "second," "third," etc., are used only to distinguish descriptions and should not be construed as indicating or implying relative importance.
[0040] This invention provides a method for correlating pulse angles from multiple radiation sources in complex environments, such as... Figure 1 , Figure 2 As shown, the steps include:
[0041] S1. Acquire signals from multiple radiation sources through multiple channels; perform AD sampling on the multi-radiation source signals acquired by each channel; perform frequency domain conversion on the AD sampling signals of each channel to obtain the amplitude and phase data of each sampling signal under each channel; sum the amplitudes of the corresponding AD sampling signals under multiple channels to obtain multiple amplitude sum signals.
[0042] Multiple channels can be understood as multiple antennas. It should be noted that each channel will simultaneously collect signals from multiple radiation sources, rather than each channel collecting signals from a single radiation source.
[0043] like Figure 3 As shown, Figure 3 Each black dot in the diagram corresponds to an amplitude and a signal. The amplitude of each amplitude and signal is equal to the sum of the amplitudes of the corresponding AD sampling signals of each channel at the same time.
[0044] S2. Threshold detection is performed on the multiple amplitudes and signals to obtain multiple pulse descriptors, wherein the pulse descriptors contain the pitch angle and azimuth angle of the corresponding radiation source signal in the inertial frame;
[0045] like Figure 3 As shown, the pulse descriptor corresponds to multiple consecutive amplitudes and signals, and the amplitude of any one of these consecutive amplitudes and signals exceeds a set system threshold. In one embodiment, the system threshold is obtained by superimposing a preset fixed threshold onto system noise.
[0046] Let P i ={t i ,pw i ,f i ,α i ,β i Let} represent the i-th pulse descriptor word, where i = 0, 1, 2…N, and N is the total number of pulse descriptor words. Among the multiple consecutive amplitude and signal values corresponding to the i-th pulse descriptor word, the acquisition times of the first and last amplitude and signal values exceeding the system threshold are t and t, respectively. i (also known as arrival time), t i (Also known as the end time). pw i pw is the width of the i-th pulse descriptor word (referred to as pulse width). i =ti ′-t i . Figure 3 The image shows two pulse descriptors, along with their corresponding pulse widths and arrival times. i Let be the frequency of the radiation source signal corresponding to the i-th pulse descriptor word.
[0047] α i β i Let α be the elevation and azimuth angles of the radiation source signal corresponding to the i-th pulse descriptor in the inertial coordinate system. i β i From θ i , It is derived from the conversion, θ i , In the coordinate system of the system (e.g., a missile or shipboard system, which can be understood as the system where the antenna is located), the i-th pulse description word P i The corresponding elevation and azimuth angles of the radiation source signal. Due to system jitter or swaying, it is necessary to convert the system's attitude angles in the inertial frame to attitude angles in the accelerometer frame.
[0048] In one embodiment, direction finding can be performed based on the amplitude and phase data corresponding to the arrival time to obtain θ. i , This can be achieved using existing direction-finding algorithms (such as spatial spectrum direction-finding algorithms). This is not the focus of this invention and will not be elaborated upon here. It should be noted that the direction-finding results for the amplitude and phase data corresponding to each channel are consistent.
[0049] S3. Based on the pitch angle and azimuth angle in the pulse descriptor, perform angular association to generate multiple association centers for multiple pulse descriptors;
[0050] The principle of this invention is that, in multidimensional parameters, the target cannot undergo abrupt changes in the spatial domain, so the identification and tracking of radiation source targets can be achieved based on angle information.
[0051] Step S3 includes:
[0052] S31. Establish a first set A1, which initially contains all pulse descriptors; let k be a count value, with an initial value of 0;
[0053] S32, Let C k ={(α) k ′,β k ′),m k B k} represents the association center; Bk represents the set of pulse descriptors associated with Ck, and the initial B kIt includes arbitrarily selecting a pulse descriptor from the first set A1, and denoting this pulse descriptor as P. k,0 ;(α k ′,β k ′) is the correlation center C in the inertial frame. k The pitch angle, azimuth angle, α k ′、β k The initial value of ′ is the pulse description word P in the inertial coordinate system. k,0 The corresponding elevation and azimuth angles of the radiation source signal, m k For set B k The number of pulse description words in m k The initial value is 1; let the count value j = 1;
[0054] S33, based on (α) k ′,β k ′), and set the second set A2 = A1 - B k The pulse descriptors in the text are respectively associated with the center C. k Perform angle correlation. If a pulse descriptor in the second set A2 satisfies the correlation condition, denote the pulse descriptor as P. k,j Update the associated center C k B k Updated to B k ∪{P k,j}; will m k Updated to m k +1; will (α) k ′,β k Updated to inertial coordinates, the pulse description word is denoted as P. k,j The elevation and azimuth angles of the corresponding radiation source signals; update j to j+1;
[0055] In one embodiment, the pulse description word P in the inertial coordinate system k,j The elevation and azimuth angles corresponding to the radiation source signal are denoted as α. k,j β k,j If (α) k,j ,β k,j ) and (α k ′,β k The Euclidean distance d of ′) k,j If the distance is less than the set distance threshold, then the pulse description word P k,j With the associated center C k The association conditions are met, where:
[0056]
[0057] Repeat step S33 until there are no pulse descriptors that satisfy the association condition in the second set A2, then proceed to S34;
[0058] S34. Update the first set A1 to A1-B k Determine whether the updated first set A1 is an empty set; if not, update k to k+1 and proceed to step S32; if yes, proceed to step S4.
[0059] S4. Based on the multiple association centers, obtain the characteristics of multiple radiation sources.
[0060] Remember the association center C k The corresponding radiation source is the k-th radiation source; the characteristics of the k-th radiation source include:
[0061] Minimum pulse width of the k-th radiation source for The minimum value in the width;
[0062] The maximum pulse width of the k-th radiation source for The maximum value in the width;
[0063] Central pulse width value of k-radiation source for and The average value;
[0064] The minimum frequency of the k-th radiation source for The minimum value in the corresponding radiation source signal frequency;
[0065] The maximum frequency of the k-th radiation source for The minimum value in the corresponding radiation source signal frequency;
[0066] The frequency center value of the kth radiation source for and The average value.
[0067] The present invention provides a method for pulse angle correlation of multiple radiation sources in complex environments. This method involves AD sampling of multi-channel acquired multi-radiation source signals, followed by time-frequency domain conversion to obtain corresponding amplitude and phase data. The amplitudes of the AD sampled signals from each channel are summed to obtain multiple amplitude sum signals. Threshold detection and amplitude-phase direction finding are then performed on these multiple amplitude sum signals to obtain pulse descriptors. Finally, using angle as a reference, Euclidean distance is used to correlate the pulse descriptors. The latest elevation and azimuth angles (in inertial coordinates) of the radiation source signal corresponding to the most recently correlated pulse descriptor are used as the elevation and azimuth angles of the correlation center. Radiation source characteristic parameters (including minimum pulse width, maximum pulse width, pulse width center value, minimum frequency, maximum frequency, and frequency center value) are calculated based on the pulse descriptors that meet the correlation requirements, thereby obtaining accurate information about the radiation source in the angle, time, and frequency domains.
[0068] It should be understood that the sequence number of each step in the above embodiments does not imply the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this application.
[0069] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any person skilled in the art can easily conceive of various equivalent modifications or substitutions within the technical scope disclosed in the present invention, and these modifications or substitutions should all be covered within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.
Claims
1. A method for correlating pulse angles from multiple radiation sources in complex environments, characterized in that, Including the following steps: S1. Acquire signals from multiple radiation sources through multiple channels; perform AD sampling on the multi-radiation source signals acquired by each channel; perform frequency domain conversion on the AD sampling signals of each channel to obtain the amplitude and phase data of each sampling signal under each channel; sum the amplitudes of the corresponding AD sampling signals under multiple channels to obtain multiple amplitude sum signals; S2. Threshold detection is performed on the multiple amplitudes and signals to obtain multiple pulse descriptors, wherein the pulse descriptors contain the pitch angle and azimuth angle of the corresponding radiation source signal in the inertial frame; S3. Based on the pitch angle and azimuth angle in the pulse descriptor, perform angular association to generate multiple association centers for multiple pulse descriptors; Step S3 includes: S31. Establish the first set The initial first set Includes all pulse descriptors; let For count values, The initial value is 0; S32, Order Indicates the associated center; Indicates and The set of associated pulse descriptors, initial Includes from the first set Choose any pulse descriptor from the list and denote it as . ; Correlation center in inertial frame Angle, The initial value is in the inertial coordinate system, pulse description word The elevation and azimuth angles of the corresponding radiation source signals, For set The number of pulse description words in the text The initial value is 1; let the count value be 1. ; S33, based on The second set The pulse descriptors in the text are respectively associated with the center. Perform angular correlation, if the second set If a pulse descriptor exists that satisfies the association condition, then this pulse descriptor is denoted as... Update the association center :Will Updated to ;Will Updated to ;Will Updated to inertial coordinate system, pulse description word The elevation and azimuth angles of the corresponding radiation source signal; Updated to ; Repeat step S33 until the second set. If no pulse descriptor that meets the association condition is found, proceed to S34; S34, Set the first set Updated to - Determine the first set after the update. Is it an empty set? If not, Updated to If yes, proceed to step S32; otherwise, proceed to step S4. S4. Based on the multiple association centers, obtain the characteristics of multiple radiation sources.
2. The method for correlation of pulse angles from multiple radiation sources in complex environments as described in claim 1, characterized in that, The pulse descriptor corresponds to multiple consecutive amplitudes and signals; the amplitude of any one of the multiple consecutive amplitudes and signals exceeds a set system threshold.
3. The method for correlation of pulse angles from multiple radiation sources in complex environments as described in claim 2, characterized in that, make Indicates the first Each pulse description word, , The total number of pulse description words; the first Among the multiple consecutive amplitudes and signals corresponding to each pulse descriptor, the acquisition times of the first and last amplitudes and signals exceeding the system threshold are respectively... , ; For the first The width of each pulse description word, ; For the first The frequency of the radiation source signal corresponding to each pulse descriptor word In the inertial coordinate system, the first The elevation and azimuth angles of the radiation source signal corresponding to each pulse description word.
4. The method for correlation of pulse angles from multiple radiation sources in complex environments as described in claim 1, characterized in that, Pulse description word in inertial coordinate system The elevation and azimuth angles of the corresponding radiation source signals are denoted as follows: ,like and Euclidean distance If the distance is less than the set distance threshold, then the pulse description word... In the associated center The association conditions are met, where: 。 5. The method for correlation of pulse angles from multiple radiation sources in complex environments as described in claim 2, characterized in that, The system threshold is obtained by superimposing a preset fixed threshold onto the system noise.
6. The method for correlation of pulse angles from multiple radiation sources in complex environments as described in claim 2, characterized in that, Record the associated center The corresponding radiation source is the first Radiation source; The characteristics of a radiation source include: No. Minimum pulse width of the radiation source , for The minimum value in the width; No. Maximum pulse width of radiation source , for The maximum value in the width; Central value of the pulse width of the radiation source , for and The average value; No. Minimum frequency of radiation source , for The minimum value in the corresponding radiation source signal frequency; No. Maximum frequency of radiation source , for The minimum value in the corresponding radiation source signal frequency; No. The center frequency value of the radiation source , for and The average value.