CN105306127B - Mobile communication timing estimation method and device - Google Patents
Mobile communication timing estimation method and device Download PDFInfo
- Publication number
- CN105306127B CN105306127B CN201410759147.9A CN201410759147A CN105306127B CN 105306127 B CN105306127 B CN 105306127B CN 201410759147 A CN201410759147 A CN 201410759147A CN 105306127 B CN105306127 B CN 105306127B
- Authority
- CN
- China
- Prior art keywords
- sequence
- main peak
- value
- mobile communication
- position corresponding
- 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
Links
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
The present invention provides a kind of mobile communication timing estimation method and devices, wherein method includes: the leader sequence for receiving mobile terminal and sending, and the leader sequence is to touch ZC sequence certainly according to first and be conjugated the first ZC retrieval；Local ZC sequence and the leader sequence are subjected to related operation, obtain the value of position corresponding to the first main peak；The conjugation of local ZC sequence and the leader sequence are subjected to related operation, obtain the value of position corresponding to the second main peak；According to the value of position corresponding to the value of position corresponding to first main peak and second main peak, the carrier wave time delay in mobile communication is obtained.The mobile communication timing estimation method and device provided through the invention can accurately measure the time delay of carrier wave, meet the needs of mobile satellite communication.
Description
Technical field
The present invention relates to network communication technology field more particularly to a kind of mobile communication timing estimation methods and device.
Background technique
Satellite communication is important one of the means of communication of today's society.Satellite communication area coverage is big, changes traditional
Communication mode can easily provide global transmission for multiple business, while satellite colleague can overcome terrestrial coverage area
The shortcomings that finiteness, and communication distance is remote, is capable of providing the advantages that transmission of high capacity, in current all kinds of new traffic business
To application and quick development.
In mobile communication system, due to relative motion between satellite and mobile terminal, it will lead to what receiving end received
Wave frequency shifts, and this phenomenon is known as Doppler effect.It is how general for the satellite mobile communication system of stationary orbit
Strangle the movement that Producing reason is mainly user terminal；For non-geo satellite, Doppler frequency shift depends primarily on satellite
Quick movement relative to terrestrial user.Due to non-geo satellite orbit altitude, classification of track, earth station latitude and defending
Situations such as position of the star area of coverage, is different, and the Doppler frequency shift effect thus caused is also different.When earth station sees satellite from ground
Plane rises or has a maximum Doppler frequency shift when disappearing, and when surface of the satellite in mobile terminal, Doppler frequency shift is
Zero.
In long term evolution (Long Term Evolution, abbreviation LTE) satellite mobile communication system, if connect at random
There are biggish Doppler frequency shifts during entering, then will affect from the correlation for touching (Zadoff-Chu, abbreviation ZC) sequence, thus
The performance of signal detection is seriously affected.In the random access procedure of LTE satellite mobile communication system, accidental access signal is used
Be ZC sequence, for ZC sequence, Doppler frequency shift will lead to it and lose zero auto-correlation, seriously affect the detection of signal
Performance.In ground system, downlink user equipment (User Equipment, abbreviation UE) carry out frequency offset correction after, uplink with
Machine accesses between usually not more than one times of frequency deviation of the subcarrier of channel (Random Access Channel, abbreviation RACH)
Every the position of main peak will not shift, and can observe apparent peak value in detection window or in neighbouring several windows.
But in satellite mobile communication system, Doppler frequency shift may will offset out completely inspection considerably beyond subcarrier spacing, main peak
Window is surveyed, the detection on ground under existing high velocity environment is unable to satisfy the communicating requirement on satellite.
Summary of the invention
For the defects in the prior art, the present invention provides a kind of mobile communication timing estimation method and device, Neng Gouzhun
The true time delay for measuring carrier wave, meets the needs of mobile satellite communication.
In a first aspect, the present invention provides a kind of mobile communication timing estimation method, comprising:
The leader sequence that mobile terminal is sent is received, the leader sequence is to touch ZC sequence and conjugation first certainly according to first
ZC retrieval；
Local ZC sequence and the leader sequence are subjected to related operation, obtain the value of position corresponding to the first main peak；
The conjugation of local ZC sequence and the leader sequence are subjected to related operation, obtain position corresponding to the second main peak
Value；
According to the value of position corresponding to the value of position corresponding to first main peak and second main peak, moved
Carrier wave time delay in dynamic communication.
Further, the leader sequence are as follows:
Wherein, ZCu,NIt (n) is the first ZC sequence,For the conjugation of the first ZC sequence, N is the length of ZC sequence, u
For the radical exponent of ZC sequence.
Further, described that local ZC sequence and the leader sequence are subjected to related operation, it is right to obtain the first main peak institute
The value for the position answered, comprising:
Local ZC sequence and the leader sequence are subjected to related operation, obtain the first main peak；
According to first main peak, the value of position corresponding to first main peak is determined；
First main peak are as follows:
Wherein, Δ fk,d=fk,d/fRA,fRAFor subcarrier spacing, fk,dFor frequency deviation, hkIt is rung for the channel impulse of k-th of user
It answers, τkFor the propagation delay of k-th of user, N is the length of ZC sequence, and u is the radical exponent of ZC sequence, and m is the time delay of ZC sequence.
Further, the value of position corresponding to first main peak are as follows:
d1=τk+Δfk,du-1+aN
Wherein, τkFor the propagation delay of k-th of user, Δ fk,dFor frequency deviation, N is the length of ZC sequence, and a is integer, u-1It is full
Foot ((uu-1))N=1.
Further, the conjugation by local ZC sequence and the leader sequence carry out related operation, and it is main to obtain second
The value of position corresponding to peak, comprising:
The conjugation of local ZC sequence and the leader sequence are subjected to related operation, obtain the second main peak；
According to second main peak, the value of position corresponding to second main peak is determined；
Second main peak are as follows:
Wherein, Δ fk,d=fk,d/fRA,fRAFor subcarrier spacing, fk,dFor frequency deviation, hkIt is rung for the channel impulse of k-th of user
It answers, τkFor the propagation delay of k-th of user, N is the length of ZC sequence, and u is the radical exponent of ZC sequence, and m is the time delay of ZC sequence.
Further, the value of position corresponding to second main peak are as follows:
d2=τk-Δfk,du-1+bN
Wherein, τkFor the propagation delay of k-th of user, Δ fk,dFor frequency deviation, N is the length of ZC sequence, and b is integer, u-1It is full
Foot ((uu-1))N=1.
Further, position corresponding to the value of the position according to corresponding to first main peak and second main peak
The value set obtains the carrier wave time delay in mobile communication, comprising:
According to the value of position corresponding to the value of position corresponding to first main peak and second main peak, use
First formula obtains the carrier wave time delay in mobile communication；
First formula are as follows:
Wherein, c=a+b, c are integer, 0 < τk<N,τkFor integer.
Further, which is characterized in that the method also includes:
If the carrier wave time delay in the mobile communication is unsatisfactory for preset subcarrier spacing integral multiple, to described leading
Sequence carries out frequency compensation, obtains compensated leader sequence；
Local ZC sequence and the compensated leader sequence are subjected to related operation, obtain position corresponding to the first main peak
The value set；
The conjugation of local ZC sequence and the compensated leader sequence are subjected to related operation, it is right to obtain the second main peak institute
The value for the position answered；
According to the value of position corresponding to the value of position corresponding to first main peak and second main peak, moved
Carrier wave time delay in dynamic communication.
It is further, described that leader sequence is compensated are as follows:
Wherein, k is positive integer.
Second aspect, the present invention provide a kind of mobile communication timing estimating apparatus, comprising:
Receiving unit, for receiving the leader sequence of mobile terminal transmission, the leader sequence is to touch ZC certainly according to first
The first ZC retrieval of sequence and conjugation；
First main peak unit obtains the first main peak for local ZC sequence and the leader sequence to be carried out related operation
The value of corresponding position；
Second main peak unit obtains for the conjugation of local ZC sequence and the leader sequence to be carried out related operation
The value of position corresponding to two main peaks；
Carrier wave time delay elements, for corresponding to the value of the position according to corresponding to first main peak and second main peak
Position value, obtain mobile communication in carrier wave time delay.
As shown from the above technical solution, a kind of mobile communication timing estimation method and device provided through the invention,
In, method includes: the leader sequence for receiving mobile terminal and sending, and the leader sequence is to touch ZC sequence and conjugation certainly according to first
First ZC retrieval；Local ZC sequence and the leader sequence are subjected to related operation, obtained corresponding to the first main peak
The value of position；The conjugation of local ZC sequence and the leader sequence are subjected to related operation, obtain position corresponding to the second main peak
The value set；According to the value of position corresponding to the value of position corresponding to first main peak and second main peak, moved
Carrier wave time delay in dynamic communication.The mobile communication timing estimation method and device provided through the invention, can accurately measure
The time delay of carrier wave, meets the needs of mobile satellite communication.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below to embodiment or existing
Attached drawing needed in technical description is briefly described, it should be apparent that, the accompanying drawings in the following description is only this hair
Bright some embodiments for those of ordinary skill in the art without creative efforts, can be with root
Other attached drawings are obtained according to these attached drawings.
Fig. 1 is a kind of flow diagram of mobile communication timing estimation method provided in an embodiment of the present invention；
Fig. 2 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=1 |
corrkk1| peak value schematic diagram；
Fig. 3 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=1.3 |
corrkk1| peak value schematic diagram；
Fig. 4 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=1.7 |
corrkk1| peak value schematic diagram；
Fig. 5 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=2 |
corrkk1| peak value schematic diagram；
Fig. 6 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=1 |
corrkk2| peak value schematic diagram；
Fig. 7 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=1.3 |
corrkk2| peak value schematic diagram；
Fig. 8 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=1.7 |
corrkk2| peak value schematic diagram；
Fig. 9 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=2 |
corrkk2| peak value schematic diagram；
In the case that Figure 10 is a kind of zero offset provided in an embodiment of the present invention | corrkk1| peak condition schematic diagram；
In the case that Figure 11 is a kind of zero offset provided in an embodiment of the present invention | corrkk2| peak condition schematic diagram；
Figure 12 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=10.45 |
corrkk1| peak Distribution schematic diagram；
Figure 13 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=10.45+ δ
| corrkk1| peak Distribution schematic diagram；
Figure 14 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=10.45+ δ
| corrkk2| peak Distribution schematic diagram；
Figure 15 is a kind of structural schematic diagram of mobile communication timing estimating apparatus provided in an embodiment of the present invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, the technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is only
It is only a part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiment of the present invention, ordinary skill people
Member's every other embodiment obtained without making creative work, shall fall within the protection scope of the present invention.
Fig. 1 is a kind of flow diagram of mobile communication timing estimation method provided in an embodiment of the present invention, such as Fig. 1 institute
Show, the mobile communication timing estimation method of the present embodiment is as described below.
101, the leader sequence that mobile terminal is sent is received, the leader sequence is to touch ZC sequence together certainly according to first
The first ZC retrieval of yoke.
It should be understood that receiving the leader sequence that mobile terminal is sent, leader sequence is the first ZC sequence and the first ZC sequence
The sum of the conjugation of column.
ZC sequence has good correlation, and in the signal of random access, generallys use ZC sequence.
Leader sequence are as follows:
Wherein, ZCu,NIt (n) is the first ZC sequence,For the conjugation of the first ZC sequence, N is the length of ZC sequence, u
For the radical exponent of ZC sequence.
102, local ZC sequence and the leader sequence are subjected to related operation, obtain position corresponding to the first main peak
Value.
It should be understood that the leader sequence received and local ZC sequence are carried out related operation, the first main peak is obtained
The value of corresponding position.
It is understood that local ZC sequence and the leader sequence are carried out related operation, the first main peak is obtained,
According to first main peak, the value of position corresponding to first main peak is determined.
First main peak are as follows:
Wherein, Δ fk,d=fk,d/fRA,fRAFor subcarrier spacing, fk,dFor frequency deviation, hkIt is rung for the channel impulse of k-th of user
It answers, τkFor the propagation delay of k-th of user, N is the length of ZC sequence, and u is the radical exponent of ZC sequence, and m is the time delay of ZC sequence.
The value of position corresponding to first main peak are as follows:
d1=τk+Δfk,du-1+aN
Wherein, τkFor the propagation delay of k-th of user, Δ fk,dFor frequency deviation, N is the length of ZC sequence, and a is integer, u-1It is full
Foot ((uu-1))N=1.
103, the conjugation of local ZC sequence and the leader sequence are subjected to related operation, obtained corresponding to the second main peak
Position value.
It should be understood that the leader sequence received and local ZC sequence are carried out related operation, the second main peak institute is obtained
The value of corresponding position.
The conjugation of local ZC sequence and the leader sequence are subjected to related operation, obtain the second main peak；
According to second main peak, the value of position corresponding to second main peak is determined；
Second main peak are as follows:
Wherein, Δ fk,d=fk,d/fRA,fRAFor subcarrier spacing, fk,dFor frequency deviation, hkIt is rung for the channel impulse of k-th of user
It answers, τkFor the propagation delay of k-th of user, N is the length of ZC sequence, and u is the radical exponent of ZC sequence, and m is the time delay of ZC sequence.
The value of position corresponding to second main peak are as follows:
d2=τk-Δfk,du-1+bN
Wherein, τkFor the propagation delay of k-th of user, Δ fk,dFor frequency deviation, N is the length of ZC sequence, and b is integer, u-1It is full
Foot ((uu-1))N=1.
104, the value of position corresponding to the value of the position according to corresponding to first main peak and second main peak, is obtained
Obtain the carrier wave time delay in mobile communication.
It is understood that the value of position corresponding to the value of the position according to corresponding to the first main peak and the second main peak,
Obtain the carrier wave time delay in mobile communication.
According to the value of position corresponding to the value of position corresponding to first main peak and second main peak, using
One formula obtains the carrier wave time delay in mobile communication；
First formula are as follows:
Wherein, c=a+b, c are integer, 0 < τk<N,τkFor integer.
If the carrier wave time delay in the mobile communication is unsatisfactory for preset subcarrier spacing integral multiple, to described leading
Sequence carries out frequency compensation, obtains compensated leader sequence；
Local ZC sequence and the compensated leader sequence are subjected to related operation, obtain position corresponding to the first main peak
The value set；
The conjugation of local ZC sequence and the compensated leader sequence are subjected to related operation, it is right to obtain the second main peak institute
The value for the position answered；
According to the value of position corresponding to the value of position corresponding to first main peak and second main peak, moved
Carrier wave time delay in dynamic communication.
Leader sequence is compensated are as follows:
Wherein, k is positive integer.
The mobile communication timing estimation method provided through this embodiment can accurately measure the time delay of carrier wave, meet
The demand of mobile satellite communication.
Fig. 2 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=1 |
corrkk1| peak value schematic diagram, Fig. 3 is in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,d
When=1.3 | corrkk1| peak value schematic diagram, Fig. 4 be a kind of mobile communication timing estimation method provided in an embodiment of the present invention
In as Δ fk,dWhen=1.7 | corrkk1| peak value schematic diagram, Fig. 5 is that a kind of mobile communication provided in an embodiment of the present invention is fixed
When estimation method in as Δ fk,dWhen=2 | corrkk1| peak value schematic diagram, Fig. 6 be a kind of shifting provided in an embodiment of the present invention
As Δ f in dynamic communication timing estimation methodk,dWhen=1 | corrkk2| peak value schematic diagram, Fig. 7 provides for the embodiment of the present invention
A kind of mobile communication timing estimation method in as Δ fk,dWhen=1.3 | corrkk2| peak value schematic diagram, Fig. 8 be the present invention
As Δ f in a kind of mobile communication timing estimation method that embodiment providesk,dWhen=1.7 | corrkk2| peak value schematic diagram, figure
9 be in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=2 | corrkk2| peak value show
It is intended to, in the case that Figure 10 is a kind of zero offset provided in an embodiment of the present invention | corrkk1| peak condition schematic diagram, Figure 11
In the case where a kind of zero offset provided in an embodiment of the present invention | corrkk2| peak condition schematic diagram, Figure 12 is that the present invention is real
It applies in a kind of mobile communication timing estimation method of example offer as Δ fk,dWhen=10.45 | corrkk1| peak Distribution signal
Figure, Figure 13 are in a kind of mobile communication timing estimation method provided in an embodiment of the present invention as Δ fk,dWhen=10.45+ δ |
corrkk1| peak Distribution schematic diagram, Figure 14 be in a kind of mobile communication timing estimation method provided in an embodiment of the present invention when
Δfk,dWhen=10.45+ δ | corrkk2| peak Distribution schematic diagram, as shown in Fig. 2-14, the mobile communication of the present embodiment is fixed
When estimation method it is as described below.
For ZC sequence, Doppler frequency shift will lead to it and lose zero auto-correlation, to seriously affect the inspection of signal
Survey performance.In the presence of frequency deviation Δ f, ZC sequence can generate the distortion such as following formula:
Wherein, TseqIt is Physical Random Access Channel (Physical Random Access Channel, abbreviation PRACH)
The sampling interval of leader sequence.With the frequency deviation Δ f=± Δ f of PRACH subcarrier same sizeRA=± 1/TseqLead to ZC sequence
xμ(n) d onμ=(± 1/ μ) modNZCCyclic shift.It is worth noting that, | Δ f |=Δ fRAPlace, leader sequence peak value exist
Desired position completely disappears.
It is as follows to establish system model,
Wherein, h is channel impulse response, and τ is Systematic Communication time delay, fdFor Doppler frequency shift, μ, ν are that the root of ZC sequence refers to
Number, w (t) are system noise function, and k, l represent different accessing users.Above system model is subjected to discretization, can be obtained:
Wherein, Δ fk,d=fk,d/fRA, fRAFor subcarrier spacing.
New random access leader sequence, i.e., it is real by the sum of original ZC sequence and ZC sequence conjugate as leader sequence
The leader sequence of the random access at current family, design is as follows:
It is utilized respectively ZC sequence and ZC conjugate sequence when signal detection and makees related operation with the leader sequence received, then will
The correlation summation being calculated, to realize the signal detection and timing estimation in unknown Doppler frequency shift.At present
The case where only considering single user, specifically:
(1) ZC is utilizedu,N(n) related operation is carried out with the leader sequence received.
Wherein corrkkIndicate use and ZCu,N(n) user and ZC of sequence same root uu,N(n) correlation, corrklTable
Show use and ZCu,N(n) user and ZC of sequence difference root uu,N(n) correlation,Indicate noise and ZCu,N(n) sequence
Correlation.According to the characteristic of ZC sequence, the correlation of different root sequences is constantIt withWhen only influencing to detect
Background noise level, therefore selective analysis corr herekkValue.
corrkk1Indicate the part ZC and local ZC sequence correlation, corr in received leader sequencekk2It indicates to receive sequence
The correlation of middle ZC conjugate moiety and local ZC sequence.Since ZC is conjugated from the cross-correlation of ZC sequence similar to different roots and ZC sequence
Column cross-correlation only influences background noise level, therefore, selective analysis corrkk1。
As (um- (u τk+Δfk,d))N=0, | corrkk1| there is maximum value；
As (um- (u τk+Δfk,d))N=q, and q be integer when, | corrkk1|=0；
As (um- (u τk+Δfk,d))N=q, and q be decimal when, | corrkk1| there is relationship with q；
When Fig. 2 to Fig. 5 shows that frequency deviation becomes larger from 1 times of subcarrier spacing to 2 times of subcarrier spacings, signal and ZC sequence are received
The changing rule of peak value after column are related.When frequency deviation is m times of subcarrier spacing of integer, peak value is located at mdμmodNZCPlace, wherein
dμ=(± 1/ μ) modNZC, μ=4, N in this simulation processZC=839, so dμ=210.When frequency deviation from 1 times of subcarrier spacing to
When 2 times of subcarrier spacing variations, peak energy is gradually from dμIt is transferred to 2dμ.Also, when frequency deviation is decimal times subcarrier spacing,
Energy will disperse dμIntegral multiple point at.
(2) it usesRelated operation is carried out with the leader sequence received.
As (- um+ (u τk-Δfk,d))N=0, | corrkk2| there is maximum value；
As (- um+ (u τk-Δfk,d))N=q, and q be integer when, | corrkk2|=0；
As (- um+ (u τk-Δfk,d))N=q, and q be decimal when, | corrkk2| there is relationship with q；
When Fig. 6 to Fig. 9 shows that frequency deviation becomes larger from 1 times of subcarrier spacing to 2 times of subcarrier spacings, signal and ZC sequence are received
The changing rule of peak value after the conjugation of column is related.As it can be seen that the changing rule for receiving signal and the correlation peak of ZC sequence conjugate is same
It receives signal and the variation of ZC sequence correlation peak is symmetrical, the raw effect of conjugation fecund of the identical frequency deviation to ZC and ZC is opposite.
(3) in detection and timing estimation, it is assumed that frequency deviation is the integral multiple of subcarrier spacing at this time
If | corrkk1| and | corrkk2| the corresponding position of main peak is d1And d2, then
d1=τk+Δfk,du-1+aN
d2=τk-Δfk,du-1+bN
Wherein u-1Meet ((uu-1))N=1；A and b is integer；
d1+d2=2 τk+ (a+b) N, i.e. d1+d2=2 τk+ cN, c are integer；
0<τk< N, and τkFor integer, can acquire:
In conjunction with Fig. 2 and Fig. 6, d1=210, d2=629,
In conjunction with Fig. 5 and Fig. 9, d1=420, d2=419,
Therefore, it is right in the case where Doppler frequency shift is the integral multiple of subcarrier spacing finally to be realized by this algorithm
The detection and timing estimation of leader sequence.
It is minimum to be likely to occur the case where frequency deviation is subcarrier spacing integral multiple but in practical communication environment.Therefore, will
To how to be carried out in the case of small several times detection and timing estimation carry out analysis verifying.By analysis above, can be realized
Doppler frequency shift be subcarrier spacing integral multiple in the case where detection and timing estimation to leader sequence, then can lead to
The mode for crossing frequency drift compensation makes decimal frequency multiplication shift to integral multiple frequency displacement to approach, and then converts integral multiple for small several times problem approximation
Problem finally realizes the detection in the case of small several times.
Assuming that leader sequence Cν=150, system Doppler frequency shift fk,d=Δ fk,d·fRA=10.45fRA。
As shown in Figure 10,11, zero-frequency move in the case where, the position of correlation peak will not shift, peak coordinate with
CνIt is equal.Due to | corrkk1| with | corrkk2| influenced to be symmetrical by frequency displacement, therefore we only consider for the time being | corrkk1|
With the situation of change of frequency drift compensation.Under the conditions of Gaussian channel, when frequency displacement is fk,dWhen, peak Distribution is as shown in figure 12:
It can be found that there are two more apparent peak values in Figure 12, the two peak value sizes are close, which can not judge
It is a be real peak value which be pseudo- peak.Moreover the two peak values are smaller, certain false-alarm probability condition may be unable to satisfy
The requirement of lower predetermined threshold value, so energy is transferred to another from one of them by way of frequency drift compensation, Ke Yitong
It is compensated after several timesMode by Δ fk,dStep wise approximation integer.
As δ=+ 0.5, | corrkk1| as shown in figure 13:
After frequency drift compensation, become a peak value of figure l3 from two peak values in Figure 12, peak value also accordingly becomes larger.
It is then possible to further decrease | δ | to improve the precision of compensation.At this point, d1=782.
Figure 14 is after overcompensation δ | corrkk2| peak Distribution situation.At this point, d2=357.
Therefore, in the case where realizing the small several times that Doppler frequency shift is subcarrier spacing by frequency compensated mode
Detection and timing estimation to leader sequence.
In conclusion passing through while sending the conjugation of ZC sequence Yu ZC sequence, and local ZC is respectively adopted in receiving end
Sequence and its conjugation are associated, to realize the signal detection and timing estimation under maximum Doppler frequency shift.Pass through two groups of correlations
Symmetric relation under identical Doppler frequency shift offsets influence of the Doppler frequency shift to timing estimation, using the progressive compensation of frequency displacement
Mode by the energy transfer to main peak at pseudo- peak, and then determine the actual position of main peak, realize the detection to leader sequence
And timing estimation.
Figure 15 is a kind of structural schematic diagram of mobile communication timing estimating apparatus provided in an embodiment of the present invention, such as Figure 15
Shown, the mobile communication timing estimating apparatus of the present embodiment is as described below.
Mobile communication timing estimating apparatus includes: receiving unit 151, the first main peak unit 152, the second main peak unit 153
With carrier wave time delay elements 154.
Receiving unit 151, for receive mobile terminal transmission leader sequence, the leader sequence be according to first from
It touches ZC sequence and is conjugated the first ZC retrieval.
First main peak unit 152 obtains first for local ZC sequence and the leader sequence to be carried out related operation
The value of position corresponding to main peak.
Second main peak unit 153 is obtained for the conjugation of local ZC sequence and the leader sequence to be carried out related operation
Obtain the value of position corresponding to the second main peak.
Carrier wave time delay elements 154, value and the second main peak institute for the position according to corresponding to first main peak
The value of corresponding position obtains the carrier wave time delay in mobile communication.
The mobile communication timing estimating apparatus provided through this embodiment can accurately measure the time delay of carrier wave, meet
The demand of mobile satellite communication.
Those of ordinary skill in the art will appreciate that: realize that all or part of step of above method embodiment can lead to
The relevant hardware of program instruction is crossed to complete, program above-mentioned can store in computer-readable storage medium, the journey
Sequence when being executed, executes step including the steps of the foregoing method embodiments；And storage medium above-mentioned include: ROM, RAM, magnetic disk or
In the various media that can store program code such as CD.
The above description is merely a specific embodiment, and still, protection scope of the present invention is not limited to this, appoints
What those familiar with the art in the technical scope disclosed by the present invention, the variation or substitution that can be readily occurred in, all
It is covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with the scope of protection of the claims
Subject to.
Claims (10)
1. a kind of mobile communication timing estimation method characterized by comprising
The leader sequence that mobile terminal is sent is received, the leader sequence is to touch the first ZC sequence of ZC sequence and conjugation certainly according to first
What column obtained；
Local ZC sequence and the leader sequence are subjected to related operation, obtain the value of position corresponding to the first main peak；
The conjugation of local ZC sequence and the leader sequence are subjected to related operation, obtain position corresponding to the second main peak
Value；
According to the value of position corresponding to the value of position corresponding to first main peak and second main peak, obtain mobile logical
Carrier wave time delay in letter.
2. mobile communication timing estimation method according to claim 1, which is characterized in that the leader sequence are as follows:
Wherein, ZCu,NIt (n) is the first ZC sequence,For the conjugation of the first ZC sequence, N is the length of ZC sequence, u ZC
The radical exponent of sequence.
3. mobile communication timing estimation method according to claim 1, which is characterized in that described by local ZC sequence and institute
It states leader sequence and carries out related operation, obtain the value of position corresponding to the first main peak, comprising:
Local ZC sequence and the leader sequence are subjected to related operation, obtain the first main peak；
According to first main peak, the value of position corresponding to first main peak is determined；
First main peak are as follows:
Wherein, Δ fk,d=fk,d/fRA,fRAFor subcarrier spacing, fk,dFor frequency deviation, hkFor the channel impulse response of k-th of user,
τkFor the propagation delay of k-th of user, N is the length of ZC sequence, and u is the radical exponent of ZC sequence, and m is the time delay of ZC sequence.
4. according to claim 1 or mobile communication timing estimation method described in 3, which is characterized in that the first main peak institute
The value of corresponding position are as follows:
d1=τk+Δfk,du-1+aN
Wherein, τkFor the propagation delay of k-th of user, Δ fk,dFor frequency deviation, N is the length of ZC sequence, and a is integer, u-1Meet
((uu-1))N=1.
5. mobile communication timing estimation method according to claim 1, which is characterized in that being total to local ZC sequence
Yoke and the leader sequence carry out related operation, obtain the value of position corresponding to the second main peak, comprising:
The conjugation of local ZC sequence and the leader sequence are subjected to related operation, obtain the second main peak；
According to second main peak, the value of position corresponding to second main peak is determined；
Second main peak are as follows:
Wherein, Δ fk,d=fk,d/fRA,fRAFor subcarrier spacing, fk,dFor frequency deviation, hkFor the channel impulse response of k-th of user,
τkFor the propagation delay of k-th of user, N is the length of ZC sequence, and u is the radical exponent of ZC sequence, and m is the time delay of ZC sequence.
6. according to claim 1 or mobile communication timing estimation method described in 5, which is characterized in that the second main peak institute
The value of corresponding position are as follows:
d2=τk-Δfk,du-1+bN
Wherein, τkFor the propagation delay of k-th of user, Δ fk,dFor frequency deviation, N is the length of ZC sequence, and b is integer, u-1Meet
((uu-1))N=1.
7. mobile communication timing estimation method according to claim 1, which is characterized in that described according to first main peak
The value of position corresponding to the value of corresponding position and second main peak obtains the carrier wave time delay in mobile communication, comprising:
According to the value of position corresponding to the value of position corresponding to first main peak and second main peak, using the first public affairs
Formula obtains the carrier wave time delay in mobile communication；
First formula are as follows:
Wherein, c=a+b, c are integer, 0 < τk<N,τkFor integer.
8. mobile communication timing estimation method according to claim 1, which is characterized in that the method also includes:
If the carrier wave time delay in the mobile communication is unsatisfactory for preset subcarrier spacing integral multiple, to the leader sequence
Frequency compensation is carried out, compensated leader sequence is obtained；
Local ZC sequence and the compensated leader sequence are subjected to related operation, obtain position corresponding to the first main peak
Value；
The conjugation of local ZC sequence and the compensated leader sequence are subjected to related operation, obtained corresponding to the second main peak
The value of position；
According to the value of position corresponding to the value of position corresponding to first main peak and second main peak, obtain mobile logical
Carrier wave time delay in letter.
9. mobile communication timing estimation method according to claim 8, which is characterized in that it is described to the leader sequence into
Line frequency compensation are as follows:
Wherein, k is positive integer.
10. a kind of mobile communication timing estimating apparatus characterized by comprising
Receiving unit, for receiving the leader sequence of mobile terminal transmission, the leader sequence is to touch ZC sequence certainly according to first
With the first ZC retrieval of conjugation；
It is right to obtain the first main peak institute for local ZC sequence and the leader sequence to be carried out related operation for first main peak unit
The value for the position answered；
It is main to obtain second for the conjugation of local ZC sequence and the leader sequence to be carried out related operation for second main peak unit
The value of position corresponding to peak；
Carrier wave time delay elements, for position corresponding to the value of the position according to corresponding to first main peak and second main peak
The value set obtains the carrier wave time delay in mobile communication.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410759147.9A CN105306127B (en) | 2014-12-10 | 2014-12-10 | Mobile communication timing estimation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410759147.9A CN105306127B (en) | 2014-12-10 | 2014-12-10 | Mobile communication timing estimation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105306127A CN105306127A (en) | 2016-02-03 |
| CN105306127B true CN105306127B (en) | 2019-01-25 |
Family
ID=55202948
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410759147.9A Active CN105306127B (en) | 2014-12-10 | 2014-12-10 | Mobile communication timing estimation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105306127B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106253944A (en) * | 2016-09-13 | 2016-12-21 | 北京捷联微芯科技有限公司 | The method of information transmission and node |
| CN110673120B (en) * | 2019-09-18 | 2022-03-29 | 四川天邑康和通信股份有限公司 | Moving rate determination method based on Zadoff-Chu sequence |
| CN111083080B (en) * | 2019-12-26 | 2022-12-30 | 北京华力创通科技股份有限公司 | Method and device for synchronizing PRACH (physical random Access channel) of broadband satellite communication system |
| CN114079606B (en) * | 2020-08-17 | 2023-10-24 | 海能达通信股份有限公司 | Air interface time alignment method and device and electronic equipment |
| CN112887241B (en) | 2021-04-29 | 2021-07-16 | 成都星联芯通科技有限公司 | Frequency offset estimation method, device, communication equipment and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102104951A (en) * | 2009-12-17 | 2011-06-22 | 中兴通讯股份有限公司 | Method and device for timing shift estimation |
| MY148707A (en) * | 2007-05-02 | 2013-05-31 | Qualcomm Inc | Selection of root indices in polyphase cazac sequences |
| JP2014127897A (en) * | 2012-12-27 | 2014-07-07 | National Institute Of Information & Communication Technology | Data transmitter and data transmission method |
| CN103929825A (en) * | 2014-04-30 | 2014-07-16 | 电子科技大学 | Multi-user detection method based on ZC sequence |
| CN104038461A (en) * | 2013-03-04 | 2014-09-10 | 联想(北京)有限公司 | Symbol synchronization and channel estimation method, apparatus and system |
-
2014
- 2014-12-10 CN CN201410759147.9A patent/CN105306127B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| MY148707A (en) * | 2007-05-02 | 2013-05-31 | Qualcomm Inc | Selection of root indices in polyphase cazac sequences |
| CN102104951A (en) * | 2009-12-17 | 2011-06-22 | 中兴通讯股份有限公司 | Method and device for timing shift estimation |
| JP2014127897A (en) * | 2012-12-27 | 2014-07-07 | National Institute Of Information & Communication Technology | Data transmitter and data transmission method |
| CN104038461A (en) * | 2013-03-04 | 2014-09-10 | 联想(北京)有限公司 | Symbol synchronization and channel estimation method, apparatus and system |
| CN103929825A (en) * | 2014-04-30 | 2014-07-16 | 电子科技大学 | Multi-user detection method based on ZC sequence |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105306127A (en) | 2016-02-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105306127B (en) | Mobile communication timing estimation method and device | |
| US8295311B2 (en) | Detection of time-frequency hopping patterns | |
| KR100804904B1 (en) | Base station identification in orthongonal frequency division multiplexing based spread spectrum multiple access systems | |
| CN103098379B (en) | Interference is eliminated with the method estimating UTDOA in detection reference signal | |
| CN101295999A (en) | Detection method for accidental precursor access | |
| CN101489238B (en) | Time difference measuring method, system and apparatus | |
| CN102298097B (en) | Method for estimating thunder impulse signal Time Difference of Arrival (TDOA) | |
| WO2014063275A1 (en) | Method for determining remote same-frequency interference source and locating method therefor | |
| EP3554138A1 (en) | Method for determining candidate line-of-sight path, and wireless local area network device | |
| CN104022996A (en) | Channel estimation-based timing synchronization method for orthogonal frequency division multiplexing (OFDM) system | |
| CN101635598B (en) | Method and device for estimating noise power | |
| CN105531600B (en) | Time analysis in wireless network for user velocity estimation | |
| CN105072063A (en) | Method and device for suppressing interference signal | |
| US8923376B2 (en) | Apparatus and method for estimating noise and interference power in communication system | |
| CN104779993B (en) | Deep space TT＆C system and method based on frequency domain equalization | |
| CN105306097B (en) | A kind of detection method of accidental access signal, device and system | |
| CN106817707B (en) | Method and device for detecting and assisting in detecting signal source in base station | |
| CN108495365A (en) | Method of locating terminal based on the estimation of narrowband Internet of Things time delay | |
| RU2582590C1 (en) | Method of estimating frequency shift for communication systems using ofdm signals | |
| CN104243373B (en) | A kind of asynchronous interference detection method and device reached under scene of multi-user | |
| Leria et al. | Timing-based location estimation for OFDM signals with application in LTE, WLAN and WIMAX | |
| Haghparast et al. | OFDM ranging in multipath channels using time reversal method | |
| CN114640410B (en) | Wireless channel type detection method and device | |
| Nausner et al. | Positioning with 5G reference signals for vehicular applications | |
| Shuaieb et al. | Radio-location techniques under adverse channel conditions |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |