153 results about "Zero padding" patented technology

The invention relates to the technical field of charging load prediction, in particular to an electric vehicle charging load prediction method and system combining sequence decomposition and feature extraction and a storage medium. Comprising the steps of performing variational mode decomposition on time sequence load data to obtain a plurality of intrinsic mode functions; performing zero padding expansion on each intrinsic mode function, and then converting the intrinsic mode function to a frequency domain; enhancing the frequency domain data through a multi-layer perceptron based on complex domain calculation; constructing a prediction model based on a complex field attention module and a residual module, inputting the enhanced spectrum data into the prediction model, extracting a data characteristic frequency combination through the complex field attention module and the residual module, and carrying out weighted summation on all intrinsic mode functions as the output of the prediction model; and calculating a difference value between real data and prediction data according to a mean square error of a time sequence and an average absolute value error of a frequency spectrum to obtain time domain and frequency domain mixed loss, and updating prediction model parameters. According to the invention, the problem of large prediction error under the influence of random load fluctuation is solved.

The invention discloses a joint inversion method for near-well characteristics of a fracturing crack based on high-frequency tube waves and water attack, and relates to the technical field of fracturing crack inversion. The method comprises the following steps: establishing a tube wave shaft crack model; carrying out segmentation, zero padding and Fourier transform on the pre-pressing pipe wave excitation signal, then substituting the pre-pressing pipe wave excitation signal into a pipe wave shaft crack model for inversion to obtain an optimal pipe wave quality factor and an optimal pipe wave velocity, and correcting the model; the depth of a fracturing main liquid inlet point is calculated based on a pump stopping water attack curve; and performing segmentation, zero padding and Fourier transform on the pressed pipe wave excitation signal, and then substituting the pressed pipe wave excitation signal into the corrected pipe wave wellbore crack model for inversion. According to the method, the echo signal obtained through calculation is closer to the real situation, the depth of the main liquid inlet point is determined in advance, accurate and rapid inversion of the tube wave is facilitated, and meanwhile the inversion time can be greatly saved.

Methods and apparatus for communicating zero-padded signals related to zero padding of a transmission symbol that is to be transmitted in a wireless communication network. Signaling indicative of such zero padding is communicated between communication devices, and a zero-padded transmission symbol to which the zero padding has been applied, as well as a cyclic prefix preceding the zero-padded transmission symbol, are transmitted by a transmitting device and received by a receiving device. In some embodiments, the zero padding includes both head zero padding at a head of the transmission symbol and tail zero padding at a tail of the transmission symbol. Zero padding may also be applied to a reference signal sequence upon which a reference signal for channel estimation is based.

The invention provides an inhaul cable self-vibration fundamental frequency automatic identification method and system based on a harmonic product spectrum, and relates to the technical field of structure monitoring, and the method comprises the steps: receiving vibration acceleration data of an inhaul cable, carrying out the windowing processing of the vibration acceleration data through employing a Hanning window, and generating a windowed signal; carrying out fast Fourier transform on the windowed signal, and improving the frequency resolution by implementing zero filling to obtain an initial frequency spectrum; smoothing the initial frequency spectrum to form a smooth frequency spectrum; performing passband clipping on the smooth frequency spectrum to zero the passband outer amplitude to form a clipped frequency spectrum; calculating the signal-to-noise ratio of the frequency spectrum after cutting, and comparing the signal-to-noise ratio with a preset threshold value; if the signal-to-noise ratio is lower than a preset threshold value, judging that the signal quality is unqualified and terminating the processing flow; and if the signal-to-noise ratio is equal to or higher than the preset threshold value, continuing the normal processing flow. According to the invention, automatic and high-precision identification of the self-vibration fundamental frequency of the inhaul cable can be realized.

The invention provides a sensor data recovery method and system based on mask perception space-time modeling, and belongs to the technical field of Internet of Things and intelligent sensing data processing. The method comprises the following steps: modeling a missing position into a trainable potential vector representation through an adaptive missing representation module, and avoiding noise introduced by traditional zero filling or mean filling; a missing perception space-time decoding module is adopted, and a dynamic weight distribution mechanism of mask constraint is combined in the decoding process, so that information is effectively prevented from being excessively smooth; designing a space-time dual-channel feature aggregation unit, and capturing spatial dependence and time sequence dependence between the sensors at the same time; and finally realizing accurate completion of a large-scale space-time sensor matrix through an output recovery unit. According to the method, the accuracy and robustness of sensor data restoration can be remarkably improved, the reliability of subsequent monitoring, prediction and anomaly detection is enhanced, and the method has wide engineering application value.

The invention provides a non-cooperative bistatic radar real-time signal processing method based on a heterogeneous platform. The non-cooperative bistatic radar real-time signal processing method comprises the steps that the heterogeneous platform is built; obtaining CPU fixed memory areas corresponding to the number of GPU blocks, and creating a data interaction thread to realize parallel computing; performing data analysis through frame header identifier retrieval to obtain transmitting signal parameter information, multi-channel echo data and reference channel data, and performing calculation to obtain echo beam data; zero padding is carried out on the echo beam data and the reference channel data, matched filtering is carried out on the echo beam data after zero padding, and a distance one-dimensional image is obtained; performing fast Fourier transform calculation in the slow time dimension to obtain data after coherent accumulation; obtaining a target trace point; and sorting the target plots, grouping the sorted plots, screening plot groups by using an M / N algorithm, fusing each screened plot group into a single plot, and outputting a final target plot array. According to the invention, signals can be processed in real time, and target trace points can be output.

The application belongs to the field of digital signal processing, and particularly relates to a linear frequency modulation signal parameter estimation method under a low signal-to-noise ratio condition, which comprises the following steps: collecting signals and discretizing the signals; dividing the discrete sequence into sub-signals by framing, performing an FFT operation after a windowing and zero padding operation, and obtaining a frame sequence vector; searching for a maximum spectral peak of the spectrum of each frame of signals, and correcting the maximum spectral peak by using a cubic spline interpolation method to obtain a corresponding frequency; performing an equal-step distance and successive rotation scanning of two groups of vectors in a plane; synchronously projecting frequency vectors at different rotation angles, dividing a frequency band, and recording a frequency quantity statistical result; eliminating a trend item by using a median filtering, obtaining a maximum value of a statistical quantity, and then obtaining a preliminary estimation value of a parameter; and obtaining a final estimation result of an initial frequency and a frequency modulation rate of the linear frequency modulation signal by using a least square method fitting. The application overcomes the problems of great difficulty, low precision and complex calculation of the linear frequency modulation signal parameter estimation under the low signal-to-noise ratio condition in the prior art, and has good engineering practicability.

The invention discloses a flexible spectral efficiency chirp multi-carrier signal generation method based on zero insertion and zero padding, and belongs to the technical field of chirp multi-carrier signal generation. The invention solves the problems that the spectrum efficiency of the existing non-orthogonal chirp multicarrier multiplexing technology is limited, and the applicability and universality of the structure of the existing non-orthogonal chirp multicarrier multiplexing system are limited. According to the invention, the operations of zero insertion and zero padding are added after the first phase rotation is carried out on the transmitting and receiving ends, and the number of points of DFT and IDFT can be adjusted by controlling the number of zero insertion and zero padding so as to adjust the bandwidth compression factor in the chirp multicarrier multiplexing system, thereby providing flexible and variable spectrum efficiency, and improving the performance of the chirp multicarrier multiplexing system. The problem that the spectrum efficiency range is limited due to the fact that the number of operation points must be an integer and is the power of 2 in the chirp multi-carrier signal generation process is solved, the structural application range of a chirp multi-carrier multiplexing system is further expanded, and the chirp multi-carrier multiplexing system has higher universality. The method provided by the invention can be applied to chirp multi-carrier signal generation.

Techniques, systems and methods are provided for performing convolution operations on real-valued input data using implicit zero-padding for dealiasing. A set of real-valued input data intended for convolution is embedded into a complex array, reorganizing the real-valued input into a format amenable to fast Fourier transform (FFT) processing. FFT operations are performed on the complex array to generate intermediate frequency-domain results. Pointwise operations on these intermediate results are then performed to integrate the intermediate results in a format that unifies the earlier-separated real-valued input data. The convolution output results matrix is extracted via inverse FFT operations performed on the results of these pointwise operations.

The disclosure provides a method performed by a UE in a wireless communication system, including receiving a channel state information (CSI) reporting configuration comprising a plurality of first sub-configurations, and determining and reporting CSI based on a second sub-configuration comprising N out of the plurality of the first sub-configurations, where N is greater than or equal to 1, wherein the CSI comprises one or two CSI parts, wherein a number of the CSI parts included in the CSI is determined based on N, and wherein at least one of the one or two CSI parts included in the CSI is determined based on a first order including at least one of an order of the N first sub-configurations in the second sub-configuration, an order of CSI parameters in the CSI, and an order between zero-padding bits and the CSI parameters in the CSI.

