[go: up one dir, main page]

CN102148636B - Antenna calibration method and system - Google Patents

Antenna calibration method and system Download PDF

Info

Publication number
CN102148636B
CN102148636B CN 201110030019 CN201110030019A CN102148636B CN 102148636 B CN102148636 B CN 102148636B CN 201110030019 CN201110030019 CN 201110030019 CN 201110030019 A CN201110030019 A CN 201110030019A CN 102148636 B CN102148636 B CN 102148636B
Authority
CN
China
Prior art keywords
time domain
domain sequences
sequence
sampling
cyclic sequence
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
CN 201110030019
Other languages
Chinese (zh)
Other versions
CN102148636A (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.)
Datang Mobile Communications Equipment Co Ltd
Original Assignee
Datang Mobile Communications Equipment 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 Datang Mobile Communications Equipment Co Ltd filed Critical Datang Mobile Communications Equipment Co Ltd
Priority to CN 201110030019 priority Critical patent/CN102148636B/en
Publication of CN102148636A publication Critical patent/CN102148636A/en
Priority to PCT/CN2012/070518 priority patent/WO2012100699A1/en
Application granted granted Critical
Publication of CN102148636B publication Critical patent/CN102148636B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/20Monitoring; Testing of receivers
    • H04B17/21Monitoring; Testing of receivers for calibration; for correcting measurements
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842Weighted combining
    • H04B7/086Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Radio Transmission System (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

The embodiment of the invention relates to the technical field of wireless communication, in particular to an antenna calibration method and an antenna calibration system, which are used for solving the problem in the prior art that a sampling signal may be deviated in antenna calibration to produce remarkable influence on baseband signal reception to further cause the heavy loss of calibration sequence signals and reduce the accuracy of a calibration coefficient, beamforming performance and direction of arrival (DOA) estimation accuracy. The method provided by the embodiment of the inventioncomprises that: a transmission side adds a periodic sequence into an original time domain sequence to obtain a first time domain sequence; the transmission side transmits the first time domain sequence to a receiving side; the receiving side determines a calibration sequence according to a second time domain sequence obtained by sampling; and the receiving side performs antenna calibration according to the calibration sequence. The transmission side adds the periodic sequence into the original time domain sequence to reduce the remarkable influence of sampling deviations on the baseband signal reception in the antenna calibration system, thereby reducing the loss of the calibration sequence signals.

Description

A kind of method and system of antenna calibration
Technical field
The present invention relates to wireless communication technology field, particularly a kind of method and system of antenna calibration.
Background technology
Antenna system suppresses advantages such as signal interference, at TD-SCDMA (Time Division Synchronized Code Division Multiple Access because it has the raising cell coverage area; The time-division synchronization CDMA; CDMA, Code Division Multiple Access, code division multiple access) be extensive use of in the mobile communication system, and at LTE (Long Term Evolution, Long Term Evolution) and in LTE-A (LongTerm Evolution-Advanced, the long-term evolution upgrading) communication system can continue use.Each bay radio-frequency channel characteristic that requires when smart antenna uses to form smart antenna is consistent, there are differences because thereby the actual parameter of the radio-frequency devices of different array elements there are differences the radio-frequency channel magnitude-phase characteristics that causes each bay, therefore adopt in the system of smart antenna all to have the antenna calibration function.
As shown in Figure 1, intelligent antenna calibration system comprises: one or more radio-frequency transmissions passage and radio frequency reception channel; The transmitting-receiving coupling channel; Calibration receive path and transmission channel; The calibrating signal processor; Baseband signal processor etc.
The antenna calibration algorithm is by receipts calibrating signal or a calibrating signal are carried out channel estimating, time domain or the frequency domain response of estimated result being used as the radio-frequency channel calculate calibration factor, because sampling deviation can influence the reception of baseband signal when signal receives, thereby cause the loss of signal of calibrating sequence, finally cause the calculating of calibration factor inaccurate, influence wave beam forming performance and DOA (Direction Of Arrival, arrival direction) estimated accuracy.
In sum, sampled signal deviation can occur and then receiving baseband signal is produced significantly influence in the antenna calibration at present, thereby causes the calibrating sequence loss of signal serious, reduces accuracy, wave beam forming performance and the DOA estimated accuracy of calibration factor.
Summary of the invention
The method and system of a kind of antenna calibration that the embodiment of the invention provides, sampled signal deviation can occur and then receiving baseband signal is produced significantly influence in the antenna calibration that exists in the prior art in order to solve, thereby cause the calibrating sequence loss of signal serious, reduce the problem of accuracy, wave beam forming performance and the DOA estimated accuracy of calibration factor.
The method of a kind of antenna calibration that the embodiment of the invention provides comprises:
Transmitter side adds cyclic sequence in original time domain sequences, obtain first time domain sequences;
Described transmitter side sends described first time domain sequences to receiver side;
Described receiver side is determined calibrating sequence according to second time domain sequences that sampling obtains;
Described receiver side carries out antenna calibration according to calibrating sequence.
The system of a kind of antenna calibration that the embodiment of the invention provides comprises:
Sending module is used for adding cyclic sequence in original time domain sequences and obtains first time domain sequences, sends described first time domain sequences;
Receiver module is used for determining calibrating sequence according to second time domain sequences that sampling obtains, and carries out antenna calibration according to calibrating sequence.
Because transmitter side adds cyclic sequence in original time domain sequences, can be in antenna calibration system, reduce owing to sampling deviation produces significantly influence to receiving baseband signal, thereby reduce the loss of calibrating sequence signal, improve the accurate and antenna calibration precision of calibration factor; Further improve performance and the DOA estimated accuracy of wave beam forming.
Description of drawings
Fig. 1 is the intelligent antenna calibration system schematic diagram;
Fig. 2 is the method flow schematic diagram of embodiment of the invention antenna calibration;
Fig. 3 sends the Calibration Method schematic flow sheet for the embodiment of the invention;
Fig. 4 receives the Calibration Method schematic flow sheet for the embodiment of the invention;
Fig. 5 is the system configuration schematic diagram of embodiment of the invention antenna calibration.
Embodiment
Embodiment of the invention transmitter side sends by adding first time domain sequences that cyclic sequence obtains in the original time domain sequences, and receiver side is determined calibrating sequence and carried out antenna calibration according to the time domain sequences that sampling obtains.If sampling shifts to an earlier date or lags behind, understanding some signal sampling point does not gather, and transmitter side adds after the cyclic sequence in original time domain sequences, can reduce the loss that the data sampling point is subjected to, can in antenna calibration system, reduce because sampling deviation produces significantly influence to receiving baseband signal, thereby reduce the loss of calibrating sequence signal.
Below in conjunction with Figure of description the embodiment of the invention is described in further detail.
As shown in Figure 2, the method for embodiment of the invention antenna calibration comprises the following steps:
Step 201, transmitter side add cyclic sequence in original time domain sequences, obtain first time domain sequences.
Step 202, transmitter side send first time domain sequences to receiver side.
Step 203, receiver side are determined calibrating sequence according to second time domain sequences that sampling obtains.
Step 204.Receiver side carries out antenna calibration according to calibrating sequence.
In the step 201, transmitter side adds cyclic sequence in original time domain sequences position can have multiple, such as the front in original time domain sequences; Back in original time domain sequences; Front and back in original time domain sequences.Below the branch situation describe.
Situation one, transmitter side add cyclic sequence in the front of original time domain sequences, and this cyclic sequence can be called CP (Cyclic Prefix).
Wherein, transmitter side need be determined cyclic sequence before the step 201.
Concrete, at first transmitter side determines that the length of cyclic sequence equals L Cp1, L wherein Cp1Be positive integer, M+N≤L Cp1≤ L, L are the length of original time domain sequences, and M is the numerical value after rounding up divided by the sampling interval the maximum sampling lag deviation time, and N is that maximum sampling shifts to an earlier date the numerical value after the deviation time rounded up divided by the sampling interval, such as 5.4, is exactly 6 after rounding up; Such as being 6, be exactly 6 after rounding up, down together);
Transmitter side is with back L in the original time domain sequences then Cp1The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence the foremost of original time domain sequences to, forms first time domain sequences, specifically can represent with formula one.
y * ( i ) = y ( L - L cp 1 + i ) , 0 ≤ i ≤ L cp 1 - 1 y * ( i ) = y ( i - L cp 1 ) , L cp 1 ≤ i ≤ L + L cp 1 - 1 . . . . . . . . . . Formula one.
Wherein, L is original time domain sequences length, and y (i) is original time domain sequences, y *(i) be first time domain sequences.
Such as: back L in the original time domain sequences Cp1The position is ABCD, and then ABCD adds the foremost to original time domain sequences to as cyclic sequence, and first time domain sequences of composition is the original time domain sequences of ABCD+.
L Cp1Concrete numerical value set according to the sampling time delay scope of system.
Preferable, if cyclic sequence length is greater than sampling deviation, then data sampling point can not incur loss, sampling deviation is equivalent to add a cyclic shift to the influence of the back data of sampling, further reduce because sampling deviation produces significantly influence to receiving baseband signal, reduce the loss of calibrating sequence signal.
Accordingly, in the step 203, receiver side is known the length of cyclic sequence, i.e. L Cp1Size.According to L Cp1Size, and second time domain sequences that sampling obtains just can be determined calibrating sequence.Specifically there is dual mode to determine calibrating sequence.
Mode one, steps A 1, receiver side be from sampling initial time sampled data, and form second time domain sequences according to the order of sampling.
Steps A 2, receiver side are from the N of second time domain sequences that samples 1Individual sampling is lighted L sampling number of intercepting according to (sampled point of initial time is called the 0th sampled point, down together), and forms sequence to be adjusted according to the order of intercepting, wherein N 1Be positive integer, and N≤N 1≤ L Cp1-M.
Steps A 3, receiver side move L to sequence big-endian circulation to be adjusted Cp1-N 1The position obtains calibrating sequence.
Concrete, in the steps A 1, receiver side begins to receive time domain data from the sampling initial time, and forms the second time domain sequences y (i), i=0,1, L, L+L according to the order of sampling CP1-1.
In the steps A 2, receiver side is from the second time domain sequences N 1L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted according to the order of intercepting
Figure BDA0000045704790000051
I=0,1, L, L-1.
In the steps A 3, receiver side moves L to sequence big-endian circulation to be adjusted Cp1-N 1The position obtains calibrating sequence, and concrete calibrating sequence can be by formula two expressions.
Figure BDA0000045704790000052
Formula two.
Mode two, step B1, receiver side are delayed N from the sampling initial time 1* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences, wherein N according to the order of sampling 1Be positive integer, and N≤N 1≤ L Cp1-M, Ts are the sampling intervals.
Circulation moves L to the second time domain sequences big-endian for step B2, receiver side Cp1-N 1The position obtains calibrating sequence.
Concrete, among the step B1, receiver side is delayed N from the sampling initial time 1* Ts rises constantly and begins to receive time domain data, and forms the second time domain sequences y (i), i=0,1, L, L-1 according to the order of sampling.
Among the step B2, circulation moves L to receiver side to the second time domain sequences big-endian Cp1-N 1The position obtains
Figure BDA0000045704790000053
As the calibrating sequence after the sampling.
Preferable, after mode one and mode two are removed cyclic sequence above adopting, when the subsequent calculations calibration factor, need use channel estimation value H, can compensate this moment to the time-domain cyclic shift operation, thereby further reduce because sampling deviation produces significantly influence to receiving baseband signal, reduce the loss of calibrating sequence signal.
Concrete, in the step 204, receiver side is determined channel estimation value according to calibrating sequence, and compensates according to formula three, carries out antenna calibration according to the road estimated value after the compensation.
H 1=H * e J2 π kn... .... formula three.
Wherein, H 1Be the channel estimation value after the compensation,
Figure BDA0000045704790000054
K is the subcarrier sequence number, and H is k subcarrier channel estimation, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion;
In force, the frequency domain channel estimated value on each subcarrier of being determined by calibrating sequence according to three pairs of formula of receiver side compensates.
Situation two, transmitter side add cyclic sequence in the back of original time domain sequences, and this cyclic sequence can be called cyclic suffix.
Wherein, transmitter side need be determined cyclic sequence before the step 201.
Concrete, at first transmitter side determines that the length of cyclic sequence equals L Cp2, L Cp2Be positive integer, M+N≤L wherein Cp2≤ L, L are the length of original time domain sequences, and M is the maximum sampling lag deviation time divided by the numerical value in sampling interval, and N is the maximum sampling numerical value that rounds up divided by the sampling interval of deviation time in advance;
Transmitter side is with preceding L in the original time domain sequences then Cp2The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence to original time domain sequences backmost, forms first time domain sequences, specifically can represent with formula four.
y * ( i ) = y ( i ) , 0 ≤ i ≤ L - 1 y * ( i ) = y ( i - L ) , L ≤ i ≤ L + L cp 2 - 1 . . . . . . . . . . Formula four.
Wherein, L is original time domain sequences length, and y (i) is original time domain sequences, y *(i) be first time domain sequences.
Such as: L before in the original time domain sequences Cp2The position is ABC, and then ABC adds to original time domain sequences backmost as cyclic sequence, and first time domain sequences of composition is original time domain sequences+ABC.
The concrete numerical value of Lcp2 is set according to the sampling time delay scope of system.
Preferable, if cyclic sequence length is greater than sampling deviation, then data sampling point can not incur loss, sampling deviation is equivalent to add a cyclic shift to the influence of the back data of sampling, further reduce because sampling deviation produces significantly influence to receiving baseband signal, reduce the loss of calibrating sequence signal.
Accordingly, in the step 203, receiver side is known the length of cyclic sequence, i.e. L Cp2Size.According to L Cp2Size, and second time domain sequences that sampling obtains just can be determined calibrating sequence.Specifically there are four kinds of modes to determine calibrating sequence.
Mode one, step C1, receiver side be from sampling initial time sampled data, and form second time domain sequences according to the order of sampling.
Step C2, receiver side are from the N of second time domain sequences 2L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted, wherein N according to the order of intercepting 2Be positive integer, and N≤N 2≤ L Cp2-M.
Step C3, receiver side move N to sequence little-endian circulation to be adjusted 2The position obtains calibrating sequence.
Concrete, among the step C1, receiver side begins to receive time domain data from the sampling initial time, and forms the second time domain sequences y (i), i=0,1, L, L+L according to the order of sampling CP2-1;
Among the step C2, receiver side is from the second time domain sequences N 2L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted according to the order of intercepting
Figure BDA0000045704790000071
I=0,1, L, L-1.
Among the step C3, receiver side is circulated to a high position by low level to sequence to be adjusted and moves N 2The position obtains calibrating sequence, and concrete calibrating sequence can be by formula five expressions.
Figure BDA0000045704790000072
Formula five.
Mode two, step D1, receiver side are delayed N from the sampling initial time 2* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences, wherein N according to the order of sampling 2Be positive integer, and N≤N 2≤ L Cp2-M, Ts are the sampling intervals.
Circulation moves N to the second time domain sequences little-endian for step D2, receiver side 2The position obtains calibrating sequence.
Concrete, among the step D1, receiver side is delayed N from the sampling initial time 2* Ts rises constantly and begins to receive time domain data, and forms the second time domain sequences y (i), i=0,1, L, L-1 according to the order of sampling.
Among the step D2, circulation moves N to receiver side to the second time domain sequences little-endian 2The position obtains
Figure BDA0000045704790000073
As the calibrating sequence after the sampling.
Preferable, after mode one and mode two are removed cyclic sequence above adopting, when the subsequent calculations calibration factor, need use channel estimation value H, can compensate this moment to the time-domain cyclic shift operation, thereby further reduce because sampling deviation produces significantly influence to receiving baseband signal, reduce the loss of calibrating sequence signal.
Concrete, in the step 204, receiver side is determined channel estimation value according to calibrating sequence, and compensates according to formula six, carries out antenna calibration according to the road estimated value after the compensation.
H 2=H * e -j2 π kn... .... formula six.
Wherein, H 2Be the channel estimation value after the compensation,
Figure BDA0000045704790000081
K is the subcarrier sequence number, and H is k subcarrier channel estimation, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion;
In force, the frequency domain channel estimated value on each subcarrier of being determined by calibrating sequence according to six pairs of formula of receiver side compensates.
Situation three, cyclic sequence comprise first cyclic sequence and second cyclic sequence, and transmitter side adds first cyclic sequence in the front of original time domain sequences, and transmitter side adds second cyclic sequence in the back of original time domain sequences.
Wherein, transmitter side need be determined cyclic sequence before the step 201.
Concrete, at first transmitter side determines that the length of first cyclic sequence equals L Cp3, and determine that the length of second cyclic sequence equals L Cp4, L wherein Cp3Be positive integer, L Cp4Be positive integer, N≤L Cp3≤ L, M≤L Cp4≤ L, L are the length of original time domain sequences, and M is the numerical value that rounds up divided by the sampling interval the maximum sampling lag deviation time, and N is that maximum sampling shifts to an earlier date the numerical value that the deviation time rounded up divided by the sampling interval;
Transmitter side is with back L in the original time domain sequences then Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, forms first time domain sequences, specifically can represent with formula seven.
y * ( i ) = y ( L - L cp 3 + i ) , 0 ≤ i ≤ L cp 3 - 1 y * ( i ) = y ( i - L cp 3 ) , L cp 3 ≤ i ≤ L + L cp 3 - 1 y * ( i ) = y ( i - L - L cp 3 ) , L + L cp 3 ≤ i ≤ L + L cp 3 + L cp 4 - 1 . . . . . . . . . . Formula seven.
Wherein, L is original time domain sequences length, and y (i) is original time domain sequences, y *(i) be first time domain sequences.
Such as: back L in the original time domain sequences Cp3The position is ABCD, and then ABCD adds the foremost to original time domain sequences to as first cyclic sequence, preceding L in the original time domain sequences Cp4The position is EFG, and then EFG adds to original time domain sequences backmost as second cyclic sequence, and first time domain sequences of composition is the original time domain sequences+EFG of ABCD+.
L Cp1Concrete numerical value set according to the sampling time delay scope of system.
Preferable, if cyclic sequence length is greater than sampling deviation, then data sampling point can not incur loss, sampling deviation is equivalent to add a cyclic shift to the influence of the back data of sampling, further reduce because sampling deviation produces significantly influence to receiving baseband signal, reduce the loss of calibrating sequence signal.
Accordingly, in the step 203, receiver side is known the length of first cyclic sequence, i.e. L Cp3Size, the length of second cyclic sequence, i.e. L Cp4Size.According to L Cp3And L Cp4Size, and second time domain sequences that sampling obtains just can be determined calibrating sequence.Specifically there is dual mode to determine calibrating sequence.
Mode one, step e 1, receiver side be from sampling initial time sampled data, and form second time domain sequences according to the order of sampling.
Step e 2, receiver side are from the L of second time domain sequences Cp3L sampling number certificate of intercepting lighted in individual sampling, and forms calibrating sequence according to the order of intercepting.
Concrete, in the step e 1, receiver side begins to receive time domain data from the sampling initial time, and forms the second time domain sequences y (i), 0≤i≤L according to the order of sampling Cp3+ L Cp4
In the step e 2, receiver side is from the second time domain sequences L Cp3L sampling number certificate of intercepting lighted in individual sampling, and according to the order composition sequence that intercepts
Figure BDA0000045704790000091
I=0,1, L, L-1; Receiver side is with sequence As the calibrating sequence after the sampling.
Mode two, step F 1, receiver side are delayed L from the sampling initial time Cp3* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences according to the order of sampling, and wherein Ts is the sampling interval.
Step F 2, receiver side with second time domain sequences as calibrating sequence.
Concrete, in the step F 1, receiver side is delayed L from the sampling initial time Cp3* Ts rises constantly and begins to receive time domain data, and forms the second time domain sequences y (i), i=0,1, L, L-1 according to the order of sampling.
In the step F 2, receiver side with sequences y (i) as the sampling after calibrating sequence.
Mode three, step G1, receiver side be from sampling initial time sampled data, and form second time domain sequences according to the order of sampling.
Step G2, receiver side are from the N of second time domain sequences 3L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted, wherein N≤N according to the order of intercepting 3≤ L Cp3+ L Cp4-M, and N 3≠ L Cp3
If step G3 is N 3<L Cp3, then receiver side moves L to sequence big-endian circulation to be adjusted Cp3-N 3The position obtains calibrating sequence;
If N 3>L Cp3, then receiver side moves N to sequence little-endian circulation to be adjusted 3-L Cp3The position obtains calibrating sequence.
Preferable, after mode three is removed cyclic sequence above adopting, when the subsequent calculations calibration factor, need use frequency domain channel estimated value H, can compensate this moment to the time-domain cyclic shift operation, thereby further reduce because sampling deviation produces significantly influence to receiving baseband signal, reduce the loss of calibrating sequence signal.
Concrete, if N 3<L Cp3, then receiver side is determined channel estimation value according to calibrating sequence, and compensates according to formula eight, carries out antenna calibration according to the road estimated value after the compensation.
H 3=H * e J2 π kn... .... formula eight.
Wherein, H 1Be the channel estimation value after the compensation,
Figure BDA0000045704790000101
K is the subcarrier sequence number, and H is k subcarrier channel estimation, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion;
In force, the frequency domain channel estimated value on each subcarrier of being determined by calibrating sequence according to eight pairs of formula of receiver side compensates.
If N 3>L Cp3, then receiver side is determined channel estimation value according to calibrating sequence, and compensates according to formula nine, carries out antenna calibration according to the road estimated value after the compensation:
H 3=H * e -j2 π kn... .... formula nine.
Wherein, H 3Be the frequency domain channel estimated value after the compensation, K is the subcarrier sequence number, and H is the frequency domain channel estimated value of k subcarrier, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion.
In force, the frequency domain channel estimated value on each subcarrier of being determined by calibrating sequence according to nine pairs of formula of receiver side compensates.
Mode four, step H1, receiver side are delayed N from the sampling initial time 3* Ts constantly rises and begins to receive time domain data, and forms second time domain sequences, wherein N≤N according to the order of sampling 3≤ L Cp3+ L Cp4-M, and N 3≠ L Cp3, Ts is the sampling interval.
If step H2 is N 3<L Cp3, then circulation moves L to receiver side to the second time domain sequences big-endian Cp3-N 3The position obtains calibrating sequence;
If N 3>L Cp3, then receiver side moves N to sequence little-endian circulation to be adjusted 3-L Cp3The position obtains calibrating sequence.
Preferable, after mode four is removed cyclic sequence above adopting, when the subsequent calculations calibration factor, need use frequency domain channel estimated value H, can compensate this moment to the time-domain cyclic shift operation, thereby further reduce because sampling deviation produces significantly influence to receiving baseband signal, reduce the loss of calibrating sequence signal.
Concrete, if N 3<L Cp3, N in the method that then specifically compensates and the situation three 3<L Cp3The method that compensates is identical, does not repeat them here.
If N 3>L Cp3, N in the method that then specifically compensates and the situation three 3>L Cp3The method that compensates is identical, does not repeat them here.
The mode of which kind of interpolation cyclic sequence and the length of cyclic sequence can be stipulated in agreement above concrete the employing, are perhaps consulted to determine by transmitter side and receiver side; Perhaps be configured in transmitter side and receiver side in advance respectively.
If what carry out is to send calibration, then transmitter side is baseband signal processor, and receiver side is the calibrating signal processor; If what carry out is to receive calibration, then transmitter side is the calibrating signal processor, and receiver side is baseband signal processor.
Being calibrated to example with transmission calibration and reception respectively below describes.
As shown in Figure 3, embodiment of the invention transmission Calibration Method comprises the following steps:
Step 301, baseband signal processor generate the radio-frequency channel calibrating sequence and are mapped on the subcarrier in frequency domain.
Step 302, baseband signal processor convert the frequency domain calibrating signal to the time domain baseband signal by the IFFT conversion.
Step 303, baseband signal processor send by each radio-frequency channel, road after time domain baseband signal (being original time domain sequences) is added cyclic sequence.
Step 304, baseband signal arrive calibrated channel through receiving coupling network.
Step 305, calibrating signal processor receive time domain data (being the time domain sequence) after each road stack by calibrated channel, and determine to remove the time domain data (being calibrating sequence) of cyclic sequence according to the time domain data of receiving.
The time domain data that step 306, calibrating signal processor will be removed behind the cyclic sequence transforms to frequency domain by FFT (FastFourier Transform, fast fourier transform).
Step 307, calibrating signal processor estimate that according to receiving signal each road sends time domain or the frequency domain characteristic of radio-frequency channel, calculates calibration factor, and sends to baseband signal processor.
Step 308, baseband signal processor carry out compensating coefficient in the signal emission process, guarantee that each road sends the amplitude-phase unanimity of radio-frequency channel.
As shown in Figure 4, the method for embodiment of the invention reception calibration flow process comprises the following steps:
Step 401, calibrating signal processor generate the radio-frequency channel calibrating sequence and are mapped on the subcarrier in frequency domain.
Step 402, calibrating signal processor convert the frequency domain calibrating signal to the time domain baseband signal by the IFFT conversion.
Step 403, calibrating signal processor send by calibrated channel after time domain baseband signal (being original time domain sequences) is added cyclic sequence.
Step 404, baseband signal arrive each road receiving RF channel through sending coupling network.
Step 405, baseband signal processor receive time domain data (being the time domain sequence) by each road receiving RF channel, and determine to remove the time domain data (being calibrating sequence) of cyclic sequence according to the time domain data of receiving
The time domain data that step 406, baseband signal processor will be removed behind the cyclic sequence transforms to frequency domain by FFT.
Step 407, baseband signal processor are estimated time domain or the frequency domain characteristic of each road receiving RF channel according to receiving signal, calculate calibration factor.
Step 408, baseband signal processor carry out compensating coefficient at the signal receiving course, guarantee the amplitude-phase unanimity of each road receiving RF channel.
Based on same inventive concept, a kind of system of antenna calibration also is provided in the embodiment of the invention, because the principle that this system deals with problems is similar to the method for antenna calibration, so the enforcement of this system can repeat part and repeat no more referring to the enforcement of method.
As shown in Figure 5, the system of embodiment of the invention antenna calibration comprises: sending module 10 and receiver module 20.
Sending module 10 is used for adding cyclic sequence in original time domain sequences and obtains first time domain sequences, sends first time domain sequences.
Receiver module 20 is used for determining calibrating sequence according to second time domain sequences that sampling obtains, and carries out antenna calibration according to calibrating sequence.
Sending module 10 adds cyclic sequence in original time domain sequences position can have multiple, such as the front in original time domain sequences; Back in original time domain sequences; Front and back in original time domain sequences.Below the branch situation describe.
Situation one, sending module 10 determine that the length of cyclic sequence equals L Cp1, M+N≤L wherein Cp1≤ L, L are the length of original time domain sequences, and M is the numerical value that rounds up divided by the sampling interval the maximum sampling lag deviation time, and N is that maximum sampling shifts to an earlier date the numerical value that the deviation time rounded up divided by the sampling interval; With back L in the original time domain sequences Cp1The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence the foremost of original time domain sequences to, forms first time domain sequences.
Accordingly, receiver module 20 is from sampling initial time sampled data, and forms second time domain sequences according to the order of sampling; N from second time domain sequences 1L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted, wherein N according to the order of intercepting 1Be positive integer, and N≤N 1≤ L Cp1-M; Sequence big-endian circulation to be adjusted is moved L Cp1-N 1The position obtains calibrating sequence.Perhaps
Receiver module 20 is delayed N from the sampling initial time 1* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences, wherein N according to the order of sampling 1Be positive integer, and N≤N 1≤ L Cp1-M, Ts are the sampling intervals; Circulation moves L to the second time domain sequences big-endian Cp1-N 1The position obtains calibrating sequence.
Preferable, receiver module 20 is determined channel estimation value according to calibrating sequence; Compensate according to formula three; Carry out antenna calibration according to the channel estimation value after the compensation.
Situation two sending modules 10 determine that the length of cyclic sequence equals L Cp2, M+N≤L wherein Cp2≤ L, L are the length of original time domain sequences, and M is the numerical value that rounds up divided by the sampling interval the maximum sampling lag deviation time, and N is that maximum sampling shifts to an earlier date the numerical value that the deviation time rounded up divided by the sampling interval; With preceding L in the original time domain sequences Cp2The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence to original time domain sequences backmost, forms first time domain sequences.
Accordingly, receiver module 20 is from sampling initial time sampled data, and forms second time domain sequences according to the order of sampling; N from second time domain sequences 2L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted, wherein N according to the order of intercepting 2Be positive integer, and N≤N 2≤ L Cp2-M; Sequence little-endian circulation to be adjusted is moved N 2The position obtains calibrating sequence.Perhaps
Receiver module 20 is delayed N from the sampling initial time 2* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences, wherein N according to the order of sampling 2Be positive integer, and N≤N 2≤ L Cp2-M, Ts are the sampling intervals; Circulation moves N to the second time domain sequences little-endian 2The position obtains calibrating sequence.
Preferable, receiver module 20 is determined channel estimation value according to calibrating sequence; Compensate according to formula six; Carry out antenna calibration according to the channel estimation value after the compensation.
Situation three, cyclic sequence comprise first cyclic sequence and second cyclic sequence; Sending module 10 determines that the length of first cyclic sequence equals L Cp3, and determine that the length of second cyclic sequence equals L Cp4, N≤L wherein Cp3≤ L, M≤L Cp4≤ L, L are the length of original time domain sequences, and M is the numerical value that rounds up divided by the sampling interval the maximum sampling lag deviation time, and N is that maximum sampling shifts to an earlier date the numerical value that the deviation time rounded up divided by the sampling interval; With back L in the original time domain sequences Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, forms first time domain sequences.
Accordingly, receiver module 20 is from sampling initial time sampled data, and forms second time domain sequences according to the order of sampling; L from the second time domain sequences time domain sequences Cp3L sampling number certificate of intercepting lighted in individual sampling, and forms calibrating sequence according to the order of intercepting, and wherein Ts is the sampling interval.Perhaps
Receiver module 20 is delayed L from the sampling initial time Cp3* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences according to the order of sampling, and wherein Ts is the sampling interval; With second time domain sequences as calibrating sequence.Perhaps
Receiver module 20 is from sampling initial time sampled data, and forms second time domain sequences according to the order of sampling; N from second time domain sequences 3L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted, wherein N≤N according to the order of intercepting 3≤ L Cp3+ L Cp4-M and N 3≠ L Cp3If N 3<L Cp3, then receiver side moves L to sequence big-endian circulation to be adjusted Cp3-N 3The position obtains calibrating sequence; If N 3>L Cp3, then receiver side moves N to sequence little-endian circulation to be adjusted 3-L Cp3The position obtains calibrating sequence.
Preferable, compensate if carry out 20 pairs of channel estimation values of cyclic shift receiver module.If N 3<L Cp3, N in the mode three of the method that specifically compensates and Fig. 2 situation three 3<L Cp3The method that compensates is identical, does not repeat them here.If N 3>L Cp3, N in the mode three of the method that then specifically compensates and Fig. 2 situation three 3>L Cp3The method that compensates is identical, does not repeat them here.
If what carry out is to send calibration, then sending module 10 is baseband signal processors, and receiver module 20 is calibrating signal processors; If what carry out is to receive calibration, then sending module 10 is calibrating signal processors, and receiver module 20 is baseband signal processors.
In concrete enforcement, the system of the antenna calibration of the embodiment of the invention can also be placed the base station, namely sending module 10 and receiver module 20 place the base station.
Those skilled in the art should understand that embodiments of the invention can be provided as method, system or computer program.Therefore, the present invention can adopt complete hardware embodiment, complete software embodiment or in conjunction with the form of the embodiment of software and hardware aspect.And the present invention can adopt the form of the computer program of implementing in one or more computer-usable storage medium (including but not limited to magnetic disc store, CD-ROM, optical memory etc.) that wherein include computer usable program code.
The present invention is that reference is described according to flow chart and/or the block diagram of method, equipment (system) and the computer program of the embodiment of the invention.Should understand can be by the flow process in each flow process in computer program instructions realization flow figure and/or the block diagram and/or square frame and flow chart and/or the block diagram and/or the combination of square frame.Can provide these computer program instructions to the processor of all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing device to produce a machine, make the instruction of carrying out by the processor of computer or other programmable data processing device produce to be used for the device of the function that is implemented in flow process of flow chart or a plurality of flow process and/or square frame of block diagram or a plurality of square frame appointments.
These computer program instructions also can be stored in energy vectoring computer or the computer-readable memory of other programmable data processing device with ad hoc fashion work, make the instruction that is stored in this computer-readable memory produce the manufacture that comprises command device, this command device is implemented in the function of appointment in flow process of flow chart or a plurality of flow process and/or square frame of block diagram or a plurality of square frame.
These computer program instructions also can be loaded on computer or other programmable data processing device, make and carry out the sequence of operations step producing computer implemented processing at computer or other programmable devices, thereby be provided for being implemented in the step of the function of appointment in flow process of flow chart or a plurality of flow process and/or square frame of block diagram or a plurality of square frame in the instruction that computer or other programmable devices are carried out.
Although described the preferred embodiments of the present invention, in a single day those skilled in the art get the basic creative concept of cicada, then can make other change and modification to these embodiment.So claims are intended to all changes and the modification that are interpreted as comprising preferred embodiment and fall into the scope of the invention.
Because transmitter side adds cyclic sequence in original time domain sequences, can be in antenna calibration system, reduce owing to sampling deviation produces significantly influence to receiving baseband signal, thereby reduce the loss of calibrating sequence signal, improve the accurate and antenna calibration precision of calibration factor; Further improve performance and the DOA estimated accuracy of wave beam forming.
Obviously, those skilled in the art can carry out various changes and modification to the present invention and not break away from the spirit and scope of the present invention.Like this, if of the present invention these are revised and modification belongs within the scope of claim of the present invention and equivalent technologies thereof, then the present invention also is intended to comprise these changes and modification interior.

Claims (27)

1. the method for an antenna calibration is characterized in that, this method comprises:
Transmitter side adds cyclic sequence in original time domain sequences, obtain first time domain sequences;
Described transmitter side sends described first time domain sequences to receiver side;
Described receiver side is determined calibrating sequence according to second time domain sequences that sampling obtains;
Described receiver side carries out antenna calibration according to calibrating sequence;
Wherein, described transmitter side obtains first time domain sequences and comprises: described transmitter side determines that the length of cyclic sequence equals L Cp1, L wherein Cp1Be positive integer, M+N≤L Cp1≤ L, L are the length of original time domain sequences, and M is the numerical value that rounds up divided by the sampling interval the maximum sampling lag deviation time, and N is that maximum sampling shifts to an earlier date the numerical value that the deviation time rounded up divided by the sampling interval; Described transmitter side is with back L in the original time domain sequences Cp1The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence the foremost of original time domain sequences to, forms first time domain sequences; Perhaps
Described transmitter side determines that the length of cyclic sequence equals L Cp2, L wherein Cp2Be positive integer, M+N≤L Cp2≤ L, L are the length of original time domain sequences, and M is the numerical value that rounds up divided by the sampling interval the maximum sampling lag deviation time, and N is that maximum sampling shifts to an earlier date the numerical value that the deviation time rounded up divided by the sampling interval; Described transmitter side is with preceding L in the original time domain sequences Cp2The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence to original time domain sequences backmost, forms first time domain sequences; Perhaps
Described transmitter side determines that the length of first cyclic sequence equals L Cp3, and determine that the length of second cyclic sequence equals L Cp4, L wherein Cp3Be positive integer, L Cp4Be positive integer, N≤L Cp3≤ L, M≤L Cp4≤ L, L are the length of original time domain sequences, and M is the numerical value that rounds up divided by the sampling interval the maximum sampling lag deviation time, and N is that maximum sampling shifts to an earlier date the numerical value that the deviation time rounded up divided by the sampling interval; Described transmitter side is with back L in the original time domain sequences Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, forms first time domain sequences.
2. the method for claim 1 is characterized in that, when described transmitter side by L after in original time domain sequences Cp1The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence the foremost of original time domain sequences to, and when forming first time domain sequences, described receiver side determines that calibrating sequence comprises:
Described receiver side is from sampling initial time sampled data, and forms second time domain sequences according to the order of sampling;
Described receiver side is from the N of described second time domain sequences 1L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted, wherein N according to the order of intercepting 1Be positive integer, and N≤N 1≤ L Cp1-M;
Described receiver side moves L to sequence big-endian circulation to be adjusted Cp1-N 1The position obtains calibrating sequence.
3. the method for claim 1 is characterized in that, when described transmitter side with original time domain sequences in after L Cp1The time domain data of position is formed cyclic sequence in order, and by cyclic sequence being added to the foremost of original time domain sequences, when forming first time domain sequences, described receiver side determines that calibrating sequence comprises:
Described receiver side is delayed N from the sampling initial time 1* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences, wherein N according to the order of sampling 1Be positive integer, and N≤N 1≤ L Cp1-M, Ts are the sampling intervals;
Circulation moves L to described receiver side to the second time domain sequences big-endian Cp1-N 1The position obtains calibrating sequence.
4. as the arbitrary described method of claim 1~3, it is characterized in that, when described transmitter side with original time domain sequences in after L Cp1The time domain data of position is formed cyclic sequence in order, and by cyclic sequence being added to the foremost of original time domain sequences, when forming first time domain sequences, described receiver side carries out antenna calibration according to calibrating sequence and comprises:
Described receiver side is determined the frequency domain channel estimated value according to described calibrating sequence;
Described receiver side compensates according to following formula:
H 1=H×e j2πkn
Wherein, H 1Be the frequency domain channel estimated value after the compensation,
Figure FDA00003285795900021
K is the subcarrier sequence number, and H is k subcarrier channel estimation, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion;
Described receiver side carries out antenna calibration according to the channel estimation value after compensating.
5. the method for claim 1 is characterized in that, preceding L in described transmitter side passes through original time domain sequences Cp2The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver side determines that calibrating sequence comprises:
Described receiver side is from sampling initial time sampled data, and forms second time domain sequences according to the order of sampling;
Described receiver side is from the N of described second time domain sequences 2L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted, wherein N according to the order of intercepting 2Be positive integer, and N≤N 2≤ L Cp2-M;
Described receiver side moves N to sequence little-endian circulation to be adjusted 2The position obtains calibrating sequence.
6. the method for claim 1 is characterized in that, preceding L in described transmitter side passes through original time domain sequences Cp2The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver side obtains calibrating sequence and comprises:
Described receiver side is delayed N from the sampling initial time 2* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences, wherein N according to the order of sampling 2Be positive integer, and N≤N 2≤ L Cp2-M, Ts are the sampling intervals;
Circulation moves N to described receiver side to the second time domain sequences little-endian 2The position obtains calibrating sequence.
7. as claim 1,5 or 6 described methods, it is characterized in that preceding L in described transmitter side passes through original time domain sequences Cp2The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver side carries out antenna calibration according to calibrating sequence and comprises:
Described receiver side is determined the frequency domain channel estimated value according to described calibrating sequence;
Described receiver side compensates according to following formula:
H 2=H×e -j2πkn
Wherein, H 2Be the frequency domain channel estimated value after the compensation, K is the subcarrier sequence number, and H is the frequency domain channel estimated value of k subcarrier, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion;
Described receiver side carries out antenna calibration according to the channel estimation value after compensating.
8. the method for claim 1 is characterized in that, when described transmitter side with original time domain sequences in after L Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver side determines that calibrating sequence comprises:
Described receiver side is from sampling initial time sampled data, and forms second time domain sequences according to the order of sampling;
Described receiver side is from the L of described second time domain sequences Cp3L sampling number certificate of intercepting lighted in individual sampling, and forms calibrating sequence according to the order of intercepting.
9. the method for claim 1 is characterized in that, when described transmitter side with original time domain sequences in after L Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver side obtains calibrating sequence and comprises:
Described receiver side is delayed L from the sampling initial time Cp3* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences according to the order of sampling, and wherein Ts is the sampling interval;
Described receiver side with described second time domain sequences as calibrating sequence.
10. the method for claim 1 is characterized in that, when described transmitter side with original time domain sequences in after L Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver side obtains calibrating sequence and comprises:
Described receiver side is from sampling initial time sampled data, and forms second time domain sequences according to the order of sampling;
Described receiver side is from the N of second time domain sequences 3L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted, wherein N≤N according to the order of intercepting 3≤ L Cp3+ L Cp4-M, and N 3≠ L Cp3
If N 3<L Cp3, described receiver side moves L to sequence big-endian circulation to be adjusted Cp3-N 3The position obtains calibrating sequence;
If N 3L Cp3, described receiver side moves N to sequence little-endian circulation to be adjusted 3-L Cp3The position obtains calibrating sequence.
11. the method for claim 1 is characterized in that, when described transmitter side with original time domain sequences in after L Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver side obtains calibrating sequence and comprises:
Described receiver side is delayed N from the sampling initial time 3* Ts constantly rises and begins to receive time domain data, and forms second time domain sequences, wherein N≤N according to the order of sampling 3≤ L Cp3+ L Cp4-M, and N 3≠ L Cp3, Ts is the sampling interval;
If N 3<L Cp3, circulation moves L to described receiver side to the second time domain sequences big-endian Cp3-N 3The position obtains calibrating sequence;
If N 3L Cp3, described receiver side moves N to sequence little-endian circulation to be adjusted 3-L Cp3The position obtains calibrating sequence.
12., it is characterized in that described receiver side carries out antenna calibration according to calibrating sequence and comprises as claim 10 or 11 described methods:
Described receiver side is determined the frequency domain channel estimated value according to described calibrating sequence;
If N 3<L Cp3, described receiver side compensates according to following formula:
H 3=H×e j2πkn
Wherein, H 3Be the frequency domain channel estimated value after the compensation,
Figure FDA00003285795900051
K is the subcarrier sequence number, and H is k subcarrier channel estimation, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion;
If N 3L Cp3, described receiver side compensates according to following formula:
H 3=H×e -j2πkn
Wherein, H 3Be the frequency domain channel estimated value after the compensation, K is the subcarrier sequence number, and H is the frequency domain channel estimated value of k subcarrier, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion;
Described receiver side carries out antenna calibration according to the channel estimation value after compensating.
13., it is characterized in that if described transmitter side is baseband signal processor, described receiver side is the calibrating signal processor as claim 1~3,5~6,8~11 arbitrary described methods;
If described transmitter side is the calibrating signal processor, described receiver side is baseband signal processor.
14. the system of an antenna calibration is characterized in that, this system comprises:
Sending module is used for adding cyclic sequence in original time domain sequences and obtains first time domain sequences, sends described first time domain sequences;
Receiver module is used for determining calibrating sequence according to second time domain sequences that sampling obtains, and carries out antenna calibration according to calibrating sequence;
Wherein, described sending module specifically is used for: the length of determining cyclic sequence equals L Cp1, L wherein Cp1Be positive integer, M+N≤L Cp1≤ L, L are the length of original time domain sequences, and M is the numerical value that rounds up divided by the sampling interval the maximum sampling lag deviation time, and N is that maximum sampling shifts to an earlier date the numerical value that the deviation time rounded up divided by the sampling interval; With back L in the original time domain sequences Cp1The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence the foremost of original time domain sequences to, forms first time domain sequences; The length of perhaps determining cyclic sequence equals L Cp2, L wherein Cp2Be positive integer, M+N≤L Cp2≤ L, L are the length of original time domain sequences, and M is the numerical value that rounds up divided by the sampling interval the maximum sampling lag deviation time, and N is that maximum sampling shifts to an earlier date the numerical value that the deviation time rounded up divided by the sampling interval; With preceding L in the original time domain sequences Cp2The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence to original time domain sequences backmost, forms first time domain sequences; Perhaps
Described cyclic sequence comprises first cyclic sequence and second cyclic sequence; Described sending module specifically is used for: the length of determining first cyclic sequence equals L Cp3, and determine that the length of second cyclic sequence equals L Cp4, L wherein Cp3Be positive integer, L Cp4Be positive integer, N≤L Cp3≤ L, M≤L Cp4≤ L, L are the length of original time domain sequences, and M is the numerical value that rounds up divided by the sampling interval the maximum sampling lag deviation time, and N is that maximum sampling shifts to an earlier date the numerical value that the deviation time rounded up divided by the sampling interval; With back L in the original time domain sequences Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, forms first time domain sequences.
15. system as claimed in claim 14 is characterized in that, L after described sending module specifically is used for original time domain sequences Cp1The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence the foremost of original time domain sequences to, and when forming first time domain sequences, described receiver module specifically is used for:
From sampling initial time sampled data, and form second time domain sequences according to the order of sampling; N from described second time domain sequences 1L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted, wherein N according to the order of intercepting 1Be positive integer, and N≤N 1≤ L Cp1-M; Sequence big-endian circulation to be adjusted is moved L Cp1-N 1The position obtains calibrating sequence.
16. system as claimed in claim 14 is characterized in that, L after described sending module specifically is used for original time domain sequences Cp1The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence the foremost of original time domain sequences to, and when forming first time domain sequences, described receiver module specifically is used for:
Delay N from the sampling initial time 1* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences, wherein N according to the order of sampling 1Be positive integer, and N≤N 1≤ L Cp1-M, Ts are the sampling intervals; Described second time domain sequences big-endian circulation is moved L Cp1-N 1The position obtains calibrating sequence.
17., it is characterized in that L after described sending module specifically is used for original time domain sequences as the arbitrary described system of claim 14~16 Cp1The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence the foremost of original time domain sequences to, and when forming first time domain sequences, described receiver module specifically is used for:
Determine channel estimation value according to described calibrating sequence; Compensate according to following formula; Carry out antenna calibration according to the channel estimation value after the compensation;
H 1=H×e j2πkn
Wherein, H 1Be the channel estimation value after the compensation,
Figure FDA00003285795900081
K is the subcarrier sequence number, and H is k subcarrier channel estimation, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion.
18. system as claimed in claim 14 is characterized in that, L before described sending module specifically is used for original time domain sequences Cp2The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver module specifically is used for:
From sampling initial time sampled data, and form second time domain sequences according to the order of sampling; N from described second time domain sequences 2L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted, wherein N according to the order of intercepting 2Be positive integer, and N≤N 2≤ L Cp2-M; Sequence little-endian circulation to be adjusted is moved N 2The position obtains calibrating sequence.
19. system as claimed in claim 14 is characterized in that, L before described sending module specifically is used for original time domain sequences Cp2The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver module specifically is used for:
Delay N from the sampling initial time 2* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences, wherein N according to the order of sampling 2Be positive integer, and N≤N 2≤ L Cp2-M, Ts are the sampling intervals; Described second time domain sequences little-endian circulation is moved N 2The position obtains calibrating sequence.
20., it is characterized in that L before described sending module specifically is used for original time domain sequences as claim 14,18 or 19 described systems Cp2The time domain data of position is formed cyclic sequence in order, and adds cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver module specifically is used for:
Determine channel estimation value according to described calibrating sequence; Compensate according to following formula; Carry out antenna calibration according to the channel estimation value after the compensation;
H 2=H×e -j2πkn
Wherein, H 2Be the channel estimation value after the compensation,
Figure FDA00003285795900091
K is the subcarrier sequence number, and H is k subcarrier channel estimation, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion.
21. system as claimed in claim 14 is characterized in that, when described cyclic sequence comprises first cyclic sequence and second cyclic sequence, described sending module specifically is used for original time domain sequences back L Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver module specifically is used for:
From sampling initial time sampled data, and form second time domain sequences according to the order of sampling; L from the described second time domain sequences time domain sequences Cp3L sampling number certificate of intercepting lighted in individual sampling, and forms calibrating sequence according to the order of intercepting, and wherein Ts is the sampling interval.
22. system as claimed in claim 14 is characterized in that, when described cyclic sequence comprises first cyclic sequence and second cyclic sequence, described sending module specifically is used for original time domain sequences back L Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver module specifically is used for:
Delay L from the sampling initial time Cp3* Ts rises constantly and begins to sample the L point data, and forms second time domain sequences according to the order of sampling, and wherein Ts is the sampling interval; With described second time domain sequences as calibrating sequence.
23. system as claimed in claim 14 is characterized in that, when described cyclic sequence comprises first cyclic sequence and second cyclic sequence, described sending module specifically is used for original time domain sequences back L Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver module specifically is used for:
From sampling initial time sampled data, and form second time domain sequences according to the order of sampling; N from second time domain sequences 3L sampling number certificate of intercepting lighted in individual sampling, and forms sequence to be adjusted, wherein N≤N according to the order of intercepting 3≤ L Cp3+ L Cp4-M, and N 3≠ L Cp3If N 3<L Cp3, sequence big-endian circulation to be adjusted is moved L Cp3-N 3The position obtains calibrating sequence; If N 3L Cp3, sequence little-endian circulation to be adjusted is moved N 3-L Cp3The position obtains calibrating sequence.
24. system as claimed in claim 14 is characterized in that, when described cyclic sequence comprises first cyclic sequence and second cyclic sequence, described sending module specifically is used for original time domain sequences back L Cp3The time domain data of position is formed first cyclic sequence in order, and adds first cyclic sequence foremost of original time domain sequences to, and with L before in the original time domain sequences Cp4The time domain data of position is formed second cyclic sequence in order, and adds second cyclic sequence to original time domain sequences backmost, and when forming first time domain sequences, described receiver module specifically is used for:
Delay N from the sampling initial time 3* Ts constantly rises and begins to receive time domain data, and forms second time domain sequences, wherein N≤N according to the order of sampling 3≤ L Cp3+ L Cp4-M, and N 3≠ L Cp3, Ts is the sampling interval; If N 3<L Cp3, circulation moves L to described receiver side to the second time domain sequences big-endian Cp3-N 3The position obtains calibrating sequence; If N 3L Cp3, described receiver side moves N to sequence little-endian circulation to be adjusted 3-L Cp3The position obtains calibrating sequence.
25. as claim 23 or 24 described systems, it is characterized in that described receiver module specifically is used for:
Determine the frequency domain channel estimated value according to described calibrating sequence;
If N 3<L Cp3, compensate according to following formula:
H 3=H×e j2πkn
Wherein, H 3Be the frequency domain channel estimated value after the compensation,
Figure FDA00003285795900101
K is the subcarrier sequence number, and H is k subcarrier channel estimation, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion;
If N 3L Cp3, compensate according to following formula:
H 3=H×e -j2πkn
Wherein, H 3Be the frequency domain channel estimated value after the compensation,
Figure FDA00003285795900111
K is the subcarrier sequence number, and H is the frequency domain channel estimated value of k subcarrier, N FFTTotal sampling number when doing the FFT conversion when being time-frequency conversion;
Carry out antenna calibration according to the channel estimation value after the compensation.
26., it is characterized in that if described sending module is baseband signal processor, described receiver module is the calibrating signal processor as claim 14~16,18~19,21~24 arbitrary described systems;
If described sending module is the calibrating signal processor, described receiver module is baseband signal processor.
27. as claim 14~16,18~19,21~24 arbitrary described systems, it is characterized in that the system of described antenna calibration is in the base station.
CN 201110030019 2011-01-27 2011-01-27 Antenna calibration method and system Active CN102148636B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN 201110030019 CN102148636B (en) 2011-01-27 2011-01-27 Antenna calibration method and system
PCT/CN2012/070518 WO2012100699A1 (en) 2011-01-27 2012-01-18 Method and system for antennas calibration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110030019 CN102148636B (en) 2011-01-27 2011-01-27 Antenna calibration method and system

Publications (2)

Publication Number Publication Date
CN102148636A CN102148636A (en) 2011-08-10
CN102148636B true CN102148636B (en) 2013-09-04

Family

ID=44422667

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110030019 Active CN102148636B (en) 2011-01-27 2011-01-27 Antenna calibration method and system

Country Status (2)

Country Link
CN (1) CN102148636B (en)
WO (1) WO2012100699A1 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102148636B (en) * 2011-01-27 2013-09-04 大唐移动通信设备有限公司 Antenna calibration method and system
WO2015085510A1 (en) * 2013-12-11 2015-06-18 Telefonaktiebolaget L M Ericsson (Publ) Methods and apparatus for antenna calibration
CN104023344B (en) * 2014-05-21 2017-12-22 北京创毅视讯科技有限公司 A kind of beam form-endowing method and user equipment

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1446000A (en) * 2002-12-25 2003-10-01 大唐移动通信设备有限公司 Method for adjusting intelligences antenna array system in real time
KR20060073257A (en) * 2004-12-24 2006-06-28 삼성전자주식회사 Fast Fourier Transform Start Position Correction Device Using Cyclic Prefix Length Estimation in Frequency Space Block Code-Orthogonal Frequency Division Multiplexing System
CN101064902A (en) * 2006-04-25 2007-10-31 大唐移动通信设备有限公司 Method for real-time calibrating intelligent antenna
CN101188448A (en) * 2006-11-15 2008-05-28 大唐移动通信设备有限公司 A smart antenna calibration method, device and system
JP2008311800A (en) * 2007-06-13 2008-12-25 Hitachi Communication Technologies Ltd OFDM cellular radio communication base station
CN101854323A (en) * 2009-04-03 2010-10-06 中兴通讯股份有限公司 Method and system for calibrating antenna
CN101958856A (en) * 2009-07-13 2011-01-26 联芯科技有限公司 Method and device for carrying out channel estimation by using cyclic prefix code

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007031691A1 (en) * 2005-09-14 2007-03-22 France Telecom Device and method for configuring a receiver in a communication device
KR101200736B1 (en) * 2006-03-14 2012-11-13 퀄컴 인코포레이티드 Apparatus and method for transmitting and receiving packet data in a smart antenna system
CN101783722B (en) * 2009-01-21 2013-09-18 普天信息技术研究院有限公司 Transmission method and device for virtual MIMO
CN102148636B (en) * 2011-01-27 2013-09-04 大唐移动通信设备有限公司 Antenna calibration method and system

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1446000A (en) * 2002-12-25 2003-10-01 大唐移动通信设备有限公司 Method for adjusting intelligences antenna array system in real time
KR20060073257A (en) * 2004-12-24 2006-06-28 삼성전자주식회사 Fast Fourier Transform Start Position Correction Device Using Cyclic Prefix Length Estimation in Frequency Space Block Code-Orthogonal Frequency Division Multiplexing System
CN101064902A (en) * 2006-04-25 2007-10-31 大唐移动通信设备有限公司 Method for real-time calibrating intelligent antenna
CN101188448A (en) * 2006-11-15 2008-05-28 大唐移动通信设备有限公司 A smart antenna calibration method, device and system
JP2008311800A (en) * 2007-06-13 2008-12-25 Hitachi Communication Technologies Ltd OFDM cellular radio communication base station
CN101854323A (en) * 2009-04-03 2010-10-06 中兴通讯股份有限公司 Method and system for calibrating antenna
CN101958856A (en) * 2009-07-13 2011-01-26 联芯科技有限公司 Method and device for carrying out channel estimation by using cyclic prefix code

Also Published As

Publication number Publication date
CN102148636A (en) 2011-08-10
WO2012100699A1 (en) 2012-08-02

Similar Documents

Publication Publication Date Title
KR102036210B1 (en) Method and apparatus for searching cell in beamforming system
CN110301123B (en) Coherent detection of large physical random access control channel (PRACH) delays
CN101753185B (en) A kind of methods, devices and systems realizing multi-cell multi-antenna calibration
CN103873396A (en) Frequency offset estimation and compensation method and device
JP2020504584A (en) Method and apparatus for transmitting and receiving uplink signals in a wireless communication system
US10178639B2 (en) Method for controlling timing of terminal in wireless communication system, and electronic device therefor
US10931482B2 (en) Method and apparatus for channel estimation in wireless communication system
JP2021536154A (en) Methods implemented on network equipment and methods implemented on terminal equipment
EP2901576B1 (en) Apparatus, system and method of communicating filter information
CN104105984A (en) Method and apparatus for estimating speed of a mobile terminal
CN102148636B (en) Antenna calibration method and system
CN103209010A (en) Antenna calibration method and base band unit
US20090232051A1 (en) Method and system for the extension of frequency offset estimation range based on correlation of complex sequences
US20210219258A1 (en) Methods and apparatuses for communicating position reference signals
EP2159927B1 (en) Method and system for the extension of frequency offset range estimation based on correlation of complex sequences
US10090893B2 (en) Frequency division multiple antenna distribution
US11683710B2 (en) Method and apparatus for interference measurement in wireless communication system
CN102783108A (en) Method and apparatus for estimating a Doppler frequency
CN103139884A (en) Method and device of enhancing calibration performance of intelligent antenna of test driven development (TDD) system
CN105656535A (en) Method and device for downlink calibration in Massive MIMO (Multiple Input Multiple Output) system
CN101321008A (en) Descending beam forming emission method and device
CN102421184A (en) Method and equipment for determining time advance
EP3711200B1 (en) Method and device for positioning communication device
CN102665269B (en) Symbol timing synchronization method and symbol timing synchronization system
US20230375655A1 (en) Dynamic antenna array angular phase deviation compensation for broad beam positioning

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant