CN111934728B - Digital multi-beam antenna processing method, device and equipment - Google Patents

Digital multi-beam antenna processing method, device and equipment Download PDF

Info

Publication number
CN111934728B
CN111934728B CN202010549747.8A CN202010549747A CN111934728B CN 111934728 B CN111934728 B CN 111934728B CN 202010549747 A CN202010549747 A CN 202010549747A CN 111934728 B CN111934728 B CN 111934728B
Authority
CN
China
Prior art keywords
frequency
incoming wave
satellite
space
signal
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.)
Active
Application number
CN202010549747.8A
Other languages
Chinese (zh)
Other versions
CN111934728A (en
Inventor
任超
屠恩源
吴红甲
郑永翔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BEIJING BDSTAR NAVIGATION Co.,Ltd.
Original Assignee
BEIJING BDSTAR NAVIGATION CO LTD
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by BEIJING BDSTAR NAVIGATION CO LTD filed Critical BEIJING BDSTAR NAVIGATION CO LTD
Priority to CN202010549747.8A priority Critical patent/CN111934728B/en
Publication of CN111934728A publication Critical patent/CN111934728A/en
Application granted granted Critical
Publication of CN111934728B publication Critical patent/CN111934728B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H04B7/0426—Power distribution
    • H04B7/043—Power distribution using best eigenmode, e.g. beam forming or beam steering
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/0003—Software-defined radio [SDR] systems, i.e. systems wherein components typically implemented in hardware, e.g. filters or modulators/demodulators, are implented using software, e.g. by involving an AD or DA conversion stage such that at least part of the signal processing is performed in the digital domain
    • H04B1/0007—Software-defined radio [SDR] systems, i.e. systems wherein components typically implemented in hardware, e.g. filters or modulators/demodulators, are implented using software, e.g. by involving an AD or DA conversion stage such that at least part of the signal processing is performed in the digital domain wherein the AD/DA conversion occurs at radiofrequency or intermediate frequency stage
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06—Receivers
    • H04B1/16—Circuits
    • 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

Abstract

A digital multi-beam antenna processing method, apparatus, device, and computer-readable storage medium, wherein the method comprises: converting a received analog intermediate frequency signal into a digital intermediate frequency signal, and converting the digital intermediate frequency signal into frequency domain data through Fourier transformation; determining a space frequency incoming wave direction constraint matrix based on an incoming wave angle set of a satellite; determining a filtering weight vector of each wave beam according to the space-frequency incoming wave direction constraint matrix and the frequency domain covariance matrix; performing space-frequency filtering on the frequency domain data by using the filtering weight vector, and performing inverse Fourier transform on the filtered frequency domain signal to obtain a time domain signal corresponding to each satellite; and carrying out correlation peak detection on the time domain signal, and updating an incoming wave angle set of the satellite based on a detection result. According to the method and the device, the incoming wave angle set of the satellite is updated based on the correlation peak detection result, the space-time filtering under the condition without the expected signal prior is realized, and the robustness of the space-time filtering is enhanced.

Description

Digital multi-beam antenna processing method, device and equipment
Technical Field
The present disclosure relates to the field of communications, and more particularly, to a digital multi-beam antenna processing method, apparatus, device, and computer-readable storage medium.
Background
In order to suppress forced-press interference in a complex electromagnetic environment, the array antenna beam forming technology is widely applied to a spread spectrum communication system.
The adaptive beamforming technology often requires a priori knowledge of an expected signal incoming direction and the like, but the priori knowledge of the expected signal, an array channel error and the like cannot be obtained in an array application environment, so that the broadband beamformer fails.
Disclosure of Invention
A digital multi-beam antenna processing method, apparatus, device and computer readable storage medium are provided to enhance robustness of wideband beamformer space-time filtering.
The embodiment of the application provides a digital multi-beam antenna processing method, which comprises the following steps:
converting a received analog intermediate frequency signal into a digital intermediate frequency signal, and converting the digital intermediate frequency signal into frequency domain data through Fourier transformation;
determining a space frequency incoming wave direction constraint matrix based on an incoming wave angle set of a satellite, wherein the incoming wave angle set is preset with an initial value;
determining a filtering weight vector of each wave beam according to the space-frequency incoming wave direction constraint matrix and the frequency domain covariance matrix;
performing space-frequency filtering on the frequency domain data by using the filtering weight vector, and performing inverse Fourier transform on the filtered frequency domain signal to obtain a time domain signal corresponding to each satellite;
and carrying out correlation peak detection on the time domain signal, and updating an incoming wave angle set of the satellite based on a detection result.
An embodiment of the present application further provides a digital multi-beam antenna processing apparatus, including:
the signal conversion module is used for converting the received analog intermediate frequency signal into a digital intermediate frequency signal and converting the digital intermediate frequency signal into frequency domain data through Fourier transformation;
the constraint matrix module is used for determining a space frequency incoming wave direction constraint matrix based on an incoming wave angle set of a satellite;
the weight vector module is used for determining a filtering weight vector of each wave beam according to the space frequency incoming wave direction constraint matrix and the frequency domain covariance matrix;
the filtering output module is used for performing space-frequency filtering on the frequency domain data by using the filtering weight vector and performing inverse Fourier transform on the filtered frequency domain signal to obtain a time domain signal corresponding to each satellite;
and the updating module is used for carrying out correlation peak detection on the time domain signal and updating the incoming wave angle set of the satellite based on the detection result.
An embodiment of the present application further provides a digital multi-beam antenna processing apparatus, including: a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the digital multi-beam antenna processing method when executing the program.
Embodiments of the present application also provide a computer-readable storage medium storing computer-executable instructions for performing the digital multi-beam antenna processing method.
Compared with the related art, the embodiment of the application comprises the following steps: converting a received analog intermediate frequency signal into a digital intermediate frequency signal, and converting the digital intermediate frequency signal into frequency domain data through Fourier transformation; determining a space frequency incoming wave direction constraint matrix based on an incoming wave angle set of a satellite, wherein the incoming wave angle set is preset with an initial value; determining a filtering weight vector of each wave beam according to the space-frequency incoming wave direction constraint matrix and the frequency domain covariance matrix; performing space-frequency filtering on the frequency domain data by using the filtering weight vector, and performing inverse Fourier transform on the filtered frequency domain signal to obtain a time domain signal corresponding to each satellite; and carrying out correlation peak detection on the time domain signal, and updating an incoming wave angle set of the satellite based on a detection result. According to the method and the device, the incoming wave angle set of the satellite is updated based on the correlation peak detection result, the space-time filtering under the condition without the expected signal prior is realized, and the robustness of the space-time filtering is enhanced.
Additional features and advantages of the application will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by the practice of the application. Other advantages of the present application may be realized and attained by the instrumentalities and combinations particularly pointed out in the specification and the drawings.
Drawings
The accompanying drawings are included to provide an understanding of the present disclosure and are incorporated in and constitute a part of this specification, illustrate embodiments of the disclosure and together with the examples serve to explain the principles of the disclosure and not to limit the disclosure.
FIG. 1 is a schematic diagram of a wideband beamformer architecture;
fig. 2 is a flowchart of a digital multi-beam antenna processing method according to an embodiment of the present application;
fig. 3 is a schematic diagram of digital multi-beam antenna processing according to an exemplary application of the present application;
fig. 4 is a schematic diagram of a digital multi-beam antenna processing apparatus according to an embodiment of the present application.
Detailed Description
The present application describes embodiments, but the description is illustrative rather than limiting and it will be apparent to those of ordinary skill in the art that many more embodiments and implementations are possible within the scope of the embodiments described herein. Although many possible combinations of features are shown in the drawings and discussed in the detailed description, many other combinations of the disclosed features are possible. Any feature or element of any embodiment may be used in combination with or instead of any other feature or element in any other embodiment, unless expressly limited otherwise.
The present application includes and contemplates combinations of features and elements known to those of ordinary skill in the art. The embodiments, features and elements disclosed in this application may also be combined with any conventional features or elements to form a unique inventive concept as defined by the claims. Any feature or element of any embodiment may also be combined with features or elements from other inventive aspects to form yet another unique inventive aspect, as defined by the claims. Thus, it should be understood that any of the features shown and/or discussed in this application may be implemented alone or in any suitable combination. Accordingly, the embodiments are not limited except as by the appended claims and their equivalents. Furthermore, various modifications and changes may be made within the scope of the appended claims.
Further, in describing representative embodiments, the specification may have presented the method and/or process as a particular sequence of steps. However, to the extent that the method or process does not rely on the particular order of steps set forth herein, the method or process should not be limited to the particular sequence of steps described. Other orders of steps are possible as will be understood by those of ordinary skill in the art. Therefore, the particular order of the steps set forth in the specification should not be construed as limitations on the claims. Further, the claims directed to the method and/or process should not be limited to the performance of their steps in the order written, and one skilled in the art can readily appreciate that the sequences may be varied and still remain within the spirit and scope of the embodiments of the present application.
As shown in fig. 1, the wideband beamformer is a schematic diagram of a wideband beamformer, which may also be referred to as a space-frequency adaptive processor, a space-frequency two-dimensional processor, or a wideband receiver.
Broadband signal s taking into account Q far fieldsiQ is equal to θ 1,2iThe direction of the Q is incident on the M-element uniform line in space, i being 1,2, and the 1 st term s is generally assumed1Corresponds to the desired signal of interest, and siQ corresponds to Q-1 interfering signals. The array unit interval is half of the wavelength corresponding to the highest frequency of the received signal. Assuming that the bandwidth of the signal is B, the received data of the mth array element
Figure BDA0002542026620000041
Wherein tau ismiRepresenting the delay of the i-th signal received at the m-th element, nmAnd (t) is the thermal noise of the mth array element channel. One broadband signal can be divided into a plurality of frequency points through Fourier transform, and each frequency point is processed according to a narrowband beamforming method. Wherein, the wideband signal can be divided into K narrow sub-bands in frequency domain by K-point Discrete Fourier Transform (DFT), and each narrow sub-band is processed with narrow-band beam forming, and then the beam output is converted into time domain output, so that the frequency is increasedThe domain processing is a processing method for dividing sub-bands.
If the observation data of one observation time subsection is subjected to discrete Fourier transform of K points, the following broadband model can be obtained:
X(fk)＝A(fk)S(fk)+N(fk) k＝1,2,...,K
A(fk)＝[a(θ1,fk)a(θ2,fk)...a(θN,fk)]
in the formula, X (f)k)，S(fk)，N(fk) The discrete fourier transform of a received data vector, a signal vector and a noise vector corresponding to a certain frequency, respectively, K means that a signal with a bandwidth of B is divided into K subbands. Wherein the content of the first and second substances,
X(fk)＝[x1k,x2k,...,xMk]T k＝1,2,...,K
using M x 1 dimensional vectors w (f)k) The representation corresponds to the received data fkA processor weight vector of a subband, then
w(fk)＝[w1k,w2k,...,wMk]T k＝1,2,...,K
The filtered output signal of the wideband beamformer is
Figure BDA0002542026620000051
As shown in fig. 2, an embodiment of the present application provides a digital multi-beam antenna processing method, including:
step 101, converting a received analog intermediate frequency signal into a digital intermediate frequency signal, and converting the digital intermediate frequency signal into frequency domain data through fourier transform.
The fourier transform may be a Fast Fourier Transform (FFT) or a Discrete Fourier Transform (DFT).
In this step, a wideband beam former structure with M array elements may be adopted to perform AD (analog to digital) sampling on the analog intermediate frequency received by each array element to obtain a digital intermediate frequency signal x1(h),...,xM(h) And h represents the discrete point sequence number of the moment, the M paths of digital intermediate frequency signals are respectively subjected to K-point discrete Fourier transform, and each path of digital intermediate frequency signal is divided into K narrow sub-bands on the frequency domain. Obtaining a frequency domain signal xmk(l) The method indicates that the mth array element digital intermediate frequency data corresponds to a frequency point f after Fourier transformationkThe ith fast beat number of (c).
Step 102, determining a space frequency incoming wave direction constraint matrix based on an incoming wave angle set of a satellite, wherein the incoming wave angle set is preset with an initial value.
In one embodiment, the initial incoming wave direction angle of each satellite is 0 degrees, i.e., zenith normal.
The set of incoming wave directions of the N visible satellites may be set as: theta ═ theta1,θ2,...,θNAnd the initial incoming wave angle set is {0, 0.
In other embodiments, the initial incoming wave direction angle of each satellite may be set to a value other than 0.
In one embodiment, step 102 comprises:
and determining a space-frequency two-dimensional steering vector of a frequency point corresponding to each satellite signal according to the incoming wave angle set of the satellite, and constructing a space-frequency incoming wave direction constraint matrix based on the space-frequency two-dimensional steering vector.
The space-frequency two-dimensional steering vector of the nth satellite signal corresponding to the frequency point k is as follows:
Figure BDA0002542026620000061
n＝1,2,...N；k＝1,2,...K
wherein d is the array element spacing, M is the array element number, and c is the light velocity constant.
Constructing a space frequency incoming wave direction constraint matrix A (f) according to a space frequency two-dimensional steering vectork)＝[a(θ1,fk)a(θ2,fk)...a(θN,fk)]H。
And 103, determining a filtering weight vector of each wave beam according to the space frequency incoming wave direction constraint matrix and the frequency domain covariance matrix.
Corresponding to frequency point fkOf the frequency domain covariance matrix R (f) of the sampled datak) From the sample covariance matrix
Figure BDA0002542026620000063
Instead, it is expressed as:
Figure BDA0002542026620000062
where L represents the number of snapshots, x, of the frequency domain datamk(l) Showing that the m-th array element sampling data corresponds to a frequency point f after being subjected to Fourier transformkThe ith fast beat number of (c).
Filter weight vector w of nth beam(n)(fk)＝μR-1(fk)a(θn,fk)，μ＝[aH(θn,fk)R-1(fk)a(θn,fk)]-1Is a constant.
And 104, performing space-frequency filtering on the frequency domain data by using the filtering weight vector, and performing inverse Fourier transform on the filtered frequency domain signal to obtain a time domain signal corresponding to each satellite.
In this step, the weight vector w(n)(fk) Performing space-frequency filtering to obtain frequency domain signal [ y ] output by the nth wave beam filtering(n)(f1) y(n)(f2) ... y(n)(fK)]N is 1, 2. The frequency domain output of each wave beam is subjected to inverse Fourier transform (IFFT), and the time domain signal [ y ] of each satellite can be obtained(1)(t) y(2)(t) ... y(N)(t)]。
And 105, performing correlation peak detection on the time domain signal, and updating an incoming wave angle set of the satellite based on a detection result.
In one embodiment, step 105 comprises:
carrying out correlation peak detection on the time domain signal to obtain a detection result which is a correlation peak function of each satellite; and correcting the incoming wave angle of each satellite based on the correlation peak function to obtain an updated incoming wave angle set.
Wherein the signal [ y ] is converted into a signal(1)(t) y(2)(t) ... y(N)(t)]Respectively carrying out related peak detection to obtain related peak function F [ theta ] of each satelliten(h)]。
Wherein the correlation peak detection may comprise: output result y for each satellite(n)(t) performing finite length FFT to obtain Yn1(t) performing finite length FFT processing on the satellite pseudo code sequence PN1 to obtain Yn2(t) reacting Yn1(t).*conj(Yn2(t)) is subjected to IFFT, and correlation peak results are obtained. Wherein ". times.a. dot product, conj () represents a conjugate.
In one embodiment, the incoming wave angle θ of the nth satellite is corrected according to the following formulan(h+1)：
Figure BDA0002542026620000071
Wherein h represents the discrete point sequence number of the time, gamma is a learning factor and is an empirical value, and gamma belongs to (0, 1)]，F[θn(h)]Denotes thetan(h) Correlation peak function of direction.
The updated set of incoming wave angles is: [ theta ] of1(h+1),θ2(h+1),...,θN(h+1)]。
And calculating a space-frequency incoming wave direction constraint matrix by using the updated incoming wave angle set, repeating the space-frequency filtering processing, and iteratively converging to a maximum value of a related peak according to the rising gradient of a related peak-to-peak function, thereby obtaining the optimal processing of the digital wave beam under the parallel single-satellite constraint.
According to the method and the device, the incoming wave angle set of the satellite is updated based on the correlation peak detection result, the space-time filtering under the condition without the expected signal prior is realized, and the robustness of the space-time filtering is enhanced.
The following is a description of an application example. Referring to fig. 3, the receiving antenna array has M array elements, and M signals pass throughObtaining M analog intermediate frequency signals after the processing of the radio frequency circuit, and obtaining M paths of digital intermediate frequency signals x by the M analog intermediate frequency signals through AD sampling1(h),...,xM(h) In that respect And respectively carrying out K-point discrete Fourier transform on the M paths of digital intermediate frequency signals, and dividing each path of digital intermediate frequency signal into K narrow sub-bands on a frequency domain. Obtaining a frequency domain signal xmk(l) The corresponding frequency point f of the mth array element digital intermediate frequency data after FFT is shownkThe ith fast beat number of (c).
The incoming wave direction set of N visible satellites is set as follows: theta ═ theta1,θ2,...,θN}. The initial incoming wave direction angle is not set to 0 degree, i.e. the zenith normal direction. Then the initial set of angles is {0, 0.., 0 }.
Calculating a space-frequency two-dimensional steering vector of a frequency point k corresponding to the nth satellite signal according to the initial angle set,
Figure BDA0002542026620000081
n＝1,2,...N；k＝1,2,...K
wherein d is the array element spacing and c is the light velocity constant.
In a wideband beamformer, xmk(l) Frequency point f for m array element receivingkUpper frequency domain signal, let R (f)k) Is a frequency point fkUpper received signal covariance matrix, in practical applications, R (f)k) Usually by a sample covariance matrix
Figure BDA0002542026620000082
Instead, corresponding to frequency point fkThe frequency-domain covariance matrix of the sampled data of (a) can be expressed as a matrix of M × M dimensions as follows:
Figure BDA0002542026620000083
where L represents the fast beat number of the frequency domain data and xmk(l) Indicating that the m-th array element sampling data corresponds to a frequency point f after FFTkThe ith fast beat number of (c).
According to the initial angle set of N visible satellites, the frequency point f corresponding to the nth satellite can be calculatedkThe digital beam constraint equation of
Figure BDA0002542026620000084
Where Min represents the minimum, s.t. represents the constraint relationship, E represents the weighted average, w(n)(fk) For the nth satellite corresponds to frequency point fkThe weight vector of (2) is obtained by obtaining N groups of optimal weights according to a Lagrange multiplier method, and [ a ] is set as muH(θn,fk)R-1(fk)a(θn,fk)]-1Is a constant, then w(n)(fk)＝μR-1(fk)a(θn,fk) And obtaining the filtering weight vector of the nth wave beam.
Corresponding to the frequency domain fkThe parallel multi-beam output expression is:
Y(fk)＝[y(1)(fk) y(2)(fk) ... y(N)(fk)]T＝A(fk)[R-1(fk)]TX(fk)
A(fk)＝[a(θ1,fk) a(θ2,fk) ... a(θN,fk)]H
in the formula, A (f)k) Defined as a space-frequency incoming wave direction constraint matrix. The process of implementing space-frequency filtering is shown in fig. 1 for a wideband beamformer. Obtaining a frequency domain signal [ y ] output by the nth wave beam after space-frequency filtering(n)(f1) y(n)(f2) ... y(n)(fK)]N is 1, 2. The frequency domain output of each wave beam is subjected to inverse Fourier transform, and the time domain signal [ y ] of each satellite can be obtained(1)(t) y(2)(t) ... y(N)(t)]。
Respectively carrying out acquisition correlation peak processing on the time domain signal of each satellite to obtain the correlation peak size [ F ] of each satellite(1)(h) F(2)(h) ... F(N)(h)]And the related peak-to-peak function is subjected to iterative convergence, and the change factor is that the incoming wave direction set is as follows: theta ═ theta1,θ2,...,θNWhen the direction of each satellite wave is updated in a step manner of Δ θ, the function of the direction of the satellite wave of the satellite 1 can be expressed as:
Figure BDA0002542026620000091
wherein gamma ∈ (0, 1)]Is a learning factor. The other N-1 stars iterate the same, and the set of incoming wave directions returned to the digital beam processing is [ theta ]1(n+1),θ2(n+1),...,θN(n+1)]And calculating a space-frequency incoming wave direction constraint matrix by using the updated incoming wave direction set, repeating the space-frequency filtering processing, and iteratively converging to a maximum value of a related peak according to the rising gradient of a related peak-to-peak function, thereby obtaining the optimal processing of the digital wave beam under the parallel single-satellite constraint.
As shown in fig. 4, an embodiment of the present application further provides a digital multi-beam antenna processing apparatus, including:
a signal conversion module 21, configured to convert a received analog intermediate frequency signal into a digital intermediate frequency signal, and convert the digital intermediate frequency signal into frequency domain data through fourier transform;
a constraint matrix module 22, configured to determine a space-frequency incoming wave direction constraint matrix based on the incoming wave angle set of the satellite;
a weight vector module 23, configured to determine a filtering weight vector of each beam according to the space-frequency incoming wave direction constraint matrix and the frequency domain covariance matrix;
a filtering output module 24, configured to perform space-frequency filtering on the frequency domain data by using the filtering weight vector, and perform inverse fourier transform on the filtered frequency domain signal to obtain a time domain signal corresponding to each satellite;
and an updating module 25, configured to perform correlation peak detection on the time domain signal, and update an incoming wave angle set of the satellite based on a detection result.
In one embodiment, the constraint matrix module 22 is configured to:
determining a space-frequency two-dimensional steering vector of a frequency point corresponding to each satellite signal according to an incoming wave angle set of the satellite;
and constructing the space-frequency incoming wave direction constraint matrix based on the space-frequency two-dimensional steering vector.
In an embodiment, in the incoming wave angle set of the satellites, an initial incoming wave direction angle of each satellite is 0 degree.
In an embodiment, the update module 25 is configured to:
carrying out correlation peak detection on the time domain signal to obtain a detection result which is a correlation peak function of each satellite;
and correcting the incoming wave angle of each satellite based on the correlation peak function to obtain an updated incoming wave angle set.
In an embodiment, the update module 25 is configured to:
correcting the incoming wave angle theta of the nth satellite according to the following formulan(h+1)：
Figure BDA0002542026620000101
Wherein h represents the discrete point sequence number of the time, gamma is a learning factor, and gamma belongs to (0, 1)]，F[θn(h)]Denotes thetan(h) Correlation peak function of direction.
The embodiment of the application provides a digital multi-beam antenna processing method and device for a broadband beam former. On the basis of analyzing a broadband beam former model, received data of each array element is converted into a frequency domain, initial constraint is carried out on normal beams through a frequency domain parallel one-star constraint equation, the sizes of related peaks after stepping pointing angles are compared and gradually converged in the direction of an expected signal, the output signal-to-noise ratio of each frequency point in the signal bandwidth is improved, and the attenuation of a traditional algorithm to a useful signal is effectively prevented. Finally, each satellite signal is inversely transformed to the time domain, so that all satellite signals after interference suppression are obtained.
An embodiment of the present application further provides a digital multi-beam antenna processing apparatus, including: a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the digital multi-beam antenna processing method when executing the program.
Embodiments of the present application also provide a computer-readable storage medium storing computer-executable instructions for performing the digital multi-beam antenna processing method.
In this embodiment, the storage medium may include, but is not limited to: a U-disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a removable hard disk, a magnetic or optical disk, and other various media capable of storing program codes.
It will be understood by those of ordinary skill in the art that all or some of the steps of the methods, systems, functional modules/units in the devices disclosed above may be implemented as software, firmware, hardware, and suitable combinations thereof. In a hardware implementation, the division between functional modules/units mentioned in the above description does not necessarily correspond to the division of physical components; for example, one physical component may have multiple functions, or one function or step may be performed by several physical components in cooperation. Some or all of the components may be implemented as software executed by a processor, such as a digital signal processor or microprocessor, or as hardware, or as an integrated circuit, such as an application specific integrated circuit. Such software may be distributed on computer readable media, which may include computer storage media (or non-transitory media) and communication media (or transitory media). The term computer storage media includes volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data, as is well known to those of ordinary skill in the art. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, Digital Versatile Disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can accessed by a computer. In addition, communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media as known to those skilled in the art.

Claims (8)

1. A digital multi-beam antenna processing method, comprising:
converting a received analog intermediate frequency signal into a digital intermediate frequency signal, and converting the digital intermediate frequency signal into frequency domain data through Fourier transformation;
determining a space frequency incoming wave direction constraint matrix based on an incoming wave angle set of a satellite, wherein the incoming wave angle set is preset with an initial value;
determining a filtering weight vector of each wave beam according to the space-frequency incoming wave direction constraint matrix and the frequency domain covariance matrix;
performing space-frequency filtering on the frequency domain data by using the filtering weight vector, and performing inverse Fourier transform on the filtered frequency domain signal to obtain a time domain signal corresponding to each satellite;
carrying out correlation peak detection on the time domain signal to obtain a detection result which is a correlation peak function of each satellite;
and correcting the incoming wave angle of each satellite based on the correlation peak function to obtain an updated incoming wave angle set.
2. The method of claim 1, wherein the satellite-based set of incoming wave angles determines a space-frequency incoming wave direction constraint matrix, comprising:
determining a space-frequency two-dimensional steering vector of a frequency point corresponding to each satellite signal according to an incoming wave angle set of the satellite;
and constructing the space-frequency incoming wave direction constraint matrix based on the space-frequency two-dimensional steering vector.
3. The method according to claim 1 or 2, wherein the preset initial values of the incoming wave angle sets comprise:
the initial incoming wave direction angle of each satellite is 0 degree.
4. The method of claim 1, wherein the modifying the angle of the incoming wave for each satellite based on the correlation peak function comprises:
correcting the incoming wave angle theta of the nth satellite according to the following formulan(h+1)：
Figure FDA0003269860500000021
Wherein h represents the discrete point sequence number of the time, gamma is a learning factor, and gamma belongs to (0, 1)]，F[θn(h)]Denotes thetan(h) Correlation peak function of direction.
5. A digital multi-beam antenna processing apparatus, comprising:
the signal conversion module is used for converting the received analog intermediate frequency signal into a digital intermediate frequency signal and converting the digital intermediate frequency signal into frequency domain data through Fourier transformation;
the constraint matrix module is used for determining a space frequency incoming wave direction constraint matrix based on an incoming wave angle set of a satellite;
the weight vector module is used for determining a filtering weight vector of each wave beam according to the space frequency incoming wave direction constraint matrix and the frequency domain covariance matrix;
the filtering output module is used for performing space-frequency filtering on the frequency domain data by using the filtering weight vector and performing inverse Fourier transform on the filtered frequency domain signal to obtain a time domain signal corresponding to each satellite;
the updating module is used for carrying out correlation peak detection on the time domain signal, and the obtained detection result is a correlation peak function of each satellite; and correcting the incoming wave angle of each satellite based on the correlation peak function to obtain an updated incoming wave angle set.
6. The apparatus of claim 5,
in the incoming wave angle set of the satellites, the initial incoming wave direction angle of each satellite is 0 degree.
7. A digital multi-beam antenna processing apparatus, comprising: memory, processor and computer program stored on the memory and executable on the processor, characterized in that the processor implements the method according to any of claims 1 to 4 when executing the program.
8. A computer-readable storage medium storing computer-executable instructions for performing the method of any one of claims 1-4.
CN202010549747.8A 2020-06-16 2020-06-16 Digital multi-beam antenna processing method, device and equipment Active CN111934728B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010549747.8A CN111934728B (en) 2020-06-16 2020-06-16 Digital multi-beam antenna processing method, device and equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010549747.8A CN111934728B (en) 2020-06-16 2020-06-16 Digital multi-beam antenna processing method, device and equipment

Publications (2)

Publication Number Publication Date
CN111934728A CN111934728A (en) 2020-11-13
CN111934728B true CN111934728B (en) 2022-02-11

Family

ID=73316180

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010549747.8A Active CN111934728B (en) 2020-06-16 2020-06-16 Digital multi-beam antenna processing method, device and equipment

Country Status (1)

Country Link
CN (1) CN111934728B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117250644B (en) * 2023-11-08 2024-02-13 中国船舶集团有限公司第七〇七研究所 Ship application-oriented anti-interference method and system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102879790A (en) * 2011-07-13 2013-01-16 北京泰豪联星技术有限公司 Anti-interference system and method based on digital beam forming and space-time zeroing cascade
CN108650007A (en) * 2018-05-03 2018-10-12 西安交通大学 A kind of highly reliable channel equalization method based on null tone adaptive-filtering

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012168157A (en) * 2011-02-11 2012-09-06 National Univ Corp Shizuoka Univ Multi-beam type radar apparatus for on-vehicle use, multi-beam type radar method and multi-beam type radar program
US10348394B1 (en) * 2014-03-14 2019-07-09 Tarana Wireless, Inc. System architecture and method for enhancing wireless networks with mini-satellites and pseudollites and adaptive antenna processing
CN105375938A (en) * 2014-10-16 2016-03-02 航天恒星科技有限公司 Signal processing method based on inertial navigation auxiliary space frequency, apparatus and receiver
CN104280745B (en) * 2014-10-28 2017-04-19 北京北斗星通导航技术股份有限公司 Improved space-frequency signal processing method and processing device
CN105158735B (en) * 2015-06-03 2017-06-16 西安电子科技大学 Null tone Two-Dimensional Spectral Estimation method based on compression sampling array
CN106646531B (en) * 2016-11-16 2019-05-17 和芯星通科技(北京)有限公司 A kind of more stars constrain steady null tone anti-interference processing method and device
CN111277310B (en) * 2020-01-21 2023-06-20 北京北斗星通导航技术股份有限公司 Blind beam pointing airspace filtering processing method, device and equipment

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102879790A (en) * 2011-07-13 2013-01-16 北京泰豪联星技术有限公司 Anti-interference system and method based on digital beam forming and space-time zeroing cascade
CN108650007A (en) * 2018-05-03 2018-10-12 西安交通大学 A kind of highly reliable channel equalization method based on null tone adaptive-filtering

Also Published As

Publication number Publication date
CN111934728A (en) 2020-11-13

Similar Documents

Publication Publication Date Title
CN103954938B (en) Multi-sub-band reception synthetic method for SAR echo signals
US7885688B2 (en) Methods and systems for signal selection
WO2007103037A2 (en) System and method for generating a separated signal
CN107431525B (en) Receiving apparatus and receiving method
KR101603630B1 (en) Method and Apparatus for suppressing jammer signals and estimating Angle Of Arrival of original signal using orthogonal of transmitting signal waveform
CN103346756A (en) Space-time self-adaptive filtering method and device
CN111934728B (en) Digital multi-beam antenna processing method, device and equipment
CN104777491A (en) Blind beam broadband interference suppression method and device
US9444558B1 (en) Synthetic robust adaptive beamforming
CN104716982A (en) Robust anti-interference processing method and device of spread frequency system
Bilgehan et al. Fast detection and DOA estimation of the unknown wideband signal sources
CN103701515A (en) Digital multi-beam forming method
US8798569B2 (en) Channel estimation method and device in a multi-antenna system
Lee et al. Large-array signal processing for deep-space applications
CN109188366A (en) Broadband emission Adaptive beamformer method based on subband maximum signal noise ratio principle
CN111446997B (en) Self-adaptive digital beam synthesis method based on deep learning
CN111277310B (en) Blind beam pointing airspace filtering processing method, device and equipment
CN113419241B (en) Signal reconstruction method and equipment
CN111934650A (en) Spatial domain adaptive filtering processing method, device and equipment
CN104483665B (en) Beam forming method and beam forming system of passive acoustic sensor array
CN113075613A (en) Interference direction finding method in satellite mobile communication system
JP7168446B2 (en) Beamforming architecture for scalable radio frequency frontends
EP3712626B1 (en) High-rate dft-based data manipulator and data manipulation method for high performance and robust signal processing
CN117452442B (en) High-precision satellite navigation anti-interference method based on space-frequency self-adaptive processing
WO2024031880A1 (en) Beam decoupling method for multiple beams, transmitting end and receiving end

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
TA01 Transfer of patent application right

Effective date of registration: 20210521

Address after: 100094 Beidou Xingtong building, 7 Fengxian East Road, Yongfeng high tech industrial base, Zhongguancun, Haidian District, Beijing

Applicant after: BEIJING BDSTAR NAVIGATION Co.,Ltd.

Address before: 100094 room A301, Beidou Xingtong building, 7 Fengxian East Road, Haidian District, Beijing

Applicant before: Hexin Xingtong Technology (Beijing) Co.,Ltd.

Applicant before: UNICORECOMM (SHANGHAI) TECHNOLOGY Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant