Frequency hopping network station signal sorting method and system

Abstract

The invention discloses a frequency hopping network station signal sorting method and system, and the method comprises the steps: inputting a plurality of frequency hopping signals to be sorted, and carrying out the clustering according to the frequency hopping frequency and the frequency hopping period, and obtaining a plurality of groups of frequency hopping signal sets; mapping the TOA of each frequency hopping signal in each group of frequency hopping signal set to 0-1 to obtain a TOA mapping value and a distribution diagram of the TOA changing along with time; performing correction processing on each TOA mapping value according to the trend slope angle to obtain a corrected TOA mapping value; and sorting each frequency hopping signal in each frequency hopping signal set again according to the corrected TOA mapping value and the frequency hopping frequency so as to sort out the frequency hopping signals with similar frequency hopping frequency and frequency hopping period in each frequency hopping signal set. According to the method, the sorting accuracy and the sorting efficiency are improved by comprehensively utilizing a plurality of parameters of the frequency hopping signals, so that the frequency hopping signals of each network station can be clearly separated, and the accurate sorting of the frequency hopping signals under limited resource cost is realized.

Description

Technical Field

[0001] The present invention relates to the technical field of frequency hopping signal processing, and in particular to a method and system for sorting frequency hopping network station signals. Background Art

[0002] Frequency hopping communication means that the carrier frequency of the transmitted signal jumps pseudo-randomly. For non-cooperative parties, they can intercept part of the information, but because the hopping pattern of the carrier frequency is unknown, it is difficult to predict and interfere with the communication. For non-cooperative parties, the reconnaissance of frequency hopping signals includes signal reception, detection, sorting and modeling and prediction of sequences. Among them, signal sorting occupies a very important position, because the received signal is often not a single signal. If the frequency hopping signal is not sorted, very effective information cannot be obtained, and the received information will be useless. Due to the advantages of strong anti-interference, low interception probability, high spectrum utilization and strong anti-fading ability, frequency hopping communication has received great attention in the civil and military fields and has developed rapidly in wireless communication technology.

[0003] In a real electromagnetic environment within a limited observation time, the frequency hopping signals intercepted by the reconnaissance equipment often come from multiple frequency hopping stations of both communication parties. In order to analyze the characteristics of a single frequency hopping station signal, it must be sorted out from multiple frequency hopping station signals. And the sorting of frequency hopping signals is mainly based on the parameters of the frequency hopping signal. In the prior art, the frequency hopping station signals are usually sorted according to a single type of frequency hopping signal parameters such as the power parameter, hopping speed, hopping time TOA and directional angle of the frequency hopping signal. However, this type of method of sorting using a single type of frequency hopping signal parameter will have the problems of high cost, poor practicality, insufficient utilization of the frequency hopping signal parameter information, complex calculation, and poor effectiveness of the sorting results. For example, when using the directional angle for sorting, since it is necessary to obtain the azimuth of the frequency hopping signal, it is necessary to perform direction finding in multiple channels to obtain the directional parameters, which will greatly increase the hardware cost, or when using a single frequency hopping signal hopping time TOA for sorting, frequency hopping signals with different hopping speeds and different hopping frequency sets will appear at the same time in a complex electromagnetic environment, and it is impossible to sort them effectively at this time. Summary of the invention

[0004] The technical problem to be solved by the present invention is as follows: In view of the above-mentioned problems in the prior art, a method and system for sorting frequency hopping network station signals are provided to reasonably utilize various frequency hopping signal parameter information to accurately sort frequency hopping signals in a complex electromagnetic environment.

[0005] In order to solve the above technical problems, the technical solution adopted by the present invention is: A method for sorting frequency hopping network station signals, comprising: Step S1, input multiple frequency hopping signals to be sorted, cluster the input frequency hopping signals according to the frequency hopping frequencies and frequency hopping periods of the frequency hopping signals, and obtain multiple groups of frequency hopping signal sets, each group of frequency hopping signal sets corresponding to a category; Step S2, mapping the TOA of each frequency hopping signal in each set of frequency hopping signals to between 0 and 1, obtaining the TOA mapping value of each frequency hopping signal and a distribution diagram of TOA changes over time, wherein the TOA is the arrival time; Step S3, calculating the trend slope angle of the distribution diagram of each frequency hopping signal, and performing correction processing on the TOA mapping value of each frequency hopping signal according to the trend slope angle to obtain a corrected TOA mapping value, wherein the trend slope angle is the angle between a straight line formed by the TOA of the frequency hopping signal in the distribution diagram changing with time and a horizontal time axis; Step S4, sorting the frequency hopping signals in each frequency hopping signal set again according to the modified TOA mapping value and the frequency hopping frequency, so as to sort out the frequency hopping signals with similar frequency hopping frequency and frequency hopping period in each frequency hopping signal set, and obtain the final sorting result output.

[0006] Furthermore, in step S2, the TOA mapping value is obtained by dividing the TOA of the frequency hopping signal by the frequency hopping period of the frequency hopping signal and taking the remainder, and then dividing by the frequency hopping period. The calculation expression is:

[0007] In the above formula, For the The TOA of a frequency hopping signal, is the average value of the frequency hopping period in each set of frequency hopping signals, To take the remainder, For the The mapping value of TOA of a frequency hopping signal on 0-1.

[0008] Further, in step S3, calculating the trend slope angle of the distribution diagram of each frequency hopping signal includes: Step S31, performing Hough transform on the distribution map to obtain Hough space mapping data; Step S32, accumulating and summing the Hough space mapping data in the angle dimension to obtain the projection cumulative sum value at each angle, and taking the angle corresponding to the maximum value of the projection cumulative sum value as the trend slope angle.

[0009] Furthermore, after step S31 and before step S32, the Hough space mapping data is further subjected to interference removal processing, including: Sort the Hough space mapping data in ascending order; Calculate the interference removal threshold value according to the data length and data length ratio index position in the Hough space mapping data sorted in ascending order; The Hough space mapping data are compared with the interference removal threshold value respectively, the value less than the interference removal threshold value is set to 0, and the value greater than or equal to the interference removal threshold value is kept unchanged, so as to obtain the Hough space mapping data after interference removal processing.

[0010] Furthermore, in step S3, the TOA mapping value of each frequency hopping signal is corrected according to the trend slope angle according to the following formula:

[0011] In the above formula, For the The normalized value of a frequency hopping signal in time, For the The mapping value of TOA of a frequency hopping signal on 0-1, is the trend slope angle.

[0012] Furthermore, the step S3 further includes determining whether the calculated corrected TOA mapping value is between 0 and 1, and if so, mapping the corrected TOA mapping value to between 0 and 1 to obtain a final corrected TOA mapping value.

[0013] Furthermore, the modified TOA mapping value is mapped to between 0 and 1 according to the following formula to obtain the final modified TOA mapping value:

[0014] In the above formula, For the The corrected mapping value of TOA of a frequency hopping signal, is the ceiling function, is the floor function, For the The TOA of the frequency hopping signal is finally corrected and mapped.

[0015] Furthermore, after step S2 and before step S3, the following steps are further included: The TOA of each frequency hopping signal is divided by the TOA value of the maximum frequency hopping signal in each set of frequency hopping signals, so as to perform normalization processing on the distribution data in the time dimension.

[0016] The present invention also provides a frequency hopping network station signal sorting system, comprising a microprocessor and a memory connected to each other, wherein the microprocessor is programmed or configured to execute the above frequency hopping network station signal sorting method.

[0017] The present invention also provides a computer-readable storage medium, in which a computer program is stored. The computer program is used to be programmed or configured by a microprocessor to execute the above-mentioned frequency hopping network station signal sorting method.

[0018] Compared with the prior art, the advantages of the present invention are: The present invention performs preliminary sorting by first using the frequency hopping frequency and the frequency hopping period of the frequency hopping signal, and then maps the TOA of the frequency hopping signal to between 0 and 1 and then uses the trend slope angle of the TOA distribution diagram to correct it, so that the TOA parameter can be effectively used and the error accumulated by the TOA value over time can be effectively eliminated, so that the corrected TOA mapping value and the frequency hopping frequency are used to perform secondary sorting on different frequency hopping network station signals, and the network station frequency hopping signals with similar frequency hopping frequencies and frequency hopping periods can be accurately sorted, so that the frequency hopping frequency, the frequency hopping period, the TOA parameter and multiple parameters are comprehensively used to effectively improve the sorting accuracy and sorting efficiency, so that the frequency hopping signal of each network station can be clearly separated, and accurate sorting of the frequency hopping signal is achieved under limited resource costs. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] Figure 1 Schematic diagram of the principles of various parameters of frequency hopping signals.

[0020] Figure 2 Schematic diagram of the implementation flow of the frequency hopping network station signal sorting method of this embodiment.

[0021] Figure 3 It is a distribution scatter plot of the frequency hopping signals of four network stations in three dimensions of frequency hopping signal power, frequency hopping signal hopping frequency and frequency hopping signal hopping period in a specific application embodiment of the present invention.

[0022] Figure 4 This is a characteristic distribution diagram of frequency hopping signals of four networks in a specific application embodiment of the present invention.

[0023] Figure 5 It is a schematic diagram of the principles of various parameters of the frequency hopping signal in a specific application embodiment of the present invention.

[0024] Figure 6 It is a TOA (0-1 mapping)-time distribution diagram of frequency hopping signals of four network stations in a specific application embodiment of the present invention.

[0025] Figure 7 It is a normalized TOA (0-1 mapping)-time distribution diagram of frequency hopping signals of four network stations in a specific application embodiment of the present invention.

[0026] Figure 8 It is a schematic diagram of the mapping result in Hough space obtained in a specific application embodiment of the present invention.

[0027] Fig. 9 It is a schematic diagram of the mapping result after the interference of the Hough space mapping graph obtained in the specific application embodiment of the present invention is removed.

[0028] Fig.10 It is a schematic diagram of the result obtained by accumulating and summing the mapping results after interference removal at various angles in a specific application embodiment of the present invention.

[0029] Fig.11 It is a TOA (0-1 mapping)-time distribution diagram of the frequency hopping signal after the TOA is corrected in the specific application embodiment of the present invention.

[0030] Fig.12 It is a frequency hopping frequency-TOA (0-1 mapping) distribution diagram of the frequency hopping signals of four network stations obtained in a specific application embodiment of the present invention. DETAILED DESCRIPTION

[0031] In order to better understand the above technical solution, the above technical solution will be described in detail below in conjunction with the accompanying drawings and specific implementation methods.

[0032] After completing the parameter measurement of the frequency hopping signal, the characteristic parameters of the frequency hopping signal can be obtained, which mainly include: hopping speed, frequency hopping period, TOA of the frequency hopping signal, dwell time, frequency hopping frequency set, frequency hopping bandwidth, etc. The schematic diagram of the frequency hopping parameters is shown in Figure 1 As shown, where: Hop rate refers to the number of carrier frequency jumps per second, and the unit is hop / s; The frequency hopping period is the inverse of the hopping speed, that is, frequency hopping period = 1 / hopping speed; The TOA (Time of Arrival) of a frequency hopping signal is the time when the frequency hopping signal appears; The dwell time is the duration of each hopping carrier frequency signal; The frequency hopping frequency set is the set of all hopping frequencies of the frequency hopping signal; The frequency hopping bandwidth is the range that the frequency hopping frequency set can cover, that is, the difference between the maximum frequency and the minimum frequency in the frequency set.

[0033] One solution for frequency hopping network station signal sorting is to use DBSCAN clustering algorithm to implement frequency hopping signal sorting. By clustering the characteristic parameters such as azimuth, power, and period of the intercepted signal, the frequency hopping network station signal is sorted, that is, azimuth, power, and period are selected as frequency hopping signal sorting parameters. Although this type of method uses multiple characteristic parameters, each characteristic parameter is used independently, and the actual sorting accuracy and efficiency are not high. There are also problems of high cost and poor practicality. First, the azimuth of the frequency hopping signal must be obtained, which requires multi-channel direction finding, which greatly increases the hardware cost. Secondly, the power parameter feature is used for sorting. In the actual complex electromagnetic environment, the frequency hopping signals sent by the frequency hopping radio station on different frequency hopping channels are quite different, and power is not reliable as a frequency hopping signal sorting parameter. It is actually difficult to separate frequency hopping signals using the frequency hopping cycle. In actual scenarios, the frequency hopping speeds that can be set by third-party communication radio stations are limited (e.g., only 203 hop / s, 609 hop / s, and 1000 hop / s are available). As a result, there are many radio stations with the same hopping speed communicating in the actual environment. In this case, the frequency hopping period of the frequency hopping signal cannot separate the frequency hopping signals of different network stations.

[0034] Another solution is to use the TOA difference statistical histogram method of the frequency hopping signal to sort the frequency hopping signal. The principle is to first count the TOA difference of the frequency hopping signal, estimate one or several possible PRIs (radar signal sorting algorithm based on pulse repetition interval) according to a certain criterion, and then use the obtained PRI to perform sequence search. However, this type of solution only applies the frequency hopping signal TOA parameter to sort the frequency hopping signal, and the actual sorting accuracy is not high, especially in a complex electromagnetic environment, frequency hopping signals with different hopping speeds and different hopping frequency sets will appear at the same time, and only using the TOA parameter of the hopping time cannot effectively distinguish them. In addition, the above method also needs to set two key parameters: K-level TOA difference and threshold for estimating possible PRI. The K-level TOA difference is the value obtained by subtracting the K adjacent frequency hopping signals, and the threshold of PRI is to judge that the signal of such PRI exists when the number of PRI is greater than the threshold value. The values ​​of the above two parameters will directly affect the sorting accuracy, but it is actually difficult to accurately determine the values ​​of these two parameters. For example, the K value is usually configured based on the number of frequency hopping signals, but the actual number of stations that exist before sorting is often unknown. If the K value is set too small, there will be missed detections, and if the K value is set too large, the amount of calculation will increase.

[0035] The TOA of the frequency hopping signal is a key parameter for frequency hopping signal sorting. In real frequency hopping communication application scenarios, there are multiple frequency hopping network station signals communicating at the same time. The radio stations in each network station communicate using synchronous orthogonal networking (radios in the same network station use the same time reference), and the frequency hopping signal time references of different networks are different. Therefore, this feature can be used to use the TOA of the frequency hopping signal to distinguish different frequency hopping network station signals. However, the TOA of the frequency hopping signal is a value that increases over time, and as the TOA value of the frequency hopping signal increases, the TOA estimate will have errors and will continue to accumulate, causing the TOA estimate to shift according to a fixed slope. Therefore, directly using the TOA of the frequency hopping signal cannot accurately achieve sorting.

[0036] The present invention performs preliminary sorting by first using the frequency hopping frequency and the frequency hopping period of the frequency hopping signal, and then maps the TOA of the frequency hopping signal to between 0 and 1 and then uses the trend slope angle of the TOA distribution diagram to correct it, so that the TOA parameter can be effectively used and the error accumulated by the TOA value over time can be effectively eliminated, so that the corrected TOA mapping value and the frequency hopping frequency are used to perform secondary sorting on different frequency hopping network station signals, and the network station frequency hopping signals with similar frequency hopping frequencies and frequency hopping periods can be accurately sorted, so that the frequency hopping frequency, the frequency hopping period, the TOA parameter and multiple parameters are comprehensively used to effectively improve the sorting accuracy and sorting efficiency, so that the frequency hopping signal of each network station can be clearly separated, and accurate sorting of the frequency hopping signal is achieved under limited resource costs.

[0037] like Figure 2 As shown, the steps of the frequency hopping network station signal sorting method in this embodiment include: Step S1, input multiple frequency hopping signals to be sorted, cluster the input frequency hopping signals according to the frequency hopping frequencies and frequency hopping periods of the frequency hopping signals, and obtain multiple groups of frequency hopping signal sets, each group of frequency hopping signal sets corresponding to a category; Step S2, mapping the TOA of each frequency hopping signal in each frequency hopping signal set to between 0 and 1, obtaining the TOA mapping value of each frequency hopping signal and a distribution diagram of TOA changing with time, where TOA is the arrival time; Step S3, calculating the trend slope angle of the distribution diagram of each frequency hopping signal, and performing correction processing on the TOA mapping value of each frequency hopping signal according to the trend slope angle to obtain a corrected TOA mapping value, wherein the trend slope angle is the angle between a straight line formed by the TOA of the frequency hopping signal in the distribution diagram changing with time and the horizontal time axis; Step S4, sorting the frequency hopping signals in each frequency hopping signal set again according to the modified TOA mapping value and the frequency hopping frequency, so as to sort out the frequency hopping signals with similar frequency hopping frequency and frequency hopping period in each frequency hopping signal set, and obtain the final sorting result output.

[0038] It can be understood that this embodiment uses a clustering method to perform preliminary sorting based on the frequency hopping frequency and the frequency hopping period, and preliminarily classifies the frequency hopping signals that are similar in two dimensions of the frequency hopping frequency and the frequency hopping period into one category to form multiple groups of frequency hopping signal sets, and then maps the TOA of each frequency hopping signal in each group of frequency hopping signal sets to between 0 and 1, and at the same time obtains a distribution diagram of TOA changes over time, and uses the distribution diagram to calculate the trend slope angle, and then uses the trend slope angle to correct the TOA mapping value, so as to adaptively correct the TOA mapping value according to the distribution state of TOA, and then uses the corrected TOA mapping value and the frequency hopping frequency to perform secondary sorting on each frequency hopping network station signal in each group of frequency hopping signal sets, so as to accurately sort out the network station frequency hopping signals with similar frequency hopping frequencies and frequency hopping periods, thereby making full use of multiple parameters of the frequency hopping signal to improve the sorting accuracy and sorting efficiency.

[0039] In a specific application embodiment, the specific implementation method of step S1 is described by taking four network station signals in a real external field electromagnetic environment as an example. Figure 3 The figure shows the distribution scatter plot of the four network station signals in the three dimensions of frequency hopping signal power, frequency hopping signal hopping frequency and frequency hopping signal hopping period in the real external field electromagnetic environment. Figure 3 It can be seen that the frequency hopping signals sent by a network station can be clustered into one category. It can also be proved that in the actual electromagnetic environment, even the power variation of the frequency hopping signals of the same network station will be large, so it is not appropriate to use the power of the frequency hopping signal for signal sorting.

[0040] Therefore, the frequency hopping frequency and frequency hopping period parameter characteristics of the network station frequency hopping signal can be used for clustering, and the frequency hopping signals with similar frequency hopping frequency and frequency hopping period can be classified into one category, so as to perform the initial clustering of the network station signal to obtain multiple groups of frequency hopping signal sets. Specifically, Figure 4 It is a distribution scatter plot of the same set of frequency hopping signal data in two dimensions: the frequency hopping frequency of the frequency hopping signal and the frequency hopping period of the frequency hopping signal. Figure 4 It can be seen that the frequency hopping periods of the four network station signals are all concentrated around 1.8, one network station frequency hopping frequency is around 25MHz, one is around 10MHz, and the other two network station frequency hopping signals have frequency hopping frequencies around 0MHz, which are integrated. It can also be seen from the above figure that some network station frequency hopping signals will use similar hopping speeds and hopping frequencies, which makes it more difficult to sort the frequency hopping signals. In order to further separate the network station frequency hopping signals with the same hopping speed and overlapping hopping frequencies, it is necessary to further use other parameter characteristics of the frequency hopping signal.

[0041] The TOA of the frequency hopping signal is a value that increases over time. In order to utilize the parameter feature of TOA, this embodiment first performs mapping processing on TOA. In step S2 of this embodiment, the TOA mapping value is obtained by dividing the TOA of the frequency hopping signal by the frequency hopping period of the frequency hopping signal and then dividing it by the frequency hopping period, thereby mapping the TOA to between 0 and 1. The calculation expression can be expressed as: (1) In the above formula, For the The TOA of a frequency hopping signal, is the average value of the frequency hopping period in each set of frequency hopping signals, To take the remainder, For the The TOA mapping value of a frequency hopping signal on 0-1.

[0042] If the frequency hopping signal is sent by the same network station signal, the TOA mapping result on 0-1 should be the same, because the stations in the same network station follow the same time base. When the frequency hopping period of the frequency hopping signal is used Parameters, and It is also estimated. Once The estimated value has errors. When the TOA value of the frequency hopping signal is extended, the calculation error will continue to accumulate, resulting in the calculated The value will shift with a fixed slope. If the estimation is accurate, then the result of mapping the TOA of the frequency hopping signal of the same network station on 0-1 will be a horizontal line as time extends.

[0043] like Figure 5 As shown in the figure, it shows two frequency hopping network station signals overlapping in the frequency hopping frequency, with the same frequency hopping period or hopping speed, and different TOA of the frequency hopping signal. Assume that the two frequency hopping network station frequency hopping signals (black for one group, green for another group) overlap in the frequency hopping frequency, and have the same frequency hopping period or hopping speed, but the TOA of the frequency hopping signal is different. Assume that the TOA of the first frequency hopping signal of the black frequency hopping network station frequency hopping signal is 0, and the TOA of the second frequency hopping signal is 1, and so on. Figure 5The black dotted line corresponds to the moment. Assuming that the estimated frequency hopping period is accurate at this time, which is 1, the TOA of the first frequency hopping signal is 0 divided by 1 and the remainder is 0, and the TOA of the second frequency hopping signal is 1 divided by 1 and the remainder is 0, and so on. As long as the TOA of the frequency hopping signal of this network station is divided by the accurate frequency hopping period and the remainder is taken, the value of 0 will be obtained, which is the mapping of the TOA of the frequency hopping signal on 0-1. Similarly, the TOA of the first frequency hopping signal of the green frequency hopping network station is 0.5, and the TOA of the second frequency hopping signal is 1.5, and so on. Assuming that the estimated frequency hopping period is accurate at this time, which is 1, the TOA of the first frequency hopping signal is 0.5 divided by 1 and the remainder is 0.5, and the TOA of the second frequency hopping signal is 1.5 divided by 1 and the remainder is 0.5, and so on. As long as the hopping time of the frequency hopping signal of the network station is divided by the accurate frequency hopping period and the remainder is taken, the value of 0.5 will be obtained.

[0044] When the estimated frequency hopping period has errors, it is equivalent to using a ruler with errors to measure data. As the measurement data lengthens, the accumulated errors will become larger and larger. Assuming that the estimated frequency hopping period is 0.9, the arrival time 0 of the first frequency hopping signal is divided by 0.9 and the remainder is 0, the arrival time 1 of the second frequency hopping signal is divided by 0.9 and the remainder is 0.1, the arrival time 2 of the third frequency hopping signal is divided by 0.9 and the remainder is 0.2, and the arrival time 3 of the fourth frequency hopping signal is divided by 0.9 and the remainder is 0.3. It can be seen that with the subsequent hopping times, the mapping of the estimated TOA on 0-1 will change according to a constant trend slope.

[0045] like Figure 6 As shown in the obtained frequency hopping signal TOA (0-1 mapping)-time distribution diagram of the four network stations, due to the error in the estimation of the frequency hopping period, as the frequency hopping signal TOA is extended, the result of calculating the TOA mapping on 0-1 will shift according to a fixed slope. The fixed slope shift actually corresponds to the angle between the straight line formed by the TOA of the frequency hopping signal changing with time in the distribution diagram and the horizontal time axis, that is, the trend slope angle.

[0046] like Figure 7 The figure shows the TOA (0-1 mapping)-time distribution (normalized) diagram of the frequency hopping signal of four stations, where the horizontal direction is defined as 0 degrees, counterclockwise rotation is positive, and clockwise rotation is negative. It can be seen that the angle range of the frequency hopping signal TOA (0-1 mapping)-time distribution data is arrive The trend slope angle is the angle between the straight line formed by the change of TOA over time in the distribution graph corresponding to the distribution data and the horizontal time axis. , so by using the trend slope angle The TOA can be modified by using the value so that the result of TOA mapping on 0-1 can be extended to a horizontal line over time.

[0047] This embodiment calculates the trend slope angle according to the frequency hopping signal TOA (0-1 mapping)-time distribution diagram, and uses the trend slope angle to correct the TOA mapping, which can effectively eliminate the above-mentioned error accumulated over time, so that the TOA parameters can be effectively used to achieve sorting, while improving the sorting accuracy.

[0048] After the mapping is completed, this embodiment needs to perform normalization processing in the time dimension. Specifically, after step S2 and before step S3, it also includes: dividing the TOA of each frequency hopping signal by the TOA value of the maximum frequency hopping signal in each set of frequency hopping signals to normalize the distribution data in the time dimension.

[0049] For example, after obtaining the frequency hopping signal TOA (0-1 mapping)-time distribution data, normalization processing can be performed in the time dimension according to the following formula: (2) In the above formula, is the maximum TOA value in each set of frequency hopping signals, For the A frequency hopping signal TOA, For the The normalized value of a frequency hopping signal in time.

[0050] In step S3 of this embodiment, calculating the trend slope angle of the distribution diagram of each frequency hopping signal includes: Step S31, performing Hough transform on the distribution data to obtain Hough space mapping data; Step S32, performing cumulative summation on the Hough space mapping data in the angle dimension to obtain the cumulative summation value of the projection at each angle, and taking the angle corresponding to the maximum value of the cumulative summation value of the projection as the trend slope angle.

[0051] In a specific application embodiment, taking the four network station signals in a real external electromagnetic environment as an example, after obtaining the following Figure 7 After the frequency hopping signal TOA (0-1 mapping)-time distribution (normalized) diagram shown, what needs to be calculated is Figure 7 The trend slope angle corresponding to the distribution graph In order to calculate Figure 7 The corresponding trend slope angle ,First, the obtained frequency hopping signal TOA (0-1 mapping) - time distribution (normalized) data is ,transformed by Hough transform to obtain the Hough space mapping ,result, as shown in the figure. Figure 8 The figure shows the Hough space mapping of the above data. It can be clearly seen from the figure that the points where the curves overlap (bright spots) are all at an angle, which is Figure 7The corresponding trend slope angle.

[0052] In order to obtain data at points where curves overlap and avoid interference from other data, this embodiment further performs interference removal processing on the Hough space mapping data, specifically including: Step T1, sort the Hough space mapping data in ascending order, the expression is: (3) In the above formula, is the sort value, For all values ​​in Hough space, It is the sort ascending function; Step T2, calculating the interference removal threshold value according to the data length and data length ratio index position in the Hough space mapping data after ascending sorting, the expression is: (4) In the above formula, is the decision threshold to be obtained, To obtain Data length, is the rounding function, is the value of the ratio index position (i.e., the threshold for interference removal, preferably ); Step T3: compare all values ​​of the Hough space mapping data with the de-interference threshold, set values ​​less than the de-interference threshold to 0, and keep values ​​greater than or equal to the de-interference threshold unchanged.

[0053] Specifically, we get After the value, all the values ​​in the Hough space and Compare, greater than or equal to The value remains the same, less than Then set it to 0.

[0054] After interference removal, we can get Fig. 9 The Hough space mapping data interference removal processing result is shown. The data after Hough space mapping interference removal is obtained, and the cumulative sum is performed in the angle dimension to obtain the cumulative sum of projections at each angle. Fig.10 The figure shows the cumulative sum of the above data from -90 to 90 degrees. It can be seen from the figure that the angle value corresponding to the maximum cumulative sum value is around -57 degrees, which is Figure 7 The corresponding trend slope angle .

[0055] In step S3 of this embodiment, according to the trend slope angle The TOA mapping of each frequency hopping signal can be corrected by using the following calculation expression: (5) In the above formula, For the The normalized value of a frequency hopping signal in time, For the The mapping value of TOA of a frequency hopping signal on 0-1, is the trend slope angle, For the A first modified mapping value of the TOA of a frequency hopping signal.

[0056] In step S3 of this embodiment, it is also included to determine whether the calculated modified TOA mapping value is between 0 and 1. If it is, the modified TOA mapping value is mapped to between 0 and 1 to obtain the final modified TOA mapping value. Summarize to the range of 0-1.

[0057] For example, the modified TOA mapping value may be mapped using the following expression: (6) In the above formula, is the ceiling function, is the floor function, For the The TOA of the frequency hopping signal is finally corrected and mapped.

[0058] In a specific application embodiment, the obtained trend slope angle is used to Figure 7 The TOA (0-1 mapping)-time distribution data of the multi-network frequency hopping signal is corrected by TOA (0-1 mapping) to obtain the corrected frequency hopping signal TOA (0-1 mapping)-time distribution data, such as Fig.11 It can be seen that at this time, the result of mapping the frequency hopping signal TOA of the same network station on 0-1 is already a horizontal line extending over time, that is, the influence caused by the frequency hopping period estimation error can be eliminated.

[0059] After obtaining the data corrected by the TOA (0-1 mapping) of the frequency hopping signal, this embodiment uses the corrected result of the TOA (0-1 mapping) of the frequency hopping signal and the information on the two dimensions of the frequency hopping frequency to perform secondary sorting on the result after the preliminary sorting in step S1, so as to sort out the network station frequency hopping signals with similar frequency hopping frequencies and frequency hopping periods after preliminary clustering (that is, the scenarios that cannot be sorted out by the preliminary clustering of the frequency hopping signals in step S1).

[0060] In a specific application embodiment, the frequency hopping signal TOA (0-1 mapping) correction result and the information on the two dimensions of frequency hopping frequency are further sorted, and the distribution diagram of the frequency hopping signal data of the four network stations in the dimensions of frequency and TOA (0-1 mapping) is as follows: Fig.12 As shown. Fig.12 It can be seen that the frequency hopping signals of the four network stations can be clearly separated and sorted out. So far, the frequency hopping network station signal sorting is completed.

[0061] It can be understood that this embodiment can complete the accurate sorting of multiple frequency hopping network station signals in a complex electromagnetic environment. By making full and reasonable use of various frequency hopping characteristic parameters of the frequency hopping signals, the frequency hopping signals can be sorted at a limited resource cost. The processing steps are convenient, and the frequency hopping network station communication scenarios in complex environments where the frequency hopping frequencies and hopping periods of the network stations are similar and it is difficult to distinguish between them are fully considered, effectively improving the accuracy and reliability of sorting.

[0062] The present invention also provides a frequency hopping network station signal sorting system, comprising a microprocessor and a memory connected to each other, wherein the microprocessor is programmed or configured to execute the above frequency hopping network station signal sorting method.

[0063] The present invention also provides a computer-readable storage medium, in which a computer program is stored. The computer program is used to be programmed or configured by a microprocessor to execute the above-mentioned frequency hopping network station signal sorting method.

[0064] The system and medium of the present invention correspond to the above method and also have the advantages described in the above method.

[0065] The present invention implements all or part of the processes in the above-mentioned embodiment method, and can also be completed by instructing related hardware through a computer program. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by the processor, the steps of the above-mentioned method embodiment can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. Computer-readable media include: any entity or device that can carry computer program code, recording medium, U disk, mobile hard disk, disk, optical disk, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, RandomAccess Memory), electric carrier signal, telecommunication signal and software distribution medium. The memory is used to store computer programs and / or modules. The processor implements various functions by running or executing computer programs and / or modules stored in the memory, and calling data stored in the memory. The memory may include a high-speed random access memory and may also include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) card, a Flash Card, at least one disk storage device, a flash memory device, or other volatile solid-state storage devices.

[0066] The above is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above embodiments. All technical solutions under the concept of the present invention belong to the protection scope of the present invention. It should be pointed out that for ordinary technicians in this technical field, some improvements and modifications without departing from the principle of the present invention should also be regarded as the protection scope of the present invention.

Claims

1. A method for sorting frequency hopping network station signals, characterized in that: include: Step S1, input multiple frequency hopping signals to be sorted, cluster the input frequency hopping signals according to the frequency hopping frequencies and frequency hopping periods of the frequency hopping signals, and obtain multiple groups of frequency hopping signal sets, each group of frequency hopping signal sets corresponding to a category; Step S2, mapping the TOA of each frequency hopping signal in each set of frequency hopping signals to between 0 and 1, obtaining the TOA mapping value of each frequency hopping signal and a distribution diagram of TOA changing with time, wherein the TOA is the arrival time, and the TOA is calculated using the frequency hopping period of the frequency hopping signal; Step S3, calculating the trend slope angle of the distribution diagram of each frequency hopping signal, and performing correction processing on the TOA mapping value of each frequency hopping signal according to the trend slope angle to obtain a corrected TOA mapping value, wherein the trend slope angle is the angle between a straight line formed by the TOA of the frequency hopping signal in the distribution diagram changing with time and a horizontal time axis; Step S4, sorting the frequency hopping signals in each frequency hopping signal set again according to the modified TOA mapping value and the frequency hopping frequency, so as to sort out the frequency hopping signals with similar frequency hopping frequency and frequency hopping period in each frequency hopping signal set, and obtain the final sorting result output.

2. The method for selecting frequency hopping network station signals according to claim 1, characterized in that: In step S2, the TOA mapping value is obtained by dividing the TOA of the frequency hopping signal by the frequency hopping period of the frequency hopping signal and taking the remainder and then dividing by the frequency hopping period. The calculation expression is: In the above formula, For the The TOA of a frequency hopping signal, is the average value of the frequency hopping period in each set of frequency hopping signals, To take the remainder, For the The mapping value of TOA of a frequency hopping signal on 0-1.

3. The method for selecting frequency hopping network station signals according to claim 1, characterized in that: In the step S3, calculating the trend slope angle of the distribution diagram of each frequency hopping signal includes: Step S31, performing Hough transform on the distribution map to obtain Hough space mapping data; Step S32, accumulating and summing the Hough space mapping data in the angle dimension to obtain the projection cumulative sum value at each angle, and taking the angle corresponding to the maximum value of the projection cumulative sum value as the trend slope angle.

4. The method for selecting frequency hopping network station signals according to claim 3, characterized in that: After step S31 and before step S32, the method further includes performing interference removal processing on the Hough space mapping data, including: Sort the Hough space mapping data in ascending order; Calculate the interference removal threshold value according to the data length and data length ratio index position in the Hough space mapping data sorted in ascending order; The Hough space mapping data are compared with the interference removal threshold value respectively, the value less than the interference removal threshold value is set to 0, and the value greater than or equal to the interference removal threshold value is kept unchanged, so as to obtain the Hough space mapping data after interference removal processing.

5. The method for selecting frequency hopping network station signals according to claim 1, characterized in that: In the step S3, the TOA mapping value of each frequency hopping signal is corrected according to the trend slope angle according to the following formula: In the above formula, For the The normalized value of a frequency hopping signal in time, For the The mapping value of TOA of a frequency hopping signal on 0-1, is the trend slope angle.

6. The method for selecting frequency hopping network station signals according to claim 5, characterized in that: The step S3 further includes determining whether the calculated corrected TOA mapping value is between 0 and 1, and if so, mapping the corrected TOA mapping value to between 0 and 1 to obtain a final corrected TOA mapping value.

7. The method for selecting frequency hopping network station signals according to claim 6, characterized in that: The corrected TOA mapping value is mapped to between 0 and 1 according to the following formula to obtain the final corrected TOA mapping value: In the above formula, For the The corrected mapping value of TOA of a frequency hopping signal, is the ceiling function, is the floor function, For the The TOA of the frequency hopping signal is finally corrected and mapped.

8. The method for selecting frequency hopping network station signals according to any one of claims 1 to 6, characterized in that: After step S2 and before step S3, the following steps are also included: The TOA of each frequency hopping signal is divided by the TOA value of the maximum frequency hopping signal in each set of frequency hopping signals, so as to perform normalization processing on the distribution data in the time dimension.

9. A frequency hopping network station signal sorting system, comprising a microprocessor and a memory connected to each other, characterized in that: The microprocessor is programmed or configured to execute the frequency hopping network station signal sorting method as described in any one of claims 1 to 8.

10. A computer-readable storage medium having a computer program stored therein, characterized in that: The computer program is used to be programmed or configured by a microprocessor to execute the frequency hopping network station signal sorting method as described in any one of claims 1 to 8.

Images