CN116566471A - GPU-based multi-antenna array satellite signal synthesis method, device and medium - Google Patents
GPU-based multi-antenna array satellite signal synthesis method, device and medium Download PDFInfo
- Publication number
- CN116566471A CN116566471A CN202310548816.7A CN202310548816A CN116566471A CN 116566471 A CN116566471 A CN 116566471A CN 202310548816 A CN202310548816 A CN 202310548816A CN 116566471 A CN116566471 A CN 116566471A
- Authority
- CN
- China
- Prior art keywords
- signal
- signals
- processing
- sub
- sampling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000001308 synthesis method Methods 0.000 title abstract description 21
- 238000012545 processing Methods 0.000 claims abstract description 114
- 238000005070 sampling Methods 0.000 claims abstract description 84
- 238000000034 method Methods 0.000 claims abstract description 62
- 238000007781 pre-processing Methods 0.000 claims abstract description 45
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 10
- 238000006243 chemical reaction Methods 0.000 claims description 15
- 125000004122 cyclic group Chemical group 0.000 claims description 13
- 238000000605 extraction Methods 0.000 claims description 13
- 238000005457 optimization Methods 0.000 claims description 11
- 239000002131 composite material Substances 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 10
- 238000013461 design Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 7
- 230000008520 organization Effects 0.000 claims description 6
- HPTJABJPZMULFH-UHFFFAOYSA-N 12-[(Cyclohexylcarbamoyl)amino]dodecanoic acid Chemical compound OC(=O)CCCCCCCCCCCNC(=O)NC1CCCCC1 HPTJABJPZMULFH-UHFFFAOYSA-N 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 abstract description 16
- 238000003786 synthesis reaction Methods 0.000 abstract description 16
- 230000008569 process Effects 0.000 abstract description 4
- 238000004891 communication Methods 0.000 description 18
- 238000001228 spectrum Methods 0.000 description 10
- 230000005540 biological transmission Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000005236 sound signal Effects 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000013475 authorization Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18513—Transmission in a satellite or space-based system
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
- G06F18/15—Statistical pre-processing, e.g. techniques for normalisation or restoring missing data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18519—Operations control, administration or maintenance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Astronomy & Astrophysics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Signal Processing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Probability & Statistics with Applications (AREA)
- Radio Relay Systems (AREA)
Abstract
The invention provides a method, a device, a medium and electronic equipment for synthesizing multi-antenna array satellite signals based on a Graphic Processing Unit (GPU). The invention discloses a GPU-based multi-antenna array satellite signal synthesis method, which comprises the following steps: acquiring sampling signals of the same signal source through multiple paths of antennas; carrying out first pretreatment on the sub-channels and determining the center frequency of the sub-channels; performing second preprocessing on the sampling signals according to the center frequency of the sub-channel to determine baseband signals; and carrying out first-stage channelization and second-stage channelization on the baseband signal, and determining a target synthesized signal, wherein the target synthesized signal comprises narrowband signals with various output sampling rates. The invention can separate the sub-channels and the narrow-band signals in the satellite signals in the signal synthesis stage, realize the matching of various target output sampling rates of the narrow-band signals, process the narrow-band signals into target synthesis signals, and improve the signal synthesis performance and the signal-to-noise ratio gain capability.
Description
Technical Field
The invention relates to the technical field of communication, in particular to a method, a device, a medium and electronic equipment for synthesizing multi-antenna array satellite signals based on a Graphic Processing Unit (GPU).
Background
With the development of satellite communication systems in recent years, the demand for high-rate satellite communication is also increasing, and the modulation order is required to be increased. Because of the increase of the modulation order, the ground receiving end needs to obtain a higher signal-to-noise ratio of the received signal so as to reach the lowest demodulation threshold, the signal synthesis method widely adopted at present is full spectrum synthesis so as to improve the signal-to-noise ratio of the received signal, and the full spectrum synthesis method comprises time domain full spectrum synthesis and frequency domain full spectrum synthesis, wherein the full spectrum synthesis method has the problems of large related operation amount between signals, large hardware cost and unsatisfactory synthesis gain caused by the fact that the time delay compensation precision is limited by the sampling rate; the frequency domain synthesis method is widely applied in the field of deep space communication, is suitable for synthesizing broadband signals, and is not suitable for synthesizing satellite signals.
The satellite communication at present widely adopts a Frequency Division Multiple Access (FDMA) access mode to allocate different frequency band resources for users, a plurality of sub-channels exist in the bandwidth in the signal, a plurality of narrow-band signals exist in the sub-channels, and the narrow-band signals are needed by the input of processing modules such as back-end demodulation decoding and the like. Under non-cooperative satellite communication, the frequency points of the sub-channels are unknown and are unevenly distributed, and different sub-channel pairs have different target output sampling rates. The receiving end has very high real-time requirements on signal processing, the traditional full-spectrum synthesis method cannot flexibly configure sub-channel frequency points needing to be separated, cannot meet the requirement of multiple output sampling rates, and is not suitable for synthesizing satellite signals under non-cooperative communication.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a method, a device, a medium and electronic equipment for synthesizing multi-antenna array satellite signals based on a GPU.
In order to achieve the above object, according to a first aspect of the present invention, a GPU-based multi-antenna array satellite signal synthesis method is provided.
Optionally, sampling signals of the same signal source are obtained through multiple paths of antennas;
performing first preprocessing on a sub-channel, and determining the center frequency of the sub-channel;
performing second preprocessing on the sampling signals according to the center frequency of the sub-channel to determine baseband signals;
and carrying out first-stage channelization and second-stage channelization on the baseband signal to determine a target synthesized signal, wherein the target synthesized signal comprises narrowband signals with various output sampling rates.
Optionally, the baseband signal performs a first stage of channelization and a second stage of channelization, to determine a target composite signal, including:
performing linear shift processing on the baseband signal through first-stage channelization, and determining a signal after the linear shift processing;
performing cyclic shift processing on the baseband signal subjected to the linear shift processing through second-stage channelization, and determining a signal subjected to the cyclic shift processing;
and performing discrete Fourier transform processing on the signal subjected to the cyclic shift processing to determine a target synthesized signal.
Optionally, the method further comprises:
after the step of carrying out first-stage channelization and second-stage channelization on the baseband signals, carrying out phase difference processing on narrowband signals at the same frequency point of the multipath antenna to determine a phase difference;
fitting the phase difference based on a least square method to determine a residual time delay and a phase deviation value;
performing decimal time delay compensation processing on the narrowband signals at the same frequency point of the multipath antenna according to the residual time delay;
and carrying out phase compensation processing on the narrowband signals at the same frequency point of the multipath antenna according to the phase offset value.
Optionally, the second preprocessing is performed on the sampling signal according to the center frequency of the sub-channel, and determining a baseband signal includes:
performing down-conversion processing on the sampling signal by taking the center frequency of the sub-channel as a target;
performing P times interpolation processing and Q times extraction processing on the sampling signal subjected to the down-conversion processing;
and performing low-pass filtering processing on the sampling signals subjected to the P times interpolation processing and the Q times extraction processing to determine baseband signals.
Optionally, the method further comprises:
and carrying out second preprocessing on the sampling signals according to the center frequency of the sub-channels, and carrying out integer time delay compensation processing on the baseband signals after the step of determining the baseband signals.
Optionally, the method further comprises:
and transplanting the design method to a CPU-GPU heterogeneous platform for execution through CUDA programming.
Optionally, the method further comprises:
when the design method is executed on the CPU-GPU heterogeneous platform, parallel processing is carried out in a mode of asynchronous copying, page locking memory optimization, on-chip programmable memory optimization and thread organization structure.
According to a second aspect of the present invention, there is provided a GPU-based multi-antenna array satellite signal synthesis apparatus, comprising:
the sampling module is used for acquiring sampling signals of the same signal source through multiple paths of antennas;
the first preprocessing module is used for carrying out first preprocessing on the sub-channels and determining the center frequency of the sub-channels;
the second preprocessing module is used for carrying out second preprocessing on the sampling signals according to the center frequency of the sub-channel to determine baseband signals;
and the determining module is used for carrying out first-stage channelization and second-stage channelization on the baseband signal to determine a target synthesized signal, wherein the target synthesized signal comprises narrowband signals with various output sampling rates.
According to a third aspect of the present invention there is provided a non-transitory computer readable storage medium having stored thereon a computer program, characterized in that the program when executed by a processor implements the steps of the method provided by the first aspect of the present invention.
According to a fourth aspect of the present invention, there is provided an electronic device comprising:
a memory having a computer program stored thereon;
a processor for executing the computer program in the memory to implement the steps of any of the methods provided in the first aspect of the present invention.
Compared with the prior art, the embodiment of the invention has at least one of the following beneficial effects:
acquiring sampling signals of the same signal source through multiple paths of antennas, performing first preprocessing on sub-channels, and determining center frequencies of the sub-channels so as to flexibly configure sub-channel frequency points according to actual requirements, and the method is suitable for application scenes in which prior information is missing on the received signal frequency points under non-cooperative satellite communication; and carrying out second preprocessing on the sampling signals according to the center frequency of the sub-channels to determine baseband signals, and carrying out first-stage channelization and second-stage channelization on the baseband signals to separate the sub-channels and narrowband signals in the satellite signals, so as to realize the matching of various target output sampling rates of the narrowband signals, and then processing the narrowband signals into target synthesized signals, thereby improving the signal synthesis performance and the signal-to-noise ratio gain capability.
Drawings
Other features, objects and advantages of the present invention will become more apparent upon reading of the detailed description of non-limiting embodiments, given with reference to the accompanying drawings in which:
fig. 1 is a flow chart illustrating a GPU-based multi-antenna array satellite signal synthesis method, according to an example embodiment.
Fig. 2 is a flow chart illustrating a method of determining a baseband signal according to an exemplary embodiment.
Fig. 3 is a flowchart illustrating a GPU-based multi-antenna array satellite signal synthesis method, according to another example embodiment.
Fig. 4 is a flowchart illustrating a method of determining a target composite signal, according to an example embodiment.
Fig. 5 is a flow chart illustrating a method of determining a target composite signal according to another exemplary embodiment.
FIG. 6 is a flowchart illustrating a GPU-based multi-antenna array satellite signal synthesis method executing on a CPU-GPU heterogeneous platform, according to an example embodiment.
FIG. 7 is a schematic diagram illustrating a GPU-based multi-antenna array satellite signal synthesizer according to an exemplary embodiment
Fig. 8 is a block diagram of an electronic device, according to an example embodiment.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the present invention, but are not intended to limit the invention in any way. It should be noted that variations and modifications could be made by those skilled in the art without departing from the inventive concept. These are all within the scope of the present invention.
It should be noted that, all actions of acquiring signals, information or data in the present invention are performed under the condition of conforming to the corresponding data protection rule policy of the country of the location and obtaining the authorization given by the owner of the corresponding device.
Fig. 1 is a flow chart illustrating a method of GPU-based multi-antenna array satellite signal synthesis, according to an example embodiment. As shown in fig. 1, the method may include the steps of:
s11, sampling signals of the same signal source are obtained through multiple paths of antennas.
Wherein the sampling signal is a broadband intermediate frequency signal, the signal type is a real signal, and the sampling rate of the sampling signal is f s 。
In some possible embodiments, after the sampling signals of the same signal source are acquired through multiple antennas, the sampling signals are packaged according to a preset duration, and a packaged signal is determined.
The preset duration may be a fixed duration, the packet signal is srl, and the data length of the packet signal is ILEN.
S12, carrying out first preprocessing on the sub-channels, and determining the center frequency of the sub-channels.
By performing the first preprocessing on the sub-channels, where the first preprocessing may use a dual sliding window energy detection, correlation detection, and wavelet transform method, those skilled in the art should understand that the present invention may also use other methods to perform the first preprocessing on the sub-channels, which all fall within the protection scope of the present invention.
Wherein the center frequency of the determined sub-channel is f c ＝{f 1 ,f 2 ,...,f n }。
S13, carrying out second preprocessing on the sampling signals according to the center frequency of the sub-channels to determine baseband signals.
The target sampling rates of the sub-channels at different frequency points in the sampling signal are different due to different modulation modes, and the sampling signal needs to be subjected to second preprocessing to determine multiple output sampling rates.
In the invention, the second preprocessing comprises down-conversion processing, P times interpolation processing and Q times extraction processing, the sampling signal is subjected to down-conversion processing by taking the center frequency of a sub-channel as a target, so that the sub-channels positioned at different frequency points are converted to a baseband, wherein the sub-channel frequency points can be flexibly configured into any center frequency point determined through preprocessing, meanwhile, the P times interpolation processing and the Q times extraction processing are adopted, and finally, the low-pass filtering processing is adopted to inhibit spectrum aliasing, so that the preliminary adjustment of various output sampling rates of the sampling signal is realized.
S14, the baseband signal is subjected to primary channelization and secondary channelization, and a target synthesized signal is determined.
Wherein the target composite signal comprises narrowband signals of a plurality of output sampling rates.
The baseband signal is first extracted D times and then first-stage channelized and second-stage channelized successively. The baseband signal is channelized at a first stage to output a first-stage narrowband signal, and the first-stage narrowband signal is channelized at a second stage to output a second-stage narrowband signal.
In a possible embodiment, the sub-channels at different frequency points of the sampled signal adopt different modulation modes, so that the requirements of the sub-channels on output sampling rates of the narrowband signals after the channelizing are different, the sampling rate required to be output finally is set to be m kinds, the preliminary processing of any multiplying power is realized through the first-stage channelizing, the phase difference is eliminated through the second-stage channelizing before the structuring (DTF) operation, and then the target synthesized signal is determined, wherein the target synthesized signal comprises m kinds of narrowband signals with output sampling rates.
Acquiring sampling signals of the same signal source through multiple paths of antennas, performing first preprocessing on sub-channels, and determining center frequencies of the sub-channels so as to flexibly configure sub-channel frequency points according to actual requirements, and the method is suitable for application scenes in which prior information is missing on the received signal frequency points under non-cooperative satellite communication; and carrying out second preprocessing on the sampling signals according to the center frequency of the sub-channels to determine baseband signals, and carrying out first-stage channelization and second-stage channelization on the baseband signals to separate the sub-channels and narrowband signals in the satellite signals, so as to realize the matching of various target output sampling rates of the narrowband signals, and then synthesizing the narrowband signals into target synthesized signals, thereby improving the signal synthesis performance and the signal-to-noise ratio gain capability.
Fig. 2 is a flow chart illustrating a method of determining a baseband signal according to an exemplary embodiment.
In a possible embodiment, the second preprocessing is performed on the sampling signal according to the center frequency of the sub-channel to determine the baseband signal, as shown in fig. 2, including S21 to S23.
S21, down-conversion processing is carried out on the sampling signals with the center frequency of the sub-channel as the target.
The frequency conversion method used in the down conversion process is as follows:
wherein f c Is the carrier frequency, T s For sampling period, n is the number of samples
S22, carrying out P times interpolation processing and Q times extraction processing on the sampling signal after the down-conversion processing.
Let the original sampling rate be f s Is adjusted to f s *P/Q
S23, carrying out low-pass filtering processing on the sampling signals subjected to P times interpolation processing and Q times extraction processing, filtering spectrum aliasing caused by interpolation and extraction, and determining a baseband signal.
The n baseband signals are sd= { sd1, sd2, & gt, sd }, and the input sampling rate is primarily adjusted through down-conversion processing, P-time interpolation processing and Q-time extraction processing in the preprocessing process.
Fig. 3 is a flowchart illustrating a GPU-based multi-antenna array satellite signal synthesis method, according to another example embodiment.
In some possible embodiments, as shown in fig. 3, the method further comprises S15.
After S13, S15 is included:
s15, integer time delay compensation processing is carried out on the baseband signal.
As an example, in the integer delay compensation process, the integer delay uses a preset length of the baseband signal sd1, so that the calculation amount can be reduced and the calculation efficiency can be improved without affecting the result of the delay. And determining the number of offset integer sample points by adopting a cross-correlation method on the baseband signal so as to carry out sample point offset on the baseband signal, and compensating the baseband signal according to the residual time delay and the phase difference value fed back by the rear end.
Fig. 4 is a flowchart illustrating a method of determining a target composite signal, according to an example embodiment.
In some possible embodiments, the baseband signal is first-stage channelized and second-stage channelized to determine a target composite signal, as shown in fig. 4, including S31 through S33.
S31, the baseband signal is subjected to linear shift processing through the first-stage channelization, and the signal subjected to the linear shift processing is determined.
S32, performing cyclic shift processing on the baseband signal subjected to linear shift processing through second-stage channelization, and determining the signal subjected to cyclic shift processing.
S33, performing discrete Fourier transform processing on the signal subjected to the cyclic shift processing to determine a target synthesized signal.
According to the technical scheme, the baseband signal is subjected to linear shift processing through the first-stage channelization, so that the linear shift processing based on the programmable memory is introduced, two problems of continuous access or repeated access are avoided in the linear shift processing, and when the GPU-based multi-antenna array satellite signal synthesis method is executed on a CPU-GPU heterogeneous platform, the frequency of global access to the memory can be reduced based on the programmable memory, the data hit rate is improved, and the access speed and the data processing performance are improved.
As an example, the embodiments of S31 to S33 described above perform matching of multiple output sampling rates by linear shift processing, polyphase filter processing, cyclic shift processing, and discrete fourier transform in order, including:
setting the sampling signal as a D path, outputting z of a kth path signal in the multiphase filtering k (m) is:
wherein y (n) is an input signal, h LP (mD-i) represents the time window of the filter, and mD represents the delay of h (mD-i).
Replacing the time window sliding according to the translation signal, and making r=i-mD, then
In one possible embodiment, the window function is relatively fixed due to the change in output caused by the change in m, and the sampling signal is taken as the input signal at the engineering implementation angle, and the sampling rate is f s It is loaded into D channels, m varies over time, and it corresponds to K output values per transform.
Setting h LP And (N) the length of the filter is N, wherein if K is less than or equal to N, in order to facilitate discrete Fourier transform, y (r) is folded once every K points, and the calculated results of each group after folding are accumulated, and the result of the accumulated operation is used as an output result, so that the number of channels K and the signal extraction multiple D are decoupled, and the channelizing is no longer satisfied with the constraint condition of K=D.
In one possible embodiment, the input register and the filter register are set to be N in size, the input data is loaded with D points as step sizes, multiplied by corresponding filter coefficients and updated, and folded with K points as a group, N/K groups of data are shared, and the weighted sum of K points is determined by calculation from the N/K groups of data.
And performing discrete Fourier transform (FFT) processing on the K points, and determining one point output updated by each path at each moment. As an example, if the preset number of output channels is K to realize K/D times of sampling rate adjustment, the input sampling rate f s Adjusted to f s Output sampling rate of/D.
As an example, the original output sampling rate is 100kHz, and the sampling rate is now adjusted to 200kHz, where let k=32, d=16, the adjustment of the output sampling rate can be achieved, and no additional sampling operation is required after the channelization, so as to save the calculation overhead.
And finally synthesizing the narrowband signal subjected to the second-stage channelization into a target synthesized signal according to the determined output sampling rate.
Fig. 5 is a flowchart illustrating a method of determining a target composite signal according to another exemplary embodiment.
In some possible embodiments, as shown in fig. 5, the method further comprises:
s34 to S37 are also included after the steps of first-stage channelization and second-stage channelization of the baseband signal, i.e., after S32.
S34, carrying out phase difference processing on the narrowband signals at the same frequency point of the multipath antenna to determine the phase difference.
In one possible embodiment, the signal determined by the baseband signal undergoing the first-stage channelization and the second-stage channelization is a narrowband signal, the narrowband signals at the same frequency point of different antennas are subjected to phase difference processing, and the phase difference is denoted as phi.
Let the first sub-band signal be s 1 (t)＝Asin[ω c t+θ 1 (t)]The second sub-band signal is s 2 (t)＝Asin[ω c t+θ 2 (t)]。
Wherein A represents the amplitude of the subband signal, ω c Represents the angular frequency, θ, of the subband signal 1 (t) represents the phase of the first sub-band signal, θ 2 (t) represents the phase of the second subband signal.
Conjugate multiplication processing is carried out on each path of subband signal in the multipath antenna, then average processing is carried out, and average signal is determinedAnd then carrying out arctangent processing on the average signal to determine the phase difference delta theta of each path of subband signal.
And S35, fitting the phase difference based on a least square method to determine the residual time delay and the phase deviation value.
And fitting the phase differences at the same frequency point in the multipath antenna based on a least square method to determine a straight line, wherein the slope of the straight line represents residual time delay which is residual decimal time delay, and the intercept represents a phase deviation value, namely an initial phase difference.
S36, performing decimal delay compensation processing on broadband signals of the multipath antennas according to the residual delay.
S37, carrying out phase compensation processing on the narrowband signals at the same frequency point of the multipath antenna according to the phase offset value.
After compensation, the signal-to-noise ratio gain of the narrowband signal reaches 2.5dB under the condition of double-antenna receiving, and is more than 4.5dB under the condition of four-antenna receiving.
FIG. 6 is a flowchart illustrating a GPU-based multi-antenna array satellite signal synthesis method executing on a CPU-GPU heterogeneous platform, according to an example embodiment.
In some possible embodiments, as shown in fig. 6, the method further comprises:
and transplanting the method to a CPU-GPU heterogeneous platform for execution through CUDA programming.
Wherein, in the CPU:
the first step: starting the treatment;
and a second step of: data transfer to host device (data transfer hosttopdevice);
and a third step of: center frequency of the pre-processing sub-channel (pre-processing positioning frequency point f c )；
Fourth step: waiting for GPU-end kernel execution (executing kernel);
fifth step: data transfer device to host (data transfer devicehost);
sixth step: and separating output according to the channels.
In the GPU:
the first step: kernel down conversion processing (down conversion kernel);
and a second step of: kernel channelization (channelizing);
and a third step of: kernel full spectrum synthesis (full spectrum synthesis kernel).
After the third step is executed in the CPU, the GPU starts to execute the first step, when the third step of the GPU is executed, the CPU executes the fifth step, the data transmission equipment is connected to the host, and finally the CPU executes the sixth step, and the data is separated and output according to the channel.
According to the technical scheme, the CPU-GPU heterogeneous platform is adopted to execute the multi-antenna array satellite signal synthesis method based on the GPU, so that the calculated amount can be effectively reduced, the data hit rate is improved, and faster access speed and better processing performance are provided.
In some possible embodiments, the method further comprises:
when the design method is executed on the CPU-GPU heterogeneous platform, parallel processing is carried out in a mode of asynchronous copying, page locking memory optimization, on-chip programmable memory optimization and thread organization structure.
Wherein, for lock page memory and asynchronous copy:
according to the length of the processed data, the memory of the host is opened up into the page-locking memory corresponding to the length of the processed data, and the data transmission method with asynchronous non-blocking property is adopted, and the time for returning the data transmission control right to the host side is utilized to execute the kernel function in parallel, so that the time for data transmission is hidden.
According to the technical scheme, the page-locked host memory is mapped to the address space of the device, so that the time for copying the page-locked host memory to the device memory or copying the page-locked host memory from the device memory is omitted, the bandwidth between the host and the device is improved, the copying time is shortened, namely, the cost of addressing time from a virtual address to a physical address and the cost of memory replacement time from a disk space to the physical space are saved, the data transmission time is hidden by using an asynchronous copying method, and the processing speed is improved.
For programmable memory and thread organization, the highly multiplexed data is stored temporarily in the memory on the GPU chip by software defining the programmable memory, in the present invention, the thread block dimension value is set as (128,4).
In the technical scheme, the access delay of the on-chip memory is lower than that of the off-chip memory, so that the cache hit rate is improved, and the access and reading time of data is shortened. The hardware self architecture determines that the size of the usable on-chip memory and the number of registers proprietary to the computing thread are fixed, the higher the utilization rate of the on-chip memory and the thread registers is, the faster the data access speed is, but the parallel thread block data on each stream processor can be reduced at the same time, the constraint factors are comprehensively considered, and the dimension value of the thread block is set, so that the GPU-CPU heterogeneous platform has better processing performance.
In conclusion, by adopting asynchronous copying, page-locked memory optimization, on-chip programmable memory optimization and thread organization structure to perform parallel optimization processing, the computational complexity can be effectively reduced, so that the signal synthesis method based on software implementation has high flexibility and high throughput and real-time processing capability.
Based on the same inventive concept, the invention also provides a multi-antenna array satellite signal synthesis device based on the GPU. Fig. 7 is a schematic diagram illustrating a GPU-based multi-antenna array satellite signal synthesis apparatus, according to an example embodiment. Referring to fig. 7, the GPU-based multi-antenna array satellite signal synthesizing apparatus 100 includes a sampling module 110, a first preprocessing module 120, a second preprocessing module 130, and a determining module 140:
the sampling module 110 is configured to obtain sampling signals of the same signal source through multiple antennas;
a first preprocessing module 120, configured to perform a first preprocessing on a subchannel, and determine a center frequency of the subchannel;
a second preprocessing module 130, configured to perform a second preprocessing on the sampling signal according to the center frequency of the subchannel, and determine a baseband signal;
a determining module 140, configured to perform a first-stage channelization and a second-stage channelization on the baseband signal, and determine a target composite signal, where the target composite signal includes narrowband signals with multiple output sampling rates.
According to the technical scheme, sampling signals of the same signal source are obtained through the multipath antennas, first preprocessing is carried out on the sub-channels, and the center frequencies of the sub-channels are determined, so that sub-channel frequency points are flexibly configured according to actual requirements, and the method is suitable for application scenes in which prior information is missing on the received signal frequency points under non-cooperative satellite communication; and carrying out second preprocessing on the sampling signals according to the center frequency of the sub-channels to determine baseband signals, and carrying out first-stage channelization and second-stage channelization on the baseband signals to separate the sub-channels and narrowband signals in the satellite signals, so as to realize the matching of various target output sampling rates of the narrowband signals, and then processing the narrowband signals into target synthesized signals, thereby improving the signal synthesis performance and the signal-to-noise ratio gain capability.
Optionally, the determining module 140 includes:
the first determining submodule is used for carrying out linear shift processing on the baseband signal through first-stage channelization and determining a signal after the linear shift processing;
a second determining submodule, configured to perform cyclic shift processing on the baseband signal after the linear shift processing through second-stage channelization, and determine a signal after the cyclic shift processing;
and the third determining submodule is used for performing discrete Fourier transform processing on the signal subjected to the cyclic shift processing and determining a target synthesized signal.
Optionally, the apparatus 100 further comprises:
the phase difference processing module is used for carrying out phase difference processing on the narrowband signals at the same frequency point of the multipath antenna after the step of carrying out first-stage channelization and second-stage channelization on the baseband signals, so as to determine a phase difference;
the fitting module is used for carrying out fitting processing on the phase difference based on a least square method and determining residual time delay and phase deviation values;
the decimal time delay compensation module is used for carrying out first decimal time delay compensation processing on the narrowband signals at the same frequency point of the multipath antenna according to the residual time delay;
and the phase compensation module is used for carrying out phase compensation processing on the narrowband signals at the same frequency point of the multipath antenna according to the phase offset value.
Optionally, the second preprocessing module 130 includes:
the first preprocessing sub-module is used for carrying out down-conversion processing on the sampling signal by taking the center frequency of the sub-channel as a target;
the second preprocessing submodule carries out P times interpolation processing and Q times extraction processing on the sampling signals after the down-conversion processing;
and the third preprocessing sub-module is used for carrying out low-pass filtering processing on the sampling signals subjected to the P times interpolation processing and the Q times extraction processing to determine baseband signals.
Optionally, the apparatus 100 further comprises:
and the integer time delay compensation module is used for carrying out second preprocessing on the sampling signals according to the center frequency of the sub-channels, and carrying out integer time delay compensation processing on the baseband signals after the step of determining the baseband signals.
Optionally, the apparatus 100 further comprises:
and the execution module is used for transplanting the design method to a CPU-GPU heterogeneous platform for execution through CUDA programming.
Optionally, the apparatus 100 further comprises:
and the parallel processing module is used for carrying out parallel processing in the modes of asynchronous copying, page locking memory optimization, on-chip programmable memory optimization and thread organization structure when the design method is executed on the CPU-GPU heterogeneous platform.
The specific manner in which the various modules perform the operations in the apparatus of the above embodiments have been described in detail in connection with the embodiments of the method, and will not be described in detail herein.
Fig. 8 is a block diagram of an electronic device, according to an example embodiment. As shown in fig. 8, the electronic device 800 may include: a processor 801, a memory 802. The electronic device 800 may also include one or more of a multimedia component 803, an input/output interface 804, and a communication component 805.
The processor 801 is configured to control overall operation of the electronic device 800 to perform all or part of the steps in the GPU-based multi-antenna array satellite signal synthesis method according to the first aspect. The memory 802 is used to store various types of data to support operation at the electronic device 800, which may include, for example, instructions for any application or method operating on the electronic device 800, as well as application-related data, such as contact data, messages sent and received, pictures, audio, video, and so forth. The Memory 802 may be implemented by any type or combination of volatile or non-volatile Memory devices, such as static random access Memory (Static Random Access Memory, SRAM for short), electrically erasable programmable Read-Only Memory (Electrically Erasable Programmable Read-Only Memory, EEPROM for short), erasable programmable Read-Only Memory (Erasable Programmable Read-Only Memory, EPROM for short), programmable Read-Only Memory (Programmable Read-Only Memory, PROM for short), read-Only Memory (ROM for short), magnetic Memory, flash Memory, magnetic disk, or optical disk. The multimedia component 803 may include a screen and an audio component. Wherein the screen may be, for example, a touch screen, the audio component being for outputting and/or inputting audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signals may be further stored in the memory 802 or transmitted through the communication component 805. The audio assembly further comprises at least one speaker for outputting audio signals. The input/output interface 804 provides an interface between the processor 801 and other interface modules, which may be a keyboard, mouse, buttons, etc. These buttons may be virtual buttons or physical buttons. The communication component 805 is used for wired or wireless communication between the electronic device 800 and other devices. Wireless communication, such as Wi-Fi, bluetooth, near field communication (Near Field Communication, NFC for short), 2G, 3G, 4G, NB-IOT, eMTC, or other 5G, etc., or one or a combination of more of them, is not limited herein. The corresponding communication component 305 may thus comprise: wi-Fi module, bluetooth module, NFC module, etc.
In another exemplary embodiment, there is also provided a non-transitory computer readable storage medium including program instructions which, when executed by a processor, implement the steps of the GPU-based multi-antenna array satellite signal synthesis method of the first aspect described above. For example, the computer readable storage medium may be the memory including program instructions described above, which are executable by a processor of an electronic device to perform the GPU-based multi-antenna array satellite signal synthesis method described above.
In another exemplary embodiment, a computer program product is also provided, the computer program product comprising a computer program executable by a programmable apparatus, the computer program having code portions for performing the GPU-based multi-antenna array satellite signal synthesis method described above when executed by the programmable apparatus.
The foregoing describes specific embodiments of the present invention. It is to be understood that the invention is not limited to the particular embodiments described above, and that various changes and modifications may be made by one skilled in the art within the scope of the claims without affecting the spirit of the invention. The above-described preferred features may be used in any combination without collision.
Claims (10)
1. The method for synthesizing the multi-antenna array satellite signals based on the GPU is characterized by comprising the following steps of:
acquiring sampling signals of the same signal source through multiple paths of antennas;
performing first preprocessing on a sub-channel, and determining the center frequency of the sub-channel;
performing second preprocessing on the sampling signals according to the center frequency of the sub-channel to determine baseband signals;
and carrying out first-stage channelization and second-stage channelization on the baseband signal to determine a target synthesized signal, wherein the target synthesized signal comprises narrowband signals with various output sampling rates.
2. The method of claim 1, wherein said subjecting the baseband signal to a first stage of channelization and a second stage of channelization to determine a target composite signal comprises:
performing linear shift processing on the baseband signal through first-stage channelization, and determining a signal after the linear shift processing;
performing cyclic shift processing on the baseband signal subjected to the linear shift processing through second-stage channelization, and determining a signal subjected to the cyclic shift processing;
and performing discrete Fourier transform processing on the signal subjected to the cyclic shift processing to determine a target synthesized signal.
3. The method according to claim 1, wherein the method further comprises:
after the step of carrying out first-stage channelization and second-stage channelization on the baseband signals, carrying out phase difference processing on narrowband signals at the same frequency point of the multipath antenna to determine a phase difference;
fitting the phase difference based on a least square method to determine a residual time delay and a phase deviation value;
performing decimal delay compensation processing on broadband signals of the multipath antennas according to the residual delay;
and carrying out phase compensation processing on the narrowband signals at the same frequency point of the multipath antenna according to the phase offset value.
4. The method of claim 1, wherein said second preprocessing the sampled signal based on the center frequency of the sub-channel to determine a baseband signal comprises:
performing down-conversion processing on the sampling signal by taking the center frequency of the sub-channel as a target;
performing P times interpolation processing and Q times extraction processing on the sampling signal subjected to the down-conversion processing;
and performing low-pass filtering processing on the sampling signals subjected to the P times interpolation processing and the Q times extraction processing to determine baseband signals.
5. The method according to claim 1, wherein the method further comprises:
and after the step of determining the baseband signal, carrying out integer delay compensation processing on the baseband signal.
6. The method according to claim 1, wherein the method further comprises:
and transplanting the design method to a CPU-GPU heterogeneous platform for execution through CUDA programming.
7. The method of claim 6, wherein the method further comprises:
when the design method is executed on the CPU-GPU heterogeneous platform, parallel processing is carried out in a mode of asynchronous copying, page locking memory optimization, on-chip programmable memory optimization and thread organization structure.
8. The utility model provides a many antennas group array satellite signal synthesizer based on GPU which characterized in that includes:
the sampling module is used for acquiring sampling signals of the same signal source through multiple paths of antennas;
the first preprocessing module is used for carrying out first preprocessing on the sub-channels and determining the center frequency of the sub-channels;
the second preprocessing module is used for carrying out second preprocessing on the sampling signals according to the center frequency of the sub-channel to determine baseband signals;
and the determining module is used for carrying out first-stage channelization and second-stage channelization on the baseband signal to determine a target synthesized signal, wherein the target synthesized signal comprises narrowband signals with various output sampling rates.
9. A non-transitory computer readable storage medium having stored thereon a computer program, characterized in that the program when executed by a processor realizes the steps of the method according to any of claims 1-7.
10. An electronic device, comprising:
a memory having a computer program stored thereon;
a processor for executing the computer program in the memory to implement the steps of the method of any one of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310548816.7A CN116566471A (en) | 2023-05-16 | 2023-05-16 | GPU-based multi-antenna array satellite signal synthesis method, device and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310548816.7A CN116566471A (en) | 2023-05-16 | 2023-05-16 | GPU-based multi-antenna array satellite signal synthesis method, device and medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116566471A true CN116566471A (en) | 2023-08-08 |
Family
ID=87487612
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310548816.7A Pending CN116566471A (en) | 2023-05-16 | 2023-05-16 | GPU-based multi-antenna array satellite signal synthesis method, device and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116566471A (en) |
-
2023
- 2023-05-16 CN CN202310548816.7A patent/CN116566471A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102480302B (en) | Method and system for reducing interference caused by clock signal in receiver | |
| US8125213B2 (en) | System for extraction and analysis of significant radioelectric signals | |
| US20020027958A1 (en) | Feedback channel signal recovery | |
| US10039020B2 (en) | Dominant signal detection method and apparatus | |
| CN102318198A (en) | The efficient sampling and the reconstruction of sparse multiband signal | |
| WO2007110750A2 (en) | Low complexity subband-domain filtering in the case of cascaded filter banks | |
| JP2009300284A (en) | Arrival electric wave direction measuring device, arrival electric wave direction measuring method, and arrival electric wave direction measuring program | |
| CN110210101A (en) | A kind of non-homogeneous narrow transition band filter group of dynamic based on CEM FRM and design method | |
| US20190036675A1 (en) | Multifunction channelizer/ddc architecture for a digital receiver/exciter | |
| US20050128937A1 (en) | Determination of correlation in the frequency domain | |
| Im et al. | Implementation of SDR-based digital IF channelizer/de-channelizer for multiple CDMA signals | |
| US20210091752A1 (en) | Resource conserving weighted overlap-add channelizer | |
| CN116566471A (en) | GPU-based multi-antenna array satellite signal synthesis method, device and medium | |
| WO2021227217A1 (en) | Signal-to-noise ratio estimation method and apparatus, and electronic device and storage medium | |
| CN111464228B (en) | Multichannel VDES receiving and dispatching preprocessing system based on satellite-borne DBF | |
| CN109687826B (en) | Multi-channel down-conversion method and device based on model design and storage medium | |
| CN108900210B (en) | Method and system for eliminating parallel interference of multimode measurement and control signals | |
| CN111934728B (en) | Digital multi-beam antenna processing method, device and equipment | |
| CN113252957A (en) | Digital down-conversion implementation method and digital oscilloscope | |
| CN113067150A (en) | Anti-interference antenna and control method for anti-interference antenna | |
| US7142607B2 (en) | Efficient subband channelizer and subband combiner with oversampled subbands | |
| JP3146093B2 (en) | Two-stage fast Fourier transform method | |
| Azza et al. | Implementation of Cognitive Radio applications on a Software Defined Radio plateform | |
| US20230379205A1 (en) | System, method, apparatus, and computer program product for providing adaptive wireless noise cancellation | |
| CN113126033B (en) | Radar signal digital quadrature down-conversion method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |