CN100367690C - Frequency deviation estimation and correction method in TD-SCDMA system - Google Patents

Frequency deviation estimation and correction method in TD-SCDMA system Download PDF

Info

Publication number
CN100367690C
CN100367690C CNB2004100091755A CN200410009175A CN100367690C CN 100367690 C CN100367690 C CN 100367690C CN B2004100091755 A CNB2004100091755 A CN B2004100091755A CN 200410009175 A CN200410009175 A CN 200410009175A CN 100367690 C CN100367690 C CN 100367690C
Authority
CN
China
Prior art keywords
frequency offset
user
estimation
frequency deviation
frequency
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.)
Expired - Fee Related
Application number
CNB2004100091755A
Other languages
Chinese (zh)
Other versions
CN1705252A (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.)
Jinjiang High And New Technology Development Office
Original Assignee
ZTE Corp
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 ZTE Corp filed Critical ZTE Corp
Priority to CNB2004100091755A priority Critical patent/CN100367690C/en
Publication of CN1705252A publication Critical patent/CN1705252A/en
Application granted granted Critical
Publication of CN100367690C publication Critical patent/CN100367690C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

The present invention relates to a frequency deviation estimation and correction method in a TD-SCDMA system. The present invention comprises the steps that the frequency deviation of each user, which is estimated by the previous sub-frame, is used as the frequency deviation estimation of the sub-frame to generate a local training sequence having the frequency deviation; the local training sequence and the receiving training sequence make channel estimation; the frequency deviation of each user, which is estimated by the previous sub-frame, and a spread spectrum scrambling composite sequence generated at a location are processed to form the spread spectrum scrambling composite sequence having the frequency deviation; the spread spectrum scrambling composite sequence is used for making joint detection to estimate a symbol; the symbol comprises a linear phase caused by a frequency deviation difference of the sub-frame and the previous sub-frame of the user; the linear phase is used for estimating the frequency deviation difference, and the frequency deviation difference is used for correcting the phase; the symbol of which the phase is corrected is output. The present invention has the advantage that the aim for improving frequency deviation correcting performance is achieved by the channel estimation and the joint detection.

Description

