CN108957411A - High-precision array signal generation method - Google Patents
High-precision array signal generation method Download PDFInfo
- Publication number
- CN108957411A CN108957411A CN201810473018.1A CN201810473018A CN108957411A CN 108957411 A CN108957411 A CN 108957411A CN 201810473018 A CN201810473018 A CN 201810473018A CN 108957411 A CN108957411 A CN 108957411A
- Authority
- CN
- China
- Prior art keywords
- array
- signal
- array element
- frequency spectrum
- incoming 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
- G01S7/4021—Means for monitoring or calibrating of parts of a radar system of receivers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/23—Testing, monitoring, correcting or calibrating of receiver elements
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/20—Monitoring; Testing of receivers
- H04B17/29—Performance testing
Abstract
The invention relates to the technical field of array signal processing, in particular to a high-precision array signal generation method, which comprises the following steps: (S1) performing fourier transform on the incident signal to obtain a frequency spectrum of the incident signal; (S2) establishing a northeast coordinate system by taking the reference receiving point as a coordinate origin, and calculating the time delay difference of the incident signal reaching each array element according to the incident direction of the incident signal and the coordinates of each array element of the antenna array; (S3) performing a phase rotation operation on the frequency spectrum of the incident signal to generate a frequency spectrum of a received signal of each array element of the antenna array; (S4) inverse fourier transform is performed on the frequency spectrum of each array element received signal of the antenna array to obtain the time domain form of each array element received signal, i.e. the final array signal. The array signal generating method can accurately generate the array signal required by the wired test of the array receiving equipment including wireless communication and satellite navigation, and is used for accurately evaluating the anti-interference performance of the array receiving equipment.
Description
Technical field
The present invention relates to the generation methods of array signal processing technology, especially array signal, more particularly relate to
And a kind of high-precision array for carrying out wired test to radar, wireless communication and satellite navigation array received equipment is believed
Number generation method.
Background technique
Currently, antenna array is widely used to the fields such as radar, wireless communication and satellite navigation, conventional reception is set
Standby to use single-antenna received signal, when receiving device is by intentional or malice electromagnetic interference, receiving device can only be from time domain
Or frequency domain inhibits interference.However, time domain and solution in frequency domain method can only inhibit narrowband to interfere, when electromagnetic interference
When bandwidth is more than the 10% of echo signal bandwidth, interference suppressioning effect is unsatisfactory.When receiving device receives letter using antenna array
Number when, as long as electromagnetic interference is different from the incident direction of echo signal, receiving device can all have interference from airspace
Effect inhibits, and insensitive to the bandwidth of interference.In addition, receiving device can also carry out Wave beam forming processing using antenna array to mention
The signal-to-noise ratio of high echo signal, to obtain better receptivity.Based on these advantages, wireless communication high-end at present with
And satellite navigation receiving equipment especially military equipment mostly uses greatly antenna array to receive signal.
The performance test means of array received equipment include wireless test and wired test, wherein wired test is in array
Proof-Of Principle and debugging stage in receiving device prototype design play an important role.Wired survey is carried out to array acceptor
Examination generates each array element received signal (i.e. array signal) of antenna array firstly the need of simulation.Since incoming signal reaches respectively
The time of a array element is different, and there are delay inequalitys for each array element received signal.
If the expression formula of incoming signal s (t) are as follows:
In formula, d (t) indicates baseband signal, f0Indicate the carrier frequency of incoming signal, j is imaginary unit, and e is nature pair
Several truth of a matter, π are pi, and t indicates the time.Incoming signal determines that can be echo signal can also according to actual test demand
To be various electromagnetic interferences.
Assuming that the array number of antenna array is N, then the signal that each array element of antenna array receives may be expressed as:
In formula, rk(t) signal received for k-th of array element, τkGinseng is reached for incoming signal k-th of array element of arrival is opposite
Examine the delay inequality of receiving point, k=1,2 ..., N, ()TIndicate transposition.R (t) is to need to simulate generation in wired test
Array signal.When the size of antenna array is smaller (such as diameter is the circle battle array of 200mm), above-mentioned delay inequality is generally measured in nanosecond
Grade, therefore the difficult point that array signal generates is how to carry out accurate timing_delay estimation to incoming signal.
Traditional array signal generation method is mostly based on narrowband it is assumed that thinking that incoming signal reaches each array element
Delay inequality be much smaller than the inverse of incoming signal bandwidth, therefore neglected the delay inequality of baseband signal d (t), array is believed at this time
It number can be approximately:
In formula,For a N-dimensional vector, in array signal
It is known as steering vector in processing.From above formula, it can be seen that, when being assumed using narrowband, the realization that array signal generates is relatively simple,
It does not need to carry out timing_delay estimation to incoming signal, it is only necessary to steering vector is calculated, then by steering vector and incoming signal s
(t) it is multiplied.
When the broader bandwidth of incoming signal, the approximate processing that narrowband is assumed will introduce biggish error.Nie Junwei etc.
Result of study in the impact analysis that narrowband assumes to GNSS antenna battle array anti-jamming performance evaluation shows (bibliography: Nie Jun
Big, Ge Rui, Li Min etc., " assuming the impact analysis to GNSS antenna battle array anti-jamming performance evaluation in narrowband ", the National University of Defense technology is learned
Report, in October, 2011, vol.33, no.5, pp.128-133.), when assessing the performance of the anti-broadband interference of array received equipment,
If generating array signal according to traditional array signal generation method, assessment result will be partially optimistic, not can accurately reflect
The actual performance of equipment.
Summary of the invention
In view of the defects existing in the prior art, the object of the present invention is to provide a kind of high-precision array signal generation method,
When for carrying out wired test to the array received equipment including wireless communication, satellite navigation, required survey is accurately generated
Trial signal.Specific technical solution is as follows:
A kind of high-precision array signal generation method, it is characterised in that include the following steps:
(S1) Fourier transformation is carried out to incoming signal s (t), obtains the frequency spectrum S (f) of incoming signal:
In formula, f indicates frequency, T1For the duration of incoming signal.
(S2), as coordinate origin, to establish northeast day coordinate system (East-North-Up, abbreviation: ENU) with reference to receiving point,
And incoming signal is calculated according to the coordinate of the incident direction of incoming signal and each array element of antenna array in ENU coordinate system and is reached
The delay inequality τ of each array elementk, k=1,2 ..., N:
In formula, (xk,yk,zk) it is coordinate of k-th of the array element in ENU coordinate system, θ is the pitch angle of incoming signal,
For the azimuth of incoming signal, c is the light velocity；
(S3) phase rotation operation is carried out to the frequency spectrum S (f) of incoming signal, generates each array element of antenna array and receives signal
Frequency spectrum, by taking k-th of (k=1,2 ..., N) array element as an example, the frequency spectrum for receiving signal is obtained by the phase rotation operation of following formula
It arrives:
In formula, Yk(f) frequency spectrum of signal is received for k-th of array element.
(S4) frequency spectrum for receiving signal to each array element of antenna array carries out Fourier inversion, obtains each array element and receives letter
Number forms of time and space, still with k-th of (k=1,2 ..., N) array element for, receive signal forms of time and space obtained by following formula
It arrives:
In formula, B is the bandwidth of incoming signal.yk(t) it is k-th of component of the array signal ultimately generated, ultimately generates
Array signal indicate are as follows:
Y (t)=[y1(t) … yk(t) … yN(t)]T。
Preferably, described with reference to the geometric center or some array element position that receiving point is antenna array.
Using the present invention obtain the utility model has the advantages that the present invention can accurately generate to include wireless communication, satellite navigation exists
Interior array received equipment carries out array signal required when wired test, for the anti-dry of accurate evaluation array received equipment
Immunity energy.
Detailed description of the invention
Fig. 1 is the flow chart of array signal generation method of the present invention；
Fig. 2 is the line array schematic diagram of four array elements in embodiment.
Specific embodiment
Below in conjunction with drawings and examples, the present invention is described further.
Fig. 1 is the flow chart of array signal generation method of the present invention, as shown, the present invention the following steps are included:
Step S1 carries out Fourier transformation to incoming signal, obtains the frequency spectrum of incoming signal.The frequency spectrum table of incoming signal
It is shown as:
In formula, f indicates frequency, T1For the duration of incoming signal, j indicates imaginary unit, and e indicates the bottom of natural Exponents, π
For pi, t indicates the time.Incoming signal s (t) can be generated by general signal source, can be public using Agilent in embodiment
The vector signal generator E4438C of department.
Step S2, as origin, to calculate the delay inequality that incoming signal reaches each array element of antenna array with reference to receiving point.First
, as origin, to establish ENU coordinate system with reference to receiving point.Then, according to the incident direction of incoming signal and each array element of antenna array
Coordinate in ENU coordinate system calculates the delay inequality that incoming signal reaches each array element, with k-th of (k=1,2 ..., N) battle array
For member, the delay inequality for receiving signal is calculate by the following formula:
In formula, (xk,yk,zk) it is coordinate of k-th of the array element in ENU coordinate system, θ is the pitch angle of incoming signal,
For the azimuth of incoming signal, c is the light velocity；
Antenna array selection quaternary line array in the present embodiment, i.e. N=4, as shown in Figure 2.It chooses array element 1 and (uses A1 table in figure
Showing array element 1, remaining array element is A2, A3, A4) position is with reference to receiving point to establish ENU coordinate system, it is calculate by the following formula kth
A (k=1,2,3,4) array element receives the delay inequality of signal:
Wherein, d is array element spacing, and θ is the pitch angle of incoming signal, the azimuth of incoming signalIt is light equal to 0, c
Speed；
Step S3 carries out phase rotation operation to the frequency spectrum of incoming signal, obtains the frequency spectrum that each array element receives signal.With
For k-th of (k=1,2 ..., N) array element, the expression formula of phase rotation operation is as follows:
In formula, Yk(f) frequency spectrum of signal is received for k-th of (k=1,2 ..., N) array element.
Step S4, the frequency spectrum for receiving signal to each array element carry out Fourier inversion, obtain final array signal.With
For k-th of (k=1,2 ..., N) array element, firstly, obtaining its time domain shape for receiving signal by the Fourier inversion of following formula
Formula:
In formula, B is the bandwidth of incoming signal.yk(t) k-th of component of the array signal as ultimately generated, most throughout one's life
At array signal are as follows:
Y (t)=[y1(t) … yk(t) … yN(t)]T；
In the present embodiment, the array signal ultimately generated are as follows:
Y (t)=[y1(t) y2(t) y3(t) y4(t)]T。
Contain the explanation of the preferred embodiment of the present invention above, this be for the technical characteristic that the present invention will be described in detail, and
Be not intended to for summary of the invention being limited in concrete form described in embodiment, according to the present invention content purport carry out its
He is also protected by this patent modifications and variations.The purport of the content of present invention is to be defined by the claims, rather than by embodiment
Specific descriptions defined.
Claims (4)
1. a kind of high-precision array signal generation method, it is characterised in that include the following steps:
(S1) Fourier transformation is carried out to incoming signal, obtains the frequency spectrum of incoming signal；
(S2), as coordinate origin, to establish northeast day coordinate system with reference to receiving point, and according to the incident direction and day of incoming signal
The coordinate of each array element of linear array calculates the delay inequality that incoming signal reaches each array element；
(S3) phase rotation operation is carried out to the frequency spectrum of incoming signal, generates the frequency spectrum that each array element of antenna array receives signal；
(S4) frequency spectrum for receiving signal to each array element of antenna array carries out Fourier inversion, obtain each array element receive signal when
Domain form, i.e., final array signal.
2. a kind of high-precision array signal generation method as described in claim 1, it is characterised in that: described to be with reference to receiving point
The geometric center of antenna array or some array element position.
3. a kind of high-precision array signal generation method as described in claim 1, it is characterised in that: the calculating of the delay inequality
Formula are as follows:
Wherein, τkIndicate that incoming signal reaches the delay inequality of k-th of array element, (xk,yk,zk) it is k-th of array element in northeast day coordinate
Coordinate in system, θ are the pitch angle of incoming signal,For the azimuth of incoming signal, c is the light velocity.
4. a kind of high-precision array signal generation method as claimed in claim 3, it is characterised in that: the tool of the step (S3)
Body process are as follows:
In formula, Yk(f) frequency spectrum of signal is received for k-th of array element, S (f) is the frequency spectrum of incoming signal, and f indicates frequency, and j indicates empty
Number unit, the value range of k are 1,2 ..., N；N is the array element total number of antenna array.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810473018.1A CN108957411A (en) | 2018-05-17 | 2018-05-17 | High-precision array signal generation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810473018.1A CN108957411A (en) | 2018-05-17 | 2018-05-17 | High-precision array signal generation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108957411A true CN108957411A (en) | 2018-12-07 |
Family
ID=64499183
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810473018.1A Pending CN108957411A (en) | 2018-05-17 | 2018-05-17 | High-precision array signal generation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108957411A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113111470A (en) * | 2021-04-13 | 2021-07-13 | 中国民用航空飞行学院 | Arraying method for identifying number features of targets of low-altitude small unmanned aerial vehicle |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0326906D0 (en) * | 2003-11-19 | 2003-12-24 | Roke Manor Research | A method of calibrating an adaptive antenna array of a satellite system |
| CN1606299A (en) * | 2003-10-10 | 2005-04-13 | 三星电子株式会社 | Apparatus, method and program for correcting common phase error of OFDM signal symbols |
| CN102062851A (en) * | 2009-11-16 | 2011-05-18 | 佘彩云 | Direction finding method based on improved L array star-carrying broadband multipurpose |
| CN102355443A (en) * | 2011-08-18 | 2012-02-15 | 广州海格通信集团股份有限公司 | Method for realizing federated filtering and timing synchronization in digital communication system |
| CN103297379A (en) * | 2013-07-01 | 2013-09-11 | 重庆邮电大学 | Time-varying orthogonal frequency division multiplexing (TV-OFDM) multi-carrier modulation system and modulation method |
| CN105549037A (en) * | 2015-12-18 | 2016-05-04 | 中国人民解放军国防科学技术大学 | High-precision satellite navigation broadband array signal generation method |
| CN106210719A (en) * | 2015-05-29 | 2016-12-07 | 美国西门子医疗解决公司 | Channel data based on Ultrasound beamforming device compresses |
| CN106443573A (en) * | 2016-09-14 | 2017-02-22 | 清华大学 | Angle-of-arrival single snapshot quick estimation method and device for massive antenna array signals |
-
2018
- 2018-05-17 CN CN201810473018.1A patent/CN108957411A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1606299A (en) * | 2003-10-10 | 2005-04-13 | 三星电子株式会社 | Apparatus, method and program for correcting common phase error of OFDM signal symbols |
| GB0326906D0 (en) * | 2003-11-19 | 2003-12-24 | Roke Manor Research | A method of calibrating an adaptive antenna array of a satellite system |
| CN102062851A (en) * | 2009-11-16 | 2011-05-18 | 佘彩云 | Direction finding method based on improved L array star-carrying broadband multipurpose |
| CN102355443A (en) * | 2011-08-18 | 2012-02-15 | 广州海格通信集团股份有限公司 | Method for realizing federated filtering and timing synchronization in digital communication system |
| CN103297379A (en) * | 2013-07-01 | 2013-09-11 | 重庆邮电大学 | Time-varying orthogonal frequency division multiplexing (TV-OFDM) multi-carrier modulation system and modulation method |
| CN106210719A (en) * | 2015-05-29 | 2016-12-07 | 美国西门子医疗解决公司 | Channel data based on Ultrasound beamforming device compresses |
| CN105549037A (en) * | 2015-12-18 | 2016-05-04 | 中国人民解放军国防科学技术大学 | High-precision satellite navigation broadband array signal generation method |
| CN106443573A (en) * | 2016-09-14 | 2017-02-22 | 清华大学 | Angle-of-arrival single snapshot quick estimation method and device for massive antenna array signals |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113111470A (en) * | 2021-04-13 | 2021-07-13 | 中国民用航空飞行学院 | Arraying method for identifying number features of targets of low-altitude small unmanned aerial vehicle |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109298388B (en) | Direct over-the-horizon target geographic coordinate estimation method based on azimuth information | |
| CN106950529B (en) | Acoustic vector near field sources ESPRIT and MUSIC method for parameter estimation | |
| CN102997988B (en) | Pool testing method of low-frequency acoustic directivity of large submerged buoy vector hydrophone | |
| CN105334266A (en) | Rock acoustic emission source positioning method | |
| WO2006068888A1 (en) | System and technique for calibrating radar arrays | |
| Park et al. | Compact HF surface wave radar data generating simulator for ship detection and tracking | |
| CN108776342A (en) | A kind of high speed platform SAR moving-target detection and speed estimation method at a slow speed | |
| CN112379342B (en) | Echo simulation and echo characteristic parameter precision estimation method for satellite-borne cloud detection radar | |
| CN103869298B (en) | A kind of distributed MIMO sky-wave OTH radar sea clutter emulation mode | |
| Li et al. | A Barankin-type bound on direction estimation using acoustic sensor arrays | |
| CN111198387A (en) | Space-time sampling navigation positioning method capable of resisting deception jamming | |
| CN108957411A (en) | High-precision array signal generation method | |
| Zhao et al. | FPGA-Based Real-Time Synchronous Parallel System for Underwater Acoustic Positioning and Navigation | |
| CN109884621B (en) | Radar altimeter echo coherent accumulation method | |
| CN104360335B (en) | Ionized layer parametric inversion method based on AIS target indication | |
| Willerton | Array auto-calibration | |
| Li et al. | A multi‐step method for passive broadband source localisation using a single vector sensor | |
| Xinghua et al. | Performance gain bounds of coherently combining multiple radars in a target-based calibration manner | |
| Shi et al. | Weighted direct position determination via the dimension reduction method for noncircular signals | |
| CN113534219A (en) | Beidou positioning outdoor target method based on multipath utilization | |
| Zhang et al. | A vortex electromagnetic imaging method based on two-dimensional matrix compensation for off-grid targets | |
| Li et al. | Underwater source localization in the presence of strong interference | |
| Wei et al. | Prediction of ship motion attitude from radial velocity of water particle using coherent S-band radar | |
| Yu et al. | Error Analysis on Ultra-Short Baseline Positioning System | |
| Cao et al. | Research on Electromagnetic Pulse Prediction for Sensitive Equipment on Ship |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181207 |