CN101257324B  Linear combined channel estimation method in TDSCDMA system  Google Patents
Linear combined channel estimation method in TDSCDMA system Download PDFInfo
 Publication number
 CN101257324B CN101257324B CN2007100377239A CN200710037723A CN101257324B CN 101257324 B CN101257324 B CN 101257324B CN 2007100377239 A CN2007100377239 A CN 2007100377239A CN 200710037723 A CN200710037723 A CN 200710037723A CN 101257324 B CN101257324 B CN 101257324B
 Authority
 CN
 China
 Prior art keywords
 channel
 estimation
 matrix
 scdma
 channel estimation
 Prior art date
Links
 239000011159 matrix materials Substances 0.000 claims description 28
 280000867207 Lambda companies 0.000 claims description 10
 230000000875 corresponding Effects 0.000 claims description 10
 239000000562 conjugates Substances 0.000 claims description 6
 239000000284 extracts Substances 0.000 claims description 6
 238000004364 calculation methods Methods 0.000 claims description 5
 280000987480 New Channel companies 0.000 claims description 3
 230000015572 biosynthetic process Effects 0.000 claims description 3
 125000004122 cyclic group Chemical group 0.000 claims description 3
 238000000605 extraction Methods 0.000 claims description 3
 238000005755 formation reactions Methods 0.000 claims description 3
 238000000638 solvent extraction Methods 0.000 claims description 3
 230000017105 transposition Effects 0.000 claims description 3
 238000000034 methods Methods 0.000 description 9
 238000010586 diagrams Methods 0.000 description 4
 238000005516 engineering processes Methods 0.000 description 4
 230000005540 biological transmission Effects 0.000 description 3
 230000000694 effects Effects 0.000 description 2
 238000000926 separation method Methods 0.000 description 2
 238000004088 simulation Methods 0.000 description 2
 102000014150 Interferons Human genes 0.000 description 1
 108010050904 Interferons Proteins 0.000 description 1
 280000874760 Third Generation companies 0.000 description 1
 238000004422 calculation algorithm Methods 0.000 description 1
 238000004891 communication Methods 0.000 description 1
 238000005562 fading Methods 0.000 description 1
 229940079322 interferon Drugs 0.000 description 1
 230000004048 modification Effects 0.000 description 1
 238000006011 modification reactions Methods 0.000 description 1
 238000001228 spectrum Methods 0.000 description 1
 230000001629 suppression Effects 0.000 description 1
 230000001360 synchronised Effects 0.000 description 1