Frequency offset estimation and correction method in TD-SCDMA system
Technical Field
The invention belongs to the technical field of mobile communication, and particularly relates to a frequency offset estimation and correction method in a TD-SCDMA system.
Background
In an actual mobile communication system, since a transmitter and a receiver use independent clocks, and there is a relative speed between a transmitting end and a receiving end, a frequency offset exists between a signal carrier received at the receiving end and a local carrier. The presence of such frequency offsets may result in a degradation of reception performance, and it is therefore necessary to correct such frequency offsets. The conventional frequency offset correction method uses a phase-locked loop circuit. US5303257 in 1994 proposes a digital frequency offset correction method based on analysis and processing of a received digital signal to perform frequency offset correction.
In TD-SCDMA system, symbol level digital frequency offset correction method is generally used, i.e. the de-spread symbol is analyzed and processed, frequency offset information is extracted from the symbol, and frequency offset correction is performed, and the system block diagram is shown in fig. 1. The received chip-level signal 101 is a superposition of multiple paths for multiple users, including a training sequence. The channel estimation 102 is performed by using the received training sequence 101 and the locally generated training sequence 103, then the joint detection 105 is performed on the received chip-level user data signal 101 by using the locally generated spreading and scrambling sequence 106 and the channel estimation value 104, a symbol 107 is estimated, then the frequency offset estimation and correction 108 are performed on the symbol 107, and the corrected symbol 109 is output. This method is effective when the frequency offset is small. The principle of the method is that when the frequency offset is small, the channel estimation value 104 is less affected by the frequency offset and is accurate, so that the phase of the jointly detected symbol 107 includes a linear phase caused by the frequency offset, and the phase of the symbol is extracted to perform frequency offset estimation and correction 108. However, when the frequency offset is large, the performance of this method becomes poor. The reason is that when the frequency offset is large, the channel estimation value 104 is affected by the frequency offset and becomes inaccurate enough, so that the phase error of the jointly detected symbol 107 is large, and the linear phase caused by the frequency offset cannot be accurately extracted from the phase of the symbol 107, so that the frequency offset cannot be effectively compensated, and the performance is deteriorated.
Disclosure of Invention
The invention aims to provide a novel frequency offset estimation and correction method in a TD-SCDMA system, which improves the frequency offset correction performance in large frequency offset.
The purpose of the invention is realized by the following steps:
the invention discloses a frequency offset estimation and correction method in a TD-SCDMA system, which is used for estimating and correcting the frequency offset existing between a received signal carrier and a local carrier, and comprises the following steps:
(1) The frequency offset of each user estimated by the previous subframe is used as the frequency offset estimation of the subframe to generate a training sequence generated locally to form a local training sequence with frequency offset, and then the local training sequence with frequency offset and the received training sequence with frequency offset are used for channel estimation;
(2) Processing the frequency offset of each user estimated by the previous subframe and a locally generated spreading and scrambling code composite sequence to form a spreading and scrambling code composite sequence with frequency offset, and performing joint detection by using the spreading and scrambling code composite sequence with frequency offset to estimate a symbol;
(3) Extracting symbol estimation of a certain user from the estimated symbols of all users, wherein the symbol estimation comprises a linear phase caused by the difference of the frequency offset of the current subframe and the previous subframe of the user, and the frequency offset difference can be estimated by utilizing the linear phase and is used for correcting the phase; the phase corrected symbol is output;
(4) And adding the estimated frequency deviation with the frequency deviation of the previous subframe, storing the frequency deviation as the frequency deviation estimation of the subframe, and using the frequency deviation estimation and correction of the next subframe.
The step (1) comprises the following steps:
(11) The training sequences with frequency offset of all users in a time slot are superposed after passing through a multipath channel to obtain a training sequence signal of the time slot:
e m =M fo h+n
wherein e m Is N m X 1 column vector, nm length of training sequence, h KW x 1 column vector formed by channel impulse response of all users in a time slot, K number of users, W length of channel impulse response of one user in chip unit, N N m Noise column vector of x 1, M fo Is a training sequence matrix with frequency offset, dimension N m ×KW;
Taking N m = KW, can obtain M fo
Figure C20041000917500111
Wherein, Δ f k And phi k Is the frequency offset and initial phase of the kth user, tc is the chip width, i is the position of the first chip of the training sequence in a time slot,is a transmitted training sequence;
(12) Using the frequency offset estimated by the last sub-frame as the frequency offset estimation of the sub-frame to obtain the local training sequence matrix estimation with frequency offset
Figure C20041000917500122
Figure C20041000917500123
Wherein
Figure C20041000917500124
Is the frequency offset 211 of the kth user estimated in the previous subframe;
(13) By using
Figure C20041000917500125
Estimated channel
Figure C20041000917500126
Figure C20041000917500127
Figure C20041000917500128
If frequency deviation of each user is delta f k If the phase difference is very small, then
Figure C20041000917500129
Approximating to form a circular matrix, thereby reducing the channel by Fourier transform
Figure C200410009175001210
The amount of computation of the term.
The step (2) comprises the following steps:
(21) Superposing the data signals with frequency offset of all users in a time slot after passing through a multipath channel to obtain the user data signal of the time slot:
e d =A fo d+n
wherein e is d Is (N) s Q + W-1) × 1 column vector, ns isThe number of user symbols, Q is the spreading factor, d is the column vector formed by all user symbols, and the dimension is N s K × 1,n is N s K × 1 noisy column vector, A fo Is a system transmission matrix with dimensions of (N) s Q+W-1)×N s K;
A fo The structure of (a) is as follows:
Figure C20041000917500131
A fo is a diagonal array of blocks, block V fo (n) Can be expressed as:
Figure C20041000917500132
n=1,…,N s
wherein, b k,fo (n) Is the convolution of the spreading scrambling code composite sequence of the k user with frequency offset and the channel impulse response of the k user, the dimension is (Q + W-1) x 1, and is expressed as:
Figure C20041000917500133
Figure C20041000917500134
wherein, c fo k Is the composite sequence of the spreading scrambling code of the k user with frequency deviation, and the dimension is (Q + W-1) xW, h k Is the channel impulse response of the kth user with dimension W × 1.
Locally generated spreading scrambling code composite sequence c of k user k Expressed as:
Figure C20041000917500141
(22) The frequency deviation estimated by the previous subframe is used as the frequency deviation estimation of the subframe, thereby generating the frequency deviation-containing spread spectrum scrambling code composite sequence estimation
Figure C20041000917500142
Wherein the content of the first and second substances,
Figure C20041000917500144
is the frequency offset estimation of the k user of the previous subframe;
(23) Calculated from step (13)Extracting the channel impulse response estimation of each user
Figure C20041000917500146
K =1, …, K, andand
Figure C20041000917500148
synthesis of
Figure C20041000917500149
V fo (n) And systematic transmission matrix estimation with frequency offset
Figure C200410009175001410
Figure C200410009175001411
k=1,…,K,n=1,…,N s
Figure C20041000917500151
Figure C20041000917500152
Minimum mean square error, MMSE, estimation of all user symbols
Figure C20041000917500153
Comprises the following steps:
Figure C20041000917500154
wherein R is n Is a noise correlation matrix, R d Is the correlation matrix of the user symbols.
The invention discloses a frequency offset estimation and correction method in a TD-SCDMA system, which is used for estimating and correcting the frequency offset existing between a received signal carrier and a local carrier, and comprises the following steps:
(1) Receiving a chip-level signal;
(2) Extracting from the chip-level signal a training sequence em and a user data part ed:
(3) Generating a local training sequence matrix M:
Figure C20041000917500155
(4) Frequency offset estimated from the last subframe
Figure C20041000917500156
K =1, …, K, which is the frequency offset estimation of the subframe:
(5) Using local training sequence matrices M and
Figure C20041000917500157
k =1, …, K, producing a local training sequence matrix estimate with frequency offset
Figure C20041000917500158
Figure C20041000917500161
(6) Using training sequences e m And
Figure C20041000917500162
performing channel estimation to obtain a channel estimation value
Figure C20041000917500163
Figure C20041000917500164
Figure C20041000917500165
(7) Generating a local spread spectrum scrambling code composite sequence matrix C k ,k=1,…,K:
Figure C20041000917500166
(8) Using c k K =1, …, K and
Figure C20041000917500167
k =1, …, K, producing a local spread spectrum scrambling code composite sequence matrix estimate with frequency offset
Figure C20041000917500168
k=1,…,K,n=1,…,N s
Figure C20041000917500171
(9) From
Figure C20041000917500172
Extracting channel estimation value of each user
Figure C20041000917500173
K =1, … K; by usingk=1,…,K,n=1,…,N s And
Figure C20041000917500175
k =1, …, K productionk=1,…,K,n=1,…,N s
Figure C20041000917500177
k=1,…,K,n=1,…,N s
(10) By using
Figure C20041000917500178
k=1,…,K,n=1,…,N s Composition V fo (n) ,n=1,…,N s
Figure C20041000917500179
n=1,…,N s
(11) By usingV fo (n) ,n=1,…,N s Forming system transmission matrix with frequency deviation
Figure C200410009175001710
Figure C200410009175001711
(12) By using
Figure C200410009175001712
And e d Performing joint detection to estimate user symbols
Figure C200410009175001713
Figure C20041000917500181
(13) By using
Figure C20041000917500182
Performing frequency offset estimation to estimate the frequency offset difference between the current subframe and the previous subframe
Figure C20041000917500183
k=1,…,K;
(14) By usingK =1, …, K pairs of user symbols
Figure C20041000917500185
Carrying out frequency offset correction and outputting a corrected symbol;
(15) Adding the obtained frequency deviation and the estimated frequency deviation of the previous subframe to be used as the estimated frequency deviation of the next subframe:
Figure C20041000917500186
k=1,…,K。
in short, the invention carries out channel estimation and joint detection by using the training sequence with frequency offset and the spreading scrambling code composite sequence with frequency offset, thereby realizing the following two points and further achieving the purpose of improving the frequency offset correction performance. First, a channel is estimated more accurately; and secondly, converting the estimated large frequency offset into a small frequency offset of two adjacent subframes.
Drawings
FIG. 1 is a basic block diagram of symbol-level frequency offset estimation and correction of a prior TD-SCDMA system;
fig. 2 is a system diagram of frequency offset estimation and correction of TD-SCDMA system discussed in the present invention.
Detailed Description
The system block diagram of the invention is shown in figure 2. The basic idea is as follows.
And taking the frequency offset 211 of each user estimated in the previous subframe as the frequency offset estimation of the subframe, substituting the frequency offset 211 into the training sequence 203 generated locally to form a local training sequence 204 with frequency offset, and performing channel estimation by using the local training sequence 204 with frequency offset and the received training sequence 201 with frequency offset. Generally, it can be considered that the frequency offset between two adjacent subframes is not changed much by 5 ms. Under such conditions, the estimated channel values 205 are more accurate than those estimated using the local training sequence without frequency offset.
Then, the frequency offset 211 of each user estimated in the previous sub-frame is substituted into the locally generated spreading and scrambling code composite sequence 207 to form a spreading and scrambling code composite sequence 208 with frequency offset, and then the spreading and scrambling code composite sequence with frequency offset is used for joint detection 206. The symbols thus estimated include the linear phase caused by the difference between the frequency offset of the present subframe and the frequency offset of the previous subframe. Since the difference between the frequency offset of the present subframe and the frequency offset of the previous subframe can be considered to be relatively small, the frequency compensation can be performed better by the frequency offset estimation and correction 210 at the following symbol level.
The present invention will be described in detail below.
First, channel estimation
Received training sequence signal e of a time slot m 201 is the superposition of the training sequences with frequency offset of all users in the time slot after passing through the multipath channel, which can be expressed as
e m =M fo h+n(1)
Wherein e m Is N m X 1 column vector, nm length of training sequence, h KW x 1 column vector formed by channel impulse response of all users in a time slot, K number of users, W length of channel impulse response of one user in chip unit, N N m Noise column vector of x 1, M fo Is a training sequence matrix with frequency deviation and dimension N m X KW. When designing a system, to simplify the calculation, generally take N m =KW。 M fo Can be expressed as
Figure C20041000917500191
Wherein Δ f k And phi k Is the frequency offset and initial phase of the kth user, tc is the chip width, i is the position of the first chip of the training sequence in a time slot,
Figure C20041000917500192
is the transmitted training sequence.
The locally generated training sequence matrix M203 is represented as
Figure C20041000917500193
It can be seen that M is a circulant matrix, i.e. each row is cyclically shifted to the right by one element, i.e. the next row. Using the frequency offset estimated from the previous sub-frame as the sub-frameFrequency offset estimation, thereby producing a local training sequence matrix estimate with frequency offset
Figure C20041000917500194
Figure C20041000917500201
Wherein
Figure C20041000917500202
The frequency offset 211 of the kth user estimated in the previous subframe.
By using
Figure C20041000917500203
Estimated channel
Figure C20041000917500204
Can be expressed as
Figure C20041000917500205
Figure C20041000917500206
And the channel estimated by the local training sequence matrix MCan be expressed as
Figure C20041000917500208
When frequency deviation is delta f k When the frequency offset is larger, under the condition that the frequency offset difference between two adjacent subframes is not large,
Figure C200410009175002010
ratio ofIs accurate.
In the formula (5), the following requirements are satisfied
Figure C200410009175002012
The calculation amount is large. If frequency offset of each user is delta f k If the phase difference is very small, then it will be
Figure C200410009175002013
An approximation process is performed to form a circulant matrix, thereby reducing the amount of calculation by using FFT (Fourier transform). Assuming frequency offset estimation for each user
Figure C200410009175002014
Average value of (2)
Figure C200410009175002015
Is composed of
Figure C200410009175002016
Will be provided with
Figure C200410009175002017
Make the following approximation
Figure C20041000917500211
Then
Figure C20041000917500212
Becomes a circulant matrix that can be computed using a Fourier transform FFT
Figure C20041000917500213
Figure C20041000917500214
Second, joint detection
Received user data signal e of one time slot d 201 is the superposition of the data signals with frequency offset of all users in the time slot after passing through the multipath channel, which can be expressed as
e d =A fo d+n
Wherein e d Is (N) s Q + W-1) x 1, ns is the number of user symbols, Q is the spreading factor, d is the column vector of all user symbols, and the dimension is N s K × 1,n is N s K × 1 noise column vector, afo is the system transmission matrix, with dimensions of (N) s Q+W-1)×N s K. Afo has the following structure:
Figure C20041000917500215
afo is a block diagonal matrix, block V fo (n) Can be expressed as
Figure C20041000917500216
n=1,…,N s (11)
Wherein b is k,fo (n) Is the convolution of the spreading scrambling code composite sequence of the k user with frequency deviation and the channel impulse response of the k user, the dimension is (Q + W-1) multiplied by 1, which is expressed as
Figure C20041000917500221
Figure C20041000917500222
Wherein c is fo k Is the composite sequence of the spreading scrambling code of the k user with frequency deviation, and the dimension is (Q + W-1) xW, h k Is the channel impulse response of the kth user with dimension W × 1.
The locally generated spreading scrambling code composite sequence ck207 of the k-th user is denoted as
Figure C20041000917500223
The frequency deviation 211 estimated by the previous subframe is used as the frequency deviation estimation of the subframe, thereby generating the frequency deviation-containing spread spectrum scrambling code composite sequence estimation
Figure C20041000917500224
Is shown as
Figure C20041000917500231
WhereinFrequency offset estimate 211 for the k user of the previous subframe. Calculated from the first step
Figure C20041000917500233
Extracting the channel impulse response estimation of each user
Figure C20041000917500234
K =1, …, K, and
Figure C20041000917500235
and
Figure C20041000917500236
synthesis of
Figure C20041000917500237
V fo (n) And systematic transmission matrix estimation with frequency offset
Figure C20041000917500238
Figure C20041000917500239
k=1,…,K,n=1,…,N s (15)
Figure C200410009175002310
n=1,…,N s (16)
MMSE estimation of all user symbols
Figure C20041000917500241
Is composed of
Figure C20041000917500242
Wherein R is n Is a noise correlation matrix, R d Is the correlation matrix of the user symbols. R n Can be estimated from measured data, R d Can be calculated by a priori knowledge.
Notably, the use of systematic transmission matrix estimation with frequency offset
Figure C20041000917500244
To perform symbol estimation, the estimated symbol
Figure C20041000917500245
The method includes a linear phase caused by the difference of the frequency offset of the current subframe and the previous subframe, but not a linear phase caused by the frequency offset of the current subframe.
Third, frequency offset estimation and correction
All user symbols estimated from the second step
Figure C20041000917500246
Extracting symbol estimation of the k user
Figure C20041000917500247
Figure C20041000917500248
The linear phase caused by the frequency offset difference between the current subframe and the previous subframe of the kth user is included, and the frequency offset difference between the current subframe and the previous subframe can be generally considered to be small (< 200 Hz). The frequency offset 210 can be estimated by using the linear phase and corrected by using the frequency offset
Figure C20041000917500249
Phase 210 of (a). The phase corrected symbol 213 is output. And meanwhile, adding the estimated frequency deviation with the frequency deviation of the previous subframe, storing the frequency deviation as the frequency deviation estimation of the subframe, and using the frequency deviation estimation and correction of the next subframe. There are many documents and patents on how to perform the frequency offset estimation and correction 210, which are not in the scope of the present invention.
The implementation steps of the technical scheme are further described in detail in the following steps with the combination of the accompanying drawings:
1. receiving a chip-level signal 201
2. The training sequence e is extracted from the received chip-level signal 201 m And a user data section e d
3. A local training sequence matrix M203 is generated, as in equation (3).
4. Frequency offset estimated from the last subframeK =1, …, K211 as the frequency offset estimate of the present subframe.
5. With M and
Figure C200410009175002411
k =1, … K, producing a local training sequence matrix estimate with frequency offsetSuch as formula 4.
6. With e m And
Figure C20041000917500251
channel estimation 202 is performed as shown in equation (5). Obtaining a channel estimate
Figure C20041000917500252
7. Generating a composite sequence matrix c of local spreading scrambling codes k K =1, …, K207, as in formula (13).
8. By c k K =1, …, K and
Figure C20041000917500253
k =1, … K, producing a local spread spectrum scrambling code composite sequence matrix estimate with frequency offset
Figure C20041000917500254
k=1,…,K,n=1,…N s 208, as in formula (14).
9. From
Figure C20041000917500255
Extracting channel estimation value of each user
Figure C20041000917500256
k =1, … k. By using
Figure C20041000917500257
k=1,…,K,n=1,…,N s And
Figure C20041000917500258
k =1, …, K productionk=1,…,K,n=1,…N s Such as formula 15.
10. By usingk=1,…,K,n=1,…,N s Composition V fo (n) ,n=1,…,N s Such as formula (16).
11. By V fo (n) ,n=1,…,N s Forming system transmission matrix with frequency deviation
Figure C200410009175002511
Such as formula (17).
12. By usingAnd e d Performing joint detection to estimate user symbols
Figure C200410009175002513
As in formula (18)
13. By using
Figure C200410009175002514
Performing frequency offset estimation 210 to estimate the frequency offset between the current subframe and the previous subframek=1,…,K
14. By using
Figure C200410009175002516
K =1, …, K pairs of user symbols
Figure C200410009175002517
Frequency offset correction 210 is performed and corrected symbol 213 is output.
15. Order to
Figure C200410009175002518
k=1,…,K
The invention achieves the purpose of improving the frequency offset correction performance in large frequency offset through the following two points:
1. and the local training sequence with frequency offset is used for channel estimation, so that the channel estimation is more accurate.
2. And performing symbol estimation by using the local spread spectrum scrambling code composite sequence with the frequency offset, so that the estimated symbol comprises a linear phase caused by the frequency offset difference between the current subframe and the previous subframe, and further converting the estimated large frequency offset into the small frequency offset of the adjacent subframe.

