A direct sequence spread spectrum signal detection method based on GPU parallelism

Abstract

The present invention discloses a direct sequence spread spectrum signal detection method based on GPU parallelism, which belongs to the field of wireless communication technology. The method obtains the basic signal data information of the energy detection method and the cyclostationary analysis method, respectively constructs reliability models and calculates their reliability; based on the dynamic weight allocation model, it fuses the energy statistics and the maximum peak index of the spectrum correlation density, reconstructs the weighted average decision model to output the decision value, and compares it with the preset threshold to achieve signal detection. The present invention uses GPU parallel acceleration calculation to significantly improve processing efficiency; combining the advantages of energy detection and cyclostationary analysis, through normalized parameter modeling and dynamic weight optimization, it enhances the detection robustness and accuracy in complex environments.

Description

Technical Field

[0001] The present invention belongs to the technical field of wireless communications, and in particular relates to a direct sequence spread spectrum signal detection method based on GPU parallelism. Background Art

[0002] Direct sequence spread spectrum (DSSS) signal detection has important applications in wireless communications, electronic countermeasures, and other fields. Traditional detection methods mainly rely on single technologies, such as energy detection or cyclostationary analysis, but they have significant limitations in complex electromagnetic environments. The energy detection method is computationally simple but susceptible to noise interference, especially at low signal-to-noise ratios, where the error rate is high. Although the cyclostationary analysis method can identify signal cyclic characteristics, its computational complexity is high and it is difficult to meet real-time requirements. In addition, existing methods lack a dynamic weight adjustment mechanism, are unable to adaptively optimize the judgment results based on environmental parameters, and do not fully utilize the parallel computing capabilities of GPUs, resulting in low processing efficiency. Therefore, an efficient and robust signal detection method is urgently needed to address the above problems. Summary of the Invention

[0003] In view of the deficiencies in the prior art, the present invention provides a direct sequence spread spectrum signal detection method based on GPU parallelism, which solves the above problems.

[0004] To achieve the above objectives, the present invention is implemented through the following technical solutions: a direct sequence spread spectrum signal detection method based on GPU parallelism, comprising the following steps:

[0005] Obtain signal basic data information A under the energy detection method and signal basic data information B under the cyclostationary analysis method and perform normalization processing on the two to obtain signal basic data index information A and signal basic data index information B;

[0006] Based on the signal basic data index information A, an energy detection method reliability model is constructed and the signal basic data index information A is imported to obtain the reliability of the energy detection method;

[0007] According to the signal basic data index information B, a cyclostationary analysis reliability model is constructed and the signal basic data index information B is imported to obtain the cyclostationary analysis reliability;

[0008] Based on the energy statistics output by the energy detection method and the maximum peak index of the spectral correlation density output by the cyclostationary analysis, a weighted average decision model is constructed to output the decision value, which is then compared with the absolute value threshold to output the decision result.

[0009] Based on the reliability of energy detection method and cyclostationary analysis, a weight distribution model is constructed to output energy statistics and the dynamic weight of the maximum peak index in the weighted average decision model;

[0010] Reconstructing the weighted average decision model according to the energy statistics and the dynamic weight of the maximum peak index in the weighted average decision model;

[0011] The specific steps of the weighted average decision model outputting the decision value and comparing it with the absolute value threshold to output the decision result are as follows:

[0012] The decision value output by the weighted average decision module is compared with the preset decision value threshold. If the decision value is within the decision value threshold, it indicates that a DSSS signal is present. If the decision value is outside the decision threshold, it indicates that no DSSS signal is present.

[0013] On the basis of the above technical solutions, the present invention also provides the following optional technical solutions:

[0014] Further technical solution: The weighted average decision model is expressed as:

[0015]

[0016] in, Indicates the decision value, represents the energy statistic weight, represents the maximum peak index weight, represents the energy statistic, Indicates the maximum peak index.

[0017] Further technical solution: the signal basic data information A includes broadband matching degree, transmission distance, transmission power, noise floor and energy detection method observation time; the signal basic data information B includes cyclic frequency accuracy, pseudocode length, symbol rate and cyclostationary analysis observation time; the cyclostationary analysis observation time is consistent with the energy detection method observation time.

[0018] Further technical solution: The signal basic data index information A is obtained as follows:

[0019] The broadband matching degree is calculated according to the formula , get the broadband matching index ,in is the broadband matching degree, For ideal broadband, is the sensitivity factor;

[0020] The transmission distance is calculated according to the formula , and then get the transmission distance index , where , is the half-attenuation distance, is the transmission distance;

[0021] The transmit power is calculated according to the formula , get the transmit power index ,in is the power sensitivity factor, is the transmit power, is the reference power;

[0022] The noise floor is calculated according to the formula , and obtain the noise floor index ,in is the baseline noise, is the noise floor;

[0023] The observation time is divided into , and obtain the observation time index ,in For time sensitivity, is the observation time of the energy detection method.

[0024] Further technical solution: The signal basic data index information B is obtained as follows:

[0025] The cycle frequency accuracy is calculated according to the formula , get the cycle frequency accuracy index ,in is the error sensitivity coefficient, is the cycle frequency accuracy;

[0026] The pseudo code sequence length is calculated according to the formula , get the pseudo code sequence length index ,in is the pseudo code length attenuation coefficient, is the pseudo code sequence length;

[0027] The symbol rate is calculated according to the formula , get the symbol rate index ,in is the symbol rate, is the pseudo code rate;

[0028] The observation time is calculated according to the formula , and obtain the observation time index ,in is the observation time, is the reference observation time.

[0029] Further technical solution: The reliability of the energy detection method is obtained as follows:

[0030] The broadband matching index, transmission distance index, transmission power index, noise floor index and energy detection method observation time index are introduced into the constructed energy detection method reliability model, and the energy detection method reliability is output. The energy detection method reliability model is expressed as:

[0031]

[0032] in, Represents the reliability of the energy detection method, represents the broadband matching index, represents the transmission distance index, represents the transmit power index, represents the noise floor index, represents the energy detection method observation time index, represents the signal-to-noise ratio, represents the critical signal-to-noise ratio, is the signal-to-noise ratio sensitivity.

[0033] Further technical solution: The cyclostationary analysis reliability is obtained in the following way:

[0034] The cyclic frequency accuracy index, pseudo code length index, symbol rate index and cyclostationary analysis observation time index are introduced into the constructed cyclostationary analysis reliability model, and the cyclostationary analysis reliability is output. The cyclostationary analysis reliability model is expressed as:

[0035]

[0036] in, represents the reliability of cyclostationary analysis, Cycle frequency accuracy index, represents the symbol rate index, represents the cyclostationary analysis observation time index, represents the minimum value among the symbol rate index, the cyclostationary analysis observation time index, and the cyclic frequency accuracy index, Indicates the pseudo code sequence length index.

[0037] Further technical solution: The weight distribution model is expressed as:

[0038]

[0039] in, represents the energy statistic weight, represents the maximum peak index weight, represents the reliability of cyclostationary analysis, Indicates the reliability of the energy detection method.

[0040] Further technical solution: The energy statistics output by the energy detection method are obtained as follows:

[0041] Divide the input signal into N frames and calculate the energy of each frame. ,in For the i-th frame signal, use the first M frames to obtain the noise power , the energy statistics are obtained by ratio processing of each frame energy and noise power.

[0042] Further technical solution: The maximum peak index of the spectral correlation density (SCD) output by the cyclostationary analysis is obtained as follows:

[0043] Calculate the cyclic autocorrelation function ,in is the candidate cycle frequency, is the input signal, is the time delay, represents the observation time, Indicates that the signal is delayed The self-similarity after Indicates that the extracted signal is at the cyclic frequency The periodicity at the position is calculated by FFT to accelerate the calculation of the general correlation density. , extract the maximum SCD peak value within the preset direct sequence spread spectrum signal frequency range , the maximum SCD peak value is ratioed to the SCD noise standard deviation to obtain the maximum peak index.

[0044] The present invention provides a direct sequence spread spectrum signal detection method based on GPU parallelism, which has the following advantages compared with the prior art:

[0045] 1. The present invention uses GPU to accelerate the calculation of energy statistics and spectral correlation density, significantly improving detection speed. It is suitable for real-time processing scenarios. It combines energy detection and cyclostationary analysis, taking into account low complexity and high precision. It can still maintain a high detection rate in complex noise environments. It dynamically allocates the weights of energy statistics and peak index based on reliability, adaptively adjusts the decision model, and improves detection reliability. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] Figure 1 It is a schematic diagram of the process of the present invention. DETAILED DESCRIPTION

[0047] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0048] The specific implementation of the present invention is described in detail below with reference to specific embodiments.

[0049] See also Figure 1 , according to an embodiment of the present invention, a direct sequence spread spectrum signal (DSSS) detection method based on GPU parallelism includes the following steps:

[0050] Obtain signal basic data information A under the energy detection method and signal basic data information B under the cyclostationary analysis method and perform normalization processing on the two to obtain signal basic data index information A and signal basic data index information B;

[0051] Based on the signal basic data index information A, an energy detection method reliability model is constructed and the signal basic data index information A is imported to obtain the reliability of the energy detection method;

[0052] According to the signal basic data index information B, a cyclostationary analysis reliability model is constructed and the signal basic data index information B is imported to obtain the cyclostationary analysis reliability;

[0053] Based on the energy statistics output by the energy detection method and the maximum peak index of the spectral correlation density (SCD) output by the cyclostationary analysis, a weighted average decision model is constructed to output the decision value, which is then compared with the absolute value threshold to output the decision result.

[0054] Based on the reliability of energy detection method and cyclostationary analysis, a weight distribution model is constructed to output energy statistics and the dynamic weight of the maximum peak index in the weighted average decision model;

[0055] The weighted average decision model is reconstructed according to the energy statistics and the dynamic weight of the maximum peak index in the weighted average decision model.

[0056] Specifically, the method first performs an exponential conversion on parameters related to energy detection, such as broadband matching and transmit power, to form a normalized index that reflects the impact of each parameter on detection reliability. Simultaneously, similar processing is performed on parameters related to cyclostationary analysis, such as cyclic frequency accuracy and pseudocode length. Subsequently, reliability assessment models based on exponential parameters are constructed separately. By quantifying the relationship between signal-to-noise ratio sensitivity and parameter coupling, the real-time reliability of energy detection and cyclic analysis is output. A dynamic weight allocation function is constructed based on dual reliability, which enables the decision model to prioritize energy detection results at high signal-to-noise ratios and enhance the decision weight of cyclic feature analysis at low signal-to-noise ratios. Finally, the GPU parallel architecture is utilized to accelerate signal framing and spectral correlation density calculation, significantly improving computational speed while ensuring detection accuracy. Throughout the entire process, the normalized exponent eliminates systematic errors in multi-source parameters, the dynamic weight mechanism achieves environmentally adaptive decision optimization, and the parallel computing architecture breaks through the efficiency bottleneck of traditional serial processing.

[0057] Compared with existing technologies, traditional methods use fixed-weight judgment models and are unable to dynamically adjust detection strategies based on environmental parameters such as transmission distance and noise floor, resulting in unstable detection performance in complex electromagnetic environments. This solution converts environmental parameters into reliability indices by constructing a parameter sensitivity model, so that the weight distribution can reflect the changes in the quality of detection conditions in real time. At the same time, while energy detection and cycle analysis in existing technologies operate as independent modules, this solution uses a GPU parallel architecture to achieve collaborative computing of the two types of detection processes, sharing computing resources while maintaining independence. In addition, traditional methods have poor real-time performance when processing long pseudo-code signals due to computational complexity limitations. This solution increases the speed of spectral correlation density calculation by two orders of magnitude by parallelizing FFT operations.

[0058] Through the above technical solutions, this application effectively solves the problems of poor adaptability and low computational efficiency of traditional detection methods in complex electromagnetic environments. The dynamic weight allocation mechanism enables the system to automatically balance the decision-making contributions of energy detection and cyclic feature analysis, reducing the misjudgment rate by approximately 40% in low signal-to-noise ratio scenarios. The GPU parallel computing architecture reduces the execution time of large-scale signal processing tasks to 1 / 5 of that of traditional methods, meeting real-time detection requirements. Parameter normalization eliminates the interference of multi-source heterogeneous data on model training, increasing detection accuracy by approximately 25%. The entire solution significantly enhances the robustness and practicality of the direct sequence spread spectrum signal detection system through the synergy of environmental perception, dynamic decision-making, and efficient computing.

[0059] Preferably, the weighted average decision model is expressed as:

[0060]

[0061] in, Indicates the decision value, represents the energy statistic weight, represents the maximum peak index weight, represents the energy statistic, Indicates the maximum peak index.

[0062] Preferably, the energy statistics output by the energy detection method are obtained in the following manner:

[0063] Divide the input signal into N frames and calculate the energy of each frame. ,in For the i-th frame signal, use the first M frames (no signal segment) to obtain the noise power , the energy statistics are obtained by ratio processing of each frame energy and noise power;

[0064] Among them, the threshold setting is based on the Neyman-Pearson criterion, setting the false alarm probability , calculate the threshold , where Q is the Gaussian complementary cumulative distribution function.

[0065] Preferably, the maximum peak index of the spectral correlation density (SCD) output by the cyclostationary analysis is obtained as follows:

[0066] Calculate the cyclic autocorrelation function ,in is the candidate cycle frequency, is the input signal, is the time delay, represents the observation time, Indicates that the signal is delayed The self-similarity after Indicates that the extracted signal is at the cyclic frequency Periodicity Accelerate the calculation of the general correlation density by FFT , extract the maximum SCD peak value within the preset direct sequence spread spectrum signal frequency range , the maximum SCD peak value is processed with the ratio of the SCD noise standard deviation to obtain the maximum peak index;

[0067] Among them, the threshold setting sets the SCD peak threshold based on the noise statistical characteristics .

[0068] Preferably, the specific steps of outputting the decision value of the weighted average decision model and comparing it with the absolute decision value threshold and then outputting the decision result are:

[0069] The decision value output by the weighted average decision module is compared with the preset decision value threshold. If the decision value is within the decision value threshold, it indicates that a DSSS signal is present. If the decision value is outside the decision threshold, it indicates that no DSSS signal is present.

[0070] Preferably, the signal basic data information A includes broadband matching, transmission distance, transmission power, noise floor and energy detection method observation time, and the signal basic data information BB includes cyclic frequency accuracy, pseudo code length, symbol rate and cyclostationary analysis observation time, and the cyclostationary analysis observation time is consistent with the energy detection method observation time;

[0071] Broadband matching refers to the degree of matching between the signal bandwidth and the preset bandwidth of the receiving device. It can be calculated by the ratio of the actual bandwidth to the ideal bandwidth. It is used to reflect the integrity of signal reception and the effectiveness of energy detection.

[0072] Transmission distance refers to the spatial distance between the signal transmitter and the receiving device. It can be measured by receiving signal strength or arrival time difference, and is used to characterize the path loss characteristics during signal propagation.

[0073] Transmitted power refers to the electromagnetic wave power output by the signal source. It can be achieved by using a power meter or inversely calculating the received signal strength. It is used to measure the degree of energy attenuation of the signal in the propagation environment.

[0074] The noise floor refers to the power level of background noise in the receiving environment. It can be achieved by measuring the output power of the receiving device when there is no signal. It is used to determine the signal-to-noise ratio threshold for energy detection.

[0075] The observation time of the energy detection method refers to the length of time the receiving device accumulates signal energy. It can be implemented using a fixed time window or an adaptive time window setting to control the detection sensitivity of the energy statistics.

[0076] Cycle frequency accuracy refers to the estimated error range of the cycle frequency in cycle spectrum analysis, which can be achieved by setting the cycle frequency search step or spectrum resolution parameters, and is used to affect the accuracy of cycle feature extraction;

[0077] The pseudo-code length refers to the sequence length of the pseudo-random code in the spread spectrum signal. It can be specifically calculated by the pseudo-code period or the number of code chips, and is used to determine the autocorrelation characteristics and cyclostationary characteristics of the spread spectrum signal.

[0078] The symbol rate refers to the transmission rate of data symbols, which can be achieved by calculating the inverse of the symbol period or the spectrum occupied bandwidth. It is used to reflect the time-frequency characteristics of the spread spectrum signal modulation structure.

[0079] Specifically, for the energy detection method, broadband matching, transmission distance, transmission power, noise floor, and observation time are selected as basic data information A. By quantifying signal bandwidth matching, propagation loss, environmental noise, and time accumulation effects, the input parameter system required for the energy detection method reliability model is established. For the cyclostationary analysis method, cyclic frequency accuracy, pseudocode length, symbol rate, and observation time are selected as basic data information B. By defining cyclic feature extraction accuracy, spread spectrum signal structure parameters, and spectral analysis time base, the parameter set of the cyclostationary analysis reliability model is constructed. The observation time of the two methods is synchronized to ensure that the energy statistics and spectral correlation density calculations are based on the same time base, avoiding parameter correlation breaks caused by time asynchrony, thereby maintaining the spatiotemporal consistency of dynamic weight allocation in the joint decision model.

[0080] Compared with existing technologies, traditional signal detection methods do not clearly define the parameter composition of basic data information, resulting in arbitrary selection of detection model input parameters. The energy detection method in existing technologies only considers a single parameter, signal-to-noise ratio, and ignores the impact of transmission distance on path loss. The cyclostationary analysis method does not establish a parameter system that associates pseudocode length with symbol rate, and the observation time of the two methods is set independently, resulting in data benchmark deviation. This solution systematically defines the specific parameters of the two types of basic data information and establishes a multi-dimensional parameter system covering signal propagation, receiving environment, signal structure, and time benchmark, thus solving the core problem of inaccurate detection model parameter selection.

[0081] Through the above technical solution, this application achieves a standardized definition of the parameters for building reliability models for the energy detection method and cyclostationary analysis method, ensuring that the model input parameters can fully reflect the signal propagation characteristics, environmental noise level, and signal structure characteristics. The synchronous setting of observation time eliminates the data correlation errors caused by time base differences in traditional methods, and provides a unified spatial and temporal data foundation for the dynamic weight calculation in the joint decision model, thereby significantly improving the parameter selection accuracy of the signal detection model and the reliability of the decision results.

[0082] Preferably, the signal basic data index information A is obtained in the following manner:

[0083] The broadband matching degree is calculated according to the formula , get the broadband matching index ,in is the broadband matching degree, For ideal broadband, is the sensitivity factor;

[0084] The transmission distance is calculated according to the formula , and then get the transmission distance index , where , is the half-attenuation distance, is the transmission distance, is the half-life coefficient;

[0085] The transmit power is calculated according to the formula , get the transmit power index ,in is the power sensitivity factor, is the transmit power, is the reference power;

[0086] The noise floor is calculated according to the formula , and obtain the noise floor index ,in is the baseline noise, is the noise floor;

[0087] The observation time is divided into , and obtain the observation time index ,in For time sensitivity, is the observation time of the energy detection method.

[0088] Specifically, for the five heterogeneous parameters involved in the energy detection method, namely, broadband matching, transmission distance, transmission power, noise floor, and observation time, nonlinear normalized models matching their physical characteristics are constructed. The broadband matching is processed by an exponential decay function. When the actual bandwidth is close to the ideal bandwidth, the exponent approaches 1. When the deviation increases, the exponent decreases nonlinearly. The sensitivity factor Control the attenuation speed. The transmission distance adopts a negative exponential model to simulate the path loss law of electromagnetic wave propagation, and the half-attenuation coefficient The transmission distance corresponds to the signal power attenuated to 50% of the initial value. The transmit power is normalized using the Sigmoid function. When the power exceeds the reference benchmark, the exponential saturates quickly to avoid nonlinear distortion caused by high power. The noise floor maps the actual base noise to the 0-1 range through the proportional function. The reference noise As a normalized benchmark. The observation time adopts a gradual saturation curve to simulate the effect of time accumulation on the signal detection probability, and the sensitivity coefficient Controls the time threshold required to reach maximum detection probability.

[0089] Compared to existing technologies, traditional methods use linear normalization for multi-source heterogeneous parameters, which fails to reflect the nonlinear characteristics of parameter influences. For example, the effect of transmission distance on signal strength follows the inverse square law, but linear normalization underestimates the degree of signal attenuation over long distances. This solution constructs a nonlinear mapping model that matches the physical laws of each parameter. This eliminates dimensional differences while preserving the nonlinear characteristics of the impact of parameter changes on detection results. This allows the normalized index information to more accurately reflect the actual environmental impact.

[0090] Through the above technical solution, this application solves the problem of insufficient normalization accuracy caused by parameter dimension differences in traditional detection methods, quantifies heterogeneous parameters such as broadband matching deviation, transmission distance change, and power fluctuation under a unified dimension, provides accurate standardized input for the subsequent construction of reliability models, and improves the robustness of direct sequence spread spectrum signal detection in complex electromagnetic environments.

[0091] Preferably, the signal basic data index information B is obtained in the following manner:

[0092] The cycle frequency accuracy is calculated according to the formula , get the cycle frequency accuracy index ,in is the error sensitivity coefficient, is the cycle frequency accuracy;

[0093] The pseudo code sequence length is calculated according to the formula , get the pseudo code sequence length index ,in is the pseudo code length attenuation coefficient, is the pseudo code sequence length;

[0094] The symbol rate is calculated according to the formula , get the symbol rate index ,in is the symbol rate, is the pseudo code rate;

[0095] The observation time is calculated according to the formula , and obtain the observation time index ,in is the observation time, is the reference observation time.

[0096] Specifically, in a complex electromagnetic environment, the detection reliability is improved by normalizing the four parameters. For the cyclic frequency accuracy, the high-frequency error is compressed to the [0,1) interval using a formula structure whose denominator includes the error sensitivity coefficient and the error product, reducing its negative impact on the detection model. For the pseudo-code sequence length, an exponential decay function is used to convert the long code sequence into an exponential value, effectively suppressing the calculation complexity of the spectral correlation density caused by the increase in code length. For the symbol rate, the ratio of the symbol rate to the pseudo-code rate is mapped to a periodic attenuation curve through the cosine square function, so that the model maintains high sensitivity when the symbol rate is close to the pseudo-code rate, and automatically reduces the weight when it deviates. For the observation time, the ratio of the sum of the numerator and denominator is used, and when the observation time is too short, the ratio is reduced by Provides benchmark compensation to avoid overfitting by increasing the denominator when the observation time is too long.

[0097] Compared with existing technologies, traditional cyclostationary analysis methods fail to normalize detection parameters, making the detection results susceptible to differences in parameter dimensions and extreme values. For example, when the length of the pseudo-code sequence exceeds a threshold, the spectral correlation density calculation time of traditional methods increases exponentially. However, this solution automatically suppresses the weight of long code sequences through an exponential decay function. When the symbol rate and pseudo-code rate mismatch occur, traditional methods cannot effectively distinguish signal characteristics from noise. However, this solution automatically reduces the contribution of abnormal symbol rates through a cosine square function.

[0098] Through the above technical solution, this application solves the problem of insufficient detection reliability caused by parameter sensitivity differences in complex electromagnetic environments using cyclostationary analysis. Through the directional normalization of four parameters, the computational burden of long pseudo-code sequences is suppressed under low signal-to-noise ratio conditions, feature extraction stability is maintained when the symbol rate fluctuates, and spectral correlation density distortion is avoided when observation time is insufficient, ultimately improving the accuracy and real-time performance of direct sequence spread spectrum signal detection.

[0099] Preferably, the reliability of the energy detection method is obtained in the following manner:

[0100] The broadband matching index, transmission distance index, transmission power index, noise floor index and energy detection method observation time index are introduced into the constructed energy detection method reliability model, and the energy detection method reliability is output. The energy detection method reliability model is expressed as:

[0101]

[0102] in, Represents the reliability of the energy detection method, represents the broadband matching index, represents the transmission distance index, represents the transmit power index, represents the noise floor index, represents the energy detection method observation time index, represents the signal-to-noise ratio, represents the critical signal-to-noise ratio, is the signal-to-noise ratio sensitivity.

[0103] Preferably, the cyclostationary analysis reliability is obtained in the following manner:

[0104] The cyclic frequency accuracy index, pseudo code length index, symbol rate index and cyclostationary analysis observation time index are introduced into the constructed cyclostationary analysis reliability model, and the cyclostationary analysis reliability is output. The cyclostationary analysis reliability model is expressed as:

[0105]

[0106] in, represents the reliability of cyclostationary analysis, Cycle frequency accuracy index, represents the symbol rate index, represents the cyclostationary analysis observation time index, represents the minimum value among the symbol rate index, the cyclostationary analysis observation time index, and the cyclic frequency accuracy index, Indicates the pseudo code sequence length index.

[0107] Specifically, in complex electromagnetic environments, when errors exist in the cyclic frequency estimation, the reliability is automatically reduced through the denominator of the cyclic frequency accuracy index. When the symbol rate and the pseudo-code rate do not match, the cosine square function causes the symbol rate exponent to decay rapidly. When the observation time is insufficient, the observation time exponent limits the reliability improvement through the reference time term in the denominator. Taking the minimum value of these three can quickly identify the weak link effect of key parameters, ensuring that reliability is significantly reduced when any parameter does not meet the standard. At the same time, the pseudo-code sequence length exponent weakens the negative impact of long pseudo-code on reliability through exponential decay. Finally, the influence of each parameter is coupled through a product operation to output the final reliability value.

[0108] Compared to existing technologies, traditional cyclostationary analysis methods lack a multi-parameter coupled reliability assessment mechanism and are unable to distinguish differences in parameter sensitivity across different scenarios. Existing methods typically use fixed thresholds or single parameters to evaluate detection results, which can easily lead to misjudgments when the symbol rate changes suddenly or when observation time is insufficient. This solution, by constructing a dynamic parameter correlation model, accurately quantifies key influencing factors while maintaining computational efficiency, addressing the reliability fluctuations caused by differences in parameter sensitivity.

[0109] Through the above technical solution, this application can automatically suppress interference caused by cyclic frequency errors in low signal-to-noise ratio environments, trigger a reliability degradation mechanism when the symbol rate and pseudo-code rate mismatch occur, and dynamically weaken the negative impact of long pseudo-code signals, effectively improving the stability of cyclostationary analysis results in complex electromagnetic environments. This model provides an accurate reliability basis for subsequent weight allocation, enabling the detection system to adaptively adjust its decision strategy based on real-time environmental parameters.

[0110] Preferably, the weight distribution model is expressed as:

[0111]

[0112] in, represents the energy statistic weight, represents the maximum peak index weight, represents the reliability of cyclostationary analysis, Indicates the reliability of the energy detection method.

[0113] Specifically, the weight allocation model emphasizes the synergistic effect of the reliability of the energy detection method and the reliability of the cyclostationary analysis through the product term. When the two reliabilities are high, the product term in the numerator dominates, and the weight allocation model will give the energy statistic a higher weight, strengthening the decision-making position of the energy detection method in the judgment; when either reliability decreases, the complementary term in the denominator will suppress the interference of unreliable detection results on the weight allocation, thereby dynamically balancing the contributions of the two methods. For example, in a low signal-to-noise ratio scenario, if the reliability of the energy detection method decreases and the reliability of the cyclostationary analysis remains stable, the model will automatically reduce the weight of the energy statistic and increase the weight of the maximum peak index, ensuring that the judgment result relies more on the feature extraction capability of the cyclostationary analysis method. This process does not require manual intervention, and mathematical constraints ensure that the sum of the weights is 1, avoiding the problem of insufficient adaptability caused by subjective weight setting.

[0114] Compared to existing technologies, traditional methods typically use fixed weights or weight allocation based on empirical rules, which cannot dynamically adjust to changes in signal characteristics and environmental parameters. For example, fixed weight schemes may overly rely on energy detection, which is susceptible to interference, in high-noise environments, while threshold-based rule-based adjustment schemes have difficulty handling situations where reliability continuously varies. This solution establishes a nonlinear mapping relationship between reliability and weights, and uses a mathematical model to automatically generate weights, addressing the lack of robustness in decisions caused by fixed weights in traditional methods.

[0115] Through the above technical solution, this application realizes the adaptive allocation of weights of energy detection method and cyclostationary analysis method, and solves the problem of misjudgment caused by the reliability fluctuation of a single detection method in a complex electromagnetic environment. In scenarios with low signal-to-noise ratio or fuzzy signal characteristics, the model can suppress the weight of unreliable detection results and improve the accuracy of the judgment results; when the signal-to-noise ratio is high or the signal characteristics are significant, the model can strengthen the synergy of the two methods and enhance the stability of detection. This solution replaces manual experience intervention with a mathematically driven approach, significantly improving the adaptability of the signal detection system to dynamic environments.

[0116] It should be noted that, in this document, relational terms such as first and second, etc., are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that includes a list of elements includes not only those elements but also other elements not explicitly listed, or elements inherent to such process, method, article, or apparatus.

[0117] While embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions, and variations may be made to these embodiments without departing from the principles and spirit of the invention, and that the scope of the invention is defined by the appended claims and their equivalents.

Claims

1. A direct sequence spread spectrum signal detection method based on GPU parallelism, characterized in that: The following steps are involved: Obtain signal basic data information A under the energy detection method and signal basic data information B under the cyclostationary analysis method and perform normalization processing on the two to obtain signal basic data index information A and signal basic data index information B; Based on the signal basic data index information A, an energy detection method reliability model is constructed and the signal basic data index information A is imported to obtain the reliability of the energy detection method; According to the signal basic data index information B, a cyclostationary analysis reliability model is constructed and the signal basic data index information B is imported to obtain the cyclostationary analysis reliability; Based on the energy statistics output by the energy detection method and the maximum peak index of the spectral correlation density output by the cyclostationary analysis, a weighted average decision model is constructed to output the decision value, which is then compared with the absolute value threshold to output the decision result. Based on the reliability of energy detection method and cyclostationary analysis, a weight distribution model is constructed to output energy statistics and the dynamic weight of the maximum peak index in the weighted average decision model; The weight distribution model is expressed as: ； in, represents the energy statistic weight, represents the maximum peak index weight, represents the reliability of cyclostationary analysis, Indicates the reliability of the energy detection method; The weighted average decision model is reconstructed according to the dynamic weights of the energy statistics and the maximum peak index in the weighted average decision model; The specific steps of the weighted average decision model outputting the decision value and comparing it with the absolute value threshold to output the decision result are as follows: The decision value output by the weighted average decision module is compared with the preset decision value threshold. If the decision value is within the decision value threshold, it indicates that a DSSS signal is present. If the decision value is outside the decision threshold, it indicates that no DSSS signal is present.

2. The direct sequence spread spectrum signal detection method based on GPU parallelism according to claim 1, characterized in that: The weighted average decision model is expressed as: ； in, Indicates the decision value, represents the energy statistic weight, represents the maximum peak index weight, represents the energy statistic, Indicates the maximum peak index.

3. The direct sequence spread spectrum signal detection method based on GPU parallelism according to claim 2, characterized in that: The signal basic data information A includes broadband matching degree, transmission distance, transmission power, noise floor and energy detection method observation time, and the signal basic data information B includes cyclic frequency accuracy, pseudo code length, symbol rate and cyclostationary analysis observation time. The cyclostationary analysis observation time is consistent with the energy detection method observation time.

4. The direct sequence spread spectrum signal detection method based on GPU parallelism according to claim 1, characterized in that: The signal basic data index information A is obtained in the following manner: The broadband matching degree is calculated according to the formula , get the broadband matching index ,in is the broadband matching degree, For ideal broadband, is the sensitivity factor; The transmission distance is calculated according to the formula , and then get the transmission distance index , where , is the half-attenuation distance, is the transmission distance; The transmit power is calculated according to the formula , get the transmit power index ,in is the power sensitivity factor, is the transmit power, is the reference power; The noise floor is calculated according to the formula , and obtain the noise floor index ,in is the baseline noise, is the noise floor; The observation time is divided into , and obtain the observation time index ,in For time sensitivity, is the observation time of the energy detection method.

5. The direct sequence spread spectrum signal detection method based on GPU parallelism according to claim 1, characterized in that: The signal basic data index information B is obtained in the following manner: The cycle frequency accuracy is calculated according to the formula , get the cycle frequency accuracy index ,in is the error sensitivity coefficient, is the cycle frequency accuracy; The pseudo code sequence length is calculated according to the formula , get the pseudo code sequence length index ,in is the pseudo code length attenuation coefficient, is the pseudo code sequence length; The symbol rate is calculated according to the formula , get the symbol rate index ,in is the symbol rate, is the pseudo code rate; The observation time is calculated according to the formula , and obtain the observation time index ,in is the observation time, is the reference observation time.

6. The direct sequence spread spectrum signal detection method based on GPU parallelism according to claim 1, characterized in that: The reliability of the energy detection method is obtained as follows: The broadband matching index, transmission distance index, transmission power index, noise floor index and energy detection method observation time index are introduced into the constructed energy detection method reliability model, and the energy detection method reliability is output. The energy detection method reliability model is expressed as: ； in, Represents the reliability of the energy detection method, represents the broadband matching index, represents the transmission distance index, represents the transmit power index, represents the noise floor index, represents the energy detection method observation time index, represents the signal-to-noise ratio, represents the critical signal-to-noise ratio, is the signal-to-noise ratio sensitivity.

7. The direct sequence spread spectrum signal detection method based on GPU parallelism according to claim 1, characterized in that: The cyclostationary analysis reliability is obtained in the following way: The cyclic frequency accuracy index, pseudo code length index, symbol rate index and cyclostationary analysis observation time index are introduced into the constructed cyclostationary analysis reliability model, and the cyclostationary analysis reliability is output. The cyclostationary analysis reliability model is expressed as: ； in, represents the reliability of cyclostationary analysis, Cycle frequency accuracy index, represents the symbol rate index, represents the cyclostationary analysis observation time index, represents the minimum value among the symbol rate index, the cyclostationary analysis observation time index, and the cyclic frequency accuracy index, Indicates the pseudo code sequence length index.

8. The direct sequence spread spectrum signal detection method based on GPU parallelism according to claim 1, characterized in that: The energy statistics output by the energy detection method are obtained as follows: Divide the input signal into N frames and calculate the energy of each frame. ,in For the i-th frame signal, use the first M frames to obtain the noise power , the energy statistics are obtained by ratio processing of each frame energy and noise power.

9. The direct sequence spread spectrum signal detection method based on GPU parallelism according to claim 1, characterized in that: The maximum peak index of the spectral correlation density output by the cyclostationary analysis is obtained as follows: Calculate the cyclic autocorrelation function ,in is the candidate cycle frequency, is the input signal, is the time delay, represents the observation time, Indicates that the signal is delayed The self-similarity after Indicates that the extracted signal is at the cyclic frequency The periodicity at the position is calculated by FFT to accelerate the calculation of the general correlation density. , extract the maximum SCD peak value within the preset direct sequence spread spectrum signal frequency range , the maximum SCD peak value is ratioed to the SCD noise standard deviation to obtain the maximum peak index.

Images