Abstract
Description
Technical field
The present invention relates to the channel estimation methods of receiver in a kind of wireless communication system, particularly the linear combined channel estimation method in a kind of TDSCDMA system.
Background technology
Next TDSCDMA system uplink data will carry out channel estimating and jointdetection after the base station receives by antenna.Associated detection technique then refers to and takes full advantage of multiple access interference (MA I), a kind of signal separation techniques of within the step all users' signal all being separated.The very suitable use associated detection technique of TDSCDMA system, because code channel number is limited and be confined on the most suitable number in this system, and all users in the same time slot are synchronous.Associated detection technique has become the focus in the 3G (Third Generation) Moblie technology.The associated detection technique that adopts in the TDSCDMA system is on the basis of traditional detection technology, take full advantage of and cause all subscriber signals that multiple access disturbs and the prior information of multipath thereof, the separation of subscriber signal being used as the joint detection process that is mutually related of a unification finishes, thereby has good interference free performance, reduced the requirement of system to the power control precision, therefore can effectively utilize the uplink spectrum resource more, improve power system capacity significantly.
Channel estimating is a vital part in joint detection process, and it will influence the testing process performance in many aspects.The accuracy of the detection information of combined detector output is directly related with the quality of channel estimating, and how obtaining accurately efficiently, each user's channel impulse response is the key that followup base band algorithm is realized.As a rule, channel estimating is more accurate, just can obtain more good systematic function.
In the TDSCDMA system under the identical networking situation, the receiver channel method of estimation generally has two big class, i.e. single cell channel estimation method and multiple cell combined channel estimation methods.
Single cell channel estimation method is exactly to continue to use Steiner estimator in the past, will be used as to disturb from the signal of other residential quarter and handle, and obviously its performance seriously is subject to interference, usually can't operate as normal in the TDSCDMA of identical networking system.
Summary of the invention
The purpose of this invention is to provide the linear combined channel estimation method in a kind of TDSCDMA system, to improve performance for estimating channel, make it to be applicable to the TDSCDMA system of identical networking.
For solving the problems of the technologies described above, this multiple cell combined channel estimation method is applied to channel estimating with multiuser detection, to be used as different users from the signal of different districts, and adopt suitable multiuser detection partly to eliminate the phase mutual interference between the different districts.Technical scheme of the present invention is:
Linear combined channel estimation method in a kind of TDSCDMA system is characterized in that this method may further comprise the steps:
(1) the realtime signal that receives K residential quarter from receiver;
(2) data splitting block is isolated the reception signal x of P chip of corresponding intermediate code part from the reception signal r of a time slot, as the realtime input of channel estimator;
(3) basic midamble code according to K the residential quarter that participates in the combined channel estimation constructs a P * KP rank complex matrix The cyclic shift of classifying k residential quarter basic midamble code as;
(4) the Fourier transform FFT{M of the basic midamble code of K residential quarter of calculating participation combined channel estimation _{k}(:, 1) }, M wherein _{k}(:, 1) representing matrix M _{k}First row, i.e. basic midamble code itself, being expressed as matrix form is Λ _{k}=[λ _{K, 1}λ _{K, P}] ^{T}=FFT{M _{k}(:, 1), p=1 ..., P, storage result of calculation Λ _{k}Fixedly input value as channel estimator;
(5) make S=MM
^{H}, characteristic value Λ=[λ of calculating S
_{1}λ
_{P}]
^{T}, wherein
^{H}The expression complexconjugate transpose,
^{T}The expression transposition, computational methods are
(6) calculate the Fourier transform FFT{x} that receives signal x;
(7) calculate the least square solution of estimation of the channel impulse response of each residential quarter
(8) based on the result of abovementioned channel estimating, selected tap is used for estimating again;
(9) selected tap forms a new channel impulse response vector h ', new matrix M of row formation of extraction correspondence from matrix M ' and, the least square solution of channel impulse response estimation is
Wherein Expression MoorePenrose pseudoinverse;
(10) h as a result with channel estimating is applied to jointdetection with the data symbol part that receives signal r, to obtain sending the estimation of data symbol.
The criterion of the selected tap when estimating again in the described step 9 is based on the hypothesis to the multidiameter delay expansion.
The condition that described M ' for estimation again is a row nonsingular matrix is that the columns of M ' is less than or equal to P.
Described P=128, set basis is referring to 3GPP TS25.221 physical channel and transmission channel mapping (TDD) the 5A.2.3 joint to physical channel.
Described reception signal x can be expressed as
x＝Mh+n
Wherein Be the channel impulse response of k residential quarter, Be noise;
Minimize the Euclidean norm of evaluated error The least square solution that can obtain the estimation of channel impulse response is
Wherein ^{1}Expression is inverted.
Described new matrix M for estimation again ' be a row nonsingular matrix, the least square solution of channel impulse response estimation is
Adopt shortcut calculation in the described estimation again: calculate M ' ^{H}In M ' time, calculate M by the method for fast fourier transform ^{H}M therefrom extracts desired element then; Calculate (M ' ^{H}M ') ^{1}The time, earlier with the partitioning of matrix, carry out Cholesky then and decompose; Calculate M ' ^{H}During x, calculate M by the method for fast fourier transform equally ^{H}X therefrom extracts desired element then.
Compared with prior art, beneficial effect of the present invention is as follows:
In the channel estimation results of technical solution of the present invention, arbitrary residential quarter accesses supression preferably to the Recombinant Interferon of other residential quarter, and the channel estimation window of the user's correspondence in the residential quarter can reflect the multipath environment that this user experiences more exactly.Therefore, compare with single cell channel estimated result, performance for estimating channel of the present invention has greatly improved.
Description of drawings
Fig. 1 is a typical TDSCDMA structure of time slot schematic diagram.
Fig. 2 is channel estimating and jointdetection schematic diagram.
Fig. 3 is the simulation result figure of single cell channel estimation method and multiple cell combined channel estimation method.
Embodiment
The present invention will be further described below in conjunction with drawings and Examples, but should not limit protection scope of the present invention with this.
As shown in Figure 1, provide a typical TDSCDMA structure of time slot schematic diagram.A time slot is made of at interval the intermediate code of the data symbol of two sections 352 chips, 144 chips and the protection of 16 chips, altogether 864 chips.
In general different basic midamble codes is used in different residential quarters, and basic midamble code generates intermediate code and saves referring to 3GPP TS25.221 physical channel and transmission channel mapping (TDD) 5A.2.3 to physical channel.
Consider one and can receive the situation that K residential quarter transmits simultaneously.From receiver, receive the signal of K residential quarter in real time.
Fig. 2 provides a channel estimating and jointdetection schematic diagram, as shown in Figure 2, the reception signal r of a time slot comprises 864 chips, pass through data splitting block, intercepting is the reception signal x of corresponding intermediate code part wherein, namely 496 chips of the 369th chip to the have P chip altogether as the realtime input of channel estimator, get P=128 here.
The reception signal x of corresponding intermediate code part can be expressed as
x＝Mh+n
Wherein Be that M is a P * KP rank complex matrix, The cyclic shift of classifying k residential quarter basic midamble code (Basic Midamble Code) as. Be the channel impulse response of k residential quarter, wherein ^{T}The expression transposition. Be noise.
The least square of the estimation of channel impulse response (LS, Least Square) solution is
Wherein ^{H}The expression complexconjugate transpose, ^{1}Expression is inverted.
Make S=MM
^{H}, M so
_{k}With S all be circular matrix (Circulant Matrices).According to the character of circular matrix, M
_{k}Characteristic value be λ
_{K, p}, p=1 ..., P is expressed as matrix form, i.e. Λ
_{k}=[λ
_{K, 1}λ
_{K, P}]
^{T}=FFT{M
_{k}(:, 1) }, Mk (:, 1) representing matrix M wherein
_{k}First row, i.e. basic midamble code itself.Similarly, the characteristic value of S can be expressed as Λ=[λ
_{1}λ
_{P}]
^{T}, wherein
Store described result of calculation Λ _{k}, Λ is as the fixedly input value of channel estimator.Calculate the Fourier transform FFT{x} that receives signal x.
Therefore the least square solution of the estimation of channel impulse response also can basis
Based on abovementioned channel estimation results, according to certain criterion, for example based on the hypothesis to the multidiameter delay expansion, select some may exist the tap of multipath signal to estimate again.Selected tap can form a new channel impulse response vector h ', new matrix M of row formation of extraction correspondence from matrix M ' and, the least square solution of channel impulse response estimation can be write as so
Wherein Expression MoorePenrose pseudoinverse.
When abovementioned new matrix M ' columns when being less than or equal to P, M ' is a row nonsingular matrix, following formula can be write as so
In the abovementioned process of estimating again, can use some shortcut calculations.For example, calculating M ' ^{H}In M ' time, can calculate M by the method for fast fourier transform ^{H}M therefrom extracts desired element then; Calculate (M ' ^{H}M ') ^{1}The time, can carry out Cholesky then and decompose earlier with the partitioning of matrix; Calculating M ' ^{H}In the time of x, can calculate M by the method for fast fourier transform equally ^{H}X therefrom extracts desired element then.
As shown in Figure 2, the result of channel estimating will be applied to jointdetection, obtain sending the estimation of data symbol with the data symbol part that receives signal.
By an emulation good effect of the technical program is described, simulated conditions as shown in Table 1 and Table 2, set basis is referring to 3GPP TS25.102 subscriber equipment wireless transmission and reception (TDD).
Table 1 simulated conditions
Table 2 multipath fading environments propagation conditions
Based on abovementioned simulated conditions, Fig. 3 has provided the Matlab simulation result figure of single cell channel estimation method and multiple cell combined channel estimation method, and abscissa is tap among the figure, and ordinate is intensity.Two figure in the left side are single cell channel results estimated, and last figure is single cell channel estimated result of residential quarter 1, and figure below is single cell channel estimated result of residential quarter 2; Two results that figure is technical solution of the present invention in the right, last figure is the linear combined channel estimation results in many residential quarters of residential quarter 1, figure below is the multiplot joint linear channel estimated result of residential quarter 2.
Obviously, in the channel estimation results of single residential quarter, residential quarter 2 has been subjected to the interference from residential quarter the last 1, and in the channel estimation results of technical solution of the present invention, interference from residential quarter 1 is restrained preferably, can see the channel estimation window of user's correspondence in the residential quarter 2 more clearlyly, reflect the multipath environment of this user's experience more exactly.Therefore, compare with single cell channel estimated result, the present invention has bigger performance improvement.
The above is preferred embodiment of the present invention only, is not to limit practical range of the present invention.Be that all equivalences of doing according to the content of the present patent application claim change and modification, all should be technology category of the present invention.
Claims (6)
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN2007100377239A CN101257324B (en)  20070228  20070228  Linear combined channel estimation method in TDSCDMA system 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN2007100377239A CN101257324B (en)  20070228  20070228  Linear combined channel estimation method in TDSCDMA system 
Publications (2)
Publication Number  Publication Date 

CN101257324A CN101257324A (en)  20080903 
CN101257324B true CN101257324B (en)  20131009 
Family
ID=39891838
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN2007100377239A CN101257324B (en)  20070228  20070228  Linear combined channel estimation method in TDSCDMA system 
Country Status (1)
Country  Link 

CN (1)  CN101257324B (en) 
Families Citing this family (2)
Publication number  Priority date  Publication date  Assignee  Title 

CN101931450B (en) *  20090619  20120905  联芯科技有限公司  Method for selecting codebooks and device thereof 
CN102739282B (en) *  20110402  20150204  联芯科技有限公司  Multiantenna combined detection system of TDSCDMA system and multiantenna combined detection method thereof 
Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

CN1533191A (en) *  20030321  20040929  乐金电子(中国)研究开发中心有限公司  Adaptive channel evaluating method and device based on expected medium code of construction 
EP1487165A1 (en) *  20030614  20041215  Lg Electronics Inc.  TDSCDMA channel estimation method and apparatus 
CN1770650A (en) *  20041015  20060510  展讯通信（上海）有限公司  Maximum ratio combining of channel estimation for joint detection in TDSCDMA systems 

2007
 20070228 CN CN2007100377239A patent/CN101257324B/en active IP Right Grant
Patent Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

CN1533191A (en) *  20030321  20040929  乐金电子(中国)研究开发中心有限公司  Adaptive channel evaluating method and device based on expected medium code of construction 
EP1487165A1 (en) *  20030614  20041215  Lg Electronics Inc.  TDSCDMA channel estimation method and apparatus 
CN1770650A (en) *  20041015  20060510  展讯通信（上海）有限公司  Maximum ratio combining of channel estimation for joint detection in TDSCDMA systems 
Also Published As
Publication number  Publication date 

CN101257324A (en)  20080903 
Similar Documents
Publication  Publication Date  Title 

Onggosanusi et al.  Performance analysis of closedloop transmit diversity in the presence of feedback delay  
US7336599B2 (en)  Apparatus and method for joint detection receiving irrespective of orthogonal code length in mobile communication system  
KR100669962B1 (en)  Single user detection  
CN100450115C (en)  Method and apparatus for estimation of phase offset between communication channels  
CN100403659C (en)  Family of linear multiuser detectors (MUDS)  
JP4213752B2 (en)  RAKE receiver  
Pätzold  Mobile radio channels  
KR100559070B1 (en)  Adaptive antenna array and method for control thereof  
TWI388135B (en)  Fast joint detection  
TWI408927B (en)  Advanced receiver with sliding window block linear equalizer  
TW386325B (en)  Selfsynchronizing equalization techniques and systems  
TW497340B (en)  Adaptive channel estimation in a wireless communication system  
KR100759297B1 (en)  Low complexity data detection using fast fourier transform of channel correlation matrix  
CN101888259B (en)  Method, apparatus and system used for CDMA communication  
KR100770875B1 (en)  Beam forming apparatus and method using estimating interference power in array antenna system  
JP2004521539A (en)  Method and apparatus for estimating reduced rank channel in a communication system  
TWI260173B (en)  Efficient joint detection  
US7088693B2 (en)  Adaptive diversity combining for Wide Band Code Division Multiple Access (WCDMA) based on iterative channel estimation  
CN100542080C (en)  Associated detecting method and device and use the system of this device  
TWI287364B (en)  Iterative and turbobased method and apparatus for equalization of spreadspectrum downlink channels  
CN101814931B (en)  Doppler frequency shift estimation and compensation method in TDSCDMA (Time DivisionSynchronization Code Division Multiple Access) system  
US20060126753A1 (en)  Apparatus and method for detecting spacetime multiuser signal of base station having array antenna  
TWI250731B (en)  Extended algorithm data estimator  
US7606295B2 (en)  Generalized rake receiver for wireless communication  
US20030142762A1 (en)  Wireless receiver method and apparatus using spacecovertime equalization 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
C06  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
GR01  Patent grant  
C14  Grant of patent or utility model  
TR01  Transfer of patent right 
Effective date of registration: 20200604 Address after: 101318 8 807, building 6, ronghuiyuan, airport economic core area, Xicheng District, Beijing Patentee after: Xin Xin finance leasing (Beijing) Co.,Ltd. Address before: 201203 Shanghai City Songtao road Pudong Zhangjiang hi tech Park No. 696 35 Patentee before: SPREADTRUM COMMUNICATIONS (SHANGHAI) Co.,Ltd. 