Claims (5)

1. A frequency deviation estimation and correction method in TD-SCDMA system is used for estimating and correcting the frequency deviation existing between the received signal carrier and the local carrier, and is characterized in that the method comprises the following steps:
step 1, the frequency offset of each user estimated by the previous subframe is used as the frequency offset estimation of the subframe to generate a local training sequence with frequency offset, and then the local training sequence with frequency offset and the received training sequence with frequency offset are used for channel estimation;
step 2, processing the frequency offset of each user estimated by the previous subframe and a locally generated spreading and scrambling code composite sequence to form a spreading and scrambling code composite sequence with frequency offset, and performing joint detection by using the spreading and scrambling code composite sequence with frequency offset to estimate a symbol;
step 3, extracting symbol estimation of a certain user from the estimated symbols of all users, wherein the symbol estimation comprises a linear phase caused by the difference of the frequency offset of the current subframe and the previous subframe of the user, the frequency offset difference can be estimated by utilizing the linear phase, and the phase is corrected by utilizing the frequency offset difference; the phase corrected symbol is output;
and 4, adding the estimated frequency deviation and the frequency deviation of the previous subframe, storing the frequency deviation as the frequency deviation estimation of the subframe, and using the frequency deviation estimation and correction of the next subframe.
2. The method of claim 1, wherein the step 1 comprises the steps of:
step 11, overlapping the training sequences with frequency offset of all users in a time slot after passing through a multipath channel to obtain a training sequence signal of the time slot:
e m =M fo h+n
wherein e m Is N m X 1 column vector, nm length of training sequence, h KW x 1 column vector formed by channel impulse response of all users in a time slot, K number of users, W length of channel impulse response of one user in chip unit, N N m Noise column vector of x 1, M fo Is a training sequence matrix with frequency deviation and dimension N m ×KW;
Get N m = KW, can obtain M fo
Figure C2004100091750003C1
Wherein, Δ f k And phi k Is the frequency offset and initial phase of the kth user, tc is the chip width, i is the position of the first chip of the training sequence in a time slot,
Figure C2004100091750003C2
is a transmitted training sequence;
step 12, using the frequency offset estimated by the previous subframe as the frequency offset estimation of the present subframe to obtain the local training sequence matrix estimation with frequency offset
Figure C2004100091750003C3
Figure C2004100091750003C4
Wherein
Figure C2004100091750003C5
Is the frequency offset 211 of the kth user estimated in the previous subframe;
step 13, utilize
Figure C2004100091750003C6
Estimated channel
Figure C2004100091750003C7
Figure C2004100091750003C8
3. The method of claim 2, wherein the frequency offset is estimated and corrected if the frequency offset of each user is Δ f k If the phase difference is very small, then it will beAs an approximationForming a circulant matrix to reduce the channel using Fourier transform
Figure C2004100091750003C10
The amount of computation of the term.
4. The method of claim 2, wherein the step 2 comprises the steps of:
step 21, overlapping the data signals with frequency offset of all users in a time slot after passing through a multipath channel to obtain the user data signal of the time slot:
e d =A fo d+n
wherein e is d Is (N) s Q + W-1) x 1, ns is the number of user symbols, Q is the spreading factor, d is the column vector of all user symbols, and the dimension is N s K × 1,n is N s K × 1 noisy column vector, A fo Is a system transmission matrix with dimensions of (N) s Q+W-1)×N s K;
A fo The structure of (a) is as follows:
Figure C2004100091750004C1
A fo is a diagonal array of blocks, block V fo (n) Can be expressed as:
Figure C2004100091750004C2
wherein, b k,fo (n) Is the convolution of the spreading scrambling code composite sequence of the k user with frequency offset and the channel impulse response of the k user, the dimension is (Q + W-1) x 1, and is expressed as:
Figure C2004100091750004C3
Figure C2004100091750004C4
wherein, c fo k Is a composite sequence of spread spectrum scrambling codes of the k user with frequency deviation, and the dimension is(Q+W-1)×W,h k Is the channel impulse response of the kth user, with dimension W x 1,
locally generated spreading scrambling code composite sequence c of k user k Expressed as:
Figure C2004100091750005C1
step 22, using the frequency offset estimated by the previous subframe as the frequency offset estimation of the present subframe, thereby generating the frequency offset-based composite sequence estimation of the spreading and scrambling codes
Figure C2004100091750005C2
Figure C2004100091750005C3
Wherein, the first and the second end of the pipe are connected with each other,
Figure C2004100091750005C4
is the frequency offset estimation of the k user of the previous subframe;
step 23, calculated from step 13Extracting the channel impulse response estimation of each user
Figure C2004100091750005C6
K =1, …, K, and
Figure C2004100091750005C7
and
Figure C2004100091750005C8
synthesis of
Figure C2004100091750005C10
And systematic transmission matrix estimation with frequency offset
Figure C2004100091750006C1
k=1,…,K,n=1,…,N s
Figure C2004100091750006C2
n=1,…,N s
Minimum mean square error, MMSE, estimation of all user symbols
Figure C2004100091750006C4
Comprises the following steps:
Figure C2004100091750006C5
wherein R is n Is a noise correlation matrix, K d Is the correlation matrix of the user symbols.
5. A frequency deviation estimation and correction method in TD-SCDMA system is used for estimating and correcting the frequency deviation existing between the received signal carrier and the local carrier, and is characterized in that the method comprises the following steps:
step 1, receiving a chip-level signal;
step 2, extracting training sequence e from the signal of the code chip level m And a user data part e d
Step 3, generating a local training sequence matrix M:
Figure C2004100091750006C6
step 4, estimating the frequency offset of the last subframe
Figure C2004100091750006C7
K =1, …, K, as the frequency offset estimation of the subframe;
step 5, utilizing local training sequence matrix M and
Figure C2004100091750006C8
k =1, …, K, producing local training with frequency offsetTraining sequence matrix estimation
Figure C2004100091750007C1
Figure C2004100091750007C2
Step 6, training sequence e is utilized m And
Figure C2004100091750007C3
performing channel estimation to obtain a channel estimation value
Figure C2004100091750007C4
Figure C2004100091750007C5
Step 7, generating a local spread spectrum scrambling code composite sequence matrix c k ,k=1,…,K:
Figure C2004100091750007C6
Step 8, utilizing c k K =1, …, K and
Figure C2004100091750007C7
k= 1 …, K, produces a local spread spectrum scrambling code composite sequence matrix estimate with frequency offset
Figure C2004100091750007C8
k=1,…,K,n=1,…,N s
Figure C2004100091750008C1
Step 9, from
Figure C2004100091750008C2
Extracting channel estimation value of each user
Figure C2004100091750008C3
K =1, …, K; by usingk=1,…,K,n=1,…,N s And
Figure C2004100091750008C5
k =1, …, K production
Figure C2004100091750008C6
k=1,…,K,n=1,…,N s
Figure C2004100091750008C7
k=1,…,K,n=1,…,N s
Step 10, utilizing
Figure C2004100091750008C8
k=1,…,K,n=1,…,N s Composition V fo (n) ,n=1,…,N s
Figure C2004100091750008C9
n=1,…,N s
Step 11, utilizing V fo (n) ,n=1,…,N s Forming system transmission matrix with frequency deviation
Figure C2004100091750008C10
Figure C2004100091750008C11
Step 12 of utilizing
Figure C2004100091750008C12
And e d Performing joint detection to estimate user symbols
Figure C2004100091750008C13
Step 13, utilize
Figure C2004100091750009C2
Performing frequency offset estimation to estimate the frequency offset difference between the current subframe and the previous subframe
Step 14, use
Figure C2004100091750009C4
K =1, …, K pairs of user symbols
Figure C2004100091750009C5
Carrying out frequency offset correction and outputting a corrected symbol;
step 15, adding the obtained frequency deviation and the estimated frequency deviation of the previous subframe to be used as the estimated frequency deviation of the next subframe:
Figure C2004100091750009C6
k=1,…,K。
CNB2004100091755A 2004-06-04 2004-06-04 Frequency deviation estimation and correction method in TD-SCDMA system Expired - Fee Related CN100367690C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2004100091755A CN100367690C (en) 2004-06-04 2004-06-04 Frequency deviation estimation and correction method in TD-SCDMA system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2004100091755A CN100367690C (en) 2004-06-04 2004-06-04 Frequency deviation estimation and correction method in TD-SCDMA system

Publications (2)

Publication Number Publication Date
CN1705252A CN1705252A (en) 2005-12-07
CN100367690C true CN100367690C (en) 2008-02-06

Family

ID=35577698

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2004100091755A Expired - Fee Related CN100367690C (en) 2004-06-04 2004-06-04 Frequency deviation estimation and correction method in TD-SCDMA system

Country Status (1)

Country Link
CN (1) CN100367690C (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101087159B (en) * 2006-06-05 2011-03-09 中兴通讯股份有限公司 A method for frequency deviation estimation
CN101304587B (en) * 2007-05-09 2012-11-21 华为技术有限公司 Association detection method, cognitive system, subscriber terminal and base station
CN101102122B (en) * 2007-07-31 2010-06-16 北京天碁科技有限公司 A frequency deviation estimation device and its method
CN101312442B (en) * 2008-05-22 2010-12-29 北京创毅视讯科技有限公司 Estimation method and estimation unit for frame head and integral frequency bias of receiver frame head
CN101958744B (en) * 2009-07-13 2013-11-20 电信科学技术研究院 Method and device for pre-calibrating frequency deviation
CN102710270B (en) * 2012-06-15 2014-04-30 深圳国人通信有限公司 Radio remote unit and automatic frequency pre-controlling method and device thereof
US9584355B2 (en) * 2013-03-25 2017-02-28 Intel Deutschland Gmbh Device and method for frequency offset estimation
CN107979547A (en) * 2016-10-24 2018-05-01 深圳超级数据链技术有限公司 Channel estimation methods and device
CN107979549A (en) * 2016-10-24 2018-05-01 深圳超级数据链技术有限公司 Channel estimation methods and device
CN109525336A (en) * 2018-10-29 2019-03-26 上海大学 Based on the radio communication channel test method of frequency deviation measurement time domain compensation under asynchronous clock

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1160981A2 (en) * 2000-05-30 2001-12-05 Nokia Mobile Phones Ltd. Method and arrangement for reducing frequency offset in a radio receiver
CN1342004A (en) * 2000-09-05 2002-03-27 华为技术有限公司 Frequency deviation estimation method applicable to code division multiple address access system
CN1486008A (en) * 2002-09-24 2004-03-31 深圳市中兴通讯股份有限公司 Method and apparatus for automatic frequency correction in CDMA system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1160981A2 (en) * 2000-05-30 2001-12-05 Nokia Mobile Phones Ltd. Method and arrangement for reducing frequency offset in a radio receiver
CN1342004A (en) * 2000-09-05 2002-03-27 华为技术有限公司 Frequency deviation estimation method applicable to code division multiple address access system
CN1486008A (en) * 2002-09-24 2004-03-31 深圳市中兴通讯股份有限公司 Method and apparatus for automatic frequency correction in CDMA system

Also Published As

Publication number Publication date
CN1705252A (en) 2005-12-07

Similar Documents

Publication Publication Date Title
CN101341663B (en) Linear turbo equalization using despread values
US20060233225A1 (en) Frequency synchronization apparatus and frequency synchronization method
KR101520362B1 (en) Symbol estimation methods and apparatuses
CN101888259B (en) Method, apparatus and system used for CDMA communication
KR100552076B1 (en) Signal receiving device in CDMA communication system
JP4271235B2 (en) Method and apparatus for frequency domain equalization in a single carrier receiver
CN100367690C (en) Frequency deviation estimation and correction method in TD-SCDMA system
US20040131109A1 (en) Bidirectional turbo ISI canceller-based DSSS receiver for high-speed wireless LAN
JP2003037526A (en) Device and method of spread code synchronizing detection of spread spectrum signals
CN103178947B (en) The relevant carrier synchronization method with coding auxiliary based on time domain
JP3377389B2 (en) Signal receiving method and apparatus in spread spectrum wireless communication system
CA2448258A1 (en) Channel estimation in cdma communications systems using both lower power pilot channel and higher power data channel
CN113141195A (en) Demodulation method for dispreading in direct sequence spread spectrum system and storage medium
US20040252793A1 (en) Channel estimation method and apparatus
US20130322507A1 (en) Phase Tracking in Communications Systems
CN102035770A (en) Method for estimating channel by means of correlation
AU2005203016A1 (en) Method of noise factor computation for a chip equaliser in a spread spectrum receiver
US6728306B1 (en) Method and apparatus for synchronizing a DS-CDMA receiver
CN101330482A (en) Method for implementing rapid channel estimation using FWT in DMB-T system
CN101116257A (en) Apparatus and method for determining a correlation value
KR100395384B1 (en) Demodulation of receiver with simple structure
JP2000196498A (en) Digitla communication system, its transmitter and receiver and frame synchronization detection circuit
JP3676986B2 (en) Radio receiving apparatus and radio receiving method
CN109561042A (en) A kind of timing frequency synchronous method of ofdm system receiver
TW200522623A (en) Reduced complexity sliding window based equalizer

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
ASS Succession or assignment of patent right

Owner name: JINJIANG HIGH-TECH DEVELOPMENT OFFICE

Free format text: FORMER OWNER: ZTE CORPORATION

Effective date: 20131111

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 518057 SHENZHEN, GUANGDONG PROVINCE TO: 362200 QUANZHOU, FUJIAN PROVINCE

TR01 Transfer of patent right

Effective date of registration: 20131111

Address after: Jinjiang City, Fujian province 362200 Chongde road Quanzhou City No. 283 friends Jade Technology Building

Patentee after: Jinjiang high and new technology development office

Address before: 518057 Nanshan District, Guangdong high tech Industrial Park, science and Technology Industrial Park, ZTE building, block A, layer 6, layer

Patentee before: ZTE Corporation

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20080206

Termination date: 20170604

CF01 Termination of patent right due to non-payment of annual fee