CN101868022A - Cluster-based distributed antenna time frequency synchronizing method in wireless communication - Google Patents

Cluster-based distributed antenna time frequency synchronizing method in wireless communication Download PDF

Info

Publication number
CN101868022A
CN101868022A CN 201010185688 CN201010185688A CN101868022A CN 101868022 A CN101868022 A CN 101868022A CN 201010185688 CN201010185688 CN 201010185688 CN 201010185688 A CN201010185688 A CN 201010185688A CN 101868022 A CN101868022 A CN 101868022A
Authority
CN
China
Prior art keywords
alpha
bunch
delta
spaced antenna
phi
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.)
Granted
Application number
CN 201010185688
Other languages
Chinese (zh)
Other versions
CN101868022B (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN2010101856887A priority Critical patent/CN101868022B/en
Publication of CN101868022A publication Critical patent/CN101868022A/en
Application granted granted Critical
Publication of CN101868022B publication Critical patent/CN101868022B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radio Transmission System (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

The invention relates to a cluster-based distributed antenna time frequency synchronizing method in wireless communication, and belongs to the technical field of wireless communication. The method comprises the following steps of: dividing mass distributed antennas into clusters, selecting cluster representative antennas, and performing distributed bidirectional coordinated synchronization on the cluster representative antennas; and then, initiating intra-cluster distributed bidirectional coordinated synchronization by the cluster representative antennas so as to finally realize time frequency synchronization of all distributed antennas in a system. The method has the following advantages that: through cluster division, a plurality of clusters are simultaneously synchronized, the total synchronization time is reduced and the synchronization speed is quickened; the synchronization between the cluster representative antennas is performed first and then the intra-cluster distributed bidirectional coordinated synchronization is performed so as to obtain higher synchronization precision and meet the requirement of a distributed wireless communication system for synchronization precision; and the synchronization between the distributed antennas is finished in advance so as to avoid multiple times of synchronization between a terminal user and a plurality of distributed antennas. The method reduces the system resource overhead of the terminal user, and improves the communication reliability.

Description

In a kind of radio communication based on bunch distributed antenna time frequency synchronizing method
Technical field
The present invention relates in a kind of radio communication based on bunch distributed antenna time frequency synchronizing method, belong to wireless communication technology field.
Background technology
In the next generation wireless communication system of identical networking,, adopt the collaborative mode in spaced antenna or many base stations to transmit for the availability of frequency spectrum of further raising system with to quality of services for users more.For example, coordinate multipoint (CoordinatedMultiple Point, hereinafter to be referred as CoMP) transmission technology is the technology of being paid close attention to by LTE-Advanced at present, it can avoid or reduce presence of intercell interference effectively on the basis that keeps system high-frequency spectrum utilance as far as possible.The networking model that is based on BBU+RRU that the CoMP model is the most frequently used, common 1 Base Band Unit (hereinafter to be referred as BBU) connects a plurality of wireless remotes unit (hereinafter to be referred as RRU) by radio frequency optical fiber (hereinafter to be referred as RoF).RRU is similar to spaced antenna in the distributing antenna system, and the BBU by RoF and similar central processing unit is connected.
Under the scene of a plurality of spaced antenna collaboration communications, time between a plurality of spaced antennas and Frequency Synchronization provide the basic premise of reliable transmission.When a plurality of spaced antennas transmit information for one or more users simultaneously,, make user's restoring signal difficulty because the frequency of the local crystal oscillator of each spaced antenna can produce deviation in time.Thereby, guarantee that the temporal frequency of a plurality of spaced antennas is synchronous, be directly connected to user's's (as cellphone subscriber, mobile device etc.) communication quality and reliability.Synchronous motor at many spaced antennas is shaped on principal and subordinate's synchronization mechanism, opens and closes ring synchronization mechanism and two-way poll synchronization mechanism etc. at present.
In principal and subordinate's synchronization mechanism, host node (a selected representative antennas or terminal use in many spaced antennas) sends reference signal to each from node (remaining spaced antenna), receive frequency and the phase information of obtaining reference signal after the reference signal by phase-locked loop (PLL) from node, and then adjust from the local crystal oscillator of node and reach synchronous purpose.But should mechanism also come with some shortcomings part:
(1) host node needs to broadcast from node to a plurality of simultaneously, than bringing more power consumption from node.When especially being host node, then need to handle a large amount of next response signals since node with terminal use.
(2) owing to each need be adjusted the local crystal oscillator of oneself according to the reference signal that host node sends from node, in case the crystal oscillator deviation of host node is excessive, occur when perhaps each is from the node processing reference signal than mistake, can cause then that main and subordinate node can't the deadline frequency synchronously.Synchronization accuracy under principal and subordinate's synchronization mechanism also can't satisfy the temporal frequency synchronisation requirement of spaced antenna cooperative system in addition.The temporal frequency synchronization accuracy that current spaced antenna/multi-base station cooperative system needs is higher than about one to two order of magnitude of principal and subordinate's synchronization mechanism precision.
In two-way poll synchronization mechanism, consider a labelled spaced antenna set { C 1, C 2..., C L, and setting antenna number, to be incremented to L from 1 be forward direction, opposite from L be decremented to 1 for the back to.Send reference signal by base station 1 to base station 2, this signal comprises the frequency and the phase information of base station 1 local crystal oscillator, and base station 2 is estimated after receiving reference signal, and adjusted local crystal oscillator according to estimated value.Continue subsequently to send reference signal, by that analogy to base station 3.Each base station is carried out the estimation of time and frequency by the mode of two-way poll.Then,, calculate time and frequency values that the local crystal oscillator of each base station need be adjusted, finally reach the purpose of a plurality of base stations time, Frequency Synchronization by the symmetry of front and back bi-directional estimation.Two-way poll synchronization mechanism can be met time, the Frequency Synchronization precision of spaced antenna/multi-base station cooperative system requirements, but also there is weak point in this synchronization mechanism:
Frequency and phase estimation are carried out according to receiving reference signal in each base station, and continue to send reference signal and all need independently time slot to the process of next base station, have L spaced antenna cooperative multicast system for one, along with the increase of spaced antenna number L, the synchronously required resource overhead of whole system can linearly increase.In DWCS/multi-base station cooperative communication system, general spaced antenna number L under a central controller controls is bigger, and synchronizing process need periodically be carried out.Though two-way poll synchronization mechanism can access the higher synchronous precision, lock in time, expense was excessive.
In a word, present existing synchronization mechanism has problems aspect system resource overhead, lock in time and synchronization accuracy, can not realize high efficiency, high precision time Frequency Synchronization to a plurality of spaced antennas in DWCS/multi-base station cooperative communication system.
Summary of the invention
The objective of the invention is to propose in a kind of radio communication based on bunch distributed antenna time frequency synchronizing method, at first, the communication system that a plurality of spaced antennas are arranged is carried out sub-clustering handle, select a bunch representative antennas, finish the high-precise synchronization between bunch representative antennas.Then, the temporal frequency of being initiated a plurality of spaced antennas in each bunch by bunch representative antennas is synchronous, and the temporal frequency of all spaced antennas is synchronous in the finally realization system.
In the radio communication that the present invention proposes based on bunch distributed antenna time frequency synchronizing method, may further comprise the steps:
(1) spaced antenna of the M=PN in the wireless communication system is divided into P bunch, p=1,2 ..., P, the spaced antenna number in each bunch are N, n=1, and 2 ..., N, n spaced antenna usefulness in p bunch
Figure GSA00000141617300021
Expression, then the spaced antenna set is in every bunch
Figure GSA00000141617300022
P=1,2 ..., P;
Setting up departments, the unified reference clock signal of spaced antenna is ψ (t)=exp{j (Ω t+ Φ) in the system }, wherein Ω is the frequency of local crystal oscillator, and t is the unified overall reference time, and Φ is the phase place of local crystal oscillator.Owing to there is the deviation of local crystal oscillator, thus in the system any one spaced antenna i local crystal oscillator signal indication be:
ψ i(t)=exp{j(Ω it+Φ i+n i(t))},i=1,2,...,M-1,M (1)
Ω in the following formula iBe local crystal oscillator frequency, Φ iThe phase place of representing local crystal oscillator, n i(t) be the noise jamming of local crystal oscillator;
According to above-mentioned (1) formula, obtain that the current local zone time of any one spaced antenna i is in the said system:
t i=(Ω it+Φ i+n i(t))/Ω=α i(t+Δ i(t)),i=1,2,...,M-1,M (2)
α wherein ii/ Ω, the ratio of the local crystal oscillator frequency of any one spaced antenna i and unified overall crystal oscillator frequency in the expression system; Δ i(t)=(Φ i+ n i(t))/Ω i, the local zone time deviation of any one spaced antenna i in the expression system considers that the local clock deviation is a fixed value in a unit propagates time slot, then the local zone time t of any one spaced antenna i in the system iWith the overall situation reference time t relational expression be t i=t+ Δ i
(2) from N the spaced antenna of above-mentioned each bunch, choose 2 spaced antennas
Figure GSA00000141617300031
With As a bunch representative antennas, then adjacent P bunch produces 2P bunch representative antennas;
(3) above-mentioned all bunch representative antennas are carried out the distributed bidirectional cooperation synchronously, detailed process is as follows:
(3-1) above-mentioned all bunch representative antennas are numbered the back ordering, form a bunch representative antennas set
Figure GSA00000141617300033
Be abbreviated as { A 1, A 2..., A 2P-1, A 2P;
(3-2) by a bunch representative antennas A 1To a bunch representative antennas A 2Send sinusoidal signal, A by that analogy 1→ A 2→ ... → A 2P-1→ A 2P, be designated as fl transmission; After the forward direction end of transmission, by a bunch representative antennas A 2PSend identical sinusoidal signal to bunch representative antennas A 2P-1, A by that analogy 2P→ A 2P-1→ ... → A 2→ A 1, being designated as reverse transfer, fl transmission and reverse transfer take 4P-2 unit altogether and propagate time slot, wherein time slot { T (1), T (2)..., T (2P-2), T (2P-1)Be that the fl transmission time slot is gathered time slot { T (2P), T (2P+1)..., T (4P-3), T (4P-2)Be that the reverse transfer time slot is gathered;
(3-3) at above-mentioned fl transmission time slot T (1)In, bunch representative antennas A 1The signal indication that sends is:
S 1 (1)(t)=exp{j(w 1t 11)} (3)
W wherein 1Be a bunch representative antennas A 1Local crystal oscillator frequency, Φ 1Phase place for local crystal oscillator;
Bunch representative antennas A 2Reception is from a bunch representative antennas A 1Signal, this signal is with bunch representative antennas A 2Local zone time t 2Be expressed as:
r 2 ( 1 ) ( t 2 ) = H 1,2 exp { j ( α 1 w 1 ( t 1 / α 2 + Δ 1 - Δ 2 ) + Φ 1 ) } + ω 2 ( 1 ) ( t 2 ) - - - ( 4 )
H wherein 1,2For from a bunch representative antennas A 1To bunch representative antennas A 2Channel Transmission gain, and set fl transmission and the back is symmetry equivalent to the channel transmitted transmission gain,
Figure GSA00000141617300035
Be a bunch representative antennas A 2At T (1)In the noise that receives, then bunch representative antennas A 2Estimate that the crystal oscillator frequency and the phase place that obtain are:
w ^ 2 ( 1 ) = ( α 1 w 1 + w ~ 2 ( 1 ) ) / α 2 - - - ( 5 )
Φ ^ 2 ( 1 ) = α 1 w 1 ( Δ 1 - Δ 2 ) + Φ 1 + Φ ~ 2 ( 1 ) - - - ( 6 )
Wherein
Figure GSA00000141617300043
With
Figure GSA00000141617300044
Be estimated frequency error and phase estimation error;
(3-4) repeat above-mentioned steps (3-2)-(3-3), obtain a bunch representative antennas A iFl transmission crystal oscillator frequency and the phase place estimated be:
w ^ i = α i - 1 w ^ i - 1 / α i = α 1 w 1 / α i - - - ( 7 )
Φ ^ i = α 1 w 1 ( Δ i - 1 - Δ i ) + Φ ^ i - 1 = α 1 w 1 ( Δ 1 - Δ i ) + Φ 1 ; - - - ( 8 )
(3-5) reverse transfer is by a bunch representative antennas A 2PBeginning is up at reverse transfer time slot T (4p-2)Intra-cluster representative antennas A 1Receive till the final signal;
(3-6) at above-mentioned reverse transfer time slot T (2P)In, bunch representative antennas A 2PThe signal indication that sends is:
s 2P (2P)(t)=exp{j(w 2Pt 2P2P)} (9)
W wherein 2PBe a bunch representative antennas A 2PLocal crystal oscillator frequency, Φ 2PPhase place for local crystal oscillator;
Bunch representative antennas A 2P-1Reception is from a bunch representative antennas A 2PSignal, this signal is with bunch representative antennas A 2P-1Local zone time t 2p-1Be expressed as:
r 2 P - 1 ( 2 P ) ( t 2 P - 1 ) = H 2 P , 2 P - 1 exp { j ( α 2 P w 2 P ( t 2 P - 1 / α 2 P - 1 + Δ 2 P - Δ 2 P - 1 ) + Φ 2 P ) } + ω 2 P - 1 ( 2 P ) ( t 2 P - 1 ) - - - ( 10 )
H wherein 2P, 2P-1For from a bunch representative antennas A 2PTo bunch representative antennas A 2P-1Transmission channel gain,
Figure GSA00000141617300048
Be a bunch representative antennas A 2P-1At T (2P)In the noise that receives, then bunch representative antennas A 2P-1Estimate that the crystal oscillator frequency and the phase place that obtain are:
w ^ 2 P - 1 ( 2 P ) = ( α 2 P w 2 P + w ~ 2 P - 1 ( 2 P ) ) / α 2 P - 1 - - - ( 11 )
Φ ^ 2 P - 1 ( 2 P ) = α 2 P w 2 P ( Δ 2 P - Δ 2 P - 1 ) + Φ 2 P + Φ ~ 2 P - 1 ( 2 P ) - - - ( 12 )
Wherein
Figure GSA000001416173000411
With
Figure GSA000001416173000412
Be estimated frequency error and phase estimation error;
(3-7) repeat above-mentioned steps (3-5)-(3-6), obtain a bunch representative antennas A iReverse transfer crystal oscillator frequency and the phase place estimated be:
w ^ i = α i + 1 w ^ i + 1 / α i = α 2 P w 2 P / α i - - - ( 13 )
Φ ^ i = α 2 P w 2 P ( Δ i + 1 - Δ i ) + Φ ^ i + 1 = α 2 P w 2 P ( Δ 2 P - Δ i ) + Φ 2 P ; - - - ( 14 )
(3-8) according to the fl transmission estimated result and the reverse transfer estimated result of above-mentioned each that obtains bunch representative antennas, bunch representative antennas A then iFinal local crystal oscillator frequency and the phase place estimated are:
w ^ i = ( α 1 w 1 + α 2 P w 2 P ) / α i - - - ( 15 )
Φ ^ i = α 1 w 1 ( Δ 1 - Δ i ) + α 2 P w 2 P ( Δ 2 P - Δ i ) + Φ 1 + Φ 2 P - - - ( 16 )
And obtain bunch representative antennas A based on unified overall reference time t iLocal crystal oscillator signal be:
ψ i(t)=exp{j((α 1w 12Pw 2P)t+λ 12P)},i=1,2,...,2P-1,2P (17)
λ wherein nnw nΔ n+ Φ n
(3-9) according to above-mentioned bunch of representative antennas A iLocal crystal oscillator frequency and the estimated value of phase place, obtain the set of bunch representative antennas
Figure GSA00000141617300054
The synchronized result of local crystal oscillator signal be:
ψ n ( p ) ( t ) = exp { j ( ( α 1 ( 1 ) w 1 ( 1 ) + α N ( P ) w N ( P ) ) t + λ 1 ( P ) + λ N ( P ) ) } , p=1,2,...,P,n=1,2,...,N; (18)
(4) according to the synchronized result of bunch representative antennas between above-mentioned each bunch, then, distributed bidirectional cooperation between a plurality of spaced antennas in being initiated bunch by bunch representative antennas of each bunch synchronously, the same in default wherein and method for synchronous and the above-mentioned steps (3) is example with p bunch N spaced antenna, p=1,2, ..., P, concrete steps are as follows:
(4-1) after finishing the set synchronously of bunch representative antennas, to bunch in N spaced antenna be numbered the back ordering, form spaced antenna and gather
Figure GSA00000141617300056
Be abbreviated as { B 1, B 2..., B N, spaced antenna wherein
Figure GSA00000141617300057
Be two bunches of representative antennas of this bunch;
(4-2) by bunch in spaced antenna B 1To spaced antenna B 2Send sinusoidal signal, B by that analogy 1→ B 2→ ... → B N-1→ B N, be designated as fl transmission; After the forward direction end of transmission, by spaced antenna B NSend identical sinusoidal signal to spaced antenna B N-1, B by that analogy N→ B N-1→ ... → B 2→ B 1, being designated as reverse transfer, fl transmission and reverse transfer take 2N-2 unit altogether and propagate time slot, wherein time slot
Figure GSA00000141617300058
Be the set of fl transmission time slot, time slot
Figure GSA00000141617300059
Be the set of reverse transfer time slot;
(4-3) at above-mentioned fl transmission time slot
Figure GSA000001416173000510
In, spaced antenna B 1The signal indication that sends is:
x 1 (1)(t)=exp{j(w 1t 11)} (19)
W wherein 1Be spaced antenna B 1Local crystal oscillator frequency, Φ 1Phase place for local crystal oscillator;
Then, spaced antenna B 2Reception is from spaced antenna B 1Signal, this signal spaced antenna B 2Local zone time t 2Be expressed as:
y 2 ( 1 ) ( t 2 ) = H 1,2 exp { j ( α 1 w 1 ( t 2 / α 2 + Δ 1 - Δ 2 ) + Φ 1 ) } + ω 2 ( 1 ) ( t 2 ) - - - ( 20 )
H wherein 1,2For from spaced antenna B 1To spaced antenna B 2Channel Transmission gain, and set fl transmission and the back is symmetry equivalent to the channel transmitted transmission gain,
Figure GSA00000141617300062
Be spaced antenna B 2 In the noise that receives, then spaced antenna B 2Estimate that the crystal oscillator frequency and the phase place that obtain are:
w ^ 2 ( 1 ) = ( α 1 w 1 + w ~ 2 ( 1 ) ) / α 2 - - - ( 21 )
Φ ^ 2 ( 1 ) = α 1 w 1 ( Δ 1 - Δ 2 ) + Φ 1 + Φ ~ 2 ( 1 ) - - - ( 22 )
Wherein
Figure GSA00000141617300066
With
Figure GSA00000141617300067
Be estimated frequency error and phase estimation error;
(4-4) repeat above-mentioned steps (4-2)-(4-3), spaced antenna B in the spaced antenna set in obtaining bunch iFl transmission crystal oscillator frequency and the phase place estimated be:
w ^ i = α i - 1 w ^ i - 1 / α i = α 1 w 1 / α i - - - ( 23 )
Φ ^ i = α 1 w 1 ( Δ i - 1 - Δ i ) + Φ ^ i - 1 = α 1 w 1 ( Δ 1 - Δ i ) + Φ 1 ; - - - ( 24 )
(4-5) reverse transfer is by spaced antenna B NInitiate, up to time slot Middle spaced antenna B 1Receive till the final signal;
(4-6) at above-mentioned reverse transfer time slot
Figure GSA000001416173000611
In, spaced antenna B NThe signal indication that sends is:
x N (N)(t)=exp{j(w Nt NN)} (25)
W wherein NBe spaced antenna B NLocal crystal oscillator frequency, Φ NPhase place for local crystal oscillator;
Then, spaced antenna B N-1Reception is from spaced antenna B NSignal, this signal spaced antenna B N-1Local zone time t N-1Be expressed as:
y N - 1 ( N ) ( t N - 1 ) = H N , N - 1 exp { j ( α N w N ( t N - 1 / α N - 1 + Δ N - Δ N - 1 ) + Φ N ) } + ω N - 1 ( N ) ( t N - 1 ) - - - ( 26 )
H wherein N, N-1For from spaced antenna B NTo spaced antenna B N-1Transmission channel gain,
Figure GSA00000141617300071
Be spaced antenna B N-1
Figure GSA00000141617300072
In the noise that receives, then spaced antenna B N-1Estimate that the frequency and the phase place that obtain are
w ^ N - 1 ( N ) = ( α N w N + w ~ N - 1 ( N ) ) / α N - 1 - - - ( 27 )
Φ ^ N - 1 ( N ) = α N w N ( Δ N - Δ N - 1 ) + Φ N + Φ N - 1 ( N ) - - - ( 28 )
Wherein
Figure GSA00000141617300075
With
Figure GSA00000141617300076
Be estimated frequency error and phase estimation error;
(4-7) repeat above-mentioned steps (4-5)-(4-6), spaced antenna B in the spaced antenna set in obtaining bunch iReverse transfer crystal oscillator frequency and the phase place estimated be:
w ^ i = α i + 1 w ^ i + 1 / α i = α N w N / α i - - - ( 29 )
Φ ^ i = α N w N ( Δ i + 1 - Δ i ) + Φ ^ i + 1 = α N w N ( Δ N - Δ i ) + Φ N ; - - - ( 30 )
(4-8) according to above-mentioned obtain bunch in the fl transmission estimated result and the reverse transfer estimated result of spaced antenna, then spaced antenna B iFinal local crystal oscillator frequency and the phase place estimated are:
w ^ i = ( α 1 w 1 + α N w N ) / α i - - - ( 31 )
Φ ^ i = α 1 w 1 ( Δ 1 - Δ i ) + α N w N ( Δ N - Δ i ) + Φ 1 + Φ N - - - ( 32 )
And obtain spaced antenna B based on unified overall reference time t iLocal crystal oscillator signal be:
ψ i(t)=exp{j((α 1w 1Nw N)t+λ 1N)},i=1,2,...,N-1,N (33)
λ wherein nnw nΔ n+ Φ n
(4-9) according to spaced antenna B in above-mentioned bunch iLocal crystal oscillator frequency and the estimated value of phase place, obtain p bunch bunch in the spaced antenna set
Figure GSA000001416173000711
The synchronized result of local crystal oscillator signal be:
ψ n ( p ) ( t ) = exp { j ( ( α 1 ( p ) w 1 ( p ) + α N ( p ) w N ( p ) ) t + λ 1 ( p ) + λ N ( p ) ) } , p=1,2,...,P,n=1,2,...,N; (34)
(4-10) because p bunch bunch representative antennas
Figure GSA000001416173000713
With a bunch representative antennas In above-mentioned (3) step, finish the synchronous of local crystal oscillator, so the P in the whole system bunch is total to M spaced antenna deadline Frequency Synchronization with bunch representative antennas of other bunches.Be about to obtain in above-mentioned formula (18) the substitution formula (34) that the final synchronized result of all spaced antennas is in the whole system:
ψ n ( p ) ( t ) = exp { j [ 2 ( α 1 ( 1 ) w 1 ( 1 ) + α N ( P ) w N ( P ) ) t + 2 ( λ 1 ( P ) + λ N ( P ) ) ] } , p=1,2,...,P,n=1,2,...,N。(35)
In the radio communication that the present invention proposes based on bunch distributed antenna time frequency synchronizing method, its advantage is:
1, time frequency synchronizing method of the present invention by with the sub-district sub-clustering, makes a plurality of bunches to carry out simultaneously under the prerequisite that guarantees higher synchronization accuracy, having reduced total lock in time synchronously, has accelerated synchronizing speed.
If total M=PN of the spaced antenna number in the whole system, a basic transmission time slot in the synchronizing process is 1.Directly the two-way poll method for synchronous of utilization realizes that the time overhead of global synchronization is T Tot=2M-2 is along with the linear growth trend of spaced antenna number M in the whole system.And expense lock in time of utilization institute of the present invention extracting method is
Figure GSA00000141617300081
Promptly the spaced antenna number when each bunch of choose reasonable is
Figure GSA00000141617300082
The time, required lock in time, expense was only with the linear growth trend of square root of the spaced antenna number M in the whole system.
2, time frequency synchronizing method of the present invention, it is synchronous earlier bunch representative antennas of each bunch to be carried out temporal frequency, and then the cooperation of the distributed bidirectional in carrying out bunch can obtain higher synchronous accuracy synchronously, can satisfy the requirement of DWCS/wireless communication systems such as multi-base station cooperative wireless communication system to synchronization accuracy.
3, time frequency synchronizing method of the present invention, utilization based on bunch distributed antenna time frequency synchronizing method, can make a plurality of spaced antennas with before terminal use's (as cellphone subscriber, mobile device etc.) communicates by letter, the temporal frequency of finishing in advance between the spaced antenna is synchronous, avoid the terminal use to carry out repeatedly synchronous with a plurality of spaced antennas respectively, and only need to finish synchronously with a spaced antenna, just can reach synchronous with the temporal frequency of the spaced antenna of all participation collaboration communications.Reduce terminal use's system resource overhead, improved the reliability of communication.
Description of drawings
Fig. 1 is the FB(flow block) of the inventive method.
Fig. 2 is a bunch division model schematic diagram in one embodiment of the present of invention.
Fig. 3 is the distributed bidirectional of bunch representative antennas set among Fig. 2 synchronous schematic diagram of cooperating.
Embodiment
Below in conjunction with Fig. 1 introduce in the radio communication that the present invention proposes based on bunch an embodiment of distributed antenna time frequency synchronizing method.The present invention sends out among the embodiment, and the central controller of wireless communication system connects 21 spaced antennas, and concrete steps are as follows:
(1) M=21 spaced antenna is divided into bunch, bunch uses the p mark, total P=3 bunch, spaced antenna n mark in bunch has N=7 spaced antenna in each bunch.N antenna in p bunch used Expression, then the set of the spaced antenna of each bunch is
Figure GSA00000141617300092
P=1,2,3, Fig. 2 has provided a model instance of bunch division, and wherein spaced antenna 1 and spaced antenna 7 are a bunch representative antennas;
(2) from 7 spaced antennas of above-mentioned each bunch, choose spaced antenna
Figure GSA00000141617300093
With
Figure GSA00000141617300094
As a bunch representative antennas, then adjacent 3 bunches produce 6 bunches of representative antennas;
(3) all bunch representative antennas are carried out the distributed bidirectional cooperation synchronously, detailed process is as follows:
(3-1) above-mentioned all bunch representative antennas are numbered the back ordering, form a bunch representative antennas set
Figure GSA00000141617300095
Be abbreviated as { A 1, A 2, A 3, A 4, A 5, A 6, as shown in Figure 3;
(3-2) by a bunch representative antennas A 1To a bunch representative antennas A 2Send sinusoidal signal, A by that analogy 1→ A 2→ ... → A 5→ A 6, be designated as fl transmission; After the forward direction end of transmission, by a bunch representative antennas A 6Send identical sinusoidal signal to bunch representative antennas A 5, A by that analogy 6→ A 5→ ... → A 2→ A 1, being designated as reverse transfer, fl transmission and reverse transfer take 10 units altogether and propagate time slot, wherein time slot { T (1), T (2), T (3), T (4), T (5)Be that the fl transmission time slot is gathered time slot { T (6), T (7), T (8), T (9), T (10)Be that the reverse transfer time slot is gathered;
(3-3) at above-mentioned fl transmission time slot T (1)In, bunch representative antennas A 1The signal indication that sends is:
s 1 (1)(t)=exp{j(w 1t 11)}
W wherein 1Be a bunch representative antennas A 1Local crystal oscillator frequency, Φ 1Phase place for local crystal oscillator;
Then, bunch representative antennas A 2Reception is from a bunch representative antennas A 1Signal, this signal is with bunch representative antennas A 2Local zone time t 2Be expressed as:
r 2 ( 1 ) ( t 2 ) = H 1,2 exp { j ( α 1 w 1 ( t 2 / α 2 + Δ 1 - Δ 2 ) + Φ 1 ) } + ω 2 ( 1 ) ( t 2 )
H wherein 1,2For from a bunch representative antennas A 1To bunch representative antennas A 2Channel Transmission gain, and set fl transmission and the back is symmetry equivalent to the channel transmitted transmission gain,
Figure GSA00000141617300097
Be a bunch representative antennas A 2At T (1)In the noise that receives, then bunch representative antennas A 2Estimate that the crystal oscillator frequency and the phase place that obtain are:
w ^ 2 ( 1 ) = ( α 1 w 1 + w ~ 2 ( 1 ) ) / α 2
Φ ^ 2 ( 1 ) = α 1 w 1 ( Δ 1 - Δ 2 ) + Φ 1 + Φ ~ 2 ( 1 )
Wherein
Figure GSA00000141617300102
With
Figure GSA00000141617300103
Be estimated frequency error and phase estimation error;
(3-4) repeat above-mentioned steps (3-2)-(3-3), obtain a bunch representative antennas A iFl transmission crystal oscillator frequency and the phase place estimated be:
w ^ i = α i - 1 w ^ i - 1 / α i = α 1 w 1 / α i
Φ ^ i = α 1 w 1 ( Δ i - 1 - Δ i ) + Φ ^ i - 1 = α 1 w 1 ( Δ 1 - Δ i ) + Φ 1
(3-5) reverse transfer is by a bunch representative antennas A 6Initiate, up to time slot T (10)In bunch representative antennas A 1Receive till the final signal;
(3-6) at above-mentioned reverse transfer time slot T (6)In, bunch representative antennas A 6The signal indication that sends is:
s 6 (6)(t)=exp{j(w 6t 66)} (36)
W wherein 6Be a bunch representative antennas A 6Local crystal oscillator frequency, Φ 6Phase place for local crystal oscillator;
Then, bunch representative antennas A 5Reception is from a bunch representative antennas A 6Signal, this signal is with bunch representative antennas A 5Local zone time t 5Be expressed as:
r 5 ( 6 ) ( t 5 ) = H 6,5 exp { j ( α 6 w 6 ( t 5 / α 5 + Δ 6 - Δ 5 ) + Φ 6 ) } + ω 5 ( 6 ) ( t 5 ) - - - ( 37 )
H wherein 6,5For from a bunch representative antennas A 6To bunch representative antennas A 5Transmission channel gain,
Figure GSA00000141617300107
Be a bunch representative antennas A 5At T (6)In the noise that receives, then bunch representative antennas A 5Estimate that the crystal oscillator frequency and the phase place that obtain are:
w ^ 5 ( 6 ) = ( α 6 w 6 + w ~ 5 ( 6 ) ) / α 5 - - - ( 38 )
Φ ^ 5 ( 6 ) = α 6 w 6 ( Δ 6 - Δ 5 ) + Φ 6 + Φ ~ 5 ( 6 ) - - - ( 39 )
Wherein
Figure GSA000001416173001010
With
Figure GSA000001416173001011
Be estimated frequency error and phase estimation error;
(3-7) repeat above-mentioned steps (3-5)-(3-6), obtain a bunch representative antennas A iReverse transfer crystal oscillator frequency and the phase place estimated be:
w ^ i = α i + 1 w ^ i + 1 / α i = α 6 w 6 / α i - - - ( 40 )
Φ ^ i = α 6 w 6 ( Δ i - 1 - Δ i ) + Φ ^ i + 1 = α 6 w 6 ( Δ 6 - Δ i ) + Φ 6 - - - ( 41 )
(3-8) according to the fl transmission estimated result and the reverse transfer estimated result of above-mentioned each that obtains bunch representative antennas, bunch representative antennas A then iFinal local crystal oscillator frequency and the phase place estimated are:
w ^ i = ( α 1 w 1 + α 6 w 6 ) / α i - - - ( 42 )
Φ ^ i = α 1 w 1 ( Δ 1 - Δ i ) + α 6 w 6 ( Δ 6 - Δ i ) + Φ 1 + Φ 6 - - - ( 43 )
And obtain bunch representative antennas A based on unified overall reference time t iLocal crystal oscillator signal be:
ψ i(t)=exp{j((α 1w 16w 6)t+λ 16)},i=1,2,3,4,5,6 (44)
λ wherein nnw nΔ n+ Φ n
(3-9) according to above-mentioned bunch of representative antennas A iLocal crystal oscillator frequency and the estimated value of phase place, obtain the set of bunch representative antennas The synchronized result of local crystal oscillator signal be:
ψ n ( p ) ( t ) = exp { j ( ( α 1 ( 1 ) w 1 ( 1 ) + α 7 ( 3 ) w 7 ( 3 ) ) t + λ 1 ( 3 ) + λ 7 ( 3 ) ) } , p=1,2,3,n=1,2,...,7 (45)
(4) according to the synchronized result of bunch representative antennas between above-mentioned each bunch, then, the distributed bidirectional cooperation between a plurality of spaced antennas in being initiated bunch by bunch representative antennas of each bunch synchronously, 7 spaced antennas with p bunch are example, p=1,2,3, concrete steps are as follows:
(4-1) after finishing the set synchronously of bunch representative antennas, to bunch in 7 spaced antennas be numbered the back ordering, form the spaced antenna set
Figure GSA00000141617300115
Be abbreviated as { B 1, B 2..., B 7, spaced antenna wherein
Figure GSA00000141617300116
Be two bunches of representative antennas of this bunch;
(4-2) by bunch in spaced antenna B 1To spaced antenna B 2Send sinusoidal signal, B by that analogy 1→ B 2→ ... → B 6→ B 7, be designated as fl transmission; After the forward direction end of transmission, by spaced antenna B 7Send identical sinusoidal signal to spaced antenna B 6, B by that analogy 7→ B 6→ ... → B 2→ B 1, being designated as reverse transfer, fl transmission and reverse transfer take 12 units altogether and propagate time slot, wherein time slot
Figure GSA00000141617300117
Be the set of fl transmission time slot, time slot
Figure GSA00000141617300118
Be the set of reverse transfer time slot;
(4-3) at above-mentioned fl transmission time slot
Figure GSA00000141617300119
In, spaced antenna B 1The signal indication that sends is:
x 1 (1)(t)=exp{j(w 1t 11)}
W wherein 1Be spaced antenna B 1Local crystal oscillator frequency, Φ 1Phase place for local crystal oscillator;
Spaced antenna B 2Reception is from spaced antenna B 1Signal, this signal spaced antenna B 2Local zone time t 2Be expressed as:
y 2 ( 1 ) ( t 2 ) = H 1,2 exp { j ( α 1 w 1 ( t 2 / α 2 + Δ 1 - Δ 2 ) + Φ 1 ) } + Φ 2 ( 1 ) ( t 2 )
H wherein 1,2For from spaced antenna B 1To spaced antenna B 2Channel Transmission gain, and set fl transmission and the back is symmetry equivalent to the channel transmitted transmission gain,
Figure GSA00000141617300122
Be spaced antenna B 2
Figure GSA00000141617300123
In the noise that receives, then spaced antenna B 2Estimate that the crystal oscillator frequency and the phase place that obtain are:
w ^ 2 ( 1 ) = ( α 1 w 1 + w ~ 2 ( 1 ) ) / α 2
Φ ^ 2 ( 1 ) = α 1 w 1 ( Δ 1 - Δ 2 ) + Φ 1 + Φ ~ 2 ( 1 )
Wherein
Figure GSA00000141617300126
With
Figure GSA00000141617300127
Be estimated frequency error and phase estimation error;
(4-4) repeat above-mentioned steps (4-2)-(4-3), spaced antenna B in the spaced antenna set in obtaining bunch iFl transmission crystal oscillator frequency and the phase place estimated be:
w ^ i = α i - 1 w ^ i - 1 / α i = α 1 w 1 / α i
Φ ^ i = α 1 w 1 ( Δ i - 1 - Δ i ) + Φ ^ i - 1 = α 1 w 1 ( Δ 1 - Δ i ) + Φ 1
(4-5) reverse transfer is by spaced antenna B 7Initiate, up to time slot
Figure GSA000001416173001210
Middle spaced antenna B 1Receive till the final signal;
(4-6) at above-mentioned reverse transfer time slot
Figure GSA000001416173001211
In, spaced antenna B 2The signal indication that sends is:
x 7 (7)(t)=exp{j(w 7t 77)} (46)
W wherein 7Be spaced antenna B 7Local crystal oscillator frequency, Φ 7Phase place for local crystal oscillator;
Spaced antenna B 6Reception is from spaced antenna B 7Signal, this signal spaced antenna B 6Local zone time t 6Be expressed as:
y 6 ( 7 ) ( t 6 ) = H 7,6 exp { j ( α 7 w 7 ( t 6 / α 6 + Δ 7 - Δ 6 ) + Φ 7 ) } + ω 6 ( 7 ) ( t 6 ) - - - ( 47 )
H wherein 7,6For from spaced antenna B 7To spaced antenna B 6Transmission channel gain,
Figure GSA000001416173001213
Be spaced antenna B 6
Figure GSA00000141617300131
In the noise that receives, then spaced antenna B 6Estimate that the frequency and the phase place that obtain are
w ^ 6 ( 7 ) = ( α 7 w 7 + w ~ 6 ( 7 ) ) / α 6 - - - ( 48 )
Φ ^ 6 ( 7 ) = α 7 w 7 ( Δ 7 - Δ 6 ) + Φ 7 + Φ ~ 6 ( 7 ) - - - ( 49 )
Wherein
Figure GSA00000141617300134
With Be estimated frequency error and phase estimation error;
(4-7) repeat above-mentioned steps (4-5)-(4-6), spaced antenna B in the spaced antenna set in obtaining bunch iReverse transfer crystal oscillator frequency and the phase place estimated be:
w ^ i = α i + 1 w ^ i + 1 / α i = α 7 w 7 / α i - - - ( 50 )
Φ ^ i = α 7 w 7 ( Δ i + 1 - Δ i ) + Φ ^ i + 1 = α 7 w 7 ( Δ 7 - Δ i ) + Φ 7 - - - ( 51 )
(4-8) according to above-mentioned obtain bunch in the fl transmission estimated result and the reverse transfer estimated result of spaced antenna, then spaced antenna B iFinal local crystal oscillator frequency and the phase place estimated are:
w ^ i = ( α 1 w 1 + α 7 w 7 ) / α i - - - ( 52 )
Φ ^ i = α 1 w 1 ( Δ 1 - Δ i ) + α 7 w 7 ( Δ 7 - Δ i ) + Φ 1 + Φ 7 - - - ( 53 )
And obtain spaced antenna B based on unified overall reference time t iLocal crystal oscillator signal be:
ψ i(t)=exp{j((α 1w 17w 7)t+λ 17)},i=1,2,...,7 (54)
λ wherein nnw nΔ n+ Φ n
(4-9) according to spaced antenna B in above-mentioned bunch iLocal crystal oscillator frequency and the estimated value of phase place, obtain p bunch bunch in the spaced antenna set The synchronized result of local crystal oscillator signal be:
ψ n ( p ) ( t ) = exp { j ( ( α 1 ( p ) w 1 ( p ) + α 7 ( p ) w 7 ( p ) ) t + λ 1 ( p ) + λ 7 ( p ) ) } , p=1,2,3,n=1,2,...,7 (55)
(4-10) because p bunch bunch representative antennas
Figure GSA000001416173001312
With a bunch representative antennas
Figure GSA000001416173001313
In above-mentioned (3) step with bunch representative antennas of other bunches finish local crystal oscillator synchronously, so in the whole system 3 bunches totally 21 spaced antenna deadline Frequency Synchronization.Be about to obtain in above-mentioned formula (45) the substitution formula (55) that the final synchronized result of all spaced antennas is in the whole system:
ψ n ( p ) ( t ) = exp { j [ 2 ( α 1 ( 1 ) w 1 ( 1 ) + α 7 ( 3 ) w 7 ( 3 ) ) t + 2 ( λ 1 ( 3 ) + λ 7 ( 3 ) ) ] } , p=1,2,3,n=1,2,...,7。(56)

Claims (1)

  1. In the radio communication based on bunch distributed antenna time frequency synchronizing method, it is characterized in that this method may further comprise the steps:
    (1) spaced antenna of the M=PN in the wireless communication system is divided into P bunch, p=1,2 ..., P, the spaced antenna number in each bunch are N, n=1, and 2 ..., N, n spaced antenna usefulness in p bunch Expression, then the spaced antenna set is in every bunch
    Figure FSA00000141617200012
    Setting up departments, the unified reference clock signal of spaced antenna is Ψ (t)=exp{j (Ω t+ Φ) in the system }, wherein Ω is the frequency of local crystal oscillator, and t is the unified overall reference time, and Φ is the phase place of local crystal oscillator.Owing to there is the deviation of local crystal oscillator, thus in the system any one spaced antenna i local crystal oscillator signal indication be:
    ψ i(t)=exp{j(Ω it+Φ i+n i(t))},i=1,2,...,M-1,M (1)
    Ω in the following formula iBe local crystal oscillator frequency, Φ iThe phase place of representing local crystal oscillator, n i(t) be the noise jamming of local crystal oscillator; According to above-mentioned (1) formula, obtain that the current local zone time of any one spaced antenna i is in the said system:
    t i=(Ω it+Φ i+n i(t))/Ω=α i(t+Δ i(t)),i=1,2,...,M-1,M (2)
    α wherein il/ Ω, the ratio of the local crystal oscillator frequency of any one spaced antenna i and unified overall crystal oscillator frequency in the expression system; Δ i(t)=(Φ i+ n i(t))/Ω i, the local zone time deviation of any one spaced antenna i in the expression system considers that the local clock deviation is a fixed value in a unit propagates time slot, then the local zone time t of any one spaced antenna i in the system iWith the overall situation reference time t relational expression be t i=t+ Δ i
    (2) from N the spaced antenna of above-mentioned each bunch, choose 2 spaced antennas
    Figure FSA00000141617200013
    With
    Figure FSA00000141617200014
    As a bunch representative antennas, then adjacent P bunch produces 2P bunch representative antennas;
    (3) above-mentioned all bunch representative antennas are carried out the distributed bidirectional cooperation synchronously, detailed process is as follows:
    (3-1) above-mentioned all bunch representative antennas are numbered the back ordering, form a bunch representative antennas set
    Figure FSA00000141617200015
    Be abbreviated as { A 1, A 2..., A 2P-1, A 2P;
    (3-2) by a bunch representative antennas A 1To a bunch representative antennas A 2Send sinusoidal signal, A by that analogy 1→ A 2→ ... → A 2P-1→ A 2P, be designated as fl transmission; After the forward direction end of transmission, by a bunch representative antennas A 2PSend identical sinusoidal signal to bunch representative antennas A 2P-1, A by that analogy 2P→ A 2P-1→ ... → A 2→ A 1, being designated as reverse transfer, fl transmission and reverse transfer take 4P-2 unit altogether and propagate time slot, wherein time slot { T (1), T (2)..., T (2P-2), T (2P-1)Be that the fl transmission time slot is gathered time slot { T (2P), T (2P+1)..., T (4P-3), T (4P-2)Be that the reverse transfer time slot is gathered; (3-3) at above-mentioned fl transmission time slot T (1)In, bunch representative antennas A 1The signal indication that sends is:
    s 1 (1)(t)=exp{j(w 1t 11)} (3)
    W wherein 1Be a bunch representative antennas A 1Local crystal oscillator frequency, Φ 1Phase place for local crystal oscillator;
    Bunch representative antennas A 2Reception is from a bunch representative antennas A 1Signal, this signal is with bunch representative antennas A 2Local zone time t 2Be expressed as:
    r 2 ( 1 ) ( t 2 ) = H 1,2 exp { j ( α 1 w 1 ( t 2 / α 2 + Δ 1 - Δ 2 ) + Φ 1 ) } + ω 2 ( 1 ) ( t 2 ) - - - ( 4 )
    H wherein 1,2For from a bunch representative antennas A 1To bunch representative antennas A 2Channel Transmission gain, and set fl transmission and the back is symmetry equivalent to the channel transmitted transmission gain,
    Figure FSA00000141617200022
    Be a bunch representative antennas A 2At T (1)In the noise that receives, then bunch representative antennas A 2Estimate that the crystal oscillator frequency and the phase place that obtain are:
    w ^ 2 ( 1 ) = ( α 1 w 1 + w ~ 2 ( 1 ) ) / α 2 - - - ( 5 )
    Φ ^ 2 ( 1 ) = α 1 w 1 ( Δ 1 - Δ 2 ) + Φ 1 + Φ ~ 2 ( 1 ) - - - ( 6 )
    Wherein
    Figure FSA00000141617200025
    With
    Figure FSA00000141617200026
    Be estimated frequency error and phase estimation error;
    (3-4) repeat above-mentioned steps (3-2)-(3-3), the crystal oscillator frequency and the phase place that obtain the fl transmission estimation of bunch representative antennas Ai are:
    w ^ i = α i - 1 w ^ i - 1 / α i = α 1 w 1 / α i - - - ( 7 )
    Φ ^ i = α 1 w 1 ( Δ i - 1 - Δ i ) + Φ ^ i - 1 = a 1 w 1 ( Δ 1 - Δ i ) + Φ 1 ; - - - ( 8 )
    (3-5) reverse transfer is by a bunch representative antennas A 2PBeginning is up at reverse transfer time slot T (4p-2)Intra-cluster representative antennas A 1Receive till the final signal;
    (3-6) at above-mentioned reverse transfer time slot T (2P)In, bunch representative antennas A 2PThe signal indication that sends is:
    s 2P (2P)(t)=exp{j(w 2Pt 2P2P)} (9)
    W wherein 2PBe a bunch representative antennas A 2PLocal crystal oscillator frequency, Φ 2PPhase place for local crystal oscillator;
    Bunch representative antennas A 2P-1Reception is from a bunch representative antennas A 2PSignal, this signal is with bunch representative antennas A 2P-1Local zone time t 2p-1Be expressed as:
    r 2 P - 1 ( 2 P ) ( t 2 P - 1 ) = H 2 P , 2 P - 1 exp { j ( α 2 P w 2 P ( t 2 P - 1 / α 2 P - 1 + Δ 2 P - Δ 2 P - 1 ) + Φ 2 P ) } + ω 2 P - 1 ( 2 P ) ( t 2 P - 1 ) - - - ( 10 )
    H wherein 2P, 2P-1For from a bunch representative antennas A 2PTo bunch representative antennas A 2P-1Transmission channel gain,
    Figure FSA00000141617200032
    Be a bunch representative antennas A 2P-1At T (2P)In the noise that receives, then bunch representative antennas A 2P-1Estimate that the crystal oscillator frequency and the phase place that obtain are:
    w ^ 2 P - 1 ( 2 P ) = ( α 2 P w 2 P + w ^ 2 P - 1 ( 2 P ) ) / α 2 P - 1 - - - ( 11 )
    Φ ^ 2 P - 1 ( 2 P ) = α 2 P w 2 P ( Δ 2 P - Δ 2 P - 1 ) + Φ 2 P + Φ ~ 2 P - 1 ( 2 P ) - - - ( 12 )
    Wherein
    Figure FSA00000141617200035
    With
    Figure FSA00000141617200036
    Be estimated frequency error and phase estimation error;
    (3-7) repeat above-mentioned steps (3-5)-(3-6), obtain a bunch representative antennas A iReverse transfer crystal oscillator frequency and the phase place estimated be:
    w ^ i = α i + 1 w ^ i + 1 / α i = α 2 P w 2 P / α i - - - ( 13 )
    Φ ^ i = α 2 P w 2 P ( Δ i + 1 - Δ i ) + Φ ^ i + 1 = α 2 P w 2 P ( Δ 2 P - Δ i ) + Φ 2 P ; - - - ( 14 )
    (3-8) according to the fl transmission estimated result and the reverse transfer estimated result of above-mentioned each that obtains bunch representative antennas, bunch representative antennas A then iFinal local crystal oscillator frequency and the phase place estimated are:
    w ^ i = ( α 1 w 1 + α 2 P w 2 P ) / α i - - - ( 15 )
    Φ ^ i = α 1 w 1 ( Δ 1 - Δ i ) + α 2 P w 2 P ( Δ 2 P - Δ i ) + Φ 1 + Φ 2 P - - - ( 16 )
    And obtain bunch representative antennas A based on unified overall reference time t iLocal crystal oscillator signal be:
    ψ i(t)=exp{j((α 1w 12Pw 2P)t+λ 12P)},i=1,2,...,2P-1,2P (17)
    λ wherein nnw nΔ n+ Φ n
    (3-9) according to above-mentioned bunch of representative antennas A iLocal crystal oscillator frequency and the estimated value of phase place, obtain the set of bunch representative antennas
    Figure FSA000001416172000311
    The synchronized result of local crystal oscillator signal be:
    ψ n ( p ) ( t ) = exp { j ( ( α 1 ( 1 ) w 1 ( 1 ) + α N ( P ) w N ( P ) ) t + λ 1 ( P ) + λ N ( P ) ) } , p = 1,2 , . . . , P , n = 1,2 , . . . , N ; - - - ( 18 )
    (4) according to the synchronized result of bunch representative antennas between above-mentioned each bunch, then, distributed bidirectional cooperation between a plurality of spaced antennas in being initiated bunch by bunch representative antennas of each bunch synchronously, the same in default wherein and method for synchronous and the above-mentioned steps (3) is example with p bunch N spaced antenna, p=1,2, ..., P, concrete steps are as follows:
    (4-1) after finishing the set synchronously of bunch representative antennas, to bunch in N spaced antenna be numbered the back ordering, form spaced antenna and gather
    Figure FSA00000141617200041
    Be abbreviated as { B 1, B 2..., B N, spaced antenna wherein
    Figure FSA00000141617200042
    Figure FSA00000141617200043
    Be two bunches of representative antennas of this bunch;
    (4-2) by bunch in spaced antenna B 1To spaced antenna B 2Send sinusoidal signal, B by that analogy 1→ B 2→ ... → B N-1→ B N, be designated as fl transmission; After the forward direction end of transmission, by spaced antenna B NSend identical sinusoidal signal to spaced antenna B N-1, B by that analogy N→ B N-1→ ... → B 2→ B 1, being designated as reverse transfer, fl transmission and reverse transfer take 2N-2 unit altogether and propagate time slot, wherein time slot
    Figure FSA00000141617200044
    Be the set of fl transmission time slot, time slot Be the set of reverse transfer time slot;
    (4-3) at above-mentioned fl transmission time slot In, the signal indication that spaced antenna B1 sends is:
    x 1 (1)(t)=exp{j(w 1t 11)} (19)
    W wherein 1Be spaced antenna B 1Local crystal oscillator frequency, Φ 1Phase place for local crystal oscillator;
    Spaced antenna B 2Reception is from spaced antenna B 1Signal, this signal spaced antenna B 2Local zone time t 2Be expressed as:
    y 2 ( 1 ) ( t 2 ) = H 1 , 2 exp { j ( α 1 w 1 ( t 2 / α 2 + Δ 1 - Δ 2 ) + Φ 1 ) } + ω 2 ( 1 ) ( t 2 ) - - - ( 20 )
    H wherein 1,2For from spaced antenna B 1To spaced antenna B 2Channel Transmission gain, and set fl transmission and the back is symmetry equivalent to the channel transmitted transmission gain,
    Figure FSA00000141617200048
    Be spaced antenna B 2
    Figure FSA00000141617200049
    In the noise that receives, then spaced antenna B 2Estimate that the crystal oscillator frequency and the phase place that obtain are:
    w ^ 2 ( 1 ) = ( α 1 w 1 + w ~ 2 ( 1 ) ) / α 2 - - - ( 21 )
    Φ ^ 2 ( 1 ) = α 1 w 1 ( Δ 1 - Δ 2 ) + Φ 1 + Φ ~ 2 ( 1 ) - - - ( 22 )
    Wherein
    Figure FSA000001416172000412
    With
    Figure FSA000001416172000413
    Be estimated frequency error and phase estimation error;
    (4-4) repeat above-mentioned steps (4-2)-(4-3), spaced antenna B in the spaced antenna set in obtaining bunch iFl transmission crystal oscillator frequency and the phase place estimated be:
    w ^ i = α i - 1 w ^ i - 1 / α i = α 1 w 1 / α i - - - ( 23 )
    Φ ^ i = α 1 w 1 ( Δ i - 1 - Δ i ) + Φ ^ i - 1 = α 1 w 1 ( Δ 1 - Δ i ) + Φ 1 ; - - - ( 24 )
    (4-5) reverse transfer is by spaced antenna B NInitiate, up to time slot
    Figure FSA00000141617200053
    Middle spaced antenna B 1Receive till the final signal;
    (4-6) at above-mentioned reverse transfer time slot
    Figure FSA00000141617200054
    In, spaced antenna B NThe signal indication that sends is:
    x N (N)(t)=exp{j(w Nt NN)}(25)
    W wherein NBe spaced antenna B NLocal crystal oscillator frequency, Φ NPhase place for local crystal oscillator;
    Spaced antenna B N-1Reception is from spaced antenna B NSignal, this signal spaced antenna B N-1Local zone time t N-1Be expressed as:
    y N - 1 ( N ) ( t N - 1 ) = H N , N - 1 exp { j ( α N w N ( t N - 1 / α N - 1 + Δ N - Δ N - 1 ) + Φ N ) } + ω N - 1 ( N ) ( t N - 1 ) - - - ( 26 )
    H wherein N, N-1For from spaced antenna B NTo spaced antenna B N-1Transmission channel gain,
    Figure FSA00000141617200056
    Be spaced antenna B N-1
    Figure FSA00000141617200057
    In the noise that receives, then spaced antenna B N-1Estimate that the frequency and the phase place that obtain are
    w ^ N - 1 ( N ) = ( α N w N + w ~ N - 1 ( N ) ) / α N - 1 - - ( 27 )
    Φ ^ N - 1 ( N ) = α N w N ( Δ N - Δ N - 1 ) + Φ N + Φ ^ N - 1 ( N ) - - - ( 28 )
    Wherein
    Figure FSA000001416172000510
    With
    Figure FSA000001416172000511
    Be estimated frequency error and phase estimation error;
    (4-7) repeat above-mentioned steps (4-5)-(4-6), spaced antenna B in the spaced antenna set in obtaining bunch iReverse transfer crystal oscillator frequency and the phase place estimated be:
    w ^ i = α i + 1 w ^ i + 1 / α i = α N w N / α i - - - ( 29 )
    Φ ^ i = α N w N ( Δ i + 1 - Δ i ) + Φ ^ i + 1 = α N w N ( Δ N - Δ i ) + Φ N ; - - - ( 30 )
    (4-8) according to above-mentioned obtain bunch in the fl transmission estimated result and the reverse transfer estimated result of spaced antenna, then spaced antenna B iFinal local crystal oscillator frequency and the phase place estimated are:
    w ^ i = ( α 1 w 1 + α N w N ) / α i - - - ( 31 )
    Φ ^ i = α 1 w 1 ( Δ 1 - Δ i ) + α N w N ( Δ N - Δ i ) + Φ 1 + Φ N - - - ( 32 )
    And obtain spaced antenna B based on unified overall reference time t iLocal crystal oscillator signal be:
    ψ i(t)=exp{j((α 1w 1Nw N)t+λ 1N)},i=1,2,...,N-1,N (33)
    λ wherein nnw nΔ n+ Φ n
    (4-9) according to spaced antenna B in above-mentioned bunch iLocal crystal oscillator frequency and the estimated value of phase place, obtain p bunch bunch in the spaced antenna set
    Figure FSA00000141617200062
    The synchronized result of local crystal oscillator signal be:
    ψ n ( p ) ( t ) = exp { j ( ( α 1 ( p ) w 1 ( p ) + α N ( p ) w N ( p ) ) t + λ 1 ( p ) + λ N ( p ) ) } , p = 1,2 , . . . , P , n = 1,2 , . . . , N ; - - - ( 34 )
    (4-10) because p bunch bunch representative antennas
    Figure FSA00000141617200064
    With a bunch representative antennas
    Figure FSA00000141617200065
    In above-mentioned (3) step, finish the synchronous of local crystal oscillator, so the P in the whole system bunch is total to M spaced antenna deadline Frequency Synchronization with bunch representative antennas of other bunches.Be about to obtain in above-mentioned formula (18) the substitution formula (34) that the final synchronized result of all spaced antennas is in the whole system:
    ψ n ( p ) ( t ) = exp { j [ 2 ( α 1 ( 1 ) w 1 ( 1 ) + α N ( P ) w N ( P ) ) t + 2 ( λ 1 ( P ) + λ N ( P ) ) ] } , p = 1,2 , . . . , P , n = 1,2 , . . . , N . - - - ( 35 )
CN2010101856887A 2010-05-21 2010-05-21 Cluster-based distributed antenna time frequency synchronizing method in wireless communication Expired - Fee Related CN101868022B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101856887A CN101868022B (en) 2010-05-21 2010-05-21 Cluster-based distributed antenna time frequency synchronizing method in wireless communication

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101856887A CN101868022B (en) 2010-05-21 2010-05-21 Cluster-based distributed antenna time frequency synchronizing method in wireless communication

Publications (2)

Publication Number Publication Date
CN101868022A true CN101868022A (en) 2010-10-20
CN101868022B CN101868022B (en) 2012-11-14

Family

ID=42959523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101856887A Expired - Fee Related CN101868022B (en) 2010-05-21 2010-05-21 Cluster-based distributed antenna time frequency synchronizing method in wireless communication

Country Status (1)

Country Link
CN (1) CN101868022B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102143572A (en) * 2011-01-24 2011-08-03 清华大学 Frequency synchronization method for cluster-based two-stage distributed antennae in wireless communication
CN102523619A (en) * 2011-12-15 2012-06-27 西安交通大学 Distributed coherent cooperation method in long distance wireless communication system
CN103327602A (en) * 2012-11-20 2013-09-25 电子科技大学 Carrier synchronization method and system for distributed wireless network
CN104022990A (en) * 2014-06-20 2014-09-03 哈尔滨工业大学 Distributed beam forming carrier phase synchronization method based on sea wireless sensor network
CN113038589A (en) * 2021-03-04 2021-06-25 重庆邮电大学 Matrix model estimation time synchronization method based on wireless network clustering topology

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101184004A (en) * 2007-03-16 2008-05-21 中科院嘉兴中心微系统所分中心 Double cluster based wireless sensor network distributed topology control method
US20080150514A1 (en) * 2006-12-21 2008-06-26 Nokia Corporation Communication method and system
CN101227242A (en) * 2008-01-31 2008-07-23 西安交通大学 Method for forming distributed aerial array beam based on channel correction

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080150514A1 (en) * 2006-12-21 2008-06-26 Nokia Corporation Communication method and system
CN101184004A (en) * 2007-03-16 2008-05-21 中科院嘉兴中心微系统所分中心 Double cluster based wireless sensor network distributed topology control method
CN101227242A (en) * 2008-01-31 2008-07-23 西安交通大学 Method for forming distributed aerial array beam based on channel correction

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《中国科学技术大学学报》 20090228 施冠超,邓单,朱近康 一种适用于分布式MIMO-OFDM系统的时间同步算法 全文 1 第39卷, 第2期 2 *
《北京邮电大学学报》 20071231 张一衡,崔琪楣,张平,陶小峰 分布式多天线系统中的TDOA估计 全文 1 第30卷, 第6期 2 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102143572A (en) * 2011-01-24 2011-08-03 清华大学 Frequency synchronization method for cluster-based two-stage distributed antennae in wireless communication
CN102143572B (en) * 2011-01-24 2013-12-11 清华大学 Frequency synchronization method for cluster-based two-stage distributed antennae in wireless communication
CN102523619A (en) * 2011-12-15 2012-06-27 西安交通大学 Distributed coherent cooperation method in long distance wireless communication system
CN102523619B (en) * 2011-12-15 2014-04-23 西安交通大学 Distributed coherent cooperation method in long distance wireless communication system
CN103327602A (en) * 2012-11-20 2013-09-25 电子科技大学 Carrier synchronization method and system for distributed wireless network
CN103327602B (en) * 2012-11-20 2016-12-28 电子科技大学 A kind of carrier synchronization method for distributed wireless networks and system
CN104022990A (en) * 2014-06-20 2014-09-03 哈尔滨工业大学 Distributed beam forming carrier phase synchronization method based on sea wireless sensor network
CN104022990B (en) * 2014-06-20 2017-06-20 哈尔滨工业大学 A kind of distributed beams based on sea wireless sense network form carrier phase synchronization method
CN113038589A (en) * 2021-03-04 2021-06-25 重庆邮电大学 Matrix model estimation time synchronization method based on wireless network clustering topology
CN113038589B (en) * 2021-03-04 2022-07-22 重庆邮电大学 Matrix model estimation time synchronization method based on wireless network clustering topology

Also Published As

Publication number Publication date
CN101868022B (en) 2012-11-14

Similar Documents

Publication Publication Date Title
CN100464612C (en) Digital cluster communication method with the control channel frequency hopping function
CN109041196B (en) Resource joint allocation method based on energy efficiency maximization in NOMA energy-carrying communication system
CN103402251B (en) Synchronizing information transmitting and receiving method, channel mapping and analyzing method and control information transmitting method
CN101868022B (en) Cluster-based distributed antenna time frequency synchronizing method in wireless communication
CN1909533B (en) Frame creating method based on crossing frequency division multiplexing in time-division duplex mode
CN104685923A (en) Method and device for allocating bandwidth, user equipment and base station
CN101951307A (en) Method for selecting cell cooperation set under CoMP
CN107580370A (en) A kind of wireless take based on OFDM amplification forward collaborations can communication means
CN103052111A (en) Method for selecting interuser communication modes under different base stations
CN109005592A (en) Power distribution method in single antenna NOMA system
CN102946295A (en) Frame structure capable of being matched with different carrier bandwidth and application thereof
CN109714806A (en) A kind of wireless power junction network optimization method of non-orthogonal multiple access
CN102143572B (en) Frequency synchronization method for cluster-based two-stage distributed antennae in wireless communication
CN103298097A (en) Combined relay selection and power distribution method implemented in bidirectional relay network
CN110113812A (en) A kind of wireless communication whole network synchronous method based on TD-LTE
CN1905405A (en) Space interface synchronous method
CN103889043A (en) Power distribution method used in cognitive wireless relay network
CN104168659A (en) Multi-cell MIMO system user scheduling method under MRT pre-coding strategy
WO2014166195A1 (en) Closely distributed wireless communication method and system therefor
CN102938754B (en) Synchronization sequence construction method in super cellular mobile communication device-to-device technology
CN113099461B (en) Symbiotic radio network design method based on non-orthogonal multiple access technology
CN103024752B (en) The system of selection of user collaboration set of cells and system
CN104066098A (en) Method used for distributed wireless network and equipment thereof
CN101136724A (en) Wireless self-organizing cooperative communication method
CN106160814B (en) A kind of multi-base station cooperative downlink transmission method based on interference alignment techniques

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20121114

Termination date: 20180521