The invention provides a downlink channel recovery method for an asymmetric large-scale MIMO system under UPA. The method comprises the following steps: S1, according to the distribution of an asymmetric array, carrying out zero padding on a received signal, defining a mixed angle and constructing a virtual array; s2, designing an antenna selection algorithm on the premise of ensuring the UPA first row array aperture, and constructing a maximum aperture virtual array with as few antennas as possible; s3, respectively executing improved Newton orthogonal matching pursuit on receiving signals of the virtual array and the first row of arrays, and obtaining all possible path gains, azimuth angles and elevation angles until the gains, azimuth angles and elevation angles respectively meet threshold values; s4, combining and pairing the estimated path gain, elevation angle and azimuth angle according to related criteria; and S5, recovering a downlink channel by using the obtained path gain and angle parameter pair and the downlink steering vector. According to the method, part of channel information contained in the virtual array is fully utilized, and the downlink channel recovery performance is effectively improved while the number of received uplink radio frequency chains and the calculation complexity are remarkably reduced.

The invention discloses a CPU end radar data flow receiving exception processing method and device, and particularly relates to the field of radar signal processing. Comprising the following steps: performing error code check on each path of optical fiber radar data; reading the verified optical fiber radar data, and carrying out cutoff and packet loss detection; when it is determined that cutoff or packet loss occurs in the radar data of any optical fiber, reading of the optical fiber radar data is stopped, zero padding operation is carried out on the optical fiber radar data, and the use value of a queue of the cutoff optical fiber is obtained; when it is determined that the use value of the queue is larger than the use threshold value or the idle value is smaller than the idle threshold value, the radar data transmission queue of the current optical fiber is null, and reading of the radar data of the current optical fiber is recovered; and rearranging and splicing the read radar data of each path of optical fiber into the form of all beam data of a single CPI. And normal operation of the program can be ensured.

Systems, methods and computer software are disclosed for providing high resolution timing advance estimation based on Physical Random Access Channel (PRACH). An example method includes receiving a preamble signal r(n) having a predetermined sampling frequency and a predetermined length; correlating a down sampled version of the received preamble with a reference preamble sequence c(n) using an FFT method to provide correlation output Ryc; using a peak value P of the correlation output Ryc to detect a preamble ID and a timing advance at a resolution of 24Ts; zero padding sequences Y(k) and C(k) so that they have a predetermined length resulting in sequences Y_hat(k) and C_hat(k), which are 1024-point FFT of y(n) and c(n); performing a maximum likelihood estimation (MLE) to estimate a timing offset; and detecting a peak value out of the R_hat(m) and using a corresponding index Q to provide a timing advance with an accuracy of 2Ts.

Systems, methods and computer software are disclosed for providing high resolution timing advance estimation based on Physical Random Access Channel (PRACH). An example method includes receiving a preamble signal r(n) having a predetermined sampling frequency and a predetermined length; correlating a down sampled version of the received preamble with a reference preamble sequence c(n) using an FFT method to provide correlation output Ryc; using a peak value P of the correlation output Ryc to detect a preamble ID and a timing advance at a resolution of 24Ts; zero padding sequences Y(k) and C(k) so that they have a predetermined length resulting in sequences Y_hat(k) and C_hat(k), which are 1024-point FFT of y(n) and c(n); performing a maximum likelihood estimation (MLE) to estimate a timing offset; and detecting a peak value out of the R_hat(m) and using a corresponding index Q to provide a timing advance with an accuracy of 2Ts.

The invention provides a function beam forming improvement method and a related device, and relates to sound source recognition based on a microphone array. An array signal receiving model is established; calculating FB output according to the array signal receiving model; according to the array signal receiving model, calculating a PSF matrix with invariant power rising space transfer; determining spatial domain convolution information according to the FB output and the PSF matrix; converting into a wavenumber domain product through FFT (Fast Fourier Transform) operation under a periodic boundary condition according to the spatial domain convolution information; constructing a difference function according to the wavenumber domain product and the FB output; a difference function is calculated through a fast iteration shrinkage threshold method, the wavenumber leakage problem caused by a zero padding boundary is solved, the operation accuracy is improved, and the imaging performance is further improved; the main lobe of the point spread function is sharpened through exponential operation, and the applicability of the point spread function space transfer invariance hypothesis is expanded. The imaging spatial resolution and the dynamic range are improved, and the effective imaging area is expanded. The technical effect of avoiding performance loss in the sound source recognition process is achieved.

An intrusion detection method performed by an intrusion detection system using CRC includes receiving in-vehicle CAN data in units of frame, generating first feature information by extracting a CAN ID from the in-vehicle CAN data in units of frame and performing zero padding, generating second feature information by extracting a CRC field from the in-vehicle CAN data in units of frame and performing the zero padding, training a learning model by setting the first feature information and the second feature information as input data and labeling, as output data, whether the CAN data is attack data or normal data, and detecting whether the CAN data is normal data or attack data by extracting a CAN ID and a CRC field from the received in-vehicle CAN data when training of the learning model is completed and inputting the CAN ID and the CRC field to the learning model.

The present invention provides a generalized constant modulus signal generation method, device and communication equipment. The method comprises: obtaining a first matrix according to a single-user signal or a multi-user signal, the first matrix including M rows and N columns, where M and N are integers greater than 1; subjecting the M signals in each column to a first transform domain processing to obtain a first sequence, the first sequence including M signals; repeating the M signals in the first sequence J times and then performing point multiplication with a filter of length M×J to obtain a second sequence of length M×J, where J is an oversampling factor and J is an integer greater than 1; subjecting the second sequence to a second transform domain processing and zero padding processing to obtain a third sequence of length M×N, and cyclically shifting the third sequence to obtain the first signal; accumulating the first signals of all columns to obtain a second signal, and performing phase modulation on the second signal to obtain a generalized constant modulus signal, the generalized constant modulus signal occupying M symbols in the time domain and N subcarriers in the frequency domain.

The invention provides a channel estimation method based on a double-domain weighting-noise threshold, which comprises the following steps of: 1, time-frequency domain transformation: after receiving signals of a receiving end are subjected to time-frequency synchronization processing, extracting a pilot frequency sequence pilot, carrying out N-point FFT (Fast Fourier Transform) on the receiving signals, and converting the receiving signals to a frequency domain to obtain PILOT; performing N-point FFT (Fast Fourier Transform) on the time domain local pilot frequency sequence pilot to generate a frequency domain local pilot frequency PILOT; 2, calculating a frequency domain weighting coefficient: calculating the frequency domain weighting coefficient Coef according to the local pilot frequency PILOT and the current linear signal-to-noise ratio sn; step 3, channel estimation calculation: carrying out signal-to-noise ratio self-adaptive adjustment on the received signal of each frequency point through the operation; 4, inverse discrete Fourier transform: performing inverse Fourier transform on the channel estimation H after frequency domain weighting, and converting a frequency domain signal into a time domain signal to obtain h; step 5, extracting useful information: calculating and determining an optimal symbol length parameter L according to the maximum time delay extension of the system, performing optimization truncation on time domain data, retaining original values of the first L sampling points, endowing a zero value to subsequent sampling data exceeding the L point, and finally obtaining an optimized time domain result hL; 6, threshold judgment: within the range of the effective symbol length L, calculating a self-adaptive noise judgment threshold, implementing intelligent noise filtering processing on the truncated L-point data by using the threshold, eliminating a noise component through the threshold, and finally outputting a purified result hLR; and 7, zero padding: setting the length of the data block needing to be balanced as M, and setting the point number of channel estimation as N. When Mgt; when M is larger than N, time domain zero padding needs to be carried out on the hLR, the number of zero padding is M-N, and hLRZ is obtained; and 8, DFT: carrying out Fourier transform on the hLRZ to obtain a final complete channel estimation HLRZ.

Embodiments disclosed herein are directed to communicating time-critical ultra-low latency (ULL) data. An access point station (AP) communicates time-critical ULL data using aggregated MAC Protocol Data Unit (A-MPDU) preemption. When time-critical ULL data for a second associated STA (STA2) becomes available at a medium access control (MAC) layer of the AP during transmission of a physical layer protocol data unit (PPDU) to a first associate station (STA1), the AP may encode the time-critical ULL data in a new A-MPDU subframe for insertion before one of the A-MPDU subframes of the PPDU that has not yet been transmitted. The new A-MPDU subframe may be encoded to include zero-padding to set a size of the new A-MPDU subframe equal to a size of the A-MPDU subframe that has been preempted. The A-MPDU subframes 606 for STA1 may be encoded include a MAC address of the STA1 and the new A-MPDU subframe 608 may be encoded include a MAC address of the STA2.

A chirp multi-carrier transmission method based on windowing and truncation belongs to the field of wireless communication technology. The purpose of the present invention is to achieve efficient spectrum utilization while improving the bit error rate performance. The present invention first performs digital modulation and serial-to-parallel conversion on a random bit stream sequence, then performs a discrete Fresnel inverse transform on the conversion result, and after the discrete Fresnel inverse transform result passes through a time domain truncation window, the intercepted signal is converted into parallel and serial, and the parallel-to-serial conversion result is then sent through a channel. The receiving end performs serial-to-parallel conversion, zero padding, discrete Fresnel transformation, and inter-subcarrier interference elimination on the received signal, and then digitally demodulates the signal after interference elimination to restore the original bit stream data. The present invention saves transmission time and improves transmission efficiency by performing windowing and truncation in the time domain, and the bit error rate performance is better than that of existing systems. The method of the present invention can be applied to the field of wireless communication technology.

To provide a matrix multiplication hardware architecture.SOLUTION: A matrix multiplication hardware architecture includes: a reduction network including a tree topology with multiple levels formed by a plurality of reduction network nodes, each of which includes a data selector and two computation paths; and a chain formed by cascading a plurality of digital signal processing units (DSPs), where output ends of adjacent DSPs are respectively connected to the two computation paths of the same reduction network node in a first level of the tree topology, and outputs of the two computation paths are connected to a reduction network node in a higher level of the tree topology through data selectors. An addition tree of the FSB is changed into an addition chain adapted to the DSP structure, so that rear adders of the DSPs are multiplexed, the hardware utilization rate is improved, and the sign extension of higher-order bits can also be changed into zero padding.SELECTED DRAWING: Figure 2

The invention provides a bee colony unmanned aerial vehicle communication system based on an enhanced OTFS frame and a hybrid topology, and relates to the technical field of unmanned aerial vehicle cluster wireless communication, and the system comprises autonomous wireless communication modules which are deployed in all unmanned aerial vehicles and are used for generating and processing orthogonal time-frequency space OTFS modulation signals; a star-mesh-linear IP hybrid topology network dynamically constructed by unmanned aerial vehicle nodes supports autonomous optimization of a communication path according to channel quality; embedding an enhanced OTFS frame structure of reduced-order cyclic prefix (RCP) and reduced-order zero filling (RZP) pilot frequency; a ground control station GSC integrating an OTFS module and a remote controller RC establish a low-delay link through UDP multicast; the invention relates to a RadCom joint module with integrated radar sensing and communication functions. According to the invention, a single-point fault risk is eliminated through decentration topology, inter-symbol and inter-carrier interference caused by multipath interference is effectively suppressed by an enhanced OTFS frame, the anti-interference capability of a complex environment is improved by a dual-band frequency hopping spread spectrum technology, and target detection and data transmission millisecond-level synchronization are realized by a RadCom module.

This disclosure relates to precoding for polar codes and provides a method comprising: determining a first precoding matrix, the size of which is determined based on the message size of a message to be encoded by polar coding and the codeword size of the polar-coded codeword; determining a second precoding matrix based on the first precoding matrix and the information positions of a zero-padding vector; and performing precoding on the message based on the second precoding matrix to generate an output vector.

Methods and apparatus for communicating zero-padded signals related to zero padding of a transmission symbol that is to be transmitted in a wireless communication network. Signaling indicative of such zero padding is communicated between communication devices, and a zero-padded transmission symbol to which the zero padding has been applied, as well as a cyclic prefix preceding the zero-padded transmission symbol, are transmitted by a transmitting device and received by a receiving device. In some embodiments, the zero padding includes both head zero padding at a head of the transmission symbol and tail zero padding at a tail of the transmission symbol. Zero padding may also be applied to a reference signal sequence upon which a reference signal for channel estimation is based.

The present invention discloses a network traffic prediction method based on local downsampling and sparse frequency grouping, the method comprising: sequentially dividing an input sequence into P regions; performing downsampling operations within the regions; performing real sequence fast Fourier transform on the downsampling results; performing frequency truncation operations on the results; performing sparse frequency grouping and sparse frequency mixing operations; performing cross-region frequency point prediction; performing zero padding and inverse real sequence fast Fourier transform irFFT; and performing sequence rearrangement operations to obtain an output sequence.

A method by a wireless device includes receiving, from the network node, a first message comprising a 3 bit indicator associated with a first number of Physical Uplink Shared Channel (PUSCH) repetitions to be used by the wireless device for a first instance of a PUSCH transmission. The wireless device receives a second message comprising a 2 bit indicator for determining a second number of PUSCH repetitions to be used by the wireless device for one or more subsequent instances of the PUSCH transmission based on said second number of PUSCH repetitions. The wireless device applies a zero-bit padding to the 2 bit indicator to obtain an updated 3 bit indicator used to communicate from physical layer to higher layers the second number of PUSCH repetitions to be used and sends the one or more subsequent PUSCH transmissions based on the second number of PUSCH repetitions.

The invention discloses a target harmonic signal complex frequency estimation method based on frequency domain least square, and belongs to the technical field of signal processing, and the method comprises the steps: obtaining a geophysical observation time sequence containing a target harmonic signal; windowing Fourier transform is adopted to obtain an initial power spectrum, and a spectrum peak of the target harmonic signal is identified; selecting power spectrum values at a plurality of frequency points near the target spectrum peak, and obtaining complex frequency parameters of the target harmonic signal by using linear least square fitting; zero-filling expansion is adopted to increase effective frequency points of a power spectrum, and the estimation precision of complex frequency parameters is improved through multiple least squares. According to the method, Fourier transform and linear least square fitting are combined, high-precision estimation of the complex frequency of the target harmonic signal can be achieved, the method is simple, efficient and easy to implement, and the method is successfully applied to scenes such as earth normal mode frequency extraction and ocean partial tide parameter inversion and can be expanded to any time sequence signal analysis such as communication and biomedicine.

This application belongs to the field of interconnection network reliability and fault diagnosis technology, and discloses a method for intermittent fault diagnosis of low-bandwidth long interconnection networks based on graph attention mechanism. Targeting the hierarchical recursive structure and high connectivity of the network, under the PMC fault diagnosis model, a multi-round testing strategy is used to obtain test symptoms within the node's neighborhood. A feature vector is constructed for each node using local statistical feature extraction methods, and preprocessed using zero-padding and noise enhancement techniques. Finally, a graph attention network model is constructed, and the importance weights of neighboring node test results are dynamically learned using the attention mechanism to achieve accurate diagnosis of node fault states in the network. This application leverages the powerful local information aggregation capability of graph attention networks to overcome the diagnostic limitations of traditional algorithms, maintaining high diagnostic accuracy and robustness even with high fault rates, incomplete test symptoms, and large network sizes.

A method and apparatus for transmitting a zero-filled signal related to zero filling of a transmission symbol to be transmitted in a wireless communication network. Signaling indicating such zero padding is transmitted between communication devices, and a zero-padding transmission symbol to which the zero padding has been applied and a cyclic prefix prior to the zero-padding transmission symbol are transmitted by a transmitting device and received by a receiving device. In some embodiments, the zero padding includes a head zero padding at a head of a transmission symbol and a tail zero padding at a tail of the transmission symbol. Zero filling may also be applied to a reference signal sequence on which a reference signal for channel estimation is based.

The present invention discloses a distributed OFDM base station cooperative sensing sequence design method, comprising the following steps: setting parameters required for the problem; generating a signal sequence, a Fourier transform matrix, a frequency domain subcarrier selection matrix, and a zero padding matrix; generating a sidelobe level coefficient matrix and an autocorrelation ratio coefficient matrix; using an alternating iterative algorithm to obtain the optimal solution to the problem, solving the waveform optimal solution and limiting the peak-to-average power ratio, updating the matched filter and performing energy limitation, updating the frequency domain waveform, updating the dual variable, and calculating various residuals, objective functions, and waveform autocorrelation characteristics. After the iteration, the optimal waveform sequence is output, and various information such as residuals, objective functions, and autocorrelation characteristics are calculated. The present invention uses an alternating iterative algorithm to obtain the optimal solution to the problem, thereby maximally suppressing the autocorrelation sidelobe characteristics of multi-station waveforms and reducing interference of clutter on the communication system.